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

Odd-even decoupling of cells

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 31, 2001, 17:45
Default Odd-even decoupling of cells
  #1
Peter
Guest
 
Posts: n/a
I have heard using CDS scheme can lead to numerical problems due to odd-even decoupling of cells. could you tell me what is the manifestation of this topic in the convergence history or the solution???

Thanks
  Reply With Quote

Old   June 1, 2001, 10:09
Default Re: Odd-even decoupling of cells
  #2
Patrick Godon
Guest
 
Posts: n/a
Hi,

odd-even decoupling of cells occurs when (for example) you use a low (second) order scheme in space and approximate

du/dx = (u(j+1)-u(j-1))/(2*delta x)

Then obviously, by jumping two indexes (from j-1 to j+1; skipping over j), you decouple the even and odd cells in the whole domain (assuming you don't have any other term).

Say, it is a first order equation, therefore, you need one boundary condition (mathematically speaking).

Next, as you impose a boundary condition on one side you actually solve two problems:

1) you solve for the even cells with one boundary condition on one side

2) you solve for the odd cells without any boundary condition (since you impose here only one boundary condition).

This problem will 'manifest' (will show) by two point oscillations in the solution, which are actually two different solution of the same problem (since for each different boundary condition you actually define - mathematically a new problem).

To overcome this problem, you might use another scheme where you do not skip over 'j'; or you might want to use an additional boundary condition, such that the odd and even cells will be coupled (for example if you give the function on the odd cells at one boundary, you might want to consider giving the derivative of the function at the same boundary for the even cells, such that the expression for the derivative includes also the function given at the odd cell boundary - this way you will have a relation between odd and even cell points). Mathematically speaking, this is OK, since the difference equations are decoupled, this means that you have actually two equations of the first order in space and need one boundary condition for each of them.

Patrick Godon
  Reply With Quote

Old   June 1, 2001, 13:19
Default Re: Odd-even decoupling of cells
  #3
Peter
Guest
 
Posts: n/a
Chees, thanks for your comments!!!
  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Import netgen mesh to OpenFOAM hsieh Open Source Meshers: Gmsh, Netgen, CGNS, ... 32 September 13, 2011 05:50
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 04:23
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43
physical boundary error!! kris CD-adapco 2 August 3, 2005 00:32


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