<p>Linear stationary reaction-convection-diffusion equations with Dirichlet boundary conditions are approximated using a simple finite difference method corresponding to central differences and the addition of a high-order stabilization term called a numerical moment. The focus is on convection-dominated equations, and the formulation for the method is motivated by various results for fully nonlinear problems. The method features higher-order local truncation errors than monotone methods consistent with the use of the central difference approximation for the gradient. Stability and rates of convergence are derived in the <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\ell ^2\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>ℓ</mi> <mn>2</mn> </msup> </math></EquationSource> </InlineEquation> norm for the constant-coefficient case. Numerical tests are provided to compare the new methods to monotone methods. The methods are also tested for stationary Hamilton-Jacobi equations where they demonstrate higher rates of convergence than the Lax-Friedrich’s method when the underlying viscosity solution is smooth and comparable performance when the underlying viscosity solution is not smooth.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Numerical Moment Stabilization of Central Difference Approximations for Linear Stationary Reaction-Convection-Diffusion Equations

  • Thomas Lewis,
  • Xiaohuan Xue

摘要

Linear stationary reaction-convection-diffusion equations with Dirichlet boundary conditions are approximated using a simple finite difference method corresponding to central differences and the addition of a high-order stabilization term called a numerical moment. The focus is on convection-dominated equations, and the formulation for the method is motivated by various results for fully nonlinear problems. The method features higher-order local truncation errors than monotone methods consistent with the use of the central difference approximation for the gradient. Stability and rates of convergence are derived in the \(\ell ^2\) 2 norm for the constant-coefficient case. Numerical tests are provided to compare the new methods to monotone methods. The methods are also tested for stationary Hamilton-Jacobi equations where they demonstrate higher rates of convergence than the Lax-Friedrich’s method when the underlying viscosity solution is smooth and comparable performance when the underlying viscosity solution is not smooth.