CFD中文网

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

生活很有规律呀！我就不行了，生物钟紊乱，感觉每天都很疲惫，特“困”生就是我！:haqi:

@veen 谢谢，我学习一下

@一朵苍凉 我也不太懂，可以看看控制领域的相频特性，也许有帮助:zoule:

块就是将你的整个计算区域分割成几个区域，每个区域有各自的网格分布规律；
模型应该是你计算的几何外形吧

kOmegaSST和kEpsilon印象中是比较好收敛的两个模型，其他的用起来似乎都不太顺手，当然最重要的是结果正确。

@veen 从源码编译of是个挺痛苦的事情，建议不把of安装在root，可以考虑试一试easybuild，它是基于module的环境管理，应该比较成熟了。我从0开始，最开始在PC上面编of-3.0，到后来在实验室的机群上home编译of-4.0, foam-extend等等，多版本的管理用module还是方便些。印象中，gcc（如果你用gcc的话）和mpi的版本很重要。

@西湖冷月
nice兄弟，d=====(￣▽￣*)b

@浪迹天大 ControlDict里面的writeControl和writeInterval似乎会覆盖前面的？

application twoPhaseEulerFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 2; deltaT 1e-4; writeControl timeStep; writeInterval 100; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable on; adjustTimeStep no; maxCo 0.9; maxDeltaT 1e-04; functions { volFieldValue1 { type volFieldValue; libs ("libfieldFunctionObjects.so"); log true; //writeControl writeTime; writeControl timeStep; writeInterval 1; //"这里的writeInterval 1 似乎会覆盖前面的100导致所有的输出都是每个时间步都输出？" writeFields true; regionType all; operation volAverage; //weightField alpha1; fields ( alpha.particles U.air U.particles p p_rgh ); } }

大家好，我在字典文件中看到

xcells #calc "round((($xmax)-($xmin))/$delta)";

想知道#calc 可以使用的标准函数包括哪些呢？哪里能查得到？

$\tau$这个雷诺应力，对角线的是normal stress，非对角线的是shear stress。不过你说的更详细，感谢分享！

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

@ibelief 在 求教fluent后处理-导入曲面（线）-输出曲面（线）参数的实现方法 中说：

或者是分割或者什么，探索一下。 祝

好的 谢谢老师

@evensun 嗯，涡粘系数是二阶张量这个按我理解其实就是非线性湍流模型了

厉害厉害了老铁 你回国了么？

对于给定的NDF，划分为$i$个$N_{pp}$，对每个$i$上定义$k$阶矩$m_k^i$，给定$m_k^i$，可以计算第$i$区间的节点$d^i_0,d^i_1$以及权重$w^i_0,w^i_1$：
\begin{equation}
\begin{split}
w^i_0&=w^i_1=0.5
\\
d^i_0&=m_1^i-\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\\
d^i_1&=m_1^i+\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\end{split}
\end{equation}
对于仅考虑破碎的PBE：
\begin{equation}\label{pbe}
\frac{\p n(d)}{\p t}=\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'-g(d)n(d)
\end{equation}
对方程\eqref{pbe}在$i$上取$k$阶矩：
\begin{equation}\label{m}
\frac{\p m_k^i}{\p t}=\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')d^k\beta(d|d')n(d')\rd d'\rd d-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}
\begin{equation}
\begin{split}
\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'\rd d&=
\int_{d_{i-1/2}}^{d_{max}}g(d')n(d')\left(\int_{d_{i-1/2}}^{d'}\beta(d|d')\rd d\right)\rd d'
\\&=
\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)
\end{split}
\end{equation}
Therefore
\begin{equation}
\frac{\p m_k^i}{\p t}=\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}

@东岳 我是按这个公式计算的，密度为1，那么u_tau=sqrt(tau_w)，但是我发现，论文里的摩擦速度刚好是我用OF算出来的1.225倍。论文是用fluent算的，fluent里面的空气密度是1.225。我现在想不通 为什么会差这1.225倍。:jingya:

https://cfd-china.com/topic/3742/结构化网格/8?_=1615050206338
倒数第二个帖子

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？ 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
3.png

@vllst 你好请问你的问题解决了吗？