PROBLEM DEFINITION
This problem is a flow in a two-dimensional (planar geometry) cavity. The inlet flow profile can either be uniform or parabolic. An expansion is encountered after the entrance. The top wall is heated, the bottom cooled and the left side wall is insulated. The component of gravity acts to generate buoyancy. The geometry is illustrated below:
EXECUTABLE AND INPUT FILE
Input file (PROB.INP)
Windows 3.x/9x/ME/2000/NT Executable (PCBUOY.EXE)
Both files compressed in a ZIP archive (PCBUOY.ZIP)
SAMPLE RESULTS
Simulations
Grid used in simulations (generated using ARROW)
Velocity vector field (generated using ARROW)
Streamfunction field contours (generated using ARROW)
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
]======================================
*
===[
PROBLEM PARAMETRES
]===============================================
Length
of duct
[ALTOT]
0.192
Height
of small duct
[HSMALL]
0.0048
Height
of large duct
[HLARGE]
0.0096
Y
Location of abrupt expansion [JSTEP]
21
Axial
velocity at inlet of flow domain [UIN]
0.1
Parabolic inlet velocity profile
[UPARAB]
.TRUE.
Temperature of bottom wall
(K)
[TIN]
293
Temperature of top wall (K)
[TWALL]
345
Gravity m/secı
[GRAV]
9.81
===[
GRID PARAMETERS
]==================================================
Cylindrical coordinate system ?
[POLAR]
.FALSE.
Number
of cells in X-direction [NI]
42
Number
of cells in Y-direction [NJ]
42
Grid
expansion factor in X-direction [FEXPX]
1.0
Grid
expansion factor in Y-direction [FEXPY]
1.0
===[
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]
.TRUE.
===[
THERMO-PHYSICAL PROPERTIES
]=======================================
Material Density
[DENSIT]
1.16
Laminar fluid viscosity
[VISCOS]
1.8E-5
Laminar
Prandtl number
[PRANDT]
0.7
Update
thermo-physical properties ? [INPRO]
.FALSE.
===[
TERMINATION CRITERIA FOR ITERATIONS
]==============================
Maximum number of iterations per
time step
[MAXIT]
500
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.1
Under-relaxation factor for
T-equation
[URFT]
0.5
===[
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]
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]
15
J-index of monitoring location
[JMON]
15
Iteration Index for field printout
[NITPRI]
50
Time
Index for field printout [NSTPRI]
1
Read
IN initial field
[INREAD]
.TRUE.
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]
.TRUE.