完全非线性孤立波的直墙反射

报道了应用边界积分方法模拟完全非线性孤立波的传播与直墙反射,给出了波形演变过程。结果表明,本模型对计算孤立波的传播与直墙反射是有效的。三阶Boussinesq方程的孤立波解比低阶方程的孤立波解更接近完全非线性的数值解.当来波波高增大时,孤立波直墙反射的相位滞后变小。若考虑大波高孤立波的直墙反射或波——波相互作用,一阶理论预报的相位滞后往往低估实际情况。     

内孤立波波-波相互作用的数值模拟

基于弱二维的KP方程，并结合南中国海东沙群岛附近内孤立波的观测资料，模拟了内孤立波的波-波相互倌用0数值结果较好的反应了内孤立波的二维特征，同时体现两个内孤立波波-波相互作用的非线性特征，即两波相交处相速随振幅的增大而变大。相比于一维的KdV方程，KP在内孤立波的仿真反演方面具有更大的优势。     

内孤立波作用下半潜平台荷载及其周围流场数值模拟研究

内孤立波对海洋平台的安全运行存在一定的威胁,基于三维数值水槽对内孤立波传播引起的半潜平台受力及其周围流场的分布进行了数值模拟研究.通过与试验对比,分析了入射波幅和分层流体深度比对半潜平台上内孤立波荷载特性的影响规律,验证了数值模拟结果是准确可靠的.研究表明,内孤立波引起的平台荷载会随着内孤立波波幅的增加而增大,随着分层流体深度比的增加而减小.基于对平台周围剖面速度场和三维涡场的演化规律研究,发现在内孤立波传播过程中,平台周围会出现明显的速度减小区,平台周围会有大量的漩涡产生并发生脱落现象.     

内孤立波作用下Spar平台的运动响应

建立了内孤立波作用下Spar平台运动响应的时域数值模型,并利用这一模型计算了内孤立波作用下Spar平台的运动响应。其中根据K-dv方程解出内孤立波的波面及水质点速度表达式,应用Morison公式求出内孤立波对Spar平台的水平作用力;由伯努力方程求解压强,通过在平台底面积分求出内孤立波对Spar平台的垂向作用力;在悬链线理论的基础上,分析锚链对平台的作用力,采用四阶Runge-Kutta方法求解运动方程。数值结果表明,当内孤立波不断接近海洋平台时,Spar平台受到的作用力随之增大,海洋平台会产生远大于表面波作用下的水平位移,从而证实了内孤立波是影响海洋平台稳定的重要因素,在海洋平台的设计和评估中内孤立波的作用是不可忽略的。     

南海东北部C型内孤立波的观测与分析

基于布放在南海东北部陆坡海域的5套潜标观测到的内孤立波波列数据和孤立波扰动KdV(PKdV)理论,研究内孤立波在趋浅陆架上的传播特征。得出如下结果:1)观测到的内孤立波属于C型内孤立波,即平均重现周期为(23.41±0.31)h。2)内孤立波在西传爬坡过程中,其振幅表现为先增大后减小再增大,与该海域温跃层深度的变化趋势一致;由观测数据和理论计算得到的孤立波振幅增长率(SAGR)数值接近,表明该海域的内孤立波的振幅变化可以采用由孤立波PKdV方程导出的趋浅温跃层理论来描述。3)随着水深变浅,内孤立波传播方向向北偏移,传播速度减小,即在A,B和D站位,传播方向分别为279°,296°和301°,偏转角度达22°;传播速度分别为2.36,2.23和1.47 m/s,减小38%。     

内孤立波与Spar型浮式风机平台耦合作用研究

针对内孤立波与系泊浮式平台的耦合作用问题,本文以单立柱Spar型浮式风机平台为研究对象,研究了平台在内孤立波作用下的耦合运动特性。利用重叠网格技术及自编耦合动力程序,以计算流体动力学求解Navier-Stokes流体控制方程,建立了Spar型浮式风机平台与内孤立波(Internal solitary waves, ISWs)流固耦合作用数值模型。通过系列数值模拟,研究了Spar型浮式风机平台在内孤立波作用下的水平动态载荷和纵荡、横荡规律。研究表明,内孤立波对Spar平台会产生较显著的水平作用力,并产生大幅的水平位移;在水深一定时,同一流体分层比下,随着内孤立波振幅提高,纵荡运动幅值增大;当振幅一定时,随着流体分层比的增大,纵荡幅值减小;通过调整系泊系统,发现非线性系泊系统下平台纵荡幅值要大于线性系泊系统下的情况,即线性系泊系统会低估纵荡运动幅值。     

风作用下孤立波在珊瑚礁上传播变形机制数值模拟研究

基于二维不可压缩两相流模型建立了数值风浪水槽,采用SST k-ω雷诺时均湍流模型,研究了风作用下孤立波在珊瑚礁上的传播变形规律。将计算结果与实验数据对比,证明了该两相流模型计算孤立波在珊瑚礁上传播的准确性,并进一步分析了不同风速对珊瑚礁上孤立波传播变形的影响。结果表明：风的作用会使波面发生随机脉动特征。当地波高随风速的增大而增大;当地波高关于风速的变化梯度随入射波高的增大而增大。风的作用会加快孤立波的传播并且使孤立波提前发生破碎;孤立波开始破碎的位置随风速的增大向远离礁坪的方向移动。反射系数随风速的增大而增大;反射系数关于风速的变化梯度随入射波高的增大而减小;透射系数随风速的增大呈增大趋势。平底区波峰剖面同一水深处的水平流速随风速的增大而增大;且一定的风速不改变水平流速沿水深的变化梯度。有风时波面上方的矢量密度和大小均明显高于无风时且与风速呈正相关,并且波峰上方气流不再循环。随着风速的增大,水气交界面附近的正涡量和负湍流剪应力减小,负涡量和正湍流剪应力增大。水体动能、势能和总能达到高值的时间随风速的增大而减少;水体动能、势能和总能随风速的增大而增大,并且风速对水体动能的相对影响大于势能。     

基于FLUENT的内孤立波质量源造波方法

为研究内孤立波质量源数值造波方法,本文采用两个点源形式的质量源,分别放置于两层流体的上下层中作为内孤立波激发源。推导源项表达式,从不可压缩流体的Navier-Stokes方程出发,结合内孤立波Kd V、e Kd V理论,基于商业软件FLUENT发展了一种内孤立波质量源数值造波方法。通过数值模拟,分析了质量源造波过程中内孤立波的生成过程,并将数值模拟结果同理论及实验作对比。结果表明:基于此方法生成的内孤立波波形、波高及波致水平速度与理论及实验吻合度较好,该方法是可行的,并且耗时短、效率高。     

基于MCC理论的内孤立波数值模拟

基于MCC(Miyata-Choi-Camassa)理论,将内孤立波诱导上下层深度平均水平速度为入口条件,结合理想流体完全非线性欧拉方程,建立了两层流体中内孤立波生成与传播的CFD(Computational Fluid Dynamics)数值水槽方法。在系列数值模拟基础上,获得了内孤立波设计振幅与数值模拟振幅之间的相关关系,实现了在振幅可控条件下的内孤立波数值模拟。结果表明,在极限振幅范围内,数值模拟所得内孤立波波形均与MCC理论解吻合良好,而KdV(Korteweg-de Vries)和mKdV(modi-fied KdV)理论解只适用于小振幅的情况。同时,利用CFD数值模拟结果,对内孤立波诱导流速场特性进行了分析,结果表明内孤立波诱导水平速度在上层及波面下方流体层中垂向衰减很小,但在内界面与波面之间流体层中垂向衰减明显。     

两层流体中内孤立波与潜体相互作用数值模拟

提出基于双推板的内孤立波数值造波方法,对两层流体中内孤立波与潜体的相互作用进行了数值模拟,分析不同潜深下潜体所受的内孤立波载荷特性.结果表明,内孤立波对水下潜体的载荷作用是不可忽视的.特别地,当潜体位于内孤立波中时,其所受的垂向力要远比潜体位于内孤立波下时的大.     

内孤立波作用下Spar平台动力响应特性

以三类内孤立波理论(Kd V、e Kd V和MCC)的适用性条件为依据,采用Morison和傅汝德-克雷洛夫公式分别计算Spar平台内孤立波水平力和垂向力,结合时域有限位移运动方程,建立了有限深两层流体中内孤立波与带分段式系泊索Spar平台相互作用的理论模型。以东沙群岛某海域实测内孤立波为对象,数值分析了在内孤立波作用下某经典式Spar平台的内孤立波动态载荷、运动响应及其系泊张力的变化特性。研究表明,内孤立波不仅会对Spar平台产生突发性冲击载荷,使其产生大幅度水平漂移运动,而且还会使其系泊张力显著增大。因此,在Spar平台等深海平台的设计应用中,内孤立波的影响不可忽视。     

内孤立波波致流场数值模拟研究

基于KdV、mKdV理论,利用Fluent计算软件,采用"平板拍击"造波方法,进行内孤立波数值模拟,并与物理实验结果进行对比验证。利用数值模拟结果,分析内孤立波波致流场变化,结果表明:上下层流体中波致水平流速方向相反,均呈现先增加后减小的变化趋势,且波谷经过时刻流速最大;在波谷经过断面处,波致水平流速在上层流体中沿垂向分布无明显变化,在波面以下的下层流体中有衰减趋势,但衰减很小;两层流体界面与波谷之间存在过渡水深范围,水平流速在该水深范围内沿垂向衰减明显,且随内孤立波振幅的增大,过渡水深范围有所增大。     

有限深两层流体内孤立波非线性数值模拟入口边界条件研究

基于有限深两层流体KdV(Korteweg-de Vries)、eKdV(extended KdV)和MCC(Miyata-Choi-Camassa)理论,以内孤立波诱导上下层深度平均水平速度为入口边界条件,采用理想流体完全非线性欧拉方程,建立了两层流体中内孤立波生成的CFD(Computational Fluid Dynamics)数值模拟方法。以系列数值模拟结果为依据,结合内孤立波非线性和色散参数的组合条件,给出了选择合适内孤立波理论解作为CFD数值模拟入口边界条件的方法,从而实现了振幅与波形可控的内孤立波完全非线性数值模拟。     

局部微分求积法的深水包络孤立波数值模拟

利用局部微分求积法(LDQ)对非线性薛定谔(Schr dinger)方程进行数值求解,分别模拟了单深水孤立波运动,同向双深水孤立波追赶碰撞耦合运动,高阶孤立波振动和孤立波的反射与透射现象,得到各情况下的数值结果。从数值模拟及图像中揭示非线性薛定谔方程的性质和特点,阐述深水孤立波形成的物理意义、运动方式和运动规律,分析在不同初值条件下波形的变化特点,验证了LDQ法对该类问题的有效性。     

内孤立波作用下竖直圆柱体的横向力研究

本文对内孤立波作用下圆柱体的水动力特性进行了数值研究。对不同直径、不同水深比、不同浸没深度的圆柱进行了水平方向受力分析,从多个因素出发,研究了内孤立波对水下竖直圆柱的水平作用情况。通过分析水平作用力和力矩的时程曲线,得到内孤立波行进过程中对于圆柱体的横向作用力和力矩与时间的关系。通过对比不同工况下内孤立波橫向作用力幅值,得到内孤立波作用力幅值与圆柱直径、水深比以及浸没深度的关系。结果表明,内孤立波横向作用力(力矩)幅值与圆柱直径呈线性相关关系,但与水深比相关关系不明显,随着圆柱浸没深度的增加,内孤立波横向作用力幅值呈现出先增大后减小的趋势。     

南海北部陆架区内孤立波向岸传播过程研究

南海北部是全球海洋中内孤立波最强和最为活跃的海域。然而,内孤立波在传入陆架区后,其形态发生显著变化,其传播演变过程表现出高度的复杂性。本研究综合卫星图像和数值模式手段研究了内孤立波在向岸传播过程中的空间变化特征。可见光卫星图像研究结果显示,南海北部陆架区存在三种形态的内孤立波,分别为第一模态下凹型内孤立波、第一模态上凸型内孤立波和第二模态内孤立波。受水深和层结变化的控制,它们的分布区域显著不同。基于MITgcm的数值模拟研究表明,上凸型内孤立波由第一模态下凹内孤立波经过极性转换过程发展而来,而第二模态内孤立波由第一模态下凹内孤立波与急剧变浅地形相互作用而产生。     

内孤立波作用下FPSO动力响应特性

以KdV、eKdV和MCC理论的适用性条件为依据,利用内孤立波产生的瞬时速度和瞬时压力沿船体表面求面积分的方法分别计算浮式生产储油卸油装置(FPSO)内孤立波摩擦阻力和Froude-krylov力,再结合时域有限位移运动方程,研究了有限深两层流体中的内孤立波与带分段式系泊缆FPSO的相互作用问题。以东沙群岛某海域实测170 m振幅的内孤立波为对象,数值分析了在内孤立波作用下FPSO的内孤立波动态载荷、运动响应及其系泊张力的变化特性。研究表明,内孤立波会对FPSO产生突发的冲击,使FPSO产生大幅度水平漂移运动,导致系泊缆张力显著增大。在该振幅内孤立波作用下,FPSO受到的水平力最大可达56 t,垂向力最大可达1 961 t;纵荡位移可达47.62 m,系泊缆张力增加最大可达106.4 t。因此,在FPSO等深海浮式结构物的设计应用中,内孤立波的影响是不可忽视的。     

内孤立波数值模拟方法研究

本文选用RNGk-ε湍流模型封闭二维不可压缩N-S方程组,选用KdV、mKdV理论作为输入条件进行内孤立波造波模拟。运用流体模拟软件Fluent建立二维内波数值水槽,用动网格手段实现网格变形模拟物理边界造波法,用VOF方法进行密度分层流体的内界面追踪,对不同非线性强度的二相流内孤立波进行了数值模拟。并将数值模拟结果与KdV、mKdV理论波形以及内波试验水槽所造内孤立波进行对比,验证了此套造波方法在Fluent中进行内孤立波模拟的可行性。     

孤立波与直墙式多孔介质结构相互作用数值分析

以连续性方程和雷诺平均方程为控制方程,用k-ε模型封闭方程组,采用VOF方法追踪自由面,建立了具有造波和消波功能的二维孤立波数值水槽。并基于改进的Darcy多孔介质模型研究了二维孤立波与直墙式多孔介质结构的相互作用问题。数值分析了孤立波与直墙式多孔结构相互作用引起的流场和波形变化特性,并与文献中的实验结果进行了比较,数值结果与实验结果吻合良好。     

基于类海啸波型的岛礁水动力特性数值模拟研究

基于非静压单相流模型NHWAVE建立了高精度二维数值波浪水槽,采用日本2011年实测真实海啸波型系统研究了海啸波在岛礁上传播变形的规律,并且分析了波高、礁坪淹没水深和礁前斜坡坡度等因素对孤立波和真实海啸传播变形的影响。结果表明,相比孤立波,类海啸波的波长明显大于孤立波波长,在测点处引起的水面变化持续时间更长,同等波高情况下真实海啸波型比孤立波能够携带更多的能量,与岛礁的相互作用也更为复杂,在礁坪上形成的淹没水深约为孤立波的两倍。礁前斜坡坡度和礁坪淹没水深均对类海啸波的反射和透射系数有显著影响。随着礁前斜坡坡度的增加,反射系数和透射系数均逐渐增加。随着礁坪淹没水深的增加,反射系数逐渐减小,而透射系数逐渐增大。但是,反射系数和透射系数均随着入射波高的增加而逐渐减小。