模拟二维水下爆炸问题的光滑粒子(SPH)方法

水下爆炸涉及两相流,是一个动边界的问题,迅速变化的气-液交界面给传统的基于网格的数值方法带来了严重的挑战。SPH法是一种拉格朗日方法,该方法可以准确且有效地捕捉气-液交界面。论文应用光滑粒子(SPH)法来模拟水下爆炸过程,特别给出了光滑长度的选取方法和优化算法。最后,通过两个算例验证该方法,数值模拟结果表明,光滑粒子法在处理水下爆炸这样的大变形和强非线性问题方面具有精度高、耗时少、数值稳定等优点。     

Smoothed particle hydrodynamics and its applications in fluid-structure interactions

In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena(e.g. water waves, impacts, splashing jets, etc.). To model these complex free surface flows is a tough and challenging task for most computational fluid dynamics(CFD) solvers which work in the Eulerian framework. As a Lagrangian and meshless method, smoothed particle hydrodynamics(SPH) offers a convenient tracking for different complex boundaries and a straightforward satisfaction for different boundary conditions. Therefore SPH is robust in modeling complex hydrodynamic problems characterized by free surface boundaries, multiphase interfaces or material discontinuities. Along with the rapid development of the SPH theory, related numerical techniques and high-performance computing technologies, SPH has not only attracted much attention in the academic community, but also gradually gained wide applications in industrial circles. This paper is dedicated to a review of the recent developments of SPH method and its typical applications in fluid-structure interactions in ocean engineering. Different numerical techniques for improving numerical accuracy, satisfying different boundary conditions, improving computational efficiency, suppressing pressure fluctuations and preventing the tensile instability, etc., are introduced. In the numerical results, various typical fluid-structure interaction problems or multiphase problems in ocean engineering are described, modeled and validated. The prospective developments of SPH in ocean engineering are also discussed.     

二维滑坡涌浪的SPH方法数值模拟

块体滑坡往往引起自由水面的剧烈变化，研究块体滑坡激发的水面波动一直是水利与海岸工程界非常感兴趣的问题。该文使用SPH方法（光滑质点水动力学）对水下块体下滑引起的自由表面水动力学过程进行了二维数值模拟，并将数值计算结果与试验数据作了对比，计算结果给出了块体下滑过程中，水体表面出现的二次卷破现象，以及在块体上方出现两个方向相反的旋涡，显示了SPH方法对处理块体下滑引起的自由表面大变形问题具有十分优越的特性。通过数值试验，文中还讨论了滑块下滑速度对自由表面破碎和形成的旋涡个数的影响。     

SPH方法在自由面问题中的应用

作为无网格粒子法,SPH法在处理大变形、自由面流动问题时具有显著的优势.介绍了SPH法的基本数值方法,并基于SPH法数值模拟了2个二维溃坝问题,将计算结果与试验数据进行了比较,结果表明:SPH法在处理自由面时具有很强的适应性.尽管水面发生了翻卷、破碎等剧烈的变化,但SPH法仍然能够较好地捕捉到这些流动现象,同时数值模拟得到的水头位置和自由面形状均能与试验结果相吻合,表明SPH法在处理自由面问题时具有较高的准确性及可靠性.     

On the modeling of viscous incompressible flows with smoothed particle hydrodynamics

Smoothed particle hydrodynamics(SPH) is a Lagrangian, meshfree particle method and has been widely applied to different areas in engineering and science. Since its original extension to modeling free surface flows by Monaghan in 1994, SPH has been gradually developed into an attractive approach for modeling viscous incompressible fluid flows. This paper presents an overview on the recent progresses of SPH in modeling viscous incompressible flows in four major aspects which are closely related to the computational accuracy of SPH simulations. The advantages and disadvantages of different SPH particle approximation schemes, pressure field solution approaches, solid boundary treatment algorithms and particle adapting algorithms are described and analyzed. Some new perspectives and future trends in SPH modeling of viscous incompressible flows are discussed.     

Restoring Ice-jam Floodwater to a Drying Delta Ecosystem

Abstract

Overbank flooding is essential to the ecological health of riparian landscapes, particularly river deltas. One of the world's largest freshwater deltas, the Peace-Athabasca Delta (PAD) in northern Canada, has experienced a series of wetting and drying cycles because of inter-annual variations in flooding. Recent research has found that most of the major floods affecting this system are produced by spring ice jams. For approximately two decades, however, the combination of climatic and flow-regulation effects precluded significant ice-jam flooding of the PAD. Resultant drying caused major changes to the ecology of the delta and led to the evaluation of a number of methods to restore water flows. Since most of delta is contained within a national park (Wood Buffalo National Park), a major goal was to employ non-structural measures. Hence, in an effort to manage the water problems of this delta, the final report of a multi-agency “Northern River Basins Study” made the recommendation that the spring flow-release strategy of the upstream hydro electric reservoir be modified to increase the probability of ice-jam flooding near the PAD. This was to be conducted in years when downstream hydrometeorological conditions (snowpack magnitude and ice-cover strength) appeared conducive to ice-jam formation. Such favourable conditions developed in the spring of 1996, a natural ice jam began to develop, and regulated flows were increased to assist in potential flooding. As a result, the PAD experienced its first major flood in over 20 years. This paper reviews the hydrometeorological conditions that led to the ice-jam formation, compares the conditions to historical events, analyzes the spatial extent of the flood, and evaluates the effectiveness of the flow release.     

孤立波浅化过程的SPH数值模拟

本文应用光滑粒子流体动力学(SPH)方法建立了二维数值波浪水槽,研究了孤立波沿斜坡的传播、变形和破碎过程。通过数值计算分析了波浪沿斜坡传播过程中自由表面的变化特征,给出了斜坡上不同位置处波面变化的历时曲线,并应用文献记录的试验结果对数模结果进行了验证。得到了波高的沿程变化趋势、破波点的位置和破碎波高以及不同破碎形态下波浪自由表面的变化特征,显示了SPH方法在处理较复杂自由表面流动问题方面的优越性。     

高速船甲板上浪的水体流动与冲击

针对高速船舶的甲板上浪,综合考虑船体外飘、船艏波面升高、波浪传播与船体运动等相关因素的影响,采用势流理论计算和CFD技术相结合的方法,建立了计入整船运动的船艏上浪计算模型,计算模拟并分析了甲板上水体流动形态及其对甲板结构的冲击作用。计算中船体运动通过势流理论给定,势流理论计算通过用户自定义函数嵌入Fluent软件,并通过二次开发实现了上浪现象的模拟。文中就迎浪状态下高速行驶的S175船甲板上浪进行了三维数值模拟研究,结果表明甲板上水体的流动及其对结构物的冲击载荷与试验吻合良好,该方法能够对上浪机理给出较为清晰的描述,能够分析预报甲板上浪对浮体结构物的破坏作用,为工程应用提供了良好的方法和途径。     

A 3-D SPH model for simulating water flooding of a damaged floating structure

With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during the transient flooding process. Therefore, an enhanced smoothed particle hydrodynamics(SPH) model is applied in this paper to investigate the response of a simplified cabin model under the condition of the transient water flooding. The enhanced SPH model is presented firstly including the governing equations, the diffusive terms, the boundary implementations and then an algorithm regarding the coupling motions of six degrees of freedom(6-DOF) between the structure and the fluid is described. In the numerical results, a non-damaged cabin floating under the rest condition is simulated. It is shown that a stable floating state can be reached and maintained by using the present SPH scheme. After that, three-dimensional(3-D) test cases of the damaged cabin with a hole at different locations are simulated. A series of model tests are also carried out for the validation. Fairly good agreements are achieved between the numerical results and the experimental data. Relevant conclusions are drawn with respect to the mechanism of the responses of the damaged cabin model under water flooding conditions.     

High-order finite difference solution for 3D nonlinear wave-structure interaction

This contribution presents our recent progress on developing an efficient fully-nonlinear potential flow model for simulating 3D wave-wave and wave-structure interaction over arbitrary depths (i.e. in coastal and offshore environment). The model is based on a high-order finite difference scheme OceanWave3D presented in [1, 2]. A nonlinear decomposition of the solution into incident and scattered fields is used to increase the efficiency of the wave-structure interaction problem resolution. Application of the method to the diffraction of nonlinear waves around a fixed, bottom mounted circular cylinder are presented and compared to the fully nonlinear potential code XWAVE as well as to experiments.

     

Numerical modeling of dam-break flood through intricate city layouts including underground spaces using GPU-based SPH method

This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas including underground spaces. Taking the advantage of GPUs parallel computing techniques, simulations involving more than 107 particles can be achieved. We use a virtual geometric plane boundary to handle the outermost solid wall in order to save considerable video card memory for the GPU computing. To evaluate the accuracy of the new GPU-based SPH model, qualitative and quantitative comparison to a real flooding experiment is performed and the results of a numerical model based on Shallow Water Equations （SWEs） is given with good accuracy. With the new GPU-based SPH model, the effects of the building layouts and underground spaces on the propagation of dambreak flood through an intricate city layout are examined.     

UNSTEADY FREE-SURFACE WAVES GENERATED BY BODIES IN A VISCOUS FLUID

BRIEFINTRODUCTIONOFTHEPAPER :　Inthefirstpartofthisthesis ,thefundamentalsolutionsforthesingularStokesandOseenflowsinanun boundedfluidarederivedinauniversalformwhichinvolvestheHamiltonian ,Hessian ,andLaplacianoperators ,andelementaryfunctions .Inthesecondpart,theinteractionofunsteadylow Reynolds numberflowswithafreesurfaceisinvestigatedanalytically .Thedisturbedflows ,generatedbysubmergedbodiesmovingverticallydownwardsawayfromthesurfaceofthefluid ,aregovernedbytheunsteadyStokeseq…     

A two‐dimensional self‐adaptive hydrodynamic scheme for the assessment of the effects of structures on flooding phenomena in river basins

In order to assess the effects of river and floodplain engineering projects on flooding, a new self‐adaptive hydrodynamic scheme for the simulation of two‐dimensional river flows is proposed. The depth‐averaged motion equations are solved numerically using a fractional step method, in which the convective terms are calculated using the inverse characteristics method and the remaining terms with an explicit method based on a finite difference method. The integration is performed on a dynamically self‐adaptive calculus grid, which allows representation of the movable boundary between wetting and drying regions of the basins to follow the effective development, in time and space, of the expansion phenomenon of flood. The proposed procedure allows the grid's dynamic refinement to avoid coordinate transformation or the use of unstructured grids. The proposed method is simple and allows the thickening of the grid to accommodate the flooding phenomena on the floodplain and to calculate the velocity in the domain regions in which a higher space resolution is required. Therefore, flows running through structures such as weirs, gates, bridges or culverts can be simulated. In the paper two different case studies, approached with the proposed self‐adaptive calculation scheme, are discussed. The studies concern the analysis of the effects of structures, such as roads or embankments, on flooding phenomena in the Tiber and Tanaro basins respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation and comparison of potential flow and viscous fluid models in near trapping of narrow gaps

Numerical simulation results of fluid oscillation in narrow gaps subjected to incident water waves are presented in this paper. Both viscous fluid model and potential flow model with artificial viscous damping are employed to investigate wave motions in a single narrow gap separated by twin bodies and double narrow gaps formed by three identical bodies. Variation of wave heights in the narrow gaps with incident wave frequency is compared with experimental data available in literatures. The numerical results demonstrate that both the viscous fluid model and the potential flow model are able to predict the resonant frequency reasonably well. However the conventional potential flow model (without introducing the artificial viscous damping) significantly over-predicts the wave height in narrow gaps around the resonant frequency. The potential flow model with artificial damping predicts resonant wave heights well with a properly tuned (or calibrated) damping coefficient. The horizontal wave force on the individual body is also examined in this work. It is found that a properly calibrated damping coefficient is also very important for estimating the horizontal wave forces on the structures. A damping coefficient of 0.4 in the potential flow model is found to result in satisfactory predictions of wave heights in the narrow gaps and the horizontal hydrodynamic force for the cases considered in this work.

     

The dynamics of the floodwater and the damaged ship in waves

The interaction between the dynamics of the floodwater and the damaged ship in waves is investigated by applying an integrated method, which couples a seakeeping solver and a Navier-Stokes solver. To reveal the effects of the water flooding and the sloshing on the damaged ship behaviour, the motion of a Ro-Ro ferry in regular beam seas is simulated, including the ship under the intact condition with and without internal water and the ship under a damaged condition. It is found that the shift of the natural roll frequency of a damaged ship and the decrease of its roll response are mainly due to the water sloshing inside the compartment. The effect of the resonant sloshing leads to the presence of a ship's second peak response at higher frequencies and it is significantly reduced by the water flooding through the damaged opening. The influence of the flooding and the sloshing on the ship behaviour is small with a further increase of the wave frequency.     

Development of a hybrid particle-mesh method for simulating free-surface flows

In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method(FLIP)a formulation for the hybrid method is presented for incompressible multiphase flows involving large density jumps and wave generating boundaries.The performance of the method is assessed for a standing wave and for the generation and propagation of a solitary wave over a flat and a sloping bed.A comparison is made with results obtained with a well-established SPH package.The tests demonstrate that the method is a promising and attractive tool for simulating the nearshore propagation of waves.     

Roll Waves in High Gradient Channels

Abstract

Vedernikov's number quantifies the allowable capacity that can be sustained as a stable uniform flow in a high gradient channel, but it does not provide any guidance to the design of freeboard when roll waves exist. Without considering roll waves, high gradient channels are undersized using the uniform flow approach and fail to maintain design integrity. This paper presents revised design curves to determine the stability of flow in high gradient channels. When the design condition exceeds the limiting condition, roll waves should be considered in design. This study applies the model of moving hydraulic jump to simulate roll waves. It provides estimations for wave heights when Froude number is greater than 1.50.     

Identifying the natural flow regime and the relationship with riparian vegetation for two contrasting western Australian rivers

The natural flow regime and the relationship between flows and riparian vegetation are described for sites on both the Blackwood River in south‐western Australia and the Ord River in north‐western Australia. Analysis of long‐term flow data showed the historic mean monthly river discharge for the Blackwood River is strongly seasonal and highly predictable with generally low variability each month. The Ord River showed a strong seasonality of flows with about 92% of the (total) yearly flow occurring between December and March. Flow variability was very high (e.g. coefficient of variation >100% for all months) but highly predictable, with this mostly attributed to low but constant dry‐season flows. Water depth, duration of flood events and the number of flood events per year show a significant correlation with aspects of the riparian vegetation within experimental vegetation plots. Results highlight the strong relationship between floristics, life form structure and population dynamics with stream hydrology. On the Blackwood River, species richness and cover of shrubs reduced with increased duration and frequency of flooding, while cover of exotic species and annual herbs increased with increased flooding. Germination of tree seedlings was not influenced by flood regime but size class of tree species increased with flooding frequency. On the Ord River, species richness was not influenced by flooding regime. However, cover of perennial grasses increased with flooding frequency whilst cover of shrubs decreased. There was no relationship between flooding and seedling establishment whilst tree size class decreased with increased flooding. The methods described here can be used to compare the response of different components of the riparian vegetation to different fluvial regimes (e.g. because of impoundment and abstraction). This technique can be expanded for the management of riparian zones and planning rehabilitation programmes. It may also be useful for improving the ecological knowledge base for setting environmental flows in regulated systems. Copyright © 2001 John Wiley & Sons, Ltd.     

An experimental study of impact loading on deck of shore-connecting jetties exposed to oblique waves and current

The impact pressure from waves is an important issue to be considered in the design of coastal structures. In this paper, the waves acting on the deck of a shore-connecting jetty on a slope exposed to oblique waves and in the presence of current are examined based on laboratory experiments. The impact pressures are measured on a 1：50 scale model of a jetty head with down-standing beams and berthing members. The relations of the impact pressure with the incident wave angle and the current velocity are examined. It is shown that the impact pressure is sensitive to the wave angle and the current velocity. A computational model for the impact load on the deck of shore-connecting jetties exposed to oblique waves and current is developed.     

Numerical simulations of viscous flows around surface ship by level set method

Numerical simulations of viscous flows around surface ships by coupling the 3D incompressible RANS equations with level set method are presented in this paper. The finite difference method is used to discretize the RANS equations with turbulent model SST k − ω. The fully nonlinear boundary condition at the free surface is satisfied at each time step and the evolution of the free surface is achieved by using the level set method. The coupled solver is applied to a benchmark case of viscous flows around an advancing Wigley ship with various Froude numbers. The computational results are in excellent agreement with experimental data. The simulations reveal clearly the generation and evolution of bow and stern waves.