CFD中文网

@evensun 求解器只负责求解数学式（控制方程），为什么要求解这个控制方程而不是别的控制方程是用户的选择，或者说取决于用户需要的是什么。
如果你需要百分百精确的描述一个非定常的物理问题，那你用稳态求解器原则上来讲是不合理的，但这个不合理不是说解方程这一步不对，而是你没有求解与你的需求相匹配的方程。如果你只需要部分描述这个物理问题，比如你关注时间平均量，那这时候用RANS去解这个问题，这一套下来是合理的。因为你需要的只是时均量，只是这个物理问题的一部分。
从这个角度去思考，你用稳态求解器计算了一个非定常的过程，求解器没有发散，甚至可能还体现出了一些有规律的震荡。这个现象本身其实是有意义的，但你在分析它的时候要非常小心。比如你认为“算得的结果呈现出了涡的发展变化过程”，这个表述是值得商榷的。

再啰嗦一点，关于RANS我一直觉得有个很有意思的角度。我们都知道雷诺数是决定流动的重要因素，我们可以认为控制方程也有一个“数值雷诺数”，它和普通的雷诺数定义类似，只是粘度项应该替换为分子粘度+数值粘度。在RANS中，数值粘度一般高出分子粘度2个量级甚至3个量级，所以尽管物理问题可能雷诺数很大，但数值雷诺数可能只在10的3次方或者4次方量级。10的3次方已经和我们通常认为的临界雷诺数接近了，导致RANS结果看起来是稳定的，就好像层流一样。
说白了，RANS模型对湍流的耗散非常强，把很多非定常的结构都给耗散掉了，留下了定常的部分。但这个耗散能力是有限度的，如果你求解的这个物理问题经过RANS模型后它的数值雷诺数正好处于临界雷诺数附近，那可能就是会有一些小的震荡。
这个湍流模型是自己开发的吗？你可以通过尝试不同的湍流粘度来验证一下，比如量级直接调大10倍，你说的这些震荡还存在吗？

@veen 谢谢，我学习一下

描述热扩散的无量纲物理量热扩散因子究竟是什么？有无针对不同情形下的具体公式？
1.png 2.png 3.png

$\alpha$

在某个时期怀疑人生很正常

找工作被拒很正常

找工作跟找对象一样

不只是我，你，也会看到很多人很流弊，工作没怎么找好。或者很多人大学谈恋爱过去了，毕业找了个年薪40万的外企。同样，隔壁老五老实巴交家穷，被富婆看上了，从此人生巅峰。你勤奋读书考研考博，现在还是单身无人问津。

以我自身举例，我总感觉自己很流弊，找工作都是没人要的，这个具体就不细说了。

所以很多时候，努力只是一部分，你需要机遇。机遇只能靠等。

做好定位再打算

如果确定要出国、读博。你英语水平还很好。可以一直找下去。要么是“实验室已满”，要么是“我没有钱如果你自费我们再考虑你”，这种太正常了。我经常遇到这种情况。这说明不了什么。换别的。硕士期间还没有定位研究方向，可以满世界找。

在这里我强调一下，就是不一定要找顶级高校顶级课题组。以我自身举例，2014年出国，我去的是一个世界排名400的学校，导师那个时候还是副教授。但人非常好，我学到了很多东西。很多，包括国内，最强产出者很多来自于讲师与博士后。相反，大教授往往集中精力与拿项目组团队，具体的技术都是博后来做。

所以扩大选择面。

人生三年一变

你回头看你3年前在做什么，6年前在做什么。所以，你才26岁，机会多的很。我2018年1月回国后，还失业了半年，也没影响我做CFD。

给你个2018年我申请北京某高校的拒信截图

捕获.JPG

论坛最近迁移服务器，才发现搜索功能丢了，我去找找，铁铁
先用 @random_ran 的方法

@yuanliangwojt 对的，是动网格的时候会有，我还以为你这边也是动网格

计算两个并排的小球在水中下降，关注两个小球的中心位置变化。物理原型为船舶紧急抛双锚。

目前在of2206中设置了两个小球，但设置sixDoFRigidBodyMotion时，发现貌似不支持为每个球单独设置运动参数。比如两个球大小、密度不一。

上传算例，哪位能帮助看看如何解决？falling.SphereTwoNew.zip

Screenshot from 2025-01-10 01-53-59.png

kOmegaSST和kEpsilon印象中是比较好收敛的两个模型，其他的用起来似乎都不太顺手，当然最重要的是结果正确。

@李东岳

李老师您好。

目前出口是采用如下方程
3.png

这里Uc尝试过：

使用上下层速度的平均值1.5。 在X方向上采集最靠近边界的三个点的速度值进行平均，作为Uc。

如果采用2的话，Uc是一个随着Y坐标以及时间变换的函数Uc(y,t)

Phi是压力泊松方程构造的中间量，这个查了好多文章都说的很模糊，我也就直接用了零梯度。

@纵浪大化 好的，常常交流联系

在利用per phase求解湍流方程时，可以用以下语句获取各相指针。
Domain *dp = Get_Domain(2); /primary phase domain pointer/
Domain *ds = Get_Domain(3) ;/secondary phase domain pointer/
Thread *tp = Lookup_Thread(dp, 7);
Thread *ts = Lookup_Thread(ds, 7);
2,3分别是各相domain_id，7为区域ID，可以从用户图形界面获得

@李东岳 李老师有没有兴趣看下风工程这块的有个商业软件WT，里面是速度-压力耦合求解器，真的非常快非常快。我看了下，软件底层是从PHOENICS来的，然后搞了一套叫MIGAL的算法。

现在CFD用的SIMPLE这种segregated solver确实太慢了，城市这种尺度下的流动本身大概率又比较简单，不可压，也没有化学反应之类的，可以暂时不考虑温度？可能会有颗粒流的问题，但是不是也可以先只考虑流动对颗粒的影响。这种情况下感觉流场本身的求解其实就看求解器的效率了。

我搜了下，OpenFOAM这边其实已经有一些coupled solver 的工作，不知道成不成熟。

@西湖冷月
nice兄弟，d=====(￣▽￣*)b

openfoam的多孔介质求解器，无法获得真实的速度

多孔介质大部分都是层流，不需要你手动去控制湍流以及层流

请问燃料燃烧时这种‘理论计算产生的平均总温’应该怎么计算啊？翻了基本燃烧学的书都没这方面的定义。温升燃烧效率.png

@jinlinna 感谢！

这种情况宁可相信Tecplot中Calculate variables和自己编程计算的结果，Fluent导出的速度场是没有问题的。

@ZZ给我趴下 在 Smagorinsky模型不输出k_sgs场 中说：

@李东岳这里提示重复定义，在59行有定义，而且和这个函数一样。 003f24ca251de51d6dbeaacef31c1ce.png

老帖子了，为了方便后面的人看，这里指出下问题。原始的代码159行这里是调用correctNut，而它的实现是在上面的59行。李老师的意思是换掉59行位置correctNut的函数实现代码里面的k_，不是把此处调用correctNut换成另一个新的实现，因此会报重复定义错误

@hxy970122 扫了一眼Fluent的组分输运方程，对于各个组分就没有重力的影响，所以如果组分分离就肯定是因为扩散率不同。但是在binary mixture里面，似乎两个组分的扩散系数是相同的？

重力的影响是在质量方程和动量方程里面，所以就是所有组分受到的重力作用是平均重力……

2e37529c-113a-4d19-a63f-e23bf85c3593-image.png