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