{"id":367502,"date":"2024-05-20T06:07:35","date_gmt":"2024-05-20T06:07:35","guid":{"rendered":"\/forum\/forums\/topic\/report-the-drag-lift-force-on-particles-in-dpm\/"},"modified":"2024-05-22T03:00:58","modified_gmt":"2024-05-22T03:00:58","slug":"report-the-drag-lift-force-on-particles-in-dpm","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/report-the-drag-lift-force-on-particles-in-dpm\/","title":{"rendered":"Report the drag\/lift force on particles in DPM"},"content":{"rendered":"

Hello, everyone!<\/p>\n

I am working on a project about particle-airflow interaction in a channel in Fluent 2023. I want to find out the magnitude of the drag force and lift force on each particle. <\/p>\n

I have opened the drag force, saffman lift force in physical models of discrete phase model (DPM). Generally, i can find particle related information at Results-Particle tracks-Particle variables. However, i can not find the drag force, saffman lift force under this option.<\/p>\n

I was told certain particle-related forces such as drag and lift force are not available directly through the standard Particle Tracks Reports interface. These specific forces require customization to view and analyze by UDF. And an example of UDF has been given in the last. However, this UDF did not work in my case and software (the content in this link also not work \/forum\/forums\/topic\/how-to-report-the-drag-force-on-particles-in-dpm\/ )<\/p>\n

Does anyone know how to report this when the solution is done, or how to write an correct UDF to obtain the forces on particles?<\/p>\n

Thanks!<\/p>\n

 {<\/p>\n

#include “udf.h”<\/div>\n
 <\/div>\n
\/* Define a UDF to track forces on particles *\/<\/div>\n
DEFINE_DPM_SCALAR_UPDATE(track_forces, cell, thread, initialize, p)<\/div>\n
{<\/div>\n
    real drag_force[3];<\/div>\n
    real lift_force[3];<\/div>\n
    real virtual_mass_force[3];<\/div>\n
 <\/div>\n
    \/* Get the drag force *\/<\/div>\n
    P_DRAG(p, drag_force);<\/div>\n
 <\/div>\n
    \/* Get the lift force *\/<\/div>\n
    P_SAFFMAN_LIFT_FORCE(p, lift_force);<\/div>\n
 <\/div>\n
    \/* Get the virtual mass force *\/<\/div>\n
    P_VIRTUAL_MASS_FORCE(p, virtual_mass_force);<\/div>\n
 <\/div>\n
    \/* Store the forces in user-defined memory (UDM) *\/<\/div>\n
    P_USER_REAL(p, 0) = drag_force[0];<\/div>\n
    P_USER_REAL(p, 1) = drag_force[1];<\/div>\n
    P_USER_REAL(p, 2) = drag_force[2];<\/div>\n
    P_USER_REAL(p, 3) = lift_force[0];<\/div>\n
    P_USER_REAL(p, 4) = lift_force[1];<\/div>\n
    P_USER_REAL(p, 5) = lift_force[2];<\/div>\n
    P_USER_REAL(p, 6) = virtual_mass_force[0];<\/div>\n
    P_USER_REAL(p, 7) = virtual_mass_force[1];<\/div>\n
    P_USER_REAL(p, 8) = virtual_mass_force[2];<\/div>\n
}<\/div>\n

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