聚焦波作用下平面网衣结构的水动力特性研究

由于沿海区域的限制以及愈加严重的环境污染,渔业养殖正从近海走向深远海。深远海海域的海况更加恶劣,给养殖装备的设计与性能评估带来新的挑战。为解决该问题,对极端波浪与养殖装备网衣结构的相互作用开展研究。基于waves2Foam建立数值波浪水池,极端波浪模拟采用基于NewWave理论的聚焦波模型,网衣结构模拟采用多孔介质模型,并通过与Morison模型计算的网衣受力等效分析,获得多孔介质模拟网衣结构阻力系数的直接估计方法。然后将多孔介质模型嵌入waves2Foam中,开展聚焦波与网衣结构相互作用的数值模拟,同时开展水槽试验,验证数值模拟的准确性。基于数值模拟结果,系统地分析了不同网衣密实度及不同波浪参数下网衣结构的升阻力特性以及网衣结构对波浪场的扰动规律。研究表明：聚焦波波峰幅值和网衣密实度对网衣结构的升阻力影响较大,且升力峰值出现在阻力为0的时刻;网衣结构对聚焦波的时空演化特性有影响,改变了聚焦波波形。     

网衣对波浪传播影响的三维数值模拟

为研究网衣对波浪传播的影响,采用多孔介质模型模拟网衣,建立了模拟网衣与波浪相互作用的三维数值波浪水槽模型。基于该数值模型研究了波浪经过网衣作用后波高衰减的变化规律,并与物理模型试验结果进行比较。通过对比数值模拟和物理模型试验的结果,证明了该三维数值模型模拟网衣对波浪传播影响的可行性。     

波浪与抛石潜堤相互作用的MPS法数值模型研究

为了研究波浪与抛石潜堤相互作用过程中大自由表面变形和堤内渗流等强非线性紊流运动问题,利用改进的MPS法,建立了模拟波浪与抛石潜堤相互作用的MPS法数值计算模型。模型将抛石潜堤假定为均质多孔介质,采用Drew的二相流运动方程描述多孔介质内外的流体运动;通过在动量方程中增加非线性阻力项,并引入亚粒子尺度紊流模型,模拟波浪与可渗结构物相互作用过程中的紊流运动。选取“U”型管中多孔介质内渗流过程和孤立波与可渗潜堤相互作用两个典型的渗流问题,通过将数值计算结果与理论解和实测值的对比分析,对所提出的MPS法紊流渗流模型的模拟精度进行验证。结果表明:基于改进的MPS法构建的垂向二维紊流渗流模型可以很好地再现“U”型管中多孔介质内渗流以及波浪作用下可渗潜堤内外的复杂流场,显著缓解流-固界面处的压力震荡与粒子分布不均匀问题,实现了较高的模拟精度。     

基于IDDES方法的冲刷地形海底管道绕流数值模拟

近年来,海底管道对于近海油气工业中碳氢化合物的运输十分重要,预测海底管道周围水动力特征对于研究管道损坏机理具有重要意义。基于Spalart-Allmaras改进延迟分离涡模拟(improved detached eddy simulation,简称IDDES)湍流模型,利用开源OpenFOAM求解器计算单管冲刷过程中4个特征阶段管道周围水动力特征。在数值模拟中先计算了无圆管时均匀来流充分发展后的流场,依据该流场中进口流速分布设定有管道时的边界层及进口边界条件。通过将最终模拟结果与前人试验数据及数值模拟结果进行对比,表明了该模型计算结果相对于k-ε模型及LES模型更加贴合试验结果,验证了定义边界层方法的可靠性及IDDES模型在预测管道周围流场的适用性。模型计算结果显示,随着冲刷坑深度的逐渐增加,管道下游回流区会经历一个大幅缩减又逐渐扩大的过程,且逐渐关于圆管中心所在水平线对称;通过对不同冲刷阶段尾涡形态的分析表明,冲刷前期并没有涡脱产生,中期已经有产生涡脱的趋势,后期出现了明显的且较为稳定的涡脱。     

波浪作用下锚链系泊浮式防波堤动力响应的数值模型研究

浮式防波堤充分利用波能在水深方向的分布特性,在满足工程消浪要求的同时对海域水沙交换影响较小,且能够快速布置,在某些实际工程有一定应用前景。为了深入了解波浪作用下浮式防波堤的动力响应,基于OpenFOAM标准求解器olaFlow,在刚体运动求解计算中植入锚链求解模块MOODY(mooring cable dynamics),实现了基于重叠网格方法的浮体运动与锚链受力耦合求解,建立了锚链系泊浮式防波堤动力响应的二维数值模型。利用该数值模型对锚链系泊单方箱浮式防波堤在波浪作用下的透射系数、运动响应、锚链张力进行了模拟,并和相关试验结果进行了比较。结果表明,模型能够准确模拟二维波浪和浮式防波堤的相互作用,并用于三维模型的改进。     

深埋海底隧洞内水压力下饱和土的动力特性

将衬砌分别视为多孔柔性材料和弹性介质,在频率域内研究内水压力作用下饱和分数导数黏弹性土-深埋圆形隧洞多孔柔性或弹性衬砌系统的耦合简谐振动。土体的宏观力学特性采用多孔介质理论来模拟。通过引入与土体体积分数相关的应力系数,利用衬砌和土体界面处位移连续,分别得到饱和黏弹性土和衬砌的位移、应力和孔隙水压力等的解析表达式。在此基础上,分析了多孔柔性衬砌和弹性衬砌结构的差异,并考察了应力系数、渗透系数、分数阶导数本构参数等对系统动力响应的影响,结果表明:多孔柔性衬砌材料条件下系统的动力响应明显大于弹性衬砌材料条件下系统的动力响应;随着分数导数阶数和材料参数比的增加,系统共振效应明显减弱;衬砌边界透水和不透水只是反映边界渗透性的两种极限状态。     

基于有限元方法对鲆鲽网箱耐流特性的数值模拟

鲆鲽网箱结构在海中受到水流的冲击作用会发生运动与变形,针对鲆鲽鱼特有的栖底习性,为确保网底结构的稳定有必要对其进行动力分析。为此利用有限元方法建立了流场中网箱受力和变形的数学计算模型,运用该数学模型对底框中加有支撑管结构并装配方形网目网衣的鲆鲽网箱整体位移进行了数值模拟。数值模拟结果表明,网箱的网衣部分在水流作用下形态变化比较大,网底的水平位移与垂直位移随流速的增加而增大,而网箱的底框架在不同流速条件下均能保持在水平位置,且未发生明显的倾斜。由此可见,此类鲆鲽网箱具有较好的耐流性能。     

柔性网衣结构变形影响因素数值模拟研究

为解决柔性网衣在水流作用下容易产生较大变形,不利于网箱内鱼类生长的问题,本文利用Orcaflex软件对网衣结构变形进行了数值计算分析,并将计算结果与Kristiansen等^[10]的水槽试验和Cristian数值模拟结果进行比较,验证了数值计算的准确性,再利用“九点空间坐标法”计算了不同流速和配重情况下网箱体积剩余系数,定量分析了流速和配重对网衣变形的影响。结果表明：在低流速情况下,体积剩余系数可达0.997;在高流速情况下,体积剩余系数可达0.652;同时,在流速为0.3 m/s的情况下,增加40 g配重可以将网箱体积剩余系数从0.612提高到0.684,对网衣变形程度的影响高达10%。由此得出,流速是影响网衣变形的主要因素,流速越大,网衣变形情况越严重,在高速情况下,配重是改善网衣变形的重要因素。     

厦门港潮汐、风暴潮耦合模型

建立一个适合台湾海峡及小港湾,特别是能够反映小港湾的漫滩及潮汐、风暴潮耦合效应的数值模型,对小港湾风暴潮进行数值预报,有重要的社会和经济意义.文中以厦门港作为小港湾研究区域,在漫滩、风暴潮的耦合、开边界条件等方面对小港湾风暴潮数值计算作了一定的尝试,提出了一种能够模拟小港湾台风增水,潮位和流场变化情况的途径.     

潮流、风暴潮耦合模型推算珠江口海域极值流速

分别采用模式气压场及模式台风场计算气压分布及风速分布,采用地形拟合的sigma坐标水动力方程组及紊流方程建立了潮流、风暴潮耦合模型。三维水动力方程组通过分裂算子法与Casulli半隐式差分格式直接求解。在珠江口的极值流速计算中,先进行了南海北部的潮流、风暴潮耦合模拟,模拟水位与实测值符合较好。以六区域模拟结果为边界,采用嵌套网格方法,在珠江口进行了高分辨率的潮流、风暴潮耦合模拟。选用1968～1999年共65个影响工程区域的台风过程,计算出珠江口的海域年极值流速序列,分析了该海域多年一遇的极值流速分布。     

Numerical simulation of the effects of fish behavior on flow dynamics around net cage

The computational fluid dynamics study is performed to analyze the impact of the cultured fish on the flow field through net cage and the deformation of net cage. The shear stress turbulent k-omega model is applied to simulate the flow field through the net cage, and the large deformation nonlinear structure model is adopted to conduct the structural analysis of the flexible net cage. To validate the net-fluid interaction model of the net cage in current, a series of physical model tests are conducted, which indicate that the numerical model can accurately simulate the flow field around the net cage and the deformation of the net cage. A fish model is used to simulate the effect of fish behavior on the flow pattern around the net cage and the deformation of the net cage. In addition, the flow fields around the net cage in current are investigated considering different fish group structures, fish swimming speeds, fish distributions and fish stocking densities. The results indicate that the circular movement of fish in the still water leads to a low pressure zone at the center of net cage, which causes a strong vertical flow along the center line of the net cage. The drag force on the net cage is significantly decreased with the increasing fish stocking density, but the most severe deformation of net cage occurred in the case of medium fish stocking density.     

Incompressible SPH flow model for wave interactions with porous media

The paper presents an Incompressible Smoothed Particle Hydrodynamics (ISPH) method to simulate wave interactions with a porous medium. The SPH method is a mesh free particle modeling approach that is capable of tracking the large deformation of free surfaces in an easy and accurate manner. The ISPH method employs a strict incompressible hydrodynamic formulation to solve the fluid pressure and the numerical solution is obtained by using a two-step semi-implicit scheme. The ISPH flow model solves the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the NS type model equations for the flows inside the porous media. The presence of porous media is considered by including additional friction forces into the equations. The developed ISPH model is first validated by the solitary and regular waves damping over a porous bed and the solitary wave interacting with a submerged porous breakwater. The convergence of the method and the sensitivity of relevant model parameters are discussed. Then the model is applied to the breaking wave interacting with a breakwater covered with a layer of porous materials. The computational results demonstrate that the ISPH flow model could provide a promising simulation tool in coastal hydrodynamic applications.     

A coupled numerical model of wave interaction with porous medium

A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.     

河口海岸三维双σ坐标斜压水流数值模型研究

在海底地形陡变、垂向密度分层明显的水域,三维σ坐标模式中会出现一种"伪"水平斜压梯度力,并会引起"伪"密度流,以至于影响模拟的精度。垂向上引入双σ坐标变换,建立河口海岸水域三维斜压水流数值模型。数值试验结果表明,在海底地形陡变水域,双σ坐标模式可以减小水平斜压梯度力处理引起的误差。     

Moving particle method for modeling wave interaction with porous structures

This paper presents a numerical model for simulating wave interaction with porous structures. Incompressible smoothed particle hydrodynamics in porous media (ISPHP) method is introduced in this study as a mesh free particle approach that is capable of efficiently tracking the large deformation of free surfaces in a Lagrangian coordinate system. The developed model solves two porous and pure fluid flows simultaneously by means of one equation that is equivalent to the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the extended Forchheimer equation for the flows inside the porous media. Interface boundary between pure fluid and porous media is effectively modeled by the SPH integration technique. A two-step semi-implicit scheme is also used to solve the fluid pressure satisfying the fluid incompressibility criterion.The developed ISPHP model is then validated via different experimental and numerical data. Fluid flow pattern through porous dam with different porosities is studied and regular wave attenuation over porous seabed is investigated. As a practical case, wave running up and overtopping on a caisson breakwater protected by a porous armor layer are modeled. The results show good agreements between numerical and laboratory data in terms of free surface displacement, overtopping rate and pressure distribution. Based on this study, ISPHP model is an efficient method for simulating the coastal applications with porous structures.     

Numerical simulation of the influences of sinker weight on the deformation and load of net of gravity sea cage in uniform flow

1IntroductionThe fish farming and aquaculture industry areexpanding and the demand for suitable locations forfish farms is increasing.In the future,more of thefish farms will be located offshore,as the number ofsuitable nearshore locations is limited.Futu…     

抛弃式探头下沉运动数值计算的网格划分方法研究

各类抛弃式探头下沉运动的数值计算研究都需要对探头表面和整个计算域进行网格划分的前处理过程,选择不同的网格划分方法和划分精度,会影响数值计算速度和计算精度,最终生成的网格质量决定了计算结果的收敛性和准确性。针对抛弃式探头的复杂结构,对比了两种网格划分方法、三种网格划分精度,完成探头计算区域的网格划分,采用k-ε湍流模型,进行抛弃式探头下沉运动的数值计算,并将计算结果与水箱和水库实验结果进行对比,验证了混合网格划分方法和普通精度网格对抛弃式探头下沉运动数值计算的适用性,研究结果对类似较为复杂结构的水下运动体的数值计算前处理过程具有一定的参考和借鉴意义。