河口、海岸水动力模拟技术解析.ppt

河口、海岸水动力模拟技术 第一章 绪论 海岸：是海陆相互作用的重要地带，也是海、陆、气交互作用的重要空间，这种表现在： 岸线演变（自然和人为） 飓风（台风）带来的灾难性破坏； 海洋潮汐环境的变化。 河口：海岸常伴随有江河湖泊的出海口，通常称为河口。 海岸河口问题： 潮流问题 波浪问题 径流、异重流（密度流）、污染物（COD)扩散。 研究海岸河口问题的方法 物理模型（水力学比尺模型） 数学模型（数值模拟） 数值模拟：一门综合性的模拟技术，它采用数学模型来模拟某中物理现象，并通过计算机用数值计算法进行近似求解，籍以复演自然演变过程的总称。 水力学、泥沙数值模拟：以水力学和泥沙动力学为理论基础，并结合具体工程的一门新型实用科学。 水动力泥沙数值模拟：以微分方程为理论，并通过微分方程的离散，变成代数方程，最后采用计算机进行近似求解。 数值模拟的特点： （1）一般以线性理论为基础，但实际自然现象和描述这些现象的微分方程均为非线性的； （2）需要丰富的经验，现场资料和一定的技巧； （3）数值模拟不仅仅是一种近似计算，可以作为一种实验或研究及预测方法。 数值模拟的优点： （1）实验费用少； （2）速度快、周期短； （3）可以模拟多种因素相互作用的复杂物理过程。如可以模拟水（潮）流、风、柯氏力等多种因素共同作用下的多种泥沙及地形演变的复杂过程。 （4）可以完全控制流体的物理性质（如密度、容重、粘度、含沙量等） （5）模型建成后，长期保存、随时调用修改。 （6）无法模拟微分方程不能描述的物理现象。 数值模拟工作的基本步骤 （1）建立数学模型和编制源程序 建立或选择的微分方程； 根据模拟域边界条件选择合适的网格； 按一定的格式离散方程，得到代数方程和采用合适的数值方法求解代数方程； 编制源程序求解代数方程。 数值模拟分析（收敛性、稳定性、相容性、误差程度等） （2）调试源程序 （3）模型验证 调整模型中有关参数（糙率、紊动动量掺混系数等），使模型有良好的稳定性和收敛性，并与现场资料有良好的吻合； （4）正式方案试验 河口、海岸水动力模拟的发展方向 1、河口模型四维资料同化 四维同化在河口中的一类主要应用是为河口数值模式提供优化的初始状态,从而校正不合理的边界条件所带来的偏差,提高模拟精度. 变分同化技术还可应用于确定河口模式中的未知参数. 例如确定了潮汐河口的摩擦系数,通常这种系数是通过经验获得的. 另外,河口模式的外部强迫场(如风场、海-气界面的热通量等)都可以通过这种技术反演出来. 资料同化将是河口数值模型发展和结合的一个新技术切入点,也是带动河口动力数字模拟技术革新的一种重要方法. 2、数字河口动力模型 河口水动力特征及外界强迫作用因子构建数字河口动力模型(模型计算要素象温、盐、流、泥沙的特征均以数字矩阵形式记载),若再与数字河网模型嵌套联结,最终可以获得区域河网———河口的动力-沉积-地貌的机制解释,揭示河口物理规律和解决工程实际问题. 在数字河口模型的构架下,可以将人类已经拥有的河口科学理论、知识与具有较强物理概念的水动力学模型集成于一体,为河口数值模型计算提供良好的平台。 数字河口动力模型具有许多优势: 首先,数字河口模型是基于数字区域地形构建而成的,地形要素可自动生成,无需手工操作,大大提高了工作效率; 其次,数字模型不仅能输出传统模型的结果,而且能够十分方便地给出河口水文要素和水文状态变量的空间分布场,这些对近岸河口动力科学研究与河口、港口、航道工程都有着广阔的应用前景. 总而言之,数字河口模型研究的最终目的就是利用已有的河口基础科学理论和知识,在数字区域地形的基础之上将观测点的水文信息拓展、同化至区域平面上乃至区域三维立体上的信息,并形成数字成品。 参考文献： Koutitar 著“Mathematical Model in Coastal Engineering” 1)模型简单易懂 2）附有Basic程序，而且有验证的算例 3）介绍各种数值处理技术 曹祖德、王运洪”水动力泥沙数值模拟 第二章 水动力数值模拟的理论基础 三、二、一维方程的定解条件 初始条件 u,v,w,ζ|t=0=u0,v0,w0,ζ0 边界条件 开边界：计算域水体与外部水体相接处。 (u,v,w)Γ=(u(t),v(t),w(t)) ζΓ=ζ(t) 固边界：计算域与陆地或建筑物接壤处 无滑动：u,v,w=0 有滑动： 垂直边界的速度为0。 2.2数值计算 在计算水动力、泥沙数值模拟时，大都将基本方程组离散成代数方程组，最后求解代数方程组，此处介绍微分方程组的离散技术——有限差分法和线性代数方程组的数