Surface element set self contact

I’m attempting to couple CaculiX with OpenFOAM via preCICE to simulate 2-way FSI for a heart valve. As the fluid pressure forces the valve open, contact is handled by a single contact pair between the downstream surface of the valve leaflets and the blood vessel wall. However, when the pressure gradient flips to close the valve the upstream surface of the valve will come in contact with itself. How can I properly model this contact? Is it as simple as defining a contact pair between the upstream surface and itself?

I’m working with CalculiX 2.16, in case that matters.

You can define only one surface set including both sides of the contact and make a contact pair calling the same surface set for master and slave, see the example

2021-02-04 15_52_33-rubber_bellow.inp - SciTE

1 Like

Hi mtree,

I am interested in such a case. Have you managed to simulate heart valve and newtonian blood behaviors in a 2 ways FSI simulation with contact thanks to OpenFOAM/CalculiX/preCICE?

Any leakage or instability issues?

Thank you in advance,