论文部分内容阅读
波浪与可渗结构的相互作用包含波浪反射、透射、越浪和破碎等水动力过程,涉及到湍流、渗流和耗散等多种非线性物理效应,是一个作用机理十分复杂的非线性渗流力学问题。近年来随着全球气候变暖和海平面上升,海岸防护结构面临的海洋灾害威胁也日益加剧,开展波浪与可渗结构的相互作用研究具有十分重要的科学意义和工程应用价值。本文在修正光滑粒子流体动力学(SPH)方法核函数连续性的基础上,通过提出-种新的固壁粒子压力计算模式,解决了原SPH方法存在的边界处粒子分布紊乱、固壁边界附近压力场震荡以及流体粒子与固边界分离或穿透固边界问题。基于MPI空间域分解和OpenMP数据分解的耦合算法,结合动态负荷平衡技术,提出了适合现代计算集群资源配置的混合MPI-OpenMP并行算法,解决了三维SPH数学模型计算效率低的瓶颈问题。提出了适合无网格粒子法的主动吸收式造波边界条件、周期性侧边界条件和人工粘性消波边界条件,建立了基于SPH算法的周期波数值波浪水槽/水池。基于所建立的SPH数值波浪水槽/水池,结合非线性渗流本构关系,建立了模拟非线性波浪在多孔介质结构中流动的二维和三维湍流模型。推导了多孔介质内流体运动的体平均/Favre平均形式的双平均湍流控制方程(VAFANS),湍流应力采用亚粒子紊流模式来模拟,引入非线性渗流阻力项以反映介质固体骨架对孔隙内流体的阻力作用。通过在计算域内设置定义孔隙率信息的背景网格和引入表征密度的概念,将自由流动区域和不同孔隙率的渗流区域内的流体运动控制方程形式统一,避免在整个计算域内使用不同形式的控制方程导致的不同介质界面附近的压力震荡及速度跳跃现象。提出了不同介质界面处的过渡区边界条件和适合的过渡区宽度,保证了多孔介质界面处流场和压力场的计算精度。数值计算结果表明:所建立的二维SPH渗流模型比其它数值模型能更精确地模拟出溃坝波在可渗坝体内的传播、波浪在不同可渗性海床上的衰减、波浪透过可渗潜堤或出水堤时的变形、翻卷、破碎和越浪过程以及可渗结构附近的流场和湍流能分布特性。所建立的三维SPH渗流模型能够精确模拟溃坝波冲击可渗方柱、孤立波在可渗突堤处的绕射,以及周期波作用下分段式离岸潜堤附近的三维水动力特性,可精确捕捉到波浪在可渗结构周围复杂的波面演化过程以及平均水位和波生流空间分布特性。