CFD中文网

https://www.dataphysics-instruments.com/us/knowledge/understanding-interfaces/surface-tension/
这里有一句

Thus, for liquids the surface tension and the surface energy are equal

至于动能，k就是湍动能啊，如果不只是要湍动能，那k怎么来的肯定要涉及动能啊。

对于给定的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}

@bestucan 谢谢您的分享，这些我之前看过一点。但是这个的定义就还是在初场中定义的湍流积分尺度，用到的是跟湍流强度等有关的计算。但是这样的公司似乎不能用在模拟结束之后的该参数的计算，具体为啥我感觉我也说不上来:135:

@Dingcy 谢谢

谢谢教授:xinxin:

BTW, you can SAVE 20% when you order online and enter Promo Code ENG19

You can order it here: https://www.crcpress.com/Multiphase-Particulate-Systems-in-Turbulent-Flows-Fluid-Liquid-and-Solid-Liquid/Podgorska/p/book/9781482235685

:baobao:

批注 2019-10-20 163534.png

这已经是比较成熟的方法了 应该aixprocess在15年前就已经成型了 做木材洗涤工艺的除尘优化 如果说无脑的办法 还是直接写碰撞列表做加减法 也可以直接设置某一方质量半径都为0 然后另外一方做简单的质量累加 但是实际上真实的情况比这复杂得多 液固的碰撞的结局以及最后的补集情况是多样性的 这有很多相关文献 但是对于UDF编写来说 因为是不同injection 还是有一定难度的

围观中

可以设计一个量纲的类，里面包含一个整形数组来记录量纲的分量，如[L, M, T]还要定义量纲的运算。然后设计一个标量的类，里面除了浮点数据之外，再加一个量纲，定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0]，速度的量纲是[1, 0, -1]。

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

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

@程迪 666啊老铁，听君一席话，胜读百篇论文:146: :146: :146: :146: :146:

@东岳 我的边界是没有网格的，但是默认的显示就是在紧挨边界的网格上，就是你说的内网格呗，但是同样的算例，Tutorials里就能把那个数据给放到壁面上，应该是经过处理了吧。

好的，谢谢

哈哈哈，没事，谢谢东岳老师！

当你的学生肯定很幸福

缺这个包：binutils

或者路径没配置对，可以用echo $PATH看看有没有/usr/bin这一项

@史浩 我不是大佬。最近在做一个VAR的项目，考虑温度场-电磁场-流场-组织耦合计算，需要自编程序，有程序验证的问题。对OpenFoam，我没做过这方面的案例，但我想考察2个方面：一是看一下OpenFoam与自编程序的结果对比，二是看一下OpenFoam与自编程序的计算速度对比。
ESR和VAR比较类似，有相通的地方，你如果做了OpenFoam算例的话可以交流一下。

入门不久，尝试着答一下。
我认为openFoam中的边界条件大致分两类：

强制性边界条件，直接规定数值，如fixedValue，fixedMean，这类边界条件的信息会一直存在于求解所得的结果文件，算很多步以后查看文件就会发现边界条件并没有变化； 参考性边界条件，这类边界条件大多是第二类边界条件，如出口定义inletOutlet， 对于流出的气体而言，该边界其实是zeroGradient，这时定义的value只是初始迭代的时候给的一个初值，该边界的数值由流场内部的点插值得到

@micro 网格数量越多，网格尺度越小，计算肯定更加不稳定，您可以查一下CFL数

:chouchou: