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

Discretization of the convective term of the N-S equation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 31, 2017, 07:06
Default Discretization of the convective term of the N-S equation
  #1
New Member
 
Donton
Join Date: Aug 2017
Posts: 6
Rep Power: 9
Donton is on a distinguished road
Hello everyone...

In most of the document explaining the Fractional Step Method I have found that the non-dimansionl for of the Navier-Stokes equation is discretized as,

(u^n+1 - u^n)/dt + (3/2)*H(u^n) - (1/2)*H(u^n-1) = ... the right hand side

And it is mentioned the H is the discrete convective operator. Now my question is how to expand this H?

I have done in the following way. Please tell me if I am incorrect.
Considering the non-dimentional form of the x-momentum equation

H(u^n) = ( (u^n * u^n)|i+1/2,j - (u^n * u^n)|i-1/2,j )/dx + ( (v^n * u^n)|i,j+1/2 - (v^n * u^n)|i,j-1/2 )/dy

here ^n denotes the superscript n for time level n
Now if we write the equation without the superscript just for simlicity,

H(u^n) = ( (u*u)|i+1/2,j - (u*u)|i-1/2,j )/dx + ( (v*u)|i,j+1/2 - (v*u)|i,j-1/2 )/dy

Now considering uniform structurred staggered grid arrangement and central difference scheme,

(u*u)|i+1/2,j = u|i+1/2,j * u|i+1/2,j
= ((u|i+1,j + u|i,j)/2) * ((u|i+1,j + u|i,j)/2)

(u*u)|i-1/2,j = u|i-1/2,j * u|i-1/2,j
= ((u|i-1,j + u|i,j)/2) * ((u|i-1,j + u|i,j)/2)

(v*u)|i,j+1/2 = v|i,j+1/2 * u|i,j+1/2
= ((v|i+1,j + v|i,j)/2) * ((u|i,j+1 + u|i,j)/2)

(v*u)|i,j-1/2 = v|i,j-1/2 * u|i,j-1/2
= ((v|i+1,j-1 + v|i,j-1)/2) * ((u|i,j-1 + u|i,j)/2)

Is this discretization correct? My next question is if I want to use upwind scheme how will it be?

Also in Fractional step method the convective term is estimated based on the previous time step's values. Therefore while solving velocity at each point from the momentum equation for the current time step the convective term will just become a known quantity and will go to the right hand side of the algebriac equations. Then how does it effect the solution wheather central difference or upwinding is used for the discretization of the convective term?
Donton is offline   Reply With Quote

Old   September 1, 2017, 13:01
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Well, what you see is the explicit multi-step Adams-Bashforth second order accurate time integration applied only for the convective term vector H=Div(vv). This notation is used in the hystorical work of Kim and Moin, (for example see https://ntrs.nasa.gov/archive/nasa/c...9840014260.pdf)
Then, you can introduce any type of spatial discretization you want to H^n and H^n-1. Of course, different final algebric terms appear if you use staggered or non-staggered grids, central or upwind schemes of various order. Any consistent discretization is correct, there is no a unique possible discretization. And that does not enter into the Fractional Step Method.
FMDenaro is offline   Reply With Quote

Reply


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
Adding diffusion term to interFoam transport equation Gearb0x OpenFOAM Programming & Development 3 February 14, 2023 05:16
Discretization of the diffusion term mseka Main CFD Forum 5 August 30, 2017 11:41
how to troduce a source term in ICOFOAM solver's equation sawyer86 OpenFOAM Running, Solving & CFD 0 July 20, 2012 12:15
Bouyancy term in Energy equation CFDtoy Main CFD Forum 2 August 4, 2008 12:04
bouyancy term in epsilon equation Michael Main CFD Forum 1 June 25, 1999 11:20


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