PROBLEM DEFINITION
This problem is a flow in a two-dimensional cavity. The inlet flow profile can either be uniform or parabolic. After length AL1, the cavity suddenly expands and the flow develops a recirculation until the boundary layer re-attaches and develops a parabolic profile. The geometry is illustrated below:
EXECUTABLE AND INPUT FILE
Input file (PROB.INP)
Windows 3.x/9x/ME/2000/NT Executable (PCSTEP.EXE)
Both files compressed in a ZIP archive (PCSTEP.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)
Sample results input file
(PROB.INP)
*
*
===[ PARAMETER'S FILE FOR PROB
]======================================
*
===[
PROBLEM PARAMETRES
]===============================================
Length
of small tube
[AL1]
0.0245
Length
of large tube
[AL2]
0.06
Height
of small tube
[HSMALL]
0.0052
Height
of large tube
[HLARGE]
0.0101
X
Location of abrupt expansion [ISTEP]
10
Y
Location of abrupt expansion [JSTEP]
11
Axial
velocity at inlet of flow domain [UIN]
0.1533488
0.2300232
Parabolic or uniform inlet
velocity profile
[UPARAB]
.f.
===[
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.1
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]
.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]
1250
Maximum acceptable value of residuals
[SORMAX]
5e-9
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]
30
===[
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]
10
J-index of monitoring location
[JMON]
10
Iteration Index for field printout
[NITPRI]
50
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.