FVM - 虚拟节点离散化

计算科学 有限体积 离散化 双曲-pde 外推
2021-12-09 08:04:48

我看到了一篇题为: 一种单片流体结构交互算法......

对于一维网格,论文在边界处使用虚拟节点x0xN+1(第 372 页),例如对于边界处的速度 v,通过外推(第 374 页)结束

x0=2x1x2

xN+1=2xNxN1

我想知道是否有人可以向我解释点 N 和 0 在哪里?在面上还是在节点上?还有什么是使用这种推断的原因。

谢谢!

1个回答

有限体积法 (FVM) 使用单元平均值来表示数值解。例如,xi表示中心i-th 网格单元,以及Uin表示网格单元内的数值解i有时tn. 单元边界处的通量由i±1/2下标,例如,使用您论文的符号,Φi+1/2是在右边界的数值通量i-th 单元格(在网格坐标处xi+1/2)。

您的网格单元索引范围从1N. 除此之外,作者还定义了两个带有索引的“幽灵”单元0N+1. 这些不是您的“真实”网格配置的一部分。相反,它们的定义允许您计算网格“真实”边界处的通量(即,Φ1/2ΦN+1/2)。在 1D 中,你可以想到 cell no0立即在网格单元的“左侧”1. 同样,细胞N+1立即位于单元格的“右侧”N,就像这里给出的插图一样。

使用标准外推公式计算重影单元内的解近似值。