{"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":[],"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