{"id":360113,"date":"2024-04-02T18:23:10","date_gmt":"2024-04-02T18:23:10","guid":{"rendered":"\/forum\/forums\/topic\/error-in-define_pb_break_up_rate_freq-udf-example-from-manual\/"},"modified":"2024-04-02T18:23:10","modified_gmt":"2024-04-02T18:23:10","slug":"error-in-define_pb_break_up_rate_freq-udf-example-from-manual","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/error-in-define_pb_break_up_rate_freq-udf-example-from-manual\/","title":{"rendered":"Error in DEFINE_PB_BREAK_UP_RATE_FREQ UDF example from manual"},"content":{"rendered":"

Hello!<\/p>\n

 <\/p>\n

I’m trying to use the break-up kernel from fluent manual:<\/p>\n

#include “udf.h”<\/div>\n
#include “sg_pb.h”<\/div>\n
#include “sg_mphase.h”<\/div>\n
 <\/div>\n
DEFINE_PB_BREAK_UP_RATE_FREQ(break_up_freq_tav, cell, thread, d_1)<\/div>\n
{<\/div>\n
real epsi, alpha, f1, f2, rho_d;<\/div>\n
real C1 = 0.00481, C2 = 0.08, sigma = 0.07;<\/div>\n
Thread *tm = THREAD_SUPER_THREAD(thread);\/*passed thread is phase*\/<\/div>\n
epsi = C_D(cell, tm);<\/div>\n
alpha = C_VOF(cell, thread);<\/div>\n
rho_d = C_R(cell, thread);<\/div>\n
f1 = pow(epsi, 1. \/ 3.) \/ ((1. + epsi)*pow(d_1, 2. \/ 3.));<\/div>\n
f2 = -(C2*sigma*(1. + alpha)*(1. + alpha)) \/ (rho_d*pow(epsi, 2. \/ 3.)*pow(d_1, 5. \/ 3.));<\/div>\n
return C1 * f1*exp(f2);<\/div>\n
}<\/div>\n
 <\/div>\n
I get the following error during the solution: <\/div>\n
\n

Node 0: Process 1660: Received signal SIGSEGV.<\/p>\n

==============================================================================<\/p>\n

==============================================================================<\/p>\n

Node 1: Process 22400: Received signal SIGSEGV.<\/p>\n

==============================================================================<\/p>\n

===================================================================================<\/p>\n

= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES<\/p>\n

= RANK 0 PID 1660 RUNNING AT DESKTOP-CF17HIG<\/p>\n

= EXIT STATUS: -1 (ffffffff)<\/p>\n<\/div>\n

 <\/div>\n
 <\/div>\n
 <\/div>\n
The error stops if I change epsi to a constant like “epsi = 1.;” Does anyone knows how to solve it?<\/div>\n","protected":false},"template":"","class_list":["post-360113","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["74761","15180"],"_bbp_author_ip":["23.206.193.146"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["235"],"_bbp_topic_status":["unanswered"],"_bbp_topic_id":["360113"],"_bbp_forum_id":["27792"],"_bbp_engagement":["15180","74761"],"_bbp_voice_count":["2"],"_bbp_reply_count":["1"],"_bbp_last_reply_id":["360420"],"_bbp_last_active_id":["360420"],"_bbp_last_active_time":["2024-04-04 07:22:40"]},"test":"jussan-jposgrad-ufsc-br"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/360113","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\/360113\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=360113"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}