舰船非稳态气流的非结构化网格法

计算了LPD17两栖登陆舰甲板船尾的非稳态气流，证明可以通过计算的方法来可靠地预测舰船上部结构的气流。采用基于非结构网格的不可压缩流体求解程序计算非稳态气流，计算结果与风洞模型试验中选定位置的结果数据进行了比较，有很好的一致性。研究了风洞墙壁和舰船前进方向之间的影响。计算得到的非稳态流场既可以作为载人航行模拟器的输入，又可以在海军研究实验室洞穴式的虚拟现实环境中操纵。     

大规模并行非结构滑移网格守恒型隐式处理方法

滑移网格方法的难点是在不引入误差的情况下实现多区域网格之间的数值耦合,提高CFD求解器的并行可扩展性,以及改进交界面处并行通信关系。论文介绍了大规模并行非结构滑移网格守恒型隐式处理方法,并通过算例对滑移网格并行算法的鲁棒性和通量守恒性进行验证,它与串行计算结果具有良好的一致性。通过螺旋桨瞬态流场工程算例检验该滑移网格并行算法应用于船舶与海洋工程领域大规模工程问题的适用性、鲁棒性,以及并行性能。     

明渠交汇口水流浅水二维数学模型研究

针对明渠交汇口缓流流动特性,建立基于高性能有限体积数值方法、具有二阶精度的非恒定浅水二维水流数学模型。基本控制方程采用浅水二维水流方程;采用能够有效捕捉激波的高性能有限体积方法离散控制方程;采用了隐式方法处理阻力项,从而保证了模型的稳定性。采用非结构三角形网格系离散计算区域。最后,应用实测资料对两种模型进行了验证比较,结果表明该模型可模拟明渠交汇口基本流动特性,显示了对天然河流广阔的应用前景。     

三维动态混合网格在AUV发射过程中的应用

为研究AUV从有界流场自航发射到无界流场的运动边界的扰流场,文章采用了三维动态混合网格方法进行数值模拟的策略.混合网格由三菱柱,四面体/六面体网格构成,当AUV运动时,靠近AUV的三菱柱网格随AUV运动,外层是静止的六面体网格,中部的四面体网格随AUV运动而变形或者重构.数值仿真结果给出了不同时刻AUV表面的压力分布、整个航程AUV的阻力系数变化,其值与理论结果吻合.同时研究了直径比对发射管航行的附加质量和阻力系数的影响,这为水下对接AUV提供了有效手段.     

分形理论在絮凝形态学中的应用与展望

介绍了絮体结构模型的发展与分形,分析了我国仿真模拟研究的现状,对我国未来分形理论在混凝应用中的研究做出展望。     

长江感潮河段二维潮流数值模拟*

通过非结构网格建立了长江感潮河段二维潮流数学模型,研究了水工建筑物在非结构网格中的处理,验证结果良好;为检验模型长河段沿程流量的守恒性,模型上边界分别采用枯水、中水、洪水3种恒定流量模拟感潮河段水流运动。计算结果表明:模型沿程流量守恒性良好,模型能够有效模拟长江感潮河段水流运动。     

带偏舵航天飞机的网格生成及流场数值模拟

以OV-102航天飞机为对象，研究了有舵面偏角的全机单区结构网格的生成方法，模拟并分析了高超声速条件下全机及舵面气动特性。结果表明，带舵偏航天飞机单区结构网格生成方法及流场数值模拟方法，可正确模拟该复杂外形在中小舵偏条件下的气动特性。     

瓯江口波流耦合作用下二维流场数值模拟研究

基于MIKE21 SW海浪模式,建立非结构网格下的瓯江口二维、双向波流耦合数学模型,对该海区波流共同存在下的二维流场进行了数值模拟,计算结果与实测数据吻合良好,且优于未耦合的结果,比较分析表明：考虑潮流作用大大改善了波高和周期的计算结果,使之更接近观测结果;波浪对潮位的影响较小,对流速的影响较大,其中落潮憩息时刻,波浪对流速的影响最大.     

包不锈钢螺栓球网格结构应用初探

对纯不锈钢焊接球网格结构、纯不绣钢螺栓球网格结构、包不锈钢螺栓球结构三种型式进行比较,探讨不锈钢网格节点腐蚀原因与处理方法等问题.     

一种适用波浪潮流计算的地形插值算法

提出了一种求解结构化和非结构化网格节点上地形高度值的插值算法。该算法采用全局搜索，局部插值的方法处理数据，可根据需要控制地形数据的光滑度；同时充分利用了原始数据中的边界线信息．避免了由于缺乏陆域地形数据而引起的水陆交界处插值结果的振荡；通过引入背景矩形网格间接地插值求解非结构网格节点上的高度值，有效地降低了非结构网格插值算法的时间复杂度，减小了计算量。     

A CIP-based method for numerical simulations of violent free-surface flows

A CFD model is proposed for numerical simulations of extremely nonlinear free-surface flows such as wave impact phenomena and violent wave–body interactions. The constrained interpolation profile (CIP) method is adopted as the base scheme for the model. The wave–body interaction is treated as a multiphase problem, which has liquid (water), gas (air), and solid (wave-maker and floating body) phases. The flow is represented by one set of governing equations, which are solved numerically on a nonuniform, staggered Cartesian grid by a finite-difference method. The free surface as well as the body boundary are immersed in the computation domain and captured by different methods. In this article, the proposed numerical model is first described. Then to validate the accuracy and demonstrate the capability, several two-dimensional numerical simulations are presented, and compared with experiments and with computations by other numerical methods. The numerical results show that the present computation model is both robust and accurate for violent free-surface flows.     

基于间断有限元方法的可压缩水高速冲击平板特性研究

基于任意欧拉拉格朗日（ALE）方法,忽略流体的粘性,建立水射流冲击刚性平板的数值模型。为了更为精确地捕捉流场中压力波,该文推导了在ALE框架下的间断有限元方法。该方法易于提高数值离散的空间精度,数值稳定性较
　　好,利于精确模拟高速水射流冲击过程。对于自由液面的变形,采用径向基函数方法确定网格单元的速度。该方法采用边界网格节点的信息去推导出内部网格节点的信息,不需要单元信息。平头水柱射流冲击的数值结果与Autodyn的数值结果吻合较好。在验证方法合理性的基础上,文中对平头及圆头水柱中压力波的分布特性进行了分析。数值结果表明：当冲击在平板表面时,会在产生一个压力波之后向水域内传播,且使得平板受到较大压力的作用。此外,圆头水柱的射流冲击将会产生一个更大的压力峰值。     

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.     

Two-dimensional numerical simulation of an elastic wedge water entry by a coupled FDM-FEM method

Hydroelastic behavior of an elastic wedge impacting on calm water surface was investigated. A partitioned approach by coupling finite difference method (FDM) and finite element method (FEM) was developed to analyze the fluid structure interaction (FSI) problem. The FDM, in which the Constraint Interpolation Profile (CIP) method was applied, was used for solving the flow field in a fixed regular Cartesian grid system. Free surface was captured by the Tangent of Hyperbola for Interface Capturing with Slope Weighting (THINC/SW) scheme. The FEM was applied for calculating the structural deformation. A volume weighted method, which was based on the immersed boundary (IB) method, was adopted for coupling the FDM and the FEM together. An elastic wedge water entry problem was calculated by the coupled FDM-FEM method. Also a comparison between the current numerical results and the published results indicate that the coupled FDM-FEM method has reasonably good accuracy in predicting the impact force.     

Higher-order boundary element method for the interaction of a floating body with both waves and slow current

A hybrid boundary element method is suggested to solve the problem of the interaction of floating structures with both waves and slow current. A pulsating source and its mirror image referring to the sea bottom are adopted as the Green's function. The velocity potentials are expanded into an eigenfunction expansion in the outer region of the fluid domain while higher order elements are used to discretize the boundary surface surrounding the inner region. The method is validated by comparing calculated results for a circular cylinder with the semi-analytical solutions. The method is then applied to ellipsoids of various breadth and draft to investigate the influence of body shape on the wave drift damping.     

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

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.     

Numerical simulation of nonlinear wave force on a quasi-ellipse caisson

A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio B/D was shown to have significant influence on the non-dimensional transverse wave force.     

能量耗散效应的多域边界元法(英文)

The wave diffraction and radiation around a floating body is considered within the framework of the linear potential theory in a fairly perfect fluid.The fluid domain extended infinitely in the horizontal directions but is limited by the sea bed,the body hull,and the part of the free surface excluding the body waterplane,and is subdivided into two subdomains according to the body geometry.The two subdomains are connected by a control surface in fluid.In each subdomain,the velocity potential is described by using the usual boundary integral representation involving Green functions.The boundary integral equations are then established by satisfying the boundary conditions and the continuous condition of the potential and the normal derivation across the control surface.This multi-domain boundary element method(MDBEM) is particularly interesting for bodies with a hull form including moonpools to which the usual BEM presents singularities and slow convergence of numerical results.The application of the MDBEM to study the resonant motion of a water column in moonpools shows that the MDBEM provides an efficient and reliable prediction method.     

Research developments in numerical methods of fluid-structure interactions in naval architecture and ocean engineering北大核心CSCD

流固耦合问题较为复杂,通常难以通过理论推导求得,而数值模拟则能提供一种有效的解决方案,并被广泛用于船舶与海洋工程领域。流固耦合数值方法根据其网格离散方式,可以分为贴体网格方法、非贴体网格方法、重叠网格方法和粒子类方法 4类,对这4类方法的特点及研究进展进行概述并总结得出:贴体网格方法和重叠网格方法均能精确捕捉界面的变形和演化,适合高雷诺数流动问题,在考虑结构变形时一般采用贴体网格方法,而考虑复杂几何形状的刚体运动时则常采用重叠网格方法;非贴体网格方法能够避免网格的更新操作,使计算较为简单,目前多用于模拟流动控制、水下柔性仿生航行器的研发以及多体运动干扰等问题;粒子类方法因其固有的拉格朗日属性,在模拟涉及自由液面剧烈变形、砰击、爆炸等强非线性流固耦合问题中发挥着重要作用。不同的流固耦合问题属性决定了不同方法的适用性,如何选取适合的数值方法,同时结合各类方法的优势开发新的计算方法以应对更为复杂的问题,是流固耦合算法开发的重要发展方向。     

水翼涡激振动的数值模拟研究

对二维水翼结构进行流固耦合运动分析,利用大涡模拟方法计算高雷诺数下水翼绕流场,流场力作用于二自由度刚体上导致周期性的垂荡和转动,龙格库塔法求解刚体水翼的运动方程,位移参数作为下一时间步流场计算的边界条件,具体通过编译自定义函数控制刚体运动和流场网格变化。探讨了初始攻角、刚心位置以及来流速度对水翼振动的影响,发现在来流速度增大至某一数值时水翼发生了颤振现象,并在时域与频域上对颤振的发生机理进行深入探讨。