Modifying the multiple-perpendicular-flaps tutorial for turbulent flow

Official adapters and tutorials
, ,
1

Hello
I was using preCICE openfoam-dealii to modify the perpendicular flap tutorial for my case which is a turbulent flow through a 2D channel with an obstruction to simulate the effect of flexibility using RANS models not laminar. So I have incurred an error when ever I try to change the mechanical property of the flap to a rubber either I will have:

[precice] WARNING: The coupling residual equals almost zero. There is maybe something wrong in your adapter. Maybe you always write the same data or you call advance without providing new data first or you do not use available read data. Or you just converge much further than actually necessary.

or

An error occurred in line <569> of file </home/habib/build/dealii-adapter/source/linear_elasticity/linear_elasticity.cc> in function
void Linear_Elasticity::ElastoDynamics::solve() [with int dim = 2]
The violated condition was:
velocity.linfty_norm() < 1e4
Additional information:
Linear system diverged

here is my file I modified the geometry to my case and all the boundary conditions and it runs fine only in openfoam.
the run time is shorter just for checking I am running the case doe 1.5sec

I am not sure what the issue is please Help. Thank you
fluid-openfoamcase1KEflexible.zip (5.0 MB)
solid-dealii.zip (4.1 MB)
precice-config.xml (2.4 KB)

2

If there is also other recommendation other than deallii which would be better please I am happy to accept.
Thank you

3

@habaw sorry for the late reply. Please open a new topic for each new question. Answering in a very old topic will not be very visible, and will get mixed with the other person’s original question.

Did you manage to make any progress in the meantime?

Without looking at your case, I assume that the case needs some tuning. There are limitations to what you can achieve with the current setup and methods, but you could go a long way by optimizing the mesh, time step sizes, the IQN parameters, etc.

Try bringing the stand-alone cases to a state where they converge well with some pre-deformed mesh or constant load. Then try enabling the coupling.