v_in v_in_O2*6/W_O2^2*(x-W_sep-W_H2)*(W_H2-(x-W_sep-W_H2)) "Inlet velocity"
m_OH we.Ilx/F_const*MOH "Ion mass flow rate"
eps_liquid max(ee.phic,1e-4) "Electrolyte volume fraction"
slipevelmag sqrt(max((udx-ucx)^2+(udy-ucy)^2,eps^2)) "Slip velocity magnitude"
prefac -ee.rhoc*K_O2_div_d*slipevelmag Prefactor
FD_BDx prefac*d(phid,x) "Bubble dispersion volume force, x-component in dispersed phase"
FD_BDy prefac*d(phid,y) "Bubble dispersion volume force, y-component in dispersed phase"
FC_BDx -prefac*d(phid,x)*phid/max(ee.phic,1e-4) "Bubble dispersion volume force, x-component in continuous phase"
FC_BDy -prefac*d(phid,y)*phid/max(ee.phic,1e-4) "Bubble dispersion volume force, y-component in continuous phase"