DPM模型中离散相如何加入颗粒的受力？

ytuym

在DPM模型中，在颗粒的运动方程中，除了考虑曳力以外，还需要考虑有浮力影响的重力和附加质量力，如何在Fluent软件中实现？如需要加入UDF，改如何编写程序？

zoeywang

QQ：2284170287

ytuym

谢谢你的答复，可以在这里简单说说吗？

sxlflgj

不知道这个解决没？我刚接触这个，也想弄清楚这个，希望指教。。谢谢。

5超7

好像不能

fty0083

可以加浮升力。DPM有一个体积力的宏，你计算好浮升力，然后根据体积力的那个宏折算成体积力就行了。其他的附加力也可以如此计算，计算一个综合体积力反馈给方程的体积力项。具体的，你看下UDF帮助吧。

fty0083

可以加浮升力。DPM有一个体积力的宏，你计算好浮升力，然后根据体积力的那个宏折算成体积力就行了。其他的附加力也可以如此计算，计算一个综合体积力反馈给方程的体积力项。具体的，你看下UDF帮助吧。

5超7

兄台DPM搞的很深刻啊！以后DPM方面的问题，还希望兄台多出手相助！

fty0083

都是逼出来的，很多东西需要自己实现的时候，就得想办法。
    fluent的dpm模型还是有很多问题的，比如粒径不能小于0.1um（那个量级不能用牛顿第二定律）；再比如它求解的颗粒是质点的运动，若是模拟颗粒和壁面的碰撞效果，只有当颗粒质心运动到壁面时才激活DPM的boundary condition；另外，最近感觉fluent的dpm对控制体依赖性也很大，计算时的载流体速度取得应该是控制体中心速度，这样若是采用较大的壁面网格的话，边界层里面的运动求解的是不对的。并且，实际上fluent的dpm迭代方式是有颗粒自身的迭代步长的，好像设置不了，至今没搞明白和哪个参数有关系。那个所谓的Particle time step只是个摆设。
    给你个建议，好好看udf，学会这个工具自己去反馈各种数值，然后搞懂fluent的dpm里面到底怎么计算的。不然又错误自己都不知道。做udf深入的时候，到src目录下面看看那些头文件，你会受益匪浅，至少你知道fluent计算的时候考虑那些参数。不过，看头文件的感觉只可意会，不可言传。呵呵。

5超7

非稳态DPM我也算过，但是对颗粒时间步长和连续相时间步长之间的关系我搞得不清楚，这样的情况下颗粒释放的次数我始终没有弄明白。现在都是用稳态计算的，稳态那个好说！

fty0083

那两个就是迭代上的问题，设置的那个particle time step是让你看到多长时间内dpm计算的统计结果。如果你做了双边耦合，那么非稳态计算先计算连续相时间不长上的解，再计算颗粒时间步长。不过颗粒的瞬态迭代上的时间步长呢实际上很短，这个和你选用的dpm迭代scheme有关。

dpm迭代都是计算颗粒运动的轨迹，也就是说你用稳态，它解dpm也是一个个小时间步长算的，不过具体轨道的长度和那个颗粒max number有关。如果你要是看有进出口的一个环境的的话，若max number很大，估计稳态结果应该是颗粒都飞走了。因为是稳态么，无穷远处。