CFD中文网

@bestucan 感谢您的回复，抱歉之前不知道为什么没看到。好几年过去了，今天刚看到这个之前问的问题，我那个卡住的原因是：（拿好笔记本记好了，这个问题值100万）我刚开始不会用openfoam，然后硬盘满了，所以算不下去了。

哥们，找到相关资料了吗，求教

sed -e 's|^mirrorlist=|#mirrorlist=|g' \ -e 's|^#baseurl=http://mirror.centos.org|baseurl=http://mirrors.tuna.tsinghua.edu.cn|g' \ -i.bak \ /etc/yum.repos.d/CentOS-*.repo yum clean all && yum makecache

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

再次感谢大家！

@李东岳 跟李老师汇报一下，pisoFoam可以使用fvOptions并且有了结果，非常感谢李老师！！！

从特征来看，是双曲形方程。是步进形式的问题。

理解了，谢谢李老师

各位大佬，想请教一下问题，最近在用fluent计算NASA stage这个压气机，正在计算特性曲线，但是计算结果相当不好，此外，在fluent中转静子交界面稳态情况下怎么处理才恰当啊（我的网格使用numeca画完导出来的）

@李东岳
已经看到了，非常感谢您:xiexie:

如果a/b 是一条曲线， b有的时候为0

要根治的话，就是让b不为0

tmp<fv::convectionScheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)") ) ); { reaction->correct(); dQ = reaction->dQ(); label inertIndex = -1; volScalarField Yt(0.0*Y[0]); forAll(Y, i) { if (Y[i].name() != inertSpecie) { volScalarField& Yi = Y[i]; fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi)

这个看起来并不会导致区别。有可能cfd-online的用户没做特别仔细的判定吧。

请问fluent和edem耦合时：
1.颗粒相可不可以带自热的化学反应？
2.可不可以使用动网格，网格边界发生位移的那种？
如果都不能，那么fluent和rocky耦合可以做到吗？

有时候网格数量多了反而还偏离了

有时候是会这样。也不是说网格数量越少结果越好。理论跟实践的区别，也不太好解释。

你是不是想把曲面展开为一个平面？https://blog.csdn.net/mrbaolong/article/details/106949340

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

你这个进出口边界有点问题 最后能使用固定值速度进口 固定值压力出口测试一下 瘦脸之后再调试别的边界条件
另外我看你的网格 最好网格质量能提高一下 进口那个地方看起来很奇怪

@卡洛 请问你是用的哪种方法？

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

可以啊，写个多计算域求解器就行了

各位老师，之前用的是compressibleInterFoam求解器来计算横射流雾化，后面以期做到燃烧，但由于水平有限，在添加组分方程时便遇到许多问题。
现在发现reactingTwoPhaseEulerFoam，其本身就有YEqns.H，并且据我现在初步了解他可以通过选定reactingPhaseModel相模型来计算反应。

This model represents a phase with multiple species and volumetric reactions.

想问一下如果用reactingTwoPhaseEulerFoam可以做液体横射流雾化燃烧嘛？