海蓝浅唱 发表于 2015-11-30 22:26:26

求教前辈，加载密度函数UDF后发散，不能计算。

做齿轮的润滑分析，UDF中包含定义齿轮啮合传动的程序，流场密度随压力变化的程序，run后迭代几部就出错。当密度为常数，只定义udf啮合传动时收敛很好,当只用udf定义密度函数，齿轮传动用稳态做时也能正常迭代。
# Divergence detected in AMG solver: pressure correction -> Turning off correction scaling!
# Divergence detected in AMG solver: pressure correction -> Increasing relaxation sweeps!

海蓝浅唱 发表于 2015-11-30 22:28:24

自己顶一下，UDF如下：
#include "udf.h"
DEFINE_PROPERTY(my_density,c,t)
{
    double P=C_P(c,t),a=856.12;
    double density;
    density=a*(1+6*P*pow(10,-10)/(1+1.7*P*pow(10,-9)));
    return density;
}

DEFINE_CG_MOTION(roting , dt, vel, omega, time, dtime)
{

omega=-1000;
}
DEFINE_CG_MOTION(roting1 , dt, vel, omega, time, dtime)
{

omega=896.55172;
}

海蓝浅唱 发表于 2015-11-30 22:28:41

自己顶一下，UDF如下：
#include "udf.h"
DEFINE_PROPERTY(my_density,c,t)
{
    double P=C_P(c,t),a=856.12;
    double density;
    density=a*(1+6*P*pow(10,-10)/(1+1.7*P*pow(10,-9)));
    return density;
}

DEFINE_CG_MOTION(roting , dt, vel, omega, time, dtime)
{

omega=-1000;
}
DEFINE_CG_MOTION(roting1 , dt, vel, omega, time, dtime)
{

omega=896.55172;
}

Jolie 发表于 2015-12-8 21:21:09

有类似的问题，同问，还请高手解答

漫舞飞沙 发表于 2016-4-28 20:41:54

可以试试不使用密度UDF，还可以适当降低压力的松弛因子

没晓晓 发表于 2016-8-21 22:31:01

同样的问题，我也遇到了，请问楼主是如何解决的:)

海蓝浅唱 发表于 2016-11-11 11:02:05

没晓晓 发表于 2016-8-21 22:31
同样的问题，我也遇到了，请问楼主是如何解决的

没有解决，你呢？

海蓝浅唱 发表于 2016-11-11 11:04:06

Jolie 发表于 2015-12-8 21:21
有类似的问题，同问，还请高手解答

你解决没?

海蓝浅唱 发表于 2016-11-11 11:05:50

漫舞飞沙 发表于 2016-4-28 20:41
可以试试不使用密度UDF，还可以适当降低压力的松弛因子

降低松弛因子也不行，导师要求用这个密度函数。
页: [1]
查看完整版本: 求教前辈，加载密度函数UDF后发散，不能计算。