CFD中文网

https://huaban.com/boards/15971106/

新买个房子，打算打造个看文献写代码舒适的办公环境！ :xiezuoye:

替代文字

上面这个灯搞一个

替代文字

背板可以用冲孔板

替代文字

这个看起来不够立体，应该在屏幕上方加上一些元素

替代文字

搁板下面可以加个灯带

https://www.zhihu.com/question/30550212

报错了，我检查不出错误，希望有前辈能够指导一二

x3d
用这个导出，导入到blender

@allvic 是不是可以在边界层内设监测点，得到频率分布，与T-S波的频率进行相关性验证呢？

在breaking wave领域，这种海浪结构，会将空气卷吸并形成很多小气泡。最近提交一个稿件，审稿人推荐我们讨论一下这方面的研究。我详细看了一下，确实有点意思。

看下图，一些实验研究已经证明，这种海浪拍打过程卷吸进的气泡具有尺寸分布。见过很多文献，都是用的DNS直接模拟研究的，结果也都能证实。

但目前我想的是，能否用PBE模型，或者其他类似模型做相关的研究。

捕获.JPG

气泡直径和分布的关系
捕获.JPG

一些参考文献

Deike, L., Melville, W. K. & Popinet, S. 2016. Air entrainment and bubble statistics in breaking waves. J. Fluid Mech. 801, 91–129.
Wang, Z., Yang, J., Stern, F., 2016. High-fidelity simulations of bubble, droplet and spray formation in breaking waves. J. Fluid Mech. 792, 307–327.
G. Soligo, A. Roccon, and A. Soldati, Breakage, coalescence and size distribution of surfactant-laden droplets in turbulent flow, J. Fluid Mech. 881, 244 (2019).
Ahmed, Z., Izbassarov, D., Costa, P., Muradoglu, M., Tammisola, O. 2020. Turbulent bubbly channel flows: Effects of soluble surfactant and viscoelasticity, Computers & Fluids 212, 104717.

@chenxiao
暂时还没有。最近事情非常多。:crying:

@东岳

just for fun:

They used to say if man could fly, he'd have wings. But he did fly. He discovered he had to. Do you wish that the first Apollo mission hadn't reached the moon, or that we hadn't gone on to Mars and then to the nearest star? That's like saying you wish that you still operated with scalpels and sewed your patients up with catgut like your great-great-great-great-grandfather used to. I'm in command. I could order this. But I'm not because, Doctor McCoy is right in pointing out the enormous danger potential in any contact with life and intelligence as fantastically advanced as this. But I must point out that the possibilities, the potential for knowledge and advancement is equally great. Risk. Risk is our business. That's what the CFD is all about. That's why we're aboard her. You may dissent without prejudice. Do I hear a negative vote? (silence around the table) Engineer, solving that imcompressible N.S. equations using PISO algorithm.

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

各位大佬，请问CFD-DEM coupling时，流体的控制方程组求解，用耦合式解法还是分离式解法？
非常感谢！

我模拟的对象是反应堆内的稳态流动传热，由于堆芯内组件数量庞大，精细建模复杂，因此采用多孔介质模型对堆芯区域进行简化。

为了验证这种等效的合理性和误差，我先对单个组件栅元进行多孔介质简化，并将多孔介质模拟结果与未简化的模型结果对比验证多孔介质模型的合理性。

单个组件栅元的截面图如下：
图片1.png

其二维轴对称模型如下图：
f510f5de-0d5e-4389-b577-ad2d3defe29c-image.png

模拟结果
(1)未采用多孔介质模型
我用了二维轴对称模型来模拟，设置固体流体内热源，入口质量流量和出口压力边界条件，流体和固体分别加载体积热源P1和P2（W/m3），流体和固体的导热系数均为常数，得到的稳态温度场和压损分布如下：
bca065f1-7cc4-4385-8113-d7e0acd6d23e-image.png
流体出口温度700℃，固体由于带内热源温度高于流体，固体最高温度为799℃，固体平均温度为756℃，组件压损为2400Pa

(2)采用多孔介质模型模拟
由于固体带内热源，流固存在较大温差，采用非热平衡模型
根据栅元几何结构指定：孔隙率γ=0.1 ，界面换热系数hfs=5000W/M2K（与未采用多孔介质模型算得的换热系数一致），界面面积密度Afs=6.71m-1

通过流速-压力实验关系确定粘性损失项系数和惯性损失项系数，只考虑x方向的流动损失，y方向系数放大100倍

别给多孔介质流体和多孔介质固体设置了内热源，分别为P1×孔隙率和P2×（1-孔隙率），初始化后检查内热源加载总量无误；流体固体的热导率输入值和上个模型一致

得到的模拟结果：
af1efa10-a95f-4e49-a69c-63b4b4357ca0-image.png
采用多孔介质非热平衡模型模拟得到的流体温度是正确的，压损也是正确的，但固体温度分布明显不对，固体域温度热点仅为725℃，平均温度为711℃，均远小于实际模型算得的固体温度值。

我自己的猜测：
多孔介质模型将单组件变为流固均匀介质后，不能模拟固体区域沿导热路径带来的固体温升？刚好我这个栅元模型中固体由于导热产生的温升非常高，如果是这样，说明多孔介质非热平衡模型不适用于我所研究的问题，因为简化后会显著低估固体温度热点。

请问我的理解是正确的么

@haining-luo
我现在用ansa处理极复杂三维模型，保存.mesh格式面网格，用fluent meshing或star ccm生成多面体网格或切割体网格，挺方便。
如果模型太复杂，搞科研可以用icem，工程上不太建议用ICEM生成非结构网格。

@东岳 嗷，我以为您的70g激光打印大杀器已经服役了ಥ‿ಥ

@yingqing 遇到了，我是根据错误提示去排错，然后搞定了。我遇到的大部分的报错似乎是缺少必要的文件。

个人建议，看这个之前，看看icoFoam，那个是绝对的经典~:146:

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。
但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。

所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的:135:
在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。
我把相关资料发给老师。

@shrine 没关系，无关紧要这些文件

@桥下卧龙 也没帮上忙，实在不好意思，我最近也在做类似的东西，不过没出现你这样的情况。

有的时候icem就是会报错，只能没做一些操作之后多保存，然后再进行。。

多谢楼主了，可以计算了，赞一个！