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

Cyclic boundary conditions in parallel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 29, 2011, 05:26
Default Cyclic boundary conditions in parallel
  #1
New Member
 
Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 7
thibault_pringuey is on a distinguished road
Hello,

I am using OpenFOAM 1.5.x

In order to check the order of a numerical scheme, I am translating a sine field ( sin(pi * (pos().x+pos().y)) ) in a 2D cyclic domain. Such that at time = 1s the calculated field should overlap the initial field.

While this case runs fine on a single processor, in parallel (simple method; 4 processors: 2x2; original 0-halo decomposition) I notice a shift of the sine (~2 cell diagonals) at the end of the calculation. The shift of the field is forward, i.e. it presents the solution as it should look at a later time.

When using the n-halo parallel processing I developed, the shift is divided by 2. It is worth mentioning that for this 2x2 decomposition, using the n-halo decomposition effectively divide by 2 the number of processorXtoY boundaries per sub domain.

I would really appreciate your thoughts on this. I have attached a contour plot of the initial field in 2D and the plot over the diagonal for serial, 0-halo and n-halo calculations.

Many thanks in advance for your help.

With best wishes,


Thibault
Attached Images
File Type: jpg serialDiagonalPlot.jpg (22.4 KB, 9 views)
File Type: jpg 0haloDiagonalPlot.jpg (30.5 KB, 10 views)
File Type: jpg nhaloDiagonalPlot.jpg (28.7 KB, 11 views)
File Type: jpg InitialField.jpg (41.6 KB, 8 views)
thibault_pringuey is offline   Reply With Quote

Old   April 2, 2011, 15:21
Default Solved
  #2
New Member
 
Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 7
thibault_pringuey is on a distinguished road
Solved it!

For some reason the boundaryField stored for processor patches (version 1,5) is the neighbourField and not the average of the internalField and the boundaryField as for the cyclic patches.
thibault_pringuey is offline   Reply With Quote

Reply

Tags
boundary conditions, cyclic, parallel

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
Domain Imbalance HMR CFX 3 March 6, 2011 20:10
CFX doesn't continue calculation... mactech001 CFX 6 November 15, 2009 21:25
Fluent accuracy and boundary conditions Paolo Lampitella FLUENT 0 June 12, 2008 06:25
periodic and cyclic boundary conditions Hadi Main CFD Forum 2 June 29, 2007 07:19
A problem about setting boundary conditions lyang Main CFD Forum 0 September 19, 1999 18:29


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