Debris flow simulation 2D(DFS 2D):Numerical modelling of debris flows and calibration of friction parameters

Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.

2D Coupled 3D:A New Numerical Model for Dual - Structured - Aquifer System

This paper points out that a successful numerical simulation is to construct a correct conceptional model which is very dose to the natural condition. A new model, two dimensional coupled three dimensional model (2D -3D ) is presented in the Present paper,which is the most suitable one for the dual - structured - aquifer system. An example of Wenyinghu area is shown.By using the 2D-3D model, a satisfied result of the simulated area is achieved.

Numerical Modeling of Sediment Transport and Its Effect on Algal Biomass Distribution in Lake Pontchartrain Due to Flood Release from Bonnet CarréSpillway

In order to protect the city of New Orleans from the Mississippi River flooding, the Bonnet Carré Spillway (BCS) was constructed from 1929 to 1936 to divert flood water from the river into Lake Pontchartrain and then into the Gulf of Mexico. During the BCS opening for flood release, large amounts of freshwater, nutrients, sediment, etc. were discharged into Lake Pontchartrain, and caused a lot of environmental problems. To evaluate the environmental impacts of the flood water on lake ecosystems, a two-dimensional numerical model was developed based on CCHE2D and applied to simulate the flow circulation, sediment transport and algal biomass distribution in Lake Pontchartrain. The effect of sediment concentration on the growth of algae was considered in the model. The numerical model was calibrated using field measured data provided by USGS, and then it was validated by the BCS Opening Event in 1997. The simulated results were generally in good agreement with filed data and satellite imagery. The field observation and numerical model show that during the spillway opening for flood release, the sediment concentration is very high, which greatly restricts the growth of algae, so there is no algal bloom observed in the lake. After the closure of BCS, the sediment concentration in the lake reduces gradually, and the nutrient concentration of the lake is still high. Under these conditions, numerical results and satellite imagery showed that the chlorophyll concentration was high and algal bloom might occur.

Numerical Analysis of Steady Smoldering of Biomass Rods

Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.

Numerical Simulation of Residual Circulation due to Bottom Roughness Variability Under Tidal Flows in A Semi-Enclosed Bay

Nowadays there are some chronic serious environmental problems, such as eutrophication, blue tide and so on, in a complicated coastal zone or a semi-enclosed bay, because the water exchanges between an inner bay and an outer sea is weak compared with the supply of contaminant. Under this situation, a method to improve the water quality by 3-dimensional small unsymmetrical structures has been proposed by Komatsu et al. In this paper, several numerical simulations of the tidal current and concentration for various arrangements of bottom roughness in a semi-enclosed model bay are carfled out with a depth-averaged 2-D numerical model. The model is solved by the hybrid finite analytic method with nonstaggered grid. And the SIMPLES algorithm with Rhie and Chow＇ s momentum interpolation technique is used for the simulation. The effect of Komatsu＇ s method for water purification is examined by numerical simulation. The result of numerical experiment indicates that it is possible to generate a new tidal residual current and to activate a tidal exchange by bottom roughness arrangement only.

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.

A 2-D Non-local Closure Model for Atmospheric Boundary Layer Simulations

In this paper a new approach for PBL simulation, the non-local closure scheme based on the transient turbulence theory has been used. It was set up as an alternative to local closure schemes which physical concept is reasonable and distinct. A 2-D non-local closure model was developed in order to study the PBL structure and simulatesome interesting atmospheric processes over non-ulliform underlying surface, especially under the convective and unique weather conditions, such as sea-land circulation and the TIBL structure. The modelled results show good agreement with field measurement.

Numerical Simulation of the Effects on Low Temperature Drainage from LNG on the Ocean Water Environment

In order to research the distribution of thermal field of low temperature discharge from LNG in the ocean, further to formulate reasonable wastewater marine disposal scheme and provide scientific support to the protection of ocean environment, a two-dimensional hydrodynamic model was used to simulate the tide current field in the radial sand ridge group of shallow sea area of Jiangsu maritime space. This model was based on the tide current characteristics near the LNG project and the seasonal design of the typical hydrological conditions in different seasons. This model was used as the hydrodynamic condition that simulated the transport and diffusion for the heat quantity in the change of time and spatial distribution. Considering the comprehensive heat radiation coefficient of water surface, the 2-D water temperature model was established to predict the spatial distribution and temporal variation of the temperature field in the receiving waters. Finally, results showed that the low temperature drainage less affected in the ocean water environment.

Effective Electromagnetic Log Data Interpretation in Realistic Reservoir Models

This paper analyzes some specific features of the numerical interpretation of high-frequency electromagnetic logging data in vertical, deviated and horizontal boreholes entering oil- and water-saturated formations. The interpretation is based on numerical modeling for signals.

黄河下游典型滩区“92·8”洪水分区滞洪沉沙数值模拟研究

采用平面二维水沙数学模型,依据"92·8"典型洪水情况,对黄河下游夹河滩至高村河段进行了无生产堤条件下的自然漫滩、现状生产堤条件下漫滩及典型滩区分洪淤滩,3种工况条件下的河道冲淤、洪水演进及滩区淹没的数值模拟.模拟成果分析表明,对黄河下游蓄滞洪区进行滞洪沉沙,有效增大了分洪区的削峰滞洪作用,保证和扩大了滩槽泥沙横向交换,维持了河槽活力,一定程度上可以缓解黄河下游"二级悬河"的状况.

NUMERICAL MODELING OF SUSPENDED SEDIMENT TRANSPORT IN CHANNEL BENDS

An algorithm to compute three-dimensional sediment transport effect was proposed in this paper to enhance the capability of depth-averaged numerical models. This algorithm took into account of non-uniform distributions of flow velocities and suspended sediment concentrations along water depth, it significantly enhanced the applicability of 2D models in simulating open channel flows, especially in channel bends. Preliminary numerical experiments in a U-shaped and a sine-generated experimental channel indicate that the proposed method performs quite well in predicting the change of bed-deformation in channel bends due to suspended sediment transport. This method provides an effective alternative for the simulations of channel morphodynamic changes.

Two-dimensional numerical modeling of the longitudinal and lateral channel deformations in alluvial rivers

Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbanks were improved. Then a two-dimensional numerical model of the channel deformation was developed, consisting of a 2D flow and sediment transport submodel and bank-erosion submodels of different soil riverbanks. In the model, a new technique for updating the bank geometry during the bed evolution was presented, which combines closely two kinds of submodels. The proposed model is capable of not only predicting the processes of flood routing and longitudinal channel deformation in natural rivers, but also simulating the processes of lateral channel deformation, especially the processes of lateral erosion and failure of cohesive, non-cohesive and composite riverbanks.

Regional-Scale Study on Sediment Processes of Khuran Strait at Persian Gulf with Implications for Engineering Design

In this article, the sediment transport processes in the Khuran Strait between the mainland Iran and Qeshm Island at North Central Persian Gulf are studied in regional scale in a comprehensive manner. The objectives of this study include defining the type and origin of the sediment, the influencing factors, and the dominant mode of sediment transport. Four possible scenarios of sediment processes in terms of sediment sources and influential environmental forces are examined. The results of the conceptual and 2D numerical model of MIKE21 applied for this purposes indicate that the sediment sources in the transport processes are mainly provided by the sediments suspended from the central and eastern zones of the strait bed. Other sources including input from rivers do not have direct influence on the processes. The results are applied to the study of morphological changes for engineering applications including the pattern and amount of deposit in the Rajaee port approach channel and harbor basin. The pattern and amount of annual sediment deposits in the approach channel predicted by the model is satisfactory, compatible with annual dredging records.

2-D NUMERICAL SIMULATION OF RADIONUCLIDE TRANSPORT IN THE LOWER YANGTZE RIVER

The assessment of the radiological impact of the liquid discharges from nuclear power plants is a major issue for the envi- ronmental protection. In this study, a numerical model for the radionuclide transport in the aquatic environment is built, based on the hydrodynamic equations, including the complete set of Saint-Venant equations, the sediment transport equations, with consideration of several different particle sizes and the deposition and erosion of the suspended sediments, and the radionuclide transport equations The exchanges of radionuclides between water, suspended matter and bed sediments are described in terms of kinetic transfer coe- fficients. The model is used to simulate the transport of the radionuclides discharged from a planned nuclear power plant project to be sited along the lower Yangtze River. From the model results, one may see the detailed temporal-spatial evolution of the radio- nuclide contamination in the solution, in the suspended matter and in the bed sediments. The model can be used as a basic tool for studying the environmental impacts of the liquid discharges from nuclear facilities on a river system.

2-D NUMERICAL SIMULATION OF CRUSH BEDROCK RIVER

In this paper, the erosion-resisting coefficient was introduced to computebed deformation in a crush bedrock river. In the case of crush bedrock, there has been no propercontrol equation to describe bed stability, which leads to difficulty in calculation of the beddeformation with conventional methods. The data from field survey were used to give thee-rosion-resisting capability with an appropriate coefficient. After the determination oflongitudinal distribution expressed by polynomial regression and transversal distribution expressedby normal distribution function, the plane distribution of erosion-resisting coefficient in a crushbedrock river was obtained. With the computational results from a 2-D horizontal flow mathematicalmodel, the erosion-resisting coefficient and controlling condition of local stability were employedto compute the values of bed deformation when riverbed is stable. The above method was applied in acase study, and the computational results of flow and bed deformations are in good a-greement withphysical model test data.

Simulation of Artificial Grain Feeding in Order to Reach Dynamical Bed Stabilisation Along the River Rhine

The Wesel-Xanten stretch of the fiver Rhine between km-812.5 and km-821.5 is one of the reaches where strong erosion leads to high maintenance efforts conceming navigability.In order to improve the naviga- tion conditions without aggravating the flood protection,but also ensuring that the ecological system of the river is not damaged,investigations of the morphodynamical processes in connection with artificial grain feeding ac- tivities have to be carried out by Federal Waterways Engineering and Research...

A FAST AND HIGH ACCURACY NUMERICAL SIMULATION FOR A FRACTIONAL BLACK-SCHOLES MODEL ON TWO ASSETS

In this paper,a two dimensional(2D)fractional Black-Scholes(FBS)model on two assets following independent geometric Lévy processes is solved numerically.A high order convergent implicit difference scheme is constructed and detailed numerical analysis is established.The fractional derivative is a quasidifferential operator,whose nonlocal nature yields a dense lower Hessenberg block coefficient matrix.In order to speed up calculation and save storage space,a fast bi-conjugate gradient stabilized(FBi-CGSTAB)method is proposed to solve the resultant linear system.Finally,one example with a known exact solution is provided to assess the effectiveness and efficiency of the presented fast numerical technique.The pricing of a European Call-on-Min option is showed in the other example,in which the influence of fractional derivative order and volatility on the 2D FBS model is revealed by comparing with the classical 2D B-S model.

Modelling and validation of flash flood inundation in drylands

In the context of climate change and human activities,flood disasters in arid mountainous areas have become increasingly frequent,and seriously threatened the safety of people's lives and property.Rapid and accurate flash flood inundation modelling is an essential foundational research area,which can aid in the reduction of casualties and the minimization of disaster losses;however,this modelling is also very difficult,and models need to be urgently developed to address flash flood forecasting and warnings.The objective of this study is to construct a numerical modelling method for flash floods in drylands.Based on a 2D high-resolution flood numerical model(Flood Map-Hydro Inundation2D),we hindcasted the dynamic process of flash flooding and show the spatio-temporal characteristics of flash flood inundation for the“8·18”flash flood disaster that occurred in Datong county,Qinghai province.The results showed that the model output effectively agreed with the observed inundation after the event in terms of both spatial extent and temporal process.Extensive flooding mainly occurred between 00:00 and 01:00 on August 18,2022.Qingshan,Hejiazhuang and Longwo villages were affected most heavily.We further conducted model sensitivity analysis and found that the model was highly sensitive to both roughness and hydraulic conductivity in drylands,and the effect of hydraulic conductivity was more pronounced.Our study confirmed the good performance of our model for the simulation of flash flooding in arid areas and provides a potential method for flash flood assessment and management in arid areas.

A 2-D HYDRODYNAMIC MODEL FOR THE RIVER,LAKE AND NETWORK SYSTEM IN THE JINGJIANG REACH ON THE UNSTRUCTURED QUADRANGLES

The main river, the Dongting Lake and river networks in the Jingjiang reach of the Yangtze River constitute a complex water system, for which a full 2-D hydrodynamic model is established instead of the traditional 1-D or compound models for simulation of such complex systems, based on the latest developments of computer technologies and numerical methods. To better handle irregular boundaries and keep the computation cost well in a reasonable limit, unstructured grids of moderate scale are used. In addition, a dynamic boundary tracking method is proposed to simulate variable flow domains at different floods, especially, when the moderate scale gird can not describe flows in narrow river-network channels at low water levels. The t9 semi-implicit method and the Eulerian-Lagrangian Method （ELM） are adopted, which make the model unconditionally stable with respect to the gravity wave speed and Courant number restrictions. Properties and efficiency of the model are discussed, and it is concluded that the new model is robust and efficient enough for the simulation of a big, complex water system. Validation tests show that the simulation results agree well with field data. It takes about 0.96 h to complete the computation of a 76 d flood, which indicates that the model is efficient enough for engineering applications.

Numerical study of the flow in the Yellow River with non-monotonous banks

A 2-D depth averaged RNG k- ε model is developed to simulate the flow in a typical reach of the Upper Yellow River with non-monotonic banks. In order to take account of the effect of the secondary flow in a bend, the momentum equations are modified by adding an additional source term. A comparison between the numerical simulation and the field measurements indicates that the improved 2-D depth averaged RNG k- ε model can improve the accuracy of the numerical simulation. An arc spline interpolation method is developed to interpolate the non-monotonic river banks. The method can also be reasonably applied for the 2-D interpolation of the river bed level. Through a comparison of the water surface gradients simulated in the seven bends of the studied reach, some analytical formulae are improved to reasonably calculate the longitudinal and transverse gradients in meandering river reaches. Furthermore, the positions of the maximum water depth and the maximum velocity in a typical bend are discussed.