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
]==================================================
Steady
state flow ?
[STEADY]
.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
Read
IN initial field
[INREAD]
.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.