浸没边界与LBM.docVIP

浸入边界法的产生，发展和逐步完善为解决外形复杂的结构在流场中运动的模拟、流固耦合和运动边界问题提供了新的途径 1 绪论 对于复杂流场、流固耦合和运动边界问题，常用的两种处理方法是:基于坐标变换和网格映射技术的有限差分法和基于非结构网格的有限体积法。有限差分法需要寻找适合的高精度变换矩阵计算方法进行坐标变换和网格映射，因此即便是生成一个高质量正交性好的网格都需要耗费大量的时间进行迭代运算。而随着人们所需处理问题复杂程度的不断提高，生成理想的结构化网格系统变得异常困难。对于复杂的几何体，采用非结构网格是可行的方法，然而非结构网格的质量会随着几何体的复杂程度而明显下降，从而导致模拟失真。 针对上述问题，目前不同研究领域都已提出相应的计算方法。浸入边界法因其具有良好的发展前景，成为新的研究热点。为了模拟血液在可收缩的心脏瓣膜中流动，Charles Peskin [1]在1972 年提出了浸入边界法思想，即通过将复杂结构的边界模化成Navier-Stokes动量方程的力源项，当这些力加在特定的一些网格点上时，可以成功的模拟出任意形状的结构边界。 浸入边界法(Immersed Boundary Method)，既是一种数学建模方法，又是一种数值离散方法。在数学方法上，它是采用欧拉变量去描述流体的动态，利用拉格朗日变量描述结构的运动边界，用光滑Delta近似函数通过分布节点力和插值速度来表示流场和结构物的交互作用。它整个流场计算都使用笛卡尔网格，而不是按照物体形状生成复杂的贴体网格，无需处理从物理平面到计算平面的坐标和网格转换问题，因而可以大大提高计算效率而且节省了网格生成所需的时间。尤其对于动态边界问题，它无需在每一时间步长上实时更新网格，详见图1。浸入边界法在模拟血液流动、湍流的直接数值模拟、多相流动等方面取得成功，模拟结果和实验数据非常吻合，是目前计算流体力学领域研究的热点。 2 浸入边界法的基本数学模型 在考虑包含曲线浸入弹性边界的粘性不可压流场的情况下，浸入边界系统的基本控制方程组可表示为： （1） （2） （3） （4） 其中，方程（1）和（2）表示一般不可压的N-S方程。式中，u是流场的速度，p是流场的压力，t是时间，p和u分别是流体的密度和动力学粘性系数。F是边界产生的力源项，反映弹性边界的存在。方程（3）表示力源项的分布函数。式中，f是浸入边界产生的单位力。单位力f的形式是根据具体研究问题所对应的力学模型来确定的。 表示浸入边界的位移，X(s,t)是Dirac Delta 函数，s用来描述流场和弹性体的交互作用。 3 浸没边界与LBM 根据Chapman-Enskog展开，上述N-S方程对应的格子方程为（以D2Q9为例） （4） （5） （6） 当然上述式子中的与的关系式还有其他多种表示方式。下面我们主要介绍怎样把边界模化为力项，主要工作是针对（3），（4）的离散，目前已有多种方法，先对（3）进行离散，可得如下式子 （7） 是Euler点，是Lagrange点，是边界上的弧长，而是二维的Dirac函数的逼近函数，定义如下 问题的关键在于如何确定边界上的力密度，目前发展了一系列方法 反弹法 反弹的办法是先对Lagrange点的粒子分布函数进行双线性插值，即先算出 可得到边界上新的分布函数 与是相反方向的分布函数，是边界的移动速度，我们可以算的流体对边界的作用 那么，边界对流体的作用力即为 直接力 直接力的着眼点是公式（6），在边界上有 其中可一看做没有力作用时速度，也可以通过插值的方式得到的值，这样也可得到一个力的表达形式