Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Beam and Warming scheme
open-in-new
<h2 id="introduction">Introduction</h2> <p>This scheme is a spatially factored, non iterative, <a href="/facts/Alternating_direction_implicit_method/uV9e0ccz">ADI</a> scheme and uses <a href="/facts/Backward_Euler_method/lccBRRYY">implicit Euler</a> to perform the time Integration. The algorithm is in delta-form, linearized through implementation of a <a href="/facts/Taylor_series/M0aTuftH">Taylor-series</a>. 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.<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a> </p><p>The delta form of the equation produced has the advantageous property of stability (if existing) independent of the size of the time step.<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a> </p> <h2 id="the-method">The method</h2> <p>Consider the inviscid <a href="/facts/Burgers%27_equation/oslyMzES">Burgers' equation</a> in one dimension </p> ∂ u ∂ t = − u ∂ u ∂ x with x ∈ R {\displaystyle {\frac {\partial u}{\partial t}}=-u{\frac {\partial u}{\partial x}}\quad {\text{with }}x\in R} <p>Burgers' equation in conservation form, </p> ∂ u ∂ t = − ∂ E ∂ x {\displaystyle {\frac {\partial u}{\partial t}}=-{\frac {\partial E}{\partial x}}} <p>where : E = u 2 2 {\displaystyle E={\frac {u^{2}}{2}}} </p> <h3>Taylor series expansion</h3> <p>The expansion of : u i n + 1 {\displaystyle u_{i}^{n+1}} </p> u i n + 1 = u i n + 1 2 [ ∂ u ∂ t | i n + ∂ u ∂ t | i n + 1 ] Δ t + O ( Δ t 3 ) {\displaystyle u_{i}^{n+1}=u_{i}^{n}+{\frac {1}{2}}\left[\left.{\frac {\partial u}{\partial t}}\right|_{i}^{n}+\left.{\frac {\partial u}{\partial t}}\right|_{i}^{n+1}\right]\,\Delta t+O(\Delta t^{3})} <p>This is also known as the <a href="/facts/Trapezoidal_rule/KEg8cLGa">trapezoidal formula</a>. </p> ∴ u i n + 1 − u i n Δ t = − 1 2 ( ∂ E ∂ x | i n + 1 + ∂ E ∂ x | i n + ∂ ∂ x [ A ( u i n + 1 − u i n ) ] ) {\displaystyle \therefore {\frac {u_{i}^{n+1}-u_{i}^{n}}{\Delta t}}=-{\frac {1}{2}}\left(\left.{\frac {\partial E}{\partial x}}\right|_{i}^{n+1}+\left.{\frac {\partial E}{\partial x}}\right|_{i}^{n}+{\frac {\partial }{\partial x}}\left[A(u_{i}^{n+1}-u_{i}^{n})\right]\right)} ∵ ∂ u ∂ t = − ∂ E ∂ x {\displaystyle \because {\frac {\partial u}{\partial t}}=-{\frac {\partial E}{\partial x}}} <p>Note that for this equation, </p> A = ∂ E ∂ u = u {\displaystyle A={\frac {\partial E}{\partial u}}=u} <h3>Tri-diagonal system</h3> <p>Resulting tri-diagonal system: </p> − Δ t 4 Δ x ( A i − 1 n u i − 1 n + 1 ) + u i n + 1 + Δ t 4 Δ x ( A i + 1 n u i + 1 n + 1 ) = u i n − 1 2 Δ t Δ x ( E i + 1 n − E i − 1 n ) + Δ t 4 Δ x ( A i + 1 n u i + 1 n − A i − 1 n u i − 1 n ) {\displaystyle {\begin{aligned}&-{\frac {\Delta t}{4\,\Delta x}}\left(A_{i-1}^{n}u_{i-1}^{n+1}\right)+u_{i}^{n+1}+{\frac {\Delta t}{4\,\Delta x}}\left(A_{i+1}^{n}u_{i+1}^{n+1}\right)\\[6pt]={}&u_{i}^{n}-{\frac {1}{2}}{\frac {\Delta t}{\Delta x}}\left(E_{i+1}^{n}-E_{i-1}^{n}\right)+{\frac {\Delta t}{4\,\Delta x}}\left(A_{i+1}^{n}u_{i+1}^{n}-A_{i-1}^{n}u_{i-1}^{n}\right)\end{aligned}}} <p>This resulted system of linear equations can be solved using the modified <a href="/facts/Tridiagonal_matrix_algorithm/NREF5x1c">tridiagonal matrix algorithm</a>, also known as the Thomas algorithm.<a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a> </p> <h2 id="dissipation-term">Dissipation term</h2> <p>Under the condition of shock wave, dissipation term is required for <a href="/facts/Hyperbolic_partial_differential_equation/znXgGuvh">nonlinear hyperbolic equations</a> such as this. This is done to keep the solution under control and maintain convergence of the solution. </p> D = − ε e ( u i + 2 n − 4 u i + 1 n + 6 u i n − 4 u i − 1 n + u i − 2 n ) {\displaystyle D=-\varepsilon _{e}(u_{i+2}^{n}-4u_{i+1}^{n}+6u_{i}^{n}-4u_{i-1}^{n}+u_{i-2}^{n})} <p>This term is added explicitly at level n {\displaystyle n} to the right hand side. This is always used for successful computation where high-frequency oscillations are observed and must be suppressed. </p> <h2 id="smoothing-term">Smoothing term</h2> <p>If only the stable solution is required, then in the equation to the right hand side a second-order <a href="/facts/Relaxation_(iterative_method)/9tsBUvxs">smoothing term</a> 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 </p> ∇ n ( U ) = 0 {\displaystyle \nabla ^{n}(U)=0} <p>The addition of smoothing term increases the number of steps required by three. </p> <h2 id="properties">Properties</h2> <p>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.<a class="footnote-ref" id="fnref:6" href="#fn:6"><sup>6</sup></a> Stability of the equation is </p> L 2 {\displaystyle L^{2}} -stable under CFL : | a | Δ t ≤ 2 Δ x {\displaystyle |a|\,\Delta t\leq 2\,\Delta x} <p>The order of Truncation error is </p> O ( ( Δ t ) 2 + ( Δ x ) 2 ) {\displaystyle O((\Delta t)^{2}+(\Delta x)^{2})} <p>The result is smooth with considerable overshoot (that does not grow much with time). </p> <h2 id="references">References</h2> <ol> <li id="fn:1"><p>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. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li> <li id="fn:2"><p>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) <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li> <li id="fn:3"><p>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. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li> <li id="fn:4"><p>Richard H. Pletcher (2012). Computational Fluid Mechanics and Heat Transfer, Third Edition. CRC Press. ISBN 978-1591690375. <a href="978-1591690375" target="_blank">978-1591690375</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li> <li id="fn:5"><p>Chung, T.J. (2010). Computational Fluid Dynamics, 2nd Edition. Cambridge University Press. ISBN 978-0521769693. <a href="978-0521769693" target="_blank">978-0521769693</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li> <li id="fn:6"><p>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. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li> </ol>