CFD中文网

@李东岳 谢谢李老师的回复，我最后通过转GMV格式，然后导入Paraview解决了，感觉是读取文件大小的限制，直接读.vtu还是会出错。

最近用pointwise中的非结构网格生成了一个风力机旋转域中的block，检查质量发现最小角度0.14，网格质量特别差，但目前只会调整domain的网格质量，block的没看到相关教程，有大神可以指导一下吗？谢谢！

@两月三年 是的，这样显示的时候是表示选上了。注意选的位置就好。

请问大佬们知道怎么用starccm实现传送带（带凸起）的运动吗。
目前能想到的就只有给凸起物一个个建模设置好运动路径，但应该有更简单的方法吧。提前感谢各位大佬。
0ae2d1b3-d3ba-4382-8cac-a470e7232877-image.png

这个是理论大坑啊。

如果压力基的求解器，求解的是速度和压力场，那么得用连续性方程来求压力。
如果是密度基的，连续性可以求密度，然后用状态方程求压力。

@Vortex
感谢学长^_^

计算域某个区域可以吸收CO，吸收的CO的量与该区域的每个网格单元本身的浓度值成正比，请问各位老师这个UDF对吗，我设置了这个之后CO浓度无变化

#include "udf.h"

DEFINE_SOURCE(co_absorption, c, t, dS, eqn)
{
real S;
real Ci;

Ci = C_YI(c, t, 0); // 获取CO的浓度，假设CO是第一个组分（编号为0） // 计算吸收区域吸收CO的量 S = -0.0064 * Ci; // 更新源项数组 dS[eqn] = S; return S;

}

在用tecplot360画速度矢量图时，箭头分布比较密，想把它调稀疏点，调point只能改变I方向上的点数，J,K方向的无法调，请问有遇到类似的吗？

@小葱小虫
是第三方的，要使用foam-extend才能安装。
transonicMRFDyMFoam可以计算稳态

基于2011年发表在Journal of Computational Physics的文章Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations，并根据您文章的网格自适应算法编写了一套计算程序，在计算双马赫反射算例过程中，遇到了一些问题，受困许久，希望大神指点迷津~ 在此诚挚感谢
遇到的问题具体问题如下：
在计算双马赫反射问题时，在Kelvin-Helmholtz rolls两侧流场出现较大的震荡，局部密度云图及网格自适应加密情况如图1所示，震荡区域存在粗细网格过渡,故怀疑在粗细网格过渡中存在算法错误或者代码错误（我已多次检查代码，未发现代码错误，但不能完全排除仍有代码错误存在）。
0_1536922532180_9e8d3f5c-b1c2-4828-bec3-8e472986a4d2-image.png
（a）局部密度云图（红框处为存在流场震荡处）
0_1536922552890_10e26979-d84c-4d89-956b-e81584a6adb5-image.png
（b）局部网格生成
0_1536922574620_04fd03c7-eca8-4851-a7e4-3c342cae3939-image.png
（c）总体密度云图
0_1536922591916_530a5e07-8d83-41c6-b1e3-338c5acb0b3a-image.png
（d）总体网格生成（黑色区域表示网格加密）
图1 局部密度云图及网格自适应加密情况
该程序采用C++语言编写，流场计算使用有限差分法，五阶WENO格式，Steger-Warming通量分裂，2D欧拉方程，时间推进采用3阶Runge-Kutta方法；网格自适应加密方法和加密准则与您的论文有所不同，采用基于原始根网格进行加密，每个网格单元生成3×3个子网格单元，加密准则使用密度梯度(∂ρ⁄(∂x+∂ρ⁄∂y))。一共使用三层网格，在计算域[0,3]×[0,1]中，原始根网格单元数量为300×100。Layer 0表示原始根网格，Layer 2表示最密的网格。Layer 2比Layer 1略大，覆盖Layer 1，两层网格的间隙为3个根网格单元的尺度，如图2所示。以此希望精准捕捉激波，并适当减少网格生成数量。子网格层的时间步长为父网格层的1/3。
父网格层每计算完一步，对父网格节守恒量（ρ，ρu，ρv，E）和守恒量对时间的导数[ρ’，（ρu）’，(ρv)’，E’]进行五阶加权(类似WENO加权)空间插值，求解出子网格层边界的虚网格节点的相应守恒量及其导数。并将前一时间步和本时间步插值出的通量及通量导数使用Hermite插值，以获得三个子时间步的相应变量，作为加密网格层的边界条件，并使用式（1）作为Runge-Kutta时间推进的中间步的边界条件。
0_1536922724975_5a02e4bb-a1eb-4fce-884b-9d73584cc27c-image.png
在子网格层三个时间步计算完成后，父网格层被子网格层覆盖的节点数值被子网格对应同位置节点更新。每一或两个根网格时间步重新生成一次加密网格。
0_1536922743488_5f7e9dc2-9e8c-44f9-88fb-73083e652ec1-image.png
（a）粗网格和细网格
0_1536922779864_ecb19b47-59e2-4210-97fb-4d0f9de2d195-image.png
b）不同层次网格
图2 网格加密示意图
其他相关的情况介绍：
在一次测试中使网格只加密，不减少，结果显示流场基本没有前述震荡，结果如图3所示。此算例应该可以证明该程序计算格式是正确的。然而，由于生成网格实在太多，不具实际应用意义。
0_1536922795688_29716dff-e9d7-4afb-a475-b7a750cebe66-image.png
（a）密度云图
0_1536922809856_e6bd76e2-af3b-4966-9f15-c3e59d18ea3b-image.png
（b）网格生成情况（黑色为加密区域）
图3 网格只加密不减少的计算结果

不知大家在调试程序过程中，是否遇到过类似的情况。
希望大家能够指出
（1）我对算法的理解是否存在错误，特别是粗细网格过渡部分的算法；
（2）解决该问题有何建议；
（3）能够告知粗细网格过渡区域，造成的流场震荡的原因可能有哪些;

再次感谢大家！

@东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了，算的也是温度梯度，也就是$q=-\kappa \frac{\partial T}{\partial n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度，这种差值不知道在of里面能不能实现

最近在用ICEM生成非结构化网格，网格数量大、八叉树生成方法慢，且开始时几何总有各种问题，这个过程中ICEM信息窗格一直在报告进程。在发现问题和解决问题的过程中，我才意识到信息窗格的重要性。想更多了解信息窗格的阅读tips，但没有找到相关资料，希望在这里能和大家交流。相关问题大家踊跃跟楼~hh

现在我有个问题是，信息窗格字体的不同颜色分别代表什么？比如蓝色字体代表什么？

Ai做精巧定位可以，流体做仿真测试也可以，但二合一还真没见过。

额，自问自答？

不知道大家感受到没有。换了一个开源的搜索引擎，叫meilisearch。简直屌爆了。。。。。

:146: :146: :146:

@东岳 如此说来，除了加密网格，没有好的解决方案啦:zoule: :zoule:

@hurricane007 确实礼貌又得体。

利用reactingFoam计算正癸烷在水平圆管中的裂解反应（18组分生成物）、热物理模型采用
0_1492507996639_upload-7ce04876-1271-4e49-b10f-e26cf46d3989
计算结果中出口速度偏小、正癸烷在圆管中的停留时间过长，导致出口处的正癸烷质量分数偏小，考虑到与其密度有关；因为正癸烷裂解反应是在超临界条件下进行的，算例中采用的热物理模型不能满足计算要求，故考虑采用rhoReactingFoam求解器，热物理模型采用
0_1492508023620_upload-b90bd2f7-b8c7-4863-8558-4089770b0ad8
各个热物性参数通过与温度之间的多项式关系进行计算。例：
0_1492507617167_upload-2786160a-cf82-4436-b6a6-f40fb8d70e2f
但计算一步就出现错误
0_1492507781160_upload-d1e18a85-91ef-448b-a634-cacd5cda8fc0
目前不知道什么原因导致的，请问哪位大神能给指点一下？

@zhanghan 你的压力场进出口不建议设置为全部为zeroGradient。可以把出口设置为固定值。