|
[Sponsors] |
July 14, 2016, 13:21 |
mesh Motion
|
#1 |
New Member
alice
Join Date: Jul 2016
Posts: 1
Rep Power: 0 |
Hi guys!
I would like to extract the shear stress of a precise region of my fluid domain at every time step of the simulation run on ANSYS Fluent. I'm used to do it defining the region of interest with an iso-surface or iso-clip and use the automatic export command to write the values on an ASCII file. In this particular simulation the mesh is moving (the motion is imposed by an UDF grid motion). So the iso-sufrace or iso-clip do not follow the movement. Does anybody know how to extract the values of stresses for moving nodes in the mesh? Is it possible to use a pointer (or something similar) to point a cell of the domain? Thanks in advance for the help! |
|
July 15, 2016, 03:12 |
|
#2 |
Member
Davoud Malekian
Join Date: Jan 2016
Posts: 53
Rep Power: 10 |
hi
i dont know if your region is a surface or not but u can use the udf below to find the shear stress on a particular surface: DEFINE_EXECUTE_AT_END(wall_shear_stress) { .............. BOUNDARY_FACE_GEOMETRY(f,tf,A,ds,es,A_by_es,dr0); C_UDSI(c, t, 0) = sqrt(NV_MAG(F_STORAGE_R_N3V(f, tf, SV_WALL_SHEAR)) / (NV_MAG(A) * C_R(c, t))); } C_UDSI(c, t, 0) = sqrt(NV_MAG(F_STORAGE_R_N3V(f, tf, SV_WALL_SHEAR)) / (NV_MAG(A) * C_R(c, t))); this part calculate the shear stress on the face element of your face thread and save them in the center of the cell that the face element is already a part of it. |
|
July 15, 2016, 09:26 |
|
#3 |
New Member
Prikane
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
Hello a.apostoli,
Since you are using UDF for dynamic mesh, somewhere in your UDF you should put the following part of code: begin_f_loop(f, tf) { f_node_loop (f, tf, n) { node_p = F_NODE(f, tf, n) ; X[f]=(NODE_X(node_p)) ; Y[f]=(NODE_Y(node_p)) ; } sumshear += NV_MAG(F_STORAGE_R_N3V(f,tf, SV_WALL_SHEAR)) ; } end_f_loop(f, tf) After you do this, you can extract ''sumshear'' variable in text file (e.g sumshear.txt) which you will create and locate in your Fluent folder, where are your case and data file located. For extraction, you can use the following code: fx=fopen("sumshear.txt","w+") ; fprintf(fx,"%E",sumshear) ; fclose(fx) ; Hope this will help you, Prikane |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
how to set periodic boundary conditions | Ganesh | FLUENT | 15 | November 18, 2020 06:09 |
Star CCM Overset Mesh Error (Rotating Turbine) | thezack | Siemens | 7 | October 12, 2016 11:14 |
Mesh motion with Translation & Rotation | Doginal | CFX | 2 | January 12, 2014 06:21 |
3D Hybrid Mesh Errors | DarrenC | ANSYS Meshing & Geometry | 11 | August 5, 2013 06:42 |
Dynamic moving mesh | Pei-Ying Hsieh (Hsieh) | OpenFOAM Running, Solving & CFD | 64 | June 7, 2012 10:04 |