【COMSOL】案例 模拟沸腾水中气泡运动 两相流 流体传热 蒸汽 冷凝

kappa (k_L-k_v)*tpf1.Vf2+k_v "Thermal conductivity" Cp (Cp_L-Cp_v)*tpf1.Vf2+Cp_v "Specific heat" k_v 8.3154e-5[W/m/K/K]*T-7.45556e-3[W/m/K] "Thermal conductivity of vapor" rho_v (p0*Mw)/(8.314[J/mol/K]*(T+eps)) "Density of vapor" q q0+0.1*q0*sin(2*pi*f*t[1/s]) "Heat flux applied to cavity" q0 "1e5 [W/m^2]" "Heat flux term" f 0.2 "Frequency of bubble creation" mdot if(T<T_sat,R*rho_v*(T-T_sat)/T_sat,R*rho_L*(T-T_sat)/T_sat) "Rate of vaporization" delta 6*tpf1.Vf2*(1-tpf1.Vf2)*0.5*sqrt(phipfx^2+phipfy^2+eps) "Interface delta function" phi_source -mdot*delta*(tpf1.Vf1/tpf1.rho1+tpf1.Vf2/tpf1.rho2) "Source term for level set equation" usource mdot*delta*(1/tpf1.rho1-1/tpf1.rho2) "Source term in continuity equation" Qs -mdot*L*delta "Source term in energy equation"