CFD中文网

可以找已有的代码对比一下，除非特别新的东西，常用的东西一般都有模板，还有各种代码的模板。

http://www.pudn.com/Download/item/id/3745439.html
https://github.com/Singh-T/TDMA-2D

这些八成是对的，把他们的代码一点一点的替换成你的代码，看看替换到哪就不对了:142:

就我所知，至少3种方法：1.别名；2.将paraview的路径添加到PATH；3.创建快捷方式

补充一下，不知道这些图中的这些数据点是如何提取的，谢谢大家

在思考一个审稿人的意见，太难了

情况描述：一个叶片模型；操作压力设置为0，来流使用Pressure-inlet入口，总压9027.6Pa，总温299K；出口使用Pressure-outlet，使用静压出口。上下壁面使用平移周期边界，在壁面上的静压结果下壁面前缘发现高于9027.6Pa，来流Ma数大概0.3左右。
总压云图：272e8ab6-c0ca-40a6-a8ae-85289baeec81-b6de291b644f1eae2029ad5d035c97c.png
壁面静压输出结果：71799acb-1afc-4de1-9b0b-92d1214bef24-eb82740d448428364bb1ed13c86ea1d.jpg
虽然差距很小，但是却比来流总压还要高；
出问题的位置：
aaeda908-0e6a-4dd9-b51c-6c1a826e3282-f3816906934847d300a0cdfbdf0f738.png
还有很好奇的一点，理想气体的总压-静压公式应该是和Ma数相关的，在壁面上Ma=0，按理说总压应该等于静压，但是实际上并不是；
计算公式：839a2402-39a1-47b7-822c-dceda0b70b88-1f796c8655497c12e0a5455dd059f99.png
壁面总静压差△P：
d4869df2-580e-4051-ac51-1341ca071771-image.png
所以现在很疑惑。

@mohui 感谢你的解释

这个主要是监控连续性误差，累积的连续性误差。你这个结果没有守恒。这个不是发散的原因，而是发散的表象。应该其他量也存在问题。

@硫酸亚铜 在 想请教一个wallfunction的y_plus的一个问题 中说：

30

$\k$-$\omega$要求在1吗？我在a introduction to cfd 那本书里面看到，都是采用壁面函数的方法啊，而且omega也有相应的壁面模型，如果全解析到1不是LES了吗？

不好意思，那个算例里面求解器我这面编译不了（openfoam-4.x and 5.x）.. 估计后来问题已经解决了吧

@simthere 好久没有画结构网格了，回归下，练练手，暂时还能画出来，哈哈0_1503487106341_l1.png 0_1503487116733_l2.png

thinkpad t480s i7 8550u 16G内存
平时写论文、做ppt用，也可以开虚拟机改求解器算算小规模问题，真要数值计算还得远程服务器:140:

ISB

https://www.isb.bj.edu.cn/admissions/life-in-beijing

不接受中国学生

BSB

https://www.nordangliaeducation.com/en/our-schools/beijing/shunyi/admissions/how-to-apply

暂时不清楚是否接受中国学生

德威

https://beijing.dulwich.org/zh/admissions

不接受中国学生

哈罗

https://www.harrowbeijing.cn/academic-cn/lower-school

貌似接受中国学生

乐成

https://bcis.openapply.cn/

接收中国学生

青苗

https://www.bibs.com.cn/chinese-home/admissions/faq1

招收中国学生

鼎石

https://www.keystoneacademy.cn/

招收中国学生

海嘉

http://www.bibachina.org/index.html?deptId=103

招收中国学生

世青

https://www.ibwya.net/

招收中国学生

@李东岳 车灯热分析。

话说，这个作者我在CFB11的时候还见过，和她还有张合影。

你这一会1250 一会1350 到底多少度的，而且你这个截面是截的多高处的呢，另外得看一下你的残差图

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

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

再次感谢大家！

谢谢老师~~

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