v_in v_in_H2*6/W_H2^2*x*(W_H2-x) "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_H2_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"