[Jazz]

Hey everyone,

I am modelling a solidification problem where we are tracking how inclusions move in a caster. They are represented by inert particles and I needed a way by which I can make a particle reach a wall/front as close as possible so that the distance between the particle and the wall is in nanometer scale. I am using*steady particle tracking and my distance comes out to be of order of 10^(-5)m. Even after changing the step length factor the distance comes out to be of the same order.

[Manathan]

Try to decrease the tolerance in tracking parameters in DPM model (numerics if I remember well)

[2424464575]

Originally Posted by Jazz

Hey everyone,

I am modelling a solidification problem where we are tracking how inclusions move in a caster. They are represented by inert particles and I needed a way by which I can make a particle reach a wall/front as close as possible so that the distance between the particle and the wall is in nanometer scale. I am using*steady particle tracking and my distance comes out to be of order of 10^(-5)m. Even after changing the step length factor the distance comes out to be of the same order.

Recently, I was doing the inclusion behavior simulation of molten steel after solidification. The solidification part has been calculated. Now I am ready to add the calculation of inclusion particles. I have encountered problems in UDF compilation. Do you have a teacher or senior brother to discuss it together? Thank you. My email: 2424464575@qq.com , thanks again! Good wishes!
My code is as follows:

1. When the volume fraction of liquid phase encountered by the particles is less than 0.6, the particles are considered to be captured by the solidified shell of liquid steel.


#include"udf.h"

DEFINE_DPM_LAW(stop_dpm_law,p,ci)
{
real fl;
cell_t c=P_CELL(p);/* Get the cell and thread that the particle is currently in */
Thread *t=P_CELL_THREAD(p);/* pointer to the thread of the cell that the particle is currently in*/
fl=C_LIQF(c,t);/* liquid fraction in a cell*/
if(fl<=0.6)
P_VEL(p)[0]=0.0,/* particle velocity component. */
P_VEL(p)[1]=0.0,
P_VEL(p)[2]=0.0;
else
return P_VEL(p);
}

What are the problems in my code? Please help me. Thank you!