基于无奇异边界元法模拟三维全非线性液舱晃荡（英文）

文章基于全非线性势流理论对三维液舱晃荡进行了数值模拟,其控制方程由无奇异边界积分方程法(Desingularized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺自由面。在微幅水平激励下,该文中得到的结果与解析解吻合较好。     

基于多次透射公式和无奇异边界元法模拟全非线性数值波浪水池

文章基于势流理论对全非线性的三维数值水池进行了模拟,其控制方程由无奇异边界积分方程法(Desingularized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺波面.同时,在远方辐射控制面上采用多次透射公式方法(Multi-transmitting Formula,MTF)来进行消波,文中得到的结果与理论解进行了比较,结果表明该方法可用来有效模拟全非线性的数值波浪水池.     

基于多次透射公式和无奇异边界元法模拟全非线性数值波浪水池（英文）

文章基于势流理论对全非线性的三维数值水池进行了模拟,其控制方程由无奇异边界积分方程法(Desingularized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺波面。同时,在远方辐射控制面上采用多次透射公式方法(Multitransmitting Formula,MTF)来进行消波,文中得到的结果与理论解进行了比较,结果表明该方法可用来有效模拟全非线性的数值波浪水池。     

直接边界元法时域全非线性兴波计算

文章在时域中采用直接边界元的方法计算模拟了以稳定航速运动的近水面潜体全非线性兴波现象.从格林定理出发,在边界面上布置Rankine源和偶极子,采用时间步进的方法计算模拟潜体近自由面运动的兴波现象.每个时步自由面上的势和节点位置采用全非线性拉格朗日自由面边界条件求解,重新生成自由面网格准备下一个时步计算的初始条件.文中给出了数值计算的结果并与试验和线性计算的结果相比,证实了该方法是可靠的,并可以利用该方法建立数值波浪水池.     

浅水中下潜物体有航速辐射问题的数值解

开发了一种基于非均匀有理B样条(NURBS)的Rankine源高阶面元法用于求解浅水中下潜物体有航速辐射问题。基于势流理论,流场中的扰动速度势用物体表面和自由面上分布的Rankine源来表达,同时用映像法计及水底的影响;采用NURBS曲面精确表达物面和自由面,在求出边界面上的未知源强后,物体表面的速度势分布用B样条表示;采用配置方法对边界积分方程进行求解,采用高斯—勒让德公式计算方程中的积分,同时开发了一种解决数值积分中奇异性问题的方法。采用文中数值方法求解了浅水中下潜圆球的有航速辐射问题,计算得到了下潜圆球水动力系数和自由面波形,所得水动力系数计算结果和他人的计算结果进行了比较,其吻合程度良好,从而验证了文中开发方法的有效性。     

基于N-S方程的航行船舶辐射问题数值模拟

基于N-S(Navier-Stokes)方程,进行了航行船舶辐射问题的数值模拟,包括强制垂荡和强制纵摇的模拟.计算了船舶垂荡和纵摇的附加质量和阻尼.数值模拟中,控制方程--BANS(Reynolds Averaged Navier-Stokes)方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理.文中给出了以不同航速前进的船舶强制垂荡和纵摇的力与力矩,以及船舶垂荡和纵摇的附加质量系数和阻尼系数,并与DUT(Delft University of Technology)的试验数据进行了比较,二者吻合良好.     

近水面椭球体非线性兴波的数值模拟

本文开发了一种去奇异面元法用于计算以定常前进速度运动的近水面椭球体引起的三维非线性自由面势流.采用求解非线性深水波的Stokes理论求解该定常非线性问题;物面用NURBS(非均匀有理B样条)表达;孤立的Rankine源分布在物体内部和自由面上方.给出了波型、波阻以及收敛特性随一些参数的变化规律.波阻计算结果和他人的其它数值结果进行了比较,其吻合程度令人满意.     

近自由面三维水翼的水动力分析及试验研究

采用面元法对近自由面三维水翼进行势流数值分析并进行了相关试验研究。在数值计算中,将Rankine源和偶极子置于边界面上,用时间步进法模拟水翼的势流场和自由表面波形。在自由面采用非线性自由面边界条件,在尾涡面上采用偶极子布置以满足Kutta条件。文中给出了数值计算模型的参数,对于不同浸深、不同航速和不同攻角下的水翼,计算了水翼表面上的压力分布,水翼的阻力和升力及自由表面波形。数值计算结果与试验结果进行了对比。结果表明,文中方法可用于水翼优化设计、近自由面振动翼运动及水翼船兴波等问题的研究。     

数值波浪水池及顶浪中船舶水动力计算

基于粘流理论建立了三维数值波浪水池,模拟了非线性波浪,并对规则波顶浪中前进的拘束船模的水动力进行了计算.数值模拟中,控制方程-RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波.给出了非线性规则波的模拟结果以及规则波顶浪中前进的拘束船模的水动力计算结果,并与理论解及DUT(Delfi University of Technology)的试验数据进行了比较,二者吻合良好.     

MPS方法在细长河道中溃坝流动的应用研究

溃坝问题是我国防灾减灾工作中的一个突出问题.溃坝发生时,会对下游的结构物产生剧烈的砰击作用,并伴随着自由液面的翻卷、破碎等强非线性现象.对这种短时间的自由面强非线性变化,基于网格的计算流体方法往往需要一些特殊的数值处理.文中采用本课题组自主开发的无网格粒子法求解器MLParticle-SJTU对二维带斜坡、三维带90°急弯的溃坝问题进行了模拟.研究了溃坝流的自由面演化,分析了溃坝流遭遇斜坡障碍物的复杂流动现象,并记录下监测点的波高时间历程.通过MPS的数值模拟结果与试验数据进行对比,两者吻合较好.     

三维完全非线性波浪水槽的数值模拟

用有限元求解拉普拉斯方程,采用满足完全非线性自由表面条件的半拉格朗日法追踪流体自由表面,离散单元采用20节点的六面体二次等参数单元,建立了三维完全非线性数值波浪水槽.把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播,为今后研究非线性波浪对船舶等非规则物体的作用提供了参考.     

Three-dimensional numerical wave tank simulations on fully nonlinear wave–current–body interactions

A finite-difference scheme and a marker-and-cell (MAC) method are used for numerical wave tank (NWT) simulations to investigate the characteristics of nonlinear wave motions and their interactions with a stationary three-dimensional body in the presence of steady uniform currents. The Navier–Stokes (NS) equation is solved in the computational domain, and the boundary values are updated at each time-step by a finite-difference time-marching scheme in the frame of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the marker–density function technique developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling the motions of a flexible flap wavemaker, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. Using the NS–MAC NWT, nonlinear wave and current interactions around a stationary vertical truncated circular cylinder are studied, and the results are compared with the experimental results of Mercier and Niedzwecki, a time-domain NWT based on linear potential theory, a fully nonlinear NWT, and a second-order diffraction computation. Received: July 3, 2001 / Accepted: September 25, 2001     

纵摇容器中液体晃荡的非线性数值模拟

在海洋工程领域,液体晃荡是一种普遍存在的物理现象。对于船舶而言,转动比平动有着更重要的影响。该文针对纵摇容器中的液体晃荡问题,采用高阶边界元方法建立自由水面满足完全非线性边界条件的时域数学模型。通过大地坐标系和随体坐标系之间的坐标变换,使得计算域仅控制在随体坐标系内。求解中采用半混合欧拉—拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验和数值结果的对比,验证了建立模型的准确性。进而开展大量数值试验研究容器纵摇运动频率、纵摇转动中心和容器中布置一竖直隔板对晃动波面与荷载的影响。     

全非线性方法对于沉浸双圆柱定常前行运动引起自由表面流的研究（英文）

The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.     

带有挡板的注液舱室地震响应分析（英文）

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace's equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method(BEM) is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

Seismic response of liquid-filled tank with baffles

In thispaper, the effects of a rigid baffle on the seismic response of liquid in a rigid cylindrical tank are evaluated. A baffle is an annular plate which supplies a kind of passive control on the effects of ground excitation. The contained liquid is assumed incompressible, inviscid and has irrotational motion. To estimate the seismic response, the method of superposition of modes has been applied. To analyze the rigid tank response, Laplace＆#39;s equation is considered as the governing equation of the fluid domain, in both time and frequency domains. The boundary element method （BEM） is employed to evaluate the natural modes of liquid in a cylindrical tank. To gain this goal, the fluid domain is divided into two upper and lower parts partitioned by the baffle. Linearized kinematic and dynamic boundary conditions of the free surface of the contained liquid have been considered.     

采用耦合有限元、边界元方法的矩形储液舱中液体晃荡数值仿真(英文)

Sloshing of liquid can increase the dynamic pressure on the storage sidewalls and bottom in tanker ships and LNG careers. Different geometric shapes were suggested for storage tank to minimize the sloshing pressure on tank perimeter. In this research, a numerical code was developed to model liquid sloshing in a rectangular partially filled tank. Assuming the fluid to be inviscid, Laplace equation and nonlinear free surface boundary conditions are solved using coupled FEM-BEM. The code performance for sloshing modeling is validated against available data. To minimize the sloshing pressure on tank perimeter, rectangular tanks with specific volumes and different aspect ratios were investigated and the best aspect ratios were suggested. The results showed that the rectangular tank with suggested aspect ratios, not only has a maximum surrounded tank volume to the constant available volume, but also reduces the sloshing pressure efficiently.