SPH simulation of free surface flow over a sharp-crested weir

In this paper the numerical simulation of a free surface flow over a sharp-crested weir is presented. Since in this case the usual shallow water assumptions are not satisfied, we propose to solve the problem using the full weakly compressible Navier–Stokes equations with the Tait equation of state for water. The numerical method used consists of the new meshless Smooth Particle Hydrodynamics (SPH) formulation proposed by Ferrari et al. (2009) [8], that accurately tracks the free surface profile and provides monotone pressure fields. Thus, the unsteady evolution of the complex moving material interface (free surface) can been properly solved. The simulations involving about half a million of fluid particles have been run in parallel on two of the most powerful High Performance Computing (HPC) facilities in Europe. The validation of the results has been carried out analysing the pressure field and comparing the free surface profiles obtained with the SPH scheme with experimental measurements available in literature [18]. A very good quantitative agreement has been obtained.     

A new weir surface repair layer with shock absorbing cushion design

A shock absorbing cushion has never been introduced into any traditional weir surface repair layer design.However,shocks induced by high discharge with heavy sediment can easily produce brittle fracture and peeling over a weir surface repair layer as it is impacted by floods accompanied by particles of different sizes.In this study,transcending traditional designs,the authors developed a composite unit designed with a shock absorbing cushion that has performed well during field tests,proving that the weir body can be effectively protected even if the composite units are directly laid on a severely uneven weir surface repair layer.     

Numerical model of local scour considering the unsteady sediment inflow and sediment sorting: Application to scour upstream of slit weir

The current study proposes a novel framework for the numerical model for estimating the temporal scour considering unsteady sediment inflow and the sediment sorting process. The framework was applied to local scour upstream of a slit weir. The scour model is based on an ordinary nonlinear differential equation derived from sediment continuity and scour rate equations. A one-dimensional(1-D)Boussinesq-type model coupled with nonequilibrium sediment transport was incorporated in the scour model to...     

Numerical modeling of local scour at a submerged weir with a downstream slope using a coupled moving-mesh and maskedelement approach

This paper presents a model for local scour at submerged weirs with downstream slopes that uses a coupled moving-mesh and masked-element approach.In the developed model,the fluid-sediment interface is tracked using a moving-mesh technique,and the effects of the structure on the hydrodynamics and bed morphology are resolved using a masked-element technique.Compared to traditional sediment scour models,based on the moving-mesh technique,the present model has the advantage of allowing for a simpler setup of the computational grids and a larger-amplitude deformation.Laboratory experiments on local scour at a submerged weir with a downstream slope were conducted,which provided bed profiles at different time instants.The results obtained by the present model are compared to the experimental data.The comparisons demonstrate the performance of the model in satisfactorily predicting local scour at a submerged weir with a downstream slope.The model was further modified and employed to carry out additional computations to investigate the influence of various parameters and sub-models.     

考虑流固耦合效应的辽宁葠窝水库溢流坝段抗震性能分析

针对辽宁葠窝水库混凝土重力坝抗震问题,采用耦合的拉格朗日-欧拉有限元分析技术,建立了可考虑库水-坝体-基岩动力耦合效应的典型溢流坝段抗震分析数值模型。模型中,采用等效一致粘弹性边界模拟基岩的人工截断边界;采用混凝土弥散裂缝本构模型模拟混凝土的动力特性。根据烈度与地震动之间的关系,确定了水库坝体抗震设计的输入加速度峰值。据此,分析了在不同季节水位变化条件下坝体地震反应的基本特性。研究表明：完好的辽宁葠窝水库混凝土重力坝溢流坝段能满足8度的抗震设防烈度要求。地震下溢流坝段峰值位移出现在胖坝和瘦坝的坝顶迎水面位置处,胖坝的动位移较瘦坝动位移大。胖坝在闸墩与溢流堰交接处出现了拉应力最大值。有库水条件下,瘦坝峰值拉应力出现在坝趾处,无库水条件下,瘦坝最大拉应力出现在溢流堰与闸墩交接处。     

On the comparison of laboratory electromagnetic analogue model measurements and finite-difference numerical calculations

This work compares experimental analogue model measurements and finite-difference numerical calculations of the electric and magnetic fields for the cases of induction in a highly conducting plate and wedge embedded in a poorly conducting host earth. The results indicate very good agreement between the experimental and theoretical results and confirm the validity of the analogue and numerical methods for studying complex induction problems.     

Two-dimensional effective stress analysis of liquefaction of irregular ground including soil–structure interaction

An effective stress method is presented for analysis of seismic response and liquefaction of irregular ground including soil–structure interaction, based on an implicit–explicit finite element method. A pore water pressure is computed with iteration from the total stress considering an undrained condition. The simulated pore water pressure is in reasonably good agreement with the experimental data. The proposed method of analysis is compared with other well-known methods for a one-dimensional model, which is in good agreement. The present effective stress method is also applied to liquefaction problems involving a two-dimensional soil–structure model. The structure is modelled by not only a rigid model but also as a multi-degree-of-freedom system with bi-linear springs. The numerical results are considered to be significant from the viewpoint of earthquake engineering.     

Estimation of debris flow discharge coefficient considering sediment concentration

A sabo dam has a purpose to block the path of debris flow. However, when overflow occurs, a sabo dam works as a weir, a vertical obstruction, where the fluid must flow over. Many empirical formulas and discharge coefficients for weirs relating flow depth to discharge have been proposed to calculate overflow discharges. However, only a few studies about overflow discharge coefficients are available in the case of debris flow. In this paper, experiments and numerical simulations were done to estimate debris flow discharge coefficients by considering the sediment concentration. In the numerical simulation, a complete overflow equation and a free overfall equation were implemented to calculate debris overflow discharges at a sabo dam. To determine the discharge coefficients for each equation, single factor regression analysis was used. Laboratory experiments were done to calibrate and to compare with the simulation. Study results showed that the discharge coefficients increase as the sediment concentration increases. This finding suggests debris flow discharge coefficients are derived to calculate the debris overflow discharges at a sabo dam.     

水位仪在水堰流量观测中的应用——以新疆乌苏33号井为例

详细介绍了用于测量水流量水堰的种类、结构、测量原理和方法,根据数字化水位仪投入式压力传感器的特点,结合地震观测的实际情况,主要介绍了利用水位仪在三角水堰中测量水流量的方法、计算公式及测量范围,并通过应用实例对比分析了人工与仪器的测量结果,为开展此项工作的地震台站提供了设计制造水堰的具体参数和观测方法。     

A comparison between artificial neural network algorithms and empirical equations applied to submerged weir scour evolution prediction

Estimating the time evolution of a local scour hole downstream of submerged weirs can help determine the maximum scour depth and length and is essential to designing submerged weir foundations.In the current study, artificial neural networks with a backpropagation learning algorithm were used to estimate the temporal variation of scour profiles downstream of submerged weirs under clear water conditions. Physical factors, such as the flow condition, weir size, and sediment characteristics, are ge...     

地球物理问题中电磁感应的实验室模拟模型研究

在地球物理问题中,有两种方法通常用于解释地磁和地电异常,这就是数字模型技术和实验室模拟模型法。实验室模拟模型对于分析不易求得数学解的问题非常有用,并且已经广泛用于研究复杂的地球物理问题。本文评述了国外特别是加拿大开展模拟模型研究的情况;描述了平面波、线电流、磁偶极等各种类型的场源模型;给出了该法在研究海岸效应、岛陆通道电流、海浪电磁效应和各向异性导体等方面的应用。模拟模型测量和数字计算的比较结果显示了非常好的一致性,这就进一步证实了模拟模型法对于研究复杂二维和三维感应问题的可靠性。该项研究工作对于了解天然电磁场源的性质、进行地球物理勘探以及研究地球地壳和地幔的电性结构都十分有益。     

Numerical analysis of seismic wave amplification in Nice (France) and comparisons with experiments

The analysis of site effects is very important since the amplification of seismic motion in some specific areas can be very strong. In this paper, the site considered is located in the centre of Nice on the French Riviera. Site effects are investigated considering a numerical approach (Boundary Element Method) and are compared with experimental results. The experimental results are obtained thanks to real earthquakes (weak motion) and microtremor measurements. The investigation of seismic site effects through numerical approaches is interesting because it shows the dependency of the amplification level on such parameters as wave velocity in surface soil layers, velocity contrast with deep layers, seismic wave type, incidence and damping.In this specific area of Nice, experimental measurements obtained for weak motion lead to strong site effects. A one-dimensional (1D) analytical analysis of amplification does not give a satisfactory estimation of the maximum reached levels. A boundary element model is then proposed considering different wave types (SH, P, SV) as the seismic loading. The alluvial basin is successively assumed as an isotropic linear elastic medium and an isotropic linear viscoelastic solid with Zener type behaviour (standard solid). The influence of frequency and incidence is analysed. The thickness of the surface layer, its mechanical properties, its general shape as well as the seismic wave type involved have a great influence on the maximum amplification and the frequency for which it occurs. For real earthquakes, the numerical results are in very good agreement with experimental measurements for each motion component. The boundary element method leads to amplification values very close to the actual ones and much larger than those obtained in the 1D case. Two-dimensional basin effects are then very strong and are well reproduced numerically.     

Numerical simulation of laboratory experiments in detention pond routing with long rainfall duration

This paper establishes a numerical detention pond volume model based on the hydrological continuity equation and the Runge-Kutta numerical method. Experiments for the conditions of both steady and unsteady flow have been used to verify the model. In unsteady flow cases, the outflow hydrograph by numerical simulation are fairly consistent with experimental value. Both experimental and numerical results indicate that wider rectangular sharp-crested weirs or larger rectangular slot tend to induce greater outflow discharges, which undesirably cut down the detention volume. Experiments show that the necessary detention volume of rectangular slot is smaller than that of the rectangular sharp-crested weir for a constant flood peak reduction. That is, the rectangular slot is the recommended outflow device when flood peak reduction is the design criteria. The study also shows that necessary detention volume of the short rainfall duration is less than that of the long rainfall duration under constant allowable maximum discharge.     

A numerical method for simulating non-Newtonian fluid flow and displacement in porous media

Flow and displacement of non-Newtonian fluids in porous media occurs in many subsurface systems, related to underground natural resource recovery and storage projects, as well as environmental remediation schemes. A thorough understanding of non-Newtonian fluid flow through porous media is of fundamental importance in these engineering applications. Considerable progress has been made in our understanding of single-phase porous flow behavior of non-Newtonian fluids through many quantitative and experimental studies over the past few decades. However, very little research can be found in the literature regarding multi-phase non-Newtonian fluid flow or numerical modeling approaches for such analyses.For non-Newtonian fluid flow through porous media, the governing equations become nonlinear, even under single-phase flow conditions, because effective viscosity for the non-Newtonian fluid is a highly nonlinear function of the shear rate, or the pore velocity. The solution for such problems can in general only be obtained by numerical methods.We have developed a three-dimensional, fully implicit, integral finite difference simulator for single- and multi-phase flow of non-Newtonian fluids in porous/fractured media. The methodology, architecture and numerical scheme of the model are based on a general multi-phase, multi-component fluid and heat flow simulator — TOUGH2. Several rheological models for power-law and Bingham non-Newtonian fluids have been incorporated into the model. In addition, the model predictions on single- and multi-phase flow of the power-law and Bingham fluids have been verified against the analytical solutions available for these problems, and in all the cases the numerical simulations are in good agreement with the analytical solutions. In this presentation, we will discuss the numerical scheme used in the treatment of non-Newtonian properties, and several benchmark problems for model verification.In an effort to demonstrate the three-dimensional modeling capability of the model, a three-dimensional, two-phase flow example is also presented to examine the model results using laboratory and simulation results existing for the three-dimensional problem with Newtonian fluid flow.     

Modelling transport dynamics of contaminated sediments in the headwater of a hydropower plant at the Upper Rhine River

Sediments are an essential habitat compartment in rivers, which is a subject to dynamic transport processes. In many rivers, the fine deposited sediments are contaminated with heavy metals and organic compounds. Contaminated deposits are considered as potential hot spots because of the risk of the mobilization under erosive hydraulic conditions. Numerical models for particulate contaminant transport are then necessary and can be applied to estimate and predict the potential impact of mobilized contaminants as an important contribution to sediment management. This paper focuses on the quantification of the amount of contaminated sediments resuspended during the extreme flood event in 1999 and the prediction of deposition one year after the flood event. To assess such erosive flood event, a 2D numerical transport model was developed to analyse the dynamics of erosion and sedimentation processes in the headwater of a cross dam at the Upper Rhine River. The dam consists of a weir, a hydropower plant, and a navigation lock. As the weir is operating only for flood management, a huge amount of sediment highly contaminated with the hexachlorobenzene (HCB) was deposited in the weir zone. Therefore, numerical simulations were performed to determine the spatial and temporal distribution of deposited contaminated sediments as depending on the river discharge and its distribution to the hydraulic structures. The numerical investigation presented here is taken as a retrospective analysis of the contaminated sediment dynamics in the headwater to improve future sediment management.     

Analytical approximations for flow in compressible,saturated, one-dimensional porous media

A nonlinear model for single-phase fluid flow in slightly compressible porous media is presented and solved approximately. The model assumes state equations for density, porosity, viscosity and permeability that are exponential functions of the fluid (either gas or liquid) pressure. The governing equation is transformed into a nonlinear diffusion equation. It is solved for a semi-infinite domain for either constant pressure or constant flux boundary conditions at the surface. The solutions obtained, although approximate, are extremely accurate as demonstrated by comparisons with numerical results. Predictions for the surface pressure resulting from a constant flux into a porous medium are compared with published experimental data.     

Modelling detention basins measured from high‐resolution light detection and ranging data

The construction of stormwater detention basins is a best management practice to effectively control floods, to provide additional surface storage for excess floodwater and to compensate for the adverse effects of urban development. Traditional field‐based levelling survey methods are very time consuming and subject to human‐induced arbitrariness and error. This article presents an approach to modelling detention basins measured from light detection and ranging remote sensing data. A case study is illustrated by using the White Oak Bayou watershed of Harris County, Texas. The storage–stage curve obtained from the volumetric analysis is used in a modified detention basins routing model, which was developed by adding the weir structure control to the traditional hydrologic reservoir routing equations. The model simulation showed that the peak flow of the synthetic 100‐year reoccurrence event was effectively reduced and delayed by the detention basins. The comparison with the simulation results from the traditional reservoir routing model suggested that previous studies using the reservoir routing model were likely to underestimate the flood reduction effect of detention basins. The sensitivity analysis of the parameters showed that the detention basin design and evaluation should pay more attention on the weir height and river channel's roughness. Copyright © 2011 John Wiley & Sons, Ltd.     

Tide‐Induced Air Pressure Fluctuations in a Coastal Unsaturated Zone: Effects of Thin Low‐Permeability Pavements

This paper presents an analytical case study to explore one‐dimensional subsurface air pressure variation in a coastal three‐layered unsaturated zone. The upper layer is thin and much less permeable than the middle layer, and water table is located in the very permeable lower layer. An analytical solution was derived to describe the air pressure variation caused by tide‐induced water table fluctuations. We revisited the case study at Hong Kong International Airport conducted by Jiao and Li (2004) who used a two‐dimensional numerical model. The analytical prediction using the parameter values equivalent to the two‐dimensional numerical model agreed very well with the observed air pressure, indicating the validity and applicability of our one‐dimensional model in approximating the actual situation in this coastal zone with adequate accuracy. The analysis revealed that the asphalt pavement played an important role in causing air pressure fluctuations below it. Abnormally high air pressure can be caused beneath the surface pavement when the air permeability decreases due to rainfall infiltration, which may lead to heaving problems during rising tides.     

Investigation of pressure variations over stepped spillways using smooth particle hydrodynamics

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.     

Practical implementation of the fractional flow approach to multi-phase flow simulation

Fractional flow formulations of the multi-phase flow equations exhibit several attractive attributes for numerical simulations. The governing equations are a saturation equation having an advection diffusion form, for which characteristic methods are suited, and a global pressure equation whose form is elliptic. The fractional flow approach to the governing equations is compared with other approaches and the implication of equation form for numerical methods discussed. The fractional flow equations are solved with a modified method of characteristics for the saturation equation and a finite element method for the pressure equation. An iterative algorithm for determination of the general boundary conditions is implemented. Comparisons are made with a numerical method based on the two-pressure formulation of the governing equations. While the fractional flow approach is attractive for model problems, the performance of numerical methods based on these equations is relatively poor when the method is applied to general boundary conditions. We expect similar difficulties with the fractional flow approach for more general problems involving heterogenous material properties and multiple spatial dimensions.