CFD中文网

@leo_chou 看过一个挺不错的视频：https://www.bilibili.com/video/BV1Yt411W7td

学OpenFOAM至今一直水在某大佬的OpenFOAM千人群中，个人认为还是比较活跃的之一(毕竟也水到50多级了)。
但是今天早上一起床，发现我被群内某管理员踢掉了：
qq_pic_merged_1629032541800.jpg

我寻思我也没发什么违禁言论，也没做什么违法群规的事情啊，且前天我还在群里解答疑惑:134:
关于这个问题我私聊了一下管理员，结果发现我问完他直接给我删了？:136:
30bfc521-cae5-4276-8a6b-89c6da167351-Screenshot_20210815_210758_com.tencent.mobileqq.jpg
145f5c52-5a8e-47e9-9341-302bf79191c7-Screenshot_20210815_210832_com.tencent.mobileqq_e.jpg

我不是很懂这波操作，按理说踢人不应该有个规范吗，还是按管理员想踢谁就踢谁吗？

PS：不知道这里是否允许对这类问题进行讨论，侵权麻烦管理员大大直接删帖吧

为什么我用的

forAll(prticles,i) { particles[i] }

这个particles[i]取不到particle

如果只是查看，emacs 有个 whitespace mode。
空格、tab、换行之类的特殊字符都可以显示
vim按说应该也有，没搜到

Screenshot from 2021-05-22 15-16-48.png

国内外有哪些流体力学经典著作？ - TurbPhys的回答 - 知乎
https://www.zhihu.com/question/335644582/answer/822877093

替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程，现在模型验证一直没弄好。首先是编写了一个简单的UDF，采用热平衡模型，导入fluent后提醒我定义比热的那一块有错误，求各位前辈帮我看一下；其次，我添加了Boussinesq项后，优化了一点点液相云图中熔融线一直平行的问题，但是与论文中的还是有差距，且我的熔融线上下厚度保持一致（实际上应该下面厚很多），UDF中还需要添加什么项来解决这个问题呢？以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图，求各位帮我解惑，不甚感激。2.png 3.pn替代文字g 4.png

如果不考虑时间项
\begin{equation}
\int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
\end{equation}
把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来，看起来是守恒的，剩下的只有边界条件

@ljpgbl
这个是全貌
0_1539505430619_Clip_20181014_162251.png

@东岳 老师，我重新画网格后，decomposePar经常会报错，但是不影响后面计算，没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

@东岳 李博是不是在植入这个算法，我等你:ok:

感谢，有一个问题，为什么不把入口直接设置为湍流入口，是因为计算的雷诺数不高嘛

a80fed28-45ef-4507-a08d-13abf19c7ccc-image.png 这个udf具体该怎么编写和使用？

通过vorticity求解涡量； 写个程序求解$\nabla w$
2.1 在这个程序里面，$\nabla w$的边界场$\cdot\bfn$再乘以粘性系数

说简单也简单，只不过繁琐了一些，好像没有几个命令就能出来的东西

@xjwang

CFDRFG的附录code还有个问题，求解$K$，不是按照论文中的原理计算的。他先用$PQ$没有代入$S(f)$的值，算了一组$K$，然后所有坐标点都用了同一组$K$。而根据理论，应该是要根据每个坐标点的非均匀$PQ$，求对应该点的$K$。DSRFG如果生成非均匀，理论上也要这么计算。
739858e5-ac4f-4736-94bd-9f4ec9f07a9d-image.png

关于CDRFG方法，他们的code里是对每一个点都生成一组PQK

最近看到一篇JFM文章《Guo, H., et al., An efficient and low-divergence method for generating inhomogeneous and anisotropic turbulence with arbitrary spectra. Journal of Fluid Mechanics, 2023. 970: p. A2.》，用的就是三维能谱，不是频谱，生成非均匀湍流场也是low-divergence。

更新科氏力源项：$\mathbf{f}\times\bfU=2\Omega\sin \lambda \times\bfU$，其中$\mathbf{f}$表示科氏力参数，一些文献会直接给出$\mathbf{f}$的值。如果绕z旋转，那么$\mathbf{f}$在$z$方向存在值。在这种情况下，$\mathbf{f}\times\bfU$的值为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
f_2u_3-f_3u_2\\
f_3u_1-f_1u_3\\
f_1u_2-f_2u_1\
\end{matrix}
\right]=\left[
\begin{matrix}
-f_3u_2\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
如果进一步考虑速度方向是x方向的，那么变为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
0\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
可以看出如果地球绕z旋转，对于x方向的速度，科氏力导致y方向的偏离。

$\Omega$表示地球旋转速率，单位为[弧度$/\mathrm{s}$]=$1/\mathrm{s}$，$\lambda$表示维度[无单位]。因此$f$的单位为[$1/\mathrm{s}$]

可以在join/intersect里面查找一下看看。

@两月三年 good ，好主意。我试试哈。

adept

2422ba2b-f9b8-402e-936e-8075325aae10-image.png 例如我生成一个part名字叫多“cavity1”，但是导入到fluent里面会多出“cavity1：007，cavity1：008，cavity1：008shadow，”这几个边界fb6ad2d6-0c7d-4266-a7b4-fa752482f0a6-image.png 。相应的网格也被分裂出好多细块。微信图片_20190126110946.jpg 微信图片_20190126110950.jpg

我想建立一个水槽，水流从水槽一边冲过去，然后又从另一边冲过来，两边的入口和出口边界要随时间变化，有没有大神指导一下，怎么实现？谢谢！