基于无结构、σ坐标网格的半隐式三维水沙数学模型Ⅰ-算法.pdfVIP

第 30卷 第 6期 水 力 发 电 学 报 V01．3O No．6 2011年 12月 JOURNALOFHYDROELECTRIC ENGINEERING Dec．，2011 基于无结构、 坐标网格的半隐式三维水 沙数学模型 I一算法 胡德超 ，张红武 ，哈岸英。，赵维阳 (1．清华大学水沙科学与水利水 电工程国家重点实验室 ，北京 100084；2．长江科学院河流研究所 ， 武汉 430010；3．西安理工大学 ，西安 710048；4．长江航道规划设计研究院，武汉 430011) 摘 要 ：针对天然河流水域边界不规则的特点，提出了平面无结构、垂 向ro坐标网格上三维 自由水面水沙数学模型的基 本框架 。模型基于压力分裂模式和 0半隐算法建立 ，求解分为两步，第一步求解一个关于 自由水面的二维 Poisson方程 ， 获得水位和临时流场 ；第二步求解一个关于动水压强的三维 Poisson方程 ，获得动水压强场并用它来校正临时流速场 。 在构造动水压强三维Poisson方程时借鉴了BOM模型方法，而在 自由水面处理上摒弃了其 内外模分裂方法 ，因而模型不 受模式分裂方法 中一致性 问题的影响。在三维水流模型的基础上建立 了悬移质泥沙模型，采用欧拉 一拉各朗 日方法 (ELM)求解控制方程 中的对流项，使模型计算的时间步长不受 CFL稳定条件限制 。 关键词 ：泥沙动力学；三维；数学模型 ；水流 ；泥沙 ；无结构网格 ；or坐标 ；自由水面 中图分类号 ：TV143 文献标识码 ：A 3-D hydrodynam icmodelforfree-surfaceflowsandsedimenttransporton unstructured q-coordinategrid．PartI．Numericalscheme HU Dechao一，ZHANG Hongwu ，HA Anying’ ZHAO W eiyang ， (1．StateKeyLaboratoryofHydroscienceandEngineering，TsinghuaUniversity，Beijing 100084； 2．ChangjiangRiverScientificResearchInstitute，Wuhan 430010；3．XianUniversityoftechnology，Xian 710048； 4．ChangfiangDesignandResearchInstituteofWaterway，Wuhan 430011) Abstract：A formulationof3一D hydrodynamicmodelisproposedforsimulationofcomplexfree—surfaceflows and sedimenttransportin naturalriversOilunstructured．盯一coordinate grids．Two computation steps are adoptedtosolve the equationswith apressure—splittingand0 semi-implicitmethod．First，freesurfaceand provisionalvelocityfiledareobtainedbysolvinga2-D Poissonequation；second，non—hydrostaticpressuresare computedbysolvinga3-D Poissonequationandtheprovisionalvelocityfilediscorrected．Inconsturctingthis 3-D Poissonequation，themethod