"" pf.epsilon_pf 2*pf.

ep_default m "Parameter controlling interface thickness"

"Domains 1–2" ""
"" pf.dfdphi 0 J/m³ "phi-derivative of external free energy" "Domains 1–2" ""
"" pf.ur model.input.u1 m/s "Velocity field, r component" "Domains 1–2" Meta
"" pf.uphi model.input.u2 m/s "Velocity field, phi component" "Domains 1–2" Meta
"" pf.uz model.input.u3 m/s "Velocity field, z component" "Domains 1–2" Meta
"" pf.sigma model.input.sigma N/m "Surface tension coefficient" "Domains 1–2" Meta
"" phipfrt phipfrTIME-phipfrr*d(r,TIME)-phipfrz*d(z,TIME) 1/(m·s) "Gradient of
phipf, r component, first time derivative" "Domains 1–2" ""
"" phipfzt phipfzTIME-phipfzr*d(r,TIME)-phipfzz*d(z,TIME) 1/(m·s) "Gradient of
phipf, z component, first time derivative" "Domains 1–2" ""
"" phipfrtt d(phipfrTIME-phipfrr*d(r,TIME)-phipfrz*d(z,TIME),TIME)-d(phipfrTIME-
phipfrr*d(r,TIME)-phipfrz*d(z,TIME),r)*d(r,TIME)-d(phipfrTIME-phipfrr*d(r,TIME)-
phipfrz*d(z,TIME),z)*d(z,TIME) 1/(m·s²) "Gradient of phipf, r component, second
time derivative" "Domains 1–2" ""
"" phipfztt d(phipfzTIME-phipfzr*d(r,TIME)-phipfzz*d(z,TIME),TIME)-d(phipfzTIME-
phipfzr*d(r,TIME)-phipfzz*d(z,TIME),r)*d(r,TIME)-d(phipfzTIME-phipfzr*d(r,TIME)-
phipfzz*d(z,TIME),z)*d(z,TIME) 1/(m·s²) "Gradient of phipf, z component, second
time derivative" "Domains 1–2" ""
"" phipft phipfTIME-phipfr*d(r,TIME)-phipfz*d(z,TIME) 1/s "Phase field variable,
first time derivative" "Domains 1–2" ""
"" phipftt d(phipfTIME-phipfr*d(r,TIME)-phipfz*d(z,TIME),TIME)-d(phipfTIME-
phipfr*d(r,TIME)-phipfz*d(z,TIME),r)*d(r,TIME)-d(phipfTIME-phipfr*d(r,TIME)-
phipfz*d(z,TIME),z)*d(z,TIME) 1/s² "Phase field variable, second time derivative"
"Domains 1–2" ""
"" psirt psirTIME-psirr*d(r,TIME)-psirz*d(z,TIME) 1/(m·s) "Gradient of psi, r
component, first time derivative" "Domains 1–2" ""
"" psizt psizTIME-psizr*d(r,TIME)-psizz*d(z,TIME) 1/(m·s) "Gradient of psi, z
component, first time derivative" "Domains 1–2" ""
"" psirtt d(psirTIME-psirr*d(r,TIME)-psirz*d(z,TIME),TIME)-d(psirTIME-
psirr*d(r,TIME)-psirz*d(z,TIME),r)*d(r,TIME)-d(psirTIME-psirr*d(r,TIME)-
psirz*d(z,TIME),z)*d(z,TIME) 1/(m·s²) "Gradient of psi, r component, second time
derivative" "Domains 1–2" ""
"" psiztt d(psizTIME-psizr*d(r,TIME)-psizz*d(z,TIME),TIME)-d(psizTIME-
psizr*d(r,TIME)-psizz*d(z,TIME),r)*d(r,TIME)-d(psizTIME-psizr*d(r,TIME)-
psizz*d(z,TIME),z)*d(z,TIME) 1/(m·s²) "Gradient of psi, z component, second time
derivative" "Domains 1–2" ""
"" psit psiTIME-psir*d(r,TIME)-psiz*d(z,TIME) 1/s "Phase field help variable, first
time derivative" "Domains 1–2" ""
"" psitt d(psiTIME-psir*d(r,TIME)-psiz*d(z,TIME),TIME)-d(psiTIME-psir*d(r,TIME)-
psiz*d(z,TIME),r)*d(r,TIME)-d(psiTIME-psir*d(r,TIME)-psiz*d(z,TIME),z)*d(z,TIME)
1/s² "Phase field help variable, second time derivative" "Domains 1–2" ""
"" pf.mobility pf.epsilon_pf^2*pf.chi m³·s/kg Mobility "Domains 1–2" ""
"" pf.lam 3*pf.sigma*pf.epsilon_pf/sqrt(8) N "Mixing energy density" "Domains 1–2"
""
"" pf.Vf1 1-pf.Vf2 1 "Volume fraction of fluid 1" "Domains 1–2" ""
"" pf.Vf2 min(1,max(0.5*(1+phipf),0)) 1 "Volume fraction of fluid 2" "Domains 1–2"
""
"" pf.phipfReg min(1,max(phipf,-1)) 1 "Regularized phase field variable" "Domains
1–2" ""
"" pf.chi 1 m·s/kg "Mobility tuning parameter" "Domains 1–2" ""
"" pf.dflux_phipfr pf.mobility*pf.lam*psir/pf.epsilon_pf^2 m/s "Help variable, r
component" "Domains 1–2" ""
"" pf.dflux_phipfphi 0 m/s "Help variable, phi component" "Domains 1–2" ""
"" pf.dflux_phipfz pf.mobility*pf.lam*psiz/pf.epsilon_pf^2 m/s "Help variable, z
component" "Domains 1–2" ""
"" pf.dflux_phipfr 0.5*(up(pf.dflux_phipfr)+down(pf.dflux_phipfr)) m/s "Help
variable, r component" "Boundary 4" ""
"" pf.dflux_phipfphi 0.5*(up(pf.dflux_phipfphi)+down(pf.dflux_phipfphi)) m/s "Help
variable, phi component" "Boundary 4" ""
"" pf.dflux_phipfz 0.5*(up(pf.dflux_phipfz)+down(pf.dflux_phipfz)) m/s "Help
variable, z component" "Boundary 4" ""
"" pf.dflux_phipfr down(pf.dflux_phipfr) m/s "Help variable, r component"
"Boundaries 1–3, 5–10" ""
"" pf.dflux_phipfphi down(pf.dflux_phipfphi) m/s "Help variable, phi component"
"Boundaries 1–3, 5–10" ""
"" pf.dflux_phipfz down(pf.dflux_phipfz) m/s "Help variable, z component"
"Boundaries 1–3, 5–10" ""
"" pf.dfluxtest_phipfr pf.mobility*pf.lam*test(psir)/pf.epsilon_pf^2 m/s "Help
variable, r component" "Domains 1–2" ""
"" pf.dfluxtest_phipfphi 0 m/s "Help variable, phi component" "Domains 1–2" ""
"" pf.dfluxtest_phipfz pf.mobility*pf.lam*test(psiz)/pf.epsilon_pf^2 m/s "Help
variable, z component" "Domains 1–2" ""
"" pf.dfluxtest_phipfr 0.5*(up(pf.dfluxtest_phipfr)+down(pf.dfluxtest_phipfr)) m/s
"Help variable, r component" "Boundary 4" ""
"" pf.dfluxtest_phipfphi 0.5*(up(pf.dfluxtest_phipfphi)
+down(pf.dfluxtest_phipfphi)) m/s "Help variable, phi component" "Boundary 4" ""
"" pf.dfluxtest_phipfz 0.5*(up(pf.dfluxtest_phipfz)+down(pf.dfluxtest_phipfz)) m/s
"Help variable, z component" "Boundary 4" ""
"" pf.dfluxtest_phipfr down(pf.dfluxtest_phipfr) m/s "Help variable, r component"
"Boundaries 1–3, 5–10" ""
"" pf.dfluxtest_phipfphi down(pf.dfluxtest_phipfphi) m/s "Help variable, phi
component" "Boundaries 1–3, 5–10" ""
"" pf.dfluxtest_phipfz down(pf.dfluxtest_phipfz) m/s "Help variable, z component"
"Boundaries 1–3, 5–10" ""
"" pf.dflux_psir pf.epsilon_pf^2*phipfr m "Help variable, r component" "Domains 1–
2" ""
"" pf.dflux_psiphi 0 m "Help variable, phi component" "Domains 1–2" ""
"" pf.dflux_psiz pf.epsilon_pf^2*phipfz m "Help variable, z component" "Domains 1–
2" ""
"" pf.dflux_psir 0.5*(up(pf.dflux_psir)+down(pf.dflux_psir)) m "Help variable, r
component" "Boundary 4" ""
"" pf.dflux_psiphi 0.5*(up(pf.dflux_psiphi)+down(pf.dflux_psiphi)) m "Help
variable, phi component" "Boundary 4" ""
"" pf.dflux_psiz 0.5*(up(pf.dflux_psiz)+down(pf.dflux_psiz)) m "Help variable, z
component" "Boundary 4" ""
"" pf.dflux_psir down(pf.dflux_psir) m "Help variable, r component" "Boundaries 1–
3, 5–10" ""
"" pf.dflux_psiphi down(pf.dflux_psiphi) m "Help variable, phi component"
"Boundaries 1–3, 5–10" ""
"" pf.dflux_psiz down(pf.dflux_psiz) m "Help variable, z component" "Boundaries 1–
3, 5–10" ""
"" pf.dfluxtest_psir pf.epsilon_pf^2*test(psir) m "Help variable, r component"
"Domains 1–2" ""
"" pf.dfluxtest_psiphi 0 m "Help variable, phi component" "Domains 1–2" ""
"" pf.dfluxtest_psiz pf.epsilon_pf^2*test(psiz) m "Help variable, z component"
"Domains 1–2" ""
"" pf.dfluxtest_psir 0.5*(up(pf.dfluxtest_psir)+down(pf.dfluxtest_psir)) m "Help
variable, r component" "Boundary 4" ""
"" pf.dfluxtest_psiphi 0.5*(up(pf.dfluxtest_psiphi)+down(pf.dfluxtest_psiphi)) m
"Help variable, phi component" "Boundary 4" ""
"" pf.dfluxtest_psiz 0.5*(up(pf.dfluxtest_psiz)+down(pf.dfluxtest_psiz)) m "Help
variable, z component" "Boundary 4" ""
"" pf.dfluxtest_psir down(pf.dfluxtest_psir) m "Help variable, r component"
"Boundaries 1–3, 5–10" ""
"" pf.dfluxtest_psiphi down(pf.dfluxtest_psiphi) m "Help variable, phi component"
"Boundaries 1–3, 5–10" ""
"" pf.dfluxtest_psiz down(pf.dfluxtest_psiz) m "Help variable, z component"
"Boundaries 1–3, 5–10" ""
"" pf.intnormr ppr(phipfr/sqrt(phipfr^2+phipfz^2+eps/h^2)) 1 "Interface normal, r
component" "Domains 1–2" ""
"" pf.intnormphi 0 1 "Interface normal, phi component" "Domains 1–2" ""
"" pf.intnormz ppr(phipfz/sqrt(phipfr^2+phipfz^2+eps/h^2)) 1 "Interface normal, z
component" "Domains 1–2" ""
"" pf.delta 0.75*max(1-phipf^2,0)*pf.gradphi 1/m "Fluid interface delta function"
"Domains 1–2" ""
"" pf.gradphi sqrt(phipfr^2+phipfz^2+eps) 1/m "Gradient of ϕ" "Domains 1–2" ""
"" pf.lref 10*pf.epsilon_pf m "Reference length for two phase flow" "Domains 1–2"
""
"" pf.G0i nojac(2/pf.lref) 1/m "Auxiliary wall variable" "Domains 1–2" ""
"" GIrt GIrTIME-GIrr*d(r,TIME)-GIrz*d(z,TIME) 1/(m²·s) "Gradient of GI, r
component, first time derivative" "Domains 1–2" ""
"" GIzt GIzTIME-GIzr*d(r,TIME)-GIzz*d(z,TIME) 1/(m²·s) "Gradient of GI, z
component, first time derivative" "Domains 1–2" ""
"" GIrtt d(GIrTIME-GIrr*d(r,TIME)-GIrz*d(z,TIME),TIME)-d(GIrTIME-GIrr*d(r,TIME)-
GIrz*d(z,TIME),r)*d(r,TIME)-d(GIrTIME-GIrr*d(r,TIME)-GIrz*d(z,TIME),z)*d(z,TIME) 1/
(m²·s²) "Gradient of GI, r component, second time derivative" "Domains 1–2" ""
"" GIztt d(GIzTIME-GIzr*d(r,TIME)-GIzz*d(z,TIME),TIME)-d(GIzTIME-GIzr*d(r,TIME)-
GIzz*d(z,TIME),r)*d(r,TIME)-d(GIzTIME-GIzr*d(r,TIME)-GIzz*d(z,TIME),z)*d(z,TIME) 1/
(m²·s²) "Gradient of GI, z component, second time derivative" "Domains 1–2" ""
"" GIt GITIME-GIr*d(r,TIME)-GIz*d(z,TIME) 1/(m·s) "Reciprocal initial interface
distance, first time derivative" "Domains 1–2" ""
"" GItt d(GITIME-GIr*d(r,TIME)-GIz*d(z,TIME),TIME)-d(GITIME-GIr*d(r,TIME)-
GIz*d(z,TIME),r)*d(r,TIME)-d(GITIME-GIr*d(r,TIME)-GIz*d(z,TIME),z)*d(z,TIME)
1/(m·s²) "Reciprocal initial interface distance, second time derivative" "Domains
1–2" ""
"" pf.Ddirir pf.DdirVarr 1 "Direction toward nearest initial interface, r
component" "Domains 1–2" ""
"" pf.Ddiriphi pf.DdirVarphi 1 "Direction toward nearest initial interface, phi
component" "Domains 1–2" ""
"" pf.Ddiriz pf.DdirVarz 1 "Direction toward nearest initial interface, z
component" "Domains 1–2" ""
"" pf.Dwi pf.DwVar m "Initial interface distance" "Domains 1–2" ""
"" pf.DdirVarr GIr/sqrt(max(eps,GIr^2+GIz^2)) 1 "Help variable, r component"
"Domains 1–2" ""
"" pf.DdirVarphi 0 1 "Help variable, phi component" "Domains 1–2" ""
"" pf.DdirVarz GIz/sqrt(max(eps,GIr^2+GIz^2)) 1 "Help variable, z component"
"Domains 1–2" ""
"" pf.DwVar 1/GI-1/pf.G0i m "Help variable" "Domains 1–2" ""