I want to implement Peskin’s classical IBM using openfoam and my own solid solver. I have implemented it in serial, however the code is not very optimum and works only in serial. I plan to use preCICE for this purpose.
I have been going through the documentation of Precice and plan to use that for the coupling.
I have a couple of questions.
- In precice the data is exchanged at the interface of the two meshes, but in my case I would like to exchange data at all the vertices of the solid mesh. Those on the surface as well as internal. Is it possible?
- Is it possible to implement a custom Dirac delta function kernel to do the mapping between the solid and fluid mesh?
Hi,
Yes, preCICE handles volume coupling. Keep in mind that this often requires tweaking your mappings.
If you want to compute your own mapping, then a good way is to use Direct access to received meshes | preCICE - The Coupling Library.
You receive the mesh on the solver you want to perform the mapping on and query preCICE for the vertex coordinates. Then you can build your own mapping and directly map to/from your solver.
Best regards
Frédéric
Thank you for your reply. I went through the documentation on direct access. It seems a region of interest needs to be defined beforehand. However, in my case, the solid components deform and move all over the domain. Will it still be possible to use this approach?
Difficult to say. Some visual reference would help to understand this. A drawing, figure, video or paper would be great.
We collect information on how to handle moving meshes, which may answer some of your questions:
Thank you. I’ll go through the provided resource. Meanwhile I am attaching a animation that is generated using my serially coupled solver in openfoam. I my research I am planning to simulate such particles over arrays of cilia in a fluid domain.
Animation