第十八章 计算水力学基础.pptxVIP

1; 18.1 水动力学基本方程与边界初始条件 18.2 微分方程的离散方法 18.3 有限差分法 18.4 对流—扩散—源汇方程的差分格式;水动力学基本方程和边界、初始条件一起构成了对一个水力学问题的完整的数学描述。水动力学基本方程一般描述水体的运动过程中基本动力学属性（如流速、水深、压强）的变化特征，目前有不同的控制方程可以描述水体的流动规律，如描述一般流体运动的最基本的N—S方程，在此基础上简化的水深平均平面二维浅水方程、一维圣维南方程、考虑多孔介质影响的地下水流动系列控制方程等。通过求解水体运动控制方程，一般可以获得水体不同位置的流速分布和空间变化。本章将从液体运动所满足的基本方程（N—S方程）出发，介绍计算流体力学中常见的基本方程。; （1）N—S方程的基本形式 实际的液体运动中，液体所占据的每个空间位置的不同时刻都有一定的流速、压强等，而N—S方程描述了这些物理量之间的量化关系。N—S方程一般将液体看作三维流动，各运动要素均是空间坐标（x,y,z）与时间（t）的函数，所以N—S方程研究的液体最普遍的运动形式，一般形式为 （18.1） （18.2）;（2） 边界条件与初始条件 N—S方程一般用于模拟复杂流动现象，包括有压流和无压流，以及这些水流和结构物的相互作用。 对于有压流动问题，需要给定计算区域各个边界上满足的条件。下面以经典平面二维圆柱绕流问题为例给出有压流动常见的边界条件，如图18.1 ;(a)入流边界：u=常数，或u=u(t),v=0 (b)出流边界：du/dx=0,dv/dx=0,dp/dx=0（假设流动已经 充分发展，出流边界处趋于稳态）。 (c)固壁边界：u=0,v=0 （液体在壁面处无滑移运动，考虑 了液体和壁面的粘性作用，包括计算区域上、下边界和 结构物的边界）。;对于带自由表面的流动问题，除了入流、出流、固壁边界 条件与有压流动一致外，还需要处理自由表面边界条件， 如图18.2所示。;(a)入流边界：给定入口随时间变化的流速和自由表面随时间变化的位移 (b)出流边界：一般设置在流动已经充分发展的区域，通常可以认为流动 已经达到自相平衡。可以假定出口处流动参数在出口方向的梯度为零，即 （18.3）式中，Q表示N-S方程组(18.1)、(18.2)中的速度分量u,v和压强P.(c)固壁边界:u=0，v=0,w=0(液体在壁面处无滑移运动，考虑了液体和壁面的 黏性作用，包括计算区域上、下边界和结构物的边界)。(d)自由表面：动力学边界条件和运动学边界条件。 对于自由表面的动力学边界条件，由应力平衡可得 （18.4）;18.1.2积分降维类方程N-S方程反映了水力学现象的三维流动特性和细节，但数值 求解难度较大。实际水利工程中大量的流动现象具有二维 (如水库、湖泊)或一维(如明渠流)特征，因此，在计算水力 学中可以针对水力学问题的流动特性対N-S方程沿坐标轴进 行积分处理，这样就可以将三维N-S方程降维到二维或一维 的积分方程。下面，将分别介绍沿垂向积分类(即水深方向)、 沿展向积分类(即河宽方向)、沿断面积分类方程组(即水深和 河宽两个方向)。;1.沿垂向积分——浅水方程 在河口、湖泊、海湾及水库等具有开阔水面的水域，水平 尺度远大于垂向尺度。流速在垂直方向的大小和变化远远 小于水平方向的大小和变化，因此垂直方向的流速可用流 速沿水深方向的平均值来表示，从而推导出水深积分的二 维浅水方程(shallow water equation)。;2.沿展向积分——二维水库水温模型沿水深积分的浅水方程适用于水深较浅的情况。当水库水深 较深时，水库水温会呈现分层现象，使得水温沿垂向具有一 定的温度梯度，这样沿水深积分的浅水方程就不再适用于深 水水库。对于水库水温分层问题，可以沿河流展向(宽度方向) 对N-S方程进行积分，得到立面二维积分方程组，其一般形式 为;3.沿截面积分——一维圣维南方程 现在考虑河道中的水流，占主导地位的流动方向是沿河道 方向(流向)，而沿河道展向(横截面方向)的流动较弱。在此 情况下，可以沿河流展向(横截面方向)对平面二维浅水方程 进行积分，河流展向的流动信息以展向积分的形