搜索
目录目录..................................................................................................................................................1第8章有限体积法.........................................................................................................................2§8.1有限体积法的基本概念...................................................................................................28.1.1几个术语..................................................................................................................28.1.2控制体积的选择.....................................................................................................38.1.3结构与非结构网格.................................................................................................5§8.2有限体积法构造..............................................................................................................78.9-1方程离散及基本格式............................................................................................78.9-2物理特性要求.......................................................................................................118.9-3迎风型通量格式..................................................................................................148.9-3TVD格式..............................................................................................................18§8.3非结构网格上的有限体积法........................................................................................238.3.1基本方程...............................................................................................................238.3.2离散基本思路........................................................................................................248.3.3数值通量近似.......................................................................................................25第9章在水力学问题中的应用...................................................................................................29§9.1渗流问题中的应用.........................................................................................................299.1.1饱和—非饱和地下水运动基本控制方程............................................................299.1.2方程的离散............................................................................................................319.1.3算例【陈扬硕士论文】......................................................................................33§9.2二维明渠非恒定流计算.................................................................................................389.2.1水流基本方程........................................................................................................389.2.2控制方程的离散....................................................................................................399.2.3计算实例...............................................................................................................53§9.3三维紊动分层流计算.....................................................................................................649.3.1紊动分层流基本方程............................................................................................649.3.2紊流模型及控制方程离散...................................................................................659.3.3压力校正法............................................................................................................669.3.3边界条件................................................................................................................699.3.4盐度引起的负浮力流动........................................................................................711第8章有限体积法有限差分方法是从描述各种物理现象的基本微分方程出发构造离散方程的，前文已经对其作了翔实、周密的论述。该部分将从基础算法入手分析介绍在计算流体力学界广为应用的有限体积法。基于有限体积法的实用算法在计算流体力学、计算传热学等领域得到了飞速发展[1-3]。在水力学诸多问题，如水流物质输运模拟，水工水力学模拟以及溃坝洪水波演进等水流模拟中也得到了广泛应用。§8.1有限体积法的基本概念有限体积法，又称为有限容积法，它正是从物理量守恒这一基本要求出发提出的。这也是其受计算流体力学界广为称道和喜欢之处。其以守恒型的方程为出发点，通过对流体运动的有限子区域的积分离散来构造离散方程。有限体积法有两种导出方式，一是控制容积积分法，另一个是控制容积平衡法。不管采用哪种方式导出的离散化方程，都描述了有限各控制容积物理量的守恒性，所以有限体积法是守恒定律的一种昀自然的表现形式。该方法适用于任意类型的单元网格，便于应用来模拟具有复杂边界形状区域的流体运动；只要单元边上相邻单元估计的通量是一致的，就能保证方法的守恒性；有限体积法各项近似都含有明确的物理意义；同时，它可以吸收有限元分片近似的思想以及有限差分方法的思想来发展高精度算法。由于物理概念清晰，容易编程；有限体积法成为了工程界昀流行的数值计算手段。8.1.1几个术语在进行数值计算时，要把计算区域划分成一系列互不重叠的离散小区域，然后在该小区域上离散控制方程求解待求物理量。在有限差分法中只涉及到网格节点的概念，而有限体积法因为物理解释需要，形成了以下几个常用几何要素的相关名词。2PEWSN单元2N4N3N1Newsn图8-1单元中心有限体积法中典型的矩形网格如图8-1所示，为单元中心格式（cell-centeredscheme）有限体积法中典型的矩形网格关系。其包含以下几个要素：控制体积（controlvolume）：图中阴影部分，方程积分离散时的小体积单元（二维为面积单元）。单元（cell）：控制体积的中心，常用形心来表示，为待求物理量的几何位置。图中用空心园来表示，如点P、W、E、N、S等。常用单元来代表整个控制体积。以下如若不作特殊说明，则用“单元”来代表控制体积。网格线（gridlines）：用来分割计算区域内各控制体积的交错曲线簇，如图中折线N1N2、N2N3、N3N4、N4N1等网格节点（nodes）：网格线之间的交点，图中用黑圆点来表示，如N1、N2、N3、N4等单元界面（cellfaces）：相邻两个控制体积间的公共面(二维则可以认为是公共边，这样看起来就和网格线一致了，但是要注意这不是同一个概念)，图中用小写字母e、n、w、s表示。通常定义e、n、w、s几何位置位于交界面的形心点，二维则认为在公共边的中心点。8.1.2控制体积的选择当你开始用有限体积法模拟流体流动时，而且划分好网格后，你必须选定控制体积的形成方式。目前，常用的有两种方法：单元中心方式（cell-centered）和顶点中心方式（vertex-centered）。另外一些学者还发展了由两种方式综合形成的混合方式。根据问题的特点和要求，不同的变量可以采用不同的控制体积，因此又产生了交错网格和同位网格的称谓，这里不再深入介绍，读者根据需要可以参考相关文献[1-3]。前面我们我们在讨论有限体积法的术语时，已经看到了单元中心控制体积的形成方式。这里再作一个说明并对两者作简单的比较。图8-2为两种控制体积选择方式示意图。阴影部3分表示单元的控制体积。空心圆点表示单元，实心圆点表示网格节点，在顶点中心方式中，单元和网格节点重合。PPEEWWSN控制体积控制体积2N4N3N1Newsn2N4N3N1NPEWe