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

Odd-Even Decoupling on 1D 7 pt. Stencil Finite Differences?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 4, 2018, 03:37
Default Odd-Even Decoupling on 1D 7 pt. Stencil Finite Differences?
  #1
New Member
 
Ali Berk Kahraman
Join Date: Dec 2015
Location: Braunschweig, Germany
Posts: 13
Rep Power: 10
abkahraman is on a distinguished road
Dear All,

I am writing my code from scratch for incompressible Navier-Stokes. I will use 7 pt. stencil Finite Differences on a given direction, so a total 13 pt stencil on 2D. What are your ideas on odd-even decoupling?

My guess is that since odd and even nodes are able to communicate with each other because the only zero entry will be the center node for which I have to calculate the first derivative. However, I cannot be sure. Any insights or advices?

Best Regards,

Ali
abkahraman is offline   Reply With Quote

Old   March 4, 2018, 07:15
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by abkahraman View Post
Dear All,

I am writing my code from scratch for incompressible Navier-Stokes. I will use 7 pt. stencil Finite Differences on a given direction, so a total 13 pt stencil on 2D. What are your ideas on odd-even decoupling?

My guess is that since odd and even nodes are able to communicate with each other because the only zero entry will be the center node for which I have to calculate the first derivative. However, I cannot be sure. Any insights or advices?

Best Regards,

Ali
But is your stencil used for the pressure equation? how do you apply the div grad product in discrete form?
FMDenaro is offline   Reply With Quote

Old   March 4, 2018, 08:10
Default
  #3
New Member
 
Ali Berk Kahraman
Join Date: Dec 2015
Location: Braunschweig, Germany
Posts: 13
Rep Power: 10
abkahraman is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
But is your stencil used for the pressure equation? how do you apply the div grad product in discrete form?
Yes, the stencil is also used for the pressure equation, no staggered grid. The paper I read suggests applying upwind discretization for div and downwind for grad or vice versa. Thus, odd-even decoupling would not occur.

My thinking is that div dot grad product is mathematically Laplace operator. So, I thought of using 7 pt. stencil on the Laplace operator.

This brings another question for which I have not looked up the answer, why bother with div dot grad when it is simply Laplacian?
abkahraman is offline   Reply With Quote

Old   March 4, 2018, 08:18
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Matemathically the product corresponds to the Laplacian. However, numericall there are several well known issues that drive to an Approximate or an Exact projection method.
Consider that you have to compute the gradient of the pressure to ensure the continuity constraint. I suggest to see by hand in 1D the differences if you discretize the div grad and the Lap operators.
You can find many references in the CFD literature
FMDenaro is offline   Reply With Quote

Old   March 4, 2018, 08:33
Default
  #5
New Member
 
Ali Berk Kahraman
Join Date: Dec 2015
Location: Braunschweig, Germany
Posts: 13
Rep Power: 10
abkahraman is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Matemathically the product corresponds to the Laplacian. However, numericall there are several well known issues that drive to an Approximate or an Exact projection method.
Consider that you have to compute the gradient of the pressure to ensure the continuity constraint. I suggest to see by hand in 1D the differences if you discretize the div grad and the Lap operators.
You can find many references in the CFD literature
Thank you for the explanation. Any ideas on odd-even decoupling?
abkahraman is offline   Reply With Quote

Old   March 4, 2018, 08:44
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by abkahraman View Post
Thank you for the explanation. Any ideas on odd-even decoupling?
It happens for some type of discretizations on non-staggered grids. A classical example is the second order discretization. Other examples can be found in the literature. I never used a 7-point stencil for the gradient operator so you have to develop the scheme. But I immagine that the implementation will be very difficult close to the boundaries..
FMDenaro 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
finite differences multilayer diffusion wilmermantilla Main CFD Forum 0 September 30, 2017 20:18
Continuous adjoint VS finite differences Haowang SU2 0 March 26, 2017 04:59
Lid Driven Cavity - Finite Differences BCs? aem0512 Main CFD Forum 3 May 1, 2016 23:50
conservative finite differences and finite volumes Joachim Main CFD Forum 7 January 23, 2014 17:32
Finite differences Zdravko Stojanovic Main CFD Forum 8 March 24, 2001 23:31


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