QUICK方案推导

计算科学 有限体积 插值
2021-12-11 01:42:52

我正在阅读有关计算未知变量值的快速方案ϕ在有限体积法中。给定一个局部一维流,我们假设ϕ计算为二阶多项式:

ϕ=k0+k1x+k2x2

受制于:

  • ϕ=ϕU (上风)x=xU
  • ϕ=ϕCx=xC
  • ϕ=ϕD (顺风)x=xD

对于均匀网格的情况,单元值减少为: ϕCf

ϕf=ϕC+ϕD2ϕD2ϕC+ϕU8=34ϕC18ϕW+38ϕE

我试图了解上述简化公式是如何从二阶多项式推导出来的,有什么指针吗?

1个回答

在此处输入图像描述

QUICK 利用两个迎风节点以及顺风节点处进行二次插值利用二次的拉格朗日多项式形式:xUxCxDxf

ϕ(xf)=[(xfxC)(xDxf)(xCxU)(xDxU)]ϕU+[(xfxU)(xDxf)(xCxU)(xDxC)]ϕC+[(xfxU)(xfxC)(xDxU)(xDxC)]ϕD

注意 ,等等,你的均匀网格模板很容易恢复。xCxU=ΔxxfxC=Δx/2