AeroFoamV1-2008 移植到OF41


  • 管理员

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


  • OpenFOAM讲师

    @李东岳

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

    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的代码,我很吃惊。


  • OpenFOAM讲师

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


  • 管理员

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


  • OpenFOAM讲师

    还好呀,你会freestyle么?

    😂



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


  • 管理员

    😂 :expressionless:



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



  • 小白试了你更新的代码,运行时出现下面的错误,请问是我编译位置错了吗0_1521630956325_error.png 图片地址)


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待