CFD中文网

of中多相流这块不可压好像都没有能量方程，在模仿着可压的能量方程添加时发现，剪切力做功是用湍流粘性为形参的库函数表示的，现在我不调用库函数，怎么表示这一项呢？谢谢

:chouchou:

各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？

这个其实是个好话题，很多大佬在sci里面也明确说了，LES之类的模拟，如果对于这种类似的圆管，施加cyclic边界条件，会对湍流产生抑制。但是我也没见过大佬们进行对比。

另一方面，感觉边界条件，直觉感觉，单单的一个patch引起的差异，不至于那么大。LES模型参数的影响更大。如果谁有类似的文章，可以发来看看。

顺便测试一下侧边栏回复

我用子系统上的OpenFOAM，总感觉并行计算没用，和我不开并行计算时间差不多。。。而且还decomposePar不了，好像是安装的时候少了啥东西。。。

@李东岳 感谢东岳老师，非常详细的步骤。我会好好看看这块

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

@dzw05
哦了，有空整理整理。

在使用Fluent 6dof+overset+隐式更新时，发现写出的motion history中的重心坐标与实际明显不符，差的很大，而在udf中通过DT_CG得到的却是正确的重心坐标，不知道有没有人遇到过这样的问题，是什么原因导致的呢？

@李东岳 您好！感谢您的提议。首先rhoCentralFoam应该是可以算湍流的，我是基于这个前提，想知道除了在constant文件下turbulenceModel选择RAS或LES以后，如何改动其system文件夹下的fvsolution和fvscheme字典以及在0文件夹下新加相应的变量字典来计算湍流。不知道我的这种想法是不是正确的。

@禾斗石开 fluent

您好，我也在研究topset和refinemesh的问题，您知道在blockmesh生成一个梯形（具有斜面）的网格时，使用refinemesh会出现报错是什么原因吗？

黑圈那些edge，都不能关联，要取消关联

各位大佬好~向大家请教一个问题。

式（1）是大涡模拟法夫尔滤波操作以后的连续性方程
\begin{equation}
\frac{\partial{\bar{\rho}}}{t}+\frac{\partial{\bar{\rho}\tilde{u}_i}}{x_i}=0
\end{equation}

如果对$\tilde{u}_i$取时均值，即$<\tilde{u}_i>$, 这样得到的是雷诺时均值的速度吗？如果想要取法夫尔时均值，是不是需要做$<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>$的处理？P.S(这里用<φ>表示对φ取时间平均)

第二个问题与第一个问题相似：基于法夫尔时均的可解尺度的湍动能的求解是否可以按照式（2）不可压缩流体的求解方式求解呢？我感觉式（2）的解更像是雷诺时均的解，因为可解尺度的湍流里面也包含了密度脉动。
\begin{equation}
k_{resolved}=\frac{1}{2}(<\tilde{u}_i\tilde{u}_i>-<\tilde{u}_i>^2)
\end{equation}

还是说需要用式（3）继续用法夫尔时均进行计算呢？
\begin{equation}
k_{resolved}=\frac{1}{2}\big(<\bar{\rho}\tilde{u}_i\tilde{u}_i>-<\bar{\rho}>(<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>)^2\big)/<\bar{\rho}>
\end{equation}

这个放在他们搞涡那面，估计能讨论一上午 :duang:

我个人也不太熟悉。在我的文章里面，也不太关心涡的情况。

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

@东岳

很难申请，不过我会努力的！

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢