相变传热与流体流动数值分析（第11-14讲）.pptVIP

// Zou and He pressure boundary on north side. for( i=0; ini; i++) { fi = ftemp[nj-1][i]; uy0 = -1. + ( fi[0] + fi[1] + fi[3] + 2.*( fi[2] + fi[5] + fi[6])) / rho0; ru = rho0*uy0; fi[4] = fi[2] - (2./3.)*ru; fi[7] = fi[5] - (1./6.)*ru + (1./2.)*( fi[1]-fi[3]); fi[8] = fi[6] - (1./6.)*ru + (1./2.)*( fi[3]-fi[1]); } 压力（密度）边界 14.2.3 非平衡态反弹格式 对于充分发展流动而言，如果边界条件不能使流体在整个流动区域内满足质量守恒要求，则数值结果的精度就会降低，而且数值计算的收敛速度以及稳定性都会受到影响； 如果在处理充分发展边界条件时，采取一定的方法，使得流体在整个流动区域内的总体质量守恒要求得到满足，那么数值结果的精度以及计算过程的收敛状况就会有较大程度的提高； 何雅玲等基于此思路，提出了格子Bolzmann方法处理充分发展流动的质量修正格式。 所谓质量修正，就是通过计算进口和出口的质量流率，确定一个质量修正系数，以此系数修正出口速度，并用修正后的速度更新边界节点上的未知粒子分布函数。 14.2.4 质量修正格式 具体做法如下： 假设出口截面上各点速度x方向分量的相对变化率为常数，即 然后，考虑整体质量守恒，也即进出口的质量流率相等，得质量修正系数为： 利用质量修正系数即可得到出口速度分布，利用该速度计算出口边界的粒子分布函数。 14.2.4 质量修正格式 启发式格式和动力学格式的局限性 依赖使用的格子Boltzmann模型， 通用性差； 需要对边界处的一些物理性质（如速度、密度、压力等）作假设，因而很难推广到更一般的边界处理上。 对一些非常规边界（包含梯度信息的边界）存在误差 借鉴传统计算流体力学方法中的边界处理方法来构造格子Bolzman方法边界处理格式； 外推格式 非平衡外推格式 14.3 外推格式 1996年Chen等人提出； 假定在物理边界外有一层虚拟边界，物理边界和虚拟边界格点上的分布函数值都按照标准的演化步骤进行演化，但虚拟边界格点上指向流场的分布函数值使用外推方法确定。 14.3.1 外推格式 流体层 物理边界 虚拟边界 在碰撞过程中在物理边界上使用给定的边界条件（速度或压力）来计算平衡态分布函数。 二阶精度 优点 普适性好，可以根据这种方法容易地设计出一般边界条件的格式； 计算量也比较小，容易实现。在确定未知分布函数时，不需要额外的假设条件，这是外推格式的技巧所在； 具有二阶精度； 缺点 数值稳定性较差（在外推格式中，因为有负的系数出现，所以数值误差有可能随时间增加，从而导致数值不稳定的发生） 14.3.1 外推格式 受外推方法和非平衡态反弹方法的启发，吸收两种方法的优点，2002年Guo等人提出了一种新的边界处理方法，即非平衡态外推方法。 其基本思想是，将边界节点上的分布函数分解为平衡态和非平衡态两部分，并根据具体的边界条件定义新的平衡态分布来近似平衡态部分，而非平衡态部分则使用非平衡态外推（一阶精度）确定，所得的分布函数的整体逼近精度为二阶精度。 14.3.2 非平衡外推格式 流场内部节点： G 边界节点： O 14.3.2 非平衡外推格式 流体内部 物理边界 流场外部 对于非平衡态部分，用邻近的流体格子节点G出的非平衡态部分近似得到，即 14.3.2 非平衡外推格式 对于速度边界：u已知，密度未知 空间精度是二阶的，时间精度也是二阶的，这与格子Boltzmann方法的整体精度是一致的。 由于低阶外推格式的数值稳定性一般要优于高阶外推格式的数值稳定性，所以非平衡态外推格式的数值稳定性要好于根据线性外推得到的分布函数外推格式的数值稳定性。 同时，与外推格式一样，非平衡态外推格式也具有适用范围广、计算简单，容易实现的优点。 14.3.2 非平衡外推格式 当所要处理的物理区域并不具有规则的几何形状时 适体网格 非结构化网格 在结构化的直角正交网格，并在适当的位置采用阶梯逼近或者差值处理，以保证满足物理边界上的条件。 阶梯逼近可与诸如反弹格式等前面提到的多种边界处理格式配合使用，实施简单，但是计算精度较低，常用的处理复杂边界处理格式为Filippova与Hanel格式、Bouzidi格式、Lallemand与