有关FDLBM中的插值格式

zhangyangbo1180

在FDLBM的迁移阶段，由于一个时间步长内分布函数不能到达与其相邻的网格节点，进而需要通过插值方法得到网格节点处的分布函数。
本贴插入的图片给出了一中插值格式，这个插值格式是根据二阶泰勒展开的相关理论推导出来的，很容易可以想到这种插值方式仅仅适用于格点密度梯度不是很大的情况。


本人用LBM模拟化学反应的相关问题，由反应产生的某种物质仅仅存在于少数的网格节点处，并且该物质的格子速度也小于网格步长，所以用LBM求解该物质的输运过程是，也要用插值方法得到网格节点上的分布函数。

现在的问题是，由于物质仅仅存在于某些少数节点处，上面插值公式的密度梯度不能太大的使用条件不再满足，因此不能继续使用。请问大家对于这种问题应该使用什么样的插值方法呢？

[ 本帖最后由 zhangyangbo1180 于 2011-6-9 19:55 编辑 ]