我正在阅读有关计算未知变量值的快速方案ϕϕ在有限体积法中。给定一个局部一维流,我们假设ϕϕ计算为二阶多项式: ϕ=k0+k1x+k2x2ϕ=k0+k1x+k2x2
受制于:
对于均匀网格的情况,单元面值减少为: ϕϕCCffϕf=ϕC+ϕD2−ϕD−2ϕC+ϕU8=34ϕC−18ϕW+38ϕEϕf=ϕC+ϕD2−ϕD−2ϕC+ϕU8=34ϕC−18ϕW+38ϕE
我试图了解上述简化公式是如何从二阶多项式推导出来的,有什么指针吗?
QUICK 利用两个迎风节点和以及顺风节点处进行二次插值。利用二次的拉格朗日多项式形式:xUxUxCxCxDxDxfxf
ϕ(xf)=[−(xf−xC)(xD−xf)(xC−xU)(xD−xU)]ϕU+[(xf−xU)(xD−xf)(xC−xU)(xD−xC)]ϕC+[(xf−xU)(xf−xC)(xD−xU)(xD−xC)]ϕDϕ(xf)=[−(xf−xC)(xD−xf)(xC−xU)(xD−xU)]ϕU+[(xf−xU)(xD−xf)(xC−xU)(xD−xC)]ϕC+[(xf−xU)(xf−xC)(xD−xU)(xD−xC)]ϕD
注意 ,等等,你的均匀网格模板很容易恢复。xC−xU=ΔxxC−xU=Δxxf−xC=Δx/2xf−xC=Δx/2