CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Impose volume mesh vertex motion with external library

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

Like Tree1Likes
  • 1 Post By daveatstyacht

LinkBack Thread Tools Search this Thread Display Modes
Old   March 11, 2019, 11:51
Default Impose volume mesh vertex motion with external library
Senior Member
Join Date: Jul 2010
Posts: 100
Rep Power: 15
daveatstyacht is on a distinguished road
Hello all,
I am currently trying to couple an external mesh deformation library to pimpleDymFoam of v1812. The external library is in Fortran90 so I pass the components of the mesh vertices for both the volume mesh and the boundary mesh to the library. The library returns the components of the mesh vertex displacements and I can convert these 1D arrays into the vectorField type of OpenFOAM. I manage to apply the deformations of the boundaries and propagate the displacements using:

pointDisplacement.boundaryFieldRef()[iPatch] == displacementValues;

where displacementValues is a vectorField of the given patch's vertex displacements. Here I am using dynamicMotionSolverFvMesh and displacementLaplacian. In contrast, I am struggling to directly deform the volume mesh along the lines of:
Manipulation of every single mesh node inside the domain (OpenFoam)

Trying to add the code of the above post into the pimpleDymFoam solver I find that it compiles fine, but that the mesh doesn't move.

vectorField newPoints = mesh.points() + volumeDisp;

Adding mesh.update() changes nothing.

Has anyone got any experience with how to directly impose the internal mesh vertex positions?

Should I continue to try using something along the lines of: pointDisplacement.somethingHere() == volumeDisp (of type vectorField)

or is the solution further afield like adding a meshSolver type that takes as input a user input vectorField? Note, re-implementing the deformation library internally to OpenFOAM is not an option for me.

Thanks in advance for any help!
atulkjoy likes this.
daveatstyacht is offline   Reply With Quote

Old   March 13, 2019, 07:37
Atul Kumar
Join Date: Dec 2015
Location: National Centre for Combustion Research and Development
Posts: 48
Rep Power: 10
atulkjoy is on a distinguished road
Can you please put whole case/changes ??
atulkjoy is offline   Reply With Quote


mesh deformation, openfoam 1812, volume mesh

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
ill defined primitive entry 'boundary' francisco.angel OpenFOAM Pre-Processing 16 September 8, 2022 06:11
[mesh manipulation] Mesh Refinement Luiz Eduardo Bittencourt Sampaio (Sampaio) OpenFOAM Meshing & Mesh Conversion 42 January 8, 2017 12:55
[CAD formats] my stl surface is seen as just a line rcastilla OpenFOAM Meshing & Mesh Conversion 2 January 6, 2010 01:30
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 14:09
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 14:00

All times are GMT -4. The time now is 08:00.