基于无网格MPS方法数值分析方舱破舱进水问题

采用移动粒子半隐式法(moving particle semi-implicit,简称MPS)对自由漂浮二维方舱的破舱瞬时进水过程进行数值模拟。首先,采用基于GPU平台自主开发的MPS软件模拟破舱进水问题,并与其他方法得到的数值模拟结果进行对比验证。然后,对该二维方舱的各种模型进行了数值模拟,其中开孔位于不同位置以表示舷侧不同高度下的损坏。此外,还研究了不同类型的挡板对破舱进水后方舱稳定性的影响。结果表明损坏的孔洞和内部挡板会影响损坏舱段的运动特性,开孔距静水面的距离越大引起舱段的横摇等运动幅度越大,垂直挡板比水平挡板对舱内洪水的影响更大。     

Numerical prediction of green water incidents

Green water loads on moored or sailing ships occur when an incoming wave significantly exceeds the freeboard and water runs onto the deck. In this paper, a Navier–Stokes solver with a free surface capturing scheme (i.e. the VOF model; [Hirt and Nichols, 1981]) is used to numerically model green water loads on a moored FPSO exposed to head sea waves. Two cases are investigated: first, green water on a fixed vessel has been analysed, where resulting waterheight on deck, and impact pressure on a deck mounted structure have been computed. These results have been compared to experimental data obtained by [Greco, 2001] and show very favourable agreement. Second, a full green water incident, including vessel motions has been modelled. In these computations, the vertical motion has been modelled by the use of transfer functions for heave and pitch, but the rotational contribution from the pitch motion has been neglected. The computed water height on deck has been compared to the experimental data obtained by [Buchner, 1995a] and it also shows very good agreement. The modelling in the second case was performed in both 2-D and 3-D with very similar results, which indicates that 3-D effects are not dominant.     

Shared-Memory parallelization of consistent particle method for violent wave impact problems

A shared-memory parallelization is implemented to the recently developed Consistent Particle. Method (CPM) for violent wave impact problems. The advantages of this relatively new particle method lie in four key aspects: (1) accurate computation of Laplacian and gradient operators based on Taylor series expansion, alleviating spurious pressure fluctuation and being able to model two-phase flows characterized by large density difference, (2) a thermodynamics-based compressible solver for modelling compressible air that eliminates the need of determining artificial sound speed, (3) seamless coupling of the compressible air solver and incompressible water solver, and (4) parallelization of the numerical model based on Open Multi-Processing (OpenMP) and a parallel direct sparse solver (Pardiso) to significantly improve computational efficiency. Strong and weak scaling analyses of the parallelized CPM are conducted, showing an efficiency speedup of 100 times or more depending on the size of simulated problem. To demonstrate the accuracy of the developed numerical model, three numerical examples are studied including the benchmark study of wave impact on seawall, and our experimental studies of violent water sloshing under rotational excitations and sloshing impact with entrapped air pocket. CPM is shown to accurately capture highly deformed breaking waves and violent wave impact pressure including pressure oscillation induced by air cushion effect.     

Theoretical and experimental study on dynamic behavior of a damaged ship in waves

Most of the large scaled casualties are caused by loss of structural strength and stability due to the progressive flooding and the effect of waves and wind. To prevent foundering and structural failure, it is necessary to predict the motion of the damaged ship in waves.This paper describes the motion of damaged ship in waves resulting from a theoretical and experimental study. A time domain theoretical model, which can be applied to any type of ship or arrangement, for the prediction of damaged ship motion and accidental flooding has been developed considering the effects of flooding of compartments. To evaluate the accuracy of the model, model tests are carried out in ship motion basin for three different damaged conditions: engine room bottom damage, side shell damage and bow visor damage of Ro–Ro ship in regular and irregular waves with different wave heights and directions.     

A finite element method for dynamic analysis of long slender marine structures under combined parametric and forcing excitations

This paper presents a numerical analysis of lateral responses of a long slender marine structure under combined parametric and forcing excitations. In the development of the 3-D numerical program, a finite element method is implemented in the time domain using the Newmark constant acceleration method. Some example studies are performed for various water depths, environmental conditions and vessel motions. The relative amplitudes of combined excitations to a conventional forcing excitation are examined. The response amplitude of a combined excitation is much greater than that of a forcing excitation in the even number of instability regions of the Mathieu stability chart. The results demonstrate that a combined excitation needs to be considered for the accurate dynamic analysis of long slender marine structures subjected to a surface vessel motion.     

Numerical modeling of perforated plates in oscillating flow

Perforated plates, relevant for several marine applications, are experimentally and numerically investigated. The numerical investigations are performed using a presently developed Navier–Stokes solver. Several comparison and sensitivity studies are presented, in order to validate and verify the solver. Forced heave experiments are performed on two perforated plates with perforation ratios 19% and 28%. Amplitude-dependent added mass and damping coefficients are presented. Good agreement is obtained between the solver and the present experiments. Consistent with existing data, the results show that the hydrodynamic coefficients of perforated plates are highly amplitude dependent. The damping force is found to dominate over added mass force. The damping force dominance increases with increasing perforation ratio. It is highlighted that plate-end flow separation has an important effect on the damping coefficient. The developed numerical solver is two-dimensional, but is found to yield reasonable estimates of hydrodynamic force coefficients when compared with a previous three-dimensional experimental investigation. This could indicate that three-dimensional effects are not dominant for the hydrodynamic forces of perforated plates, and that a two-dimensional viscous flow solver could have relevance as a tool for estimating hydrodynamic forces on three-dimensional perforated structures.     

A Numerical Study on Wave-Mud Interaction

1 .Introduction The dynamics of soft mud under surface water waves is of great importance to the sedimentationprocesses in approach channels and harbors ,and has long been drawing attention. Advancements innumerous engineering applications inthe shallowco…     

Comparison of numerical solution and semi-empirical formulas to predict the effects of important design parameters on porpoising region of a planing vessel

Stability of the marine vessels in different conditions is one of the most important problems in the design of a planing vessel. In this research, the effects of some important design parameters (mass, longitudinal center of mass, deadrise angle, and length) of DTMB 62 model 4667-1 planing hull on the drag and also on the longitudinal dynamic stability (porpoising) are investigated in the velocity range of 2.12–8.486 m/s in calm water. In this paper, both numerical simulation of Reynolds Average Naiver Stokes (RANS) equations and semi-empirical formulas of Savitsky are used to analyze the motion of a 4667 planing vessel in calm water with two degrees of freedom (2DOF). For this purpose a finite volume, ANSYS-FLUENT, code is used to solve the Navier-Stokes equations for the simulation of the flow field around the vessel. In addition, an explicit VOF scheme and SST- Kω model is used with dynamic mesh scheme to capture the interface of a two-phase flow and to model the turbulence respectively, in 2DOF model (heave and pitch). Also, the results of both methods are compared with each other. According to the present results, changing the aspect ratio of the vessel and also the longitudinal center of gravity have the most effect on the porpoising region.     

Motion response simulation of a twin-hulled semi-submersible

The purpose of the study was to develop a prediction technique to simulate the motion response of a damaged platform under wave, wind and current forces. The equations of motion were obtained using Newton's second law and the numerical solution technique of non-linear equations of motion is explained for intact and damaged cases. The analysis technique employs large displacement non-linear equations of motion. Solutions were obtained in the time-domain to predict the motion characteristics. In this study, analysis procedures were developed to calculate: (a) wave loading on asymmetrical structural configurations; (b) hydrodynamic reaction forces (inertia or moment of inertia, damping and restoring forces) on asymmetrical shapes. During the damage simulation, change in the mass of the structure as well as wave and hydrodynamic reaction forces, were taken into account. The computer program developed for the time-domain simulation is introduced. In order to avoid slowly decaying transient motions of the structure due to wave excitation forces, an exponential ramp function is used. The application of a ramp function enables a quick convergence in the time-domain solution of equations of motion. Results of a numerical motion simulation program and the experimental studies are also presented in order to make comparisons. Comparison of the test results with the numerical simulations shows good agreement for heave, roll and pitch motions. The formulations and the computational procedures given in this paper provide useful tools for the investigation of the non-linear dynamic stability characteristics of floating structures in waves for intact, damaged and post-flooding conditions in six-degrees of freedom.     

The phase difference effects on 3-D structure of wave pressure acting on a composite breakwater

In designing the coastal structures, the accurate estimation of the wave forces on them is of great importance. In this paper, the influences of the phase difference on wave pressure acting on a composite breakwater installed in the three-dimensional (3-D) wave field are studied numerically. We extend the earlier model [Hur, D.S., Mizutani, N., 2003. Coastal Engineering 47, 329–345] to simulate 3-D wave fields by introducing 3-D Navier–Stokes solver with the Smagorinsky's sub-grid scale (SGS) model. For the validation of the model, the wave field around a 3-D asymmetrical structure installed on a submerged breakwater, in which the complex wave deformations generate, is simulated, and the numerical solutions are compared to the experimental data reported by Hur, Mizutani, Kim [2004. Coastal Engineering (51, 407–420)]. The model is then adopted to investigate 3-D characteristics of wave pressure and force on a caisson of composite breakwater, and the numerical solutions were discussed with respect to the phase difference between harbor and seaward sides induced by the transmitted wave through the rubble mound or the diffraction. The numerical results reveal that wave forces acting on the composite breakwater are significantly different at each cross-section under influence of wave diffraction that is important parameter on 3-D wave interaction with coastal structures.     

岛礁附近浮式平台运动响应特性数值分析

利用基于开源平台Open FOAM自主开发的船舶与海洋工程水动力性能求解器naoe-FOAM-SJTU,数值模拟了近岛礁环境下规则波的演化特性以及带有系泊系统的浮式平台在相应波浪作用下的水动力性能。对于平台的水动力性能的研究发现,仿真结果与试验结果在平台自由衰减运动固有周期及RAO(response amplitude operator)方面吻合良好。对于波浪在近岛礁地形下的演化现象的研究,分析了波浪演化不对称性特性的成因,并分别给出了不同参数下波浪在地形上爬升时演化的具体特性,对于波高变化及波浪演化的频率成分进行量化的探究。研究发现,波浪周期越大,波高变化越明显,演化的非线性现象越明显,且波浪随着传播距离的增大演化出的阶数也在增大。     

Numerical and Experimental Research on the Global Performances of Cell-Truss Spar Platform

Spar technology has been applied to the deep-sea oil and gas exploitation for several years.From the first generation of classic spar,the spar platform has developed into the second generation of truss spar and the latest cell spar.Owing to its favorable adaptability to wide range of water depth and benign motion performances,spar has aroused quite a lot of interests from oil companies,universities and research institutes.In the present paper,a new cell-truss spar concept,put forward by the State Key Laboratory of Ocean Engineering(SKLOE)at Shanghai Jiao Tong University,is studied both numerically and experimentally.The numerical simulation was conducted by means of nonlinear time-domain fully coupled analysis,and its results were compared to the experimental data.Whereafter,detailed analysis was carried out to obtain the global performances of the new spar concept.Proposals for the improvement of numerical calculation and experimental technique were tabled meanwhile.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

Wave Refraction-Diffraction in the Presence of a Current

-Wave refraction-diffraction due to a large ocean structure and topography in the presence of a 'current are studied numerically. The mathematical model is the mild-slope equation developed by Kirby (1984). This equation is solved using a finite and boundary element method. The physical domain is devid-ed into two regions: a slowly varying topography region and a constant water depth region. For waves propagating in the constant water depth region, without current interfering, the mild- slope equation is then reduced to the Helmholtz equation which is solved by boundary element method. In varying topography region, this equation will be solved by finite element method. Conservation of mass and energy flux of the fluid between these two regions is required for composition of these two numerical methods. The numerical scheme proposed here is capable of dealing with water wave problems of different water depths with the main characters of these two methods.     

Liquid moment amplitude assessment in sloshing type problems with smooth particle hydrodynamics

Sloshing moment amplitudes in a rectangular tank for a wide range of rolling frequencies are investigated both experimentally and numerically. In a previous paper, Souto et al. [2004. Simulation of anti-roll tanks and sloshing type problems with smoothed particle hydrodynamics. Ocean Eng. 31 (8–9), 1169–1192] numerical results obtained with a 3-D Smooth Particle Hydrodynamics (SPH) formulation were presented. These only corresponded to the phase lag between the tank motion and the liquid response moment. This paper is aimed at improving those results by obtaining accurate values for the moment amplitudes. We present the corrections with respect to the aforementioned implementation that focus on the time integration scheme and on the treatment of the boundary conditions. In addition better quality experimental results are presented.     

Simulation for the propulsion of a micro-hydro robot with an unstructured grid

The flow mechanism of contractive and dilative motion was numerically investigated to obtain a propulsive force in a highly viscous fluid. The computing program for the analysis of complicated motions was numerically developed with a cell-centered, unstructured grid scheme. The developed program was validated by the well-known equation of an oscillating plane below viscous fluid for an unsteady problem, which is known as Stokes’ second problem. Validation has continued through comparison with the experimental results.In this case, sinusoidal motion was applied to the validation, instead of trochoidal motion, because it was very difficult to actually simulate trochoidal motion in this experiment. Finally, the validation and comparison with the nodal-point scheme was accomplished by Stokes’ problem, which is the famous problem at a low Reynolds number. The validated code was applied to contractive and dilative motion in a narrow tube, whose motion was embodied by trochoidal movement. In a highly viscous fluid, such as a very sticky honey or a swamp, the computed results show that a viscous force can be used for propulsion instead of a dynamic force.From the present results, it was found that a propulsive force can be obtained by contractive and dilative motion at a low Reynolds number, which can be applied to the propulsion of micro-robots in a highly viscous fluid such as a blood vessel or a swamp. This research could also be considered fundamental research for the propulsion of micro-hydro robots, which are expected to be actively studied in the future in accord with further development of nanotechnology.     

新型直立式透空堤消浪性能数值研究

为使防波堤同时具有良好的掩护效果和水体交换能力,提出了两种带有透浪通道的新型直立式防波堤。基于Fluent求解器建立了三维数值波浪水槽,通过与试验结果对比,验证了该数值水槽求解波浪与透空堤作用具有较高的精度。对两种防波堤在规则波作用下的透浪特性进行了研究,结果表明:透射系数K_t与透空率呈正线性相关,且可通过调整透浪通道间距,使相同透空率下K_t降低20%～30%。对同一结构,K_t随相对波长的增大而显著增大,但受相对波高的影响较小。在透空率大于0.16后,异型沉箱防波堤的消浪性能明显优于错位沉箱。基于数值计算结果,给出了以上两种透空堤波浪透射系数的经验公式。     

A combined wind and wave energy-converter concept in survival mode: Numerical and experimental study in regular waves with a focus on water entry and exit

A combined wind and wave energy converter concept, named STC concept was proposed. Model tests were performed in terms of operational and survival modes. Water entry and exit phenomena as well as green water on deck were observed during the survivability model tests. In this paper, a nonlinear numerical model based on a blended station-keeping potential-flow solver with a local impact solution for bottom slamming events and an approximated model for the water shipped on the deck is proposed to simulate these nonlinear phenomena. Physical investigation of the water entry and exit process was firstly carried out and uncertainty analysis of the model test results were performed. Numerical comparisons between the nonlinear solver and model test results are then performed in terms of mean, wave frequency and double wave frequency motion response components. The slamming and green water involved in the water entry process are specially investigated, in terms of the physical evolution and the effects on the dynamic motion responses. The validation work on the occurrence of slamming and water on deck as well as the slamming pressure are performed.     

Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

非线性波浪时域计算的三维耦合模型

将计算区域Ω划分为内域Ω₁和外域Ω₂(Ω₂=Ω-Ω₁),外域控制方程采用改进线性频散特性的二维Boussinesq方程,用预报一校正法数值求解;结构物附近的内域控制方程为三维Navier-Stokes方程,由VOF方法数值求解。通过在外域和内域相匹配的交界面上设置合适的速度和波面边界条件,建立了三维非线性波浪时域计算的耦合模型。模拟试验表明:(1)耦合模型数值波浪水池可以产生稳定的、重复性较好的波动过程;(2)用耦合模型数值波浪水池求解较大浅水区域上的非线性波浪数值计算问题可以取得较高的计算效率,同时又能得出结构物附近的复杂流场。