LIBRARY:
CEL:
EXPRESSIONS:
AD = (1-Phi morph)*ADBubb+Phi morph*ADSurf
ADBubb = 6*max(Gas3.Volume Fraction, minVF)/DBubb
ADDrop = 10e-4[m^-1]
ADSX = Gas3.Volume Fraction.Gradient X/ADSurf
ADSY = Gas3.Volume Fraction.Gradient Y/ADSurf
ADSZ = Gas3.Volume Fraction.Gradient Z/ADSurf
ADSurf = max(ADDrop, (2*Gas3 VF Grad*Liquid VF Grad)/(Gas3 VF \
Grad+Liquid VF Grad))
Ain = 0.000488723 [m^2]
AlphaCritBubb = 0.3
AlphaCritDrop = 0.3
AlphaGasSchlange = min( max(Gas3.Volume Fraction, AlphaCritBubb - \
DeltaAlphaBubb), AlphaCritBubb + DeltaAlphaBubb)
AlphaLiqSchlange = min( max(Liquid.Volume Fraction, AlphaCritDrop - \
DeltaAlphaDrop), AlphaCritDrop + DeltaAlphaDrop)
CD = (1-Phi morph)*CDBubb+Phi morph*CDSurf
CD fs = min(1.0E+5, max(1.0E-2, (Liquid.Volume Fraction*tlfs + \
Gas3.Volume Fraction*tgfs)/max(1.0E-3 [kg m^-1 s^-2], DenLiq*uslip^2)))
CD1 = 1
CD2 = 1
CDBubb = if(Cdsphere>=Cdellipse,Cdsphere,min(Cdellipse,Cdcap))
CDDrop = 0.44
CDSurf = max(CD fs*Phi surf, CDDrop*(1-Phi surf))
CFX = (1[m^2 s^-2])*(Liquid.Volume Fraction.Gradient X)*(DenLiq)*(Phi \
clust)
CFY = (1[m^2 s^-2])*(Liquid.Volume Fraction.Gradient Y)*(DenLiq)*(Phi \
clust)
CFZ = (1[m^2 s^-2])*(Liquid.Volume Fraction.Gradient Z)*(DenLiq)*(Phi \
clust)
CL = (CLcg)*(-0.27)
CLcg = (if(Phi surf < 0.5, 1, if(Phi surf >= 0.5 && Phi surf < 1, \
(0.5*tanh(10*(0.6-Phi surf))+0.5), 0)))*(if(Phi surf ==1 || \
Gas3.Volume Fraction > 0.5, 0, 1))
Cdcap = (8/3)
Cdellipse = ((2/3)*(Eo^(1/2)))
Cdsphere = (24/Nre)*(1+(0.15*((Nre)^0.687)))
Cent = 0.1
CoalescCoeffwall = max(1,aq)
Coeff = 1.0
CondGas = 0.24780326e-01 [W m^-1 K^-1]
CondLiq = 0.67767141e+00 [W m^-1 K^-1]
ContactAngle = 120
CpGas = 0.20321079e+04 [J kg^-1 K^-1]
CpLiq = 0.42177068e+04 [J kg^-1 K^-1]
DBubb = Gas3.Mean Particle Diameter
DP = 25e-3 [m]
De = Lamb*DP
DeltaAlphaBubb = 0.05
DeltaAlphaDrop = 0.05
DeltaGradGas = 1/(n*delta x)*0.1
DenGas = 0.59579337e+00 [kg m^-3]
DenLiq = 0.95842977e+03 [kg m^-3]
DenMix = Liquid.Volume Fraction*Liquid.density + Gas3.Volume \
Fraction*Gas3.density
Dmean = dp1
EnthGas = 0.26755776e+07 [J kg^-1]
EnthLiq = 0.41867484e+06 [J kg^-1]
Eo = ((DenLiq-DenGas)*g*(DBubb^2))/SurfTen
EpsIN = 0.001
Gas3 VF Grad = max(sqrt((Gas3.Volume Fraction.Gradient \
X)^2+(Gas3.Volume Fraction.Gradient Y)^2+(Gas3.Volume \
Fraction.Gradient Z)^2),minVFGrad)
Gas3VEL = sqrt((Gas3.Velocity u)^2+(Gas3.Velocity v)^2+(Gas3.Velocity \
w)^2)
GasDTot = Gas1.Volume Fraction + Gas2.Volume Fraction
GasTot = GasDTot+Gas3.Volume Fraction
Ginm = VIN
GradAlphaGasMag = sqrt((Gas3.Volume Fraction.Gradient X)^2 + \
(Gas3.Volume Fraction.Gradient Y)^2 + (Gas3.Volume Fraction.Gradient \
Z)^2)
GradAlphaGasMagMin = max(sqrt((Gas3.Volume Fraction.Gradient X)^2 + \
(Gas3.Volume Fraction.Gradient Y)^2 + (Gas3.Volume Fraction.Gradient \
Z)^2), MinVFGrad)
GradAlphaLiqMag = sqrt((Liquid.Volume Fraction.Gradient X)^2 + \
(Liquid.Volume Fraction.Gradient Y)^2 + (Liquid.Volume \
Fraction.Gradient Z)^2)
GradCritGas = 1/(n*delta x)
GradGasSchlange = min( max(GradAlphaGasMag, GradCritGas - \
DeltaGradGas), GradCritGas + DeltaGradGas)
GradLiqSchlange = min( max(GradAlphaLiqMag, GradCritLiq - \
DeltaGradLiq), GradCritLiq + DeltaGradLiq)
GridCellSize = delta x
HighCoefTurb = 1.0E+2
InTempProf = (510-z*10/(0.5 [m])) [K]
Lamb = 0.6
Liquid VF Grad = max(sqrt((Liquid.Volume Fraction.Gradient \
X)^2+(Liquid.Volume Fraction.Gradient Y)^2+(Liquid.Volume \
Fraction.Gradient Z)^2), minVFGrad)
LiquidHoldup = volumeAve(Liquid.Volume Fraction)@tube
LiquidTemperaturepoint = Liquid.Temperature
LiquidVEL = sqrt((Liquid.Velocity u)^2+(Liquid.Velocity \
v)^2+(Liquid.Velocity w)^2)
Liquidout = areaAve(Liquid.Volume Fraction)@Out
MEntrain = PhiFS*Gas3.Volume Fraction*DenGas*Qg
MaxCN = maxVal(Courant Number)@tube
MinVFArea = 1.0E-7
MinVFforArea = 1.0E-7
MinVolFrac = 1.0E-10
NSD = NSD6
NSD1 = max(1 [m^-2], NSDRef*((min(Tsup,9[K])/NSDTRef)^NSDPow))
NSD2 = max(1 [m^-2], NSDRef*((min(Tsup,19[K])/NSDTRef)^NSDPow))
NSD3 = max(1 [m^-2], NSDRef*((min(Tsup,29[K])/NSDTRef)^NSDPow))
NSD4 = max(1 [m^-2], NSDRef*((min(Tsup,39[K])/NSDTRef)^NSDPow))
NSD5 = max(1 [m^-2], NSDRef*((min(Tsup,49[K])/NSDTRef)^NSDPow))
NSD6 = max(1 [m^-2], NSDRef*((Tsup/NSDTRef)^NSDPow))
NSDPoint1 = Gas1 | Liquid.Nucleation Site Density
NSDPoint2 = Gas2 | Liquid.Nucleation Site Density
NSDPow = 1.805
NSDRef = 1.5e+5 [m^-2]
NSDTRef = 10.0 [K]
Nre = (DenLiq *(max((uslip),1E-7[m s^-1]))*DBubb)/(viscl)
OmegaDampLiq = \
HighCoefTurb*6.0*(Liquid.viscosity/Liquid.density)/(0.075*delta x^2)
PHydro = (hout-z)*DenLiq*g
PRef = 101325 [Pa]
PRef1 = PRef-(DenLiq*g*hout)
PabsMax = 5.5e+6 [Pa]
PabsMin = 3.5e+6 [Pa]
Phi clust = (if(Phi surf < 0.5, 1, if(Phi surf >= 0.5 && Phi surf < 1, \
(0.5*tanh(10*(0.6-Phi surf))+0.5), 0)))*(if(Phi surf ==1 || \
Gas3.Volume Fraction > 0.5, 0, 1))
Phi coal = (if(Phi surf<0.5 && Gas3.Volume Fraction < 0.5,0,1))*(if(Phi \
surf>=0.5 && Phi surf<=1.0,0,1))*(if(Gas3.Volume \
Fraction>0.5,(5*tanh(-20*(0.5-Gas3.Volume Fraction))+0.5),0))
Phi morph = (if(Phi surf<0.5 && Gas3.Volume Fraction < \
0.5,0,1))*(if(Phi surf>=0.5 && Phi surf<=1.0 || Gas3.Volume Fraction \
> 0.5,(0.5*tanh(-5*(0.5-Gas3.Volume Fraction))+0.5),1))
Phi surf = PhiFS
PhiEnt = n*GridCellSize
PhiFS = 1 - PsiSurfMag
PhiSurf = 0.5*(1 + cos(pi*(GradGasSchlange - (GradCritGas - \
DeltaGradGas))/(2*DeltaGradGas)))
PsiSurf = PhiSurf*(WeightBubb - WeightDrop)
PsiSurfMag = if(PsiSurf>=0, PsiSurf, -PsiSurf)
QLiquidPoint = Gas1 | Liquid.Heat Flux to Liquid
QVapourPoint = Gas1 | Liquid.Heat Flux to Vapour
QavgLiquid = areaAve( Gas1 | Liquid.Heat Flux to Liquid)@HeatSurf
QavgVapour = areaAve( Gas1 | Liquid.Heat Flux to Vapour)@HeatSurf
Qconvg = Gas1 | Liquid.Heat Flux to Vapour -Qeva
Qconvl = Gas1 | Liquid.Heat Flux to Liquid - Qque1
Qeva = pi*1/6*dw^3*fr/4*NSD*DenGas*hfg
Qeva2 = min(1,a2)*2/3*dw*fr/4*DenGas*hfg
Qg = Cent/(g*PhiEnt)*Liquid.Turbulence Kinetic Energy*dun dn
Qque = Qque1+Qque2
Qque1 = hquenching * (dtl) *min(1,max(aq,0))
Qque2 = qwall*min(1,a2)-Qeva2
Qw = 350E+3 [W m^-2]
Qwliquid = Gas1 | Liquid.Heat Flux to Liquid
Qwtotal = (Qwvapor+Qwliquid)/2
Qwvapor = Gas1 | Liquid.Heat Flux to Vapour
RSLX1 = if(SLX1 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw1s/2))/SLX1 ,0 )))
RSLX2 = if(SLX2 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw2s/2))/SLX2 ,0 )))
RSLX3 = if(SLX3 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw3s/2))/SLX3 ,0 )))
RSLX4 = if(SLX4 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw4s/2))/SLX4 ,0 )))
RSLX5 = if(SLX5 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw5s/2))/SLX5 ,0 )))
RSLX6 = if(SLX6 < 0.00000001 [m], 0, min \
(1,max(((1/4*((NSD))^(-1/2))-(dw6s/2))/SLX6 ,0 )))
SLX1 = max(0.0000000000001 \
[m],dw1b*(SlidingL(Ginm,TsubL,T1a)-1/2)-dw1s/2)
SLX2 = max(0.0000000000001 \
[m],dw2b*(SlidingL(Ginm,TsubL,T2a)-1/2)-dw2s/2)
SLX3 = max(0.0000000000001 \
[m],dw3b*(SlidingL(Ginm,TsubL,T3a)-1/2)-dw3s/2)
SLX4 = max(0.0000000000001 \
[m],dw4b*(SlidingL(Ginm,TsubL,T4a)-1/2)-dw4s/2)
SLX5 = max(0.0000000000001 \
[m],dw5b*(SlidingL(Ginm,TsubL,T5a)-1/2)-dw5s/2)
SLX6 = max(0.0000000000001 \
[m],dw6b*(SlidingL(Ginm,TsubL,T6a)-1/2)-dw6s/2)
STX = ((SurfTen)*(kappa))*(Gas3.Volume Fraction.Gradient X)
STY = ((SurfTen)*(kappa))*(Gas3.Volume Fraction.Gradient Y)
STZ = ((SurfTen)*(kappa))*(Gas3.Volume Fraction.Gradient Z)
SmWaTurb = abs(PhiFS*2/3*DenLiq*VELBET*(0.5*(ubound^2-lbound^2)))
SurfTen = 0.072 [N m^-1]
T1 = 10 [K]
T1a = 5 [K]
T2 = 20 [K]
T2a = 15 [K]
T3 = 30 [K]
T3a = 25 [K]
T4 = 40 [K]
T4a = 35 [K]
T5 = 50 [K]
T5a = 45 [K]
T6 = 60 [K]
T6a = 55 [K]
TIN = Tsat - 0 [K]
TLiqsup = Liquid.Temperature-Tsat
TMax = 600 [K]
TMin = 450 [K]
TWall = Tsat + Tsup
TempLiquidDomain = volumeAve(Liquid.Temperature)@tube
TempLiquidOut = areaAve(Liquid.Temperature)@Out
Tsat = 0.37303432e+03 [K]
TsubL = max(0 [K],(Tsat - Liquid.Temperature))
TsubLc = Tsat - Liquid.Temperature
Tsup = 10 [K]
TsupPoint1 = Gas1 | Liquid.Temperature Superheating
TsupPoint2 = Gas2 | Liquid.Temperature Superheating
Tsupeff = if (aq<1, Tsupeff1, if(aq1+aq2+aq3+aq4+aq5<1,aq1*T1a + \
aq2*T2a + aq3*T3a +aq4*T4a +aq5*T5a + \
(1-(aq1+aq2+aq3+aq4+aq5))*T6a,if(aq1+aq2+aq3+aq4<1,aq1*T1a + aq2*T2a \
+ aq3*T3a +aq4*T4a + \
(1-(aq1+aq2+aq3+aq4))*T5a,if(aq1+aq2+aq3<1,aq1*T1a + aq2*T2a + \
aq3*T3a + (1-(aq1+aq2+aq3))*T4a,if(aq1+aq2<1,aq1*T1a + aq2*T2a + \
(1-(aq1+aq2))*T3a,if(aq1<1,aq1*T1a + (1-(aq1))*T2a,T1a)) ))))
Tsupeff1 = min(1 [K],if(Tsup < T1, a21*T1a+ \
(aq-a2)*((Tsup+T1a)/2),if(Tsup < T2,(a21*T1a + a22*T2a)+ \
(aq-a2)*((Tsup+(a21*T1a + a22*T2a)/a2)/2), if(Tsup < T3,(a21*T1a + \
a22*T2a+ a23*T3a)+(aq-a2)*((Tsup+(a21*T1a + \
a22*T2a+a23*T3a)/a2)/2),if(Tsup < T4,(a21*T1a+ a22*T2a+ \
a23*T3a+a24*T4a)+(aq-a2)*((Tsup+(a21*T1a + \
a22*T2a+a23*T3a+a24*T4a)/a2)/2),if(Tsup < T5,(a21*T1a + a22*T2a + \
a23*T3a + a24*T4a+a25*T5a)+(aq-a2)*((Tsup+(a21*T1a + \
a22*T2a+a23*T3a+a24*T4a+a25*T5a)/a2)/2),(a21*T1a + a22*T2a + a23*T3a \
+ a24*T4a +a25*T5a+a26*T6a)+(aq-a2)*((Tsup+(a21*T1a + \
a22*T2a+a23*T3a+a24*T4a+a25*T5a+a26*T6a)/a2)/2))))))) +(1-aq)*Tsup
TurbDampLiq = AD*delta x*0.075*Liquid.density*OmegaDampLiq^2
VELBEGas3 = sqrt((Gas3.Velocity u.Gradient X*ADSX)^2+(Gas3.Velocity \
v.Gradient Y*ADSY)^2+(Gas3.Velocity w.Gradient Z*ADSZ)^2)
VELBELIQUID = sqrt((Liquid.Velocity u.Gradient \
X*ADSX)^2+(Liquid.Velocity v.Gradient Y*ADSY)^2+(Liquid.Velocity \
w.Gradient Z*ADSZ)^2)
VELBET = Liquid.Velocity u.Gradient X + Liquid.Velocity v.Gradient Y + \
Liquid.Velocity w.Gradient Z
VF Gas tot = Gas1.Volume Fraction+Gas3.Volume Fraction
VF Gas3 crit = 0.2
VF Liquid min = 0.1
VIN = 0.2 [m s^-1]
VisGas = 0.12264347e-04 [kg m^-1 s^-1]
VisLiq = 0.28214892e-03 [kg m^-1 s^-1]
WD = Gas3.Wall Distance
WSSGas3 = Gas3.vfc*sqrt(WSSxGas3^2+WSSyGas3^2+WSSzGas3^2)
WSSLIQUID = Liquid.vfc*sqrt(WSSxLiquid^2+WSSyLiquid^2+WSSzLiquid^2)
WSSLiquid = Liquid.vfc*sqrt(WSSxLiquid^2+WSSyLiquid^2+WSSzLiquid^2)
WSSxGas3 = (Gas3.Volume Fraction.Gradient \
X/ADSurf*viscg*(2*Gas3.Velocity u.Gradient X))+(Gas3.Volume \
Fraction.Gradient Y/ADSurf*viscg*(Gas3.Velocity u.Gradient \
Y+Gas3.Velocity v.Gradient X))+(Gas3.Volume Fraction.Gradient \
Z/ADSurf*viscg*(Gas3.Velocity u.Gradient Z+Gas3.Velocity w.Gradient X))
WSSxLiquid = (Liquid.Volume Fraction.Gradient \
X/ADSurf*viscl*(2*Liquid.Velocity u.Gradient X))+(Liquid.Volume \
Fraction.Gradient Y/ADSurf*viscl*(Liquid.Velocity u.Gradient \
Y+Liquid.Velocity v.Gradient X))+(Liquid.Volume Fraction.Gradient \
Z/ADSurf*viscl*(Liquid.Velocity u.Gradient Z+Liquid.Velocity \
w.Gradient X))
WSSyGas3 = (Gas3.Volume Fraction.Gradient X/ADSurf*viscg*(Gas3.Velocity \
v.Gradient X+Gas3.Velocity u.Gradient Y))+(Gas3.Volume \
Fraction.Gradient Y/ADSurf*viscg*(2*Gas3.Velocity v.Gradient \
X))+(Gas3.Volume Fraction.Gradient Z/ADSurf*viscg*(Gas3.Velocity \
v.Gradient Z+Gas3.Velocity w.Gradient Y))
WSSyLiquid = (Liquid.Volume Fraction.Gradient \
X/ADSurf*viscl*(Liquid.Velocity v.Gradient X+Liquid.Velocity \
u.Gradient Y))+(Liquid.Volume Fraction.Gradient \
Y/ADSurf*viscl*(2*Liquid.Velocity v.Gradient X))+(Liquid.Volume \
Fraction.Gradient Z/ADSurf*viscl*(Liquid.Velocity v.Gradient \
Z+Liquid.Velocity w.Gradient Y))
WSSzGas3 = (Gas3.Volume Fraction.Gradient X/ADSurf*viscg*(Gas3.Velocity \
w.Gradient X+Gas3.Velocity u.Gradient Z))+(Gas3.Volume \
Fraction.Gradient Y/ADSurf*viscg*(Gas3.Velocity w.Gradient \
Y+Gas3.Velocity v.Gradient Z))+(Gas3.Volume Fraction.Gradient \
Z/ADSurf*viscg*(2*Gas3.Velocity w.Gradient Z))
WSSzLiquid = (Liquid.Volume Fraction.Gradient \
X/ADSurf*viscl*(Liquid.Velocity w.Gradient X+Liquid.Velocity \
u.Gradient Z))+(Liquid.Volume Fraction.Gradient \
Y/ADSurf*viscl*(Liquid.Velocity w.Gradient Y+Liquid.Velocity \
v.Gradient Z))+(Liquid.Volume Fraction.Gradient \
Z/ADSurf*viscl*(2*Liquid.Velocity w.Gradient Z))
Wall Distance Grad = max(sqrt((Gas3.WDV.Gradient \
X)^2+(Gas3.WDV.Gradient Y)^2+(Gas3.WDV.Gradient Z)^2),minWDGrad)
WeightBubb = 0.5*(1 + cos(pi*(AlphaGasSchlange - (AlphaCritBubb - \
DeltaAlphaBubb))/(2*DeltaAlphaBubb)))
WeightDrop = 0.5*(1 + cos(pi*(AlphaLiqSchlange - (AlphaCritDrop - \
DeltaAlphaDrop))/(2*DeltaAlphaDrop)))
a1 = max(1-aq, 0.0)
a2 = min(1,if(Tsup < T1, a21,if(Tsup < T2,(a21 + a22), if(Tsup < \
T3,(a21 + a22 + a23 ),if(Tsup < T4,(a21 + a22 + a23 +a24),if(Tsup < \
T5,(a21 + a22 + a23 + a24+a25),(a21 + a22 + a23 + a24 +a25+a26)))))))
a21 = min(min(NSD,NSD1)*1/NSD,((pi*(dw1^2)/4*min(NSD,NSD1))))
a21s = min(2,((pi*(dw1s^2)/4*min(NSD,NSD1))))
a22 = \
min((min(NSD,NSD2)-NSD1)*1/NSD,((pi*(dw2^2)/4*(min(NSD,NSD2)-NSD1))))
a22s = min(2-a21s,((pi*(dw2s^2)/4*(min(NSD,NSD2)-NSD1))))
a23 = \
min((min(NSD,NSD3)-NSD2)*1/NSD,((pi*(dw3^2)/4*(min(NSD,NSD3)-NSD2))))
a23s = min(2-a21s-a22s,((pi*(dw3s^2)/4*(min(NSD,NSD3)-NSD2))))
a24 = \
min((min(NSD,NSD4)-NSD3)*1/NSD,((pi*(dw4^2)/4*(min(NSD,NSD4)-NSD3))))
a24s = min(2-a21s-a22s-a23s,((pi*(dw4s^2)/4*(min(NSD,NSD4)-NSD3))))
a25 = \
min((min(NSD,NSD5)-NSD4)*1/NSD,((pi*(dw5^2)/4*(min(NSD,NSD5)-NSD4))))
a25s = min(2-a21s-a22s-a23s-a24s,((pi*(dw5s^2)/4*(min(NSD,NSD5)-NSD4))))
a26 = min((min(NSD,NSD6)-NSD5)/NSD,((pi*(dw6^2)/4*(min(NSD,NSD6)-NSD5))))
a26s = \
min(2-a21s-a22s-a23s-a24s-a25s,((pi*(dw6s^2)/4*(min(NSD,NSD6)-NSD5))))
a2df = if(Tsup < T1, a21*dw1*1/(tg1+tw1),if(Tsup < \
T2,(a21*dw1*1/(tg1+tw1) + a22*dw2*1/(tg2+tw2)), if(Tsup < \
T3,(a21*dw1*1/(tg1+tw1) + a22*dw2*1/(tg2+tw2) + a23*dw3*1/(tg3+tw3) \
),if(Tsup < T4,(a21*dw1*1/(tg1+tw1) + a22 *dw2*1/(tg2+tw2)+ \
a23*dw3*1/(tg3+tw3) +a24*dw4*1/(tg4+tw4)),if(Tsup < \
T5,(a21*dw1*1/(tg1+tw1) + a22*dw2*1/(tg2+tw2) + a23*dw3*1/(tg3+tw3) + \
a24*dw4*1/(tg4+tw4)+a25*dw5*1/(tg5+tw5)),(a21*dw1*1/(tg1+tw1) + \
a22*dw2*1/(tg2+tw2) + a23*dw3*1/(tg3+tw3) + a24 \
*dw4*1/(tg4+tw4)+a25*dw5*1/(tg5+tw5)+a26*dw6*1/(tg6+tw6)))))))
a2e = pi/4*if(Tsup < 0[K],0[m],((if(Tsup < T1, \
dw1^3*min(NSD,NSD1)*f1,if(Tsup < T2,dw1^3*min(NSD,NSD1)*f1 + \
dw2^3*(min(NSD,NSD2)-NSD1)*f2, if(Tsup < \
T3,dw1^3*min(NSD,NSD1)*f1+dw2^3*(NSD2-NSD1)*f2 + \
dw3^3*(min(NSD,NSD3)-NSD2)*f3,if(Tsup < \
T4,dw1^3*min(NSD,NSD1)*f1+dw2^3*(NSD2-NSD1)*f2 +dw3^3*(NSD3-NSD2)*f3+ \
dw4^3*(min(NSD,NSD4)-NSD3)*f4,if(Tsup < \
T5,dw1^3*min(NSD,NSD1)*f1+dw2^3*(NSD2-NSD1)*f2 \
+dw3^3*(NSD3-NSD2)*f3+dw4^3*(NSD4-NSD3)*f4+ \
dw5^3*(min(NSD,NSD5)-NSD4)*f5,dw1^3*min(NSD,NSD1)*f1+dw2^3*(NSD2-NSD1)\
*f2 +dw3^3*(NSD3-NSD2)*f3+dw4^3*(NSD4-NSD3)*f4+dw5^3*(NSD5-NSD4)*f5+ \
dw6^3*(min(NSD,NSD6)-NSD5)*f6))))))/(fr/4*NSD))^(1/3))/dw
a2exp = ((0.0001+(pi*(dw^2)/4*NSD))^-4 +1)^-0.25
absslip1 = max(sqrt(uslip1*uslip1+vslip1*vslip1+wslip1*wslip1),VSlipMin)
absslip2 = ((uslip2*uslip2)+(vslip2*vslip2)+(wslip2*wslip2))^(1./2.)
absslip3 = ((uslip3*uslip3)+(vslip3*vslip3)+(wslip3*wslip3))^(1./2.)
ag = if(aq>= 1 && (tw1 <= mintw1 ||tw2 <= mintw2||tw3 <= mintw3||tw4 <= \
mintw4||tw5 <= mintw5||tw6 <= mintw6) && \
qwall-min(1,a2)*1/6*dw*fr*DenGas*hfg > 0 [W m^-2],1,0)
alphamax = 0.8
alphdm = 1.0
aq = if(Tsup < T1, aq1,if(Tsup < T2,(aq1 + aq2), if(Tsup < T3,(aq1 + \
aq2 + aq3 ),if(Tsup < T4,(aq1 + aq2 + aq3 +aq4),if(Tsup < T5,(aq1 + \
aq2 + aq3 + aq4+ aq5),(aq1 + aq2 + aq3 + aq4 +aq5+aq6))))))
aq1 = min(1,((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,5[K])-1)*dw1))*dw1+pi*(dw1^2)/4)*min(NSD,NSD1)))
aq2 = ((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,15[K])-1)*dw2))*dw2+pi*(dw2^2)/4)*(min(NSD,NSD2)-NSD1))
aq3 = ((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,25[K])-1)*dw3))*dw3s+pi*(dw3^2)/4)*(min(NSD,NSD3)-NSD2))
aq4 = ((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,35[K])-1)*dw4))*dw4s+pi*(dw4^2)/4)*(min(NSD,NSD4)-NSD3))
aq5 = ((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,45[K])-1)*dw5))*dw5s+pi*(dw5^2)/4)*(min(NSD,NSD5)-NSD4))
aq6 = ((min(1/2*((NSD))^(-1/2),max(0 [m],(SlidingL(Ginm, \
TsubL,55[K])-1)*dw6))*dw6s+pi*(dw6^2)/4)*(min(NSD,NSD6)-NSD5))
asign = if(PsiSurf>=0,1,0)
betastar = 0.09
db1 = 0.5 [mm]
db10 = 9.5 [mm]
db2 = 1.5 [mm]
db3 = 2.5 [mm]
db4 = 3.5 [mm]
db5 = 4.5 [mm]
db6 = 5.5 [mm]
db7 = 6.5 [mm]
db8 = 7.5 [mm]
db9 = 8.5 [mm]
dcrit = 4.0*sqrt(SurfTen/(g*(DenLiq-DenGas)))
ddt = 0.0015 [s]
delta x = 0.0005 [m]
dp1 = Gas1.mean particle diameter
dp2 = Gas2.mean particle diameter
dp3 = Gas3.mean particle diameter
dtp = 0.005 [s]
dtpf = 0.25 [s]
dun dn = max(Liquid.Velocity u.Gradient X*inorm x + Liquid.Velocity \
v.Gradient Y*inorm y + Liquid.Velocity w.Gradient Z*inorm z, 0[s^-1])
dw = if(Tsup < 0[K],0[m],((if(Tsup < T1, a21*dw1*f1,if(Tsup < \
T2,dw1*a21*f1 + if(a21+a22<1,dw2*a22*f2,0 [m s^-1]), if(Tsup < \
T3,dw1*a21*f1+if(a21+a22<1,dw2*a22*f2,0[m s^-1]) + \
if(a21+a22+a23<1,dw3*a23*f3,0[m s^-1]),if(Tsup < \
T4,dw1*a21*f1+if(a21+a22<1,dw2*a22*f2,0[m s^-1]) + \
if(a21+a22+a23<1,dw3*a23*f3,0[m s^-1])+ \
if(a21+a22+a23+a24<1,dw4*a24*f4,0[m s^-1]),if(Tsup < \
T5,dw1*a21*f1+if(a21+a22<1,dw2*a22*f2,0[m s^-1]) + \
if(a21+a22+a23<1,dw3*a23*f3,0[m s^-1])+ \
if(a21+a22+a23+a24<1,dw4*a24*f4,0[m s^-1])+ \
if(a21+a22+a23+a24+a25<1,dw5*a25*f5,0[m \
s^-1]),dw1*a21*f1+if(a21+a22<1,dw2*a22*f2,0[m s^-1]) + \
if(a21+a22+a23<1,dw3*a23*f3,0[m s^-1])+ \
if(a21+a22+a23+a24<1,dw4*a24*f4,0[m s^-1])+ \
if(a21+a22+a23+a24+a25<1,dw5*a25*f5,0[m s^-1]) + \
if(a21+a22+a23+a24+a25<1,dw6*a26*f6,0[m s^-1]))))))))/(fr/4*min(1,a2)))
dw1 = dw1s+(dw1b-dw1s)*RSLX1
dw1b = dBubD(Ginm, TsubL, T1a)*max(1,SlidingA(Ginm, TsubL,T1a)*2-1)
dw1s = dBubD(Ginm, TsubL, T1a)
dw2 = dw2s+(dw2b-dw2s)*RSLX2
dw2b = dBubD(Ginm, TsubL, T2a)*max(1,SlidingA(Ginm, TsubL,T2a)*2-1)
dw2s = dBubD(Ginm, TsubL, T2a)
dw3 = dw3s+(dw3b-dw3s)*RSLX3
dw3b = dBubD(Ginm, TsubL, T3a)*max(1,SlidingA(Ginm, TsubL,T3a)*2-1)
dw3s = dBubD(Ginm, TsubL, T3a)
dw4 = dw4s+(dw4b-dw4s)*RSLX4
dw4b = dBubD(Ginm, TsubL, T4a)*max(1,SlidingA(Ginm, TsubL,T4a)*2-1)
dw4s = dBubD(Ginm, TsubL, T4a)
dw5 = dw5s+(dw5b-dw5s)*RSLX5
dw5b = dBubD(Ginm, TsubL, T5a)*max(1,SlidingA(Ginm, TsubL,T5a)*2-1)
dw5s = dBubD(Ginm, TsubL, T5a)
dw6 = dw6s+(dw6b-dw6s)*RSLX6
dw6b = dBubD(Ginm, TsubL, T6a)*max(1,SlidingA(Ginm, TsubL,T6a)*2-1)
dw6s = dBubD(Ginm, TsubL, T6a)
dwlr = Liquid.Velocity w.Gradient X
epsintsource1 = min(uC3*(sqrt(Liquid.ke)/dp1)*kintsource1,epsmax)
epsintsource2 = min(uC3*(sqrt(Liquid.ke)/dp2)*kintsource2,epsmax)
epsmax = 10000[W m^-3 s^-1]
f = max(1/m*(1-VF Gas tot)^(0.5), 10e-07)
f1 = 1/(tw1 +tg1)
f2 = 1/(tw2 +tg2)
f23 = max(0,Gas3.Volume Fraction/VF Gas3 crit-1)
f3 = 1/(tw3 +tg3)
f32 = max(0,1-Gas3.Volume Fraction/VF Gas3 crit)
f4 = 1/(tw4 +tg4)
f5 = 1/(tw5 +tg5)
f6 = 1/(tw6 +tg6)
filmhf = 1.E-8 [m]
filmho = 1.E-4 [m]
fr = 1/(tw +tg)/CoalescCoeffwall
frd = (viscl/viscm)*((1-(Gas3.Volume Fraction))^(1/2))
gasholdup1 = volumeAve(Gas1.vf)@tube
gasholdup2 = volumeAve(Gas2.vf)@tube
gasholdup3 = volumeAve(Gas3.vf)@tube
gasout1 = areaAve(Gas1.vf)@Out
gasout2 = areaAve(Gas2.vf)@Out
gasout3 = areaAve(Gas3.vf)@Out
gn = g*inorm x
hfg = EnthGas-EnthLiq
hin = 0 [m]
hlg = EnthGas - EnthLiq
hout = 0.5 [m]
hquenching = 0 [W m^-2. K^-1]
htube = hout-hin
hwtotal = Qwtotal/Tsup
inorm x = Liquid.Volume Fraction.Gradient X/max(sqrt(Liquid.Volume \
Fraction.Gradient X^2+Liquid.Volume Fraction.Gradient \
Y^2+Liquid.Volume Fraction.Gradient Z^2),0.00005[m^-1])
inorm y = Liquid.Volume Fraction.Gradient Y/max(sqrt(Liquid.Volume \
Fraction.Gradient X^2+Liquid.Volume Fraction.Gradient \
Y^2+Liquid.Volume Fraction.Gradient Z^2),0.00005[m^-1])
inorm z = Liquid.Volume Fraction.Gradient Z/max(sqrt(Liquid.Volume \
Fraction.Gradient X^2+Liquid.Volume Fraction.Gradient \
Y^2+Liquid.Volume Fraction.Gradient Z^2),0.00005[m^-1])
jg = 0.0 [m s^-1]
jl = 0.0 [m s^-1]
kappa = -1*(Gas3.nclwXV.Gradient X+Gas3.nclwYV.Gradient \
Y+Gas3.nclwZV.Gradient Z)
kintsource1 = min(0.75*(CD1/dp1)*Gas1.vf*DenLiq*absslip1^3, kmax)
kintsource2 = min(0.75*(CD2/dp2)*Gas2.vf*DenLiq*absslip2^3, kmax)
kmax = 10000 [W m^-3]
lbound = sqrt(abs(7.7e-4*gn*delta x + 0.22*SurfTen/DenLiq/delta x))
minVF = 1e-7
minVFArea = 1.0E-7
minVFGrad = 0.00005[m^-1]
minWDGrad = 0.0000005
mintw1 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T1a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mintw2 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T2a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mintw3 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T3a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mintw4 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T4a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mintw5 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T5a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mintw6 = \
(pi*CondLiq/(2*qwall))^2*max(0[K],T6a+TsubLc)^2*(DenLiq*CpGas)/(pi*Con\
dLiq)
mix visc = max((1-VF Gas tot)^(-2.5*m), 10e-07)
n = 6
ncWX = (Gas3.WDV.Gradient X)/(Wall Distance Grad)
ncWY = (Gas3.WDV.Gradient Y)/(Wall Distance Grad)
ncWZ = (Gas3.WDV.Gradient Z)/(Wall Distance Grad)
nclX = (Gas3.Volume Fraction.Gradient X)/(Gas3 VF Grad)
nclY = (Gas3.Volume Fraction.Gradient Y)/(Gas3 VF Grad)
nclZ = (Gas3.Volume Fraction.Gradient Z)/(Gas3 VF Grad)
nclwX = if((Gas3.Wall Distance <= 1*(delta \
x)),((ncWX*(cos((((ContactAngle-180)/360)*2*pi)[rad])))+(nclX*(sin((((\
ContactAngle-180)/360)*2*pi)[rad])))),nclX)
nclwY = if((Gas3.Wall Distance <= 1*(delta \
x)),((ncWY*(cos((((ContactAngle-180)/360)*2*pi)[rad])))+(nclY*(sin((((\
ContactAngle-180)/360)*2*pi)[rad])))),nclY)
nclwZ = if((Gas3.Wall Distance <= 1*(delta \
x)),((ncWZ*(cos((((ContactAngle-180)/360)*2*pi)[rad])))+(nclZ*(sin((((\
ContactAngle-180)/360)*2*pi)[rad])))),nclZ)
ngabs = max(sqrt((Gas3.Volume Fraction.Gradient X)^2 + (Gas3.Volume \
Fraction.Gradient Y)^2 + (Gas3.Volume Fraction.Gradient Z)^2), \
1.0E-3[m^-1])
ngx = -Gas3.Volume Fraction.Gradient X/ngabs
ngy = -Gas3.Volume Fraction.Gradient Y/ngabs
ngz = -Gas3.Volume Fraction.Gradient Z/ngabs
nlabs = max(sqrt((Liquid.Volume Fraction.Gradient X)^2 + (Liquid.Volume \
Fraction.Gradient Y)^2 + (Liquid.Volume Fraction.Gradient Z)^2), \
1.0E-3[m^-1])
nlx = -Liquid.Volume Fraction.Gradient X/nlabs
nly = -Liquid.Volume Fraction.Gradient Y/nlabs
nlz = -Liquid.Volume Fraction.Gradient Z/nlabs
qwall = 100E+3 [W m^-2]
radius = if(De>=DP,((DP-(2e-3[m]))/2),De/2)
rate13 = DenGas*Gas1.Volume Fraction*(1/tau23)
rate23 = DenGas*Gas2.Volume Fraction*(1/tau23)
rdm = 1
rezBub = Gas1.Volume Fraction/Gas1.Mean Particle Diameter+Gas2.Volume \
Fraction/Gas2.Mean Particle Diameter
rpipe = 0.0125 [m]
tau = 1e-4 [s]
tau23 = 1e-4 [s]
tefintsource1 = (1.0/betastar/Liquid.ke*epsintsource1 - \
Liquid.tef/Liquid.ke*kintsource1)
tefintsource2 = (1.0/betastar/Liquid.ke*epsintsource2 - \
Liquid.tef/Liquid.ke*kintsource2)
tende = 5.0 [s]
tg = if(Tsup < T1, tg1,if(Tsup < T2,(tg1+tg2)/2, if(Tsup < \
T3,(tg1+tg2+tg3)/3,if(Tsup < T4,(tg1+tg2+tg3+tg4)/4,if(Tsup < \
T5,(tg1+tg2+tg3+tg4+tg5)/5,(tg1+tg2+tg3+tg4+tg5+tg6)/6)))))
tg1 = GTime(Ginm, TsubL, T1a) + (LTime(Ginm, TsubL, T1a) - GTime(Ginm, \
TsubL, T1a))* RSLX1
tg2 = GTime(Ginm, TsubL, T2a)+ (LTime(Ginm, TsubL, T2a)- GTime(Ginm, \
TsubL, T2a))* RSLX2
tg3 = GTime(Ginm, TsubL, T3a) + (LTime(Ginm, TsubL, T3a)- GTime(Ginm, \
TsubL, T3a))* RSLX3
tg4 = GTime(Ginm, TsubL, T4a) + (LTime(Ginm, TsubL, T4a)- GTime(Ginm, \
TsubL, T4a))* RSLX4
tg5 = GTime(Ginm, TsubL, T5a) + (LTime(Ginm, TsubL, T5a)- GTime(Ginm, \
TsubL, T5a))* RSLX5
tg6 = GTime(Ginm, TsubL, T6a) + (LTime(Ginm, TsubL, T6a)- GTime(Ginm, \
TsubL, T6a))* RSLX6
tgfs = sqrt(tgfsx^2 + tgfsy^2 + tgfsz^2)
tgfsx = tgx*(1-ngx^2) - tgy*ngx*ngy - tgz*ngx*ngz
tgfsy = -tgx*ngy*ngx + tgy*(1-ngy^2) - tgz*ngy*ngz
tgfsz = -tgx*ngz*ngx - tgy*ngz*ngy + tgz*(1-ngz^2)
tgx = viscg*(2*Gas3.u.Gradient X*ngx + (Gas3.u.Gradient Y + \
Gas3.v.Gradient X)*ngy + (Gas3.u.Gradient Z + Gas3.w.Gradient X)*ngz)
tgy = viscg*((Gas3.v.Gradient X + Gas3.u.Gradient Y)*ngx + \
2*Gas3.v.Gradient Y*ngy + (Gas3.v.Gradient Z + Gas3.w.Gradient Y)*ngz)
tgz = viscg*((Gas3.w.Gradient X + Gas3.u.Gradient Z)*ngx + \
(Gas3.w.Gradient Y + Gas3.v.Gradient Z)*ngy + 2*Gas3.w.Gradient Z*ngz)
tlfs = sqrt(tlfsx^2 + tlfsy^2 + tlfsz^2)
tlfsx = tlx*(1-nlx^2) - tly*nlx*nly - tlz*nlx*nlz
tlfsy = -tlx*nly*nlx + tly*(1-nly^2) - tlz*nly*nlz
tlfsz = -tlx*nlz*nlx - tly*nlz*nly + tlz*(1-nlz^2)
tlx = Liquid.viscosity*(2*Liquid.u.Gradient X*nlx + (Liquid.u.Gradient \
Y + Liquid.v.Gradient X)*nly + (Liquid.u.Gradient Z + \
Liquid.w.Gradient X)*nlz)
tly = Liquid.viscosity*((Liquid.v.Gradient X + Liquid.u.Gradient Y)*nlx \
+ 2*Liquid.v.Gradient Y*nly + (Liquid.v.Gradient Z + \
Liquid.w.Gradient Y)*nlz)
tlz = Liquid.viscosity*((Liquid.w.Gradient X + Liquid.u.Gradient Z)*nlx \
+ (Liquid.w.Gradient Y + Liquid.v.Gradient Z)*nly + \
2*Liquid.w.Gradient Z*nlz)
tw = if(a2 <1, if(Tsup < T1, 4*tg1,if(Tsup < T2,2*(tg1 +tg2), if(Tsup < \
T3,4/3*(tg1 +tg2 + tg3),if(Tsup < T4,4/4*(tg1 +tg2 + tg3 + \
tg4),if(Tsup < T5,4/5*(tg1 +tg2 + tg3 + tg4+tg5),4/6*(tg1 +tg2 + tg3 \
+ tg4+tg5 +tg6)))))),if(Tsup < T1, mintw1,if(Tsup < T2,(mintw1 + \
mintw2)/2, if(Tsup < T3,(mintw1+mintw2+mintw3)/3,if(Tsup < \
T4,(mintw1+mintw2+mintw3+mintw4)/4,if(Tsup < \
T5,(mintw1+mintw2+mintw3+mintw4+mintw5)/5,(mintw1+mintw2+mintw3+mintw6\
)/6))))))
tw1 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T1a),if(qwall*a21- \
a21*2/3*dw1*1/(tg1+WTime(Ginm, TsubL, T1a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T1a),(2/3*dw1*DenGas*hfg)/qwall-tg1)),max(mintw1,if(qwall*a21- \
a21*2/3*dw1*1/(tg1+mintw1)*DenGas*hfg > 0 [W \
m^-2],mintw1,(2/3*dw1*DenGas*hfg)/qwall-tg1)))
tw2 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T2a),if(qwall*a22- \
a22*2/3*dw2*1/(tg2+WTime(Ginm, TsubL, T2a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T2a),(2/3*dw2*DenGas*hfg)/qwall-tg2)),max(mintw2,if(qwall*a22- \
a22*2/3*dw2*1/(tg2+mintw2)*DenGas*hfg > 0 [W \
m^-2],mintw2,(2/3*dw2*DenGas*hfg)/qwall-tg2)))
tw3 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T3a),if(qwall*a23- \
a23*2/3*dw3*1/(tg3+WTime(Ginm, TsubL, T3a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T3a),(2/3*dw3*DenGas*hfg)/qwall-tg3)),max(mintw3,if(qwall*a23- \
a23*2/3*dw3*1/(tg3+mintw3)*DenGas*hfg > 0 [W \
m^-2],mintw3,(2/3*dw3*DenGas*hfg)/qwall-tg3)))
tw4 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T4a),if(qwall*a24- \
a24*2/3*dw4*1/(tg4+WTime(Ginm, TsubL, T4a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T4a),(2/3*dw4*DenGas*hfg)/qwall-tg4)),max(mintw4,if(qwall*a24- \
a24*2/3*dw4*1/(tg4+mintw4)*DenGas*hfg > 0 [W \
m^-2],mintw4,(2/3*dw4*DenGas*hfg)/qwall-tg4)))
tw5 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T5a),if(qwall*a25- \
a25*2/3*dw5*1/(tg5+WTime(Ginm, TsubL, T5a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T5a),(2/3*dw5*DenGas*hfg)/qwall-tg5)),max(mintw5,if(qwall*a25- \
a25*2/3*dw5*1/(tg5+mintw5)*DenGas*hfg > 0 [W \
m^-2],mintw5,(2/3*dw5*DenGas*hfg)/qwall-tg5)))
tw6 = if(TsubLc>0[K],max(WTime(Ginm, TsubL, T6a),if(qwall*a26- \
a26*2/3*dw6*1/(tg6+WTime(Ginm, TsubL, T6a))*DenGas*hfg > 0 [W \
m^-2],WTime(Ginm, TsubL, \
T6a),(2/3*dw6*DenGas*hfg)/qwall-tg6)),max(mintw6,if(qwall*a26- \
a26*2/3*dw6*1/(tg6+mintw6)*DenGas*hfg > 0 [W \
m^-2],mintw6,(2/3*dw6*DenGas*hfg)/qwall-tg6)))
uC3 = 1.0
ubound = sqrt(abs(0.13*gn*delta x + 1.57*SurfTen/DenLiq/delta x))
uslip = sqrt((Liquid.u-Gas3.u)^2+(Liquid.v-Gas3.v)^2+(Liquid.w-Gas3.w)^2)
uslip2 = Liquid.u - Gas1.u
uslip3 = Liquid.u - Gas3.u
vfair = if(De<=DP,(0.5*tanh(((radius)-sqrt(((x)^2 + (y)^2 + \
(z-(radius+(3e-3[m])))^2)))/(0.0012[m]))+0.5),(0.5*tanh(((1)-((x/(radi\
us))^2 + (y/(radius))^2 + \
(((z)-(zmax+(3e-3[m])))/(zmax))^2))/(1e-3))+0.5))
vfwater = 1-vfair
viscg = 1.34e-5 [kg m^-1 s^-1]
viscl = 0.0009 [kg m^-1 s^-1]
viscm = ((1-((Gas3.Volume \
Fraction)/rdm))^((-2.5)*(rdm)*((viscg+(0.4*viscl))/(viscg+viscl))))*(v\
iscl)
vofg01 = 0.9895979772497
vofg02 = 0.0104020227502996
vofg1 = vofg01 * vofg
vofg2 = vofg02 * vofg
vslip2 = Liquid.v - Gas1.v
vslip3 = Liquid.v - Gas3.v
wslip2 = Liquid.w - Gas1.w
wslip3 = Liquid.w - Gas3.w
zmax = (((De/2)^3)/(radius^2))
END
FUNCTION: GTime
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = [s]
INTERPOLATION DATA:
Data = \
0.1,0,5,0.00437351,0.1,0,6,0.004094463,0.1,0,7,0.003716653,0.1,0,8,0\
.003277104,0.1,0,9,0.002993625,0.1,0,10,0.002692336,0.1,0,11,0.00240\
3399,0.1,0,12,0.002208736,0.1,0,13,0.001912386,0.1,0,14,0.001746907,\
0.1,0.75,5,0.00437351,0.1,0.75,6,0.004096932,0.1,0.75,7,0.003719122,\
0.1,0.75,8,0.003430193,0.1,0.75,9,0.003060292,0.1,0.75,10,0.00269480\
5,0.1,0.75,11,0.002403399,0.1,0.75,12,0.002159353,0.1,0.75,13,0.0019\
12386,0.1,0.75,14,0.001746907,0.1,1.5,5,0.004380917,0.1,1.5,6,0.0041\
0187,0.1,1.5,7,0.003719122,0.1,1.5,8,0.003341303,0.1,1.5,9,0.0029936\
25,0.1,1.5,10,0.002694805,0.1,1.5,11,0.002403399,0.1,1.5,12,0.002211\
205,0.1,1.5,13,0.001914855,0.1,1.5,14,0.001746907,0.1,2.25,5,0.00437\
8448,0.1,2.25,6,0.00410187,0.1,2.25,7,0.003719122,0.1,2.25,8,0.00334\
3772,0.1,2.25,9,0.003050416,0.1,2.25,10,0.002697274,0.1,2.25,11,0.00\
2403399,0.1,2.25,12,0.002213675,0.1,2.25,13,0.001919794,0.1,2.25,14,\
0.001746907,0.1,3,5,0.004378448,0.1,3,6,0.004096932,0.1,3,7,0.003719\
122,0.1,3,8,0.003346241,0.1,3,9,0.002996094,0.1,3,10,0.002697274,0.1\
,3,11,0.002403399,0.1,3,12,0.002213675,0.1,3,13,0.001919794,0.1,3,14\
,0.001749376,0.2,0,5,0.001590801,0.2,0,6,0.001440124,0.2,0,7,0.00126\
7237,0.2,0,8,0.001143742,0.2,0,9,0.001038043,0.2,0,10,0.0009515723,0\
.2,0,11,0.0008577026,0.2,0,12,0.0007692157,0.2,0,13,0.0006975615,0.2\
,0,14,0.0006357892,0.2,0.75,5,0.001590801,0.2,0.75,6,0.001440124,0.2\
,0.75,7,0.001267237,0.2,0.75,8,0.001143742,0.2,0.75,9,0.001038043,0.\
2,0.75,10,0.0009515723,0.2,0.75,11,0.0008601717,0.2,0.75,12,0.000769\
2157,0.2,0.75,13,0.0006975615,0.2,0.75,14,0.0006357892,0.2,1.5,5,0.0\
01590801,0.2,1.5,6,0.001442593,0.2,1.5,7,0.001267237,0.2,1.5,8,0.001\
143742,0.2,1.5,9,0.001038043,0.2,1.5,10,0.0009515723,0.2,1.5,11,0.00\
08601717,0.2,1.5,12,0.0007346481,0.2,1.5,13,0.0006975615,0.2,1.5,14,\
0.0006357892,0.2,2.25,5,0.001590801,0.2,2.25,6,0.001442593,0.2,2.25,\
7,0.001267237,0.2,2.25,8,0.001143742,0.2,2.25,9,0.001038043,0.2,2.25\
,10,0.0009515723,0.2,2.25,11,0.0008601717,0.2,2.25,12,0.0007346481,0\
.2,2.25,13,0.0006975615,0.2,2.25,14,0.0006382583,0.2,3,5,0.00159327,\
0.2,3,6,0.001442593,0.2,3,7,0.001267237,0.2,3,8,0.001143742,0.2,3,9,\
0.001040512,0.2,3,10,0.0009515723,0.2,3,11,0.0008601717,0.2,3,12,0.0\
007346481,0.2,3,13,0.0006975615
Option = Three Dimensional
END
END
FUNCTION: LTime
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = [s]
INTERPOLATION DATA:
Data = \
0.1,0,5,0.00437351,0.1,0,6,0.004094463,0.1,0,7,0.003716653,0.1,0,8,0\
.003277104,0.1,0,9,0.002993625,0.1,0,10,0.002692336,0.1,0,11,0.00240\
3399,0.1,0,12,0.002208736,0.1,0,13,0.001912386,0.1,0,14,0.001746907,\
0.1,0.75,5,0.00437351,0.1,0.75,6,0.004096932,0.1,0.75,7,0.003719122,\
0.1,0.75,8,0.003430193,0.1,0.75,9,0.003060292,0.1,0.75,10,0.00269480\
5,0.1,0.75,11,0.002403399,0.1,0.75,12,0.002159353,0.1,0.75,13,0.0019\
12386,0.1,0.75,14,0.001746907,0.1,1.5,5,0.004380917,0.1,1.5,6,0.0041\
0187,0.1,1.5,7,0.003719122,0.1,1.5,8,0.003341303,0.1,1.5,9,0.0029936\
25,0.1,1.5,10,0.002694805,0.1,1.5,11,0.002403399,0.1,1.5,12,0.002211\
205,0.1,1.5,13,0.001914855,0.1,1.5,14,0.001746907,0.1,2.25,5,0.00437\
8448,0.1,2.25,6,0.00410187,0.1,2.25,7,0.003719122,0.1,2.25,8,0.00334\
3772,0.1,2.25,9,0.003050416,0.1,2.25,10,0.002697274,0.1,2.25,11,0.00\
2403399,0.1,2.25,12,0.002213675,0.1,2.25,13,0.001919794,0.1,2.25,14,\
0.001746907,0.1,3,5,0.004378448,0.1,3,6,0.004096932,0.1,3,7,0.003719\
122,0.1,3,8,0.003346241,0.1,3,9,0.002996094,0.1,3,10,0.002697274,0.1\
,3,11,0.002403399,0.1,3,12,0.002213675,0.1,3,13,0.001919794,0.1,3,14\
,0.001749376,0.2,0,5,0.001590801,0.2,0,6,0.001440124,0.2,0,7,0.00126\
7237,0.2,0,8,0.001143742,0.2,0,9,0.001038043,0.2,0,10,0.0009515723,0\
.2,0,11,0.0008577026,0.2,0,12,0.0007692157,0.2,0,13,0.0006975615,0.2\
,0,14,0.0006357892,0.2,0.75,5,0.001590801,0.2,0.75,6,0.001440124,0.2\
,0.75,7,0.001267237,0.2,0.75,8,0.001143742,0.2,0.75,9,0.001038043,0.\
2,0.75,10,0.0009515723,0.2,0.75,11,0.0008601717,0.2,0.75,12,0.000769\
2157,0.2,0.75,13,0.0006975615,0.2,0.75,14,0.0006357892,0.2,1.5,5,0.0\
01590801,0.2,1.5,6,0.001442593,0.2,1.5,7,0.001267237,0.2,1.5,8,0.001\
143742,0.2,1.5,9,0.001038043,0.2,1.5,10,0.0009515723,0.2,1.5,11,0.00\
08601717,0.2,1.5,12,0.0007346481,0.2,1.5,13,0.0006975615,0.2,1.5,14,\
0.0006357892,0.2,2.25,5,0.001590801,0.2,2.25,6,0.001442593,0.2,2.25,\
7,0.001267237,0.2,2.25,8,0.001143742,0.2,2.25,9,0.001038043,0.2,2.25\
,10,0.0009515723,0.2,2.25,11,0.0008601717,0.2,2.25,12,0.0007346481,0\
.2,2.25,13,0.0006975615,0.2,2.25,14,0.0006382583,0.2,3,5,0.00159327,\
0.2,3,6,0.001442593,0.2,3,7,0.001267237,0.2,3,8,0.001143742,0.2,3,9,\
0.001040512,0.2,3,10,0.0009515723,0.2,3,11,0.0008601717,0.2,3,12,0.0\
007346481,0.2,3,13,0.0006975615
Option = Three Dimensional
END
END
FUNCTION: SlidingA
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = []
INTERPOLATION DATA:
Data = \
0.1,0,5,1,0.1,0,6,1,0.1,0,7,1,0.1,0,8,1,0.1,0,9,1,0.1,0,10,1,0.1,0,1\
1,1,0.1,0,12,1,0.1,0,13,1,0.1,0,14,1,0.1,0.75,5,1,0.1,0.75,6,1,0.1,0\
.75,7,1,0.1,0.75,8,1,0.1,0.75,9,1,0.1,0.75,10,1,0.1,0.75,11,1,0.1,0.\
75,12,1,0.1,0.75,13,1,0.1,0.75,14,1,0.1,1.5,5,1,0.1,1.5,6,1,0.1,1.5,\
7,1,0.1,1.5,8,1,0.1,1.5,9,1,0.1,1.5,10,1,0.1,1.5,11,1,0.1,1.5,12,1,0\
.1,1.5,13,1,0.1,1.5,14,1,0.1,2.25,5,1,0.1,2.25,6,1,0.1,2.25,7,1,0.1,\
2.25,8,1,0.1,2.25,9,1,0.1,2.25,10,1,0.1,2.25,11,1,0.1,2.25,12,1,0.1,\
2.25,13,1,0.1,2.25,14,1,0.1,3,5,1,0.1,3,6,1,0.1,3,7,1,0.1,3,8,1,0.1,\
3,9,1,0.1,3,10,1,0.1,3,11,1,0.1,3,12,1,0.1,3,13,1,0.1,3,14,1,0.2,0,5\
,1,0.2,0,6,1,0.2,0,7,1,0.2,0,8,1,0.2,0,9,1,0.2,0,10,1,0.2,0,11,1,0.2\
,0,12,1,0.2,0,13,1,0.2,0,14,1,0.2,0.75,5,1,0.2,0.75,6,1,0.2,0.75,7,1\
,0.2,0.75,8,1,0.2,0.75,9,1,0.2,0.75,10,1,0.2,0.75,11,1,0.2,0.75,12,1\
,0.2,0.75,13,1,0.2,0.75,14,1,0.2,1.5,5,1,0.2,1.5,6,1,0.2,1.5,7,1,0.2\
,1.5,8,1,0.2,1.5,9,1,0.2,1.5,10,1,0.2,1.5,11,1,0.2,1.5,12,1,0.2,1.5,\
13,1,0.2,1.5,14,1,0.2,2.25,5,1,0.2,2.25,6,1,0.2,2.25,7,1,0.2,2.25,8,\
1,0.2,2.25,9,1,0.2,2.25,10,1,0.2,2.25,11,1,0.2,2.25,12,1,0.2,2.25,13\
,1,0.2,2.25,14,1,0.2,3,5,1,0.2,3,6,1,0.2,3,7,1,0.2,3,8,1,0.2,3,9,1,0\
.2,3,10,1,0.2,3,11,1,0.2,3,12,1,0.2,3,13,1
Option = Three Dimensional
END
END
FUNCTION: SlidingL
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = []
INTERPOLATION DATA:
Data = \
0.1,0,5,1,0.1,0,6,1,0.1,0,7,1,0.1,0,8,1,0.1,0,9,1,0.1,0,10,1,0.1,0,1\
1,1,0.1,0,12,1,0.1,0,13,1,0.1,0,14,1,0.1,0.75,5,1,0.1,0.75,6,1,0.1,0\
.75,7,1,0.1,0.75,8,1,0.1,0.75,9,1,0.1,0.75,10,1,0.1,0.75,11,1,0.1,0.\
75,12,1,0.1,0.75,13,1,0.1,0.75,14,1,0.1,1.5,5,1,0.1,1.5,6,1,0.1,1.5,\
7,1,0.1,1.5,8,1,0.1,1.5,9,1,0.1,1.5,10,1,0.1,1.5,11,1,0.1,1.5,12,1,0\
.1,1.5,13,1,0.1,1.5,14,1,0.1,2.25,5,1,0.1,2.25,6,1,0.1,2.25,7,1,0.1,\
2.25,8,1,0.1,2.25,9,1,0.1,2.25,10,1,0.1,2.25,11,1,0.1,2.25,12,1,0.1,\
2.25,13,1,0.1,2.25,14,1,0.1,3,5,1,0.1,3,6,1,0.1,3,7,1,0.1,3,8,1,0.1,\
3,9,1,0.1,3,10,1,0.1,3,11,1,0.1,3,12,1,0.1,3,13,1,0.1,3,14,1,0.2,0,5\
,1,0.2,0,6,1,0.2,0,7,1,0.2,0,8,1,0.2,0,9,1,0.2,0,10,1,0.2,0,11,1,0.2\
,0,12,1,0.2,0,13,1,0.2,0,14,1,0.2,0.75,5,1,0.2,0.75,6,1,0.2,0.75,7,1\
,0.2,0.75,8,1,0.2,0.75,9,1,0.2,0.75,10,1,0.2,0.75,11,1,0.2,0.75,12,1\
,0.2,0.75,13,1,0.2,0.75,14,1,0.2,1.5,5,1,0.2,1.5,6,1,0.2,1.5,7,1,0.2\
,1.5,8,1,0.2,1.5,9,1,0.2,1.5,10,1,0.2,1.5,11,1,0.2,1.5,12,1,0.2,1.5,\
13,1,0.2,1.5,14,1,0.2,2.25,5,1,0.2,2.25,6,1,0.2,2.25,7,1,0.2,2.25,8,\
1,0.2,2.25,9,1,0.2,2.25,10,1,0.2,2.25,11,1,0.2,2.25,12,1,0.2,2.25,13\
,1,0.2,2.25,14,1,0.2,3,5,1,0.2,3,6,1,0.2,3,7,1,0.2,3,8,1,0.2,3,9,1,0\
.2,3,10,1,0.2,3,11,1,0.2,3,12,1,0.2,3,13,1
Option = Three Dimensional
END
END
FUNCTION: WTime
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = [s]
INTERPOLATION DATA:
Data = \
0.1,0,5,0.03262673,0.1,0,6,0.03627384,0.1,0,7,0.03943728,0.1,0,8,0.0\
42438,0.1,0,9,0.04632295,0.1,0,10,0.05006,0.1,0,11,0.05422448,0.1,0,\
12,0.05938181,0.1,0,13,0.06387521,0.1,0,14,0.06996145,0.1,0.75,5,0.0\
3493002,0.1,0.75,6,0.03897178,0.1,0.75,7,0.04265578,0.1,0.75,8,0.046\
91781,0.1,0.75,9,0.05038999,0.1,0.75,10,0.05424761,0.1,0.75,11,0.058\
75734,0.1,0.75,12,0.06385005,0.1,0.75,13,0.06913562,0.1,0.75,14,0.07\
558935,0.1,1.5,5,0.03754041,0.1,1.5,6,0.04196446,0.1,1.5,7,0.0460282\
8,0.1,1.5,8,0.04984545,0.1,1.5,9,0.05425054,0.1,1.5,10,0.05870774,0.\
1,1.5,11,0.06357507,0.1,1.5,12,0.06967971,0.1,1.5,13,0.07469572,0.1,\
1.5,14,0.08150459,0.1,2.25,5,0.04041097,0.1,2.25,6,0.04523954,0.1,2.\
25,7,0.04966837,0.1,2.25,8,0.05405055,0.1,2.25,9,0.05883794,0.1,2.25\
,10,0.06346015,0.1,2.25,11,0.06868015,0.1,2.25,12,0.07516462,0.1,2.2\
5,13,0.08075058,0.1,2.25,14,0.08770223,0.1,3,5,0.04358862,0.1,3,6,0.\
04879702,0.1,3,7,0.05366,0.1,3,8,0.05835287,0.1,3,9,0.06329292,0.1,3\
,10,0.06848508,0.1,3,11,0.07406762,0.1,3,12,0.08092698,0.1,3,13,0.08\
686313,0.1,3,14,0.09419955,0.2,0,5,0.02305883,0.2,0,6,0.02598002,0.2\
,0,7,0.02877802,0.2,0,8,0.03212689,0.2,0,9,0.03626369,0.2,0,10,0.040\
48716,0.2,0,11,0.04520719,0.2,0,12,0.05022872,0.2,0,13,0.05565792,0.\
2,0,14,0.06192935,0.2,0.75,5,0.02530286,0.2,0.75,6,0.02856931,0.2,0.\
75,7,0.03175207,0.2,0.75,8,0.03547584,0.2,0.75,9,0.03998999,0.2,0.75\
,10,0.04458094,0.2,0.75,11,0.04969561,0.2,0.75,12,0.0550624,0.2,0.75\
,13,0.0608492,0.2,0.75,14,0.06747824,0.2,1.5,5,0.0278046,0.2,1.5,6,0\
.03146569,0.2,1.5,7,0.03501347,0.2,1.5,8,0.03911212,0.2,1.5,9,0.0439\
9869,0.2,1.5,10,0.04895959,0.2,1.5,11,0.05444421,0.2,1.5,12,0.059390\
82,0.2,1.5,13,0.06635252,0.2,1.5,14,0.0733441,0.2,2.25,5,0.03058875,\
0.2,2.25,6,0.03462967,0.2,2.25,7,0.03855726,0.2,2.25,8,0.04303574,0.\
2,2.25,9,0.04828979,0.2,2.25,10,0.05362311,0.2,2.25,11,0.05947768,0.\
2,2.25,12,0.06478684,0.2,2.25,13,0.07212835,0.2,2.25,14,0.07948741,0\
.2,3,5,0.03367011,0.2,3,6,0.03808838,0.2,3,7,0.04238592,0.2,3,8,0.04\
723929,0.2,3,9,0.05289292,0.2,3,10,0.0585715,0.2,3,11,0.06479602,0.2\
,3,12,0.07047018,0.2,3,13,0.07818659
Option = Three Dimensional
END
END
FUNCTION: dBubD
Argument Units = [m s^-1],[K],[K]
Option = Interpolation
Result Units = [m]
INTERPOLATION DATA:
Data = \
0.1,0,5,0.001860352,0.1,0,6,0.001969655,0.1,0,7,0.002032295,0.1,0,8,\
0.002055065,0.1,0,9,0.002106164,0.1,0,10,0.002126805,0.1,0,11,0.0021\
40436,0.1,0,12,0.002176565,0.1,0,13,0.002151925,0.1,0,14,0.002171197\
,0.1,0.75,5,0.001859473,0.1,0.75,6,0.001969623,0.1,0.75,7,0.00203248\
5,0.1,0.75,8,0.0020937,0.1,0.75,9,0.002125637,0.1,0.75,10,0.00212726\
4,0.1,0.75,11,0.002140163,0.1,0.75,12,0.002157002,0.1,0.75,13,0.0021\
51702,0.1,0.75,14,0.002171041,0.1,1.5,5,0.001860643,0.1,1.5,6,0.0019\
7014,0.1,1.5,7,0.002032108,0.1,1.5,8,0.002071845,0.1,1.5,9,0.0021053\
77,0.1,1.5,10,0.002126886,0.1,1.5,11,0.002139871,0.1,1.5,12,0.002176\
999,0.1,1.5,13,0.002152195,0.1,1.5,14,0.002170905,0.1,2.25,5,0.00185\
936,0.1,2.25,6,0.001969632,0.1,2.25,7,0.002031592,0.1,2.25,8,0.00207\
2055,0.1,2.25,9,0.002122278,0.1,2.25,10,0.002127318,0.1,2.25,11,0.00\
2139647,0.1,2.25,12,0.002177889,0.1,2.25,13,0.002154101,0.1,2.25,14,\
0.002170726,0.1,3,5,0.001858754,0.1,3,6,0.001967662,0.1,3,7,0.002030\
859,0.1,3,8,0.00207189,0.1,3,9,0.002105498,0.1,3,10,0.002126873,0.1,\
3,11,0.002139255,0.1,3,12,0.002177806,0.1,3,13,0.002153895,0.1,3,14,\
0.002171692,0.2,0,5,0.001092886,0.2,0,6,0.001140039,0.2,0,7,0.001166\
593,0.2,0,8,0.001189794,0.2,0,9,0.001219732,0.2,0,10,0.001242678,0.2\
,0,11,0.001254581,0.2,0,12,0.001267706,0.2,0,13,0.001276965,0.2,0,14\
,0.001285804,0.2,0.75,5,0.001092634,0.2,0.75,6,0.001139861,0.2,0.75,\
7,0.001166421,0.2,0.75,8,0.001189758,0.2,0.75,9,0.001219753,0.2,0.75\
,10,0.0012426,0.2,0.75,11,0.001255825,0.2,0.75,12,0.001267605,0.2,0.\
75,13,0.0012769,0.2,0.75,14,0.00128575,0.2,1.5,5,0.001092379,0.2,1.5\
,6,0.001140547,0.2,1.5,7,0.001166275,0.2,1.5,8,0.001189638,0.2,1.5,9\
,0.001219636,0.2,1.5,10,0.001242509,0.2,1.5,11,0.001255724,0.2,1.5,1\
2,0.001244078,0.2,1.5,13,0.001276904,0.2,1.5,14,0.00128629,0.2,2.25,\
5,0.001092186,0.2,2.25,6,0.001140271,0.2,2.25,7,0.001166113,0.2,2.25\
,8,0.001189477,0.2,2.25,9,0.001219515,0.2,2.25,10,0.001242325,0.2,2.\
25,11,0.00125563,0.2,2.25,12,0.001244013,0.2,2.25,13,0.001276762,0.2\
,2.25,14,0.001288102,0.2,3,5,0.001092796,0.2,3,6,0.001140095,0.2,3,7\
,0.001165934,0.2,3,8,0.001189361,0.2,3,9,0.001220448,0.2,3,10,0.0012\
42236,0.2,3,11,0.001255437,0.2,3,12,0.001243901,0.2,3,13,0.001276687
Option = Three Dimensional
END
END
END
ADDITIONAL VARIABLE: AD V
Option = Definition
Tensor Type = SCALAR
Units = [m^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: ADBubb V
Option = Definition
Tensor Type = SCALAR
Units = [m^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: ADSurf V
Option = Definition
Tensor Type = SCALAR
Units = [m^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Bubdia V
Option = Definition
Tensor Type = SCALAR
Units = [m ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: CD V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: CDSurf V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: CF V
Option = Definition
Tensor Type = VECTOR
Units = [N m^-3]
Variable Type = Specific
END
ADDITIONAL VARIABLE: CLIFT
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Ctd V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: GasDTot V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: GasTot V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Growtime
Option = Definition
Tensor Type = SCALAR
Units = [s]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Lifttime
Option = Definition
Tensor Type = SCALAR
Units = [s]
Variable Type = Specific
END
ADDITIONAL VARIABLE: MEnt
Option = Definition
Tensor Type = SCALAR
Units = [kg m^-3 s^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Phi clust V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Phi coal V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Phi morph V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Phi surf V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: QwLiquidV
Option = Definition
Tensor Type = SCALAR
Units = [W m^-2 ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: QwtotalV
Option = Definition
Tensor Type = SCALAR
Units = [W m^-2 ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: QwvaporV
Option = Definition
Tensor Type = SCALAR
Units = [W m^-2 ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: ST V
Option = Definition
Tensor Type = VECTOR
Units = [N m^-3]
Variable Type = Specific
END
ADDITIONAL VARIABLE: SlidingA
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: SlidingL
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: TsupLiquid
Option = Definition
Tensor Type = SCALAR
Units = [K ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: WDV
Option = Definition
Tensor Type = SCALAR
Units = [m]
Variable Type = Specific
END
ADDITIONAL VARIABLE: Waitingtime
Option = Definition
Tensor Type = SCALAR
Units = [s ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: a2V
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: aqV
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: dBubD
Option = Definition
Tensor Type = SCALAR
Units = [m]
Variable Type = Specific
END
ADDITIONAL VARIABLE: dBubD2
Option = Definition
Tensor Type = SCALAR
Units = [m ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: hwalltotalV
Option = Definition
Tensor Type = SCALAR
Units = [W m^-2 K^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: nclwXV
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: nclwYV
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: nclwZV
Option = Definition
Tensor Type = SCALAR
Units = [ ]
Variable Type = Specific
END
ADDITIONAL VARIABLE: uslip V
Option = Definition
Tensor Type = SCALAR
Units = [m s^-1]
Variable Type = Specific
END
ADDITIONAL VARIABLE: waitingtime2
Option = Definition
Tensor Type = SCALAR
Units = [s]
Variable Type = Specific
END
MATERIAL: DefLiquid
Material Group = User
Option = Pure Substance
Thermodynamic State = Liquid
PROPERTIES:
Option = General Material
EQUATION OF STATE:
Density = DenLiq
Molar Mass = 18.02 [kg kmol^-1]
Option = Value
END
SPECIFIC HEAT CAPACITY:
Option = Value
Specific Heat Capacity = CpLiq
Specific Heat Type = Constant Pressure
END
REFERENCE STATE:
Option = Specified Point
Reference Pressure = 1 [atm]
Reference Specific Enthalpy = EnthLiq
Reference Temperature = Tsat
END
DYNAMIC VISCOSITY:
Dynamic Viscosity = VisLiq
Option = Value
END
THERMAL CONDUCTIVITY:
Option = Value
Thermal Conductivity = CondLiq
END
THERMAL EXPANSIVITY:
Option = Value
Thermal Expansivity = 1/Tsat
END
END
END
MATERIAL: DefVapour
Material Group = User
Option = Pure Substance
Thermodynamic State = Gas
PROPERTIES:
Option = General Material
EQUATION OF STATE:
Density = DenGas
Molar Mass = 18 [kg kmol^-1]
Option = Value
END
SPECIFIC HEAT CAPACITY:
Option = Value
Specific Heat Capacity = CpGas
Specific Heat Type = Constant Pressure
END
REFERENCE STATE:
Option = Specified Point
Reference Pressure = 1 [atm]
Reference Specific Enthalpy = EnthGas
Reference Temperature = Tsat
END
DYNAMIC VISCOSITY:
Dynamic Viscosity = VisGas
Option = Value
END
THERMAL CONDUCTIVITY:
Option = Value
Thermal Conductivity = CondGas
END
THERMAL EXPANSIVITY:
Option = Value
Thermal Expansivity = 1/Tsat
END
END
END
END
FLOW: Flow Analysis 1
SOLUTION UNITS:
Angle Units = [rad]
Length Units = [m]
Mass Units = [kg]
Solid Angle Units = [sr]
Temperature Units = [K]
Time Units = [s]
END
ANALYSIS TYPE:
Option = Transient
EXTERNAL SOLVER COUPLING:
Option = None
END
INITIAL TIME:
Option = Automatic with Value
Time = 0 [s]
END
TIME DURATION:
Option = Total Time
Total Time = 7 [s]
END
TIME STEPS:
Option = Timesteps
Timesteps = 0.0015 [s]
END
END
DOMAIN: tube
Coord Frame = Coord 0
Domain Type = Fluid
Location = FLUID
BOUNDARY: HeatSurf
Boundary Type = WALL
Location = MANTEL
BOUNDARY CONDITIONS:
HEAT TRANSFER:
Fixed Temperature = TWall
Option = Fixed Temperature
END
MASS AND MOMENTUM:
Option = Fluid Dependent
END
WALL CONTACT MODEL:
Option = Specify Area Fraction
END
WALL ROUGHNESS:
Option = Smooth Wall
END
END
FLUID PAIR: Gas1 | Liquid
BOUNDARY CONDITIONS:
WALL BOILING MODEL:
Option = RPI Model
END
END
END
FLUID: Gas1
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL CONTACT AREA:
Area Fraction = 0
Option = Area Fraction
END
END
END
FLUID: Gas2
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL CONTACT AREA:
Area Fraction = 0
Option = Area Fraction
END
END
END
FLUID: Gas3
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL CONTACT AREA:
Area Fraction = 0
Option = Area Fraction
END
END
END
FLUID: Liquid
BOUNDARY CONDITIONS:
MASS AND MOMENTUM:
Option = No Slip Wall
END
WALL CONTACT AREA:
Area Fraction = 1
Option = Area Fraction
END
END
END
END
BOUNDARY: IN
Boundary Type = INLET
Location = INLET
BOUNDARY CONDITIONS:
FLOW REGIME:
Option = Subsonic
END
HEAT TRANSFER:
Option = Static Temperature
Static Temperature = TIN
END
MASS AND MOMENTUM:
Normal Speed = VIN
Option = Normal Speed
END
TURBULENCE:
Option = Medium Intensity and Eddy Viscosity Ratio
END
END
FLUID: Gas1
BOUNDARY CONDITIONS:
SIZE GROUP: Group 1
Option = Automatic
END
SIZE GROUP: Group 2
Option = Automatic
END
SIZE GROUP: Group 3
Option = Automatic
END
SIZE GROUP: Group 4
Option = Automatic
END
SIZE GROUP: Group 5
Option = Automatic
END
SIZE GROUP: Group 6
Option = Automatic
END
VOLUME FRACTION:
Option = Value
Volume Fraction = 0.0
END
END
END
FLUID: Gas2
BOUNDARY CONDITIONS:
SIZE GROUP: Group 7
Option = Automatic
END
SIZE GROUP: Group 8
Option = Automatic
END
SIZE GROUP: Group 9
Option = Automatic
END
VOLUME FRACTION:
Option = Value
Volume Fraction = 0.0
END
END
END
FLUID: Gas3
BOUNDARY CONDITIONS:
SIZE GROUP: Group 10
Option = Automatic
END
VOLUME FRACTION:
Option = Value
Volume Fraction = if(x^2+y^2 < (0.01 [m])^2,0.99,0.01)
END
END
END
FLUID: Liquid
BOUNDARY CONDITIONS:
VOLUME FRACTION:
Option = Value
Volume Fraction = if(x^2+y^2 < (0.01 [m])^2,0.01,0.99)
END
END
END
END
BOUNDARY: Out
Boundary Type = OUTLET
Location = OUTLET
BOUNDARY CONDITIONS:
FLOW REGIME:
Option = Subsonic
END
MASS AND MOMENTUM:
Implicit Pressure Averaging = Off
Option = Average Static Pressure
Pressure Profile Blend = 0.05
Relative Pressure = 0 [Pa]
END
PRESSURE AVERAGING:
Option = Average Over Whole Outlet
END
END
END
DOMAIN MODELS:
BUOYANCY MODEL:
Buoyancy Reference Density = DenLiq
Gravity X Component = 0 [m s^-2]
Gravity Y Component = 0 [m s^-2]
Gravity Z Component = -g
Option = Buoyant
BUOYANCY REFERENCE LOCATION:
Option = Automatic
END
END
DOMAIN MOTION:
Option = Stationary
END
MESH DEFORMATION:
Option = None
END
REFERENCE PRESSURE:
Reference Pressure = PRef1
END
END
FLUID DEFINITION: Gas1
Material = DefVapour
Option = Material Library
MORPHOLOGY:
Minimum Volume Fraction = MinVolFrac
Option = Polydispersed Fluid
END
END
FLUID DEFINITION: Gas2
Material = DefVapour
Option = Material Library
MORPHOLOGY:
Minimum Volume Fraction = MinVolFrac
Option = Polydispersed Fluid
END
END
FLUID DEFINITION: Gas3
Material = DefVapour
Option = Material Library
MORPHOLOGY:
Minimum Volume Fraction = MinVolFrac
Option = Polydispersed Fluid
END
END
FLUID DEFINITION: Liquid
Material = DefLiquid
Option = Material Library
MORPHOLOGY:
Minimum Volume Fraction = MinVolFrac
Option = Continuous Fluid
END
END
FLUID MODELS:
ADDITIONAL VARIABLE: AD V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: ADBubb V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: ADSurf V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Bubdia V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: CD V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: CDSurf V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: CF V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: CLIFT
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: GasDTot V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: GasTot V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Growtime
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Lifttime
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: MEnt
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Phi clust V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Phi coal V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Phi morph V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Phi surf V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: QwLiquidV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: QwtotalV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: QwvaporV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: ST V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: SlidingA
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: SlidingL
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: TsupLiquid
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: WDV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: Waitingtime
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: a2V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: aqV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: dBubD
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: dBubD2
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: hwalltotalV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: nclwXV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: nclwYV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: nclwZV
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: uslip V
Option = Fluid Dependent
END
ADDITIONAL VARIABLE: waitingtime2
Option = Fluid Dependent
END
COMBUSTION MODEL:
Option = None
END
FLUID: Gas1
ADDITIONAL VARIABLE: Growtime
Additional Variable Value = GTime(Ginm, TsubL, T1a)
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Lifttime
Additional Variable Value = LTime(Ginm, TsubL, T1a)
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: MEnt
Additional Variable Value = MEntrain
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: QwLiquidV
Additional Variable Value = Qwliquid
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: QwtotalV
Additional Variable Value = Qwtotal
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: QwvaporV
Additional Variable Value = Qwvapor
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: SlidingA
Additional Variable Value = SlidingA(Ginm,TsubL,T3a)
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: SlidingL
Additional Variable Value = SlidingL(Ginm,TsubL,T3a)
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Waitingtime
Additional Variable Value = tw
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: a2V
Additional Variable Value = a2
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: aqV
Additional Variable Value = aq
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: dBubD
Additional Variable Value = dBubD(Ginm, TsubL, Tsup)
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: dBubD2
Additional Variable Value = dw
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: hwalltotalV
Additional Variable Value = hwtotal
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: waitingtime2
Additional Variable Value = WTime(Ginm, TsubL, T1a)
Option = Algebraic Equation
END
FLUID BUOYANCY MODEL:
Option = Density Difference
END
HEAT TRANSFER MODEL:
Fluid Temperature = Tsat
Option = Isothermal
END
TURBULENCE MODEL:
Option = Dispersed Phase Zero Equation
END
END
FLUID: Gas2
FLUID BUOYANCY MODEL:
Option = Density Difference
END
HEAT TRANSFER MODEL:
Fluid Temperature = Tsat
Option = Isothermal
END
TURBULENCE MODEL:
Option = Dispersed Phase Zero Equation
END
END
FLUID: Gas3
ADDITIONAL VARIABLE: AD V
Additional Variable Value = AD
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: ADBubb V
Additional Variable Value = ADBubb
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: ADSurf V
Additional Variable Value = ADSurf
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: CD V
Additional Variable Value = CD
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: CDSurf V
Additional Variable Value = CDSurf
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: CF V
Option = Vector Algebraic Equation
Vector xValue = CFX*(-1)
Vector yValue = CFY*(-1)
Vector zValue = CFZ*(-1)
END
ADDITIONAL VARIABLE: CLIFT
Additional Variable Value = CL
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Phi clust V
Additional Variable Value = Phi clust
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Phi coal V
Additional Variable Value = Phi coal
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Phi morph V
Additional Variable Value = Phi morph
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Phi surf V
Additional Variable Value = Phi surf
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: ST V
Option = Vector Algebraic Equation
Vector xValue = STX*(Gas3.Volume Fraction)
Vector yValue = STY*(Gas3.Volume Fraction)
Vector zValue = STZ*(Gas3.Volume Fraction)
END
ADDITIONAL VARIABLE: WDV
Additional Variable Value = WD
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: nclwXV
Additional Variable Value = nclX
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: nclwYV
Additional Variable Value = nclY
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: nclwZV
Additional Variable Value = nclZ
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: uslip V
Additional Variable Value = uslip
Option = Algebraic Equation
END
FLUID BUOYANCY MODEL:
Option = Density Difference
END
HEAT TRANSFER MODEL:
Fluid Temperature = Tsat
Option = Isothermal
END
TURBULENCE MODEL:
Option = Dispersed Phase Zero Equation
END
END
FLUID: Liquid
ADDITIONAL VARIABLE: ADBubb V
Additional Variable Value = ADBubb
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: Bubdia V
Additional Variable Value = GasDTot/rezBub
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: GasDTot V
Additional Variable Value = GasDTot
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: GasTot V
Additional Variable Value = GasTot
Option = Algebraic Equation
END
ADDITIONAL VARIABLE: TsupLiquid
Additional Variable Value = TLiqsup
Option = Algebraic Equation
END
FLUID BUOYANCY MODEL:
Option = Density Difference
END
HEAT TRANSFER MODEL:
Include Viscous Work Term = Off
Option = Total Energy
END
TURBULENCE MODEL:
Option = SST
BUOYANCY TURBULENCE:
Option = None
END
END
TURBULENT WALL FUNCTIONS:
High Speed Model = Off
Option = Automatic
END
END
HEAT TRANSFER MODEL:
Homogeneous Model = False
Option = Fluid Dependent
END
THERMAL RADIATION MODEL:
Option = None
END
TURBULENCE MODEL:
Homogeneous Model = False
Option = Fluid Dependent
END
END
FLUID PAIR: Gas1 | Liquid
Surface Tension Coefficient = SurfTen
INTERPHASE HEAT TRANSFER:
Option = Two Resistance
FLUID1 INTERPHASE HEAT TRANSFER:
Option = Zero Resistance
END
FLUID2 INTERPHASE HEAT TRANSFER:
Option = Ranz Marshall
END
END
INTERPHASE TRANSFER MODEL:
Minimum Volume Fraction for Area Density = MinVFArea
Option = Particle Model
END
MASS TRANSFER:
Option = Phase Change
PHASE CHANGE MODEL:
Heat Transfer Coefficient Under Relaxation Factor = 0.1
Option = Thermal Phase Change
Saturation Temperature = Tsat
WALL BOILING MODEL:
Bubble Diameter Influence Factor = 2.0
Fixed Yplus for Liquid Subcooling = 250.0
Mass Source Under Relaxation = 0.1
Maximum Area Fraction of Bubble Influence = 0.5
Onset of Boiling Superheating = 0.05 [K]
Option = RPI Model
Wall Heat Flux Partitioning Under Relaxation = 0.1
Wall Superheating Under Relaxation = 0.1
BUBBLE DEPARTURE DIAMETER:
Bubble Departure Diameter = dw
Option = User Defined
END
BUBBLE DETACHMENT FREQUENCY:
Bubble Detachment Frequency = fr
Option = User Defined
END
BUBBLE WAITING TIME:
Bubble Waiting Time = tw
Option = User Defined
END
LIQUID QUENCHING HEAT TRANSFER COEFFICIENT:
Option = User Defined
Quenching Heat Transfer Coefficient = hquenching
END
PARTITIONING AREA FRACTIONS:
Evaporative Area Fraction = min(1,a2)*(1-Gas3.Volume Fraction)
Gas Convective Area Fraction = 0
Liquid Convective Area Fraction = \
(1-min(0.999,max(0.001,aq*2)))*(1-Gas3.Volume Fraction)
Option = User Defined
Quenching Area Fraction = min(1,max(aq,0.001))*(1-Gas3.Volume \
Fraction)
END
USER WALL HEAT PARTITIONING TERM:
Option = To Liquid
User Partitioning Term = \
(qwall*min(0.999,aq*2)-min(1,a2)*1/6*dw*fr*DenGas*hfg)*(1-Gas3\
.Volume Fraction)
END
WALL NUCLEATION SITE DENSITY:
Option = User Defined
Wall Nucleation Site Density = NSD
END
END
END
END
MOMENTUM TRANSFER:
DRAG FORCE:
Option = Ishii Zuber
END
LIFT FORCE:
Option = Tomiyama
END
TURBULENT DISPERSION FORCE:
Option = Favre Averaged Drag Force
Turbulent Dispersion Coefficient = 1.0
END
VIRTUAL MASS FORCE:
Option = None
END
WALL LUBRICATION FORCE:
Lubrication Coefficient C1 = -0.01
Lubrication Coefficient C2 = 0.05
Option = Antal
END
END
TURBULENCE TRANSFER:
ENHANCED TURBULENCE PRODUCTION MODEL:
Option = Sato Enhanced Eddy Viscosity
END
END
END
FLUID PAIR: Gas2 | Liquid
Surface Tension Coefficient = SurfTen
INTERPHASE HEAT TRANSFER:
Option = Two Resistance
FLUID1 INTERPHASE HEAT TRANSFER:
Option = Zero Resistance
END
FLUID2 INTERPHASE HEAT TRANSFER:
Option = Ranz Marshall
END
END
INTERPHASE TRANSFER MODEL:
Minimum Volume Fraction for Area Density = MinVFArea
Option = Particle Model
END
MASS TRANSFER:
Option = Phase Change
PHASE CHANGE MODEL:
Heat Transfer Coefficient Under Relaxation Factor = 0.1
Option = Thermal Phase Change
Saturation Temperature = Tsat
END
END
MOMENTUM TRANSFER:
DRAG FORCE:
Option = Ishii Zuber
END
LIFT FORCE:
Option = Tomiyama
END
TURBULENT DISPERSION FORCE:
Option = Favre Averaged Drag Force
Turbulent Dispersion Coefficient = 1.0
END
VIRTUAL MASS FORCE:
Option = None
END
WALL LUBRICATION FORCE:
Lubrication Coefficient C1 = -0.01
Lubrication Coefficient C2 = 0.05
Option = Antal
END
END
TURBULENCE TRANSFER:
ENHANCED TURBULENCE PRODUCTION MODEL:
Option = Sato Enhanced Eddy Viscosity
END
END
END
FLUID PAIR: Gas3 | Liquid
Surface Tension Coefficient = SurfTen
INTERPHASE HEAT TRANSFER:
Option = Two Resistance
FLUID1 INTERPHASE HEAT TRANSFER:
Option = Zero Resistance
END
FLUID2 INTERPHASE HEAT TRANSFER:
Option = Hughmark
END
END
INTERPHASE TRANSFER MODEL:
Interfacial Area Density = AD
Minimum Volume Fraction for Area Density = MinVFArea
Option = Particle Model
END
MASS TRANSFER:
Option = Phase Change
PHASE CHANGE MODEL:
Heat Transfer Coefficient Under Relaxation Factor = 0.1
Option = Thermal Phase Change
Saturation Temperature = Tsat
WALL BOILING MODEL:
Bubble Diameter Influence Factor = 2.0
Fixed Yplus for Liquid Subcooling = 250.0
Mass Source Under Relaxation = 0.1
Maximum Area Fraction of Bubble Influence = 0.5
Onset of Boiling Superheating = 0.05 [K]
Option = RPI Model
Wall Heat Flux Partitioning Under Relaxation = 0.1
Wall Superheating Under Relaxation = 0.1
BUBBLE DEPARTURE DIAMETER:
Bubble Departure Diameter = dw
Option = User Defined
END
BUBBLE DETACHMENT FREQUENCY:
Bubble Detachment Frequency = fr
Option = User Defined
END
BUBBLE WAITING TIME:
Bubble Waiting Time = tw
Option = User Defined
END
LIQUID QUENCHING HEAT TRANSFER COEFFICIENT:
Option = User Defined
Quenching Heat Transfer Coefficient = 0 [W m^-2 K^-1]
END
PARTITIONING AREA FRACTIONS:
Evaporative Area Fraction = 0
Gas Convective Area Fraction = Gas3.Volume Fraction
Liquid Convective Area Fraction = 0
Option = User Defined
Quenching Area Fraction = 0
END
USER WALL HEAT PARTITIONING TERM:
Option = To Liquid
User Partitioning Term = 0 [W m^-2]
END
WALL NUCLEATION SITE DENSITY:
Option = User Defined
Wall Nucleation Site Density = NSD
END
END
END
END
MOMENTUM TRANSFER:
DRAG FORCE:
Drag Coefficient = CD
Option = Drag Coefficient
END
LIFT FORCE:
Lift Coefficient = CL
Option = Lift Coefficient
END
TURBULENT DISPERSION FORCE:
Option = None
END
VIRTUAL MASS FORCE:
Option = None
END
WALL LUBRICATION FORCE:
Option = None
END
END
TURBULENCE TRANSFER:
ENHANCED TURBULENCE PRODUCTION MODEL:
Option = None
END
END
END
INITIALISATION:
Option = Automatic
FLUID: Gas1
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = VIN
END
SIZE GROUP: Group 1
Option = Automatic
END
SIZE GROUP: Group 2
Option = Automatic
END
SIZE GROUP: Group 3
Option = Automatic
END
SIZE GROUP: Group 4
Option = Automatic
END
SIZE GROUP: Group 5
Option = Automatic
END
SIZE GROUP: Group 6
Option = Automatic
END
VOLUME FRACTION:
Option = Automatic with Value
Volume Fraction = 0.01
END
END
END
FLUID: Gas2
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = VIN
END
SIZE GROUP: Group 7
Option = Automatic
END
SIZE GROUP: Group 8
Option = Automatic
END
SIZE GROUP: Group 9
Option = Automatic
END
VOLUME FRACTION:
Option = Automatic with Value
Volume Fraction = 0.01
END
END
END
FLUID: Gas3
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = VIN
END
SIZE GROUP: Group 10
Option = Automatic
END
VOLUME FRACTION:
Option = Automatic with Value
Volume Fraction = if(x^2+y^2 < (0.01 [m])^2,0.97,0.01)
END
END
END
FLUID: Liquid
INITIAL CONDITIONS:
Velocity Type = Cartesian
CARTESIAN VELOCITY COMPONENTS:
Option = Automatic with Value
U = 0 [m s^-1]
V = 0 [m s^-1]
W = VIN
END
TEMPERATURE:
Option = Automatic with Value
Temperature = TIN
END
TURBULENCE INITIAL CONDITIONS:
Option = Medium Intensity and Eddy Viscosity Ratio
END
VOLUME FRACTION:
Option = Automatic with Value
Volume Fraction = if(x^2+y^2 < (0.01 [m])^2,0.01,0.97)
END
END
END
INITIAL CONDITIONS:
STATIC PRESSURE:
Option = Automatic with Value
Relative Pressure = PHydro
END
END
END
MULTIPHASE MODELS:
Homogeneous Model = False
FREE SURFACE MODEL:
Option = None
END
POLYDISPERSED FLUID: Gas
Option = Inhomogeneous MUSIG
Polydispersed Fluids List = Gas1,Gas2,Gas3
BREAKUP MODEL:
Breakup Coefficient = 0.01
Option = Luo and Svendsen
END
COALESCENCE MODEL:
Option = Prince and Blanch
Turbulence Coalescence Coefficient = 2
END
SIZE GROUP DISTRIBUTION:
Number of Size Groups = 10
Option = User Defined
SIZE GROUP: Group 1
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db1
END
SIZE GROUP: Group 10
Option = Definition
Polydispersed Fluid = Gas3
Size Group Diameter = db10
END
SIZE GROUP: Group 2
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db2
END
SIZE GROUP: Group 3
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db3
END
SIZE GROUP: Group 4
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db4
END
SIZE GROUP: Group 5
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db5
END
SIZE GROUP: Group 6
Option = Definition
Polydispersed Fluid = Gas1
Size Group Diameter = db6
END
SIZE GROUP: Group 7
Option = Definition
Polydispersed Fluid = Gas2
Size Group Diameter = db7
END
SIZE GROUP: Group 8
Option = Definition
Polydispersed Fluid = Gas2
Size Group Diameter = db8
END
SIZE GROUP: Group 9
Option = Definition
Polydispersed Fluid = Gas2
Size Group Diameter = db9
END
END
END
END
SUBDOMAIN: CF
Coord Frame = Coord 0
Location = FLUID
FLUID: Gas3
SOURCES:
MOMENTUM SOURCE:
GENERAL MOMENTUM SOURCE:
Momentum Source X Component = CFX*(-1)
Momentum Source Y Component = CFY*(-1)
Momentum Source Z Component = CFZ*(-1)
Option = Cartesian Components
END
END
END
END
FLUID: Liquid
SOURCES:
MOMENTUM SOURCE:
GENERAL MOMENTUM SOURCE:
Momentum Source X Component = CFX
Momentum Source Y Component = CFY
Momentum Source Z Component = CFZ
Option = Cartesian Components
END
END
END
END
END
SUBDOMAIN: Complete Coalescence
Coord Frame = Coord 0
Location = FLUID
FLUID: Gas1
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = -rate13*Phi coal
VARIABLE: Group 1.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 2.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 3.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 4.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 5.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 6.sf
Option = Value
Value = 1 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
FLUID: Gas2
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = -rate23*Phi coal
VARIABLE: Group 7.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 8.sf
Option = Value
Value = 1 []
END
VARIABLE: Group 9.sf
Option = Value
Value = 1 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
FLUID: Gas3
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = (rate13+rate23)*Phi coal
VARIABLE: Group 10.sf
Option = Value
Value = 1 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
END
SUBDOMAIN: Entrainment
Coord Frame = Coord 0
Location = FLUID
FLUID: Gas1
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = MEntrain*vofg01
VARIABLE: Group 1.sf
Option = Value
Value = 0.16 []
END
VARIABLE: Group 2.sf
Option = Value
Value = 0.16 []
END
VARIABLE: Group 3.sf
Option = Value
Value = 0.16 []
END
VARIABLE: Group 4.sf
Option = Value
Value = 0.16 []
END
VARIABLE: Group 5.sf
Option = Value
Value = 0.16 []
END
VARIABLE: Group 6.sf
Option = Value
Value = 0.16 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
FLUID: Gas2
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = MEntrain*vofg02
VARIABLE: Group 7.sf
Option = Value
Value = 0.33 []
END
VARIABLE: Group 8.sf
Option = Value
Value = 0.33 []
END
VARIABLE: Group 9.sf
Option = Value
Value = 0.33 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
FLUID: Gas3
SOURCES:
EQUATION SOURCE: continuity
Option = Fluid Mass Source
Source = -MEntrain
VARIABLE: Group 10.sf
Option = Value
Value = 1 []
END
VARIABLE: vel
Option = Cartesian Vector Components
xValue = 0 [m s^-1]
yValue = 0 [m s^-1]
zValue = 0 [m s^-1]
END
END
END
END
END
SUBDOMAIN: GasC_Turbulence
Coord Frame = Coord 0
Location = FLUID
FLUID: Liquid
SOURCES:
EQUATION SOURCE: ke
Option = Source
Source = SmWaTurb
END
EQUATION SOURCE: tef
Option = Source
Source = TurbDampLiq
END
END
END
END
SUBDOMAIN: Surface Tension
Coord Frame = Coord 0
Location = FLUID
FLUID: Gas3
SOURCES:
MOMENTUM SOURCE:
GENERAL MOMENTUM SOURCE:
Momentum Source X Component = STX*(Gas3.Volume Fraction)
Momentum Source Y Component = STY*(Gas3.Volume Fraction)
Momentum Source Z Component = STZ*(Gas3.Volume Fraction)
Option = Cartesian Components
END
END
END
END
FLUID: Liquid
SOURCES:
MOMENTUM SOURCE:
GENERAL MOMENTUM SOURCE:
Momentum Source X Component = STX*(Liquid.Volume Fraction)
Momentum Source Y Component = STY*(Liquid.Volume Fraction)
Momentum Source Z Component = STZ*(Liquid.Volume Fraction)
Option = Cartesian Components
END
END
END
END
END
END
OUTPUT CONTROL:
BACKUP DATA RETENTION:
Option = Delete Old Files
END
BACKUP RESULTS: Backup Results 1
File Compression Level = Default
Option = Standard
OUTPUT FREQUENCY:
Option = Timestep Interval
Timestep Interval = 10
END
END
MONITOR OBJECTS:
Monitor Coefficient Loop Convergence = On
MONITOR BALANCES:
Option = Full
END
MONITOR FORCES:
Option = Full
END
MONITOR PARTICLES:
Option = Full
END
MONITOR POINT: GasHoldup1
Coord Frame = Coord 0
Expression Value = gasholdup1
Option = Expression
END
MONITOR POINT: GasHoldup2
Coord Frame = Coord 0
Expression Value = gasholdup2
Option = Expression
END
MONITOR POINT: GasHoldup3
Coord Frame = Coord 0
Expression Value = gasholdup3
Option = Expression
END
MONITOR POINT: GasHoldupDis
Coord Frame = Coord 0
Expression Value = gasholdup1+gasholdup2
Option = Expression
END
MONITOR POINT: GasHoldupTot
Coord Frame = Coord 0
Expression Value = gasholdup1+gasholdup2+gasholdup3
Option = Expression
END
MONITOR POINT: GasOut1
Coord Frame = Coord 0
Expression Value = gasout1
Option = Expression
END
MONITOR POINT: GasOut2
Coord Frame = Coord 0
Expression Value = gasout2
Option = Expression
END
MONITOR POINT: GasOut3
Coord Frame = Coord 0
Expression Value = gasout3
Option = Expression
END
MONITOR POINT: LiquidTempAvg
Coord Frame = Coord 0
Expression Value = TempLiquidDomain
Option = Expression
END
MONITOR POINT: LiquidTempOUT
Coord Frame = Coord 0
Expression Value = TempLiquidOut
Option = Expression
END
MONITOR POINT: Liquidholdup
Coord Frame = Coord 0
Expression Value = LiquidHoldup
Option = Expression
END
MONITOR POINT: MaxCourant
Coord Frame = Coord 0
Expression Value = MaxCN
Option = Expression
END
MONITOR POINT: Q_Liquid
Coord Frame = Coord 0
Expression Value = QavgLiquid
Option = Expression
END
MONITOR POINT: Q_Vapour
Coord Frame = Coord 0
Expression Value = QavgVapour
Option = Expression
END
MONITOR POINT: Tw_point
Cartesian Coordinates = 0.0125 [m], 0 [m], 0.45 [m]
Coord Frame = Coord 0
Option = Cartesian Coordinates
Output Variables List = Gas1.Waitingtime
MONITOR LOCATION CONTROL:
Interpolation Type = Nearest Vertex
END
POSITION UPDATE FREQUENCY:
Option = Initial Mesh Only
END
END
MONITOR POINT: dw_Point
Cartesian Coordinates = 0.0125 [m], 0 [m], 0.45 [m]
Coord Frame = Coord 0
Option = Cartesian Coordinates
Output Variables List = Gas1.dBubD
MONITOR LOCATION CONTROL:
Interpolation Type = Nearest Vertex
END
POSITION UPDATE FREQUENCY:
Option = Initial Mesh Only
END
END
MONITOR RESIDUALS:
Option = Full
END
MONITOR TOTALS:
Option = Full
END
END
RESULTS:
File Compression Level = Default
Option = Standard
Output Equation Residuals = All
END
TRANSIENT RESULTS: Transient Results 1
File Compression Level = Default
Include Mesh = On
Option = Selected Variables
Output Variables List = Gas1.Volume Fraction,Gas1.Velocity,Gas2.Volume \
Fraction,Gas2.Velocity,Gas3.Volume \
Fraction,Gas3.Velocity,Liquid.Temperature,Courant \
Number,Gas1.hwalltotalV,Gas1.QwLiquidV,Gas1.QwtotalV,Gas1.QwvaporV
OUTPUT FREQUENCY:
Option = Timestep Interval
Timestep Interval = 2
END
END
TRANSIENT RESULTS: Transient Results 2
File Compression Level = Default
Option = Standard
Output Equation Residuals = All
OUTPUT FREQUENCY:
Option = Time Interval
Time Interval = dtpf
END
END
TRANSIENT STATISTICS: Transient Statistics 1
Option = Arithmetic Average
Output Variables List = Gas1.Volume Fraction,Gas2.Volume \
Fraction,Gas3.Volume Fraction,GasDTot V,GasTot \
V,Liquid.Temperature,Liquid.Volume \
Fraction,TsupLiquid,Gas1.QwLiquidV,Gas1.QwtotalV,Gas1.QwvaporV,Gas1.hw\
alltotalV
END
END
SOLVER CONTROL:
Turbulence Numerics = First Order
ADVECTION SCHEME:
Option = Upwind
END
CONVERGENCE CONTROL:
Maximum Number of Coefficient Loops = 40
Minimum Number of Coefficient Loops = 4
Timescale Control = Coefficient Loops
END
CONVERGENCE CRITERIA:
Residual Target = 1.E-4
Residual Type = RMS
END
MULTIPHASE CONTROL:
Volume Fraction Coupling = Segregated
END
TRANSIENT SCHEME:
Option = First Order Backward Euler
END
END
EXPERT PARAMETERS:
build artificial wall = f
relax mass = 0.1
END
END
COMMAND FILE:
Version = 18.2
END