CFD中文网

网格加密到Y+<1，计算结果变化不大，对流通量用二阶MUSCL, AUSM格式

你不觉得大材小用了，浪费那么多的时间搞科研，不说为人类的进步献一份力把，至少让其学有所用啊，他的那些东西中学根本用不上。 人生价值何在啊。何必当年读博士。

才发现是Fluent板块。。。不过这功能应该不难实现

如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

@李东岳 谢谢东哥，就是我把resolveFeatrueAngle调大之后画网格没有捕捉这条特征线，画几何的时候因为在两个不同平面上，这条特征线是存在的

可以试试screen，tmux，byobu

苹果retina屏那个本也不错。就是能不能装win系统？

http://item.jd.com/2819597.html

@tidedrinker 在 关于fvSchemes字典中wallDist的疑问 中说：

哪些模型最简单的方法是空置这一项，或者故意填个错的，然后报错的时候会提示你有那些可选
确实k-omega SST需要，没有就会报错

对，然后不计算边界的话就是0

请问大家为什么我用meshing划分流体域网格时,明明两个面没有重叠，模型树上还会出现图片上这种接触区域呢？![6392L5}X9$J}C{GTK_ZZ7F.jpg

结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说，结构化和非结构化区分是文件保存格式的区分，结构化是使用索引，非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格，因为我最近也在学openfoam和icem，icem可以通过自动化分成四面体非结构网格，也可以分块化成6面体网格，但是我通常画完6面体网格后，icem有可以转化成非结构网格按钮，将画好的网格转化成非结构体网格（转化后依然是6面体）。
所以如果你要话非结构化网格也可以分块（画了6面体转化一下），如果你说的是画四面体非结构化网格分块的话，这方面知识我就不知道了

多谢老师解答，问题已经解决了

@李东岳 在 MULES::correct() 和 MULES::explicitSolve() 有什么区别？ 中说：

MULES要看基金会的版本，这是基金会那面提出的。我记得CMULS已经被删了。我在无痛苦ns方程笔记里面说的很齐全。你看看参考一下

好的，谢谢李老师

请教一下你是每一层加密都对三个方向进行加密吗，shm中能能不能只对z方向进行加密呢。@一朵苍凉

加密前网格由blockMesh生成blockMeshDict

FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // scale 1; vertices ( (-0.15242 0 -0.005) ( 0 0 -0.005) ( 0.3048 0.081670913853 -0.005) (-0.15242 0.1524 -0.005) ( 0 0.1524 -0.005) ( 0.3048 0.1524 -0.005) (-0.15242 0 0.005) ( 0 0 0.005) ( 0.3048 0.081670913853 0.005) (-0.15242 0.1524 0.005) ( 0 0.1524 0.005) ( 0.3048 0.1524 0.005) ); blocks ( hex (0 1 4 3 6 7 10 9 ) (80 80 1) simpleGrading (1 1 1) hex (1 2 5 4 7 8 11 10) (160 80 1) simpleGrading (1 1 1) ); edges ( ); boundary ( inlet { type patch; faces ( (0 6 9 3) ); } outlet { type patch; faces ( (2 5 11 8) ); } bottom { type symmetryPlane; faces ( (0 1 7 6) ); } top { type symmetryPlane; faces ( (3 9 10 4) (4 10 11 5) ); } obstacle { type patch; faces ( (1 2 8 7) ); } ); mergePatchPairs ( ); // ************************************************************************* //

加密后是由fluent.msh转换的fluent.msh
加密前：controlDict
加密后：controlDict
其他设置基本是一样的
物性参数：thermophysicalProperties
湍流模型laminarturbulenceProperties
离散格式：fvSchemes
求解器控制：fvSolution
来流条件：p T U

这个是理论大坑啊。

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

绘图1.png
李老师，上面算的算例因为网格问题要重新计算了，我参照您发我的算例在ICEM中画网格。如图所示A的严格对称点应该是B点，您画在了C点，所以要在openfoam里面用cyclicAMI来周期上下两个面。出现了面不匹配的老问题，我尝试在ICEM中设定周期，但是好像这种叶珊上下面应该没法设置节点的周期把？随后我仔细观察了您的算例，发现在boundary文件中的cyclicAMI的容差设置为1，随后我在createPatchDict采用了相同的容差。这一步算是可以通过了，但是在decomposePar中采用scotch方法进行分块时，还是在cyclicAMI这里报错了，麻烦李老师给看看我上述的操作是漏了什么嘛？
另附一下 网格处理命令：
fluentMeshToFoam ICEM.mesh
autoPatch 45 -overwrite
createPatch -overwrite

屏幕截图 2023-12-30 235244.png

求教岳神，我想做反应熔渗的仿真，涉及及多孔介质毛细现象，碳与硅液生成碳化硅的化学反应，反应后又会填充多孔区域孔隙，结果需要填充过程趋势，有好的方法做这方面的仿真吗

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

再次感谢大家！