A stabilized numerical method with high accuracy is proposed to solve the stability problems
which are generated by spectral element method for two-dimensional transient reaction-convection-diffusion equation. The method combines Chebyshev spectral element method with consistent approximate upwind method in space and adopts fractional-step θ-pattern in time. The numerical example is solved and compared with analytical solution to verify the accuracy and numerical stability. Reaction-convection-diffusion problems with different kinds of boundary layers are also solved. This approach reveals that the stability domain of the spectral element method for reaction-convection-diffusion equation is enlarged with the addition of the consistent approximate upwind term
and the high accuracy of the numerical solution is maintained when convection and reaction term are dominated. For the complicated reaction-convection-diffusion problems with boundary layers
the numerical solution is able to obtain uniform convergence in the whole computational domain.
关键词
Keywords
references
ZIENKIEWICZ O C, HEINRICH J C. The finite element method and convection problems in fluid mechanics [M]. New York, USA: John Wiley & Sons Inc., 1978: 1-22.
LIMA R C, MESQUITA A L A, BLANCO C J C, et al. On the application of SUPG/h-method in 2D advection-diffusion-reaction simulation [J]. J Braz Soc Mech Sci Eng, 2014, 36: 591-603
BOCHEV P, GUNZBURGER M, SHADID J. Stability of the SUPG finite element method for transient advection-diffusion problems [J]. Comput Methods Appl Mech Engrg, 2004, 193: 2301-2323.
LIN R C. A discontinuous Galerkin least-squares finite element method for solving coupled singularly perturbed reaction-diffusion equations [J]. J Comput Appl Math, 2016, 307: 134-142.
DOCARMO E G D, ALVAREZ G B. A new stabilized finite element formulation for scalar convection-diffusion problems: the streamline and approximate upwind/Petrov-Galerkin method [J]. Comput Methods Appl Mech Engrg, 2003, 192: 3379-3396.
GALE(~overA)O A C, ALMEIDA R C, MALTA S M C, et al. Finite element analysis of convection dominated reaction-diffusion problems [J]. Appl Numer Math, 2004, 48: 205-222.
PETERA A T. A spectral element method for fluid dynamics: laminar flow in a channel expansion [J]. J Comput Phys, 1984, 54(3): 468-488.
QIN Guoliang, XU Zhong. A spectral element method for incompressible Navier-Stokes equations [J]Chin J Appl Mech, 2000, 17(4): 20-26.
LUBE G, RAPIN G. Residual-based stabilized higher-order FEM for advection dominated problems [J]. Comput Methods Appl Mech Engrg, 2006, 195: 4124-4138.
JOHN V, KNOBLOCH P. On spurious oscillations at layers diminishing(SOLD)methods for convection-diffusion equations: part I A review [J]. Comput Methods Appl Mech Engrg, 2007, 196: 2197-2215.