~17 分钟阅读
在曲线网格 FVM 里，把单元中心雅可比当作体积会让 ∂x/∂x 偏离 1。只有对角线叉积面积与直边面向量围成的多边形一致。
在曲线坐标下把应变写成协变分量，本构张量也必须一起变换。漏掉这一步，扭曲单元的刚度就整个算错。
在默认的波浪网格上，J落在0.62到1.38之间，符号一直为正。把通量直接写成y_η之后，自由来流残差会降到约1e-16。
2P/r 项离散不当,静止气体就会从恒星中心漏出。
相关标签