这方面工作目前读到的主要有Stanford的Matringe和University of Bergen的Hagland。
第一种方法源自于Pollock方法，最初Pollock方法严格说来只能用于rectangular elements和diagonal permeability tensors，经过几何影射可以用于quadrilateral。这一方法的核心在于velocity interpolation inside the element。简单的一阶线性的interpolation满足divergence=0的性质，但会有较大误差。
第二种方法来自Mixed FEM。这应该是一种更为自然的方法。MFEM同时解出了Pressure和velocity。使用Raviart-Thomas (RT0) 或 Brezzi-Douglas-Marini (BDM1) element可以直接进行streamline tracing。对三维问题理论上会有难度。
就现行的商业FD simulator来说，只有第一种方法可行。第二种方法必须要基于FEM的simulator才行。但实际上对第一种方法而言，velocity interpolation的根基应该还是RT0和BDM1。