Caculix diverged when solving multiphase FSI problem with compressibleInterDyMFoam

Official adapters and tutorials
, ,
1

diverge_solide_1e-8ce
diverge_solide_1e-8ce1269×1378 48.2 KB

1
11497×970 159 KB

2
21234×786 184 KB

Hello, everyone!
I’m using precice and caculix for the first time. I modified the vertical-flap example by eliminating the incoming flow and adding a bubble near the elastic plate,hoping that bubble collapse would cause displacement and deformation of the plate.
But the solid solution diverges at coupling. There is no displacement of the plate and only a small deformation in the result.

I wonder if anyone understand why the solid solver diverges? Thank you very much!
I uploaded pictures of caculix divergence and geometric picturesprecice-config.zip (996 Bytes)
solid-calculix.zip (10.8 KB)

I also uploaded the entire case.
I hope you can help me see how to fix it

2

fluid-openfoam.zip (165.9 KB)

3

log.openfoam.log (52.6 KB)
This is the openfoam log file (omit the middle output otherwise the file is too big)

4

@zyx is the solid supposed to compress next to the bubble? Are the material parameters for the solid correct?

If everything is correct, I could imagine that the adapter computes and couples the wrong force (unreasonably high) on the region of the interface touched by the bubble. Could you check the force acting on the plate?

You could also use the preCICE exports: Export configuration | preCICE - The Coupling Library

5

Thank you very much for your reply! I want to see deflection in a two-dimensional plate rather than distortion. I set the plate as linear elastic material with an elastic modulus of 20MPa.

Now that I’ve eliminated the subcycle, it can count more time steps. But all you see is the distortion of the plate.I uploaded the initial and divergent calculations.

1
11924×1030 335 KB

2
21920×1019 171 KB

I also think that the force on the board is acting too much, maybe the adapter calculated and coupled the wrong force.

I was wondering if you have any good advices about the adapter calculation and coupling error?
I appreciate very much !

6

Let’s take a step back: are you running the official adapter, or a modified version? There are a few discussions already regarding this: Multiphase FSI case with the OpenFOAM adapter - #2 by Makis

See also the multiphase tag: Topics tagged multiphase

7

Thank you very much for your reply!
I’m using openfoam-v2206 and the official openfoam-adapter.
My fluid solver is compressibleInterFoam (but with the modification of removing T),which may be different from interFoam in terms of the reason for the error in calculating the force.
I’m still confused about the deformation of the interface mesh.I need further checks.
Any suggestions will be appreciated! :blush:

8

Have a look at the links I posted above, I understand that other users have faced similar issues.