曲线坐标下波浪抛物型缓坡方程及近岸流数值模型研究

建立了曲线坐标系下的波浪抛物型缓坡方程及近岸流数值模型,对不规则岸线地形下的波浪传播变形及近岸流运动进行了数值模拟研究.首先基于曲线化抛物型缓坡方程建立了正交曲线坐标系下的近岸波浪传播数值模型;然后通过正交曲线变换建立了正交曲线坐标系下的近岸流数值模型;进而基于波浪辐射应力的概念,耦合上述数值模型,对Gourlay模型实验和Hamm模型实验中不规则岸线地形条件下的波浪传播变形及近岸流运动进行了数值模拟研究,并结合实测资料对数值模拟结果进行了验证分析;最后基于上述耦合数值模型对美国San Francisco海湾的Ocean Beach海岸波浪斜向入射所产生的波浪场及波浪破碎所形成的近岸流进行了数值模拟研究.对比分析表明,本文数值结果与实验结果或其它学者所得数值结果基本一致,从而验证了本文数值模型的有效性.     

基于曲线坐标系下缓坡方程模拟近岸波浪及波生流

为了更好地模拟曲线边界等不规则岸线附近,波浪传播及其破碎后形成的波生流问题,该文建立了正交曲线坐标下近岸波浪和波生流耦合数值模型;即采用可综合考虑近岸区波浪反射、折射和绕射等效应的双曲型缓坡方程计算波浪场,采用近岸深度平均的二维浅水方程计算波生流场,应用辐射应力的概念,使波浪模型为流场模型提供驱动力.波浪和波流模型基于相同的正交曲线网格,应用高效的ADI有限差分格式进行数值求解.将曲线坐标系下的耦合模型应用到半圆形水槽和凹型水槽两个物理模型中,数值模拟的结果分别和解析解、实验值及已发表的数值结果做了比较.模型实验验证表明,该文建立的数值模型比较准确、有效,可以应用到防波堤和河口等具有复杂边界和复杂地形的实际海域.     

抛物形缓坡方程的数值研究

对2种典型的抛物形缓坡方程进行比较系统的数值研究。通过对4种典型地形上的波浪变形的数值模拟，详细讨论了网格节点数对数值解精度的影响、模型对初始入射角的敏感程度、非线性项对数值结果的影响等。研究结果可为实际应用抛物形缓坡方程研究大区域复杂地形上的波浪传播提供一定的理论指导。     

应用时间关联型缓坡方程数值模拟波流相互作用

从含水流项的时间关联型缓坡方程出发,采用有限差分法,建立缓变水深和水流水域波浪传播的数值模型.模型对空间导数项采用三点差分格式离散,对时间导数项采用Euler预测-校正格式离散.基于开边界条件与不同反射特性的固壁边界条件相统一的边界条件表达式,对边界条件进行处理.该模型数值解与椭圆型缓坡方程有限元模型的数值解较为吻合,表明所得数值结果合理可信.     

波浪Boussinesq方程与缓坡方程的比较

分析、比较了近岸浅水波浪传播变形的Boussinesq方程和缓坡方程的形式与特点,建立了包含底摩擦能耗效应、波浪破碎效应和子网格效应的Boussinesq方程波浪数学模型,并介绍了处理动边界问题的窄缝法以及处理消波边界的海绵层技术．采用经验非线性色散关系,结合含非线性项的缓坡方程,得到考虑非线性作用影响的缓坡方程模型．用物理试验结果对两种模型进行验证,并用相关性分析方法对两模型的计算精度进行了分析与说明．     

黄茅海潮流波浪计算分析

采用海岸工程水动力学数学模型对黄茅海潮流场及波浪场进行了计算.潮流计算中利用水平方向上的曲线正交坐标与垂直方向上的Sigma坐标相结合以及三维模型二维化的方法,波浪计算采用双曲型缓坡方程,并引入改进型的非线性弥散关系,从而能反映波浪的非线性效应.计算结果与实测结果比较表明：模型能较好地复演黄茅海潮流及波浪场,可以用于工程实际中沿岸及河口湾潮流波浪场问题的计算.     

基于非线性弥散关系的缓坡方程波浪传播变形模拟研究

【目的】波浪由外海传播到近海时,由于受到地形、建筑物等影响,波浪非线性增强,线性弥散关系不能够很好的描述波浪弱非线性效应。为了对比研究非线性弥散关系的缓坡方程在波浪传播变形的作用。【方法】采用改进型缓坡方程数值模式,并结合Li提出的非线性弥散关系,对Berkhoff椭圆经典地形进行波浪传播变形模拟研究,探讨线性和非线性弥散关系的数值模拟计算结果与实验值的关系,并对两种计算结果进行了比较分析。【结果】非线性弥散关系的计算结果与实验值的误差较线性弥散关系的结果小,非线性模型要优于线性模型。【结论】非线性模型更适合近海海域弱非线性波浪传播变形的研究。     

基于分层Boussinesq类模型的孤立波爬坡模拟

 波浪爬坡是自然界中普遍存在的现象,是近岸防波堤等结构物设计需要考虑的因素。本文采用分层Boussinesq类方程对近海岸波浪爬坡开展数值模拟研究,采用有限差分方法数值离散二层Boussinesq类方程,建立数值波浪模型。根据Synolakis 1986年的物理实验,数值模拟孤立波在梯形水道上的爬坡,验证了分层Boussinesq方程波浪模型在该研究中具有很好的适用性和较高的精度。模型较好地预报了波浪爬高的变化趋势。预报的波浪爬坡高度大于实验测量值,主要原因是实际的底摩擦大于模型中的底摩擦。     

缓坡方程计算分辨率选取的数值研究

自适应曲线网格技术是近年来国际上近岸海洋计算中采用较多的先进技术,而如何减少计算量一直是缓坡方程的一个研究热点,如将自适应曲线网格技术引入缓坡方程,则既提高局部海域分辨率又减少计算量,为了在缓坡方程计算模型中很好地引入自适应曲线网格技术作理论准备,绎缓坡方程计算分辨率的迭取问题作较深入的探讨,首先建立一个Ebersol型缓坡方程数值模型,然后设计了一组网格分辨率对计算结果影响的数值实验,讨论分辨率的取值范围,并且对缓坡方程大网格计算的可行性作了初步的理论分析,最后得出结论,若仅考虑地形的影响Ebersol型缓坡方程计算分辨率的选取主要决定于地形,即网络大小可突破波浪波长限制。     

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场