第八章边界条件.doc

第八章 边界条件 任何数值模拟都可以认为仅仅是在物理区域或系统的一部分中进行的。区域的断层产生了人工边界，在这个断层中有我们处理的物理量。此外，还有暴露在流体中的自然边界。边界条件的数值处理需要特别注意。在实际的系统中处理不当模拟就会出现偏差。与此同时，稳定性和求解方案中的合成速度同样对数值模拟有消极的影响。 下边边界条件的类型是我们在欧拉方程和N-S方程中数值计算最常见的几种： ·固体壁面 ·外表面的远场和流体内部流出或流出的表面 ·对称面 ·平整切割和周期性边界。 ·平板间的边界 这些边界条件的处理问题在以后几节中会进行详细的介绍。对于文献中进一步涉及的边界条件，比如壁面上的热辐射或者是自由表面上的（热辐射），读者可以在3.4节中了解。 8.1 虚拟单元的概念 在我们讨论边界条件时，我们需要提到虚拟单元（也可以被称作虚拟点）这个概念。在规则的网格中这种方法非常的流行。然而，在不规则网格中，虚拟单元仍然有很多的优点。虚拟单元是在物理区域外部附加层上的一些网格点。这个可以由图8.1中的二维规则网格中看到。正如我们看到的，整个计算区域被两层虚拟单元包围着（由虚线标出），虚拟单元（点）通常不会像区域内的网格一样产生（除过多平板的网格）。尽管它仍然有几何形状，比如体积或者表面的矢量，但是它仅仅是虚拟的。 利用虚拟单元可以简化计算沿边界的通量，梯度，散度等等。这是由于在边界上可以将空间离散的模型进行扩展。正如图8.1中我们看到的，在物理区域内同样可以进行离散。因此，我们可以在所有的“物理”网格点中求解控制方程。这种方法可以使离散工作非常简单。此外，所有规则的网格点可以存在在一个单独的区域内，这在矢量计算中非常很有用。虚拟单元不但包含有守恒变量，同时也有几何量。很明显的是，虚拟单元层必须完全覆盖物理区域外。几何量通常由边界的控制体积来求得。在多网格平板中（3.1节），所有的流体变量和几何变量可以从相邻的平板求得。 图8.1中灰色的虚拟单元产生了一个问题，由于不清楚怎么设置它的值（如果没有相邻网格板），通过环绕形的离散模型，我们不需要知道它的值。但是，在计算梯度（粘性通量—见4.4节）或者在多网格中计算转移量时就很有比必要知道它的值。通常，如图8.1中用箭头表示的一样，相邻的“规则的”虚拟单元的平均值是很有用的。 8．2 固体壁面 8.2.1非粘性流动 在非粘性流动中，流体流过表面，由于没有摩擦力，速度矢量与壁面相切，与就是说在壁面的法线方向没有作用力。即： SKIPIF 1 0 在壁面上， （8.1） n表示壁面上的单位法相矢量。因此，相对应的速度V等于0。即，对流通量的矢量简化成只有压力项： SKIPIF 1 0 （8.2） 其中Pw为壁面压力。 规则的单元中心方案 在以单元中心的方案中，在单元的质心可以求出压力。但是，在等式（8.2）中，边界表面单元的Pw需要求解出。我们可以很容易的通过内部区域推导出壁面的压力。考虑到图8.2，我们可以简单的设Pw=P2。通过任一个两两相连点，我们可以取得较高的精确度。 SKIPIF 1 0 （8.3） 或者一个三点相连的推导公式是 SKIPIF 1 0 （8.4） 为了解释网格的延伸，到壁面的距离可以替代常数系数。【1】 上式的推导公式8.3和8.4没有涉及到网格和表面的几何形状。一个可供选择的方法—Rizzi发明的被称作为法向动量的关系。它是基于在非粘性流动中，壁面是流线型的。在式8.1中沿着流线型壁面的法线流动方向的导数为0。并且动量方程可以替换成如下的公式： SKIPIF 1 0 （8.5） 等式（8.5）包括密度，速度和壁面法线方向的压力。法向动量在【1】中给出了精确的结果。但是在复杂的几何表面，法向动量的数值求解存在着问题。推导中详细的描述和精确的比较可以在文献【1】中找到。 虚拟单元中守恒变量的值在壁面内部能够被线性的推导出： SKIPIF 1 0 （8.6） 公式（8.6）中的参数和图8.2中的相对应。如果虚拟单元在空间上离散后可以被利用，那么对流通量的计算就可以和公式（8.2）一致。 规则的单元顶点的方案 通过重叠的控制体积（见4.2.2），对单元顶点的离散可以直接由方程（8.1）在边界条件上应用。壁面的对流通量可