Three phase flow problem

I'm getting the following errors when I try simulating droplet separation process. I've attached my file for reference.

Failed to find consistent initial values.
Attempt to evaluate real square root of negative number.
Function: sqrt
Failed to evaluate expression.
Expression: nojac(1/sqrt((54*emetric2)+((4*comp1.spf.rho)*emetric(comp1.u,comp1.v))))
Failed to evaluate variable.
Variable: comp1.tfpf1.tau mCross
Failed to evaluate variable.
Variable: comp1.tfpf1.vhNS lin, Defined as: (sqrt(max((((comp1.tfpf1.xi m*comp1.tfpf1.res u)^2)+((comp1.tfpf1.xi m*comp1.tfpf1.res v)^2))*comp1.tfpf1.tau mCross,0)*comp1.tfpf1.B0YigijtauNorm2)*0.5)/(comp1.tfpf1.B0YigijtauNorm2+eps)
Failed to evaluate expression.
Expression: nojac((1.4142135623730951^(1-shapeorder(comp1.u)))*max(comp1.tfpf1.vhNS lin-comp1.tfpf1.vhNS quad,0))
Failed to evaluate variable.
Variable: comp1.spf.vhNS
Failed to evaluate temporary symbolic derivative variable.
Variable: comp1.spf.crosswindnsVDN${test24}
Failed to evaluate expression.
Expression: d((comp1.spf.crosswindns)*(dvol),{test24})
Failed to evaluate Jacobian of expression.
Expression: (comp1.spf.crosswindns)*(dvol)
Last time step is not converged.