CFD中文网

@东岳 谢谢李老师的回复！
我看了下跟公式（7）到（8）的过程很相似，但是为什么公式7到公式8的过程中将公式7左边第四项给略去呢？

有这么一个公式，应该是出自这篇文章里（Spalding, D.B., 1974. A method for computing steady and unsteady flows possessing discontinuities of density. CHAM report 910/2.），但是这篇文章找不到，所以不是很理解这个公式的含义，这个公式也是摘自（Physical and Numerical Modeling of Unsteady Cavitation Dynamics）这篇文章里的，应该是由这个公式可以推导出Schnerr-Sauer空化模型的。

但是从这个公式里感觉是把密度对xyz方向的偏导数跟省略掉了，所以把密度rho提到了对流项的括号外面了，请问大家，这个公式的具体含义是什么呢？运用的条件是什么呢？或者有谁能找到（Spalding, D.B., 1974. A method for computing steady and unsteady flows possessing discontinuities of density. CHAM report 910/2.）这篇文献的，不胜感激！

@bart 慢慢增大进出口压差计算，在看看结果

@东岳 准备尝试下，因为非定常计算感觉挺耗时间的，所以我基本用的定常计算的，不过之前有试过在最小压力到2bar的时候，切换非定常计算，还是容易出现问题，模型尺度大概300到1000微米，我用的10e-9 s,还是出现发散，是不是得再低些呢，这会儿试下从一开始就用非定常计算，慢慢减小进出口压差，算算看下

李老师，那一般松弛因子的范围又大概是多少呢？我的case里密度是随压力变化的分段函数，当绝对压力压力降到2000Pa后，整个密度会随着压力降很多，从1000kg/m3到0.01kg/m3，我觉得case发散的原因是，每次迭代压力变化的太大，整个流场更新的太快，可能忽略中间的某些流场情况，所以我是想将压力的梯度变化变得缓慢些，慢慢捕捉流场变化，那如果松弛因子不能不能变得很小的话，还有其他方法可以减缓梯度变化吗？

请问，通常认为松弛因子的范围在0到1之间，但是一般大家减小松弛因子，最小大概在0.1左右吧，那么对于收敛性比较难的case，是否可以将松弛因子设置到0.01，甚至0.001呢？在我的case里压力梯度比较大，全局的绝对压力大概在10Pa到60MPa之间（单流体形式的空化模型），在绝对压力到饱和蒸汽压之下后，整个收敛就很差，经常出现以下错误：（附：为了避免网格的原因，整个模型为简单的2D喷嘴，网格均为正交性网格，纵横比也很小）
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: k
Divergence detected in AMG solver: epsilon
因此为了保证计算稳定，我在松弛因子部分的设置如下：

甚至有的时候得到UDF可能到0.001，请问对于我的这种设置是否有问题呢？另外对于以上出现的发散错误，除了网格方面，且边界条件也没有问题的情况下，是否有其他的解决方法？（或者得在Multigrid中进行AMG设置吗？）

#include "udf.h"
DEFINE_ADJUST(my_adjust,d)
{
Thread *t;
real p,dp;
cell_t c;
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
C_P(c,t) = MAX(0,C_P(c,t));
}
end_c_loop(c,t)
}
}

对于这个程序，是否可以对整个流场的压力最小值进行调节？目前加载之后，计算出的绝对压力依旧有负值（操作压力为0），请问是怎么回事呢？

如果气体的物性参数要求是等熵气体状态方程，也就是如下公式：

请问在fluent中如何实现呢？对于常数C和比热容比该在哪里进行设置？（在不使用UDRGM的情况下）

谢谢您的意见，提取的肯定就是面上各个点的速度数据了，那如果这个时候求液相平均速度，是否就是根据液相体积分数做速度的加权平均呢？

在出口截面上有液相和气相同时通过（mixture气液两相），如何单独的求出液相的平均速度呢？目前尝试了三种，一是面积加权平均，二是mass flow weighted average，三是tecplot里以液相体积分数为加权函数的速度平均（velocity+phase-liquid）,目前三种得出的结果无论数值还是趋势都不太一样且目前也没有实验数据，请各位给予指导！

tecplot里面的，看上去感觉上面是等值线，下面是流线

对于我目前所计算的模型所涉及是湍流模型以及多相流空化模型，在湍流模型尝试使用RSM（quadratic pressure-strain ）时，总是提示错误，而对于linear pressure strain或者Stress-Omega Model 都能够正常计算，并且每次计算前都使用了RNG k-e计算的结果求得初始值，网格部分没有问题，ICEM结构化网格质量0.5以上，skew在0.65一下，请问各位是怎么原因导致的呢？

在cfd-post官方教程里指出可以把cfd-post导出的cvf格式的文件内嵌到ppt中，请问这个操作该如何实现呢？

好像是编译的方法才能正常使用

补充一下，目前用profile成功将动网格算出来了，说明我的网格以及设置是没有问题的，可能就是UDF编写或者是编译方面出现问题了，请遇到这情况的朋友指导下

在fluent中使用动网格的时候，在TUI下面显示incorrect cg motion UDF piston on zone 26 (assuming no motion)，并且预览网格毫无变化，请问这个如何解决呢？

#include "udf.h"
#include "dynamesh_tools.h"
DEFINE_CG_MOTION(piston, dt, vel, omega, time, dtime)
{
/* reset velocities */
NV_S (vel, =, 0.0);
NV_S (omega, =, 0.0);
if (time < 0.00001) 
   {
     vel[1] = 15;
   }
  else if (time < 0.000015)
   {
     vel[1] = 0;
   }  
  else if (time < 0.000025)
   {
      vel[1] = 15;
   } 
  else
   {
     vel[1] = 0;
   }
}

我的理解是，absolute格式表示的是你当前计算所得的残差与你设置的残差进行比较，relative格式表示的是你在当前时间步所在的迭代数的残差与你这个时间步一开始的残差间的比较，比如你设置的是10-5，而你这个时间步的开始残差是10-2，那么你最后收敛的要求应该是10-7。我的理解是这样，可以参考下！

scalable wall functions应该跟 Standard Wall Functions差不多，只是多了个，无论边界层网格怎么细化，它都能得出一致的结果

另外，如果EWT与MLT是求解近壁面流动的，那么对于边界层本身很薄的情况，是否是一半的壁面函数就能满足要求了

（1）help文件中所提及的 y+-independent与y+-insensitive究竟是什么意思呢？
（2）Enhanced Wall Treatment 与Menter-Lechner near-wall Treatment得到的结果就肯定比标准壁面函数或者非平衡要好吗？（满足网格条件下，网格的问题第三问）
（3）help文件提及EWT与MLT时属于two-layer model，其实在网格粗或者细的时候都是可以使用的，那么网格的粗细对其结果有影响吗？另外如果需要求解低雷诺数近壁面流动，那么对EWT于MLT就必须使得y+小于等于1吗？