CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[ImmersedBoundary] About the moving immersed boundary tutorial: icoDyMIbFoam+movingCylinderInChannelIco

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

Like Tree9Likes

LinkBack Thread Tools Search this Thread Display Modes
Old   September 14, 2016, 05:29
New Member
Daniel Brunner
Join Date: Sep 2015
Posts: 8
Rep Power: 7
DaniB1989 is on a distinguished road

Any motion will be described by the solidBodyMotionFunctions-class. This is the same as one usually uses for moving meshes.
If you have the desire to add an additional function you can do that by copying an existing function e.g. translation and modify it as you please. After you did that, add the path into the files (Make folder) and compile the entire class. However it's not the best solution since you're changing the source code however it works.

Best wishes
ShinyRiver likes this.
DaniB1989 is offline   Reply With Quote

Old   July 4, 2017, 13:40
Default Make immersed boundarys compressible
New Member
Robert Grandl
Join Date: Aug 2015
Posts: 16
Rep Power: 7
RobertG is on a distinguished road
Hello everybody,
I am using ImmersedBounderys in foamExtend-3.2 and I want to adapt the compressible rhoPisoFoam-Solver to calculate use ImmersedBoundaries.
So far, I changed the "immersedBoundaryAdjustPhi.H" from "fvc :: makeRelative (phi, U);" to "fvc :: makeRelative (phi, rho, U);"

Using rhoPisoFoam as a basement, I've set up the solver as the solver_coockbook suggested.

But every time I try to run it, an error occurs, caused by a division by zero.

Responsible is basically cellIbMaskExt at 0 positions (solid body). The values for rho, h ans psi are at the same position also 0.
Unfortunately, rho and h are poor implemented. So that correctBoundaryConditions() does not help at all.

I hope someone can help me to make the solver run.

I read, that the used implementation of immersedBoundaries, are of "Discrete Forcing Approach" and uses "direct imposition of boundary conditions". Are there more infomations about the used algorithm?


Last edited by RobertG; July 5, 2017 at 17:56.
RobertG is offline   Reply With Quote

Old   October 8, 2018, 05:39
New Member
Alessandra Bavo
Join Date: Aug 2015
Location: Belgium
Posts: 2
Rep Power: 0
abavo is on a distinguished road
Hi all,

I am using foam-extend 4.0 and currently trying to run the tutorial movingCylinderInChannelIco. The case seems to run, but as soon as I want to display the results in paraFoam I get this error:

size of field refValue (96) is not the same size as the patch (0)
on patch ibCylinder of field U in file "/home/username/foam/foam-extend-4.0/tutorials/immersedBoundary/movingCylinderInChannelIco/0.2/U"

By inspecting the velocity files, in 0/U the refValue is defined as uniform, while in 0.2/U refvalue becomes nonUniform and defined with a List of vectors. Is the error related to this difference? How can I fix it? Is it a paraFoam issue or there is something wrong with the case? I did set up a simple case very similar to this tutorial, and modified the U file (e.g. by specifying a nonUniform List also in 0/U) but the error is the same.

Thank you in advance for your feedback
abavo is offline   Reply With Quote

Old   July 16, 2020, 08:45
Default foam-extend-4.0/4.1 tutorial blockMesh
New Member
Max Zwicker
Join Date: Apr 2019
Posts: 2
Rep Power: 0
mfrz is on a distinguished road
Dear Foamers,

I am running a bunch of the Immersed Boundary method tutorials, specifically the movingCylinderInChannelIco tut. I played around extensively and a lot of stuff works quite well but I cannot refine the blockmesh. When changing the mesh in the blockmeshdict the solvers and particularly the meshInitialisation complain. I attached the error below. It's seems to be the same for both 4.0 and 4.1. Can someone help with it?

Thanks in advance,


Create dynamic mesh for time = 0

Selecting dynamicFvMesh immersedBoundarySolidBodyMotionFvMesh

Error in point ordering: mixed used and unused points at the end of point list.
Last used point: 0 (-1 -0.5 -0.04)
First unused point: 1 (-0.982857 -0.5 -0.04)
and point 1 (-0.982857 -0.5 -0.04) is used by a live face.
Face 0 4(1 177 22353 22177) with points 4((-0.982857 -0.5 -0.04) (-0.982857 -0.492 -0.04) (-0.982857 -0.492 0.04) (-0.982857 -0.5 0.04))

From function void polyMesh::initMesh()
in file meshes/polyMesh/polyMeshInitMesh.C at line 166.

FOAM aborting
mfrz is offline   Reply With Quote


immersed boundary method

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
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 01:44
Domain Imbalance HMR CFX 5 October 10, 2016 05:57
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 04:39
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 10:53
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00

All times are GMT -4. The time now is 20:19.