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.