search

precice / openfoam-adapter

OpenFOAM-preCICE adapter
https://precice.org/adapter-openfoam-overview.html
GNU General Public License v3.0
135 stars 80 forks source link

FSI: Add an option for stress mapping #41

Closed MakisH closed 4 years ago

MakisH commented 5 years ago

At the moment, we map forces. However, it would be very nice to also be able to map stresses, which would be easy to implement, having already implemented Force.C.

A limitation here is that we cannot use stresses in CalculiX at the moment. So, a prerequisite for this would be to either:

A. Implement stresses reading also in the calculix-adapter, B. Add connectivity information in the openfoam-adapter, so that we use a preCICE action to convert the stresses to forces (see also #2 ). C. Use another structure solver that already supports reading stresses (is there one already?).

This is relevant to #7 .

derekrisseeuw commented 5 years ago

In my branch FSI_devStress I added the stress class. I plan to merge this branch soon with the FSI_dev branch on my fork.

You will find that the code is very similar to the force class with a lot of code duplication. The difference between the classes is that the force is multiplied by the surface face areas to get the conservative quantity, while this is not done for the stresses.

If I am correct, deal-ii can handle stresses?

davidscn commented 5 years ago

In general, you can choose in the deal.ii code, whether you want to calculate with stresses or with forces. But consider, that the deal.ii-precice version is currently a prototype. It could be useful for testing of course.

uekerman commented 5 years ago

Just for completion: Option D: Use connectivity in the structure solver to apply the action.

MakisH commented 4 years ago

This was closed by #125.