Premium
Oscillation limits for weighted residual methods applied to convective diffusion equations
Author(s) -
Jensen Ole Krogh,
Finlayson Bruce A.
Publication year - 1980
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620151109
Subject(s) - mathematics , galerkin method , hermite interpolation , mathematical analysis , quartic function , hermite polynomials , collocation method , quadratic equation , finite element method , geometry , physics , differential equation , ordinary differential equation , pure mathematics , thermodynamics
Convection–diffusion equations are difficult to solve when the convection term dominates because most solution methods give solutions which oscillate in space. Previous criteria based on the one‐dimensional convection–diffusion equation have shown that finite difference and Galerkin (linear or quadratic basis functions) will not give oscillatory solutions provided the Peclet number times the mesh size (Pe Δ x ) is below a critical value. These criteria are based on the solution at the nodes, and ensure that the nodal values are monotone. Similar criteria are developed here for other methods: quadratic Galerkin with upwind weighting, cubic Galerkin, orthogonal collocation on finite elements with quadratic, cubic or quartic polynomials using Lagrangian interpolation, cubic or quartic polynominals using Hermite interpolation, and the method of moments. The nodal values do not oscillate for collocation or moments methods with Hermite cubic polynomials regardless of the value of Pe Δ x . A new criterion is developed for all methods based on the monotonicity of the solutions throughout the domain. This criterion is more restrictive than one based only on the nodal values. All methods that are second order (Δ x 2 ) or better in truncation error give oscillatory solutions (based on the entire domain) unless Pe Δ x is below a critical value. This value ranges from 2 for finite difference methods to 4·6 for Hermite, quartic, collocation methods.