code_saturne User's Forum

Dear Code Saturne Support Team,

My name is Valeria and I am a student in FH Aachen University. I have a geometry file to be computed, however I am not sure if it is possible to be done with Code Saturne.
It is a medical device, in a shape of cylinder with a half-schere extrusion on a side going around the cylinder. I have included screenshots of the geometry file for a better understanding.
The problem is that half-scherical extrusion must pulsatory change the form, ejecting the blood of the volume to the outside. There are three pictures named "State1 -> State2 -> State3", which are illustrating the way the geometry is changing. Red arrows illustrate the force from the outside that compressing the geometry and the blue arrows are approximation of the blood flow. The inlet is on the bollom of the cylinder and it is equiped with a valve which is closed once the spreare is pressed and the blood is ejected to the outlet. I hope it is possible to understand from the given pictures. Please feel free to ask me for additional information or clarification.
I have never dealed with such cases and would like to ask you for support.
Do you what much be the strategy for such types of geometry problems?

Thank you very much!

Best Regards,
Valeria

Hello,

The pictures are not completely clear, but if I understand correctly, you have both a deformation and a time-varying inlet.

If you are using the rotor-stator coupling, you cannot change the topology of the mesh in addition to the changes induced by the rotor-stator aspect, but you should be able to use the deforming mesh (ALE) features as long as the deforming boundary is not at the rotor-stator interface (this has never been tested before, so if you do run into issues, please forward them so we can provide patches to test, but the fixes should not be too complex).

For the time-varying inlet, you can use a time-varying inlet velocity for the fluid defining your inlet profile in cs_user_boundary_conditions.f90 (it might even be possible in the GUI, depending on how you can express the velocity).

If you want to avoid both deforming the mesh and defining varying boundary conditions, you migh also start simply with a time-varying (fluid-)mass source-term, in the area near the inlet.

Regards,

Yvan