Tangjiashan landslide is a typical high-speed consequent landslide of medium-steep dip angle. This landslide triggered by earthquake took place in about semi-minute. The relative sliding displacement is 900 meters, so...Tangjiashan landslide is a typical high-speed consequent landslide of medium-steep dip angle. This landslide triggered by earthquake took place in about semi-minute. The relative sliding displacement is 900 meters, so average sliding speed is about 30 meters per second. The longitudinal length of barrier dam which is formed by high-speed landslide along river is 803.4 meters; and maximum width crossing river is 611.8 meters. And its volume is estimated about 20.37 million steres. Through detailed geological investigation of the barrier dam, together with early geological information before earthquake, geological structures of the barrier dam and its stability of upstream and downstream slopes are studied when water level reaches different elevations in condition of continual after shocks with seismic intensity of 7 or 8 Richter scale. On this basis, dam-breaking mode of barrier dam is discussed deeply. Thereby, analytic results provide significant guidance and advices to front headquarters of Tangjiashan barrier dam, so that some proper engineering measures can be implemented and flood discharge can be carried out well.展开更多
Past investigations of the hydrodynamic forces on vertical columns have generally been based on rigid structure assumptions.The effects of structural flexibility and geometry characteristics on the hydrodynamic force ...Past investigations of the hydrodynamic forces on vertical columns have generally been based on rigid structure assumptions.The effects of structural flexibility and geometry characteristics on the hydrodynamic force distribution are not well understood.In this study,fluid-structure interaction models are developed for numerical analyses.This modeling technique is verified with an experimental test in the literature using both circular and rectangular cross-sections.A series of material elasticities that present structural properties ranging from rigid to flexible is then used to conduct analyses.This finding indicates that an increase in structural flexibility can decrease the impact force to some extent,but this effect is limited.A concrete bridge pier with fluid flow impact can be considered rigid when it is fixed at the bottom.After that,the effects of the initial downstream water height and the width of water tank on the hydrodynamic force are thoroughly investigated.The results demonstrate that the increase in the downstream water height with a constant upstream water height corresponds to a decreased force.Moreover,the vertical column results in a blockage effect on the fluid flow.The greater the blockage effect,the higher the hydrodynamic force.The blockage effect from the vertical column can be neglected when the tank width is greater than eight times the structural cross-section diameter.展开更多
One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approache...One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.展开更多
Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformabl...Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.展开更多
The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to t...The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area I is flooded at 2 500 s, and after 5 500 s, most of industrial area II is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area I and II. Based on the above analysis, the optimal evacuation scheme is determined.展开更多
Dam-break analysis is of great importance in mountain environment,especially where reservoirs are located upstream of densely populated areas and hydraulic hazard should be assessed for land planning purposes.Accordin...Dam-break analysis is of great importance in mountain environment,especially where reservoirs are located upstream of densely populated areas and hydraulic hazard should be assessed for land planning purposes.Accordingly,there is a need to identify suitable operative tools which may differ from the ones used in flat flood-prone areas.This paper shows the results provided by a 1D and a 2D model based on the Shallow Water Equations(SWE) for dam-break wave propagation in alpine regions.The 1D model takes advantage of a topographic toolkit that includes an algorithm for pre-processing the Digital Elevation Model(DEM) and of a novel criterion for the automatic cross-section space refinement.The 2D model is FLO-2D,a commercial software widely used for flood routing in mountain areas.In order to verify the predictive effectiveness of these numerical models,the test case of the Cancano dam-break has been recovered from the historical study of De Marchi(1945),which provides a unique laboratory data set concerning the consequences of the potential collapse of the former Cancano dam(Northern Italy).The measured discharge hydrograph at the dam also provides the data to test a simplified method recently proposed for the characterization of the hydrograph following a sudden dam-break.展开更多
A two-dimensional （2D） dam-break flow numerical model was developed based on the finite-volume total variation diminishing （TVD） and monotone upstream-centered scheme for conservation laws （MUSCL）-Hancock scheme...A two-dimensional （2D） dam-break flow numerical model was developed based on the finite-volume total variation diminishing （TVD） and monotone upstream-centered scheme for conservation laws （MUSCL）-Hancock scheme, which has second-order accuracy in both time and space. A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver was used to evaluate fluxes. The TVD MUSCL-Hancock numerical scheme utilizes slope limiters, such as the minmod, double minmod, superbee, van Albada, and van Leer limiters, to prevent spurious oscillations and maintain monotonicity near discontinuities. A comparative study of the impact of various slope limiters on the accuracy of the numerical flow model was conducted with several dam-break examples including wet and dry bed cases. The numerical results of the superbee and double minmod limiters agree better with the theoretical solution and have higher accuracy than other limiters in one-dimensional （1D） space. The ratio of the downstream water depth to the upstream water depth was used to select the proper slope limiter. For the 2D numerical model, the superbee limiter should not be used, owing to significant numerical dispersion.展开更多
When studying the dam-break flow phenomenon,the basic hydrodynamic features of the dam-break flow at the gate location should be verified primarily.In this study,laboratory experiments were performed in a rectangular ...When studying the dam-break flow phenomenon,the basic hydrodynamic features of the dam-break flow at the gate location should be verified primarily.In this study,laboratory experiments were performed in a rectangular and horizontal flume with the same initial water head setting on the dry and wet downstream bed conditions.Water surface elevation was extracted through image analysis and validated by comparing with the data measured using a wave gauge.Temporal variation of the water surface elevation at the gate location,quantified in terms of high-speed video recorded images,can be divided into three stages,the sharp decreasing stage,the relatively steady stage,and the gradually decreasing stage.Applicability of several classic analytical solutions of the dam-break problem at the gate location was validated using present experimental data.Ritter's solution is effective for the dry bed condition while Stoker's solution could be applied to the wet bed case,and both are only applicable during the steady stage.Lin' solution reproduces the gate-site hydrographs well during both the relatively steady and the gradually decreasing stages,especially for the condition under which the down-upstream water depth ratio is smaller than 0.138.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
Active tectonic movements and geological disasters frequently occur in the upper reaches of the Jinsha River,increasing the likelihood of landslides obstructing the river.Taking the ancient Rongcharong landslide dam f...Active tectonic movements and geological disasters frequently occur in the upper reaches of the Jinsha River,increasing the likelihood of landslides obstructing the river.Taking the ancient Rongcharong landslide dam failure events in the Suwalong reach as an example,this paper first analyzes the accuracy and applicability of the commonly used methods in calculating the peak flow of the dam failure,such as the empirical formula,the numerical method based on the physical mechanism,and the computational fluid dynamics(CFD)method.Then,the peak flood flow of the Rongcharong-dammed lake when it overflows the dam is determined to be 28393-64272 m~3/s.At the same time,the failure process of landslide dam due to flood erosion was elucidated using the CFD method,which can be divided into three stages:gradual erosion in the initial stage,rapid development in the middle stage,and gradual expansion in the final stage.Finally,the factors that affect the peak flow of floods are analyzed,and suggestions for emergency treatment of landslide dams are put forward.The findings of this research can serve as a valuable reference for disaster prevention and mitigation strategies to adapt to the increasing frequency of landslide-induced river blockages.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
In the process of storm surge, the seawater often overflows and even destroys the seawall. The buildings near the shore are usually inundated by the seawater through the breach. However, at present, there is little st...In the process of storm surge, the seawater often overflows and even destroys the seawall. The buildings near the shore are usually inundated by the seawater through the breach. However, at present, there is little study focusing on the effects of buildings and breach on the seawall-break flows. In this paper, the lattice Boltzmann （LB） model with nine velocities in two dimensions （1920,9） for the shallow water equations is adopted to simulate the seawall-break flows. The flow patterns and water depth distributions for the seawall-break flows under various densities, layouts and shapes of buildings and different breach discharges, sizes and locations are investigated. It is found that when buildings with a high enough density are perpendicular to the main flow direction, an obvious backwater phenomenon appears near buildings while this phenomenon does not occur when buildings with the same density are paraJlel to the main flow direction. Moreover, it is observed that the occurrence of backwater phenomenon is independent of the building shape. As to the effects of breach on the seawall-break flows, it is found that only when the breach discharge is large enough or the breach size is small enough, the effects of asymmetric distribution of buildings on the seawail-break flows become important. The breach location only changes the flow pattern in the upstream area of the first building that seawater meets, but has little impact on the global water depth distribution.展开更多
In the current study,the treatment of air/water interface has been made on dam-break induced tsunami-like wave by the Coupled Level Set and Volume of Fluid(CLSVOF)three-dimensional modelling.The overall CLSVOF method ...In the current study,the treatment of air/water interface has been made on dam-break induced tsunami-like wave by the Coupled Level Set and Volume of Fluid(CLSVOF)three-dimensional modelling.The overall CLSVOF method adopts a Tangent of Hy-perbola for INterface Capturing(THINC)scheme with the Weighted Linear Interface Calculation(WLIC)and Level Set(LS)function for capturing interface and calculating normal vector,respectively.As far as THINC/WLIC scheme is concerned,since the convection problem of the VOF function can be solved well,the numerical diffusion can be avoided.The spatial terms in the LS equation were discretized by the Optimized Compact Reconstruction Weighted Essentially Non-Oscillatory(OCRWENO)scheme with fourth-order accuracy,which can avoid false oscillation of LS solution.By combining CLSVOF method with Immersed Boundary(IB)method,the simulation of dam-break induced tsunami-like wave impacting on a stationary breakwater can be carried out.Grid sensitivity,mass error and free-surface profile are first calculated for the tsunami-like wave problem to validate the proposed numerical algorithm,which shows excellent agreement between the numerical results and experimental data.Tsunami-like waves with varied tailgater levels are then investigated.Calculations of velocity magnitude,free-surface profile and wave elevation of the tsunami-like wave are conducted to investigate its dynamics and kinematics.展开更多
Overtopping is one of the main reasons for the breaching of concrete-face sand-gravel dams(CFSGDs).In this study,a refined mathematical model was established based on the characteristics of the overtopping breaching o...Overtopping is one of the main reasons for the breaching of concrete-face sand-gravel dams(CFSGDs).In this study,a refined mathematical model was established based on the characteristics of the overtopping breaching of CFSGDs.The model characteristics were as follows:(1)Based on the Renormailzation Group(RNG)k-εturbulence theory and volume of fluid(VOF)method,the turbulent characteristics of the dam-break flow were simulated,and the erosion surface of the water and soil was tracked;(2)In consideration of the influence of the change in the sediment content on the dam-break flow,the dam material transport equation,which could reflect the characteristics of particle settlement and entrainment motion,was used to simulate the erosion process of the sand gravels;(3)Based on the bending moment balance method,a failure equation of the concrete face slab under dead weight and water load was established.The proposed model was verified through a case study on the failure of the Gouhou CFSGD.The results showed that the proposed model could well simulate the erosion mode of the special vortex flow of the CFSGD scouring the support body of the concrete face slab inward and reflect the mutual coupling relationship between the dam-break flow,sand gravels,and concrete face slabs.Compared with the measured values,the relative errors of the peak discharge,final breach average width,dam breaching duration,and maximum failure length of the face slab calculated using the proposed model were all less than 12%,thus verifying the rationality of the model.The proposed model was demonstrated to perform better and provide more detailed results than three selected parametric models and three simplified mathematical models.The study results can aid in establishing the risk level and devising early warning strategies for CFSGDs.展开更多
In order to model the movement progress in case of risks such as dam collapse and coastal inundation, particle-based simulation methods, including the discrete-element method and smoothed particle hydrodynamics, which...In order to model the movement progress in case of risks such as dam collapse and coastal inundation, particle-based simulation methods, including the discrete-element method and smoothed particle hydrodynamics, which have specific advantages in modeling complex three-dimensional environmental fluid and particulate flows, are adopted as an effective way to illustrate environmental applications possibly happening in the real world. The theory of these methods and their relative advantages compared with traditional methods are discussed. Examples of 3-D flows on realistic topography including the flooding of a river valley as a result of a dam collapse and coastal inundation by a tsunami are introduced. Issues related to validation and quality data availability are also discussed. The results show that the simulations provide a valuable insight in a given situation for the security management of reservoir dams. Validation can only be performed where both the initial and final states can be very well characterized.展开更多
A dam failure analysis was performed for the City of Laiyang to determine the flood inundation extents and the breach velocities of a hypothetical failure of the Muyu Reservoir Dam in Shandong Province, China. The stu...A dam failure analysis was performed for the City of Laiyang to determine the flood inundation extents and the breach velocities of a hypothetical failure of the Muyu Reservoir Dam in Shandong Province, China. The study was carried out in four steps. First, the possible reasons of the dam failure were analyzed, and the extreme hydrological event was identified as the most critical cause Second, a numerical simulation of the dam break flood was performed. The key factors considered in this simulation include the water level at the dam （the check flood water level）, the inflow to the reservoir （the check flood hydrograph）, and two typical structural scenarios （with/without-levees） along the river reach downstream of the reservoir. Third, based on flood simulations, flood hazard risk analysis was conducted and the results indicate that the peak discharge would be over 40,000 m3/s at the breachand over 22,000 m3/s near Laiyang City. As a result, most of the urban area would be flooded; the duration from the beginning of the dam break to the arrival of peak discharge at Laiyang City is less than half an hour. Finally, a plan for integrated dam break flood management was presented to mitigate the flood risk; both available structural and nonstructural measures were proposed according to the situations in the river basin. The main strategic advises include flood emergency planning, flood risk mapping, flood monitoring and early warning system construction, reservoir regulation, and flood risk public awareness.展开更多
On July 23,a dam-break occurred to the Auxiliary Dam of XePian-XeNamnoy Hydropower Station located in Southern Laos,and 5 billion cubic meters of water flowed out in a few hours,which has caused the most serious flood...On July 23,a dam-break occurred to the Auxiliary Dam of XePian-XeNamnoy Hydropower Station located in Southern Laos,and 5 billion cubic meters of water flowed out in a few hours,which has caused the most serious flood of the past ten years in Laos and incurred the loss of many lives.展开更多
Large-scale landslide dams may block the river flow and cause inundation upstream, and subsequently fail and result in severe flooding and damage in the downstream. The need for enhanced understanding of the inundatio...Large-scale landslide dams may block the river flow and cause inundation upstream, and subsequently fail and result in severe flooding and damage in the downstream. The need for enhanced understanding of the inundation and flooding is evident. This article presents an experimental study of the inundation and landslide dam-break flooding over erodible bed in open channels. A set of automatic water-level probes is deployed to record the highly transient stage, and the post-flooding channel bed elevation is measured. New experimental data resources are provided for understanding the processes of landslide-induced flooding and for testing mathematical rivers models.展开更多
Dam break can cause a significant disaster in the downstream, especially, in a valley with cascade reservoirs, which would aggravate the disaster extent. The experimental studies of the dam-break flow of cascade reser...Dam break can cause a significant disaster in the downstream, especially, in a valley with cascade reservoirs, which would aggravate the disaster extent. The experimental studies of the dam-break flow of cascade reservoirs are few and far between at the present, Most of related studies concern the failure of a single dam.. This article presents an experimental study of the characteristics of an instantly filled dam-break flow of cascade reservoirs in a rectangular glass flume with a steep bottom slope. A new method was used to simulate the sudden collapse of the dam. A series of sensors for automatic water-levels were deployed to record the rapid water depth fluctuation. The experimental results show that, the ratio of the initial water depth of the downstream reservoir to that of the upstream reservoir would greatly affect the flood peak water depth in the downstream reservoir area and in the stream channel behind the downstream dam, while the influence of the dam spacing is insignificant. In addition, the comparison between the single reservoir and the cascade reservoirs shows some difference in the dam-break flow pattern and the stage hydrograph at the corresponding gauging points.展开更多
基金funding from the National Natural Science Foundation Project (Grant No. 40772175, 40972175)the Scientific research fund of Southwest Jiaotong University (Grant No.2008-A01)+1 种基金the Opening Fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Chengdu University of Technology) (Grant No. DZKJ–08012)the National Natural Science Foundation Project-mutual fund of Yunnan Province (Grant No.U1033601)
文摘Tangjiashan landslide is a typical high-speed consequent landslide of medium-steep dip angle. This landslide triggered by earthquake took place in about semi-minute. The relative sliding displacement is 900 meters, so average sliding speed is about 30 meters per second. The longitudinal length of barrier dam which is formed by high-speed landslide along river is 803.4 meters; and maximum width crossing river is 611.8 meters. And its volume is estimated about 20.37 million steres. Through detailed geological investigation of the barrier dam, together with early geological information before earthquake, geological structures of the barrier dam and its stability of upstream and downstream slopes are studied when water level reaches different elevations in condition of continual after shocks with seismic intensity of 7 or 8 Richter scale. On this basis, dam-breaking mode of barrier dam is discussed deeply. Thereby, analytic results provide significant guidance and advices to front headquarters of Tangjiashan barrier dam, so that some proper engineering measures can be implemented and flood discharge can be carried out well.
基金The National Natural Science Foundation of China(No.52222804,U21A20154).
文摘Past investigations of the hydrodynamic forces on vertical columns have generally been based on rigid structure assumptions.The effects of structural flexibility and geometry characteristics on the hydrodynamic force distribution are not well understood.In this study,fluid-structure interaction models are developed for numerical analyses.This modeling technique is verified with an experimental test in the literature using both circular and rectangular cross-sections.A series of material elasticities that present structural properties ranging from rigid to flexible is then used to conduct analyses.This finding indicates that an increase in structural flexibility can decrease the impact force to some extent,but this effect is limited.A concrete bridge pier with fluid flow impact can be considered rigid when it is fixed at the bottom.After that,the effects of the initial downstream water height and the width of water tank on the hydrodynamic force are thoroughly investigated.The results demonstrate that the increase in the downstream water height with a constant upstream water height corresponds to a decreased force.Moreover,the vertical column results in a blockage effect on the fluid flow.The greater the blockage effect,the higher the hydrodynamic force.The blockage effect from the vertical column can be neglected when the tank width is greater than eight times the structural cross-section diameter.
基金funded by Natural Science Foundation of China (Grants No. 11172217 and 11432015)National Key Basic Research and Development Program (i.e., 973 Program) of China (Grant No. 2007CB714106)
文摘One of the largest known megafloods on earth resulted from a glacier dam-break,which occurred during the Late Quaternary in the Altai Mountains in Southern Siberia.Computational modeling is one of the viable approaches to enhancing the understanding of the flood events.The computational domain of this flood is over 9460 km2 and about 3.784 × 106 cells are involved as a 50 m × 50 m mesh is used,which necessitates a computationally efficient model.Here the Open MP(Open Multiprocessing) technique is adopted to parallelize the code of a coupled 2D hydrodynamic and sediment transport model.It is shown that the computational efficiency is enhanced by over 80% due to the parallelization.The floods over both fixed and mobile beds are well reproduced with specified discharge hydrographs at the dam site.Qualitatively,backwater effects during the flood are resolved at the bifurcation between the Chuja and Katun rivers.Quantitatively,the computed maximum stage and thalweg are physically consistent with the field data of the bars and deposits.The effects of sediment transport and morphological evolution on the flood are considerable.Sensitivity analyses indicate that the impact of the peak discharge is significant,whilst those of the Manningroughness,medium sediment size and shape of the inlet discharge hydrograph are marginal.
基金National Natural Science Foundation of China under Grant No.90510017
文摘Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.
基金Supported by Natural Science Foundation of Tianjin (No.09JCYBJC08700)the Foundation for Innovative Research Groups of National Natural Science Foundation of China (No.51021004)National Natural Science Foundation of China (No.90815019)
文摘The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area I is flooded at 2 500 s, and after 5 500 s, most of industrial area II is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area I and II. Based on the above analysis, the optimal evacuation scheme is determined.
基金developed within the European Project Kulturisk (Grant agreement 265280)
文摘Dam-break analysis is of great importance in mountain environment,especially where reservoirs are located upstream of densely populated areas and hydraulic hazard should be assessed for land planning purposes.Accordingly,there is a need to identify suitable operative tools which may differ from the ones used in flat flood-prone areas.This paper shows the results provided by a 1D and a 2D model based on the Shallow Water Equations(SWE) for dam-break wave propagation in alpine regions.The 1D model takes advantage of a topographic toolkit that includes an algorithm for pre-processing the Digital Elevation Model(DEM) and of a novel criterion for the automatic cross-section space refinement.The 2D model is FLO-2D,a commercial software widely used for flood routing in mountain areas.In order to verify the predictive effectiveness of these numerical models,the test case of the Cancano dam-break has been recovered from the historical study of De Marchi(1945),which provides a unique laboratory data set concerning the consequences of the potential collapse of the former Cancano dam(Northern Italy).The measured discharge hydrograph at the dam also provides the data to test a simplified method recently proposed for the characterization of the hydrograph following a sudden dam-break.
基金supported by the National Natural Science Foundation of China(Grants No.51679170,51379157,and 51439007)
文摘A two-dimensional （2D） dam-break flow numerical model was developed based on the finite-volume total variation diminishing （TVD） and monotone upstream-centered scheme for conservation laws （MUSCL）-Hancock scheme, which has second-order accuracy in both time and space. A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver was used to evaluate fluxes. The TVD MUSCL-Hancock numerical scheme utilizes slope limiters, such as the minmod, double minmod, superbee, van Albada, and van Leer limiters, to prevent spurious oscillations and maintain monotonicity near discontinuities. A comparative study of the impact of various slope limiters on the accuracy of the numerical flow model was conducted with several dam-break examples including wet and dry bed cases. The numerical results of the superbee and double minmod limiters agree better with the theoretical solution and have higher accuracy than other limiters in one-dimensional （1D） space. The ratio of the downstream water depth to the upstream water depth was used to select the proper slope limiter. For the 2D numerical model, the superbee limiter should not be used, owing to significant numerical dispersion.
基金financially supported by the Natural Science Foundation of Zhejiang Province,China(No.LR14E090002)the National Natural Science Foundation of China(No.11632012)the Open Research Fund Program of State key Laboratory of Hydroscience and Engineering(No.sklhse-2016-B-02)
文摘When studying the dam-break flow phenomenon,the basic hydrodynamic features of the dam-break flow at the gate location should be verified primarily.In this study,laboratory experiments were performed in a rectangular and horizontal flume with the same initial water head setting on the dry and wet downstream bed conditions.Water surface elevation was extracted through image analysis and validated by comparing with the data measured using a wave gauge.Temporal variation of the water surface elevation at the gate location,quantified in terms of high-speed video recorded images,can be divided into three stages,the sharp decreasing stage,the relatively steady stage,and the gradually decreasing stage.Applicability of several classic analytical solutions of the dam-break problem at the gate location was validated using present experimental data.Ritter's solution is effective for the dry bed condition while Stoker's solution could be applied to the wet bed case,and both are only applicable during the steady stage.Lin' solution reproduces the gate-site hydrographs well during both the relatively steady and the gradually decreasing stages,especially for the condition under which the down-upstream water depth ratio is smaller than 0.138.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
基金supported by the Second Tibetan Plateau Scientific Expedition and Research Program(STEP)(No.2019QZKK0904)the National Natural Science Foundation of China(No.42307247)。
文摘Active tectonic movements and geological disasters frequently occur in the upper reaches of the Jinsha River,increasing the likelihood of landslides obstructing the river.Taking the ancient Rongcharong landslide dam failure events in the Suwalong reach as an example,this paper first analyzes the accuracy and applicability of the commonly used methods in calculating the peak flow of the dam failure,such as the empirical formula,the numerical method based on the physical mechanism,and the computational fluid dynamics(CFD)method.Then,the peak flood flow of the Rongcharong-dammed lake when it overflows the dam is determined to be 28393-64272 m~3/s.At the same time,the failure process of landslide dam due to flood erosion was elucidated using the CFD method,which can be divided into three stages:gradual erosion in the initial stage,rapid development in the middle stage,and gradual expansion in the final stage.Finally,the factors that affect the peak flow of floods are analyzed,and suggestions for emergency treatment of landslide dams are put forward.The findings of this research can serve as a valuable reference for disaster prevention and mitigation strategies to adapt to the increasing frequency of landslide-induced river blockages.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
基金Supported by the National Natural Science Foundation of China under Grant No.11502124the Natural Science Foundation of Zhejiang Province under Grant No.LQ16A020001+2 种基金the Scientific Research Fund of Zhejiang Provincial Education Department under Grant No.Y201533808the Natural Science Foundation of Ningbo under Grant No.2016A610075sponsored by K.C.Wong Magna Fund in Ningbo University
文摘In the process of storm surge, the seawater often overflows and even destroys the seawall. The buildings near the shore are usually inundated by the seawater through the breach. However, at present, there is little study focusing on the effects of buildings and breach on the seawall-break flows. In this paper, the lattice Boltzmann （LB） model with nine velocities in two dimensions （1920,9） for the shallow water equations is adopted to simulate the seawall-break flows. The flow patterns and water depth distributions for the seawall-break flows under various densities, layouts and shapes of buildings and different breach discharges, sizes and locations are investigated. It is found that when buildings with a high enough density are perpendicular to the main flow direction, an obvious backwater phenomenon appears near buildings while this phenomenon does not occur when buildings with the same density are paraJlel to the main flow direction. Moreover, it is observed that the occurrence of backwater phenomenon is independent of the building shape. As to the effects of breach on the seawall-break flows, it is found that only when the breach discharge is large enough or the breach size is small enough, the effects of asymmetric distribution of buildings on the seawail-break flows become important. The breach location only changes the flow pattern in the upstream area of the first building that seawater meets, but has little impact on the global water depth distribution.
基金financially supported by the National Key R&D Program of China(Grant No.2016YFC0503562207)the National Natural Science Foundation of China(Grant No.51979178)+4 种基金the Department of Science and Technology of Sichuan Province(Grant No.2019YJ0118)the Fundamental Research Funds for the Central Universities(Grant No.YJ201837)the Innovation Spark Project(Grant No.SCUH0049)the Opening Foundation of Sichuan UniversityState Key Laboratory of Hydraulics and Mountain River Engineering(Grant No.Skhl1820)。
文摘In the current study,the treatment of air/water interface has been made on dam-break induced tsunami-like wave by the Coupled Level Set and Volume of Fluid(CLSVOF)three-dimensional modelling.The overall CLSVOF method adopts a Tangent of Hy-perbola for INterface Capturing(THINC)scheme with the Weighted Linear Interface Calculation(WLIC)and Level Set(LS)function for capturing interface and calculating normal vector,respectively.As far as THINC/WLIC scheme is concerned,since the convection problem of the VOF function can be solved well,the numerical diffusion can be avoided.The spatial terms in the LS equation were discretized by the Optimized Compact Reconstruction Weighted Essentially Non-Oscillatory(OCRWENO)scheme with fourth-order accuracy,which can avoid false oscillation of LS solution.By combining CLSVOF method with Immersed Boundary(IB)method,the simulation of dam-break induced tsunami-like wave impacting on a stationary breakwater can be carried out.Grid sensitivity,mass error and free-surface profile are first calculated for the tsunami-like wave problem to validate the proposed numerical algorithm,which shows excellent agreement between the numerical results and experimental data.Tsunami-like waves with varied tailgater levels are then investigated.Calculations of velocity magnitude,free-surface profile and wave elevation of the tsunami-like wave are conducted to investigate its dynamics and kinematics.
基金supported by the National Science Fund for Distinguished Young Scholars(Grant No.52125904)the National Natural Science Foundation of China(Grant No.51979224)the Program 2022TD-01 for Shaanxi Provincial Innovative Research Team(Grant No.2022TD-01)。
文摘Overtopping is one of the main reasons for the breaching of concrete-face sand-gravel dams(CFSGDs).In this study,a refined mathematical model was established based on the characteristics of the overtopping breaching of CFSGDs.The model characteristics were as follows:(1)Based on the Renormailzation Group(RNG)k-εturbulence theory and volume of fluid(VOF)method,the turbulent characteristics of the dam-break flow were simulated,and the erosion surface of the water and soil was tracked;(2)In consideration of the influence of the change in the sediment content on the dam-break flow,the dam material transport equation,which could reflect the characteristics of particle settlement and entrainment motion,was used to simulate the erosion process of the sand gravels;(3)Based on the bending moment balance method,a failure equation of the concrete face slab under dead weight and water load was established.The proposed model was verified through a case study on the failure of the Gouhou CFSGD.The results showed that the proposed model could well simulate the erosion mode of the special vortex flow of the CFSGD scouring the support body of the concrete face slab inward and reflect the mutual coupling relationship between the dam-break flow,sand gravels,and concrete face slabs.Compared with the measured values,the relative errors of the peak discharge,final breach average width,dam breaching duration,and maximum failure length of the face slab calculated using the proposed model were all less than 12%,thus verifying the rationality of the model.The proposed model was demonstrated to perform better and provide more detailed results than three selected parametric models and three simplified mathematical models.The study results can aid in establishing the risk level and devising early warning strategies for CFSGDs.
基金The National High Technology Research and Development Program of China (863Program) (No.2006AA12Z214)the National Fundamental Surveying and Mapping Project of China (No.1460130524207)the Key Laboratory of Geo-Informatics of the State Bureau ofSurveying and Mapping Project (No.A1717)
文摘In order to model the movement progress in case of risks such as dam collapse and coastal inundation, particle-based simulation methods, including the discrete-element method and smoothed particle hydrodynamics, which have specific advantages in modeling complex three-dimensional environmental fluid and particulate flows, are adopted as an effective way to illustrate environmental applications possibly happening in the real world. The theory of these methods and their relative advantages compared with traditional methods are discussed. Examples of 3-D flows on realistic topography including the flooding of a river valley as a result of a dam collapse and coastal inundation by a tsunami are introduced. Issues related to validation and quality data availability are also discussed. The results show that the simulations provide a valuable insight in a given situation for the security management of reservoir dams. Validation can only be performed where both the initial and final states can be very well characterized.
文摘A dam failure analysis was performed for the City of Laiyang to determine the flood inundation extents and the breach velocities of a hypothetical failure of the Muyu Reservoir Dam in Shandong Province, China. The study was carried out in four steps. First, the possible reasons of the dam failure were analyzed, and the extreme hydrological event was identified as the most critical cause Second, a numerical simulation of the dam break flood was performed. The key factors considered in this simulation include the water level at the dam （the check flood water level）, the inflow to the reservoir （the check flood hydrograph）, and two typical structural scenarios （with/without-levees） along the river reach downstream of the reservoir. Third, based on flood simulations, flood hazard risk analysis was conducted and the results indicate that the peak discharge would be over 40,000 m3/s at the breachand over 22,000 m3/s near Laiyang City. As a result, most of the urban area would be flooded; the duration from the beginning of the dam break to the arrival of peak discharge at Laiyang City is less than half an hour. Finally, a plan for integrated dam break flood management was presented to mitigate the flood risk; both available structural and nonstructural measures were proposed according to the situations in the river basin. The main strategic advises include flood emergency planning, flood risk mapping, flood monitoring and early warning system construction, reservoir regulation, and flood risk public awareness.
文摘On July 23,a dam-break occurred to the Auxiliary Dam of XePian-XeNamnoy Hydropower Station located in Southern Laos,and 5 billion cubic meters of water flowed out in a few hours,which has caused the most serious flood of the past ten years in Laos and incurred the loss of many lives.
基金Project supported by the National Natural Science Foundation of China (Grant No 50459001)
文摘Large-scale landslide dams may block the river flow and cause inundation upstream, and subsequently fail and result in severe flooding and damage in the downstream. The need for enhanced understanding of the inundation and flooding is evident. This article presents an experimental study of the inundation and landslide dam-break flooding over erodible bed in open channels. A set of automatic water-level probes is deployed to record the highly transient stage, and the post-flooding channel bed elevation is measured. New experimental data resources are provided for understanding the processes of landslide-induced flooding and for testing mathematical rivers models.
基金Project supported by the National Basic Research Program of China(973 Program,Grant No.2007CB714105)the National Natural Science Foundation of China(Grant No.50909067)
文摘Dam break can cause a significant disaster in the downstream, especially, in a valley with cascade reservoirs, which would aggravate the disaster extent. The experimental studies of the dam-break flow of cascade reservoirs are few and far between at the present, Most of related studies concern the failure of a single dam.. This article presents an experimental study of the characteristics of an instantly filled dam-break flow of cascade reservoirs in a rectangular glass flume with a steep bottom slope. A new method was used to simulate the sudden collapse of the dam. A series of sensors for automatic water-levels were deployed to record the rapid water depth fluctuation. The experimental results show that, the ratio of the initial water depth of the downstream reservoir to that of the upstream reservoir would greatly affect the flood peak water depth in the downstream reservoir area and in the stream channel behind the downstream dam, while the influence of the dam spacing is insignificant. In addition, the comparison between the single reservoir and the cascade reservoirs shows some difference in the dam-break flow pattern and the stage hydrograph at the corresponding gauging points.
