CFD中文网

同理，下述网页中都需要处理$\sum_f p_f^{t}\bfS_f$为$\frac{1}{V_\rP}\sum_f p_f^{t}\bfS_f$，均已更正，谢谢

http://dyfluid.com/rhoPimpleFoam.html
http://dyfluid.com/rhoSimpleFoam.html
http://dyfluid.com/simplefoam.html

@李东岳 在 HyperMesh为STARCCM+生成面网格 中说：

感谢分享！有益于后来人

共同进步！

@tyloryu 明白了，感谢！！！

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。
但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。

所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的:135:
在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。
我把相关资料发给老师。

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

@李东岳 好的。以前在CFD-DEM的时候要几千个颗粒，DNS的话应该得缩短一下管道长度，算少一点颗粒量。

有人用过Microsoft Azure吗，可以分享一些入门学习经验吗？

Sep 25, 2018 , Jan 08, 2019, Required Reviews Completed

@LYHua 太感谢了，解决了，给力

@李东岳 感谢李老师的回复，一下子就通透了

:baobao: :baobao:

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。

密度就是单位体积的质量，速度的物质导数就是加速度。
所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？

力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。

物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。

要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。

另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。

我碰巧最近在看这方面的工作，打算在OKS课进行植入，中文还对不上，但是应该就是楼上说的 转捩模型。我参考的是Menter的文章A One-Equation Local Correlation-Based Transition Model，他们的参考案例也很有意思，T3A系列，单独用层流、湍流都不行，只能用附加转变的模型。下图：

捕获.JPG

@东岳 多谢东岳老师的回复。VOF+自适应在OpenFOAM里面有植入吗？
目前我是做一个通电导电流体中液滴的运动，在电流不大的情况下，液滴很少，就一两个；但是在大电流的情况下，液滴在电磁力作用下会破碎，这个估计会有小液滴产生。
整个属于工程问题，不用特别的精确，但是算的要快一些，需要模拟液滴600s的运动情况。目前暂时在interFoam的基础上开发，但是感觉还是不够快，这才有了上面自适应网格的想法，其余的还有没有其他计算效率比较高的求解器？

用的windAroundBuildings算例，840万网格的二维模型，网格质量检查过了，计算域是14000*4000米。epsilonInlet 0.0003; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $epsilonInlet;
}

OUTLET { type zeroGradient; } WALL { type epsilonWallFunction; value uniform $epsilonInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
INLET
{
type fixedValue;
value uniform $kInlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type kqRWallFunction; value uniform $kInlet; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
nut
internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}

OUTLET { type calculated; value uniform 0; } WALL { type nutUSpaldingWallFunction; value uniform 0; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
p
internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
p0 uniform 0;
}
WALL
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}

#includeEtc "caseDicts/setConstraintTypes"

}
Uinlet (10 0 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform $Uinlet;
}
OUTLET
{
type zeroGradient;
}

WALL { type noSlip; } frontAndBackPlanes { type empty; } #includeEtc "caseDicts/setConstraintTypes"

}
时间步长试过1 0.1 0.01 其他的都没动

可以相信这一家：https://www.phoronix.com/review/amd-ryzen9-7950x3d-linux/

@maoyanjun_dut

对于 Vim 用户 Spacemacs 了解一下？

这个是自己写的求解器，所以有可能差别确实很大

感谢东岳老师，我尝试着改小了，问题得到了部分解决。

考虑最简单的动量交换，液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA
\end{equation}
乍一看，如果$\bfA$是负的，那么会导致$\bfU_\rd$向下走（为负），$\bfU_\rc$（为正）。但是这并不符合物理，考虑一个管子的颗粒，如果颗粒向下走，必然会带动空气同时向下走。但是方程缺不是这种体现，为什么呢？

实际上，上述方程并没有写完整，完整形式应该是这样：
液相：
\begin{equation}
\frac{\rd \bfU_\rc}{\rd t}=-\bfA=-\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
气相：
\begin{equation}
\frac{\rd \bfU_\rd}{\rd t}=\bfA=\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
\end{equation}
如果开始的时候$\bfU_\rc=0$，则为