Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure

As a gridless particle method, the MPS (Moving Particle Semi-implicit) method has proven useful in a wide variety of engineering applications including free-surface hydrodynamic flows. Despite its wide range of applicability, the MPS method suffers from some shortcomings such as non-conservation of momentum and spurious pressure fluctuation. By introducing new formulations for the pressure gradient and a new formulation of the source term of the Poisson Pressure Equation (PPE), and by allowing a slight compressibility, we have proposed modified MPS methods for the prediction of wave impact pressure on a coastal structure. The improved performance of the modified methods is shown through the simulation of numerous wave impact problems (including the impacts by a dam break flow, a flip-through and two cases of slightly-breaking waves) in comparison with the experimental data.     

New modified gradient models for MPS method applied to free-surface flow simulations

In this paper, two modified pressure gradient models based on Taylor series expansion are proposed to enhance the higher order source term MPS (MPS-HS) method. The modified models consist of gradient correction matrices applied to the existing (base) pressure gradient models. To validate the modified pressure gradient models first hydrostatic pressure test is simulated and compared to both the base and modified MPS methods. Using the modified models are shown to reduce unphysical pressure oscillations observed in the base models. Second, an evolution of an elliptical drop in a 2D flow field is examined and shown to verify the models. Third, the proposed models illustrated appropriate stability and consistency properties against analytical solutions when an altered gravitational acceleration was superimposed to the hydrostatic pressure test. In addition, an improved performance is observed when Higher order Laplacian (HL) and Error-Compensating Source (ECS) of the Poisson Pressure Equation (PPE) schemes are coupled with the modified pressure gradient models compared to coupling them with the base gradient models. Finally, the modified MPS methods enhanced performances are validated in a free-surface flow simulation for a dam break problem with impact pressure, and a violent sloshing flow in a rectangular tank when compared to the base MPS methods against an existing experimental data.     

Derivation of Reliability Index Vector Formula for Series System and Its Application

In this study, a reliability index vector formula is proposed for series system with two failure modes in term of the concept of reliability index vector and equivalent failure modes. Firstly, the reliability index vector is introduced to determine the correlation coefficient between two failure modes, and then, the reliability index vector of a series system can be obtained. Several numerical cases and an analysis on offshore platform are performed, and the results show that this scheme provided here has better computational accuracy, and its calculation process is simpler for the series systems reliability calculations compared with the other methods. Also this scheme is more convenient for the engineering applications.     

Derivation of reliability index vector formula for series system and its application

In this study, a reliability index vector formula is proposed for series system with two failure modes in term of the concept of reliability index vector and equivalent failure modes. Firstly, the reliability index vector is introduced to determine the correlation coefficient between two failure modes, and then, the reliability index vector of a series system can be obtained. Several numerical cases and an analysis on offshore platform are performed, and the results show that this scheme provided here has better computational accuracy, and its calculation process is simpler for the series systems reliability calculations compared with the other methods. Also this scheme is more convenient for the engineering applications.     

Assimilation of simulated float data in Lagrangian coordinates

We implement an approach for the accurate assimilation of Lagrangian data into regional general ocean circulation models. The forward model is expressed in Lagrangian coordinates and simulated float data are incorporated into the model via four-dimensional variational data assimilation. We show that forward solutions computed in Lagrangian coordinates are reliable for time periods of up to 100 days with phase speeds of 1 m/s and deformation radius of 35 km. The position and depth of simulated floats are assimilated into the viscous, Lagrangian shallow water equations. The weights for the errors in the model and data are varied and the assimilation results react appropriately. We show the effect of different spatial and temporal samplings of float data on all Lagrangian trajectories in the computational domain. At the end of the assimilation period, results from the Lagrangian shallow water equations could be interpolated and used as initial and boundary conditions in an Eulerian general ocean circulation model.     

Multi-resolution MPS for incompressible fluid-elastic structure interactions in ocean engineering

The paper presents a multi-resolution MPS (Moving Particle Semi-implicit)-based FSI (Fluid-Structure Interaction) solver for efficient and accurate simulations of incompressible fluid flows interacting with elastic structures. The fluid model is founded on the projection-based MPS solution of continuity and Navier-Stokes equations. The structure model is set based on MPS-based discretization of linear and angular momenta corresponding to an isotropic elastic solid. Fluid-structure coupling is conducted in a mathematically-physically consistent manner along with implementation of a multi-resolution scheme comprising of common radius of influence, revised weight function, revised number density and potential number density concepts to enhance i) consistency of particle-based discretizations, ii) imposition of boundary conditions and iii) volume conservation at fluid-structure interface. A set of previously developed enhanced schemes are also adopted for the fluid model. The robustness and efficiency of proposed Enhanced Multi-resolution MPS-MPS FSI solver are investigated through a set of ocean engineering-related benchmark tests. To the best knowledge of authors, this study presents the first multi-resolution particle method for FSI corresponding to incompressible fluid and elastic structures.     

渤海三维潮余流的数值计算(英文)

作者根据冯士筰的三维弱非线性理论,计算了渤海三维潮致拉格朗日余流,其中包括由冯等人引进的作为二阶量的拉格朗日漂移椭圆。文中不仅给出了三维结果,也给出了深度平均的结果,最后还与实际图象进行了比较。     

An accurate and efficient SPH modeling of the water entry of circular cylinders

The present work is dedicated to the enhancement and application of the recently developed δ⁺-SPH scheme for two and three dimensional water entry problems. The cylinder surface presents a circular shape and therefore the position where the flow separation happens is tightly related to the implementation of the boundary condition. A special treatment for the particle shifting between fluid and solid wall is highlighted. Adaptive particle refinement (APR) is applied in this work to reduce the computational cost. It is found that, APR and the δ⁺-SPH scheme benefit to each other. That is to say the former reduces the computing cost of the latter while in return the latter solves the problem of particle disorder of the inactive particles of the former. Thanks to the combination of APR and δ⁺-SPH especially in three dimensional (3D) cases, the overall computational cost is significantly reduced while sufficiently fine particle resolution can be obtained in the flow region characterized by large pressure gradient close to the structure surface. The fairly good agreements between the SPH results and the experimental data prove the present SPH model to be a reliable tool in accurately solving the fluid–structure interacting problems in ocean engineering.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

MPS法高精度自由表面识别方法研究

为了提高粒子法即移动粒子半隐式法(Moving Particle Semi-implicit method, MPS法)中自由表面粒子的识别精度,降低由粒子误判引起的非物理压力振荡,对MPS法的自由表面识别方法进行了改进。在原始自由表面判别标准的基础上增加辅助判别条件,提出2种新的自由表面判别法即压力判别法和填充率判别法。利用对静水问题和溃坝流问题的模拟计算,对比分析选用不同自由表面判别法得到的数值计算结果,揭示粒子识别精度对压力计算的重要影响。研究结果表明:新提出的压力判别法和填充率判别法可以有效地提高自由表面粒子的识别精度,减轻压力计算中的非物理压力振荡现象,从而提高压力计算的稳定性以及整体数值计算的模拟精度。     

渤海三维污染物迁移轨迹的数值计算

根据一种三维潮致拉格朗日余流的弱非线性理论,考虑单一频率潮波系统,在零阶天文潮和一阶拉格朗日余流数值计算的基础上,计算了二阶拉格朗日余流,然后对8个具有代表性的标识质点运动轨迹进行了数值跟踪。对水平及垂直运动性质进行了分析,并用卫片资料对结果进行了验证。提出一些有启发性的初步结果。     

Civil 3D在港区工程量计算中的应用

港区陆域形成的填挖方量对工程经济性影响非常大,寻找精确地计算港区陆域工程量使纵向设计最优化的方法是港区工程全寿命周期各个环节的关键。基于传统方法的计算原理与应用现状分析了现有方法在港区土方计算中存在的问题,借鉴建筑信息模型BIM(Building Information Modeling)思想,对Civil3D的土方计算原理与操作流程进行研究,最后以某港区石化仓储库区场平方案为例,参考传统计算软件HTCAD的计算结果,验证了Civil3D在港区工程量计算中的便捷性与可靠性。结果显示:Civil3D的计算精度比HTCAD高7%~20%,计算速度是HTCAD的数十倍,在修改过程中计算速度达数百倍。研究成果为采用Civil3D进行港区陆域土方计算提供理论和技术依据。     

A conservative nonlinear filter for the high-frequency range of wind wave spectra

Filtering of the high-frequency part of a wind wave spectrum may be useful in a numerical wind wave model for various reasons. First, it can be used to augment (or be part of) a parameterization of the resonant nonlinear interactions, that are essential to third-generation wind wave models. Second, when combined with a dynamic time stepping scheme for source term integration, it may result in smoother (and hence faster) wave model integration. In this study, such a filter is proposed, based on the traditional Discrete Interaction Approximation (DIA) for the resonant four-wave nonlinear interactions. This filter retains all conservative properties of the interactions. For small time steps and/or smooth spectra, it is formulated as a traditional source term. For larger time steps and/or non-smooth spectra it is formulated as a filter. This formulation guarantees stability of the filter itself and will enhance overall computational stability in a full wave model. The stability properties of this filter are illustrated using traditional wave growth computations. Examples are given where the filter improves model economy, and where it is shown to remove spurious high-frequency noise from a wave model.     

Generation of waves in Boussinesq models using a source function method

A method for generating waves in Boussinesq-type wave models is described. The method employs a source term added to the governing equations, either in the form of a mass source in the continuity equation or an applied pressure forcing in the momentum equations. Assuming linearity, we derive a transfer function which relates source amplitude to surface wave characteristics. We then test the model for generation of desired incident waves, including regular and random waves, for both one and two dimensions. We also compare some model results with analytical solution and available experiment data.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

近海海上风电场水下噪声传播模型适用性研究

通过现场采集近海海上风电场工程区运营期风机水下噪声和背景噪声数据,计算了噪声信号的倍频带声压级,功率谱级和峰值声压级,确定了海上风电场水下噪声总声源级为148.3 d B,以此开展近海海上风电工程风机水下噪声频域特性、功率密度谱特性等研究。在此基础上使用Kraken简正波模型和Bellhop射线模型对风电场运营期风机水下噪声在水平与垂直方向上的传播进行模拟,模拟了噪声在不同频带内的衰减程度,结果显示模型模拟结果在不同频率下的衰减趋势有着很大差异,产生了明显的多途干涉现象,通过实测数据对建立的噪声传播模型进行验证,发现Kraken简正波模型在500 Hz以下,Bellhop射线模型在500 Hz以上适合模拟实际水下噪声传播情形,同时海区本身背景噪声的存在会对预测的准确性产生影响。这些结论可用于进一步对近海海上风电场水下噪声传播的研究。     

Simplified Method for Calculating Standing Wave Pressure on Vertical Breakwater

Through numerical modeling, a kind of simplified calculating method for standing wave pressure on vertical face breakwater have been put forward. Not only the formulas proposed in this paper are simple in form and very easy in use, but also they possess continuity on the full range of standing wave. And more, the precision requiremennts of calculation can be satisfied to a certain extent in engineering practice.     

On the internal generation of waves: control volume approach

The literature contains theoretical derivation of the internal wave generation techniques by using control volume approach. The source term in three-dimensional continuity equation is derived by Reynolds transport theorem. The source term in two-dimensional wave model is derived by work-energy equation. The results confirmed that the appropriate velocity to be included in the source term is physically the group velocity which is the energy velocity for a wave model.     

Radiative transfer modeling and analysis of spatially variant and coherent illumination for undersea object detection

Increasing the optical range of target detection and recognition continues to be an area of great interest in the ocean environment. Light attenuation limits radiative and information transfer for image formation in water. In this paper, the authors briefly review current methods of imaging and then describe a variation of the spatial interferometric technique that relies upon projected spatial gratings with subsequent detection against a coherent return signal for the purpose of noise reduction and image enhancement. A model is developed that simulates the projected structured illumination through turbid water to a target and its return to a detector. The model shows an unstructured backscatter superimposed upon a structured return signal. The model has been extended to predict what a camera would actually see, so that various noise-reduction schemes can be modeled. Finally, some water-tank tests are presented, validating original hypothesis and model predictions. The method is advantageous in not requiring temporal synchronization between reference and signal beams and may use a continuous illumination source. Spatial coherency of the beam allows for the detection of the direct return, while scattered light appears as a noncoherent noise term.     

A Lagrangian particle random walk model for simulating a deep-sea hydrothermal plume with both buoyant and non-buoyant features

This paper presents a computational model of simulating a deep-sea hydrothermal plume based on a Lagrangian particle random walk algorithm. This model achieves the efficient process to calculate a numerical plume developed in a fluid-advected environment with the characteristics such as significant filament intermittency and significant plume meander due to flow variation with both time and location. Especially, this model addresses both non-buoyant and buoyant features of a deep-sea hydrothermal plume in three dimensions, which significantly challenge a strategy for tracing the deep-sea hydrothermal plume and localizing its source. This paper also systematically discusses stochastic initial and boundary conditions that are critical to generate a proper numerical plume. The developed model is a powerful tool to evaluate and optimize strategies for the tracking of a deep-sea hydrothermal plume via an autonomous underwater vehicle (AUV).