PROBLEM DEFINITION

This problem is a flow in a two-dimensional (planar geometry) cavity. The inlet flow profile can either be uniform or parabolic. An expansion is encountered after the entrance. The top wall is heated, the bottom cooled and the left side wall is insulated. The component of gravity acts to generate buoyancy. The geometry is illustrated below:

EXECUTABLE AND INPUT FILE

Input file (PROB.INP)

Windows 3.x/9x/ME/2000/NT Executable (PCBUOY.EXE)

Both files compressed in a ZIP archive (PCBUOY.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)

Temperature field contours (generated using CONTOUR)

Sample results input file (PROB.INP)

*

* ===[ PARAMETER'S FILE FOR PROB ]======================================

*

===[ PROBLEM PARAMETRES ]===============================================

Length of duct  [ALTOT]

0.192

Height of small duct  [HSMALL]

0.0048

Height of large duct  [HLARGE]

0.0096

Y Location of abrupt expansion  [JSTEP]

21

Axial velocity at inlet of flow domain  [UIN]

0.1

Parabolic inlet velocity profile  [UPARAB]

.TRUE.

Temperature of bottom wall (K)  [TIN]

293

Temperature of top wall (K)  [TWALL]

345

Gravity m/secı  [GRAV]

9.81

===[ 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.0

Grid expansion factor in Y-direction  [FEXPY]

1.0

===[ TIME PARAMETERS ]==================================================

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

.TRUE.

===[ THERMO-PHYSICAL PROPERTIES ]=======================================

Material Density  [DENSIT]

1.16

Laminar fluid viscosity  [VISCOS]

1.8E-5

Laminar Prandtl number  [PRANDT]

0.7

Update thermo-physical properties ?  [INPRO]

.FALSE.

===[ TERMINATION CRITERIA FOR ITERATIONS ]==============================

Maximum number of iterations per time step  [MAXIT]

500

Maximum acceptable value of residuals  [SORMAX]

5e-3

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.1

Under-relaxation factor for T-equation  [URFT]

0.5

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

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]

2

Number of application of line iteration for V-equation  [LITERV]

2

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]

15

J-index of monitoring location  [JMON]

15

Iteration Index for field printout  [NITPRI]

50

Time Index for field printout  [NSTPRI]

1

.TRUE.

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]

.TRUE.