Numerical Simulation of Dam-Break Flows Using Radial Basis Functions: Application to Urban Flood Inundation

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.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Coupled 2D Hydrodynamic and Sediment Transport Modeling of Megaflood due to Glacier Dam-break in Altai Mountains,Southern Siberia

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.

Numerical simulation of pore-scale flow in chemical flooding process

Chemical flooding is one of the effective technologies to increase oil recovery of petroleum reservoirs after water flooding.Above the scale of representative elementary volume(REV), phenomenological modeling and numerical simulations of chemical flooding have been reported in literatures,but the studies alike are rarely conducted at the pore-scale,at which the effects of physicochemical hydrodynamics are hardly resolved either by experimental observations or by traditional continuum-based simulations.In this paper,dissipative particle dynamics(DPD),one of mesoscopic fluid particle methods,is introduced to simulate the pore-scale flow in chemical flooding processes.The theoretical background,mathematical formulation and numerical approach of DPD are presented.The plane Poiseuille flow is used to illustrate the accuracy of the DPD simulation,and then the processes of polymer flooding through an oil-wet throat and a water-wet throat are studies, respectively.The selected parameters of those simulations are given in details.These preliminary results show the potential of this novel method for modeling the physicochemical hydrodynamics at the pore scale in the area of chemical enhanced oil recovery.

Numerical Simulation of Dam-Break Flooding of Cascade Reservoirs

Previous studies at home and abroad have mainly focused on single dam-break, and little attention has been paid so far to the dam-break of cascade reservoirs. Multi-source flooding, which can lead to three-dimensional turbulent phenomena and superposition effects, is the main difference between the dam-break of cascade and single reservoirs. Detailed descriptions of the coupled numerical simulation of multi-source flooding have little been reported, and the initial wet riverbed is rarely considered in current models. Therefore, in this paper, a method based on the three-dimensional (Formula presented.) turbulence model coupled with the volume-of-fluid method is proposed to simulate the dam-break flooding of cascade reservoirs. The upstream river, reservoir, and downstream river are connected by the internal boundary method, and the initial conditions, including river flow and reservoir water, are determined according to the results of the numerical simulation. Coupled numerical simulation of different dam-break flooding is then achieved. The present work solves the challenges presented by the enhancement and superposition of natural river flow, upstream flooding, and downstream flooding. This paper provides a theoretical basis for future studies on the dam-break flood routing of cascade reservoirs. © 2017 Tianjin University and Springer-Verlag GmbH Germany

Study of total variation diminishing (TVD) slope limiters in dam-break flow simulation

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.

Determination of minimum flood flow for regeneration of floodplain forest from inundated forest width-stage curve

Floods are essential for the regeneration and growth of floodplain forests in arid and semiarid regions. However, river flows, and especially flood flows, have decreased greatly with the increase of water diversion from rivers and/or reservoir regulation, resulting in severe deterioration of floodplain ecosystems. Estimation of the flood stage that will inundate the floodplain forest is necessary for the forest＇s restoration or protection. To balance water use for economic purposes and floodplain forest protection, the inundated forest width method is proposed for estimating the minimum flood stage for floodplain forests from the inundated forest width-stage curve. The minimum flood stage is defined as the breakpoint of the inundated forest width-stage curve, and is determined directly or analytically from the curve. For the analytical approach, the problem under consideration is described by a multi-objective optimization model, which can be solved by the ideal point method. Then, the flood flow at the minimum flood stage （minimum flood flow）, which is useful for flow regulation, can be calculated from the stage-discharge curve. In order to protect the forest in a river floodplain in a semiarid area in Xinjiang subject to reservoir regulation upstream, the proposed method was used to determine the minimum flood stage and flow for the forest. Field survey of hydrology, topography, and forest distribution was carried out at typical cross sections in the floodplain. Based on the survey results, minimum flood flows for six typical cross sections were estimated to be between 306 m3/s and 393 m3/s. Their maximum, 393 m3/s, was considered the minimum flood flow for the study river reach. This provides an appropriate flood flow for the protection of floodplain forest and can be used in the regulation of the upstream reservoir.

Three Dimensional Numerical Modeling of Flow and Pollutant Transport in a Flooding Area of 2008 US Midwest Flood

This paper presents the development and application of a three-dimensional numerical model for simulating the flow field and pollutant transport in a flood zone near the confluence of the Mississippi River and Iowa River during the US Midwest Flood in 2008. Due to a prolonged precipitation event, a levee along the Iowa River just upstream of Oakville, Iowa broke, and the small town was completely flooded for a couple of weeks. During this period, the high water level in the flood zone reached about 2.5 metersabove the ground, and wind was the major force for the flow circulation. It was observed that some pollutants were leaked from the residential and farming facilities and transported into the flood zone. Leaking of pollutants from these facilities was reported by different news media during the flood and was identified using high resolution satellite imagery. The developed 3D numerical model was first validated using experimental measurements, and then applied to the flood inundated zone in Oakville for simulating the unsteady hydrodynamics and pollutant transport. The simulated pollutant distributions were generally in good agreement with the observed data obtained from satellite imagery.

Water flooding flowing area identification for oil reservoirs based on the method of streamline clustering artificial intelligence

For the case of carbonate reservoir water flooding development, the flow field identification method based on streamline modeling result was proposed. The Ocean for Petrel platform was used to build the plug-in that exported the streamline data, and the subsequent data was processed and clustered through Python programming, to display the flow field with different water flooding efficiencies at different time in the reservoir. We used density peak clustering as primary streamline cluster algorithm, and Silhouette algorithm as the cluster validation algorithm to select reasonable cluster number, and the results of different clustering algorithms were compared. The results showed that the density peak clustering algorithm could provide better identified capacity and higher Silhouette coefficient than K-means, hierachical clustering and spectral clustering algorithms when clustering coefficients are the same. Based on the results of streamline clustering method, the reservoir engineers can easily identify the flow area with quantification treatment, the inefficient water injection channels and area with developing potential in reservoirs can be identified. Meanwhile, streamlines between the same injector and producer can be subdivided to describe driving capacity distribution in water phase, providing useful information for the decision making of water flooding optimization, well pattern adjustment and deep profile modification.

Transient pressure analysis of polymer flooding fractured wells with oil-water two-phase flow

The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.

Debris Fan Produced by Failure of Canyon-Blocking Pyroclastic Flows

Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km3 to 0.08 km3. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.

一种融合Priority-Flood算法与D8算法特点的河网提取方法

基于数字高程模型(Digital Elevation Models)进行河网提取是数字地形分析的重要研究内容和水文分析的关键步骤,而流向算法是其中最重要的部分。然而,传统的流向算法在平坦地形区域容易形成平行河网。为解决平行河网问题,该文提出了一种结合Priority-Flood算法漫水思想和改进D8算法的流向算法。借助Priority-Flood算法的漫水思想,可以保证流向符合宏观地形特征;同时将地形划分为上坡、平地和下坡3种类型,辅以PriorityFlood算法可以实现在不填洼的基础上计算流向,而且针对不同地形,采用相应的改进D8流向算法,从而保证流向符合微观地形特征。通过这两种算法的结合,可以直接使用原始DEM提取河网,消除了填洼过程,提高了流向计算效率,同时消除了填洼造成的大量平地,有效地解决了平地区域平行河网问题。最后,选取全球不同区域SRTM数据对J&D算法与新算法进行对比实验,结果表明,新算法能够有效地改善平行河网问题。

一种可抗TCP Flooding攻击的网络流量监测机制

DDoS攻击是互联网的主要安全威胁之一,而大部分DDoS攻击工具都使用TCP Flooding攻击方式,基于大量研究相关技术的基础上,提出了一种可用于局域网的网络流量监测机制,可以有效的检测出TCP Flooding攻击,解决当前各种网络安全设备在此方面存在的问题.

基于MIKE FLOOD耦合模型的新建城区防洪排涝模拟研究

为有效应对城区暴雨洪涝灾害,基于MIKE FLOOD平台,将MIKE11、MIKE21和MIKE URBAN模型耦合,构建了城区内涝耦合模型。选取新建城区东山湖产业园进行模拟。由于产业园新区无实测水文资料,采用管道流量自检验法验证了所构建的模型的合理性,模拟结果与实际情况较吻合。结果表明,在遭遇强降水时,区内出现多个溢流点,道路出现不同程度积水,部分管道处于有压流状态,检查井也多数发生溢流,无法满足排水需求。针对园区防洪排涝存在的问题,提出内涝防治方案,以MIKE FLOOD为基础,实现低影响开发(Low Impact Development,LID)措施组合布设。措施实施后,园区内地表淹没面积减少,积水深度普遍降低,积水点出现积水的时间推迟且积水总历时减小,管道排水口流量峰值明显减小,排水流量峰值的出现延迟。研究为园区防灾减灾工作提供技术支撑,也为类似城区的建设提供参考。

Effects of physical parameter range on dimensionless variable sensitivity in water flooding reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.

基于流指纹的DDoS flooding攻击检测方法

针对现有基于流量特征的DDoS泛洪攻击检测方法易被攻陷,且难以区分攻击与突发访问事件(flash crowds)等问题,提出了一种泛化的攻击检测策略。在此基础上,分别构建泛洪行为和泛洪攻击两种流指纹,通过目标超点聚合度定位泛洪行为,采用一种滑动判别算法区分泛洪攻击与flash crowds事件。实验结果表明,该检测方法可以有效检测泛洪攻击。由于该方法只需维护流到达信息,对于实现高速网络环境下的攻击检测具有一定的实际意义。

Flood susceptibility modelling using advanced ensemble machine learning models

Floods are one of nature's most destructive disasters because of the immense damage to land,buildings,and human fatalities.It is difficult to forecast the areas that are vulnerable to flash flooding due to the dynamic and complex nature of the flash floods.Therefore,earlier identification of flash flood susceptible sites can be performed using advanced machine learning models for managing flood disasters.In this study,we applied and assessed two new hybrid ensemble models,namely Dagging and Random Subspace(RS)coupled with Artificial Neural Network(ANN),Random Forest(RF),and Support Vector Machine(SVM)which are the other three state-of-the-art machine learning models for modelling flood susceptibility maps at the Teesta River basin,the northern region of Bangladesh.The application of these models includes twelve flood influencing factors with 413 current and former flooding points,which were transferred in a GIS environment.The information gain ratio,the multicollinearity diagnostics tests were employed to determine the association between the occurrences and flood influential factors.For the validation and the comparison of these models,for the ability to predict the statistical appraisal measures such as Freidman,Wilcoxon signed-rank,and t-paired tests and Receiver Operating Characteristic Curve(ROC)were employed.The value of the Area Under the Curve(AUC)of ROC was above 0.80 for all models.For flood susceptibility modelling,the Dagging model performs superior,followed by RF,the ANN,the SVM,and the RS,then the several benchmark models.The approach and solution-oriented outcomes outlined in this paper will assist state and local authorities as well as policy makers in reducing flood-related threats and will also assist in the implementation of effective mitigation strategies to mitigate future damage.

Analysis of dynamic wave model for flood routing in natural rivers

Flooding is a common natural disaster that causes enormous economic, social, and human losses. Of various flood routing methods, the dynamic wave model is one of the best approaches for the prediction of the characteristics of floods during their propagations in natural rivers because all of the terms of the momentum equation are considered in the model. However, no significant research has been conducted on how the model sensitivity affects the accuracy of the downstream hydrograph. In this study, a comprehensive analysis of the input parameters 9f the dynamic wave model was performed through field applications in natural rivers and routing experiments in artificial channels using the graphical multi-parametric sensitivity analysis （GMPSA）. The results indicate that the effects of input parameter errors on the output results are more significant in special situations, such as lower values of Manning＇s roughness coefficient and/or a steeper bed slope on the characteristics of a design hydrograph, larger values of the skewness factor and/or time to peak on the channel characteristics, larger values of Manning＇s roughness coefficient and/or the bed slope on the space step, and lower values of Manning＇s roughness coefficient and/or a steeper bed slope on the time step and weighting factor.

Identifying the driving factors of sediment delivery ratio on individual flood events in a long-term monitoring headwater basin

The Sediment Delivery Ratio(SDR) has multi-fold environmental implications both in evaluating the soil and water losses and the effectiveness of conservation measures in watersheds. Various factors, including hydrological regime and watershed properties, may influence the SDR at interannual timescales. However, the effect of certain important dynamic factors, such as rainfall peak distribution, runoff erosion power and sediment bulk density, on the sediment delivery ratio of single flood events(SDRe) has received little attention. The Qiaogou headwater basin is in the hilly-gully region of the Chinese Loess Plateau, and it encompasses a 0.45 km^2 catchment. Three large-scale field runoff plots at different geomorphological positions were chosen to obtain the observation data, and the 20-year period between 1986 and 2005 is presented. The results showed that the SDRe of the Qiaogou headwaters varied from 0.49 to 2.77. Among the numerous influential factors, rainfall and runoff were the driving factors causing slope erosion and sediment transport. The rainfall erosivity had a significant positive relationship with the sediment transport modulus(R^2=0.85, P<0.01) but had no significant relationship with SDRe. The rainfall peak coefficient was significantly positively correlated with the SDRe(R^2=0.64, P<0.05), indicating the influence of rainfall energy distribution on the SDRe. The runoff erosion power index was not only significantly related to the sediment transport modulus(R^2=0.84, P<0.01) but also significantly related to the SDRe(R^2=0.57, P<0.01). In addition, the relative bulk density was significantly related to the SDRe, indicating that hyper-concentrated flow characteristics contributed to more transported sediment in the catchment. Thus, the rainfall peak coefficient, runoff erosion power and sediment relative bulk density could be used as dynamic indexes to predict the SDRe in the hilly areas of the Chinese Loess Plateau.

Flow Behavior of Clay-Silt to Sand-Silt Water-Rich Suspensions at Low to High Shear Rates: Implications for Slurries, Transitional Flows, and Submarine Debris-Flows

Water-rich clay to sand suspensions show a shear rate dependent flow behavior and knowledge of the appropriate rheological model is relevant for sedimentological, industrial and hydraulic studies. We present experimental rheological measurements of water-rich(40 to 60 wt%) clay to silt(population A) and silt to sand(population B) suspensions mixed in different proportions. The data evidence a shear rate dependent shear thinning-shear thickening transition. At lower shear rates, the suspensions organize in chains of particles, whereas at higher shear rates, these chains disrupt so increasing the viscosity. The viscosity, consistency and yield stress decrease as the A+B fraction decreases as the content of B particles increases. This behavior reflects the competing effects of the lubrication and frictional processes as a function of particle size and water content. Transitional flows form by the incorporation of small amounts of the finer fraction while ‘oceanic floods’ form at the estuary of rivers and the submarine debris-flows increase their velocity by incorporating water. The critical Reynolds number of the studied suspensions is ～2000±100 suggesting that the grainsize plays a major role in the laminar to turbulent transition. Our results have implications for the modeling of sediment flows and the hazard related to floods.