Premium
New numerical schemes based on a criterion for constructing essentially stable and accurate numerical schemes for convection‐dominated equations
Author(s) -
Sakai Katsuhiro
Publication year - 1995
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1650211103
Subject(s) - mathematics , scheme (mathematics) , stability (learning theory) , order (exchange) , convection , numerical analysis , polynomial , upwind scheme , numerical diffusion , mathematical analysis , convection–diffusion equation , physics , mechanics , computer science , economics , finance , machine learning , discretization
In order to obtain stable and accurate numerical solutions for the convection‐dominated steady transport equations, we propose a criterion for constructing numerical schemes for the convection term that the roots of the characteristic equation of the resulting difference equation have poles. By imposing this criterion on the difference coefficients of the convection term, we construct two numerical schemes for the convection‐dominated equations. One is based on polynomial differencing and the other on locally exact differencing. The former scheme coincides with the QUICK scheme when the mesh Reynolds number ( Rm ) is \documentclass{article}\pagestyle{empty}\begin{document}$\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $\end{document} , which is the critical value for its stability, while it approaches the second‐order upwind scheme as Rm goes to infinity. Hence the former scheme interpolates a stable scheme between the QUICK scheme at Rm = \documentclass{article}\pagestyle{empty}\begin{document}$\mathop \[{\textstyle{{\rm 8} \over {\rm 3}}}\] $\end{document} and the second‐order upwind scheme at Rm = ∞. Numerical solutions with the present new schemes for the one‐dimensional, linear, steady convection‐diffusion equations showed good results.