PROBLEM DEFINITION

This problem examines the laminar 2D flow of a Newtonian fluid around a solid cube placed within a free stream. 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 (PCBOXF.EXE)

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

 

 

SAMPLE RESULTS

Simulations

 

Computational Grid (generated using CONTOUR):

 

 

 

Streamlines focus on region close to body (generated using ARROW):

 

 

 

U-Velocity field contours focus on region close to body (generated using CONTOUR):

 

 

 

V-Velocity field contours focus on region close to body (generated using CONTOUR):

 

 

 

Pressure field contours focus on region close to body (generated using CONTOUR):

 

 

 

 

 

Sample results input file (PROB.INP)

*

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

*

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

Start pos of box [XL1,YL1]

0.5,4.5

Dimension of field [XL2,YL2]

10.0,10.0

Box dimension [BXD]

1.0

Number of cells in box (X- and Y- direction) [NSTX,NSTY]

5,5

Fluid velocity and angle againt horizontal axis (degrees) [UIN,UANGLE]

0.5,0.0

===[ GRID PARAMETERS ]==================================================

Number of cells in X-direction [NI]

48

Number of cells in Y-direction [NJ]

48

Grid expansion factor in X-direction BEFORE and AFTER box [FEXPX1,FEXPX2]

0.7 1.075

Grid expansion factor in Y-direction BEFORE and AFTER box [FEXPY1,FEXPY2]

0.75 1.15

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

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

1.0

Laminar fluid viscosity [VISCOS]

1.85E-3

Laminar Prandtl number [PRANDT]

1.0

Update thermo-physical properties ? [INPRO]

.FALSE.

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

Maximum number of iterations per time step [MAXIT]

2500

Maximum acceptable value of residuals [SORMAX]

2.5e-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

Normalization quantity for T-equation residual [SNORT]

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

Under-relaxation factor for T-equation [URFT]

1.0

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

Solver selection for U-equation [ISLVU]

2

Solver selection for V-equation [ISLVV]

2

Solver selection for PP-equation [ISLVP]

2

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]

2

J-index of monitoring location [JMON]

25

Iteration Index for field printout [NITPRI]

500

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]

.TRUE.

Generate pressure field file ? [INPRES]

.TRUE.

Generate temperature field file ? [INTEMP]

.FALSE.