|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Rakesh
Join Date: May 2011
Posts: 33
Rep Power: 15 ![]() |
Dear all!
Please help me to write UDF for wall shear stress equation: tau=constant*viscosity*velocitygradient. I have tried with define profile UDF using cell loop but it is resulting in acess violation error? thanks |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
|
Hi,
when you want to use gradients, you have to prevent the solver from freeing up memory by issuing the text command solve/set/expert and then answering yes to the question Keep temporary solver memory from being freed? if it doesn't work, post your UDF and I'll correct it for you. |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Rakesh
Join Date: May 2011
Posts: 33
Rep Power: 15 ![]() |
Thanks for your suggestion! The UDF which I had written for calculating wall shear is like this
![]() #include "udf.h" DEFINE_PROFILE(stressx,t,i) { cell_t c; real dvdy_r; real alpha=0.7; real dvdy=C_DVDY(c,t); real mu_a=C_MU_L(c,t); dvdy_r=alpha*dvdy; begin_c_loop(c,t) { F_PROFILE(c,t,i)=mu_a*dvdy_r ; } end_c_loop(c,t) } |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
|
It seems that you need to focus on UDF concepts more carefully.
![]() Code:
#include "udf.h" DEFINE_PROFILE(stressx,t,i) { face_t f; cell_t c; Thread *tc; real dvdy_r; real alpha=0.7; real dvdy; real mu_a; begin_f_loop(f,t) { c=F_C0(f,t); tc=THREAD_T0(t); mu_a=C_MU_L(c,tc); dvdy=C_DVDY(c,tc); dvdy_r=alpha*dvdy; F_PROFILE(f,t,i)=mu_a*dvdy_r; } end_f_loop(f,t) } Amir |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Member
Rakesh
Join Date: May 2011
Posts: 33
Rep Power: 15 ![]() |
Dear Mr. Amir,
Thanks for correcting UDF.Yes,you are right I am new to UDF writing .I will soon give u feedback on corrected UDF. with regards, Rakesh |
|
![]() |
![]() |
![]() |
![]() |
#6 |
Member
Rakesh
Join Date: May 2011
Posts: 33
Rep Power: 15 ![]() |
Dear Mr. Amir,
Sorry for late response.I had gone to some remote place.I am working with the UDF.it is working fine.Results are coming but correctness of results is to be confirmed. Thanks again! Rakesh |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Commercial meshers] Fluent3DMeshToFoam | simvun | OpenFOAM Meshing & Mesh Conversion | 50 | January 19, 2020 15:33 |
[Paraview] Wall Shear Stress | -mAx- | OpenFOAM Post-Processing | 5 | August 30, 2019 10:30 |
Wall Shear stress in 3D | nikosb | Main CFD Forum | 1 | April 15, 2010 03:52 |
Re: About wall shear stress | Mike | FLUENT | 9 | November 17, 2003 14:41 |
WALL SHEAR STRESS | Min-Hua Wang | CFX | 0 | April 3, 2002 17:40 |