一维溃坝涌波的特征线-激波装配法

基于特征线上的相容方程和激波点处的间断条件,构造了一维溃坝涌波的一种拟合计算的显格式,并进行了一维矩形断面水坝瞬间全溃洪水波演进的数值模拟.数值模拟结果与理论解的比较表明,这种方法精度高,很好地模拟了溃坝涌波和展示了溃坝波的运动特性,是研究溃坝涌波一种高精度的数值模拟方法.     

基于特征线理论的算子分裂算法综述

文中论述了基于特征线理论的算子分裂算法的研究进展,从偏微分方程特征线理论及算子分裂算法出发,构建了其理论框架,对于CBS算法的基本方程,文中给出了推导过程,并编制了相应的三维计算程序,通过数值算例进一步验证了这一算法的有效性和准确性.     

一维溃坝洪水波的数值模拟-时空守恒法

应用改进的时空守恒法对一维溃坝洪水波进行了数值模拟,从计算结果与分析解比较可知,该法能够很好地模拟溃坝洪水波,甚至对于下游为干河的情况,依然能得到满意的结果.而且,此方法精度高,有较强的溃坝涌波捕捉能力,可展示溃坝波的运动特性,为溃坝引起的灾害预测和分析及防灾减灾提供一种有效的数值模拟.     

二维溃坝洪水波演进的数值模拟及动态显示

采用基于MacCormack预测—校正技术的隐式数值格式求解控制水流运动的二维浅水方程，建立了模拟大坝瞬间全溃或局部溃倒所致的洪水演进过程的数学模型；应用该模型对矩形明渠中缓、急流过渡的水面曲线进行了数值计算，并与理论解进行了比较；最后用该模型预测了矩形河道中大坝瞬间局部溃倒的洪水演进过程，并对模拟结果进行了定性分析。比较与分析结果表明，该模型对模拟溃坝洪水波很有效。     

二维溃坝波传播的数值模拟

采用附加人工粘性的显式Maccormack格式求解二维浅水方程,并对人工粘性加上开关函数,建立了模拟大坝瞬间局部溃倒所致的洪水演进过程的数学模型。用该模型预测了矩形河道中堤坝瞬间局部溃倒后8 min的洪水演进过程,得出了水面线和流场分布,并对溃坝洪峰到达时间和最大淹没水位的沿程分布作了分析。     

二维溃坝问题的高分辨率数值模拟

对二维浅水方程空间方向的差分分裂,引入ＭＵＳＣＬ型的两参数限量函数,采用组合型的二阶ＴＶＤ格式,建立了浅水方程的高分辨率计算模型．以Ｆｒａｃｃｒｏｌｌｏ 等的部分溃坝实验模型为物理模型,模拟了该模型长历时的溃坝洪水演进过程,给出了典型位置的水位和速度计算并与其实验结果的比较．研究表明,在间断波附近,没有产生数值振荡,数值结果与实验数据相比除溃坝初期在口门附近的水位存在一定的差异外,其余位置和时刻吻合较好,特别是计算速度与实验测定值吻合更好．应用二维浅水方程和本文的高分辨率数值方法能有效地模拟溃坝波的演进过程．     

溃坝洪水演进的数值模拟

开发了溃坝洪水演进数值模型,对模型控制方程采用有限体积法进行空间离散,采用具有总变差减小特性的MacCormack预测-校正格式进行时间离散.采用非规则局部不连续网格划分计算区域,增强了模型对复杂地形的适应能力;具有总变差减小特性的数值格式使模型具备了强大的间断波捕捉能力,并消除了模型稳定性对人工黏性的依赖;基于单元属性与基于界面属性的动边界处理技术使模型能合理地模拟非恒定性强、干湿变动剧烈的水陆边界.实际算例表明：模型在强非恒定流状态下具有稳定性及守恒性,可作为溃坝洪水预报及其灾害防范的有力工具.     

基于二维浅水方程的溃坝水流数值模拟

针对溃坝水流的特点, 提出了一种基于二维浅水方程的数学模型, 对溃坝水流进行数值模拟研究. 模型采用二阶leap-frog格式和一阶迎风格式相结合的有限差分格式进行数值离散, 计算网格采用时间和空间交错网格. 在模型建立完善的基础上, 对模型进行一维和二维的数值验证, 主要算例包括:下游有水的溃坝洪水理论解, 带有三角形障碍物的一维溃坝模型试验以及二维溃坝数值模型试验. 通过与理论解和实测值比较, 结果表明模型具有良好的稳定性以及模拟间断水流和处理动边界问题的能力.     

不可压缩非等温非牛顿流广义特征线基分裂算法

发展了用于不可压缩非等温非牛顿流模拟的广义特征线基分裂算法,推广了控制方程沿粒子特征线的离散格式,经典的特征线伽辽金法和时域的Crank-Nicolson离散格式可分别被归纳为离散格式的两个特例.为消除流体的不可压缩性导致的虚假数值振荡,采用压力稳定型分步算法.利用粘性系数关于等效应变速率的指数模型来模型化流体的非牛顿特性.平面Poisseuille流动和4∶1平面收缩流动的数值模拟结果突出地显示了该方法处理不可压缩非等温非牛顿流动问题的有效性及良好性能.     

WC-CMPS法对溃坝溢流问题的数值模拟

无网格粒子法(拉格朗日法)是计算流体力学领域最新一代的计算方法,由于不需要生成网格,无网格粒子法能够模拟任意形状的界面及其变形和碎裂,尤其适用于计算如溃坝溢流及波浪破碎等具有复杂自由表面运动的流体力学问题。其中,光滑粒子流体动力学方法(SPH)和移动粒子半隐式法(MPS)是无网格计算方法的典型代表。本文采用基于MPS法的修正弱可压缩移动粒子半隐式(Weakly Compressible Corrected Moving Particle Semi-implicit Method)法构建了用于求解自由表面流问题的数值计算模型,并用于模拟干床面与湿床面上的溃坝溢流问题。计算结果与实验结果的对比分析表明,采用基于WC-CMPS无网格法建立的数值计算模型对干、湿床面上溃坝溢流的溢流形态、水柱高度以及波前速度等均给出了良好的计算结果。     

狭窄库区和开阔溃坝区溃坝洪水物模试验研究

提出的沿水发段加（减）糙方法，使狭窄库区模型在水深大幅度变动情况下，为满足水流对模型糙率的要求而同时实现重力相似和阻力阻似，提供了有效的解决途径，阐述了溃坝洪水试验模型设计布置，列举了实例，所采用的自动控制技术和量测试器设备，可模拟各种溃决过程并能获得所需的各项试验参数，指出了开阔溃区溃坝洪水演变的基本特性。     

劈裂灌浆改善土坝坝体应力状态研究

山东黄前水库粘土心墙土坝采用劈裂灌浆技术加固，本文根据Ｂｉｏｔ固结理论，利用有限元方法对坝体的变形和应力状态及坝轴线附近劈裂缝的变形进行了计算分析，并与现场实测资料比较，计算中对坝体土料的应力应变关系采用椭圆－抛物双屈服而弹塑性模型，以反映坝体劈裂灌浆加荷的复杂性。     

溃坝问题数学模型研究

本文研究溃坝问题的数学模型,探讨了有限差分法和谱方法用于计算,提出了改进意见,数值模拟表明结果是令人满意和符合物理特性的。     

基于扩展有限元法的混凝土重力坝多裂纹数值模拟

针对含有多个初始裂纹的混凝土重力坝的裂纹扩展问题,在ABAQUS有限元分析软件上,应用扩展有限元法和黏聚裂纹模型来进行数值模拟。分析了多裂纹的扩展路径以及主导裂纹的荷载-裂缝口张开位移。结果表明,混凝土重力坝多裂纹的扩展路径与文献结果几乎相同,主导裂纹的荷载-裂缝口张开位移曲线与文献结果有差异,但差异很小。可见采用扩展有限元法能够精确地模拟混凝土重力坝多裂纹从加载到失稳破坏的全过程,无需重新划分网格,减少了计算量,为解决工程实际问题提供了一种有效的应用途径。     

基于LES方法仿生翼型气动特性数值模拟

采用大涡模拟方法,研究在翼型不同位置添加脊状结构对翼型流场及气动性能的影响,讨论了添加脊状结构后翼型流场的流动特性和涡结构特性。研究发现：1)在α=6°攻角条件下,无论riblet-Q翼型模型或riblet-H翼型模型均可改善边界层分离情况,但riblet- H翼型模型表现出更好的控制效果。2）后段布置脊状结构能够有效推迟翼型边界层分离点,抑制边界层大涡形成,控制分离涡的发展和脱落。3）riblet-H翼型模型使翼型的升力系数增大,同时也使其阻力系数降低,升阻比较原翼型有了较大提高。     

复杂地形上坝堤溃决洪水演进的非结构有限体积模型

针对复杂地形上坝堤溃决洪水演进的二维数值模拟,综合考虑底坡项、摩阻项和干湿界面对水流的影响,基于浅水二维水动力学理论建立了非结构有限体积模型.围绕传统底坡项处理技术应用于复杂地形时的不足,通过计算底坡项对物质通量和动量通量的影响,并进行相应的分解,提出了集成底坡项的近似Riemann求解器.分析并采用改进的隐式格式解决了摩阻项可能引起的刚性问题.通过求解干湿单元界面处"镜像湿单元"Riemann问题,结合变量重构技术,有效处理了动边界问题.算例研究结果表明,模型稳定健全,和谐性好,具有较好的推广应用价值.     

基于LASSO方法的结构突变理论研究综述

结构突变是统计学、经济学、信号处理和生物信息学等学科领域中的研究热点之一.Z.Harchaoui 等提出了基于LASSO的结构突变点检测方法,是近几年结构突变问题的最新研究方法.为了在国内推行该方法,系统介绍了国外基于LASSO方法的几种变点模型中的变点检测问题,其核心是把变点检测问题转化成模型选择问题来解决,并阐述了相应的算法.最后探讨该方法在不同学科领域的应用和前景展望.     

铀尾矿库周围地下水运动数值模拟

铀尾矿库中的放射性核素及重金属等有害物质随着渗水进入含水层后,将随地下水的运动而进行迁移。为了研究污染物在地下水中的迁移规律,采用二维稳定流数学模型对铀尾矿库周围地下水流场进行概化,并用Galerkin有限元法对模型进行了求解。将数值模拟的等水位线同实测的等水位线进行对比可以发现,两者拟合较好,将流场模拟的结果作为水质模型的输入参数是比较可靠的。     

三泉水库坝址工程地质特征评价

介绍了三泉水库的工程概况和地质情况,对坝址的工程地质特征及主要存在的工程地质问题进行了合理的工程地质评价。