||April 28, 2010 08:14
Additional terms in N-S equs
I am working with 2phase solid-gas flow and I want to add a term to N-S equations as follows:
div(rhos*epsilons*Us*Us) = -epsilons*grad(p) - beta*(Us-Ug) + div(taus) + rhos*epsilons*g
as you can see I have two additional terms in N-S equs as beta*(Us-Ug) and div(taus). beta in the first term is the momentum transfer, taus is the stress tensor and epsilons is the volume fraction of solid phase.
I am wondering how I can add these terms to N-S equs in PHOENICS. mostly the term beta*(Us-Ug) is a simple argument and must be very easy to add but I do not know how and I would appreciate if anyone can help me with this.