High-Order Finite Element Methods for Singularly Perturbed Elliptic and Parabolic Problems
High-Order Finite Element Methods for Singularly Perturbed Elliptic and Parabolic Problems
复制标题
DOI:
10.1137/s0036139993269345
复制
发表时间:
1995-04
期刊:
影响因子:
--
通讯作者:
S. Adjerid;M. Aiffa;J. Flaherty
中科院分区:
文献类型:
--
作者:
S. Adjerid;M. Aiffa;J. Flaherty
We develop a framework for applying high-order finite element methods to singularly-perturbed elliptic and parabolic differential systems that utilizes special quadrature rules to confine spurious effects, such as excess diffusion and nonphysical oscillations, to boundary and interior layers. This approach is more suited for use with adaptive mesh-refinement and order-variation techniques than other problem-dependent methods. Quadrature rules, developed for two-point convection-diffusion and reaction-diffusion problems, are used with finite element software to solve examples involving ordinary and partial differential equations. Numerical artifacts are confined to layers for all combinations of meshes, orders and singular perturbation parameters that were tested. Radau or Lobatto quadrature used with the finite element method to solve, respectively, convection-and reaction-diffusion problems provide the benefits of the specialized quadrature formulas and are simpler to implement.