CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Mapped vs Cyclic Boundaries in OpenFOAM

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

Like Tree4Likes
  • 4 Post By arsenis

LinkBack Thread Tools Search this Thread Display Modes
Old   July 31, 2017, 08:17
Default Mapped vs Cyclic Boundaries in OpenFOAM
New Member
Arsenis Chatzimichailidis
Join Date: Oct 2014
Posts: 6
Rep Power: 11
arsenis is on a distinguished road
Hi to all,

I use OpenFOAM 2.3.1 to simulate with LES the pollutant dispersion inside street canyons. A modified pisoFoam solver (to include the convection-diffusion equation) is used.

At the beginning I used cyclic boundaries, but changed them to mapped, because I don't want to recycle and thus capture the pollutants inside the domain. I know that there are tools in swak4Foam and in newer versions (e.g. fixedValueConstraint in 4.1), to use cyclic and don't recycle the pollutants. Unfortunately, I am not familiar with the first (do I need to install it or just compile it) and any newer OpenFOAM versions won't be installed in our cluster for the near future.

I have a question about the mapped boundary conditions in OpenFOAM:

Should all the used pisoFoam properties (i.e. U, p and nuSgs) be recycled through mapping, or just the velocity U?

Reading in the forum and searching in Google didn't give a specific answer, so I present some of my results, so anyone may observe something that I have not.

In the pitzDailyMapped tutorial in the OpenFOAM installations, the U, k and nuTilda are recycled, while p and nuSgs are not. So, I have tested both cases:
a) recycle only U and
b) recycle U, p and nuSgs

For some geometries both cases have worked, with nearly identical results and both can be accepted as validated.

However, for the case b, an irregular flow develops at the inlet. The following images present this irregularity. The pressure is shown, as it is more visible than velocity, from random time steps of the studied cases.

Street canyon with aspect ratio 1/1. The irregularity is barely visible and no problems are caused in the simulation (over 2000 seconds).

Street canyon with aspect ratio 1/2. The irregularity is fully visible, but no problems have been caused in the simulation (after 25 seconds).

Street canyon with aspect ratio 1/3. The irregularity is fully visible and now the simulation crashes after ~6 seconds. A weird recycling is also observed at the inlet patch and the velocity increases constantly in the same location.
the pressure p

and the velocity U at the same time step

Another image of the same problem, in some of my attempts.

All the above flow fields have been obtained with cyclic boundaries in the beginning and then the solution was mapped (mapFields utility) in a domain with mapped boundaries. In all the above cases, the setAverage option was set to false. I have also tested variations on the offset of the mapped patches, with no change in the results (e.g. for the width of 40m, the 39, 39.5 and 41m were tested).

For the case of recycling only the velocity U, when the initial flow field is mapped to a domain with the setAverage option set to off, a significant drop in U is observed with no other problems. When mapped to domain with the setAverage option set to on, the solution crashed. I suspect that for this option, the mapped patch expects a certain velocity profile.

On the other hand, I am currently obtaining a new solution (recycling only U and the setAverage option on) from the beginning and I don't have any problems.
I must clarify, that I used this method, when started to use mapped boundaries with no problems. I temporarily abandoned it, because I cannot fully understand from the code what is the role of setAverage (not enough programming knowledge!). Which of course will be one the uncertainties again, but in any case the results obtained in this way are also validated.

Thank you,
Arsenis Chatzimichailidis
kiski, randolph, calf.Z and 1 others like this.
arsenis is offline   Reply With Quote

Old   February 8, 2023, 21:13
Senior Member
abdikerim kurbanaliev
Join Date: Jun 2010
Location: Kyrgyzstan, Osh
Posts: 120
Rep Power: 15
kerim is on a distinguished road
Dear Arsenis,
I know it's been a long time, but I would like to know about your progress on this issue. Could you share with your successes? Any kind of assistance is highly appreciated.
kerim is offline   Reply With Quote


boundary conditions, cyclic, mapped, openfoam

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
Getting Started with OpenFOAM wyldckat OpenFOAM 25 August 14, 2022 13:55
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 10 September 2, 2021 05:29
No pressure drop in twoPhaseEulerFoam pipe with cyclic boundaries rdbisme OpenFOAM Running, Solving & CFD 4 December 25, 2016 10:29
[mesh manipulation] Problem with using createPatchDict to define cyclic boundaries kaszt OpenFOAM Meshing & Mesh Conversion 0 April 1, 2016 21:18
Wrong wall distance with cyclic boundaries sebastian OpenFOAM Bugs 4 October 31, 2012 10:24

All times are GMT -4. The time now is 04:05.