CFD中文网

想问下楼主最后成功实现仿真了，我目前已经完成了UDF的编写，并且在fluent中也能正常计算，但是液膜厚度一直保持在第一层网格内，不知道是什么原因。希望可以交流一下。

谢谢李老师，确实是alphatWallFunction出了问题。重新修改了这个函数后，就可以了。再次感谢。:xinxin:

@thegame 感谢老铁，将在下一版本更新

@lys 用的已有模型，专门针对非牛顿流体的曳力模型还没用过，希望能进一步交流啊

提醒读者： 因近期图床地址变动，以上贴子中的图片可能失效，这里附上我个人博客中的相应页面（可以正常浏览），以供后来者参考：写给工程技术人员的 LaTeX 安装指南

@wwxdli 谢谢，我试一试

@东岳 内网格，尝试着接着划吧

我也会这样子，我现在也不知道是什么原因

是OpenFOAM吗？
如果是的话参考
https://www.cfd-online.com/Forums/openfoam-meshing/61596-cyclic-bcs-pointwiseopenfoam-export.html

@lyc 你好，我很少遇到。即便遇到都是重启就能解决的那种。

@Qing 感谢老师！您说的特别清楚，这对我特别有帮助。再次谢谢您！

工科末流985硕士毕业，想在武汉找CFD相关的工作，家人们有没有推荐的公司。

就是下面这个
$$
\frac{\p T_{ij}}{\p x_i\p x_j}
$$
其实$T_{ij}=u_iu_j$，所以演变为
$$
\frac{\p u_iu_j}{\p x_i\p x_j}
$$
这个公式是可以写成$\nabla\cdot(\nabla\cdot(\bfU\bfU))$这种形式，$\frac{\p T_{ij}}{\p x_i\p x_j}$好像写不出来 :134:

大佬们看看？

@李东岳 谢谢李老师回复，前几天关了一下，没起作用。我用其他的服务器跑了一下Fluent。

@李东岳 实在是有点困难~请问老师如何让两个场（如图片中的a和vDotAlpha）对应网格的值相乘，这不是内积也不是外积，用forAll循环的话，会报assignment of read-only location错误。

forAll(a,celli) { vDotcAlphal[celli] = a[celli]*vDotcAlphal0[celli]; vDotvAlphal[celli] = a[celli]*vDotvAlphal0[celli]; }

图片2.png

很有帮助。

你好，请问现在有basilisk的交流群吗，

在opemfoam中有一个内聚力，可以使粒子之间相互黏结。它是通过粒子之间的碰撞的能量转换实现的。其中通过一个叫cohesionEnergyDensity的系数来决定它的强度。这个密度越大，黏结越明显。希望我说的这个是你想问的。

@Jianwen 已解决

测试一下能否随机显示主题