[daveatstyacht]

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.


ie:
vectorField newPoints = mesh.points() + volumeDisp;
mesh.movePoints(newPoints);



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!
-Dave

[atulkjoy]

Can you please put whole case/changes ??