PROBLEM DEFINITION

This problem is a flow in a two-dimensional (planar geometry) cavity. The cavity has static walls, one of which is kept hot, and the other, cold. A buoyancy-driven flow develops (gravity acts in the -y direction). 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 (PCCAV2.EXE)

Obtain both files in a compressed package (PCCAV2.ZIP)

SAMPLE RESULTS

Simulations

Computational Grid (generated using CONTOUR):

Velocity vector field (generated using ARROW):

Velocity vector field – zoom at top right corner (generated using ARROW):

Pressure field contours (generated using CONTOUR):

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

*

===[ CARTESIAN / CYLINDRICAL COORDINATE SYSTEM ]=======================

Define coordinate system [ POLAR ]

.FALSE.

===[ GRID POINTS ]=====================================================

Number of cells in x-direction [ NI ]

42

Number of cells in y-direction [ NJ ]

42

===[ DIMENSIONS OF CAVITY ]============================================

Width and heigth of cavity [ H ]

1e-1

===[ EXPANSION FACTORS FOR GRID ]======================================

Expansion factor in X-direction [ FEXPX ]

1.05

Expansion factor in Y-direction [ FEXPY ]

1.05

===[ COORDINATES OF MONITORING POINT ]=================================

Axial monitoring location [ IMON ]

20

Radial or Vertical monitoring location [ JMON ]

20

===[ PROGRAM CONTROL PARAMETERS ]======================================

Maximum number of iterations [ MAXIT ]

1000

Maximum residual [ SORMAX ]

1E-3

Index for field printout [ INDPRI ]

1000

===[ DEPENDENT VARIABLES ]=============================================

Solve for u-equation [ INCALU ]

.TRUE.

Solve for v-equation [ INCALV ]

.TRUE.

Solve for p'-equation [ INCALP ]

.TRUE.

Solve for t-equation [ INCALT ]

.TRUE.

Update fluid parameters [ INPRO ]

.FALSE.

===[ NUMBER OF SWEEPS FOR ALL VARIABLES ]==============================

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 p'-equation [ LITERP ]

4

Number of application of line iteration for T-equation [ LITERT ]

2

===[ UNDER-RELAXATION ]================================================

Under relaxation factor for u-velocity [ URFU ]

0.5

Under relaxation factor for v-velocity [ URFV ]

0.5

Under relaxation factor for pressure correction [ URFP ]

0.1

Under relaxation factor for temperature [ URFT ]

1

===[ COORDINATES OF FIXED PRESSURE POINT ]=============================

Axial fixed-pressure location index [ IPREF ]

2

Radial or Vertical fixed-pressure location index [ JPREF ]

2

===[ FLUID PROPERTIES+PROBLEM PARAMETERS ]=============================

Fluid density at calculation domain [ DENSIT ]

1.207

Fluid kinematic viscosity [ VISCOS ]

1.505E-05

Prandl number [ PRANDT ]

0.71

Thermal expansion coef. [ BETA ]

3.66E-03

Gravit. accel. [ GI ]

9.812

Hot wall temp. [ THOT ]

303

Cold wall temp. [ TCOLD ]

283

===[ SOLVER SELECTION ]================================================

Solver for u-velocity equations [ ISLVU ]

1

Solver for v-velocity equations [ ISLVV ]

1

Solver for pressure correction equationn [ ISLVP ]

1

Solver for temperature equations [ ISLVT ]

1

Constand used by S.I.M. Solver [ ALPHA ]

0.8

===[ FILING SYSTEM ]===================================================

Name of main output file [ OUTFIL ]

CAVITY.OUT

.FALSE.

Name of restart file  [RSTFIL]

PHIDA.

Copy OUT flow field [ INWRITE ]

.TRUE.

Name of SAVE file  [SAVFIL]

PHIDA.

Generate PLOT file [ INPLOT ]

.TRUE.

Generate Pressure Contour Lines [ INPRES ]

.TRUE.

Generate Temperature Contour lines [ INTEMP ]

.TRUE.