I am trying to model a simple porous media flow through a certain type of wood or a foam. I have a generalized capillary - saturation relationship from sources: P_c (S_w) = AS_w ^ (-B)
with initial conditions in space as S_w = 0, when I model the transient flow, the term dP_c/dS_w becomes infinite because the computation becomes 0^-(b+1).
Is there a way to get around this inconsistency. I only option in my mind is to apply a very small insignificant value of S_w (1e-100) as initial condition.
Any suggestion will be highly helpful. Thank you
Regards
Sai Kiran