{"id":250629,"date":"2023-02-20T16:22:42","date_gmt":"2023-02-20T16:22:42","guid":{"rendered":"\/forum\/?post_type=topic&p=250629"},"modified":"2023-02-20T16:37:28","modified_gmt":"2023-02-20T16:37:28","slug":"fluid-velocity-in-dpm_drag","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/fluid-velocity-in-dpm_drag\/","title":{"rendered":"Fluid velocity in dpm_drag"},"content":{"rendered":"

Hi,<\/p>\n

to implement a custom drag model I need to calculate the difference in velocity between fluid and particle. Getting the particle velocity is no problem however the fluid velocity is causing me trouble. I am trying to access it using the following UDF:<\/p>\n

 #include “udf.h”
 #include “surf.h”
 #include “dpm_types.h”
 #include “dpm.h”
 #include “dpm_mem.h”
 #include “mem.h”
 #include “materials.h”
 #include “threads.h”<\/p>\n

DEFINE_DPM_DRAG (test_drag, Re, tp)
{
    cell_t c;
    Thread *t0;
    real uvel;
    
    t0 = TP_CELL_THREAD(tp);
    c = TP_CELL(tp);
    
    uvel = C_U(c,t0);<\/p>\n

}<\/p>\n

When executing this UDF I am getting the following error:<\/p>\n

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

 <\/p>\n

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

 <\/p>\n

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

 <\/p>\n

999999: mpt_accept: error: accept failed: No error<\/p>\n

 <\/p>\n

Can you help me, what am I doing wrong?<\/p>\n

 <\/p>\n

Thank you!<\/p>\n

Chris<\/p>\n","protected":false},"template":"","class_list":["post-250629","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