[Nek5000-users] Calculation of Divergence Term

[nek5000-users at lists.mcs.anl.gov]

Hi,

I want to solve an advection-diffusion equation in conjunction with  
the momentum and mass conservation equations. But my equation has an  
extra divergence term which I don't know how to calculate it in  
Nekton. So, denoting the scalar concentration field with c and  
velocity field with u, the equation turns out to be:

dc/dt + u.grad c = Laplacian c - div (c V)

which by d/dt I mean partial derivative with respect to time. V is a  
known vector field calculated from the vorticity tensor. I wonder if  
you could show me how to evaluate this additional term (the last term)  
in the user file.


Best Regards,
Alireza Karimi