CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Coefficients discretized momentum equation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 15, 2012, 09:20
Default Coefficients discretized momentum equation
  #1
Senior Member
 
Join Date: Dec 2011
Location: Madrid, Spain
Posts: 134
Rep Power: 15
michujo is on a distinguished road
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 \phi, the discretized form takes the form:
a_P \phi_P=a_E \phi_E+a_W \phi_W+a_N \phi_N+a_S \phi_S+b

, 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 \phi 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 .
michujo is offline   Reply With Quote

Old   June 18, 2012, 16:07
Default
  #2
Senior Member
 
cdegroot's Avatar
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17
cdegroot is on a distinguished road
Quote:
Originally Posted by michujo View Post
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 \phi, the discretized form takes the form:
a_P \phi_P=a_E \phi_E+a_W \phi_W+a_N \phi_N+a_S \phi_S+b

, 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 \phi 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.
cdegroot is offline   Reply With Quote

Old   June 18, 2012, 18:49
Default
  #3
Senior Member
 
Join Date: Dec 2011
Location: Madrid, Spain
Posts: 134
Rep Power: 15
michujo is on a distinguished road
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.
michujo is offline   Reply With Quote

Old   June 19, 2012, 20:16
Default
  #4
Senior Member
 
Join Date: Aug 2011
Posts: 272
Rep Power: 15
leflix is on a distinguished road
Quote:
Originally Posted by cdegroot View Post
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.
leflix is offline   Reply With Quote

Old   June 20, 2012, 02:33
Default
  #5
Senior Member
 
cdegroot's Avatar
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17
cdegroot is on a distinguished road
Quote:
Originally Posted by leflix View Post
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.
cdegroot is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
solving the momentum equation in UEqn. callahance OpenFOAM 2 October 18, 2012 10:38
Question on the discretization of momentum equation in icoFoam MPJ OpenFOAM 3 October 4, 2011 10:44
Momentum equation of Darcy's law sambatra OpenFOAM Running, Solving & CFD 2 April 17, 2009 04:27
Momentum Equation Munir Ahmed Khan FLUENT 0 August 27, 2008 08:57
Question on momentum equation in VOF or Level Set Kai Yan Main CFD Forum 0 January 11, 2008 09:47


All times are GMT -4. The time now is 07:25.