PROBLEM DEFINITION
This problem is a flow in a two-dimensional (planar geometry) cavity. The cavity has static walls, one of which is kept hot, and the other, cold. A buoyancy-driven flow develops (gravity acts in the -y direction). 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 (PCCAV2.EXE)
Obtain both files in a compressed package (PCCAV2.ZIP)
SAMPLE RESULTS
Simulations
Computational Grid (generated using CONTOUR):
Velocity vector field (generated using ARROW):
Velocity vector field – zoom at top 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):
Temperature field contours (generated using CONTOUR):
Sample results input file
(PROB.INP)
*
*
===[ PARAMETER'S FILE FOR PROB
]=====================================
*
===[
CARTESIAN / CYLINDRICAL COORDINATE SYSTEM
]=======================
Define
coordinate system [ POLAR ]
.FALSE.
===[
GRID POINTS
]=====================================================
Number
of cells in x-direction [ NI
]
42
Number
of cells in y-direction [ NJ
]
42
===[
DIMENSIONS OF CAVITY
]============================================
Width
and heigth of cavity [
H ]
1e-1
===[
EXPANSION FACTORS FOR GRID
]======================================
Expansion factor in X-direction
[ FEXPX ]
1.05
Expansion factor in Y-direction
[ FEXPY ]
1.05
===[
COORDINATES OF MONITORING POINT
]=================================
Axial
monitoring location [ IMON ]
20
Radial
or Vertical monitoring location [ JMON
]
20
===[
PROGRAM CONTROL PARAMETERS
]======================================
Maximum number of iterations
[ MAXIT ]
1000
Maximum residual [ SORMAX ]
1E-3
Index
for field printout [ INDPRI ]
1000
===[
DEPENDENT VARIABLES
]=============================================
Solve
for u-equation [ INCALU ]
.TRUE.
Solve
for v-equation [ INCALV ]
.TRUE.
Solve
for p'-equation [ INCALP ]
.TRUE.
Solve
for t-equation [ INCALT ]
.TRUE.
Update
fluid parameters [ INPRO ]
.FALSE.
===[
NUMBER OF SWEEPS FOR ALL VARIABLES
]==============================
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 p'-equation [
LITERP ]
4
Number
of application of line iteration for T-equation [
LITERT ]
2
===[
UNDER-RELAXATION
]================================================
Under
relaxation factor for u-velocity [ URFU
]
0.5
Under
relaxation factor for v-velocity [ URFV
]
0.5
Under
relaxation factor for pressure correction [ URFP
]
0.1
Under
relaxation factor for temperature [ URFT
]
1
===[
COORDINATES OF FIXED PRESSURE POINT
]=============================
Axial
fixed-pressure location index [ IPREF
]
2
Radial
or Vertical fixed-pressure location index [ JPREF
]
2
===[
FLUID PROPERTIES+PROBLEM PARAMETERS
]=============================
Fluid
density at calculation domain [ DENSIT
]
1.207
Fluid
kinematic viscosity [
VISCOS ]
1.505E-05
Prandl number [
PRANDT ]
0.71
Thermal expansion coef. [ BETA
]
3.66E-03
Gravit.
accel. [ GI ]
9.812
Hot
wall temp. [ THOT ]
303
Cold wall
temp. [ TCOLD ]
283
===[
SOLVER SELECTION
]================================================
Solver
for u-velocity equations [ ISLVU
]
1
Solver
for v-velocity equations [ ISLVV
]
1
Solver
for pressure correction equationn [ ISLVP ]
1
Solver
for temperature equations [ ISLVT
]
1
Constand used by S.I.M. Solver [ ALPHA ]
0.8
===[
FILING SYSTEM
]===================================================
Name
of main output file [ OUTFIL
]
CAVITY.OUT
Read
IN flow field [ INREAD ]
.FALSE.
Name
of restart file
[RSTFIL]
PHIDA.
Copy
OUT flow field [ INWRITE ]
.TRUE.
Name
of SAVE file
[SAVFIL]
PHIDA.
Generate PLOT file [ INPLOT ]
.TRUE.
Generate Pressure Contour Lines
[ INPRES ]
.TRUE.
Generate Temperature Contour
lines [ INTEMP ]
.TRUE.