Page 394 - Computational Fluid Dynamics for Engineers

Appendix B 387

Pressure: R100outP.txt
Temperature: R100outT.txt
Energy Et: R100outEt.txt
Reattachment location: R100reAttach.txt
Restart file: R100restart.txt

12.6.3.2 Beam-Warming for Compressible Navier—Stokes Equations

BeamWarming.f contains the main program which reads input data and controls
the iterations and computations. The PC executable file is BeamWarming.exe
Initial.f contains two subroutines initialization and restart
Subroutine initialization sets up initial conditions according to subsection 12.6.2
Subroutine restart sets up restarting calculations.
Residual.f contains subroutine righthandside which calculates the right hand-
side according to Eq. (12.4.19).
solverb4x4.f has subroutines solver, invmatrix^ atimeb, tribandsyssolver,
thomas and gauss
Subroutine solver controls calculation of linear system described by Eqs.(12.4.19)
and (12.4.20)
Subroutine invmatrix4 computes the inverse of a 4 by 4 matrix.
Subroutine atimeb calculates 4 by 4 matrix multiplication. C = AB
Subroutine tribandsyssolver solves block tridiagonal system
Subroutines thomas and gauss solve tridiagonal system described in Chapter 4
ainvrs. f contains subroutine ainvrs which calculates the inverse of a general
matrix
FormMatrix.f contains subroutines x.directm and yAirectm
Subroutines xAireetm and y-directm form a block tridiagonal system according
to Eqs. (12.4.19) and (12.4.20), respectively.
Update. f contains subroutine update which updates the values interior points
and boundary points according to the boundary conditions discussed in subsec-
tion 12.6.1
Test cases for sudden expansion laminar duct flow with various values of
Reynolds number, Mach number, grid size and dissipation parameters.

Test Cases for Sudden Expansion Laminar Duct Flow:
Reynolds Number Effect
1. RL = 25, L = 20 a grid of 200 x 50 Mach = 0.15
Input file: R25Mpl5Ep5G200X50Inp.txt
Output files:
t
General: R25Mpl5Ep5G200X50outG. x t
Mass Integration: R25Mpl5Ep5G200X50outMass.txt