[michujo]

Dear CFD-onliners,

I came across a question when trying to solve the momentum equation following Patankar's book (1980) indications.
When solving a general transport equation for a certain variable , the discretized form takes the form:


, where P,E,W,N,S stand for point, east face, west face, north face and south face, respectively. The coefficients depend on the type of numerical scheme used (central differences, upwind, exponential, power law...).

My question: how do we calculate the discretization coefficients for the momentum equations u and v? can we just assume that the momentum equations are another set of convection-diffusion equations, where the variable are u and v? Then we could just use the same method as before (upwind, power law, ...).

Any help will be welcome.

Thanks in advance .

[cdegroot]

Originally Posted by michujo

Dear CFD-onliners,

I came across a question when trying to solve the momentum equation following Patankar's book (1980) indications.
When solving a general transport equation for a certain variable , the discretized form takes the form:


, where P,E,W,N,S stand for point, east face, west face, north face and south face, respectively. The coefficients depend on the type of numerical scheme used (central differences, upwind, exponential, power law...).

My question: how do we calculate the discretization coefficients for the momentum equations u and v? can we just assume that the momentum equations are another set of convection-diffusion equations, where the variable are u and v? Then we could just use the same method as before (upwind, power law, ...).

Any help will be welcome.

Thanks in advance .

The u and v momentum equations are not exactly like a convection-diffusion equation because of the pressure term. The convection and diffusion terms are discretized in the same way as described for a general convection-diffusion equation, however. Don't forget to think about pressure-velocity coupling as well.

[michujo]

Hi Chris, I found later in Patankar's book that we could just use any differentiation scheme for the diffusive and convective terms in the momentum equation, as it is done for general convection-diffusion equations (I also found the coefficients I was looking for in Versteeg&Malalasekera book).
Yes, I will treat the pressure gradient term differently don't worry.

I appreciate your reply anyway.

Cheers.

[leflix]

Originally Posted by cdegroot

The u and v momentum equations are not exactly like a convection-diffusion equation because of the pressure term. The convection and diffusion terms are discretized in the same way as described for a general convection-diffusion equation, however. Don't forget to think about pressure-velocity coupling as well.

Hi cdegroot,

The momentum equations are EXACTLY convection-diffusion equations with a source term. In this case the source term is the pressure gradient (dp/dx_i) + cross derivative of velocity (du_j/dx_i). (here d is partial derivative)

It can be discretized using the same schemes as the ones used for any scalar transport equation. The scalars here are the componnents of the velocity vector u and v.

[cdegroot]

Originally Posted by leflix

Hi cdegroot,

The momentum equations are EXACTLY convection-diffusion equations with a source term. In this case the source term is the pressure gradient (dp/dx_i) + cross derivative of velocity (du_j/dx_i). (here d is partial derivative)

It can be discretized using the same schemes as the ones used for any scalar transport equation. The scalars here are the componnents of the velocity vector u and v.

Yes, I agree.