{"id":318360,"date":"2023-11-19T21:55:53","date_gmt":"2023-11-19T21:55:53","guid":{"rendered":"\/forum\/forums\/topic\/error-with-mesh-deformation\/"},"modified":"2023-11-19T21:55:53","modified_gmt":"2023-11-19T21:55:53","slug":"error-with-mesh-deformation","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/error-with-mesh-deformation\/","title":{"rendered":"Error with Mesh Deformation"},"content":{"rendered":"

Hi,<\/p>\n

I’m new to Ansys Fluent and am trying to learn how to deform a mesh. Specifically, I have a chamber of the heart that I need to contract with each heartbeat. I’m trying to start simple by just having the mesh undergo deformations in all 3 axes using a sinusoidal function (see the code below). The code compiles, loads, and runs without issue but when I try to Preview Mesh Motion, there is absolutely no change to the mesh at all and I’m stumped, can anyone help me out? <\/p>\n

UDF:<\/p>\n

#include “udf.h”<\/div>\n
 <\/div>\n
DEFINE_GRID_MOTION(deform_mesh, domain, dt, time, dtime)<\/div>\n
{<\/div>\n
    real amplitude = 10.0;<\/div>\n
    real frequency = 60;<\/div>\n
    real displacement_x = amplitude * sin(2.0 * M_PI * frequency * time);<\/div>\n
    real displacement_y = amplitude * cos(2.0 * M_PI * frequency * time);<\/div>\n
    real displacement_z = amplitude * sin(2.0 * M_PI * frequency * time);<\/div>\n
 <\/div>\n
    Thread *t;<\/div>\n
    Node *v;<\/div>\n
 <\/div>\n
    \/* Loop over all nodes in the domain *\/<\/div>\n
    thread_loop_nodes(t, domain)<\/div>\n
    {<\/div>\n
        begin_node_loop(v, t)<\/div>\n
        {<\/div>\n
            real coord[ND_ND];<\/div>\n
            NV_V(coord, =, NODE_COORD(v));<\/div>\n
 <\/div>\n
            coord[0] += displacement_x;<\/div>\n
            coord[1] += displacement_y;<\/div>\n
            coord[2] += displacement_z;<\/div>\n
 <\/div>\n
            NV_S(NODE_COORD(v), =, coord);<\/div>\n
        }<\/div>\n
        end_node_loop(v, t)<\/div>\n
    }<\/div>\n
}<\/div>\n
 <\/div>\n
I imported this UDF and assigned it to the interior zone of the heart (I also tried using it with the wall, just to make sure that wasn’t my issue). From there, I tried to Preview and have tried multiple time steps, frequencies, amplitudes, etc to no avail. I’ve also tried making the mesh more\/less refined. Here is an example of what I see when I run the Preview, there is no change in any of the numbers except for the time: <\/div>\n
 <\/div>\n
\n

Updating mesh to time 4.67700e+03 (step = 00180)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.70700e+03 (step = 00181)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.73700e+03 (step = 00182)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.76700e+03 (step = 00183)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.79700e+03 (step = 00184)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.82700e+03 (step = 00185)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.85700e+03 (step = 00186)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.88700e+03 (step = 00187)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.91700e+03 (step = 00188)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n

Updating mesh to time 4.94700e+03 (step = 00189)…<\/p>\n

Dynamic Mesh Statistics:<\/p>\n

Minimum Volume = 5.44652e-15<\/p>\n

Maximum Volume = 3.02548e-08<\/p>\n

Maximum Cell Skew = 9.76395e-01 (cell zone 3)<\/p>\n

Warning: maximum cell skewness exceeds 0.95.<\/p>\n

Minimum Orthogonal Quality = 2.36051e-02 (cell zone 3)<\/p>\n

Warning: minimum orthogonal quality less than 0.05.<\/p>\n

Maximum Face Skew = 8.78838e-01 (face zone 1)<\/p>\n

 <\/p>\n

done.<\/p>\n<\/div>\n","protected":false},"template":"","class_list":["post-318360","topic","type-topic","status-closed","hentry","topic-tag-adaptivemesh-2","topic-tag-fluent-cfd-ansys-1","topic-tag-dynamic-mesh","topic-tag-mesh","topic-tag-udf"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["302528","13659"],"_bbp_author_ip":["23.206.193.146"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["396"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["318360"],"_bbp_forum_id":["27792"],"_bbp_engagement":["13659","302528"],"_bbp_voice_count":["2"],"_bbp_reply_count":["1"],"_bbp_last_reply_id":["318496"],"_bbp_last_active_id":["318496"],"_bbp_last_active_time":["2023-11-20 13:20:56"]},"test":"attanay1msu-edu"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/318360","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/318360\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=318360"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}