Hi,
Hope you are well.
I would like to ask you how can I perform a time average of the wall shear stress on the vertical wall of a square cavity i.e from the tutorial examples. I know the tangential stress can be called from post_util.f90 and pass the information in the cs_user_extra_operations.f90 in order to print it in a file but I do not know how the indexing corresponds to each component.
Thank you in advance,
Costas
Wall_shear_stress_from_post_util
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
-
- Posts: 284
- Joined: Fri Dec 04, 2015 1:42 pm
Re: Wall_shear_stress_from_post_util
Hello,
In the file cs_user_extra_operations-global_efforts.f90 there is an exmple how to compute global efforts on a subset of faces. Did you check it?
Regards,
Luciano
In the file cs_user_extra_operations-global_efforts.f90 there is an exmple how to compute global efforts on a subset of faces. Did you check it?
Regards,
Luciano
Re: Wall_shear_stress_from_post_util
Dear Luciano,
Thanks for you reply. I will check it though I do not know how valid is that approach if you use standard k-epsilon.
Regards,
Costas
Thanks for you reply. I will check it though I do not know how valid is that approach if you use standard k-epsilon.
Regards,
Costas
-
- Posts: 4077
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Wall_shear_stress_from_post_util
Hello,
The efforts in the "global_efforts" example are based on accessing precomputed efforts/stresses, that are computed based on the different term contributions during a time step, so it should be correct with any turbulent model.
Note that when computation of efforts/stresses is activated (for example when requesting their postprocessing output in the GUI), the corresponding values are stored in a field, and you can simply use the field with standard time moment definitions (in the user subroutines, as this may not be accessible yet for boundary time moments).
Best regards,
Yvan
The efforts in the "global_efforts" example are based on accessing precomputed efforts/stresses, that are computed based on the different term contributions during a time step, so it should be correct with any turbulent model.
Note that when computation of efforts/stresses is activated (for example when requesting their postprocessing output in the GUI), the corresponding values are stored in a field, and you can simply use the field with standard time moment definitions (in the user subroutines, as this may not be accessible yet for boundary time moments).
Best regards,
Yvan
Re: Wall_shear_stress_from_post_util
Dear Yvan,
Thank you for your reply.
I thought of finding a way to call directly the friction velocities from clptur.f90 i.e rho uk*uk and compute the wall shear stress. I supposed this is from where the code compute the efforts.
Though, is it possible to make a call for the variable uk in that routine and print it out in cs_user_extra_operations.f90 or do I have to declare a user field in the cs_user_parameters.f90?
Thanks,
Costas
Thank you for your reply.
I thought of finding a way to call directly the friction velocities from clptur.f90 i.e rho uk*uk and compute the wall shear stress. I supposed this is from where the code compute the efforts.
Though, is it possible to make a call for the variable uk in that routine and print it out in cs_user_extra_operations.f90 or do I have to declare a user field in the cs_user_parameters.f90?
Thanks,
Costas