PROBLEM DEFINITION
This problem examines the laminar 2D flow of a Newtonian
fluid around a solid cube placed within a free stream. We are examining the
flow at steady state. The geometry is illustrated below:
EXECUTABLE AND INPUT FILE
Input file (PROB.INP)
Windows 3.x/9x/ME/2000/NT Executable (PCBOXF.EXE)
Obtain both files in a compressed package (PCBOXF.ZIP)
SAMPLE RESULTS
Simulations
Computational Grid (generated using CONTOUR):
Streamlines – focus on region close to body (generated using
ARROW):
U-Velocity field contours – focus on region close to body
(generated using CONTOUR):
V-Velocity field contours – focus on region close to body
(generated using CONTOUR):
Pressure field contours – focus on region close to body
(generated using CONTOUR):
Sample results input file (PROB.INP)
*
* ===[
PARAMETER'S FILE FOR PROB ]======================================
*
===[
PROBLEM PARAMETRES ]===============================================
Start pos
of box [XL1,YL1]
0.5,4.5
Dimension
of field [XL2,YL2]
10.0,10.0
Box
dimension [BXD]
1.0
Number of
cells in box (X- and Y- direction)
[NSTX,NSTY]
5,5
Fluid
velocity and angle againt horizontal axis (degrees) [UIN,UANGLE]
0.5,0.0
===[ GRID
PARAMETERS ]==================================================
Number of
cells in X-direction [NI]
48
Number of
cells in Y-direction [NJ]
48
Grid expansion
factor in X-direction BEFORE and AFTER box [FEXPX1,FEXPX2]
0.7 1.075
Grid
expansion factor in Y-direction BEFORE and AFTER box [FEXPY1,FEXPY2]
0.75 1.15
===[ TIME
PARAMETERS ]==================================================
Steady
state flow ? [STEADY]
.TRUE.
Number of
time steps [NT]
1
Time step
size [DT]
1.0
===[
DEPENDENT VARIABLES ]==============================================
Solve for
U-equation ? [INCALU]
.TRUE.
Solve for
V-equation ? [INCALV]
.TRUE.
Solve for
PP-equation ? [INCALP]
.TRUE.
Solve for
T-equation ? [INCALT]
.FALSE.
===[
THERMO-PHYSICAL PROPERTIES ]=======================================
Material
Density [DENSIT]
1.0
Laminar
fluid viscosity [VISCOS]
1.85E-3
Laminar Prandtl number
[PRANDT]
1.0
Update thermo-physical
properties ? [INPRO]
.FALSE.
===[
TERMINATION CRITERIA FOR ITERATIONS ]==============================
Maximum
number of iterations per time step
[MAXIT]
2500
Maximum
acceptable value of residuals [SORMAX]
2.5e-4
Normalization
quantity for mass residual [SNORM]
1.0
Normalization
quantity for U-momentum residual [SNORU]
1.0
Normalization
quantity for V-momentum residual [SNORV]
1.0
Normalization
quantity for T-equation residual [SNORT]
1.0
===[
UNDER-RELAXATION ]=================================================
Under-relaxation
factor for U-equation [URFU]
0.4
Under-relaxation
factor for V-equation [URFV]
0.4
Under-relaxation
factor for PP-equation [URFP]
0.1
Under-relaxation
factor for T-equation [URFT]
1.0
===[ EQUATION
SOLVER SELECTION ]========================================
Solver
selection for U-equation [ISLVU]
2
Solver
selection for V-equation [ISLVV]
2
Solver
selection for PP-equation [ISLVP]
2
Solver
selection for T-equation [ISLVT]
1
Constant used
by SIM solver [ALPHA]
0.8
===[
TERMINATION CRITERIA FOR EQUATION SOLVER ]=========================
Number of
application of line iteration for U-equation
[LITERU]
2
Number of
application of line iteration for V-equation
[LITERV]
2
Number of
application of line iteration for PP-equation
[LITERP]
4
Number of
application of line iteration for T-equation
[LITERT]
1
Maximum
allowable value of solver residual sums for U
[SORU]
1.0E-6
Maximum
allowable value of solver residual sums for V
[SORV]
1.0E-6
Maximum
allowable value of solver residual sums for PP
[SORP]
1.0E-6
Maximum
allowable value of solver residual sums for T
[SORT]
1.0E-6
===[
INPUT-OUTPUT PARAMETERS ]==========================================
I-index of
monitoring location [IMON]
2
J-index of
monitoring location [JMON]
25
Iteration
Index for field printout [NITPRI]
500
Time Index
for field printout [NSTPRI]
1
Read IN
initial field [INREAD]
.FALSE.
Name of
restart file [RSTFIL]
PHIDA.
Copy OUT
FINAL field [INWRIT]
.TRUE.
Name of SAVE
file [SAVFIL]
PHIDA.
Name of
main output file [OUTFIL]
RUN.OUT
Generate
velocity field file ? [INPLOT]
.TRUE.
Generate
pressure field file ? [INPRES]
.TRUE.
Generate
temperature field file ? [INTEMP]
.FALSE.