Page 395 - Computational Fluid Dynamics for Engineers
P. 395
388 Appendix B
Velocity in x-direction: R25Mpl5Ep5G200X50outU.txt
Velocity in y-direction: R25Mpl5Ep5G200X50outV.txt
t
Density: R25Mpl5Ep5G200X50outRho. x t
Pressure: R25Mpl5Ep5G200X50outP. x t
t
Temperature: R25Mpl5Ep5G200X50outT. x t
t
Energy Et: R25Mpl5Ep5G200X50outEt.txt
2. RL - 50, L = 20 a grid of 200 x 50 Mach = 0.15.
Input file: R50Mpl5Ep5G200X50Inp.txt
3. RL = 100, L = 30 a grid of 200 x 50 Mach = 0.15.
Input file: R100Mpl5G200X50Inp.txt
4. RL = 400, L - 40 a grid of 200 x 50 Mach = 0.15.
Input file: R400Mpl5G200X50Inp.txt
Mach Number Effect
1. Mach = 0.15, RL = 50, L = 20 a grid of 300 x 50.
Input file: R50Mpl5Ep5G200X50Inp.txt
2. Mach = 0.7, RL = 50, L = 20 a grid of 300 x 50.
Input file: R50Mp7Ep5G300X50Inp.txt
3. Mach = 1.0, RL = 50, L = 20 a grid of 300 x 50.
Input file: R50MlEp5G300X50Inp.txt
12.6.3.3 Finite-Volume Method for Compressible
Navier-Stokes Equations
FiniteVolume.f contains the main program which reads the input data and
controls the iterations and computations.
The PC executable file is FiniteVolume.exe
InitialFV.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.
ResidualFV.f contains subroutine righthandside which performs the integral
calculation according to Section 12.5
RKFV.f contains subroutines RKStep, updating and Boundary Update
Subroutine RKStep performs calculations according to the fourth-order Runge-
Kutta method, see Eq. (12.5.7)
Subroutines updating and Boundary Update perform boundary calculations ac-
cording to boundary condition and update the solution.
OutputFV.f contains subroutines output-g and output^ and saves the results
into files.
Test Case for sudden expansion laminar duct flow:
RL = 100, L = 30 a grid of 300 x 50 Mach = 0.15
Input file: R100input.txt