PHOTON USE
  p
 
 
 
 
  text
  1
  1
  Particles are bounced off a heated solid surface
   0.77542E+03 0.26315E+04 CR
  text
  1
  1
  in this CONJUGATE HEAT TRANSFER case.
   0.77542E+03 0.25012E+04 CR
  upause 5
  gr out x 1
  gr out x 1 y 5 6 z 4 5 col 15
  text
  1
  15
  SOLID
   0.12753E+04 0.11268E+04 CR
  GEOMETRY READ
  T00001
  upause 1
  GEOMETRY READ
  T00002
  upause 1
  GEOMETRY READ
  T00003
  upause 1
  GEOMETRY READ
  T00004
  upause 1
  GEOMETRY READ
  T00005
  upause 1
  GEOMETRY READ
  T00006
  upause 1
  GEOMETRY READ
  T00007
  upause 1
  GEOMETRY READ
  T00008
  upause 1
  GEOMETRY READ
  T00009
  upause 1
  GEOMETRY READ
  T00010
  upause 1
  GEOMETRY READ
  T00011
  upause 1
  GEOMETRY READ
  T00012
  upause 1
  GEOMETRY READ
  T00013
  upause 5
  set vec ref 4.0
  vec x 1 y 1 m z 1 3
  vec x 1 y 1 4 z 4 6
  vec x 1 y 7 m z 4 6
  vec x 1 y 1 m z 6 m
  upause 2
  text
  3
 
  Press e to end
   0.77542E+03 0.20265E+04 CR
  enduse
 
    GROUP 1. Run title
TEXT(PARTICLES IN 2D CURVED DUCT: G506
TITLE
  DISPLAY
    This case deals with conjugate heat transfer between a metal
    solid which stands in the middle of a curved 2-D duct, and a
    cold air stream which flows through the duct,passing
    the 'obstacle'. The geometry is represented by means
    of body-fitted coordinates; Both air and liquid droplets
    are heated by the solid.
  ENDDIS
    GROUP 2. Transience; time-step specification
    GROUP 3. X-direction grid specification
    GROUP 4. Y-direction grid specification
NY=10;YVLAST=.05
    GROUP 5. Z-direction grid specification
NZ=10;ZWLAST=.20
    GROUP 6. Body-fitted coordinates or grid distortion
REAL(PI,YY0,ZZ0,ARC)
PI=3.14159;ARC=PI/3.
YY0=0.05*(3.0)**0.5
ZZ0=0.05
BFC=T
DOMAIN(1,2,1,1,1,NZ+1)
SETLIN(YC,.1*SIN(ARC+ARC*LNK)-YY0)
SETLIN(ZC,ZZ0-.1*COS(ARC+ARC*LNK))
DOMAIN(1,2,NY+1,NY+1,1,NZ+1)
SETLIN(YC,.15*SIN(ARC+ARC*LNK)-YY0)
SETLIN(ZC,ZZ0-.15*COS(ARC+ARC*LNK))
DOMAIN(1,2,1,NY+1,1,1)
SETLIN(YC,LNJ*(YL-YF)+YF)
SETLIN(ZC,LNJ*(ZL-ZF)+ZF)
DOMAIN(1,2,1,NY+1,NZ+1,NZ+1)
SETLIN(YC,LNJ*(YL-YF)+YF)
SETLIN(ZC,LNJ*(ZL-ZF)+ZF)
DOMAIN(1,NX+1,1,NY+1,1,NZ+1)
MAGIC(T)
    GROUP 7. Variables stored, solved & named
NAME(C6)=VAPO
PRESS0=1.E+05
SOLVE(P1,V1,W1,VAPO,TEM1)
SOLUTN(P1,Y,Y,Y,N,N,N)
STORE(PRPS,BLOK,DEN1,VISL,VCRT,WCRT)
    GROUP 8. Terms (in differential equations) & devices
IVARBK=50;ISOLBK=1
TERMS(TEM1,N,Y,Y,Y,P,P)
REAL(CPGAS);CPGAS=4005
    GROUP 9. Properties of the medium (or media)
    GROUP 10. Inter-phase-transfer processes and properties
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=1.0;FIINIT(PRPS)=0.0
FIINIT(BLOK)=1.0;FIINIT(VAPO)=0.0
FIINIT(TEM1)=293.0
PATCH(SOLIDS,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(SOLIDS,P1,0.0,0.0)
COVAL(SOLIDS,V1,0.0,0.0)
COVAL(SOLIDS,W1,0.0,1.0)
COVAL(SOLIDS,TEM1,0.0,293.0)
PATCH(SOLID,INIVAL,1,1,5,6,4,5,1,1)
COVAL(SOLID,PRPS,0.0,103)
COVAL(SOLID,BLOK,0.0,2.0)
COVAL(SOLID,TEM1,0.0,793.0)
CONPOR(SOLID,-1,CELL,1,1,-5,-6,-4,-5)
    GROUP 12. Unused
    GROUP 13. Boundary conditions and special sources
INLET(GXINL,LOW,1,NX,1,NY,1,1,1,1)
VALUE(GXINL,P1,1.29);VALUE(GXINL,W1,1.00)
VALUE(GXINL,TEM1,293.00);VALUE(GXINL,VAPO,0.0)
PATCH(HOTHOT,CELL,1,NX,5,6,4,5,1,1)
COVAL(HOTHOT,TEM1,FIXFLU,793.0)
PATCH(GXOUT,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(GXOUT,P1,1.0E5,0.0)
COVAL(GXOUT,W1,ONLYMS,SAME)
WALL(WALL1,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(WALL1,W1,1.0,0.0)
WALL(WALL2,SOUTH,1,NX,1,1,1,NZ,1,1)
COVAL(WALL2,W1,1.0,0.0)
    GROUP 14. Downstream pressure for PARAB=.TRUE.
    GROUP 15. Termination of sweeps
LSWEEP=70;LITER(TEM1)=100
RELAX(P1,LINRLX,0.70);RELAX(V1,LINRLX,0.5)
RELAX(W1,LINRLX,0.5);RELAX(TEM1,LINRLX,0.5)
RELAX(VAPO,LINRLX,0.7)
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
    GROUP 18. Limits on variables or increments to them
    GROUP 19. Data communicated by satellite to GROUND
    GROUP 20. Preliminary print-out
    GROUP 21. Print-out of variables
    GROUP 22. Spot-value print-out
IYMON=8;IZMON=7
    GROUP 23. Field print-out and plot control
TRACE=F
#GENTRA
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 60 - VAPORISING_DROPLETS
GPTYPE=60
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=-9.80;GGRAZ=0;GBUOYA=T
  * Stochastic model of turbulence
GSTOCH=F
  * Data for vaporising particles
GDRAG=GRND1;GNUSS=GRND1
GKONV=GRND1;GCPCON=CPGAS;GCPVAP=GRND1
GMWCON=29.0;GMWVAP=18.0;GLIQST=GRND1
GCPLIQ=10000;GHLIQD=GRND1;GVAPST=GRND1
GLATVP=100000;GDTRCT=0
  *-------------------------------------------------------
  * 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_LIQDEN_FRATE_TEMP_(NUMB)
   0.0001 0.0001 3.2 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 3.1 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 3.0 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.8 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.6 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.4 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.2 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.15 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.1 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.05 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 2.0 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 1.95 1.0 9.e-5 1000.0 1.e-6 293.0
   0.0001 0.0001 1.90 1.0 9.e-5 1000.0 1.e-6 293.0
  
  * Wall treatment, and rest coefficient if appropriate
GWALLC=3;GWREST=1
  * Porosity threshold
GPOROS=0
  *-------------------------------------------------------
  * GENTRA GROUP 3: Numerical controls
  *-------------------------------------------------------
  * 1st GENTRA sweep; frequency of calls
GSWEP1=1;GSWEPF=10
  * Maximum Lagrangian time-step; time step size multplier
GDTMAX=0.002;GRTFRC=50000.0
  * Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=1000;GTIMMX=10
  *-------------------------------------------------
  * GENTRA GROUP 4: Output controls
  *-------------------------------------------------------
  * Restart-file, history-file and frequency for output
GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1
GOUTFR=4
  * Particle number for residence-time calculation
GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: Provisions for the GENTRA-Earth run
  *-------------------------------------------------------
  * Storage of Cartesian components in BFC cases
NCRT=1
  * Storage for interphase sources of momentum
STORE(MOMY)
STORE(MOMZ)
  * Storage for interphase source of heat
STORE(HEAT)
  * Storage for interphase source of mass and for
    vapour mass-fraction
STORE(MASS)
PRNDTL(VAPO)=1;PRT(VAPO)=1
  * Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,V1,FIXFLU,GRND)
COVAL(GENPAT,W1,FIXFLU,GRND)
COVAL(GENPAT,TEM1,FIXFLU,GRND)
COVAL(GENPAT,VAPO,FIXFLU,GRND)
  * Interphase sources of mass
PATCH(GENMAS,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENMAS,P1,FIXFLU,GRND)
  * Relaxation of interphase sources
GLNRLX=1.0
RELAX(MOMY,LINRLX,GLNRLX)
RELAX(MOMZ,LINRLX,GLNRLX)
RELAX(HEAT,LINRLX,GLNRLX)
RELAX(MASS,LINRLX,GLNRLX)
  * Activation of GENTRA GROUND-station
L(GENSET
  * ----------------------------------------------------------
  * GENTRA GROUP 7: GENTRA test
  * ----------------------------------------------------------
   * GENTRA auto test  RG(51)--RG(58) are reserved
RG(51)=1.572E+00;RG(52)=1.165E+03;RG(53)=1.748E-02
RG(54)=6.500E+00;RG(55)=2.902E-01;RG(56)=1.468E+00
RG(57)=3.888E+00;RG(58)=1.045E+00