李新亮老师opencfd程序交流

gyzhangqm

最近开始使用李老师程序计算多块网格问题，发现网格连接处物理量不连续，求助！

---------------------------------自己更新----------------------------------------
原因：是后处理的问题，体心有限体积法变量存在体心，在网格块棱角附近，是靠插值得到的数据，所以tecplot显示的时候会出现不连续。
解决办法：后处理时对于棱角处并非采用八个网格格心平均得到，而是棱处采用六个网格平均，顶点处采用四个网格平均。
---------------------------------自己更新----------------------------------------

发现这样处理也不连续
---------------------------------自己更新----------------------------------------
我在李新亮老师三维opencfd-ec的基础上，添加了平衡流，化学非平衡流，现在正在尝试热化学非平衡流。李老师opencfd-comb程序使用刘君解耦的思想，易于在已有的程序上通过gamma的变化将反应的影响加入流场。我使用的是耦合的方法，也算对李老师程序比较熟悉了，相关领域的同学可以交流~
---------------------------------自己更新----------------------------------------
我对CFD学习时间比较短，感觉李老师的程序使用fortran90 ,差不多是有模块化的编程。我不知道大家有没有用Fortarn03的（或者C++）面向对象编程的。感觉这种面向过程编程的思想，只能在小作坊里发展，不容易形成团队合作。我现在闲时在学一个开源代码OFF，有同样在学习的可以交流。

gyzhangqm

另外，李老师opencfd-ec程序早期版本处理边界条件采用双层ghost cell，现在改成单层，在边界处计算无粘通量使用中心格式，我不懂为什么要这么做？

Briareos

gyzhangqm 发表于 2017-3-17 20:19
另外，李老师opencfd-ec程序早期版本处理边界条件采用双层ghost cell，现在改成单层，在边界处计算无粘通量 ...

在边界处通常网格要进行加密，尤其是壁面边界处网格雷诺数会非常小，这时采用中心格式会避免因为限制器的使用而造成的数值振荡，当然这也是我从别的文献上看来的，说的有不对的地方还请指正。

plyutjuzjuksu

gyzhangqm 发表于 2017-3-17 20:19
另外，李老师opencfd-ec程序早期版本处理边界条件采用双层ghost cell，现在改成单层，在边界处计算无粘通量 ...

你好，请问李老师的代码三维的怎么获得？

ch06

后处理问题，是因为程序是有限体积法，数据存在网格中心，而tecplot是基于节点显示的，网格中心的值跑到了节点上，看着自然就有问题了，尤其是在有曲率的边界上。

ch06

这个问题是硫酸亚铜网友发现的：
“fluent的计算数据是存储在每个单元的中心位置，而tecplot的云图显示是按照节点来显示的。如果是直接导入fluent的.cas文件和.dat文件的话，那么在tecplot中，节点上的值会通过单元中心点的值插值得到，从而造成壁面速度不为0的情况。”

boyzone

请问这个多块网格文件格式是怎样的？网格格式包含哪些量，表示啥含义？

匡之舒云

请问现在还有老师在用OpenCFD-Comb这套程序吗，可以互相交流讨论一下吗？