Closed MakisH closed 3 years ago
This is due to a new check introduced in https://github.com/precice/precice/pull/787
I rebuild preCICE yesterday in the evening and it worked with the deal.II adapter (in 2D). I think this is a problem of the initial ramp, where not very much is happening and we just write the same data due to nothing is happening in the first time step. I get also a warning in the beginning, but not an error. In 3d, it might be heavier, since the out of plane direction writes always zero data.
So, what do we do then?
Remains the error, when you remove the ramp?
Options are 1) raise a warning instead of an error in preCICE 2) remove the 3d deal.II examples 3) modify the ramp
Could you please try again with the release candidate?
This error I get with the 2D case, so I don't understand how removing the 3D examples would help.
This error I get with the 2D case, so I don't understand how removing the 3D examples would help.
My bad, i will check it again, though I'm working with the parabolic profile and groovy.
This looks really seriously, I can reproduce it. It fails in the first iteration (for some strange reason not with groovyBC, but with the original one): The problem is (my best guess), that the displacement doesn't converge according to the prescribed limits, but the initial data (zero) is passed to the solid solver again due to reloading.
Quick workaround/fix: Switch to a serial-implicit
scheme
The problem is the missing data entry in the acceleration
<acceleration:IQN-ILS>
<data name="Displacement" mesh="Solid_mesh"/>
+ <data name="Stress" mesh="Fluid-Mesh-Centers"/>
<preconditioner type="residual-sum"/>
<filter type="QR1" limit="1e-6"/>
<initial-relaxation value="0.1"/>
<max-used-iterations value="50"/>
<time-windows-reused value="10"/>
</acceleration:IQN-ILS>
However, then, dealii destroys the SolverInterface early in the second timestep
---[precice] Implicitly finalizing in destructor
I guess due to convergence problems?
All in all, I think the preCICE error message here is meaningful.
Hm, I'm currently running parallel-implicit
without acceleration in both data sets, should not cause problems.
So yes, looking at the solver data, the error is right, deal.II receives all the time the same data leading to all the time the same output..
I guess due to convergence problems?
This should happend only, when the simulation is terminated abnormal.
We used to have a check for this:
But we removed it at some point as it does not always make sense.
For the nonlinear dealii solver it also crashes with serial-implicit
for me (have not tested the linear one yet).
Yes this is related to #77 . It runs properly with a better config.
In your case, the problem should be different from this though.
linear works
Suggestion: let's remove this tutorial from the documentation until we have #77 ?
linear works
parallel or serial?
I think we can not resolve this in a few hours without looking more in detail what's going wrong, so yes.
parallel or serial?
with the linear solver both coupling schemes work
I ran all the cases (flap/turek, 2D/3D, and linear/nonlinear) in serial. Here are some results: https://github.com/precice/precice/pull/835#issuecomment-661824557
Resolved in https://github.com/precice/tutorials/pull/116
With preCICE v2.1.0 (release candidate, 295d27), either something changed in behavior, or some check became stricter:
The FSI/flapPerp_2D/OpenFOAM-deal.II is not affected.