AeroFoamV1-2008 移植到OF41
-
说得跟白马非马一样。。。
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时刻的守恒量和原始变量是不对应的。