PROBLEM DEFINITION
This problem examines the laminar axisymmetric flow of a Newtonian fluid in the entry region of a pipe. The tube wall is heated in such a way so that either the temperature or the heat flux per unit area is kept constant. 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 (PCHEAT.EXE)
Obtain both files in a compressed package (PCHEAT.ZIP)
SAMPLE RESULTS
Simulations
Computational Grid (generated using CONTOUR):
Velocity vector field (generated using ARROW):
Streamlines (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 LAMINAR FLOW IN A TUBE
OF A PL FLUID]=======
*
===[
PROBLEM PARAMETRES
]===============================================
Length
of large tube
[WIDTH]
0.8
Radius
of large tube
[RADIUS]
0.005
Axial
velocity at inlet of flow domain [UIN]
3.428
===[
GRID AND GEOMETRY VARIABLES
]=======================================
Cylindrical coordinate system ?
[POLAR]
.TRUE.
Number
of cells in X-direction [NI]
40
Number
of cells in Y-direction [NJ]
40
Grid
expansion factor in X-direction [FEXPX]
1.05
Grid
expansion factor in Y-direction [FEXPY]
0.95
===[
TIME VARIABLES
]===================================================
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]
1000.
Consistency Index
[CI]
0.30636
Power
Law Index
[PLI]
0.716
Laminar Prandtl number [PRANDT]
1.0
Update
thermo-physical properties ? [INPRO]
.TRUE.
===[
ITERATION PARAMETERS
]=============================================
Maximum number of iterations per
time step
[MAXIT]
400
Maximum acceptable value of residuals
[SORMAX]
0.0005
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 factor for
U-equation
[URFU]
0.6
Under-relaxation factor for
V-equation
[URFV]
0.6
Under-relaxation factor for
PP-equation
[URFP]
0.1
Under-relaxation factor for
T-equation
[URFT]
1.0
Under-relaxation factor for VISCOSITY(GREEK ETA)-equation [URFH]
0.7
===[
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 PARAMETERS
]======================================
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
Number
of application of line iteration for U-equation [LITERU]
1
Number
of application of line iteration for V-equation [LITERV]
1
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]
3
J-index of monitoring location
[JMON]
38
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]
RESTART.BIN
Copy
OUT FINAL field
[INWRIT]
.TRUE.
Name
of SAVE file
[SAVFIL]
RESTART.BIN
Name
of main output file
[OUTFIL]
RESULT.DAT