|
[Sponsors] |
Discretization of the convective term of the N-S equation |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 31, 2017, 07:06 |
Discretization of the convective term of the N-S equation
|
#1 |
New Member
Donton
Join Date: Aug 2017
Posts: 6
Rep Power: 9 |
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? |
|
September 1, 2017, 13:01 |
|
#2 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
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. |
|
|
|
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 |