Beam and Warming scheme

In numerical mathematics, Beam and Warming scheme or Beam–Warming implicit scheme introduced in 1978 by Richard M. Beam and R. F. Warming,[1][2] is a second order accurate implicit scheme, mainly used for solving non-linear hyperbolic equations. It is not used much nowadays.

Introduction

edit

This scheme is a spatially factored, non iterative, ADI scheme and uses implicit Euler to perform the time Integration. The algorithm is in delta-form, linearized through implementation of a Taylor-series. Hence observed as increments of the conserved variables. In this an efficient factored algorithm is obtained by evaluating the spatial cross derivatives explicitly. This allows for direct derivation of scheme and efficient solution using this computational algorithm. The efficiency is because although it is three-time-level scheme, but requires only two time levels of data storage. This results in unconditional stability. It is centered and needs the artificial dissipation operator to guarantee numerical stability.[1]

The delta form of the equation produced has the advantageous property of stability (if existing) independent of the size of the time step.[3]

The method

edit
 

Consider the inviscid Burgers' equation in one dimension

 

Burgers' equation in conservation form,

 

where : 

Taylor series expansion

edit
 

The expansion of : 

 

This is also known as the trapezoidal formula.

 
 

Note that for this equation,

 

Tri-diagonal system

edit

Resulting tri-diagonal system:

 

This resulted system of linear equations can be solved using the modified tridiagonal matrix algorithm, also known as the Thomas algorithm.[4]

Dissipation term

edit

Under the condition of shock wave, dissipation term is required for nonlinear hyperbolic equations such as this. This is done to keep the solution under control and maintain convergence of the solution.

 

This term is added explicitly at level   to the right hand side. This is always used for successful computation where high-frequency oscillations are observed and must be suppressed.

Smoothing term

edit

If only the stable solution is required, then in the equation to the right hand side a second-order smoothing term is added on the implicit layer. The other term in the same equation can be second-order because it has no influence on the stable solution if

 

The addition of smoothing term increases the number of steps required by three.

Properties

edit

This scheme is produced by combining the trapezoidal formula, linearization, factoring, Padt spatial differencing, the homogeneous property of the flux vectors (where applicable), and hybrid spatial differencing and is most suitable for nonlinear systems in conservation-law form. ADI algorithm retains the order of accuracy and the steady-state property while reducing the bandwidth of the system of equations.[5] Stability of the equation is

 -stable under CFL :  

The order of Truncation error is

 

The result is smooth with considerable overshoot (that does not grow much with time).

References

edit
  1. ^ a b Richard M. Beam, R.F. Warming (September 1976). "An Implicit Finite-Difference Algorithm for Hyperbolic Systems in Conservation-Law Form". Journal of Computational Physics. 22 (1): 87–110. Bibcode:1976JCoPh..22...87B. doi:10.1016/0021-9991(76)90110-8.
  2. ^ Richard M. Beam; R. F. Warming (April 1978). "An Implicit Factored Scheme for the Compressible Navier–Stokes Equations". AIAA Journal. 16 (4): 393–402. Bibcode:1978AIAAJ..16..393B. doi:10.2514/3.60901.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  3. ^ Richard H. Pletcher (2012). Computational Fluid Mechanics and Heat Transfer, Third Edition. CRC Press. ISBN 978-1591690375.
  4. ^ Chung, T.J. (2010). Computational Fluid Dynamics, 2nd Edition. Cambridge University Press. ISBN 978-0521769693.
  5. ^ Lee, Jon (January 1992). "Simplification of Beam and Warming's implicit scheme for two-dimensional compressible flows". AIAA Journal. 30: 266–268. Bibcode:1992AIAAJ..30..266L. doi:10.2514/3.10908.