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

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 ]===================================================

.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

.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