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

More Sai Kiran Hota's questions See All
Similar questions and discussions