PROBLEM DEFINITION
This problem is a flow in a two-dimensional (planar
geometry) cavity. The cavity has 3 static walls, and the top is moving at a
EXECUTABLE AND INPUT FILE
Input file (PROB.INP)
Windows 3.x/9x/ME/2000/NT Executable (PCCAV1.EXE)
Both files compressed in a ZIP package (PCCAV1.ZIP)
SAMPLE RESULTS
Simulations
Computational Grid (generated using CONTOUR):
Velocity vector field (generated using ARROW):
Streamlines (generated using ARROW):
Velocity vector field – zoom at bottom right corner (generated using ARROW):
Streamlines – zoom at bottom right corner (generated using ARROW):
Pressure field contours (generated using CONTOUR):
U-velocity field contours (generated using CONTOUR):
V-Velocity field contours (generated using CONTOUR):
Sample results input file
(PROB.INP)
*
* ===[
PARAMETER'S FILE FOR PROB
]======================================
*
===[
PROBLEM PARAMETRES
]===============================================
Width
of cavity
[WIDTH]
1.0
Height
of cavity
[HEIGHT]
1.0
Velocity at lid top [UTOP]
5.0e-2
===[
GRID PARAMETERS
]==================================================
Cylindrical coordinate system
? [POLAR]
.FALSE.
Number
of cells in X-direction
[NI]
48
Number
of cells in Y-direction
[NJ]
48
Grid
expansion factor in X-direction
[FEXPX]
1.05
Grid
expansion factor in Y-direction
[FEXPY]
1.05
===[
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.16
Laminar fluid viscosity [VISCOS]
0.185E-4
Laminar
Prandtl number
[PRANDT]
1.0
Update
thermo-physical properties ?
[INPRO]
.FALSE.
===[
TERMINATION CRITERIA FOR ITERATIONS
]==============================
Maximum number of iterations per
time step
[MAXIT]
250
Maximum acceptable value of
residuals
[SORMAX]
5e-3
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.7
Under-relaxation factor for
V-equation
[URFV]
0.7
Under-relaxation factor for
PP-equation
[URFP]
0.3
Under-relaxation factor for
T-equation
[URFT]
1.0
===[
COORDINATES OF FIXED PRESSURE POINT
]=============================
I-index of location where
pressure is fixed
[IPREF]
2
J-index of location where
pressure is fixed
[JPREF]
2
===[
EQUATION SOLVER SELECTION
]========================================
Solver
selection for U-equation
[ISLVU]
1
Solver
selection for V-equation
[ISLVV]
1
Solver
selection for PP-equation
[ISLVP]
1
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]
3
Number
of application of line iteration for V-equation [LITERV]
3
Number
of application of line iteration for PP-equation [LITERP]
5
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]
25
J-index of monitoring
location
[JMON]
25
Iteration Index for field
printout
[NITPRI]
1000
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.