CFD中文网

今天空闲时间又看了一遍笔记发现有些笔误：
P9-第三章第一节第二段第五行“将手竖直的防止在水中” 应为“将手竖直的放置在水中”

可以把正确的UDF发来参考一下吗？我也做类似的复现研究，谢谢！

例如冬天拎起一桶水旋转，在旋转过程中水慢慢结冰，想要采用fluent中Solidification&Melting模型和Moving Reference Frames模型来模拟这一过程。
跟仅开启Solidification&Melting模型的模拟不同，但似乎收敛性很有问题。凝固/融化模型与MRF是否可以联合使用？有谁做过类似的模拟吗？

这个啊，先排错，
看看是“只有20可以”，还是“数字小了才可以”。
比如试试19，21，试试200。

错误定的越准，描述才越准，也越容易搜到前人的经验

想以后从事结构方面的工作结合着CFD，基本就是风工程这块了。想找找这方面的职位，计算，研究，工作都可，想问下有这方面的同学吗，可以请教学习一下，谢谢！
不知道这里发帖是否合适，先谢谢大家！

dos2unix命令是用来将DOS格式的文本文件转换成UNIX格式下的文本文件，这两个系统下的文本文件的差异在于换行标志，DOS下是用\r\n断行，UNIX下采用\n作为断行。将windows下面写好的一个文本文件转移到linux下，无法直接运行，打开会发现出现很多^M的标志，可以采用dos2unix命令进行清除。dos2unix命令起且仅起这个作用。因此，对于windows下生成的mesh文件，如果转移到linux下能够正常运行跑通算例，则完全不需要dos2unix命令。
同样，checkMesh也是，如果算例都能跑，则不需要这一步骤hhhh:papa:

想了解一下，fluent中做的湍流生成程序类似CDRFG的方法，在OpenFOAM中文件改动会很大吗。本人做的是风荷载模拟，基于课题组基础上做的课题，课题组有师姐在fluent上做的湍流风入口，对课题了解还处于小白阶段。

在安装完Windows版ParaView之后, 也可以去WSL的.bashrc文件里面添加以下内容

alias paraview='"/mnt/ParaView/安装/路径/bin/paraview.exe"'

然后重新Source一下就可以使用paraview命令启动Windows版ParaView了, 可以直接打开WSL里面的文件夹, 无需在Windows目录下运行OpenFOAM算例

我觉得STARCCM+的网格功能挺6的，笛卡尔/多面体等，控制的参数有限，基本上零基础跟着教程做几遍就差不多了，网格生成速度也很快，比较复杂的几何还是很容易解决的；
其他的用point wise多一点，复杂的几何自己划分块有点麻烦，如果不能分层很容易混乱，简单几何还是非常便利的；ICEM复杂几何切起来还是比较舒心的；^_^

感谢分享。不过虚拟机的话用UTM（免费）很方便，直接加载官方镜像安装即可。我在UTM里的Ubuntu20.04(Focal)里成功编译了OF8，9和OF v2012。

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

@知行 算例文件[链接文本](链接: https://pan.baidu.com/s/1gAkf3WA5S1GNmoDPnvMF0g?pwd=spm6)百度网盘链接: https://pan.baidu.com/s/1gAkf3WA5S1GNmoDPnvMF0g?pwd=spm6 提取码: spm6

@bleach2012 这两个链接现在都，失效了:chouchou:

@simthere 可以看看流量计压降信号转流量信号公式怎么算的，如果不修正悬臂压降的影响，用二维算是最经济的。如果里面有对悬臂压降的修正，那用二维是不合适的。说到底就是这个悬臂的压降能不能被忽略，如果不能，必须用三维

@liujunCFD Snipaste_2021-03-04_11-27-23.png

@李东岳 谢谢李老师回复

在气液欧拉多相流的大涡模拟中，有没有什么模型是把liquid或者particle相的SGS viscosity和主相的SGS viscosity联系起来的？
就像RANS中的一种模型：
$\begin{equation}
\frac{\nu_{t,particle}}{\nu_{t,gas}}=(\frac{k_{particle}}{k_{gas}})^2=\frac{1}{1+\frac{t_p}{t_{fl}}}
\end{equation}$
其中$\nu_t$是湍流粘度，k是湍动能，$t_p$是particle relaxation time，$t_{fl}$是Lagarangian fluid time scale。
现在像找一个类似的模型用在LES上，但是翻了好久也没翻到 :135:

@chon01 这也只是我的猜测，因为我按了一下计算器，lg和ln函数算出来的结果差距确实比较明显，希望能帮到你吧。

@Leon-Zhao 感谢大兄弟

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