Geomorphic dam-break flows. Part I: conceptual model

Proceedings of the Institution of Civil Engineers - Water Management 163 Issue WM6

Solutions of a two-dimensional dam-break problem

Solutions of a two-dimensional dam break problem are presented for two tailwater/reservoir height ratios. The numerical scheme used is an extension of one previously given by the author [J. Hyd. Res. 26(3), 293–306 (1988)], and is based on numerical characteristic decomposition. Thus approximate solutions are obtained via linearised problems, and the method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations.

Numerical Experimental Comparison of Dam Break Flows with non-Newtonian Fluids

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Experimental investigation of dynamic pressure loads during dam break

The objectives of this work are to revisit the experimental measurements on dam break flow over a dry horizontal bed and to provide a detailed insight into the dynamics of the dam break wave impacting a vertical wall downstream the dam, with emphasis on the pressure loads. The measured data are statistically analyzed and critically discussed. As a result, an extensive set of data for validation of computational tools is provided.

Dam break wave of thixotropic fluid

Thixotropy is the characteristic of a fluid to form a gelled structure over time when it is not subjected to shearing, and to liquefy when agitated. Thixotropic fluids are commonly used in the construction industry (e.g., liquid concrete and drilling fluids), and related applications include some forms of mud flows and debris flows. This paper describes a basic study of dam break wave with thixotropic fluid. Theoretical considerations were developed based upon a kinematic wave approximation of the Saint-Venant equations down a prismatic sloping channel. A very simple thixotropic model, which predicts the basic theological trends of such fluids, was used. It describes the instantaneous state of fluid structure by a single parameter. The analytical solution of the basic flow motion and theology equations predicts three basic flow regimes depending upon the fluid properties and flow conditions, including the initial degree of jamming of the fluid (related to its time of restructuration at rest). These findings were successfully compared with systematic bentonite suspension experiments. The present work is the first theoretical analysis combining the basic principles of unsteady flow motion with a thixotropic fluid model and systematic laboratory experiments.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

溃坝流的光滑粒子法模拟

用光滑粒子法模拟几种情况的溃坝流动,对坝外无水和有水的情况进行二维模拟,对坝外有立柱的情况进行三维模拟.流动控制方程采用雷诺平均方程模拟溃坝流动的湍流效应,采用混合长度形式的涡粘模式对控制方程进行封闭,推导其相应的光滑粒子形式的方程.模拟结果表明,数值模拟的流动特征与实验结果符合得非常好,说明发展的光滑粒子法有效.

基于SPH方法的溃坝流动数值模拟

应用光滑粒子流体动力学方法对溃坝问题进行了数值模拟,在现有SPH方法的理论基础上对密度近似方程进行了重新初始化处理,分析了密度重新初始化对溃坝流动问题的影响,并对SPH数值模拟所得到的结果与试验结果以及移动粒子半隐式方法得到的计算结果进行了比较和分析.结果表明,对密度近似方程进行重新初始化保持了流场内的质量守恒,同时整个计算域内的压力分布更加规则,SPH方法数值模拟得到的结果与实验结果和MPS方法得到的结果非常吻合,验证了改进方案及所编程序的可靠性和计算结果的准确性

藏东南巨厚松散堆积体发育条件、成灾模式及防治对策——以川藏公路然乌至鲁朗段为例

The unique geologic, geomorphic and climatic conditions of southeast Tibet have made the region to develop the multi-style and frequently occurring geologic hazards, especially the collapses and landslides and debris flows along the section of Ranwu-Lulang in Sichuan-Tibet highway. However, most of those geologic hazards have close relationship with the loose accumulations. That is, the loose accumulations are the main carrier of most geologic hazards. Thereof, the huge-thick accumulations along the highway is regarded as the objective in the thesis to study the geologic background, hazarding model and mitigation methods comprehensively, based on the multi-disciplinary theories and former materials. First of all, in the paper, based on field engineering geologic investigations, the genetic type and the characteristics of spatiotemporal distribution of the huge-thick loose accumulations along the highway, have been analysized from the factors of regional geology and geomorphy and climate, as well as the coupling acting of those factors with inoculation and eruption of the loose accumulations geologic hazards. The huge-thick loose accumulations has complex genetic types and specific regulations of spatiotemporal distribution, closely controlled by the outer environment of the region. The accumulations are composed of earth and boulder, with disorder structure and poor sorting, specific forming environments and depositing conditions. And its physical and mechanic properties are greatly distinguished from rock and common earth inland. When Sichuan-Tibet highway was firstly constructed along the north bank of Purlung Tsangpo River, the huge-thick loose accumulations was cut into many high and steep slopes. Through the survey to the cut-slopes and systematic investigation to their failures, the combination of height and angle of the accumulations slope has been obtained. At the same time, the types of genetic structure of those cut-slopes are also analysized and concluded, as well as their failure models. It is studied in the paper that there are piaster, duality, multielement and complexity types in genetic structure, and rip-dump-repose, rip-shear-slip and weathering-flake types in failure models. Moreover, it is briefly introduced present engineering performance methods and techniques dealing with the deformation and failure of the accumulations cut-slope. It is also suggested that several new techniques of slope enforcement and the method of landslide and rockfall avoiding should be applied. The research of high and steep cut-slope along the highway has broadened the acknowledgement of the combination of cut-slope height and angle. Especially, the dissertation also has made the monographic studies about the geologic background and hazarding models and prevention methods of some classic but difficult accumulations geologic hazards. They are: (1) Research of the engineering geologic background of the 102 landslide group and key problems about the project of tunnel. The 102 landslide group is a famous accumulational one composed of glacial tills and glaciofuvial deposit. The tunnel project is a feasible and optional one which can solve the present plight of “sliding after just harnessing” in the 102 section. Based on the glacial geomorphy and its depositing character, distribution of seepage line, a few drillhole materials and some surveying data, the position of contact surface between gneiss and accumulations has been recognized, and the retreating velocities of three different time scales (short, medium and long term) have been approximately calculated, and the weathering thickness of gneiss has also been estimated in the paper. On the basis of above acknowledgement, new engineering geomechnic mode is established. Numerical analysis about the stability of the No.2 landslide is done by way of FLAC program, which supplies the conclusion that the landslide there develops periodically. Thereof, 4 projects of tunnel going through the landslide have been put forwards. Safety distance of the tunnel from clinohefron has been numerically analysized. (2) Research of the geologic setting and disaster model and hazard mitigation of sliding-sand-slope. From the geologic setting of talus cone, it is indicated that the sliding-sand-slope is the process of the re-transportation and re-deposit of sand under the gravity action and from the talus cone. It is the failure of the talus cone essentially. The layering structure of the sliding-sand-slope is discovered. The models of movement and failure of the sliding-sand-slope has been put forwards. The technique, “abamurus+grass-bush fence+degradable culture pan”, is suggested to enforcement and green the sliding-sand-slope. (3) Characteristics and hazarding model and disaster mitigation of debris flow. The sources of solid material of three oversize debris flows have been analysized. It is found that a large amount of moraine existing in the glacial valley and large landslide dam-break are the two important features for oversize debris flow to be taken place. The disaster models of oversize and common debris flows have been generalized respectively. The former model better interpret the event of the Yigong super-large landslide-dam breaking. The features of common debris flow along the highway section, scouring and silting and burying and impacting, are formulated carefully. It is suggested that check dam is a better engineering structure to prevent valley from steeply scouring by debris flow. Moreover, the function of check dam in enforcing the slope is numerically calculated by FLAC program. (4) Songzong ancient ice-dammed lake and its slope stability. The lacustrine profile in Songzong landslide, more than 88 meters thick, is carefully described and measured. The Optical Simulated Luminescence (OSL) ages in the bottom and top of the silty clay layer are 22.5±3.3 kaB.P., 16.1±1.7 kaB.P., respectively. It is indicated by the ages that the lacustrine deposits formed during the Last Glacial Maximum ranging from 25ka B.P. to 15ka B.P. The special characteristics of the lacustrine sediment and the ancient lake line in Songzong basin indicated that the lacustrine sediment is related to the blocking of the Purlung Tsangpo River by the glacier in Last Glacial Maximum from Dongqu valley. The characteristics of the lacustrine profile also indicate that the Songzong ice-dammed lake might run through the Last Glacial Maximum. Two dimensional numerical modeling and analysis are done to simulate the slope stability under the conditions of nature and earthquake by FLAC program. The factor of safety of the lacusrtine slope is 1.04, but it will take place horizontal flow under earthquake activity due to the liquefaction of the 18.33 m silt layer. The realign to prevent the road from landslide is suggested.

Approximate Riemann solutions of the two-dimensional shallow-water equations

A finite-difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow-water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gasdynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. An extension to the two-dimensional equations with source terms, is included. The scheme is applied to a dam-break problem with cylindrical symmetry.

Efficient shock capturing for isentropic flows using arithmetic averaging

A shock capturing scheme is presented for the equations of isentropic flow based on upwind differencing applied to a locally linearized set of Riemann problems. This includes the two-dimensional shallow water equations using the familiar gas dynamics analogy. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency, leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver where the computational expense can be prohibitive. The scheme is applied to a two-dimensional dam-break problem and the approximate solution compares well with those given by other authors.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.