Adding transport equation for tangential velocity in simpleFoam
 

Hello everyone.

I am simulating an 2D-axisymmetric case with simpleFoam. The simulation worked fine with simpleFoam.

But now I want to add a new equation for tangential component to the solver.

the transport equation for "w" is,

div(phi, w) = Laplacian(nu, w) + nu/y*(dw/dy) -nu/y*(w/y)

Can anybody give me an idea how to put the second and third
term in the code....
fvScalarMatrix UwEqn
(
fvm::div(phi, w)
- fvm::laplacian(nu, w)
.............. second term
.............. third term
);
I have attached a file for convenience. any suggestion would be welcomed!

Best regards,

Sohail Iqbal