共查询到20条相似文献，搜索用时 31 毫秒
1.
Smoothed particle hydrodynamics(SPH) is a mesh-free adaptive Lagrangian particle method with attractive features for dealing with the free surface flow.This paper applies the SPH method to simulate the large-amplitude lateral sloshing both with and without a floating body,and the vertical parametrically-excited sloshing in a two-dimensional tank.The numerical results show that the SPH approach has an obvious advantage over conventional mesh-based methods in handling nonlinear sloshing problems such as violent fluid-solid interaction,and flow separation and wave-breaking on the free fluid surface.The SPH method provides a new alternative and an effective way to solve these special strong nonlinear sloshing problems. 相似文献
2.
A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained. 相似文献
3.
In this paper the smoothed particle hydrodynamics, (SPH), technique is used to investigate the pressure distribution on steps located in the non-aerated flow region of a stepped spillway for different discharges typical of skimming flow conditions. The open source code 2D SPHysics has been employed after being validated against the laboratory model studies of flow over broad crested weirs and flow over stepped spillways. The numerical results, in terms of the water surface and velocity profiles at different sections, are in good agreement with the corresponding experimental results. The code is then applied to determine the pressure distribution on the vertical and horizontal step faces. Also, the aspects of the pressure pattern are described and the positions/magnitudes of the maximum and minimum pressure values are presented. 相似文献
4.
A two-phase numerical model using Smoothed Particle Hydrodynamics (SPH) is applied to two-phase liquid-sediments flows. The absence of a mesh in SPH is ideal for interfacial and highly non-linear flows with changing fragmentation of the interface, mixing and resuspension. The rheology of sediment induced under rapid flows undergoes several states which are only partially described by previous research in SPH. This paper attempts to bridge the gap between the geotechnics, non-Newtonian and Newtonian flows by proposing a model that combines the yielding, shear and suspension layer which are needed to predict accurately the global erosion phenomena, from a hydrodynamics prospective. The numerical SPH scheme is based on the explicit treatment of both phases using Newtonian and the non-Newtonian Bingham-type Herschel-Bulkley-Papanastasiou constitutive model. This is supplemented by the Drucker-Prager yield criterion to predict the onset of yielding of the sediment surface and a concentration suspension model. The multi-phase model has been compared with experimental and 2-D reference numerical models for scour following a dry-bed dam break yielding satisfactory results and improvements over well-known SPH multi-phase models. With 3-D simulations requiring a large number of particles, the code is accelerated with a graphics processing unit (GPU) in the open-source DualSPHysics code. The implementation and optimisation of the code achieved a speed up of x58 over an optimised single thread serial code. A 3-D dam break over a non-cohesive erodible bed simulation with over 4 million particles yields close agreement with experimental scour and water surface profiles. 相似文献
5.
A 3D numerical modeling of the wave generated by the Vajont slide, one of the most destructive ever occurred, is presented in this paper. A meshless Lagrangian Smoothed Particle Hydrodynamics (SPH) technique was adopted to simulate the highly fragmented violent flow generated by the falling slide in the artificial reservoir. The speed-up achievable via General Purpose Graphic Processing Units (GP-GPU) allowed to adopt the adequate resolution to describe the phenomenon. The comparison with the data available in literature showed that the results of the numerical simulation reproduce satisfactorily the maximum run-up, also the water surface elevation in the residual lake after the event.Moreover, the 3D velocity field of the flow during the event and the discharge hydrograph which overtopped the dam, were obtained. 相似文献
6.
This paper presents the applications of a newly developed free surface flow model to the practical, while challenging overflow problems for weirs. Since the model takes advantage of the strengths of both the level set and volume of fluid methods and solves the Navier-Stokes equations on an unstructured mesh, it is capable of resolving the time evolution of very complex vortical motions, air entrainment and pressure variations due to violent deformations following overflow of the weir crest. In the present study, two different types of vertical weir, namely broad-crested and sharp-crested, are considered for validation purposes. The calculated overflow parameters such as pressure head distributions, velocity distributions, and water surface profiles are compared against experimental data as well as numerical results available in literature. A very good quantitative agreement has been obtained. The numerical model, thus, offers a good alternative to traditional experimental methods in the study of weir problems. 相似文献
7.
A three‐dimensional, numerical free‐surface‐flow simulation model is developed to investigate hydrodynamics of a lake and the predictive capabilities of the numerical model are validated by comparing them with field results obtained from Lake Naini, Uttar Pradesh, using environmental isotope (δ18O, δD and 3H) techniques. This has been carried out in order to understand the hydrodynamics of a lake (i.e. circulations, mixing and transport of water within the lake). Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
8.
9.
This paper describes the application of the Smoothed Particle Hydrodynamics（SPH） method for modeling two dimensional waves caused by dam break over a movable bed in two dimensions.The two phase SPH method is developed to solve the Navier-Stokes equations.Both fluid and sediment phases are described by particles as weakly compressible fluids and the incompressibility is achieved by the equation of state.The sediment phase is modeled as a non-Newtonian fluid using three alternative approaches of artificial viscosity and Bingham Model.In this paper,the new formulations for two-phase flows are proposed.The numerical results obtained from the developed SPH model show acceptable accuracy with comparison to experimental data. 相似文献
10.
Flow on fracture surfaces has been identified by many authors as an important flow process in unsaturated fractured rock formations. Given the complexity of flow dynamics on such small scales, robust numerical methods have to be employed in order to capture the highly dynamic interfaces and flow intermittency. In this work we use a three-dimensional multiphase Smoothed Particle Hydrodynamics (SPH) model to simulate surface tension dominated flow on smooth fracture surfaces. We model droplet and film flow over a wide range of contact angles and Reynolds numbers encountered in such flows on rock surfaces. We validate our model via comparison with existing empirical and semi-analytical solutions for droplet flow. We use the SPH model to investigate the occurrence of adsorbed trailing films left behind droplets under various flow conditions and its importance for the flow dynamics when films and droplets coexist. It is shown that flow velocities are higher on prewetted surfaces covered by a thin film which is qualitatively attributed to the enhanced dynamic wetting and dewetting at the trailing and advancing contact lines. Finally, we demonstrate that the SPH model can be used to study flow on rough surfaces. 相似文献
11.
《Advances in water resources》2007,30(4):730-741
In this work the numerical integration of 1D shallow water equations (SWE) over movable bed is performed using a well-balanced central weighted essentially non-oscillatory (CWENO) scheme, fourth-order accurate in space and in time. Time accuracy is obtained following a Runge–Kutta (RK) procedure, coupled with its natural continuous extension (NCE). Spatial accuracy is obtained using WENO reconstructions of conservative variables and of flux and bed derivatives. An original treatment for bed slope source term, which maintains the established order of accuracy and satisfies the property of exactly preserving the quiescent flow (C-property), is introduced in the scheme. This treatment consists of two procedures. The former involves the evaluation of the point-values of the flux derivative, considered as a whole with the bed slope source term. The latter involves the spatial integration of the source term, analytically manipulated to take advantage from the expected regularity of the free surface elevation. The high accuracy of the scheme allows to obtain good results using coarse grids, with consequent gain in terms of computational effort. The well-balancing of the scheme allows to reproduce small perturbations of the free surface and of the bottom otherwise of the same order of magnitude of the numerical errors induced by the non-balancing. The accuracy, the well-balancing and the good resolution of the model in reproducing free surface flow over movable bed are tested over analytical solutions and over numerical results available in literature. 相似文献
12.
Most numerical schemes applied to solve the advection–diffusion equation are affected by numerical diffusion. Moreover, unphysical results, such as oscillations and negative concentrations, may emerge when an anisotropic dispersion tensor is used, which induces even more severe errors in the solution of multispecies reactive transport. To cope with this long standing problem we propose a modified version of the standard Smoothed Particle Hydrodynamics (SPH) method based on a Moving-Least-Squares-Weighted-Essentially-Non-Oscillatory (MLS-WENO) reconstruction of concentrations. This scheme formulation (called MWSPH) approximates the diffusive fluxes with a Rusanov-type Riemann solver based on high order WENO scheme. We compare the standard SPH with the MWSPH for different a few test cases, considering both homogeneous and heterogeneous flow fields and different anisotropic ratios of the dispersion tensor. We show that, MWSPH is stable and accurate and that it reduces the occurrence of negative concentrations compared to standard SPH. When negative concentrations are observed, their absolute values are several orders of magnitude smaller compared to standard SPH. In addition, MWSPH limits spurious oscillations in the numerical solution more effectively than classical SPH. Convergence analysis shows that MWSPH is computationally more demanding than SPH, but with the payoff a more accurate solution, which in addition is less sensitive to particles position. The latter property simplifies the time consuming and often user dependent procedure to define the initial dislocation of the particles. 相似文献
13.
《国际泥沙研究》2021,36(6):711-722
The current study tries a new approach to simulating interactions between waves and seagrass through Smoothed Particle Hydrodynamics (SPH). In this model, the plants are defined as a solid that respects Hooke's law, and are assumed to have direct interaction with the fluid. Given the characteristics of the SPH method, especially in terms of computational time, the dimensions of the simulations were limited. The first goal of the current study was to optimize the approach to avoid reaching certain limits such as the rupture of the simulated plant. Plant movements under waves and/or currents have been studied by several authors in various in-situ, physical, and numerical experiments concerning various vegetation species, thus proving that plant movements can be successfully reproduced by SPH 2D/3D. Manning's roughness coefficient, n, was calculated to confirm that the results were in accordance with what had been measured in flume studies. Even though there is still room for improvement, it is shown that this method can be used to estimate Manning's coefficient for coastal vegetation (seagrass and saltmarsh vegetation) and to greatly improve the modeling and forecasting of coastal erosion and storm surge risks by including the effects of vegetation in integrated models. 相似文献
14.
During the last four decades, several numerical formulations and specialized software have been developed in response to studies about dam break (DB) wave propagation and its hydraulic and environmental impacts on downstream hydraulic structures and valleys. These methods cannot, however, be used to predict wave propagation within partially covered channels or reservoirs located upstream of hydraulic structures. In fact, such problems require the modelling of the complex transition from a free surface flow into a pressurized one. Because rivers or channels partially covered with ice sheets are typical examples commonly met in winter in northern climates, it is vitally important to assess ice-cover effects on the DB wave propagation and develop appropriate tools to predict resulting hydrodynamic loads on hydraulic structures downstream. This paper proposes an original numerical formulation to model wave propagation and hydrodynamic pressure in partially covered channels. The proposed formulation uses one-dimensional St. Venant equations to simulate open-water flow and water hammer equations to simulate pressure flow within the partially covered channel. To illustrate the use of the hydrodynamic pressures obtained, a case study is presented where a channel cover and a dam located downstream are modelled using finite elements to investigate their dynamic structural response. 相似文献
15.
Arabinda Mukhopadhyay 《Pure and Applied Geophysics》1966,65(1):12-28
Summary The object of this paper is to study the stresses produced in a porous liquid-saturated solid of the type considered byBiot [1]2) by a normal load moving along the surface of such a solid. Two cases have been considered. In the first case, the porous solid is taken to be semi-infinite, while in the second case we consider a layer of the liquid-saturated porous solid resting on a smooth rigid foundation. Expressions have been obtained for the stresses at any point in the medium due to a normal load moving along the free surface of the porous solid. Detailed numerical calculations have been carried out for liquid-saturated sandstone, using the data given byFatt [9] for such a medium. These calculations show that the normal stress zz in the solid falls off approximately as 1/z with increase in the depthz below the free surface, and that the stresses in the solid and in the liquid fall off rapidly as we move away from the point directly below the moving load. 相似文献
16.
This paper proposes a coupled fluid layer–foundation–poroelastic half-space vibration model to study how still water affects foundations operating underwater. As an example, we consider the problem of the vertical vibration of a rigid disk on a poroelastic half-space covered by a fluid layer having a finite depth. The solution of the disk vibration problem is obtained using the boundary conditions at the free surface of the fluid layer and the boundary conditions at the fluid layer–poroelastic medium interface. The solution is expressed in terms of dual integral equations that are converted into Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the vertical dynamic impedance coefficient are examined based on different water depths, poroelastic materials, disk permeabilities and frequencies of excitation. Based on the numerical results, it is proposed that the hydrodynamic pressure caused by the foundation vibration is the intrinsic reason that the existence of a fluid layer has such a great effect on the dynamic characteristics of the foundation. In many cases, the hydrodynamic pressure caused by the foundation vibration cannot be ignored when designing dynamic underwater foundations. These results are helpful in understanding the dynamic response of foundations under still water without water waves, such as foundations in pools, lakes and reservoirs. 相似文献
17.
P.W. France 《Advances in water resources》1977,1(1):49-51
This paper describes an investigation into the motion of the water surface in a simple model surge tank, and the relevant factors governing its behaviour. The oscillation of the free water surface is an unsteady flow problem, which is amenable to a mathematical solution using a finite difference step-by-step integration procedure. For comparison, two such methods are presented: (i) a simple initial value method and (ii) a predictor-corrector technique. Computer programs have been developed linked to a graph plotter to give a visual presentation of the numerical solutions together with the experimental results of the damped oscillation of the water surface in the model surge tank. 相似文献
18.
Details are given of the development and application of a numerical model for predicting free-surface flows in estuarine and coastal basins using the finite volume method. Both second- and third-order accurate and oscillation free explicit numerical schemes have been used to solve the shallow water equations. The model deploys an unstructured triangular mesh and incorporates two types of mesh layouts, namely the ‘cell centred’ and ‘mesh vertex’ layouts, and provides a powerful mesh generator in which a user can adjust the mesh-size distribution interactively to create a desirable mesh. The quality of mesh has been shown to have a major impact on the overall performance of the numerical model.The model has been applied to simulate two-dimensional dam break flows for which transient water level distributions measured within a laboratory flume were available. In total 12 model runs were undertaken to test the model for various flow conditions. These conditions include: (1) different bed slopes (ranging from zero to 0.8%), (2) different upstream and downstream water level conditions, and (3) initially wet and dry bed conditions, downstream of the dam. Detailed comparisons have been made between model predicted and measured water levels and good agreement achieved between both sets of results. The model was then used to predict water level and velocity distributions in a real estuary, i.e. the Ribble Estuary, where the bed level varies rapidly at certain locations. In order to model the whole estuary, a 1-D numerical model has also been used to model the upper part of the estuary and this model was linked dynamically to the 2-D model. Findings from this application are given in detail. 相似文献
19.
Nils Reidar B. Olsen 《地球表面变化过程与地形》2017,42(14):2393-2401
A numerical model is presented that compute the geometrical dimensions and movement of downstream migrating antidunes. The model solves the Navier–Stokes equations together with the k‐epsilon turbulence model to find the water flow field over the bedforms. A two‐dimensional width‐averaged grid is used. The bed elevation changes are computed by solving the convection–diffusion equation for suspended sediments and bedload, together with the Engelund–Hansen sediment transport formula. The free surface is computed with an algorithm based on water continuity in the surface cells. Non‐orthogonal adaptive grids were used, moving vertically with the computed location of the bed and the free water surface. The numerical model was tested on data from a physical model study where regular downstream migrating antidunes had been observed. The numerical model started out with a flat bed and the trains of antidunes formed over time. Many of the physical processes observed in earlier studies were replicated by the numerical model. Four dune parameters were computed in the current tests: The antidune wavelength, height and celerity, together with the average water depth. The antidune wavelengths were best predicted with an accuracy of 3 to 8% compared with the measurements. The antidune heights were computed with a deviation of 11 to 25% compared with an empirical formula. The water depths over the antidunes were predicted with an accuracy of 3 to 9% related to the measured values. The average antidune celerity was the parameter with largest deviation: For the coarsest grid it was overpredicted with 37%. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
20.
1920年海原8.5级大地震诱发的石碑黄土塬滑坡,因其规模大、坡度小、滑距远的特点成为业界关注焦点,目前对该滑坡的物理力学过程仍无统一认识。基于此,文章通过理论分析和数值计算方法研究该滑坡初始状态、地震动力响应和流滑发展过程。为分析震前斜坡初始状态,建立考虑斜坡表面拉裂缝中侧向水压力作用的力学计算模型。结果表明,拉裂缝中侧向水压力削弱了斜坡整体稳定状态;为研究地震时斜坡动力响应过程,建立数值计算模型,获得地震时斜坡饱和土层的孔隙水压力比和斜坡位移的变化特性;为解释远程滑坡,将液化后土层等效为流体,借鉴泥流滑距估算理论,求算石碑塬滑坡的滑距,计算结果与前人现场考察结果吻合。 相似文献