PROBLEM DEFINITION

This problem examines the flow of a two-phase mixture of Newtonian fluids within a closed rectangular inclined cavity. The two phases are initially placed in an unstable configuration, where the denser phase lies above a lighter phase. Under the influence of gravity, the two phases overturn. The geometry and initial/boundary conditions are illustrated below: EXECUTABLE AND INPUT FILE

Input file (PROB.INP)

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

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

SAMPLE RESULTS

Movies for dense phase

Volume fraction contours, where dense phase is dark and light phase is pale (movie).

Pictures

Volume fraction contours for the dense phase at t=1s (generated using CONTOUR): Velocity vector field for the dense phase at t=1s (generated using ARROW): Streamlines for the light phase at t=1s (generated using ARROW): Sample results input file (PROB.INP)

*

* ===[ PARAMETER'S FILE FOR TFLOW ; TRANSIENT 2D SEDIMENTATION-PART2 ]==

*

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

Width of cavity  [WIDTH]

1.0

Height of cavity  [HEIGHT]

1.0

Angle of flow field against horizontal axis (degrees) [THETA]

-5.0

Volume fraction of FIRST phase

0.5

Gravity Constant [GRAV]

10.0

Interphase friction coefficient [FCOEF]

10.0

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

Cylindrical coordinate system ?  [POLAR]

.FALSE.

Number of cells in X-direction  [NI]

35

Number of cells in Y-direction  [NJ]

35

Grid expansion factor in X-direction  [FEXPX]

15.0

Grid expansion factor in Y-direction  [FEXPY]

15.0

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

.F.

Number of time steps  [NT]

50

Time step size  [DT]

0.1

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

Solve for U-equations ?  [INCALU]

.T.

Solve for V-equations ?  [INCALV]

.T.

Solve for R-equations ?  [INCALR]

.T.

Solve for PP-equation ?  [INCALP]

.T.

Solve for T-equations ?  [INCALT]

.F.

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

Material Density for both phases  [DENS1,DENS2]

1,1000

Laminar fluid viscosity for both phases  [VISC1,VISC2]

1E-5, 1E-3

Laminar Prandtl number  [PRAND1,PRAND2]

1.0,1.0

Update thermo-physical properties ?  [INPRO]

.FALSE.

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

Maximum number of iterations per time step  [MAXIT]

150

Maximum acceptable value of residuals  [SORMAX]

5E-3

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

Under-relaxation factor for U-equation  [URFU]

0.3

Under-relaxation factor for V-equation  [URFV]

0.3

Under-relaxation factor for R-equation  [URFR]

0.5

Under-relaxation factor for PP-equation  [URFP]

0.1

Under-relaxation factor for T-equation  [URFT]

1.0

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

2

Solver selection for V-equations [ISLVV]

2

Solver selection for PP-equation [ISLVP]

2

Solver selection for R-equations [ISLVP]

2

Solver selection for T-equations [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 R-equation  [LITERR]

2

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

2

===[ INPUT-OUTPUT PARAMETERS ]==========================================

I-index of monitoring location  [IMON]

6

J-index of monitoring location  [JMON]

6

Iteration Index for residual printout  [NITRES]

5

Iteration Index for field printout  [NITPRI]

1000

Time Index for field printout  [NSTPRI]

10

.F.

Name of restart file  [RSTFIL]

RESTART.

Copy OUT FINAL field  [INWRIT]

.T.

Name of SAVE file  [SAVFIL]

RESTART.

Name of main output file  [OUTFIL]

RUN.OUT

Generate velocity field file ?  [INPLOT]

.T.

Generate pressure field file ?  [INPRES]

.T.

Generate volume fraction field file ?  [INVFRC]

.T.

Generate temperature field file ?  [INTEMP]

.F.

===[ INTERPHASE TRACKING PARAMETERS ]===================================

Create interphase tracking reports ? [INTRAC]

.T.

Limit for zero volume fraction (or interphase tracking) [ZEROLIM]

0.01