PROBLEM DEFINITION
This problem examines the laminar axisymmetric flow of a Newtonian fluid through a sudden contraction (orifice). 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 (PCORFC_L.EXE)
Obtain both files in a compressed package (PCORFC_L.ZIP)
SAMPLE RESULTS
Simulations
Computational Grid (generated using CONTOUR):
Velocity vector field – focus on recirculation downstream of the orifice (generated using ARROW):
Streamlines – focus on recirculation downstream of the orifice (generated using ARROW):
Streamlines for whole flow field (generated using ARROW):
U-Velocity field contours (generated using CONTOUR):
Sample results input file
(PROB.INP)
*
*
===[ PARAMETER'S FILE FOR PROB
]======================================
*
===[
PROBLEM PARAMETRES
]===============================================
Length
of tube
[ LENGTH ]
.25
Radius
of tube
[ RADIUS ]
0.01
Radius
of orifice
[ RORIF ]
0.005
Location of orifice [ XBLCK ]
.015
Axial
velocity at inlet of flow domain [ UIN ]
0.01
Parabolic or uniform inlet
velocity profile
[ UPARAB ]
.FALSE.
===[
GRID PARAMETERS
]==================================================
Cylindrical coordinate system ?
[POLAR]
.TRUE.
Number
of cells in X-direction [NI]
48
Number
of cells in Y-direction [NJ]
48
Grid
expansion factor in X-direction [FEXPX]
6
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]
.T.
Solve
for V-equation ? [INCALV]
.T.
Solve
for PP-equation ? [INCALP]
.T.
===[
THERMO-PHYSICAL PROPERTIES
]=======================================
Material Density
[DENSIT]
995.3
Laminar fluid viscosity
[VISCOS]
0.8E-3
Laminar
Prandtl number
[PRANDT]
5.43
===[
TERMINATION CRITERIA FOR ITERATIONS
]==============================
Maximum number of iterations per
time step
[MAXIT]
5000
Maximum acceptable value of residuals
[SORMAX]
5.0E-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
===[
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
===[
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]
2
Solver
selection for V-equation [ISLVV]
2
Solver
selection for PP-equation [ISLVP]
2
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
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
===[
INPUT-OUTPUT PARAMETERS
]==========================================
I-index of monitoring location
[IMON]
10
J-index of monitoring location
[JMON]
19
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]
.T.
Generate pressure field file ?
[INPRES]
.T.