Page 391 - Computational Fluid Dynamics for Engineers
P. 391
384 Appendix B
t
10.14 Model Problem for he Implicit Method:
Quasi-ID Nozzle (nozzle_imp.c)
Generate executable nozzle.imp on unix by using
cc -o nozzle_imp nozzle_imp.c b l o c k t r i . c -lm
Test Cases with PC executable nozzleImpPC.exe and NMAX = 100
1. CFL = 0.5, e e = 1, Output: CFL0p5results.txt and CFL0p5residuals.txt
2. CFL = 1.0, e e = 1, Output: CFLlresults.txt and CFLlresiduals.txt
3. CFL = 3.0, e e = 1, Output: CFL3results.txt and CFL3residuals.txt
4. CFL = 1.0, s e = 10, Output: CFLlEelOresults.txt and
CFLlEelOresiduals.txt
Chapter 11
Incompressible Navier-Stokes Equations
11.2. INS2D Method for Incompressible Navier—Stokes Equations
ins2d. f contains main program which reads input data and controls the itera-
tions and computations.
Residual.f contains subroutines flux.de, phLih, residualx, flux-df, phLjh and
residualy.
Subroutines fluxAE and fluxAF are for numerical flux computations, see Eq.
(11.2.15).
Subroutines phLih and phLjh are defined according to Eq. (11.2.17).
Subroutines residualx and residualy calculates righthandside according to Eq.
(11.2.22).
solver. f contains subroutines solver, invmatrixS, tribandsyssolver and thomas
which use the block elimination method to solve block-tridiagonal system.
Test cases for sudden expansion laminar duct flow with various values of /?, grid
size and Reynolds number.
Test Cases for Sudden Expansion Laminar Duct Flow:
1. (3 — 1, grid size nx = 200 and ny = 100,
Input: blg200X100Inp.txt
Output: blg200X1000ut_g.txt; blg200X1000ut_Mass.txt;
blg200X1000ut_p.txt; blg200X1000ut_u.txt; blg200X1000ut_v.txt
2. /3 = 50, grid size nx = 200 and ny = 100,
Input: b50g200X100Inp.txt
Output: b50g200X1000ut_g.txt; b50g200X1000ut_Mass.txt;
t
b50g200X1000ut_p.txt; b50g200X1000ut_u. x t ; b50g200X1000ut_v.txt
3. (5 = 100, grid size nx = 200 and ny = 100,
Input: bl00g200X100Inp.txt