Posts made by zousiyu

@freeazma 如果我想在 fluent 里面同时查看/获取/输出固体和流体温度，能实现吗？

@freeazma 您好！又是我，打扰了。您知道如何同时获取多孔介质区域内的固体温度和流体温度吗？
在后处理时，在多孔介质区域仅能看到一个温度，但实际上该区域应该存在两个温度，即固体温度和流体温度。

@freeazma 好的，谢谢您！

@freeazma 您好！我最近在做一个多孔介质局部非热平衡（LTNE）的模型。模型如下。我发现，开启 LTNE 后，多孔介质区域变成了两个，一个固体，一个液体。自由流动与多孔介质区域的交界面（红色线宽）变成了 wall，对流动不可设置，看结果这条边界应该是可渗透的，对能量默认是通量为0。感觉这里的边界条件不太合适。前辈有什么意见吗？

本科生已经放假了。硕士博士可能二月份放假。

@zhangshengchang 谢谢！

@zhangshengchang 您有没有相关的文献呢？就是计算换热系数的文献。

@zhangshengchang 谢谢您！这个解释很详细！

@zhangshengchang 那我能不能直接用这个热通量 q 去计算传热系数呢？

@zhangshengchang 那您指的这个“换热热流”是扩散热通量，还是扩散和对流热通量之和？

你可以看看这个意大利小组的文献，说不定有帮助。他们自己写了个求解器，专门用于求解固定床的非均相催化现象。系统内同时存在自由流动区域和多孔介质区域。
[1] T. Maffei, G. Gentile, S. Rebughini, M. Bracconi, F. Manelli, S. Lipp, A. Cuoci, M. Maestri, A multiregion operator-splitting CFD approach for coupling microkinetic modeling with internal porous transport in heterogeneous catalytic reactors, Chemical Engineering Journal 283 (2016) 1392-1404.

我在 Fluent 的后处理菜单中看到 Total Surface Heat Flux，不知道这个是不是 the magnitude of heat flux？

从 heat flux 的定义来看，

$ q=-k\nabla T $

那 the magnitude of heat flux 是不是，

$ \sqrt{ (-k\frac{\partial T}{\partial x})^2 + (-k\frac{\partial T}{\partial y})^2 + (-k\frac{\partial T}{\partial z})^2 } $

如果在 Fluent 里面要计算壁面的热通量的大小，是不是要计算温度在各个方向的导数

@李东岳 哇！谢谢！
您收到了信息后可以删这条吗？不想暴露个人信息

@bestucan 谢谢您！学习了

经常在东岳老师的朋友圈看到他自己做的 CFD 周边，比如 NS 方程的鼠标垫，湍流家居壁纸。请问东岳老师是找谁制作的，这么精致，可不可以推荐给我们。

@mobaq 那多半没有现成的模型吧。但你植入一个对流方程应该很简单。

1. 满足第一条；
2. 想要 Turbulence modelling for CFD；
3. 邮箱 zousiyu1995@outlook.com 或者 zousiyu1995@163.com
   谢谢东岳老师！期待能中

@mobaq 看这个方程好像是添加了一个标量传递方程，把扩散项去掉，就可以实现染色了。

回帖涨个币

如果存在纳米尺度，那网格也要能分别纳米尺度，应该要不少网格吧
而且，纳米尺度的现象，能不能用 NS 方程也不好说

@cccrrryyy 牛，分析得透彻，谢谢您！

@李东岳 明白了，谢谢东岳老师

20-30万，这个薪资是做应用还是算法呢？如果是应用，好像不低了哎

最近投了篇文章，有个审稿人给我意见，说，

我的理解是，实验是一种获取数据的方法，CFD应该也算一种获取数据的方法，那我们能说 CFD 是一种 method 吗？比如，The CFD method has been widely used.

你看看这几篇文献，如果不会，我可以有偿

1. Lee, Y.K., et al., Surface mesh generation for dirty geometries by the Cartesian shrink-wrapping technique. Engineering with Computers, 2010. 26(4): p. 377-390.
2. Bu, S.S., et al., On contact point modifications for forced convective heat transfer analysis in a structured packed bed of spheres. Nuclear Engineering and Design, 2014. 270(Supplement C): p. 21-33.
3. Dyck, N.J. and A.G. Straatman, A new approach to digital generation of spherical void phase porous media microstructures. International Journal of Heat and Mass Transfer, 2015. 81: p. 470-477.
4. Rebughini, S., A. Cuoci, and M. Maestri, Handling contact points in reactive CFD simulations of heterogeneous catalytic fixed bed reactors. Chemical Engineering Science, 2016. 141: p. 240-249.
5. Wehinger, G.D., C. Fütterer, and M. Kraume, Contact Modifications for CFD Simulations of Fixed-Bed Reactors: Cylindrical Particles. Industrial & Engineering Chemistry Research, 2017. 56(1): p. 87-99.

这个工具，https://github.com/BurntSushi/ripgrep ，可以全文搜索。比如搜一下 fluent 的 src 文件夹，可以看到这两个宏的定义：
CUB

TP_N

百度云下载不了呀，太慢了。能不能上传坚果云啥的，我可以跑个测试

有一个命令行工具，叫 rg，用来在指定路径内作文本全文搜索的，https://github.com/BurntSushi/ripgrep
比方说，你这么用，是不是很快就能找到文件的定义在哪

@walten excel当然也可以，但是如果你有很多数据，一般都是编程处理比较方便。

@东岳 那老师对这个问题有什么看法呢？standard k omega 在 Re=3000 时并不能表现出网格无关性，但相同的条件，在 Re=10000 下可以表现出网格无关性。

@Juanyong 图是 MATLAB 绘制的，因为可编程才是最灵活的。参见我的文章，https://zhuanlan.zhihu.com/p/85068748

@西湖冷月 太棒了，老哥！

发散了。
一般来说，都是你 return 的那一项太大了。可以先估算 return 的那一项的大小，然后按数量级缩小，直至不发散。

吓一跳，我还以为手算数值解

@random_ran 感谢老师的意见！

@random_ran 我还有一个想法。书上都说 standard k-omega 模型是一个高雷诺数模型。我图上测试的是 Re=3000 的结果，无法达到网格无关解。我另外测试了一个 Re=16500 的结果，相同的网格，居然可以达到网格无关解。

@random_ran 老师您可以给一些意见吗？

@东岳 $l=0.07L/C_{\mu}^{3/4}$。$L$是射流喷嘴的高度。

@东岳 改成什么呀？

只有在可以设置 flux 的边界条件上，才可以指定梯度边界。
比如，温度的边界条件可以指定 heat flux，这时你就可以用 DEFINE_PROFILEDEFINE_PROFILE 定义一个值，挂载在 flux 的边界条件上。
UDS 也可以指定 flux，同理可以用 UDF 指定梯度边界条件。

各位老师，大家好！我用标准的 k-omega 模型计算了一个平面自由射流（plane free jet）的案例，分别采用了 4 种不同的网格计算。
网格信息如下，


边界条件使用 fluent 常用的设置方式。

离散格式如图，

计算的残差，mesh 1 和 mesh 2 的 continuity 在 3e-3 左右，其他物理量的残差各异，且无法小于 1e-6。mesh 3 和 mesh 4 的 continuity 在 3e-6 到 6e-6 左右，其他物理量的残差小于 1e-6。
但是这 4 种网格的计算结果并没有显示出网格无关性。结果如图所示，

继续加密的 mesh 5 也没有呈现网格无关性。而 k-epsilon 模型非常容易达到网格无关性，mesh 1 到 mesh 4 的结果几乎没有区别。
请教各位熟悉湍流模型的老师，k-omega 模型如何达到网格无关性。

用 DEFINE_SOURCE 宏，根据需求挂载在 k 或者 epsilon 的方程后面，试试

@翼型优化 抱歉。我知道怎么划分 mesh，我是想问如何绘制 mesh 的矢量图。

您好！最近我在一篇文章上看到一张网格图，是矢量图。我非常好奇他是如何绘制的，不知您的这个脚本是否有帮助，能应用到这个问题上？

• Mondal T, Das M K, Guha A. Periodic vortex shedding phenomenon for various separation distances between two plane turbulent parallel jets [J]. International Journal of Heat and Mass Transfer, 2016,99: 576-88.

@东岳 我用 1 阶精度算湍流变量。发现 Re=3000 时匹配挺好的，Re=16500 不太好，前端还是有个突起。我也不是太懂，有待进一步研究了。

@东岳 k-omega 如下图所示。我测试过网格的影响，对湍动能和速度分布的影响很小，离散方法应该只影响计算精度，不会改变趋势吧。还有待测试。感谢李老师。

@东岳 用局部速度作为参考速度，结果大概是这样。

@东岳 速度的结果很不错，能和实验对上。湍流强度呢，我用湍动能来计算的，结果偏小。
I 是湍流强度，k 是湍动能。参考速度选的入口速度，v，w为0，所以下式被我简化成了 $\sqrt{2k/3}/u_in$，不知道是否正确。

@东岳 经您这么一提醒，我搜了一下 fluent 里面是如何计算湍流强度的。原来用到了 参考速度。感谢，问题大概是已经解决了。
https://www.cfd-online.com/Forums/fluent/37012-turbulence-intensity.html

大家好！我用标准的 k-epsilon 模型模拟了一个 plane free jet 的案例，Re=3000

结果的速度场看起来很不错。

但是，中心线的湍流强度分布和湍动能却非常离谱。比如，实验的湍流强度是先上升，然后下降并趋于一个渐近值。

k-epsilon 模型所得的湍流强度的趋势也是先上升后下降，但大小就很离谱了。

我用的边界条件如下：
进口：速度入口，按下式指定湍动能 k 和耗散率 epsilon。

侧边和右侧出口指定为压力出口。喷嘴两侧的边界为无滑移壁。