CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

cyclic boundary conditions and moveMesh

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

LinkBack Thread Tools Display Modes
Old   September 21, 2012, 03:54
Cool cyclic boundary conditions and moveMesh
New Member
Join Date: Jun 2011
Posts: 13
Rep Power: 8
gtg258f is on a distinguished road
Dear all,

I running a simulation with a moving mesh an periodic boundary conditions. For the mesh motion I am using the moveMesh method from elasticNonLinULSolidFoam.

After a few iterations, the code crashes because the periodic patches have moved by a different amount. While I understand that this mismatch makes the code crash, I was wondering why it happens, since I thought that cyclic patches would have the same displacement field.

face 14 area does not match neighbour 29 by 0.0145797% -- possible face ordering problem.
patch:XPERIODIC my area:0.0209228 neighbour area:0.0209197 matching tolerance:0.0001
Looking into the code, I was wondering if the problem could come from performEdgeCorrectedVolPointInterpolation.H
I see that there is a part in there about coupled boundaries that is commented out.

From what I understand (but I could be really wrong), the displacement field at the boundary points is extrapolated from the internal field gradient (gradDU). For cyclic boundaries, should the displacement field in coupled patches (say patch 1 and 2) be computed as a mean of the extrapolations on patch 1 and patch 2? Should I uncomment the following lines?

// Update coupled boundaries
forAll (pfCorr.boundaryField(), patchI)
if (pfCorr.boundaryField()[patchI].coupled())
forAll (pfCorr.boundaryField(), patchI)
if (pfCorr.boundaryField()[patchI].coupled())

gtg258f is offline   Reply With Quote

Old   September 21, 2012, 04:26
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,759
Rep Power: 29
ngj will become famous soon enoughngj will become famous soon enough
Hi Diane,

I am facing the same problem, and as far as I have found out, the reason is the following:

1. If you are using the mesh motion, which solves in the finite volume sense (the methods that come with the OpenCFD releases), then the motion is computed for the cell centres, and the motion is subsequently interpolated to the internal points and extrapolated to the boundary points. This extrapolation, however, does not take the nature of the boundaries into consideration, and the mesh motion on the cyclic patches are not enforced to be identical.

2. If you want to try to solve for the mesh motion using the tetFem-decomposition technique, then you will find that no such thing as a cyclic boundary condition exists (at least up to the most recent update I have of the OF-1.6-ext framework).

Unfortunately, I do not have any quick solution for you, and it seems to require a substantial programming effort to make a solution, which work inside the OF-framework, but you might be able to come up with a hack, but the hack needs to be placed in the mesh motion library before the patch-check for cyclic boundaries is conducted.

Good luck,

ngj is offline   Reply With Quote

Old   November 21, 2013, 10:13
Default cyclicng BC and moving mesh
Join Date: Dec 2010
Posts: 31
Rep Power: 9
levka is on a distinguished road
hi guys,
i am trying to apply cyclic BC to the domain with mesh motion and can not succeed.
Mesh solver (Laplaciancomponent or Displacement) works wonderful with the same domain but when inlet and outlet are some fixed values or zero gradients. But when i change them to cyclic with help of CreatePatch utility mesh solver starts working not corrected. I explain correct - when only object inside domain moves; non correct - when whole domain moves...

Any hit is much appreciated
levka is offline   Reply With Quote


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

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