Explict vs Implict FInite Difference Scheme

Explicit methods calculate the state of a system at a later time from the state of the system at the current time, whileimplicit methods find a solution by solving an equation involving both the current state of the system and the later one. Mathematically, if is the current system state and is the state at the later time ( is a small time step), then, for an explicit method while for an implicit method one solves an equation to find http://earthscience.stackexchange.com/questions/2253/what-are-the-differences-between-implicit-and-explicit-methods-when-applied-to-m

http://www.cfd-online.com/Forums/main/5982-implicit-vs-explicit-method.html

du/dt = f(u)
where u is a vector and f is a known function. Then we may approximate the time derivative over a time step h as
(implicit) u(n+1) = u(n) + h.f(u(n+1))

if f(u) = Au where A is a constant, then u(n+1）=u(n)+Ahu(n+1),then, u(n+1)=u(n)/(1-Ah)

https://www.flow3d.com/home/resources/cfd-101/numerical-issues/implicit-versus-explicit-numerical-methods

y i+1 = f(xi, yi), yi+2 = f(xi,xi+1,y i,yi+1), etc.In other words, your solution proceeds by solving explicitly for a new unknown value inthe solution array, given all previous values in the array. On the other hand, implicit
methods imply the simultaneous solution of n linear algebraic equations that provide, atonce, the elements of the solution array. 在上面的式子中，我们并不能直接求解出yi,这时候需要更多的方程来求解，这些方程可以写成矩阵来求解。 1-s2.0-S0096300304004680-main.pdf role_959843.pdf tcs-30_web.pdf

In finite-difference methods, calculations are performed on a grid placed over the interesting flow domain into the x-tplane. This computational grid is defined by some equal or variable space and time steps, Δx and Δt, respectively. A network of discrete points is thus obtained and the flow variables Q, z (or Q, h; V, z; V, h) are derived only for this finite number of grid points. Figure 6.2 shows a typical computational grid. The spatial positions of the grid points are denoted by index jand the time moments by index i. Figure 6.2. Finite-difference computational grid.

As principle, finite-difference methods transform governing partial differential equations into a set of algebraic (linear or nonlinear) finite-difference equations, which are solved to allow the values of flow variables in a grid point or in all grid points on a time line. These finite-difference equations are derived by approximating the time and space derivatives with some finite-difference expressions. Not only the derivatives, but also the other terms into the Saint-Vénant equations must be defined in a certain manner.

Concerning the space derivative of a continuous function f(x,t) at time moment ti and space position xj on the grid in Figure 6.2, this can be approximated as: a forward difference approximation; a backward difference relationship; a central difference expression, (6.29)

in which fji represents the value of f at grid point (xj ,ti).

In an analogous manner, the time derivative may be defined in several different ways as, for example: ;  (6.30)

Usually, the value of the function f(x,t) is accepted at grid point (xj,ti) as fji , but some different approximations can also be decided.

A finite difference method must employ a certain type of finite-difference scheme. These schemes are grouped in two major classes: explicit and implicit finite-difference schemes.

Explicit schemes are those in which the flow variables at any point j and time level i+1 may be computed using only known data at a few adjacent points on the time line i. These schemes do not lead to a system of algebraic equations, but rather to only two finite difference equations for each grid point (xj,ti+1). By solving the two equations the unknown values of flow variables are obtained and then, the computation proceeds to the next grid point along the time line i+1.

In implicit schemes, finite-difference expressions used to approximate the space and/or time derivatives at grid point (xj, ti+1) include the unknown values of flow variables at a few adjacent points on the time line i+1. Consequently, a system of algebraic equations is produced for a given time line i+1 and by solving this system, all unknown values are simultaneously determined at time level i+1.

Replacing the continuous original problem with an integration over a discrete computational grid introduces numerical errors into the results. A finite-difference scheme is stable if such errors are not amplified during computation from one time level to the next. The numerical stability depends on the size of the time and space steps and on some flow characteristics. The Courant condition: (6.31)

is a necessary but insufficient condition for stability of an explicit scheme.

Therefore, any explicit scheme is conditionally stable, the Courant condition requiring to work with small time steps as compared with the physical phenomena evolution. Despite their computational simplicity, the explicit methods are seldom used in river modelling for reasons of this stability restriction.

On the other hand, the implicit schemes may appear more complicated, but can generally be made unconditionally stable for large computation steps and with little loss of accuracy.

Another distinction among schemes belonging to the same class is related to the way in which the non derivative terms (such as Sf(Q,h), B(h), A(h) etc.) are discretized. Because these terms are functions of dependent variables, their treatment induces a linear or nonlinear feature of the finite-difference equations.

A lot of finite-difference schemes will be presented in more details within the next sections.

Numerical solution of kinematic wave model

First form (6.21) of the kinematic wave model is used to illustrate an explicit linear scheme on a finite-difference cell as show in Figure 6.4.

A backward finite-difference is used to approximate the space derivative, while the time derivative is usually expressed at the same xj+1 spatial position. In order to obtain a linear scheme, the nonderivative term αβQβ-1 must be evaluated for a known Q value, which here is accepted as mean between the two diagonal valuesQji+1 and Qij+1. Figure 6.4. Finite-difference cell for linear kinematic wave solution.

Consequently, the finite-difference form of equation 6.21 is: (6.40)

and the unknown Qj+1i+1 results as: (6.41)

Starting with the inflow hydrograph value Q at the time level ti+1, the computation sequentially proceeds from upstream towards downstream grid points, along the current time line i+1.

However, this scheme supposes that within the coefficient α, the wetted perimeter P remains constant. On the other hand, by using Q instead of A into the time derivative ∂A/t, relative computational errors are decreased. Indeed, taking the logarithm of equation 6.20, i.e. ,

and differentiating, one obtains: ,

where β is 0.6 for the Manning’s equation. It follows that the estimation error in Q would be magnified by 1/0.6 ≈ 1.67 if the cross-sectional area A is used as dependent variable.

A finite-difference method may employ either an explicit scheme or an implicit scheme for solution. The main difference between the two is that in the explicit method, the unknown values are solved sequentially along a time line from one distance point to the next, while in the implicit method the unknown valueson a given time line are all determined simultaneously. The explicit method is simpler but can be unstable, which means that small values of Ax and A^are required for convergence of the numerical procedure.

Updated: 2016-08-16 — am11:34