CFD中文网

二变量高斯分布：
\begin{equation}
f(u,v)=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left(-\frac{1}{2(1-\rho^2)}\left[\frac{(u-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(u-\mu_1)(v-\mu_2)}{\sigma_1\sigma_2}+\frac{(v-\mu_2)^2}{\sigma_2^2}\right]\right)
\end{equation}
MGF为：
\begin{equation}
m_{i,j}=\exp\left(i\mu_1+j\mu_2+0.5(\sigma_1^2i^2+2\rho\sigma_1\sigma_2ij+\sigma_2^2j^2)\right)
\end{equation}
纯矩计算方法为
\begin{equation}
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=\mu\\
m_{2,0}&=\mu^2+\sigma^2\\
m_{3,0}&=\mu^3+3\mu\sigma^2\\
\end{split}
\end{equation}
假设$\mu_1=10,\mu_2=20,\sigma_1=\sigma_2=2,\rho=0.5$，有纯矩：
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=10\\
m_{2,0}&=104\\
m_{3,0}&=1120\\
m_{0,1}&=20\\
m_{0,2}&=404\\
m_{0,3}&=8240\\
\end{split}
同时有混合矩
\begin{equation}
m_{i,j}=\exp\left(36\right)
\end{equation}
不行，混合矩计算方法不对

:136: 好吧找到原因了，是这个版本的link edge factor默认是0

@星星星星晴 大佬你好，我试了下，但是不知道具体提供应该在哪添加2024-10-29 23-45-50 的屏幕截图.png 2024-10-29 23-45-55 的屏幕截图.png 2024-10-29 23-46-04 的屏幕截图.png

请教各位：
目前有一个6自由度运动的sloshing两相流模型，通过udf文件给定运动速度。

需要测出壁面位置某些点的压强时程变化，以及某些区域的水位（波高）时程变化。请问该如何设置这些测点，使其与动网格区域一起运动。

\begin{equation}
\zeta =x-ct,\eta=x+ct
\end{equation}
\begin{equation}
\frac{\p\eta}{\p x}=1,\frac{\p\zeta}{\p x}=1
\end{equation}
\begin{equation}
\frac{\p\eta}{\p t}=c,\frac{\p\zeta}{\p t}=-c
\end{equation}
对任意变量关于$\zeta,\eta$的函数
\begin{equation}
\frac{\p}{\p x}=\frac{\p}{\p\eta}\frac{\p\eta}{\p x}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p x}=\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}
\end{equation}
\begin{equation}
\frac{\p}{\p t}=\frac{\p}{\p\eta}\frac{\p\eta}{\p t}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p t}=c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}
\end{equation}
\begin{equation}
\frac{\p^2}{\p x^2}=\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)=\frac{\p^2}{\p\eta^2}+\frac{\p^2}{\p\zeta^2}+2\frac{\p^2}{\p\eta\p\zeta}
\end{equation}
\begin{equation}
\frac{\p^2}{\p t^2}=\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)=c^2\frac{\p^2}{\p\eta^2}+c^2\frac{\p^2}{\p\zeta^2}-2c^2\frac{\p^2}{\p\eta\p\zeta}
\end{equation}

请问有前辈知道gmesh如何能划分如图所示得网格吗，图1是图2的block

block-tuya.jpg model-tuya.jpg

多次调试发现就错在mesh.C().internalField()这个调用上

@haozhibin @youhaoyu @wangfei9088 感谢三位大佬

搬运大佬🅲🅵🅳.ⓒⓞⓐⓒⓗ的解决办法，实测好用：

“根据你的描述，估计是 system/controlDict 里面有 functionObjects，可能是一个，也有可能两个都有，针对你的 case，应该是 forces 的部分，如果把这部分屏蔽了，应该可以了。”

如题，Paraviw如何求同一个场在不同时间步下相同的值（没有发生任何变化）的那部分网格，即做交集处理？例如用interFoam时，alpha=1代表水，我想获得两个时间步下alpha始终等于1的那部分网格。
提前谢谢各位回答！

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

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

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

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

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

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

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

@winsway_zero 前辈，请问gmesh怎么划分block

针对一个高度为3mm，宽度为1mm的矩形微通道，用fluent mesh 画网格后，在入口流速为0.3 0.4 m/s时，计算都可以收敛。但是当入口流速为0.5m/s时，计算收敛不了，这可能是什么原因导致的，尝试过增加网格数量也无效

请教各位老师，外力驱动的流动，常用哪些算例做benchmark，比如Poiseuille流动，非定常的Womersley流动，还有吗？

@东岳 总压减小是因为速度变小，还是静压变小，或二者都有？无粘不可压的时候就是动能与压力势能的转化，是守恒的。

@东岳 那请问东岳老师，在t=0，和t=1时刻，这个函数返回的值是什么意义呢？

完全 - > 几乎

@maoyanjun_dut

对于 Vim 用户 Spacemacs 了解一下？

以后可以有积分制度了，回帖、在线等都有积分，积分可交易等，也和人民币挂钩，分越多，能换的人民币越多

:high: :high:

主要为了刺激大家积极性，也为了给这些一直帮助大家解决问题的人实质性奖励