safavinick_reza 
April 28, 2010 08:14 
Additional terms in NS equs
Hi
I am working with 2phase solidgas flow and I want to add a term to NS equations as follows:
div(rhos*epsilons*Us*Us) = epsilons*grad(p)  beta*(UsUg) + div(taus) + rhos*epsilons*g
as you can see I have two additional terms in NS equs as beta*(UsUg) 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 NS equs in PHOENICS. mostly the term beta*(UsUg) 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.
