AeroFoamV1-2008 移植到OF41

[程迪]

如题，从https://www.aero.polimi.it/freecase/上搞来的，所有算例都至少可以正常串行运行了。

https://github.com/chengdi123000/AeroFoamV1-2008/releases/tag/v1.0-beta-of41

感觉：of1x和of41的api差异太大了。意呆利人的代码写得，奇奇怪怪的。虽然他们号称能和code-aster一起算气弹问题，但狗x的code-aster里面是法语的关键字和函数名。

[李东岳]

AeroFoam主要用来算什么的？

[程迪]

@李东岳
不知道，看算例空动问题和气弹问题为主吧。

比较有趣的是它实现了Riemann边界条件和扩展的fvm插值模板（每个面的插值用到了四个cell）。

[赵一铭]

2008年的。。。我还在上高中！。 :big_mouth:

[李东岳]

https://www.google.de/url?sa=t&rct=j&q=&esrc=s&source=web&cd=2&cad=rja&uact=8&ved=0ahUKEwjlx4z68IfVAhUFaVAKHSfbDfwQFgg2MAE&url=https%3A%2F%2Fhome.aero.polimi.it%2Ffreecase%2F%3Fdownload%3DRomanelliSerioliMantegazza_3OFWorkshop_Presentation.pdf&usg=AFQjCNFrV8c7mCWi0NQnG9WojVYUsOD75w

简单看了一下，看起来像是另一个基于Godunov-type的求解器。

[李东岳]

http://www.enu.kz/repository/2010/AIAA-2010-176.pdf

[程迪]

@李东岳 可压有几个不是基于godunov的。

[李东岳]

是的，大可不用godunov，比如常规求解。

[程迪]

@李东岳 啥是常规解法，

[李东岳]

就是求解能量，动量，连续性方程，没啥附加的了。

[程迪]

@李东岳

不都是解的这些么？

可能我所认识的Godunov有两种东西吧：

• Godunov格式，通量取激波管问题精确解
• Godunov方法，重构-进化-投影的流程搞法。

[李东岳]

以rhoCentralFoam举例，他的描述是Density-based compressible flow solver based on central-upwind schemes of Kurganov and Tadmor，应该是你说的重构-进化-投影的流程搞法，常规解法就是Density-based compressible flow solver，不需要通量分裂，或许是你说的重构 :joking:

[程迪]

@李东岳

说得跟白马非马一样。。。

FVM里我的理解是重构=插值，FVM等价于二阶Modal DG

rCF 的KT/KNP格式其实也就是进化evolution步，总之是要算出数值通量的那一步。

而且rCF的其实不是教科书上那种密度基，可以算是分离的，segregated的密度基，非常interesting。

你看rCF的代码，丫的先算出密度，在算动量，然后用动量除以密度算速度，最后算总能，用总能除以密度减去动能得比内能，再和密度一起算压力。

而AeroFoam，DBT求解器都是先算出守恒量再统一算出原始变量的。

其实二者在密度和动量方程上差不多，但是在能量方程那一项和粘性项差别就很大了。AF和DBT的代码n+1时刻的变量只直接依赖于n时刻的变量。但是rCF的n+1时刻的密度和速度只依赖于n时刻的值，n+1时刻的能量和粘性项直接依赖的是n+1时刻的速度和密度。

这种搞法不知道对不对，但是似乎大概信息传递能快一点点。代价是不一致性，也就是rCF在n+1时刻的守恒量和原始变量是不对应的。

[李东岳]

感觉：of1x和of41的api差异太大了。意呆利人的代码写得，奇奇怪怪的。虽然他们号称能和code-aster一起算气弹问题，但狗x的code-aster里面是法语的关键字和函数名。

我刚下载了aeroFoam的代码，我很吃惊。

[程迪]

@李东岳
其实AF这种才是正常的可压代码。rCF感觉能算，但是它的求解顺序不是很正统，不过也不能算太错，除非变物性大时间步。

[李东岳]

更吃惊的是aeroFoam的代码风格。

[程迪]

还好呀，你会freestyle么？

😂

[赵一铭]

😂😂
这个表情怎么打出来的 :-) :l

[李东岳]

😂 :expressionless:

[流卜]

请教下要把这个文件放在applications文件夹里编译吗