Premium
A FULLY COUPLED, IMPLICIT, NUMERICAL SCHEME FOR LAMINAR AND TURBULENT PARABOLIC FLOWS
Author(s) -
VILLASENOR RAFAEL
Publication year - 1997
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/(sici)1097-0207(19970530)40:10<1821::aid-nme140>3.0.co;2-f
Subject(s) - laminar flow , turbulence , mathematics , discretization , numerical methods for ordinary differential equations , partial differential equation , differential equation , mechanics , ordinary differential equation , mathematical analysis , physics , differential algebraic equation
A fully coupled, implicit, numerical scheme has been developed for solving highly stiff systems of parabolic conservation equations. The finite‐domain equations are formed by integration of the governing conservation equations, expressed in vector notation, over control volumes. The central idea is to put the discretized vector equations as Differential/Algebraic equations (DAE) in the context of numerically solving a system of stiff ordinary differential equations . One of the benefits of the present numerical method is that the Jacobian matrix retains its banded property and thus the problems related to computer storage are eliminated. The mathematical interface that the code employs is a computer program that is most efficient for solving the stiff equations usually found in chemical kinetics. The method has shown good convergence rate and numerical stability. The mathematical formulation is capable of modelling laminar and turbulent flows. Two study cases are considered in this work to illustrate the applicability of the numerical method. Case one corresponds to a numerical simulation of a laminar non‐premixed methane–air flame. The chemical processes are described with a four‐step reduced mechanism that derives from a larger mechanism with twenty‐five reaction steps. The predicted velocity, temperature, and mole fractions of the confined laminar round‐jet flame investigated here are compared with experimental data taken from the literature. The versatility of the present solver for predicting high Reynolds number flows is demonstrated by simulating an isothermal turbulent air jet, labelled as case two. A standard two‐equation turbulence model is used to simulate the turbulence processes. The analytical velocity distribution for a turbulent submerged jet is used as a benchmark to test the performance of the model for turbulent jets. © 1997 by John Wiley & Sons, Ltd.