CFD中文网

@ycx1042341217 @huningdong 请问两位，壁面第一层网格高度可以做到什么程度？我发现最小只能设置到大约5e-5m:136:

比如调成2的话，这个时候的时间步长是多少？

你的时间步的2倍。

根据这幅图的说法，隐式耦合计算时，非稳态的courant我也可以设置的比1大吗？

隐式理论上可以使用任意大的时间步，因此没有库狼数限制。所以可以。

@him 谢谢解答，另外一般商业软件的定常求解器是用的时间推进来求解的吗？如果是时间推进的定常求解器，这种不稳定的振荡可否理解为流场随时间的变化？

修改自https://www.xsim.info/articles/OpenFOAM/Switch-multiple-versions.html，致谢大牛。

要添加其他版本只需仿照5-7行进行修改即可。

将该脚本拷贝到/usr/local/bin并赋予执行权限：

sudo cp ./switchFoamVersion /usr/local/bin/ && sudo chmod +x /usr/local/bin/switchFoamVersion

切换OF版本时一定要在命令前加"."或"source"，不然环境变量不能正确地传递到当前shell。

3967c521-2e7a-46bf-81e8-d0bb3be3fe33-屏幕截图 2023-09-12 115521.png

switchFoamVersion

微信截图_20201227163302.png
如图，这个代码是怎么实现的，原理是什么

Gr的计算在G之后，那么G之前的Gr是怎么算的呢？

如果代码就是这样的。那这个方程就是显性离散的。针对他的问题Equation (14) is implicit; how is this evaluated in practice? ，你可以回答这个方程是隐性的，但是我们显性计算了。

谢谢老师，我找出来看一下

@李东岳 对的，东岳老师，应该不是版本问题，我是在Foam5版本运行的不可以

@shrine

我上面说过的，字典文件中有一个坐标，这个坐标落在几何里面，生成的网格就在里面，落在几何外面，生成的网格就在外面
66148774-6b5c-42db-a783-f365dc9e147d-image.png
2.在生成几何时就对每一个面进行分part并命名，然后就可以对part分别命名，我是用ICEM画的stl几何，里面可以进行这个操作，其他软件我不了解

5右边多写了个$\rho$吧，谢谢！回北京我更新一下

貌似罗技的鼠标相关功能做的不错，还是硬件上实现比较的靠谱。

@李东岳 李老师，可能是您python版本的问题，我用的是python3.8，您应该也有python3但是系统调用的是python2.7，可以试试python3 plot.py运行。
坐等大佬改我的代码:chouchou:

@李东岳 除了phase的图像，其他压力图，密度图等都没有问题。0_1504141697200_ea5fe613-9544-43ae-acce-f76b5a6f726e-image.png 这是密度图。我在网上有人提到会不会是由于相与相之间的界限不明显，而产生这个问题。老师，请问您觉得这种说法可靠么？谢谢

明白了，感谢各位的回复

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

数值计算方法

具体是什么，解PDE么，还是什么

基于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）能够告知粗细网格过渡区域，造成的流场震荡的原因可能有哪些;

再次感谢大家！

@wangfei9088 谢谢大佬的详细解释，很全面。很感谢！本来想着用Refprop软件来求参数了（虽然也不知怎么求），幸好大佬您出手了，再次感谢。

加油！

@bestucan :mihu: 确实是直接删除block就行，俺发现好多操作都跟人家一样。得到得结果却不一样，可能因为我用的不是正版