Provide initial data for perpendicular flap case #540
There are no files selected for viewing
| Original file line number | Diff line number | Diff line change |
|---|---|---|
|
|
@@ -46,8 +46,8 @@ | |
| <time-window-size value="0.01" /> | ||
| <max-time value="5" /> | ||
| <participants first="Fluid" second="Solid" /> | ||
| <exchange data="Force" mesh="Solid-Mesh" from="Fluid" to="Solid" initialize="true" /> | ||
| <exchange data="Displacement" mesh="Solid-Mesh" from="Solid" to="Fluid" initialize="true" /> | ||
|
There was a problem hiding this comment. Why initializing the displacement? The solid is initially at rest, no? Maybe the fake solid does not mimic this behavior? There was a problem hiding this comment. The displacement is not relevant for the fake solid solver (it ignores data it reads). But I'm also suggesting to generally initialize here. Generally, I think explicitly initializing with zero is much safer than implicitly initializing with zero like we do it currently. Therefore, I would prefer to initialize all exchanged data (potentially with zero). See precice/precice#2033 There was a problem hiding this comment. The problem with this that will bite us is that then the first mapped and communicated data is zero. I am afraid that users will be very confused and wonder why we waste the resources. There was a problem hiding this comment. In the end this is a question of priorities: If we want the default case to offer maximum performance at the risk of being less flexible while maintaining accuracy then we should not initialize. If we want to allow for customization of the case to some degree without introducing errors we should initialize. I'm leaning towards prioritizing accuracy and flexibility over performance. I would also argue that potentially wasting resources once in the beginning of the simulation is not as bad as introducing an error at the beginning from which you cannot recover. There was a problem hiding this comment. I don't agree on the flexibility argument. Our current solution is flexible, user can initialize data if they want. In the end, the tutorial setup fixes whether non-zero data needs to be initialized or not. For the perpendicular flap, this is not the case. If a user wants to alter the case, they have to alter the configuration as well. Besides performance, my main concern is that users will wonder why initial data is explicitly initialized if it is zero anyway. There was a problem hiding this comment. Then let's maybe phrase what I mean with flexibility differently: Our current scenario requires you to do (and remember) more additional steps when you want to use non-zero data because the default case assumes zero data. We can avoid having to do these additional steps, but this costs us performance and might look unnecessary to the users. I think the key question is how we want to resolve precice/precice#2033. |
||
| <max-iterations value="50" /> | ||
| <relative-convergence-measure limit="5e-3" data="Displacement" mesh="Solid-Mesh" /> | ||
| <relative-convergence-measure limit="5e-3" data="Force" mesh="Solid-Mesh" /> | ||
|
|
||
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Initializing the force could be tricky. How would you do it for a non-fake fluid solver? I guess an example could help here, e.g. OpenFOAM.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
OpenFOAM would still initialize with zero forces. That's also working out of the box with this configuration.