OpenFOAM中颗粒数据的插值

在欧拉-拉格朗日方法中，离散相被认为是质点，颗粒当地对应连续相的状态量一般需要通过插值的方法来获得。以OpenFOAM自带的拉格朗日类的算例simplifiedSiwek为例，在constant/coalCloud1Properties中指定不同物理量的插值方法

solution
{
......
    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        thermo:mu       cell;
        T               cell;
        Cp              cell;
        kappa           cell;
        p               cell;
        G               cell;
    }
......

关于离散相当地参数的插值方法，OpenFOAM中常用的有以下几种选项

1. cell
   Uses the cell value for any point in the cell.
   对网格内所有的离散相均采用网格当地的参数
2. cellPatchConstrained
   Uses the cell value for any point in the cell apart from a boundary face where it uses the boundary value directly. Note: will not work on an empty patch.
   这个是对cell方法的修正版，即对边界处的网格内的离散相采用边界条件的值，对其他网格内的颗粒仍采用网格当地参数
3. cellPoint
   Given cell centre values and point (vertex) values decompose into tetrahedra and linear interpolate within them.
   将网格中心的值和网格顶点的值分解为“四面体”，并在其中线性插值来得到离散相对应的参数。
4. cellPointWallModified
   Same as interpolationCellPoint, but if interpolating a wall face, uses cell centre value instead.
   是cellPoint的修正版，即对wall边界处的网格不采用插值方法，而是直接采用网格中心的参数
5. pointMVC
   Given cell centre values interpolates to vertices and uses these to do a Mean Value Coordinates interpolation.
   将网格中心的值插值到顶点处，并利用这些值对离散相对应的物理量进行平均线性插值。

我们一般常用的是cell和cellPoint方法，至于第5种pointMVC方法和cellPoint方法有什么不同，感兴趣的小伙伴可以测试一下

明天就是五一劳动节了，提前祝大家节日愉快，代码一遍work！