PHOTON USE
p
0.20443E+04 0.15633E+04 CR
use patgeo
use genuse
@
0.51268E+03 0.55365E+03 CR
Trajectories of the first 20 rain droplets.@
@
0.20379E+04 0.29962E+03 R
RETURN to continue.@
pause
vec x 1
set vec ref
50
red
vec x 10
vie x
@
0.41268E+03 0.25365E+04 CR
View from X direction@
@
0.41268E+03 0.20365E+04 CR
Press e to END.@
ENDUSE
GROUP 1. Run title
TEXT(Sample cup in rain: G207
TITLE
DISPLAY
The continous-phase set-up is that in the PHOENICS Library
Case 237.
A simulation is made of the flow in and around a shallow
open-topped hollow cylinder with a crosswind blowing
across its top, normal to the cylinder axis. The height
of the cup to its radius is 0.4 . This geometry is a design
considered for a sampling cup used for measuring atmospheric
precipitation.
Cylindrical polar coordinates are used in 3 dimensions.
The polar axis points vertically upwards from the ground
surface, and the vertical extent of the domain extends well
beyond the top of the sampling cup. The radial extent of the
domain extends beyond the rim of the cup.
Particles (say rain droplets) enter the computational domain
at the higher end, falling in a gravitational field with an
initial velocity of 0.5 m/s, and impinging on the sample cup
and the surrounding ground. The particles are removed from
the domain as they impinge.
ENDDIS
=============================================================
REAL(PI,DY,AA);PI=3.14159
INTEGER(XFLOW1,XFLOW2,JJ1);XFLOW1=6;XFLOW2=15
GROUP 2. Transience; time-step specification
GROUP 3. X-direction grid specification
** The domain is extended over 360 degrees in this example.
CARTES=F;XCYCLE=T;GRDPWR(X,20,2.0*PI,1.0)
GROUP 4. Y-direction grid specification
YVLAST=0.35;NY=21
**Inside the cup 10 radial layers of cells are used, and
beyond the rim of the cup 11 radial layers of cells are
used. A geometric-progression grid spacing distribution
is employed both sides of the cup wall to ensure adequate
resolution of the flow at the exterior- and interior-side
of the cup wall.
YFRAC(10)=1.0;DY=.03;AA=1.2
DO JJ=9,1,-1
+ JJ1=JJ+1
+ YFRAC(JJ)=YFRAC(JJ1)-DY;DY=DY*AA
ENDDO
DY=.03;AA=1.4
DO JJ=11,19
+ JJ1=JJ-1;YFRAC(JJ)=YFRAC(JJ1)+DY;DY=DY*AA
ENDDO
YFRAC(20)=3.5;YFRAC(21)=4.0
GROUP 5. Z-direction grid specification
GRDPWR(Z,9,0.12,1.)
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1);SOLUTN( P1,Y,Y,Y,N,N,N )
GROUP 9. Properties of the medium (or media)
RHO1=1.205;ENUL=1.8E-5/RHO1;ENUT=ENUL*10.
GROUP 11. Initialization of variable or porosity fields
**The wall of the cup is represented by setting zero
porosities at the cylindrical surface IY=10 . The minus
signs preceeding the number 10 flags the activation of wall
friction at the inner and outer surface of the cup wall.
CONPOR(0.0,NORTH,1,NX,-10,-10,1,NZ/3)
GROUP 13. Boundary conditions and special sources
** The inlet region extends from pi/2 to 3pi/2. The wind-
speed at the cylindrical inlet surface is set in POLRA
which is resolved in subroutine GXPLOR into the
tangential and radial directions.
POLRA=4.0
PATCH(UPOL,CELL,XFLOW1-1,XFLOW2,NY,NY,1,NZ,1,1)
COVAL(UPOL,U1,FIXVAL,GRND1)
PATCH(VPOL,CELL,XFLOW1,XFLOW2,NY-1,NY-1,1,NZ,1,1)
COVAL(VPOL,V1,FIXVAL,GRND1)
** The pressure is fixed to zero around the entire outer
boundary of the domain...
PATCH(PINF,CELL,1,NX,NY,NY,1,NZ,1,1)
COVAL(PINF,P1,1.0E6,0.0 )
PATCH(GX1,HIGH,1,NX,1,NY,NZ,NZ,1,1)
PATCH(GX2,NORTH,1,NX,NY,NY,1,NZ,1,1)
GROUP 15. Termination of sweeps
LSWEEP=20
GROUP 16. Termination of iterations
LITER(P1) = 12;LITER(U1)=1;LITER(V1)=1
GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,0.003);RELAX(V1,FALSDT,0.003)
RELAX(W1,FALSDT,0.003)
GROUP 22. Spot-value print-out
IXMON=2;IYMON=12;IZMON=4
GROUP 23. Field print-out and plot control
NPLT=2;NXPRIN=NX/5;NYPRIN=NY/5
PATCH(MAP,CONTUR,1,NX,1,NY,3,3,1,1)
PLOT(MAP,P1,1.0,10.0);PLOT(MAP,U1,1.0,10.0)
PLOT(MAP,V1,1.0,10.0);PLOT(MAP,W1,1.0,10.0)
NOWIPE=T
TRACE=F
#GENTRA
*-------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*-------------------------------------------------------
* Particle type - 30 - ISOTHERMAL_PARTICLES
GPTYPE=30
* Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=0;GGRAZ=-9.8000E+00;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
*__XP__YP__ZP__UP__VP__WP__DIAM_DENSTY_FRATE
0.0 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
0.0 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
0.0 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
0.0 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
0.0 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
0.0 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
1.5 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.0 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
2.5 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.0 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
3.5 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.0 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.1 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.15 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.2 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.25 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.3 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
4.5 0.35 0.12 0.0 0.0 -0.5 0.001 1000.0 1.0E-5
* Wall treatment, and rest coefficient if appropriate
GWALLC=1
* Porosity threshold
GPOROS=0
*-------------------------------------------------------
* GENTRA GROUP 3: Numerical controls
*-------------------------------------------------------
* 1st GENTRA sweep; frequency of calls
GSWEP1=LSWEEP-5;GSWEPF=1
* Maximum Lagrangian time-step
GDTMAX=1
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=10000;GTIMMX=0
*-------------------------------------------------------
* GENTRA GROUP 4: Output controls
*-------------------------------------------------------
* Trajectory-file identifier, history-file identifier
and frequency for output
GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1
* Particle number for residence-time calculation
GRESTI=0
*-------------------------------------------------------
* GENTRA GROUP 5: Provisions for the GENTRA-Earth run
*-------------------------------------------------------
* Storage for interphase sources of momentum
STORE(MOMX)
STORE(MOMY)
STORE(MOMZ)
* Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,3,FIXFLU,GRND)
COVAL(GENPAT,5,FIXFLU,GRND)
COVAL(GENPAT,7,FIXFLU,GRND)
* Relaxation of interphase sources
RELAX(MOMX,LINRLX,7.0000E-01)
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)=5.817E+00;RG(52)=2.129E+03;RG(53)=1.313E-01
RG(54)=1.450E+01;RG(55)=1.758E+01;RG(56)=3.997E-02
RG(57)=4.800E+17;RG(58)=4.800E+01
EX(P1)=3.035E+01;EX(U1)=1.823;EX(V1)=2.010;EX(W1)=2.477E-01
EX(NPOR)=9.841E-01