N_n "(210*(max(T_w, T_sat)-T_sat)*1[1/K])^1.805*1[1/m^2]" "Nucleation site density"
q_q A_q*h_q*(T_w-T_const) "Thermal flux from quenching"
q_e m_e*Delta_Hlg "Thermal flux from evaporation"
V_bw pi*D_bw^3/6 "Bubble departure volume"
A_b min(K*pi*N_n*D_bw^2/4,1) "Fraction of wall area influence by bubbles"
f_n sqrt(4*g_const*(phtr.rho_s1-phtr.rho_s2)/(3*phtr.rho_s1*D_bw)) "Bubble release frequency"
m_e A_e*N_n*f_n*V_bw*phtr.rho_s2 "Vapor mass flux at wall"
h_q 2/sqrt(pi)*f_n*sqrt(t_q*nitmm1.kc_iso*phtr.rho_s1*nitmm1.Cp_s1) "Heat transfer coefficient for quenching"
t_q 0.8/f_n "Waiting time elapsed between the departure of the bubble and the nucleation of a new one"
f_alpha if(alpha_l>=alpha_crit,1-0.5*exp(-20*(alpha_l-alpha_crit)),0.5*(alpha_l/alpha_crit)^(20*alpha_crit)) "Function for RPI extension"
A_cpl (1-A_b)*f_alpha "Weight factor for convection with pure liquid"
A_q A_b*f_alpha "Weight factor for quenching"
A_e f_alpha "Weight factor for evaporation"
A_cpg (1-f_alpha) "Weight factor for convection with pure gas"
Yp_f spf.d_w_plus "y+ where properties are evaluated for convective heat flux"
Tp_f nitf1.Tplus_d "T+ near the wall used for evaluation of convective heat flux"
T_f down(T) "Temperature of liquid near the wall"
T_w up(T) "Temperature of the wall"
Tp_const nif(Yp_const<10/down(nitf1.Pr)^(1/3),nojac(down(nitf1.Pr)*Yp_const),Yp_const<sqrt(1000*spf.kappav/down(nitf1.Prt)),nojac(15*down(nitf1.Pr)^(2/3)-500/Yp_const^2),nojac(down(nitf1.Prt)*log(Yp_const)/spf.kappav+15*down(nitf1.Pr)^(2/3)-0.5*down(nitf1.Prt)*(1+log(1000*spf.kappav/down(nitf1.Prt)))/spf.kappav)) "T+ at y+=Yp_const"
T_const T_w-Tp_const/Tp_f*(T_w-T_f) "Temperature of the liquid at y+=Yp_const"
Y_f Yp_f*spf.nu/spf.uStar_d "y at y+=Yp_f"
Y_const Yp_const*spf.nu/spf.uStar_d "y at y+=Yp_const"
D_bw "min(1.4[mm], Cbw*exp(-(T_sat-T_const)/45[K]))" "Bubble departure diameter"
q_cpl down(phtr.rho_s1)*nojac(spf.uStar_d)*down(nitmm1.Cp_s1)*(ht.Td-nitf1.TdWF)/sqrt(1.0E-6+max(nitf1.Tplus_d,0)^2) "Convective heat flux with pure liquid"
q_cpg down(phtr.rho_s2)*nojac(spf.uStar_d)*down(nitmm1.Cp_s2)*(ht.Td-nitf1.TdWF)/sqrt(1.0E-6+max(nitf1.Tplus_d,0)^2) "Convective heat flux with pure gas"
q_cl A_cpl*q_cpl ""
q_cg A_cpg*q_cpg ""
qwf_d -(q_q+q_cl+q_cg) "Wall function flux, downside"