一种柔性网衣对波浪传播影响的模拟方法技术

本发明专利技术公开了一种柔性网衣对波浪传播影响的模拟方法，包括以下步骤：采用推板造波的方法，按照造波方程形式造波，结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程，通过流体体积函数方法求解波面的方式建立数值波浪水槽；数值模拟中引入多孔介质模型模拟网衣，多孔介质模型是在定义为多孔介质的区域内耦合一项由经验假设确定的流动阻力，利用人工附加的流动阻力模拟现实中网衣的阻尼效应。根据本发明专利技术的物理实验验证发现数值模拟得到的消浪结果与实际吻合，数值模拟的计算时间在10个小时以内；本发明专利技术较好地解决了现有技术中成本高、用时长、无法给出网衣附近的波浪传播变化情况的问题。

【技术实现步骤摘要】

本专利技术属于水产养殖工程、水利工程和海洋工程
，特别涉及到一种波浪作用过程中柔性网衣对波浪传播影响的数值方法。
技术介绍
浮式防波堤不仅能防浪，而且对周围的生态环境影响很小，基本上能保持原有海洋生态环境的平衡，也不会影响潮流来往及浮堤内外的海水交换。目前在浮式防波堤设计中，利用浮板和网衣组合形成浮式防波堤来消波被广泛使用，其中利用网衣结构破坏水体表层波浪运动，吸收波浪能量来消减波浪传播是板网式防波堤消浪工作的主要工作原理。目前，国内外利用物理模型实验方法作为开展网衣消浪行为的主要研究手段。由于网衣是特种海洋工程结构物，它完全不同于通常用于港口、海上油田等海工结构物，具有大变形、柔性、多孔的特点，物理模型实验模拟成本高、用时长，同时由于目前观测技术的限制，很难给出网衣附近的波浪传播变化情况。
技术实现思路
为解决现有技术存在的上述问题，本专利技术的目的是，提供一种既可模拟网衣周围波浪形态和网衣对波浪传播影响，又可以降低成本、缩短时间的柔性网衣对波浪传播影响的模拟方法。 为实现上述目的，本专利技术的技术方案如下:，包括以下步骤:A、数值波浪水槽的建立采用推板造波的方法，按照造波方程形式造波，结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程，通过流体体积函数方法求解波面的方式建立数值波浪水槽；为了消除波浪反射，数值波浪水槽引入动量源项消波方法；Al、控制方程在笛卡尔坐标系下，定义X为波浪传播方向，I为铅垂向上；流体运动的控制方程采用粘性不可压缩流体的Navier-Stokes方程:连续方程:权利要求1.，其特征在于:包括以下步骤: A、数值波浪水槽的建立 采用推板造波的方法，按照造波方程形式造波，结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程，通过流体体积函数方法求解波面的方式建立数值波浪水槽；为了消除波浪反射，数值波浪水槽引入动量源项消波方法； Al、控制方程 在笛卡尔坐标系下，定义X为波浪传播方向，I为铅垂向上；流体运动的控制方程采用粘性不可压缩流体的Navier-Stokes方程: 连续方程:全文摘要本专利技术公开了，包括以下步骤采用推板造波的方法，按照造波方程形式造波，结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程，通过流体体积函数方法求解波面的方式建立数值波浪水槽；数值模拟中引入多孔介质模型模拟网衣，多孔介质模型是在定义为多孔介质的区域内耦合一项由经验假设确定的流动阻力，利用人工附加的流动阻力模拟现实中网衣的阻尼效应。根据本专利技术的物理实验验证发现数值模拟得到的消浪结果与实际吻合，数值模拟的计算时间在10个小时以内；本专利技术较好地解决了现有技术中成本高、用时长、无法给出网衣附近的波浪传播变化情况的问题。文档编号G06F19/00GK103235878SQ201310130029公开日2013年8月7日 申请日期2013年4月15日 优先权日2013年4月15日专利技术者赵云鹏, 毕春伟, 董国海 申请人:大连理工大学本文档来自技高网...

【技术保护点】
一种柔性网衣对波浪传播影响的模拟方法，其特征在于：包括以下步骤：A、数值波浪水槽的建立采用推板造波的方法，按照造波方程形式造波，结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程，通过流体体积函数方法求解波面的方式建立数值波浪水槽；为了消除波浪反射，数值波浪水槽引入动量源项消波方法；A1、控制方程在笛卡尔坐标系下，定义x为波浪传播方向，y为铅垂向上；流体运动的控制方程采用粘性不可压缩流体的Navier?Stokes方程：连续方程：∂u∂x+∂v∂y=0---(1)动量方程：∂u∂t+u∂u∂x+v∂v∂y=-1ρ∂p∂x+fx+υ(∂2u∂x2+∂2u∂y2)+Sx???????????????????????????(2)∂v∂t+u∂v∂x+v∂v∂y=-1ρ∂p∂y+fy+υ(∂2v∂x2+∂2v∂y2)+Sy式中，t为时间，ρ为流体的密度，υ为流体的运动粘度，u、v分别为x、y方向的速度分量，p为压力的时间平均值，fx、fy分别为x、y方向的单位质量力分量，Sx、Sy分别为x、y方向的附加动量源项；A2、造波方程采用摄动展开法求解推板运动产生的波面，并对正弦运动的推板附加二阶运动项，以消除数值波浪水槽中的谐波，造波方程形式为：ξ=-an1tanhkh0[cosωt+a2h0n1(34sinh2kh0-n12)sin2ωt]---(3)式中，n1=12(1+2kh0sinh2kh0)方程中，a为波幅，k为波数，ω为波浪圆频率，t为时间，h0为水深；由推板运动产生波面形式为：η=acos(kx-ωt)+a2k4coshkh0sinh3kh0(cosh2kh0+2)cos2(kx-ωt)---(4)A3、边界条件及网格设置边界条件设置：数值波浪水槽左端为动边界，利用动网格技术按照造波方程实现推板造波，中间段为计算区域，长度大致为10倍波长，右端为消波区域，消波区域长度l为2倍波长，上边界为压力边界；波长L＝2π/k，式中波数k由色散关系ω2＝gktanhkh0确定，其中g为重力加速度；数值波浪水槽的消波方法采用动量源项消波方法：Sx=-(L′-L0l)2u,Sy=-(L′-L0l)2v---(5)式中，Sx与Sy分别为动量方程的附加源项，u、v分别为x、y方向的速度分量，l为消波段的长度，L0、L′分别为消波段起始位置的x方向坐标与附加源项处的x方向坐标；网格设置：水槽长度方向网格步长Δx为波长的1/50，高度方向对静水面附近网格进行加密且步长Δy小于1/10的波高，同时保证网格的长宽比小于2，时间步长Δt为波浪周期的1/2000；A4、数值计算数值水槽中气液界面的求解采用流体体积函数方法，利用有限体积方法对计算域进行离散，对流项的差分选用精度较高的二阶迎风格式，时间项的差分选用二阶隐式差分方法，对于压力?速度耦合方程的求解采用适于非稳定状态的压力隐式算子分割算法；B、网衣模型的建立数值模拟中引入多孔介质模型模拟网衣，多孔介质模型是在定义为多孔介质的区域内耦合一项由经验假设确定的流动阻力，利用人工附加的流动阻力模拟现实中网衣的阻尼效应；本质上，多孔介质模型是在动量方程中叠加基于达西定律的动量源项Si，在多孔介质边界外的流体区域Si=0，在多孔介质边界内部，Si=-(Dijμu+Cij12ρ|u|u),?????????????????????????????????????(6)Dij=Dn00Dt,Cij=Cn00Ct式中，μ为流体的动力粘度，Dij、Cij分别为多孔介质系数矩阵，Dn表示法向粘性阻力系数，Dt表示切向粘性阻力系数，Cn表示法向惯性阻力系数，Ct表示切向惯性阻力系数；多孔介质模型的流动阻力值F由下式计算：F＝SxλA???(7)式中，Sx为x方向的动量方程...

【技术特征摘要】

【专利技术属性】
技术研发人员：赵云鹏，毕春伟，董国海，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：