The group-preserving scheme developed by Liu (2001) for calculating the solutions of k-dimensional differential equations system adopted the Cayley transform to formulate the Lie group from its Lie algebra A ∈ so(k,1). In this paper we consider a more effective exponential mapping to derive exp(hA). In order to overcome the difficulty of numerical instabilities encountered by employing group-preserving schemes on stiff differential equations, we further combine the nonstandard finite difference method into the group-preserving schemes to obtain unconditional stable numerical methods. They provide single-step explicit time integrators for stiff differential equations. Several numerical examples are examined, some of which are compared with exact solutions showing that the nonstandard group-preserving schemes have good computational efficiency and certain accuracy.
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.