Shear rate and pulsatile flow

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
Salvador Jimenez

Shear rate and pulsatile flow

Post by Salvador Jimenez »

Hi everyone,
I am working on a flow in pipe-like
structure (artery). 
I would like to visualize the shear stress. Anyone knows how to do
this using Satune/Salome?
Thanks
Yvan Fournier

Re: Shear rate and pulsatile flow

Post by Yvan Fournier »

Hello,
You may use the usvpst.f90 user subroutine to add additional user-computed fields to the visualizable output, but you need to extract the shear stress first. This would mean using cs_user_modules.f90 (for version 2.1) or the rtuser array (version 2.0) to save the field from where it may be computed (the tricky part), then outputting it in usvpst.f90.
You may check examples in usproj.f90, but I believe you will need to navigate the code's non-user subroutines to extract the shear where computed (I am not sure where, but others than me may help here), and add one of those to your user subroutines to do this.
You may also try to recompute the shear stress from available variables (the user and theory guides may help here).
Also, if you are using a Reynolds-stress (Rij), turbulence model, the components of the Reynolds stresses are available by default for postprocessing, though that might not be what you are looking for.
Best regards,
  Yvan
Salvador Jimenez

Re: Shear rate and pulsatile flow

Post by Salvador Jimenez »

Hi Yvan,
thanks for your answer, I am actually interested for the laminar case. I'm still trying to figure it out.
I will have a look on the cs_user_modules.f90
 
Best regards,
Salvador
Post Reply