#cls
text(In-Form Equivalent Of PLANT Case Z616
title
libref=781
DISPLAY
This In-Form case does what Z616 does, but without fortran coding
or re-compilation.
In case z616, PLANT was used to introduce a turbulence
model for concentration fluctuations which employs no
conservation equations for statistical properties of
the fluctuations. It is a variant of Multi-Fluid
concept of Brian Spalding as employed by Sergei Zhubrin.
The 17-fluid model is considered here to simulate the
turbulent mixing resulting from the admission of two
separate, isothermal coaxial jets of different
composition into a concentric duct as depicted
diagrammatically below.
ENDDIS
Entrance /////////////// Duct wall /////////////////Exit
-> :::::::::::: ->
surrounding -> :::::::::: -->
jet ------->:::::Mixing layer spreading to wall--->
" -------> ---->
central jet ------->-.-.-.-.- Symmetry axis .-.-.-.-.-.----->
In parallel with above model, the conservation equation
for the transport of the square of fluctuating
concentration component is also PLANTed for comparative
purposes.
At the final stage of computations the results are
processed to get the Fluid Population Distribution
histogram which may be viewed by PHOTON. The use file
for the latter is supplied.
The example employed the following PLANT features:
- specific sources introduction,
- reference residuals calculations,
- intervention in calculation,
- processing the results and
- special print-out preparations,
PLANT information :
* Data input groups used: 13, 19
* Ground groups planted : 13, 19-6
* Headings used : SORC??, SC06??
* Functions used : SUM
* Commands used : IF, REGION
The following In-Form settings are equivalent to those of PLANT
library. They make use of the In-Form functions:
MAKE, STORE1, STORED and SOURCE.
Please inspect the bottom of the Q1 file so see how this is done.
ENDDIS
load(z616)
text(17-fluid turbulence model
libref=781
INFORM13BEGIN
** The next line de-activates PLANT **
namsat=none
SORG=SKIP
WG=SKIP
MIX=SKIP
** In-Form formulae which are equivalent to the
PLANT formulae of z616 **
** Source term for g
PATCH(ISORG,VOLUME,1,NX,1,NY,1,NZ,1,1)
(SOURCE of G at ISORG is 2.0*:RHO1:*EPKE*(GENG/(2.0*:RHO1:*EPKE+$
TINY)-G))
PATCH(IWG,VOLUME,1,NX,NY,NY,1,NZ,1,1)
(SOURCE of G at IWG is GENG/(2.0*:RHO1:*EPKE+TINY) with FIXVAL)
** Provide the re-calculation of reference residuals for G
(MAKE RES)
(STORE1 of RES is SUM(VOL*(GENG-2.*:RHO1:*EPKE*G)/(NY*NZ)))
(STORE1 of RESREF(G) is RES)
**Calculation of GENG
PATCH(PAT1,CELL,1,NX,1,NY,1,NZ-1,1,1)
(STORED of DFZ at PAT1 is ((H1[,,+1]-H1)/DZG)^2)
PATCH(PAT2,CELL,1,1,1,NY-1,1,NZ,1,1)
(STORED of DFY at PAT2 is ((H1[,+1,]-H1)/DYG)^2)
PATCH(PAT3,CELL,1,NX,1,NY,NZ,NZ,1,1)
(STORED of DFZH at PAT3 is ((H1-H1[,,-1])/DZG[,,-1])^2)
PATCH(PAT4,CELL,1,NX,NY,NY,1,NZ,1,1)
(STORED of DFYN at PAT4 is ((H1-H1[,-1,])/DYG[,-1,])^2)
(STORED of GENG is 2.8*:RHO1:*ENUT*(DFZ+DFY+DFZH+DFYN))
!!!!! note the use of the IF condition !!!!!
**Output calculations
(STORED of GG is SQRT(G) with IF(ISWEEP.EQ.LSWEEP))
(STORED of GGF is GG/(H1+TINY) with IF(ISWEEP.EQ.LSWEEP))
** Coupling/splitting rates
PATCH(iMIX,PHASEM,1,NX,1,NY,1,NZ,1,1)
* Fluid 1
(SOURCE of F1 at iMIX is :MMC:*EPKE*$
(F3+F5+F7+F9+F11+F13+F15+F17)*(0.-F1) with LINE)
* Fluid 2
(SOURCE of F2 at iMIX is 2.*:MMC:*EPKE*(F1*F3)-$
:MMC:*EPKE*(F4+F6+F8+F10+F12+F14+F16)*F2)
* Fluid 3
(SOURCE of F3 at iMIX is 2.*:MMC:*EPKE*(F2*F4+F1*F5)-$
:MMC:*EPKE*(F1+F17+F5+F7+F9+F11+F13+F15)*F3)
* Fluid 4
(SOURCE of F4 at iMIX is 2.*:MMC:*EPKE*(F3*F5+F2*F6+F1*F7)-$
:MMC:*EPKE*(F2+F6+F8+F10+F12+F14+F16)*F4)
* Fluid 5
(SOURCE of F5 at iMIX is 2.*:MMC:*EPKE*(F4*F6+F3*F7+F2*F8+F1*F9)-$
:MMC:*EPKE*(F1+F3+F17+F7+F9+F11+F13+F15)*F5)
* Fluid 6
(SOURCE of F6 at iMIX is 2.*:MMC:*EPKE*(F5*F7+F4*F8+F3*F9+$
F2*F10+F1*F11)-:MMC:*EPKE*(F2+F4+F8+F10+F12+F14+F16)*F6)
* Fluid 7
(SOURCE of F7 at iMIX is 2.*:MMC:*EPKE*(F6*F8+F5*F9+F4*F10+F3*F11$
+F2*F12+F1*F13)-:MMC:*EPKE*(F1+F3+F5+F17+F9+F11+F13+F15)*F7)
* Fluid 8
(SOURCE of F8 at iMIX is 2.*:MMC:*EPKE*(F7*F9+F6*F10+F5*F11+F4*F12+$
F3*F13+F2*F14+F1*F15)-:MMC:*EPKE*(F2+F4+F6+F10+F12+F14+F16)*F8)
* Fluid 9
CHAR(SUM1,SUM2)
SUM1=(F8*F10+F7*F11+F6*F12+F5*F13+F4*F14+F3*F15+F2*F16+F1*F17)
SUM2=(F1+F3+F5+F7+F17+F11+F13+F15)
(SOURCE of F9 at iMIX is 2.*:MMC:*EPKE*:SUM1:-:MMC:*EPKE*:SUM2:*$
F9)
* Fluid 10
SUM1=(F9*F11+F8*F12+F7*F13+F6*F14+F5*F15+F4*F16+F3*F17)
SUM2=(F2+F4+F6+F8+F12+F14+F16)
(SOURCE of F10 at iMIX is 2.*:MMC:*EPKE*:SUM1:-$
:MMC:*EPKE*:SUM2:*F10)
* Fluid 11
(SOURCE of F11 at iMIX is 2.*:MMC:*EPKE*(F10*F12+F9*F13+F8*F14+F7*F$
15+F6*F16+F5*F17)-:MMC:*EPKE*(F1+F3+F5+F7+F9+F17+F13+F15)*F11)
* Fluid 12
(SOURCE of F12 at iMIX is 2.*:MMC:*EPKE*(F11*F13+F10*F14+F9*F15+$
F8*F16+F7*F17)-:MMC:*EPKE*(F2+F4+F6+F8+F10+F14+F16)*F12)
* Fluid 13
(SOURCE of F13 at iMIX is 2.*:MMC:*EPKE*(F12*F14+F11*F15+$
F10*F16+F9*F17)-:MMC:*EPKE*(F1+F3+F5+F7+F9+F11+F15+F17)*F13)
* Fluid 14
(SOURCE of F14 at iMIX is 2.*:MMC:*EPKE*(F13*F15+F12*F16+F11*F17)-$
:MMC:*EPKE*(F2+F4+F6+F8+F10+F12+F16)*F14)
* Fluid 15
(SOURCE of F15 at iMIX is 2.*:MMC:*EPKE*(F14*F16+F13*F17)-$
:MMC:*EPKE*(F1+F3+F5+F7+F9+F11+F13+F17)*F15)
* Fluid 16
(SOURCE of F16 at iMIX is 2.*:MMC:*EPKE*(F15*F17)-$
:MMC:*EPKE*(F2+F4+F6+F8+F10+F12+F14)*F16)
* Fluid 17
(SOURCE of F17 at iMIX is :MMC:*EPKE*(F1+F3+F5+F7+F9+F11+F13+F15$
)*(0.-F17) with LINE)
INFORM13END
INFORM7BEGIN
** Output calculations
SUM1=16./16.*F1+15./16.*F2+14./16.*F3+13./16.*F4
SUM2=12./16.*F5+11./16.*F6+10./16.*F7+9./16.*F8+8./16.*F9
(STORED of CAV is :SUM1:+:SUM2: with IF(ISWEEP.EQ.LSWEEP))
SUM1=7./16.*F10+6./16.*F11+5./16.*F12+4./16.*F13
SUM2=3./16.*F14+2./16.*F15+1./16.*F16+0./16.*F17
(STORED of CAV is CAV+:SUM1:+:SUM2: with IF(ISWEEP.EQ.LSWEEP))
(STORED of MAS is F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+$
F13+F14+F15+F16+F17 with IF(ISWEEP.EQ.LSWEEP))
CHAR(CH1,CH2,CH3,CH4,CH5,CH6,CH7,CH8,CH9,CH10)
CHAR(CH11,CH12,CH13,CH14,CH15,CH16,CH17)
CH1=ABS(CAV-16./16)*F1
CH2=ABS(CAV-15./16.)*F2
CH3=ABS(CAV-14./16)*F3
CH4=ABS(CAV-13./16.)*F4
CH5=ABS(CAV-12./16)*F5
CH6=ABS(CAV-11./16.)*F6
CH7=ABS(CAV-10./16)*F7
CH8=ABS(CAV- 9./16.)*F8
CH9=ABS(CAV- 8./16)*F9
CH10=ABS(CAV- 7./16.)*F10
CH11=ABS(CAV- 6./16)*F11
CH12=ABS(CAV- 5./16.)*F12
CH13=ABS(CAV- 4./16)*F13
CH14=ABS(CAV- 3./16.)*F14
CH15=ABS(CAV- 2./16)*F15
CH16=ABS(CAV- 1./16.)*F16
CH17=ABS(CAV- 0./16)*F17
(STORED GAV is :CH1:+:CH2:+:CH3:+:CH4:+:CH5:+:CH6:)
(STORED GAV is GAV+:CH7:+:CH8:+:CH9:+:CH10:+:CH11:+:CH12:)
(STORED GAV is GAV+:CH13:+:CH14:+:CH15:+:CH16:+:CH17:)
(STORED of GF is GAV/(CAV+TINY) with IF(ISWEEP.EQ.LSWEEP))
** Output data processing for plotting PDF
PATCH(PATCH1,CELL,1,1,1,NY,1,1,1,1)
SUM1=F1[1,4,4]*MAX(ABS(F1[1,4,4]-YV)/(F1[1,4,4]-YV),0.)
(STORED of FPD at PATCH1 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH2,CELL,1,1,1,NY,2,2,1,1)
SUM1=F2[1,4,4]*MAX(ABS(F2[1,4,4]-YV)/(F2[1,4,4]-YV),0.0)
(STORED of FPD at PATCH2 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH3,CELL,1,1,1,NY,3,3,1,1)
SUM1=F3[1,4,4]*MAX(ABS(F3[1,4,4]-YV)/(F3[1,4,4]-YV),0.0)
(STORED of FPD at PATCH3 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH4,CELL,1,1,1,NY,4,4,1,1)
SUM1=F4[1,4,4]*MAX(ABS(F4[1,4,4]-YV)/(F4[1,4,4]-YV),0.0)
(STORED of FPD at PATCH4 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH5,CELL,1,1,1,NY,5,5,1,1)
SUM1=F5[1,4,4]*MAX(ABS(F5[1,4,4]-YV)/(F5[1,4,4]-YV),0.0)
(STORED of FPD at PATCH5 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH6,CELL,1,1,1,NY,6,6,1,1)
SUM1=F6[1,4,4]*MAX(ABS(F6[1,4,4]-YV)/(F6[1,4,4]-YV),0.0)
(STORED of FPD at PATCH6 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH7,CELL,1,1,1,NY,7,7,1,1)
SUM1=F7[1,4,4]*MAX(ABS(F7[1,4,4]-YV)/(F7[1,4,4]-YV),0.0)
(STORED of FPD at PATCH7 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH8,CELL,1,1,1,NY,8,8,1,1)
SUM1=F8[1,4,4]*MAX(ABS(F8[1,4,4]-YV)/(F8[1,4,4]-YV),0.0)
(STORED of FPD at PATCH8 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH9,CELL,1,1,1,NY,9,9,1,1)
SUM1=F9[1,4,4]*MAX(ABS(F9[1,4,4]-YV)/(F9[1,4,4]-YV),0.0)
(STORED of FPD at PATCH9 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH10,CELL,1,1,1,NY,10,10,1,1)
SUM1=F10[1,4,4]*MAX(ABS(F10[1,4,4]-YV)/(F10[1,4,4]-YV),0.0)
(STORED of FPD at PATCH10 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH11,CELL,1,1,1,NY,11,11,1,1)
SUM1=F11[1,4,4]*MAX(ABS(F11[1,4,4]-YV)/(F11[1,4,4]-YV),0.0)
(STORED of FPD at PATCH11 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH12,CELL,1,1,1,NY,12,12,1,1)
SUM1=F12[1,4,4]*MAX(ABS(F12[1,4,4]-YV)/(F12[1,4,4]-YV),0.0)
(STORED of FPD at PATCH12 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH13,CELL,1,1,1,NY,13,13,1,1)
SUM1=F13[1,4,4]*MAX(ABS(F13[1,4,4]-YV)/(F13[1,4,4]-YV),0.0)
(STORED of FPD at PATCH13 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH14,CELL,1,1,1,NY,14,14,1,1)
SUM1=F14[1,4,4]*MAX(ABS(F14[1,4,4]-YV)/( F14[1,4,4]-YV),0.0)
(STORED of FPD at PATCH14 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH15,CELL,1,1,1,NY,15,15,1,1)
SUM1=F15[1,4,4]*MAX(ABS(F15[1,4,4]-YV)/(F15[1,4,4]-YV),0.0)
(STORED of FPD at PATCH15 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH16,CELL,1,1,1,NY,16,16,1,1)
SUM1=F16[1,4,4]*MAX(ABS(F16[1,4,4]-YV)/(F16[1,4,4]-YV),0.0)
(STORED of FPD at PATCH16 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
PATCH(PATCH17,CELL,1,1,1,NY,17,17,1,1)
SUM1=F17[1,4,4]*MAX(ABS(F17[1,4,4]-YV)/(F17[1,4,4]-YV),0.0)
(STORED of FPD at PATCH17 is :SUM1: with IF(ISWEEP.EQ.LSWEEP))
INFORM7END
DISTIL=T
EX(P1)=5.613E-01; EX(V1)=2.532E-02; EX(W1)=2.840E+00
EX(KE)=4.442E-01; EX(EP)=6.611E+00; EX(H1)=1.939E-01
EX(FPD)=2.724E-02; EX(GF)=4.936E-01; EX(GAV)=5.976E-02
EX(MAS)=1.000E+00; EX(CAV)=1.939E-01; EX(F17)=1.402E-01
EX(F16)=7.228E-02; EX(F15)=2.095E-01; EX(F14)=3.191E-01
EX(F13)=1.075E-01; EX(F12)=4.737E-02; EX(F11)=2.690E-02
EX(F10)=1.746E-02; EX(F9)=1.188E-02; EX(F8)=8.836E-03
EX(F7)=6.076E-03; EX(F6)=5.269E-03; EX(F5)=3.765E-03
EX(F4)=3.483E-03; EX(F3)=2.355E-03; EX(F2)=3.138E-03
EX(F1)=1.499E-02; EX(GGF)=2.830E+13; EX(GG)=5.768E-02
EX(DFYN)=3.957E-01; EX(DFZH)=3.238E-06; EX(DFY)=1.765E+01
EX(DFZ)=4.226E-02; EX(GENG)=2.216E-01; EX(EPKE)=8.359E+00
EX(GEN1)=1.724E+03; EX(LEN1)=1.567E-02; EX(ENUT)=4.963E-03
EX(G)=5.591E-03
libref=781
PHOTON USE
P;phi; 10 1;;
msg( Single fluid central-jet concentration contours
con h1 x 1 fil;.001
pause
con cl; red
msg( Averaged 17-fluid contours
con cav x 1 fil;.001
pause
con cl; red
msg( Concentration fluctuations by transport equation
con gg x 1 fil;.001
pause
con cl; red
msg( Averaged 17-fluid concentration fluctuation
con gav x 1 fil;.001
msg
msg Hit Enter for FPD hystogram
msg
pause
p
20 1
con fpd x 1 fil;10.
msg Hit Enter to continue
ENDUSE
STOP