[发明专利]一种网箱对水流影响的模拟方法

权利要求书

1.一种网箱及其内部养殖物对水流影响的模拟方法，其特征在于：包括以下步骤：

A、建立数值水槽

采用右手笛卡尔坐标系，坐标原点为网箱中心在水面上的投影，x轴在水面上，定义水流方向为x轴正方向，y轴在水平面上与水流方向垂直，z轴铅直向上；数值水槽左端定义为速度入口边界；右端定义为自由出流边界；水槽侧壁、水槽底面以及水面定义为固壁边界，剪切应力为0；结合连续方程和动量方程作为控制方程来描述流体的运动，基于有限体积法离散控制方程；

流体运动的控制方程如下：

连续方程：

∂ρ∂t+∂(ρui)∂xi=0---(1)]]>

动量方程：

∂(ρui)∂t+∂(ρuiuj)∂xj=-∂ρ∂xi+∂∂xj(μ∂ui∂xj-ρui'uj'‾)+Si---(2)]]>

式中，t为时间，ρ为流体的密度，μ为流体的动力粘度，u_i、u_j为流体速度分量的时均值，u_i’、u_j’为速度分量的脉动值，p为压力的时间平均值，i、j等于1、2、3，1、2、3分别表示坐标分量x、y、z；上划线“-”表示对物理量取时间平均；

S_i是动量方程源项，在多孔介质边界外的流体区域S_i=0，在多孔介质边界内部，

Si=-(Dijμu+Cij12ρ|u|u),]]>

Dij=Dn000Dt000Dt,Cij=Cn000Ct000Ct---(3)]]>

式中，D_ij、C_ij为多孔介质系数矩阵，D_n表示法向粘性阻力系数，D_t表示切向粘性阻力系数，C_n表示法向惯性阻力系数，C_t表示切向惯性阻力系数，u为流体速度；

B、建立网箱模型

根据网箱形状，将网箱划分为若干片平面网衣，引入多孔介质模型模拟网衣；多孔介质模型是在多孔介质区域耦合了一个根据经验假设为主的流动阻力，当水流以一定的流速流过多孔介质区域时，作用在该区域的阻力值F由下式计算：

F＝S_xλA   (4)

式中，S_x为x方向的动量方程源项，λ为多孔介质厚度，A为多孔介质面积；F方向与水质点速度方向相反，公式(4)中不考虑阻力值F的方向；

将方程(3)代入方程(4)得到水阻力F_d与升力F_l的表达式，

Fd=(Dnμu+Cn12ρ|u|u)λA---(5)]]>

Fl=(Dtμu+Ct12ρ|u|u)λA---(6)]]>

多孔介质系数通过物理模型实验数据计算得到；当缺少实验数据时，对于与水流方向垂直的平面网衣，其多孔介质系数通过以下经验公式计算得到：

Dn=2.281×107Sn2-2.432×106Sn]]>

Dt=1.129×106Sn2+3.194×105Sn---(7)]]>Cn=188.375Sn2+33.163Sn]]>

Ct=135.55Sn2+4.083Sn]]>

式中，S_n为网衣密实度；

对于与水流方向成一定冲角α的平面网衣，多孔介质系数由公式(7)计算后，还须再按照以下公式进行计算，得到对应冲角网衣的多孔介质系数：

Dn'=Dn+Dt2+Dn-Dt2cos(2α')---(8)]]>Cn'=Cn+Ct2+Cn-Ct2cos(2α')]]>

Dt'=Dn-Dt2sin(2α')---(9)]]>

Ct'=Cn-Ct2sin(2α')]]>

其中，α′＝90°-α，α为冲角，是水流方向与平面网衣之间的夹角；

通过对不同冲角下平面网衣模型多孔介质系数的设置，实现对圆形网箱对水流影响的数值模拟；

C、建立养殖物模型

对网箱及其内部养殖物周围流场的数值模拟引入多孔介质模型模拟养殖区域，采用有限体积法求解控制方程，对网箱及其内部养殖物周围的流场特性进行数值计算；

C1、基本假设

为了建立网箱内养殖物的数学模型，根据养殖鱼类的特征，作如下三个基本假设：

a、养殖鱼类同等大小，并把鱼体假设成圆球；

b、养殖鱼类在网箱内均匀分布，养殖区域各向同性；

c、鱼体密度等于水的密度，因为正常情况鱼体在水中处于悬浮状态；

C2、养殖物模型

将网箱内的整个养殖体积定义为一个柱形多孔介质模型，多孔介质系数需要由以下公式计算得到；

Dn=150Dp(1-ϵ)2ϵ3---(10)]]>

Cn=3.5Dp1-ϵϵ3---(11)]]>

其中，

Dp=(6Mnπρw)13---(12)]]>

ϵ=1-ρfρw---(13)]]>

式中，D_p为鱼体模型的直径，ε为养殖区域孔隙率，M为网箱内养殖鱼类的总重，n为养殖鱼类的数目，ρ_f为网箱内鱼类的养殖密度，ρ_w为水的密度；

因为假设养殖区域各向同性，所以D_t=D_n，C_t=C_n；

将网箱模型和养殖物模型加入数值水槽中，分别设置这两个模型的多孔介质系数，即实现网箱及其内部养殖物对水流影响的模拟。