[otq]

I'm trying to solve only certain directions of the continuity equation, but still in a 3D environment.

I understand that I can take the component of a vector with the following:
U.component(2) would return the z velocity.

Now I want to solve the z component of continuity:

rho(dw/dz)+w(drho/dz)=0
Which fvm:: function should I use to solve for w? I know the equation will look like:

fvm::ddz(w) == w/rho*fvc::grad(rho)&mesh.C().component(2)

or do I need to copy the ddt template and make my own ddz function?

[haze_1986]

Hi, do you have any progress on this issue? Would like to know too.