Explicit and implicit methods: Difference between revisions

Revision as of 12:29, 14 January 2021

Latest revision as of 20:31, 4 January 2025

(6 intermediate revisions by 5 users not shown)

Line 1:

{{Short description|Approaches for approximating solutions to differential equations}}

{{more citations needed|date=December 2009}}

'''Explicit and implicit methods''' are approaches used in [[numerical analysis]] for obtaining numerical approximations to the solutions of time-dependent [[ordinary differential equation|ordinary]] and [[partial differential equation]]s, as is required in [[computer simulation]]s of [[Process (science)|physical processes]]. ''Explicit methods'' calculate the state of a system at a later time from the state of the system at the current time, while ''implicit methods'' find a solution by solving an equation involving both the current state of the system and the later one. Mathematically, if <math>Y(t)</math> is the current system state and <math>Y(t+\Delta t)</math> is the state at the later time (<math>\Delta t</math> is a small time step), then, for an explicit method

Line 6 ⟶ 7:

to find <math>Y(t+\Delta t).</math>

== Computation ==

Implicit methods require an extra computation (solving the above equation), and they can be much harder to implement. Implicit methods are used because many problems arising in practice are [[Stiff equation|stiff]], for which the use of an explicit method requires impractically small time steps <math>\Delta t</math> to keep the error in the result bounded (see [[numerical stability]]). For such problems, to achieve given accuracy, it takes much less computational time to use an implicit method with larger time steps, even taking into account that one needs to solve an equation of the form (1) at each time step. That said, whether one should use an explicit or implicit method depends upon the problem to be solved.

Line 11 ⟶ 13:

:<math>Y(t+\Delta t) = F(Y(t+\Delta t))+G(Y(t)),\,</math>

while one is treated explicitly and the other implicitly.

For usual applications the implicit term is chosen to be linear while the explicit term can be nonlinear. This combination of the former method is called ''Implicit-Explicit Method'' (short IMEX,<ref>U.M. Ascher, S.J. Ruuth, R.J. Spiteri: ''[http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.48.1525&rep=rep1&type=pdf Implicit-Explicit Runge-Kutta Methods for Time-Dependent Partial Differential Equations]'', Appl Numer Math, vol. 25(2-3), 1997</ref><ref>L.Pareschi, G.Russo: ''[https://www.researchgate.net/profile/Lorenzo_Pareschi/publication/230865813_Implicit-Explicit_Runge-Kutta_schemes_for_stiff_systems_of_differential_equations/links/0046352a03ba3ee92a000000.pdf Implicit-Explicit Runge-Kutta schemes for stiff systems of differential equations]'', Recent Trends in Numerical Analysis, Vol. 3, 269-289, 2000</ref><ref>Sebastiano Boscarino, Lorenzo Pareschi, and Giovanni Russo: ''Implicit-Explicit Methods for Evolutionary Partial Differential Equations'', SIAM, ISBN 978-1-61197-819-3 (2024).</ref>).

==Illustration using the forward and backward Euler methods==

Line 21 ⟶ 23:

;Forward Euler method:

[[File:Result of applying integration schemes.png|thumb|The result of applying different integration methods to odethe ODE: <math> y'=-y^2, \; t\in[0, 5], \; y_0=1 </math> with <math>\Delta t = 5/10</math>.]]

The forward [[Euler method]]

:<math>\left(\frac{dy}{dt}\right)_k \approx \frac{y_{k+1}-y_k}{\Delta t} = - y_k^2</math>

Line 41 ⟶ 43:

In the vast majority of cases, the equation to be solved when using an implicit scheme is much more complicated than a quadratic equation, and no analytical solution exists. Then one uses [[root-finding algorithm]]s, such as [[Newton's method]], to find the numerical solution.

;Crank -Nicolson method:

With the [[Crank-Nicolson method]]

:<math>\frac{y_{k+1}-y_k}{\Delta t} = -\frac{1}{2}y_{k+1}^2 -\frac{1}{2}y_{k}^2</math>

one finds the implicit equation

: <math>y_{k+1}+\frac{1}{2}{\Delta t} y_y^{2}_{k+1}^2 = y_k - \frac{1}{2}\Delta t y_{k}^2</math>

for <math>y_{k+1}</math> (compare this with formula (3) where <math>y_{k+1}</math> was given explicitly rather than as an unknown in an equation). This can be numerically solved using [[root-finding algorithm]]s, such as [[Newton's method]], to obtain <math>y_{k+1}</math>.

Crank -Nicolson can be viewed as a form of more general IMEX (''Im''plicit-''Ex''plicit) schemes.

;Forward-Backward Euler method:

[[File:Comparison_between_Foward-Backward-Euler_and_Foward-Euler.png|thumb|400px|The result of applying the both, the Forward Euler method as well asand the Forward-Backward Euler method for <math>a = 5</math> and <math>n = 30</math>.]]

In order to apply the IMEX-scheme, consider a slightly different differential equation:

: <math>\frac{dy}{dt} = y-y^2, \ t\in [0, a]\quad \quad (5)</math>