×

An almost symmetric Strang splitting scheme for nonlinear evolution equations. (English) Zbl 1368.65074

Summary: In this paper we consider splitting methods for the time integration of parabolic and certain classes of hyperbolic partial differential equations, where one partial flow cannot be computed exactly. Instead, we use a numerical approximation based on the linearization of the vector field. This is of interest in applications as it allows us to apply splitting methods to a wider class of problems from the sciences.
However, in the situation described, the classic Strang splitting scheme, while still being a method of second order, is not longer symmetric. This, in turn, implies that the construction of higher order methods by composition is limited to order three only. To remedy this situation, based on previous work in the context of ordinary differential equations, we construct a class of Strang splitting schemes that are symmetric up to a desired order.
We show rigorously that, under suitable assumptions on the nonlinearity, these methods are of second order and can then be used to construct higher order methods by composition. In addition, we illustrate the theoretical results by conducting numerical experiments for the Brusselator system and the KdV equation.

MSC:

65J08 Numerical solutions to abstract evolution equations
35F20 Nonlinear first-order PDEs
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs

Software:

RODAS

References:

[1] McLachlan, R. I., On the numerical integration of ordinary differential equations by symmetric composition methods, SIAM J. Sci. Comput., 16, 1, 151-168 (1995) · Zbl 0821.65048
[2] Cheng, C. Z.; Knorr, G., The integration of the Vlasov equation in configuration space, J. Comput. Phys., 22, 3, 330-351 (1976)
[3] Einkemmer, L.; Ostermann, A., Convergence analysis of Strang splitting for Vlasov-type equations, SIAM J. Numer. Anal., 52, 140-155 (2014) · Zbl 1297.65106
[4] Klein, C.; Roidot, K., Fourth order time-stepping for Kadomtsev-Petviashvili and Davey-Stewartson equations, SIAM J. Sci. Comput., 33, 6, 3333-3356 (2011) · Zbl 1298.65141
[5] Hairer, E.; Lubich, C.; Wanner, G., Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations (2006), Springer Verlag · Zbl 1094.65125
[7] Jahnke, T.; Lubich, C., Error bounds for exponential operator splittings, BIT Numer. Math., 40, 4, 735-744 (2000) · Zbl 0972.65061
[8] Holden, H.; Lubich, C.; Risebro, N. H., Operator splitting for partial differential equations with Burgers nonlinearity, Math. Comp., 82, 281, 173-185 (2013) · Zbl 1260.35184
[10] Lubich, C., On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations, Math. Comp., 77, 2141-2153 (2008) · Zbl 1198.65186
[11] Blanes, S.; Casas, F., On the necessity of negative coefficients for operator splitting schemes of order higher than two, Appl. Numer. Math., 54, 1, 23-37 (2005) · Zbl 1071.65097
[12] Hansen, E.; Ostermann, A., High order splitting methods for analytic semigroups exist, BIT Numer. Math., 49, 3, 527-542 (2009) · Zbl 1176.65066
[13] Castella, F.; Chartier, P.; Descombes, S.; Vilmart, G., Splitting methods with complex times for parabolic equations, BIT Numer. Math., 49, 3, 487-508 (2009) · Zbl 1180.65106
[14] Hairer, E.; Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems (1996), Springer Verlag · Zbl 0859.65067
[15] Kametaka, Y., Korteweg-de Vries equation, I. Global existence of smooth solutions, Proc. Japan Acad. Ser. A Math. Sci., 45, 7, 552-555 (1969) · Zbl 0194.41801
[16] Bona, J. L.; Smith, R., The initial-value problem for the Korteweg-de Vries equation, Phil. Trans. R. Soc. A, 278, 1287, 555-601 (1975) · Zbl 0306.35027
[17] Klein, C.; Sparber, C.; Markowich, P., Numerical study of oscillatory regimes in the Kadomtsev-Petviashvili equation, J. Nonlinear Sci., 17, 5, 429-470 (2007) · Zbl 1128.37043
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.