# Cyclic boundary conditions in parallel

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

March 29, 2011, 05:26
Cyclic boundary conditions in parallel
#1
New Member

Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 9
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.

With best wishes,

Thibault
Attached Images
 serialDiagonalPlot.jpg (22.4 KB, 10 views) 0haloDiagonalPlot.jpg (30.5 KB, 11 views) nhaloDiagonalPlot.jpg (28.7 KB, 12 views) InitialField.jpg (41.6 KB, 10 views)

 April 2, 2011, 15:21 Solved #2 New Member   Thibault Pringuey Join Date: Mar 2009 Posts: 17 Rep Power: 9 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.

 Tags boundary conditions, cyclic, parallel

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post HMR CFX 5 October 10, 2016 05:57 mactech001 CFX 6 November 15, 2009 22:25 Paolo Lampitella FLUENT 0 June 12, 2008 06:25 Hadi Main CFD Forum 2 June 29, 2007 07:19 lyang Main CFD Forum 0 September 19, 1999 18:29

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