河口、海岸水动力模拟技术--12讲义讲课讲稿.ppt

河口、海岸水动力模拟技术-2014-12讲义;海岸河口问题： 潮流问题 波浪问题 径流、异重流（密度流）、污染物（COD)扩散。 研究海岸河口问题的方法 物理模型（水力学比尺模型） 数学模型（数值模拟） ;沿岸过程;数值模拟：一门综合性的模拟技术，它采用数学模型来模拟某中物理现象，并通过计算机用数值计算法进行近似求解，籍以复演自然演变过程的总称。 水力学、泥沙数值模拟：以水力学和泥沙动力学为理论基础，并结合具体工程的一门新型实用科学。 ;水动力泥沙数值模拟：以微分方程为理论，并通过微分方程的离散，变成代数方程，最后采用计算机进行近似求解。 数值模拟的特点： （1）一般以线性理论为基础，但实际自然现象和描述这些现象的微分方程均为非线性的； （2）需要丰富的经验，现场资料和一定的技巧； （3）数值模拟不仅仅是一种近似计算，可以作为一种实验或研究及预测方法。 ;数值模拟的优点： （1）实验费用少； （2）速度快、周期短； （3）可以模拟多种因素相互作用的复杂物理过程。如可以模拟水（潮）流、风、柯氏力等多种因素共同作用下的多种泥沙及地形演变的复杂过程。 （4）可以完全控制流体的物理性质（如密度、容重、粘度、含沙量等） （5）模型建成后，长期保存、随时调用修改。 （6）无法模拟微分方程不能描述的物理现象。;数值模拟工作的基本步骤 （1）建立数学模型和编制源程序 建立或选择的微分方程； 根据模拟域边界条件选择合适的网格； 按一定的格式离散方程，得到代数方程和采用合适的数值方法求解代数方程； 编制源程序求解代数方程。 数值模拟分析（收敛性、稳定性、相容性、误差程度等）;（2）调试源程序 （3）模型验证 调整模型中有关参数（糙率、紊动动量掺混系数等），使模型有良好的稳定性和收敛性，并与现场资料有良好的吻合； （4）正式方案试验 ;河口、海岸水动力模拟的发展方向 1、资料同化将是河口数值模型发展和结合的一个新技术切入点,也是带动河口动力数字模拟技术革新的一种重要方法. 2、数字河口动力模型 ;;;2、数字河口动力模型 河口水动力特征及外界强迫作用因子构建数字河口动力模型(模型计算要素象温、盐、流、泥沙的特征均以数字矩阵形式记载),若再与数字河网模型嵌套联结,最终可以获得区域河网———河口的动力-沉积-地貌的机制解释,揭示河口物理规律和解决工程实际问题. 在数字河口模型的构架下,可以将人类已经拥有的河口科学理论、知识与具有较强物理概念的水动力学模型集成于一体,为河口数值模型计算提供良好的平台。;;;第二章 水动力数值模拟的理论基础;U,V,W为x,y,z 方向上的流速分量。 η（ｘ，ｙ，ｔ）为距平均海平面的自由表面水位。 ｈ（ｘ，ｙ）为平均海平面距底部边界的水深。 ＡＨ为水平扩散系数。 ＡＺ为垂直涡动系数。　;初始条件;Saint Venant 方程;三、二、一维方程的定解条件;2.2数值计算;2.2.1有限差分法;;（3）一阶中心差分;;;;（3）六点格式（Crank-Nicolson），双层六点隐式格式 在x点和n+△n/2时层，对t和x均采用中心差分;（4）Richardson格式，三层显式格式 在x点和n时层，对t和x均采用中心差分;（5）加权六点格式，隐式格式 在x点和n+θ△n时层，0≤θ≤1，对t和x均采用中心差分;2.2.2线性方程组的数值解;1、解线性方程组的两种方法: 直接法:通过有限步算术运算直接求出方程组的精确解,最常用的是消元结合代入的方法. 实际上除非是采用无穷位精度计算,一般都得不到精确解. 直接方法适用于解低阶稠密矩阵方程组.;迭代法 类似于方程求根的迭代法,用一个迭代过程逐步逼近方程组的解. 迭代有可能不收敛,或虽然收敛,但收敛速度慢. 迭代法适用于求解高阶稀疏矩阵方程组. 稀疏矩阵:矩阵非零元素较少,且在固定的位置上. 稀疏矩阵一般是人为构造的,例如36页三转角插值时方程组(8.12),(8.15)的系数矩阵.;Gauss消去法(第一次消元);Gauss消去法(第一次消元);;Gauss消去法(第k次消元);Gauss消去法(第k次消元);Gauss消去法;Gauss消去法(算法);追赶法求解三对角方程组;比较两边对应元素可以得到;因此有;在得到系数矩阵的分解后,原方程组转化为 LUx=f.;;迭代法;Jacobi迭代法与Gauss-Seidel迭代法;Jacobi迭代法;Jacobi迭代法;Jacobi迭代法;Jacobi迭代法;Jacobi迭代法;Gauss-Seidel迭代法;Gauss-Seidel迭代法;Gauss-Seidel迭代法(算例);Gauss-Seidel迭代法(算例);超