{"id":374823,"date":"2024-08-02T07:01:58","date_gmt":"2024-08-02T07:01:58","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/after-importing-the-udf-it-crashes-as-soon-as-the-calculation-is-made\/"},"modified":"2024-08-02T07:08:45","modified_gmt":"2024-08-02T07:08:45","slug":"after-importing-the-udf-it-crashes-as-soon-as-the-calculation-is-made","status":"publish","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/after-importing-the-udf-it-crashes-as-soon-as-the-calculation-is-made\/","title":{"rendered":"After importing the UDF, it crashes as soon as the calculation is made"},"content":{"rendered":"

We currently use porous media to make capillary cores, and UDFs to add capillary force to porous media so that they can be aspirated, but the code we write can be imported normally but crashes as soon as we calculate it<\/p>\n

#include “udf.h”<\/div>\n
#define PARAMETER_P 1.0 \/\/ \u66ff\u6362\u4e3a\u5b9e\u9645\u7684\u53c2\u6570\u503c<\/div>\n
#define UDS_INDEX 0<\/div>\n
#define UDM_INDEX 0<\/div>\n
#define USR_SOURCE_INDEX 0<\/div>\n
 <\/div>\n
DEFINE_ADJUST(adjust_and_store_gradient, domain)<\/div>\n
{<\/div>\n
    Thread *t;<\/div>\n
    cell_t c;<\/div>\n
    face_t f;<\/div>\n
    real grad[ND_ND]; \/\/ \u5b58\u50a8\u68af\u5ea6\u7684\u6570\u7ec4<\/div>\n
    domain = Get_Domain(2); \/\/ \u5047\u8bbedomain_ID\u4e3a2<\/div>\n
    thread_loop_c(t, domain)<\/div>\n
    {<\/div>\n
        begin_c_loop(c, t)<\/div>\n
        {<\/div>\n
            C_UDMI(c, t, UDM_INDEX) = 0.0; \/\/ \u521d\u59cb\u5316UDM<\/div>\n
            C_UDSI(c, t, UDS_INDEX) = C_VOF(c, t); \/\/ \u8bbe\u7f6eUDS\u4e3a\u5f53\u524d\u4f53\u79ef\u5206\u6570<\/div>\n
            if (C_UDSI(c, t, UDS_INDEX) > 0.0 && C_UDSI(c, t, UDS_INDEX) < 1.0)<\/div>\n
            {<\/div>\n
                F_GRAD_CELL_THREAD(c, t, &C_UDSI(c, t, UDS_INDEX), grad);<\/div>\n
                C_UDMI(c, t, UDM_INDEX) = NV_MAG(grad); \/\/ \u8ba1\u7b97\u68af\u5ea6\u7684\u6a21\u5e76\u5b58\u50a8\u5230UDM<\/div>\n
            }<\/div>\n
        }<\/div>\n
        end_c_loop(c, t)<\/div>\n
    }<\/div>\n
}<\/div>\n
 <\/div>\n
DEFINE_SOURCE(apply_momentum_source, c, t_cell, dS, eqn)<\/div>\n
{<\/div>\n
    real p = PARAMETER_P; \/\/ \u7ed9\u5b9a\u7684\u53c2\u6570p<\/div>\n
    real vol_frac = C_UDSI(c, t_cell, UDS_INDEX); \/\/ \u4f53\u79ef\u5206\u6570<\/div>\n
    real grad_mag = C_UDMI(c, t_cell, UDM_INDEX); \/\/ \u4f53\u79ef\u5206\u6570\u68af\u5ea6\u7684\u6a21<\/div>\n
    real F = 2.0 * p * vol_frac * grad_mag; \/\/ \u6839\u636e\u516c\u5f0f F = 2 * p * vol_frac * grad_mag \u8ba1\u7b97\u52a8\u91cf\u6e90\u9879<\/div>\n
    if (vol_frac > 0.0 && vol_frac < 1.0)<\/div>\n
    {<\/div>\n
        dS[2] += F; \/\/ \u5728\u6240\u6709\u65b9\u5411\u4e0a\u6dfb\u52a0\u52a8\u91cf\u6e90\u9879<\/div>\n
    }<\/div>\n
    return 0;<\/div>\n
}<\/div>\n

This is our UDF,The following issues occur\"\"<\/p>\n

Here’s the model we use\"\"<\/p>\n","protected":false},"template":"","class_list":["post-374823","topic","type-topic","status-publish","hentry","topic-tag-capillary-force","topic-tag-fluent"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t