Strang splitting

In applied mathematics Strang splitting is a numerical method for solving differential equations that are decomposable into a sum of differential operators. It is named after Gilbert Strang. It is used to speed up calculation for problems involving operators on very different time scales, for example, chemical reactions in fluid dynamics, and to solve multidimensional partial differential equations by reducing them to a sum of one-dimensional problems.

Fractional step methods

As a precursor to Strang splitting, consider a differential equation of the form

d

y

d

t

=

L

1

(

y

)

+

L

2

(

y

)

{\displaystyle {\frac {d{y}}{dt}}=L_{1}({y})+L_{2}({y})}

where

L

1

{\displaystyle L_{1}}

,

L

2

{\displaystyle L_{2}}

are differential operators. If

L

1

{\displaystyle L_{1}}

and

L

2

{\displaystyle L_{2}}

were constant coefficient matrices, then the exact solution to the associated initial value problem would be

y

(

t

)

=

e

(

L

1

+

L

2

)

t

y

0

{\displaystyle y(t)=e^{(L_{1}+L_{2})t}y_{0}}

.

If

L

1

{\displaystyle L_{1}}

and

L

2

{\displaystyle L_{2}}

commute, then by the exponential laws this is equivalent to

y

(

t

)

=

e

L

1

t

e

L

2

t

y

0

{\displaystyle y(t)=e^{L_{1}t}e^{L_{2}t}y_{0}}

.

If they do not, then by the Baker–Campbell–Hausdorff formula it is still possible to replace the exponential of the sum by a product of exponentials at the cost of a second order error:

e

(

L

1

+

L

2

)

t

y

0

=

e

L

1

t

e

L

2

t

y

0

+

O

(

t

2

)

{\displaystyle e^{(L_{1}+L_{2})t}y_{0}=e^{L_{1}t}e^{L_{2}t}y_{0}+{\mathcal {O}}(t^{2})}

.

This gives rise to a numerical scheme where one, instead of solving the original initial problem, solves both subproblems alternating:

y

~

1

=

e

L

1

Δ

t

y

0

{\displaystyle {\tilde {y}}_{1}=e^{L_{1}\Delta t}y_{0}}

y

1

=

e

L

2

Δ

t

y

~

1

{\displaystyle y_{1}=e^{L_{2}\Delta t}{\tilde {y}}_{1}}

y

~

2

=

e

L

1

Δ

t

y

1

{\displaystyle {\tilde {y}}_{2}=e^{L_{1}\Delta t}y_{1}}

y

2

=

e

L

2

Δ

t

y

~

2

{\displaystyle y_{2}=e^{L_{2}\Delta t}{\tilde {y}}_{2}}

etc.

In this context,

e

L

1

Δ

t

{\displaystyle e^{L_{1}\Delta t}}

is a numerical scheme solving the subproblem

d

y

d

t

=

L

1

(

y

)

{\displaystyle {\frac {d{y}}{dt}}=L_{1}({y})}

to first order. The approach is not restricted to linear problems, that is,

L

1

{\displaystyle L_{1}}

can be any differential operator.