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 well mixed, and then separate under the influence of gravity. Inside the cavity is an internal obstruction that captures some of the dense phase that settles from above. We are examining non-steady state flow. 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 (PCTrgh.EXE)

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

SAMPLE RESULTS

Movies for dense phase

Volume fraction contours, for a 10 degree inclination (movie).

Volume fraction contours, for a 45 degree inclination (movie).

Pictures

Volume fraction contours for the dense phase (generated using CONTOUR): Velocity vector field for the dense phase (generated using ARROW): U-velocity field contours for the dense phase (generated using CONTOUR): V-velocity field contours for the dense phase (generated using CONTOUR): Sample results input file (PROB.INP)

*

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

*

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

Width of cavity  [WIDTH]

0.5

Height of cavity  [HEIGHT]

1.0

X-Position of expansion [ITROUGH]

22

Y-Position of expansion [JTROUGH]

14

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

-45.0

Volume fraction of FIRST phase (WATER) [R1INIT]

0.5

Gravity Constant [GRAV]

9.81

Interphase friction coefficient [FCOEF]

10.

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

5.0

Grid expansion factor in Y-direction  [FEXPY]

5.0

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

.F.

Number of time steps  [NT]

100

Time step size  [DT]

1.0

===[ 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  [RHO1,RHO2]

1.0, 0.001

Laminar fluid viscosity for both phases  [VISC1,VISC2]

1e-3, 0.0

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]

25

Maximum acceptable value of residuals  [SORMAX]

5.E-03

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

Under-relaxation factor for U-equation  [URFU]

0.2 0.5

Under-relaxation factor for V-equation  [URFV]

0.2 0.5

Under-relaxation factor for R-equation  [URFR]

0.7 0.1

Under-relaxation factor for PP-equation  [URFP]

0.1 0.5

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]

4

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

Solver selection for U-equations [ISLVU]

1

Solver selection for V-equations [ISLVV]

1

Solver selection for PP-equation [ISLVP]

1

Solver selection for R-equations [ISLVP]

1

Solver selection for T-equations [ISLVT]

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

Number of application of line iteration for R-equation  [LITERR]

2

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

1

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

I-index of monitoring location  [IMON]

2

J-index of monitoring location  [JMON]

6

Iteration Index for residual printout  [NITRES]

1

Iteration Index for field printout  [NITPRI]

1000

Time Index for field printout  [NSTPRI]

10

.t.

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.