Dear Code-Saturne team,
The restart of the calculation in ALE module does not appear to be continuous. The possible reason could be the values of force, velocity, displacement and acceleration of structure are not saved and used from previous calculations. Mainly the displacement is not continuous and the velocity is taken as zero. In an implicit calculations the displacement at a next time step uses the values of displacement and velocity at current timestep. We have some snapshots of the results with this issue. Please find .pdf file attached.
Thanks a lot.
Vilas
continuation of the fluid-structure interaction calculations
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
continuation of the fluid-structure interaction calculations
- Attachments
-
- results.pdf
- displacement, velocity and force plot at a restart of calculation.
- (11.4 KiB) Downloaded 224 times
-
- Posts: 4208
- Joined: Mon Feb 20, 2012 3:25 pm
Re: continuation of the fluid-structure interaction calculat
Hello,
This issue has been reported recently on this forum (http://code-saturne.org/forum/viewtopic.php?f=2&t=1520).
I just fixed it today. The issue was due to initialization of some structure values (to zero) after reading restart files. Using user subroutines (usstr1), you can work around this, by "protecting" the initialization, for example (based on the current example) :
With the GUI, there was a bug.
I just released 3.0.5, which fixes this bug with the GUI (and a few others, see the NEWS file), as well as the user subroutine example.
Regards,
Yvan
This issue has been reported recently on this forum (http://code-saturne.org/forum/viewtopic.php?f=2&t=1520).
I just fixed it today. The issue was due to initialization of some structure values (to zero) after reading restart files. Using user subroutines (usstr1), you can work around this, by "protecting" the initialization, for example (based on the current example) :
Code: Select all
if (ntpabs.le.1) then ! Avoid resetting in case of restart
xstr0(2,1) = 2.d0
xstreq(2,1) = 1.d0
vstr0(3,2) =-0.5d0
endif
I just released 3.0.5, which fixes this bug with the GUI (and a few others, see the NEWS file), as well as the user subroutine example.
Regards,
Yvan