PHOTON USE
p
0.20443E+04 0.15633E+04 CR
use patgeo
vec x 1
@
0.73697E+03 0.10031E+04 CR
Velocity vectors for the continuous phase.
Press RETURN to continue.@
pause
red
use genuse
@
0.73697E+03 0.10031E+04
Particle trajectories.
Press e to exit.@
ENDUSE
GROUP 1. Run title
TEXT(Particles through a ball valve: G205
TITLE
DISPLAY
The problem specification for the contiuous phase is the same
as in PHOENICS Library Case 534.
The incompressible, single-phase flow of water through a
fully-open axi-symmetric ball valve is solved. The pipe-work
considered is 2 m in length; the radius at inlet is 0.15 m; the
diminishes to a minimum of 0.09 m in front of the ball; and the
maximum radius is 0.16 m (the command VIEW(I,1) displays the
geometry and grid).
A plug profile of axial velocity is prescribed at the inlet,
and a fixed-pressure condition is employed at outlet. Wall
friction is activated along the boundary of the ball, and along
the pipe wall. A fixed turbulent kinematic viscosity is set to
100 times the laminar value (the Reynolds number is of order
1.0E5).
Particles of different sizes and velocities are introduced
at the valve inlet with several initial velocities.
Gravity is acting along the axial (z) direction.
ENDDIS
GROUP 4. Y-direction grid specification
NY=12
GROUP 5. Z-direction grid specification
NZ=30
GROUP 6. Body-fitted coordinates or grid distortion
BFC=T;NONORT=T;READCO(GRID4)
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);SOLUTN(P1,Y,Y,Y,N,N,N)
GROUP 8. Terms (in differential equations) & devices
** Activate the pure upwind formulation of the finite-domain
equations
DIFCUT=0.0
GROUP 9. Properties of the medium (or media)
RHO1=1.0E3;ENUL=1.0E-6;ENUT=100.0*ENUL
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=2.0
**The cell row which passes through the ball is blocked out.
Note that the minus sign in the 6th argument of CONPOR
activates wall friction on the axial velocity.
CONPOR(0.0,CELL,1,1,1,-1,11,27)
GROUP 13. Boundary conditions and special sources
** Inlet
INLET(INLET,LOW,1,1,1,NY,1,1,1,1)
VALUE(INLET,P1,2.0*RHO1);VALUE(INLET,W1,2.0)
** Outlet
PATCH(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,1000.,0);COVAL(OUTLET,W1,ONLYMS,2.0)
** North wall
WALL(WFUN,NORTH,1,NX,NY,NY,1,NZ,1,1)
GROUP 15. Termination of sweeps
LSWEEP=50
GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,1.0);RELAX(W1,FALSDT,1.0);RELAX(P1,LINRLX,0.2)
GROUP 22. Spot-value print-out
IYMON=2;IZMON=20
GROUP 23. Field print-out and plot control
NYPRIN=2;NZPRIN=NZ/6;NPLT=2
PATCH(DOMAIN,CONTUR,1,1,1,NY,1,NZ,1,1);PLOT(DOMAIN,P1,0.,15.0)
PATCH(INNER,PROFIL,1,1,2,2,1,NZ,1,1);PLOT(INNER,W1,0.0,0.0)
PATCH(OUTER,PROFIL,1,1,NY,NY,1,NZ,1,1);PLOT(OUTER,W1,0.0,0.0)
PATCH(FRONT,PROFIL,1,1,2,NY,10,10,1,1);PLOT(FRONT,W1,0.0,0.0)
PLOT(FRONT,V1,0.0,0.0);PLOT(FRONT,P1,0.0,0.0)
PATCH(BACK,PROFIL,1,1,2,NY,28,28,1,1);PLOT(BACK,W1,0.0,0.0)
PLOT(BACK,V1,0.0,0.0);PLOT(BACK,P1,0.0,0.0)
* GENTRA exit PATCH GXOUTLET
PATCH(GXOUTLET,HIGH,1,1,1,12,30,30,1,1)
* GENTRA exit PATCH GXINLET
PATCH(GXINLET,LOW,1,1,1,12,1,1,1,1)
#GENTRA
*-------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*-------------------------------------------------------
* Particle type - 30 - ISOTHERMAL_PARTICLES
GPTYPE=30
* Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=0.0;GGRAZ=9.8;GBUOYA=F
* Stochastic model of turbulence
GSTOCH=F
* Data for isothermal particles
GDRAG=GRND1
*-------------------------------------------------------
* GENTRA GROUP 2: Boundary conditions for particles
*-------------------------------------------------------
* Inlet-data file-name
GINFIL=Q1
* Use this space for inlet-data table if GINFIL=Q1
*__YP__ZP__VP__WP__DIAM_DENSTY_FRATE
0.01 0.0 0 1 0.001 500.0 1.0E-5
0.04 0.0 0 2 0.0001 1000.0 1.0E-5
0.07 0.0 0 3 0.0015 1000.0 1.0E-5
0.10 0.0 0 0.5 0.002 1000.0 1.0E-5
0.13 0.0 0 1 0.001 1500.0 1.0E-5
* Wall treatment, and rest coefficient if appropriate
GWALLC=3;GWREST=0.75
* Porosity threshold
GPOROS=0
*-------------------------------------------------------
* GENTRA GROUP 3: Numerical controls
*-------------------------------------------------------
* 1st GENTRA sweep; frequency of calls
GSWEP1=40;GSWEPF=1
* Maximum Lagrangian time-step
GDTMAX=1
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=-100;GTIMMX=10
*-------------------------------------------------------
* GENTRA GROUP 4: Output controls
*-------------------------------------------------------
* Trajectory-file identifier, history-file identifier
and frequency for output
* Particle to compute residence time
GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1
GRESTI=2
*-------------------------------------------------------
* GENTRA GROUP 5: Provisions for the GENTRA-Earth run
*-------------------------------------------------------
* Storage of Cartesian components in BFC cases
STORE(UCRT,VCRT,WCRT)
NCRT=1
* Storage for interphase sources of momentum
STORE(MOMY);STORE(MOMZ)
* Storage for cell residence-time
STORE(REST)
* Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,5,FIXFLU,GRND)
COVAL(GENPAT,7,FIXFLU,GRND)
* Relaxation of interphase sources
RELAX(MOMY,LINRLX,7.0000E-01)
RELAX(MOMZ,LINRLX,7.0000E-01)
* Activation of GENTRA GROUND-station
L(GENSET
* ----------------------------------------------------------
* GENTRA GROUP 7: GENTRA test
* ----------------------------------------------------------
* GENTRA auto test RG(51)--RG(56) are reserved
RG(51)=1.129E+00;RG(52)=1.259E+03;RG(53)=4.518E-03
RG(54)=2.500E+00;RG(55)=5.020E-01;RG(56)=3.125E+00
RG(57)=5.000E+16;RG(58)=5.600E+00
EX(P1)=5.056E+03;EX(V1)=1.001;EX(W1)=3.533
EX(WCRT)=3.247
EX(VCRT)= 1.003E+00