{"id":362110,"date":"2024-04-15T04:11:31","date_gmt":"2024-04-15T04:11:31","guid":{"rendered":"\/forum\/forums\/topic\/need-help-with-udf-999999-mpt_accept-error-accept-failed-no-error\/"},"modified":"2024-04-15T04:11:31","modified_gmt":"2024-04-15T04:11:31","slug":"need-help-with-udf-999999-mpt_accept-error-accept-failed-no-error","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/need-help-with-udf-999999-mpt_accept-error-accept-failed-no-error\/","title":{"rendered":"Need help with UDF. “999999: mpt_accept: error: accept failed: No error.”"},"content":{"rendered":"

I want my udf to update the surface tempurature to T(f), and the mass flow inlet to m_dot at the end. When try and use the DEFINE_ADJUST macro it crashes and says “999999: mpt_accept: error: accept failed: No error.” The current code below compiles and works at printing the values (although they aren’t really what I expected). I have also attached the pseudo code of what I am trying to do.<\/p>\n

#include “udf.h”<\/div>\n
#include <math.h><\/div>\n
DEFINE_EXECUTE_AT_END(print_heat_flux)<\/div>\n
{<\/div>\n
    Domain *d;<\/div>\n
    Thread *t;<\/div>\n
    face_t f;<\/div>\n
    real Ts, heat_flux, Hv, r, Ts_old, Ts_new, m_dot, T;  \/\/ Declare all variables at the top of the block<\/div>\n
    real error;<\/div>\n
    int iter;<\/div>\n
    d = Get_Domain(1);  \/\/ Adjust the domain ID as per your case<\/div>\n
 <\/div>\n
    \/* Loop over all face threads in the domain *\/<\/div>\n
    thread_loop_f(t, d)<\/div>\n
    {<\/div>\n
        \/* Check if this is the boundary with the ID you’re interested in *\/<\/div>\n
        if (THREAD_ID(t) == 6)  \/\/ Adjust the surface ID as per your case<\/div>\n
        {<\/div>\n
            \/* Loop over all faces in this thread *\/<\/div>\n
            begin_f_loop(f, t)<\/div>\n
            {<\/div>\n
                Ts = F_T(f, t);           \/\/ Get the temperature at face f<\/div>\n
                heat_flux = BOUNDARY_HEAT_FLUX(f, t); \/\/ Calculate the heat flux<\/div>\n
                printf(“Heat flux on face %d: %g W\/m^2\\n”, f, heat_flux);<\/div>\n
Ts_old = Ts;<\/div>\n
iter = 0;<\/div>\n
error = 1;<\/div>\n
while (error>1e-6 && iter <10)<\/div>\n
{<\/div>\n
Hv = 2033631 + 310000 + 1500;<\/div>\n
r = heat_flux\/(Hv*900);<\/div>\n
Ts_new = 20557.188 \/ (8.314 * (log(0.01104) – log(r)));<\/div>\n
error=fabs(Ts_new-Ts_old);<\/div>\n
iter++;<\/div>\n
}<\/div>\n
T(f)=Ts+0.1*(Ts_new-Ts);<\/div>\n
r = heat_flux\/(Hv*900);<\/div>\n
printf(“Regression Rate: %g\\n”, r);<\/div>\n
m_dot=900*r;<\/div>\n
printf(“M_dot: %g\\n”, m_dot);<\/div>\n
            }<\/div>\n
            end_f_loop(f, t)<\/div>\n
        }<\/div>\n
    }<\/div>\n
}<\/div>\n

 <\/p>\n

\"\"<\/div>\n","protected":false},"template":"","class_list":["post-362110","topic","type-topic","status-closed","hentry","topic-tag-2d","topic-tag-fluent","topic-tag-udf"],"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