Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

A typhoon-induced storm surge numerical model with GPU acceleration based on an unstructured spherical centroidal Voronoi tessellation grid

Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.

Numerical Models and Methods of Atmospheric Parameters Originating in the Formation of the Earth’s Climatic Cycle

Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in2/s for humidity, 2.993 in2/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .

Numerical modeling of destress blasting for strata separation

Destress blasting(DB)implemented along the perimeter of safety pillars is a special application of destressing in coal longwall mining.The goal is to separate relatively more deformed mined areas from safety pillars,such as shaft pillars or cross-cut pillars,to reduce the transfer of high stresses to the protective pillar.This case study aims to numerically simulate selected destress blasts in the Czech part of the Upper Silesian Coal Basin and examine its impact on stress transfer to the safety pillar area.To separate the area between the protective pillar and the longwall(LW),two fans of five 93-mm blast holes(length of 93e100 m)were drilled from the gate roads into the overburden strata.Each set of blast holes was fired separately in two stages without time delay.The explosive charge(gelatin-type of explosive)of each stage is 3450 kg.The two DB stages were fired when the longwall face was approximately 158 m and 152 m away from the blast.A 3D mine-wide model is built and validated with in situ stress measured with hydrofracturing.Mining and destressing in three 5-m thick coal seams are simulated in the region.Numerical modeling of DB is successfully conducted using a rock fragmentation factor a of 0.05 and a stress reduction/dissipation factor β of 0.95.Buffering of transfer of additional stress from the mining area into the safety pillar is evaluated by comparison of yielding volume before and after DB.It is shown that yielding volume drops after DB by nearly 80%in the area of the destressing panel and near the safety shaft pillar.

Fracture of films caused by uniaxial tensions:a numerical model

Surface cracks are commonly observed in coatings and films.When structures with coatings are subject to stretching,opening mode cracks are likely to form on the surface,which may further lead to other forms of damage,such as interfacial delamination and substrate damage.Possible crack forms include cracks extending towards the interface and channeling across the film.In this paper,a two-dimensional numerical model is proposed to obtain the structural strain energy at arbitrary crack lengths for bilayer structures under uniaxial tension.The energy release rate and structural stress intensity factors can be obtained accordingly,and the effects of geometry and material features on fracture characteristics are investigated,with most crack patterns being confirmed as unstable.The proposed model can also facilitate the analysis of the stress distribution in periodic crack patterns of films.The results from the numerical model are compared with those obtained by the finite element method(FEM),and the accuracy of the theoretical results is demonstrated.

Stiffness and Shear Stress Distribution of Glulam Beams in Elastic-Plastic Stage:Theory,Experiments and Numerical Modelling

Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to predict the stiffness and shear distribution of glulam beams in elastic-plastic stage,and consequently,the failure of such glulam beams cannot be predicted effectively.To address these issues,an analytical method considering material nonlinearity was proposed for glulam beams,and the calculating equations of deflection and shear stress distribution for different failure modes were established.The proposed method was verified by experiments and numerical models under the corresponding conditions.Results showed that the theoretical calculations were in good agreement with experimental and numerical results,indicating that the equations proposed in this paper were reliable and accurate for such glulam beams with wood material in the elastic-plastic stage ignoring the influence of mechanic properties in radial and tangential directions of wood.Furthermore,the experimental results reported by the previous studies indicated that the method was applicable and could be used as a theoretical reference for predicting the failure of glulam beams.

Zoning Evaluation of Hourly Precipitation in High-resolution Regional Numerical Models over Hainan Island

This study assesses the performance of three high-resolution regional numerical models in predicting hourly rainfall over Hainan Island from April to October for the years from 2020 to 2022.The rainfall amount,frequency,intensity,duration,and diurnal cycle are examined through zoning evaluation.The results show that the China Meteor-ological Administration Guangdong Rapid Update Assimilation Numerical Forecast System(CMA-GD)tends to forecast a higher occurrence of light precipitation.It underestimates the late afternoon precipitation and the occurrence of short-duration events.The China Meteorological Administration Shanghai Numerical Forecast Model System(CMA-SH9)reproduces excessive precipitation at a higher frequency and intensity throughout the island.It overestimates rainfall during the late afternoon and midnight periods.The simulated most frequent peak times of rainfall in CMA-SH9 are 0-1 hour deviations from the observed data.The China Meteorological Administration Mesoscale Weather Numerical Forecasting System(CMA-MESO)displays a similar pattern to rainfall observations but fails to replicate reasonable structure and diurnal variation of frequency-intensity.It underestimates the occurrence of long-duration events and overestimates related rainfall amounts from midnight to early morning.Notably,significant discrepancies are observed in the predictions of the three models for areas with complex terrain,such as the central,southeastern,and southwestern regions of Hainan Island.

Flow-Induced Clogging in Microfiltration Membranes: Numerical Modeling and Parametric Study

Microfiltration membrane technology has been widely used in various industries for solid-liquid separation. However, pore clogging remains a persistent challenge. This study employs (CFD) and discrete element method (DEM) models to enhance our understanding of microfiltration membrane clogging. The models were validated by comparing them to experimental data, demonstrating reasonable consistency. Subsequently, a parametric study was conducted on a cross-flow model, exploring the influence of key parameters on clogging. Findings show that clogging is a complex phenomenon affected by various factors. The mean inlet velocity and transmembrane flux were found to directly impact clogging, while the confinement ratio and cosine of the membrane pore entrance angle had an inverse relationship with it. Two clog types were identified: internal (inside the pore) and external (arching at the pore entrance), with the confinement ratio determining the type. This study introduced a dimensionless number as a quantitative clogging indicator based on transmembrane flux, Reynolds number, filtration time, entrance angle cosine, and confinement ratio. While this hypothesis held true in simulations, future studies should explore variations in clogging indicators, and improved modeling of clogging characteristics. Calibration between numerical and physical times and consideration of particle volume fraction will enhance understanding.

Field and numerical modelling on the stability of underground strata in longwall workings

The longwall method is one of the feasible and efficient underground mining methods for greater depth of workings.Recently,Singareni Collieries Company Limited(SCCL),the largest coal-producing government company in the Southern part of India,has deployed a high capacity(2×1152T capacity)power support system in its Adriyala Longwall Project(ALP)at a depth of 375 m.There was a concern about the stability of the longwall workings and the protective pillars.Thus,research work was carried out with the geotechnical instruments and numerical modelling tools to analyze the stability while retreating the longwall panels.The outcomes revealed that the convergence in the gate roads increased with the longwall face advancement and the area of exposure.The pressure of the powered support legs on the dip side was less than the rise side legs,which implies a stable roof condition over the longwall face.An abutment zone was identified ahead of the line of extraction up to 10e25 m and presumed from the Factor of Safety(FoS)criterion that within this zone,the longwall pillar possibly will be of better stability at every stage of extraction.The results obtained from this study would be helpful for the mining engineers to understand the behaviour of underground strata in longwall workings.

A Heavy Sea Fog Event over the Yellow Sea in March 2005: Analysis and Numerical Modeling

In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility along the Shandong Peninsula coast to 100 m or much less at some sites. Satellite images, surface observations and soundings at islands and coasts, and analyses from the Japan Meteorology Agency （JMA） axe used to describe and analyze this event. The analysis indicates that this sea fog can be categorized as advection cooling fog. The main features of this sea fog including fog area and its movement axe reasonably reproduced by the Fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model （MM5）. Model results suggest that the formation and evolution of this event can be outlined as： （1） southerly warm/moist advection of low-level air resulted in a strong sea-surface-based inversion with a thickness of about 600 m; （2） when the inversion moved from the warmer East Sea to the colder Yellow Sea, a thermal internal boundary layer （TIBL） gradually formed at the base of the inversion while the sea fog grew in response to cooling and moistening by turbulence mixing; （3） the sea fog developed as the TIBL moved northward and （4） strong northerly cold and dry wind destroyed the TIBL and dissipated the sea fog. The principal findings of this study axe that sea fog forms in response to relatively persistent southerly waxm/moist wind and a cold sea surface, and that turbulence mixing by wind shear is the primary mechanism for the cooling and moistening the marine layer. In addition, the study of sensitivity experiments indicates that deterministic numerical modeling offers a promising approach to the prediction of sea fog over the Yellow Sea but it may be more efficient to consider ensemble numerical modeling because of the extreme sensitivity to model input.

Numerical modelling of structural controls on fluid flow and mineralization

This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization. We first give an overview of natural phenomena of structurally-controlled ore formation and the background theory and mechanisms for such controls. We then provide the results of a group of simple 2D numerical models validated through comparison with Cu-vein structure observed near the Shilu Copper deposit （Yangchun, Guangdong Province, China） and finally a case study of 3D numerical modelling applied to the Hodgkinson Province in North Queensland （Australia）. Two modelling approaches, discrete deformation modelling and continuum coupled deformation and fluid flow modelling, are involved. The 2D model-derived patterns are remarkably consistent with the Cu-vein structure from the Shilu Copper deposit, and show that both modelling approaches can realistically simulate the mechanical behaviours of shear and dilatant fractures. The continuum coupled deformation and fluid flow model indicates that pattern of the Cu- veins near the Shilu deposit is the result of shear strain localization, development of dilation and fluid focussing into the dilatant fracture segments. The 3D case-study models （with deformation and fluid flow coupling） on the Hodgkinson Province generated a number of potential gold mineralization

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced, with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.

Numerical Modeling on Suspended Sediment Transportation in the Hangzhou Bay

In this paper the characteristics of tidal flow and seasonal variation of seidment content in the Hangzhou Bay and their affecting factors are studied. Field investigations and data analysis indicate that the sediment movement is mainly influenced by the Yangtze estuary and the sediment of the Yangtze estuary is induced by wind wave and tidal flow. Owing to the variation of dynamic conditions, the instantaneous sediment content is controlled by tidal flow, wind wave, depth of water and tidal range synthetically. A sediment content relationship formula is established with related factors. A non-equilibrium 2-dimensional numerical model of suspended sediment transportation is set up, and the finite element method is applied. The computation results of the model is in accordance with field data.

Three Dimensional Baroclinic Numerical Model for Simulating Fresh and Salt Water Mixing in the Yangtze Estuary

For simulating fresh and salt water mixing in estuaries, a three dimensional nonlinear baroclinic numerical model is developed, in which the gradients of horizontal pressure contain die gradient of barotropic pressure arising from the gradient of tidal level and the gradient of baroclinic pressure due to the gradient of salinity. The Eulerian-Lagrangian method is employed to descretize both the momentum equations of tidal motion and the equation of salt water diffusion so as to improve the computational stability and accuracy. The methods to provide the boundary conditions and the initial conditions are proposed, and the criterion for computational stability of the salinity fields is presented. The present model is used for modeling fresh and salt water mixing in the Yangtze Estuary. Computations show that the salinity distribution has the characteristics of partial mixing pattern, and that the present model is suitable for simulation of fresh and salt water mixing in the Yangtze Estuary.

Examining the effect of adverse geological conditions on jamming of a single shielded TBM in Uluabat tunnel using numerical modeling

Severe shield jamming events have been reported during excavation of Uluabat tunnel through adverse geological conditions, which resulted in several stoppages at advancing a single shielded tunnel boring machine(TBM). To study the jamming mechanism, three-dimensional(3D) simulation of the machine and surrounding ground was implemented using the finite difference code FLAC3D. Numerical analyses were performed for three sections along the tunnel with a higher risk for entrapment due to the combination of overburden and geological conditions. The computational results including longitudinal displacement contours and ground pressure profiles around the shield allow a better understanding of ground behavior within the excavation. Furthermore, they allow realistically assessing the impact of adverse geological conditions on shield jamming. The calculated thrust forces, which are required to move the machine forward, are in good agreement with field observations and measurements. It also proves that the numerical analysis can effectively be used for evaluating the effect of adverse geological environment on TBM entrapments and can be applied to prediction of loads on the shield and preestimating of the required thrust force during excavation through adverse ground conditions.

Overhanging rock slope by design:An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.

Numerical Models of Higher-Order Boussinesq Equations and Comparisons with Laboratory Measurement

Nonlinear water wave propagation passing a submerged shelf is studied experimentally and numerically. The applicability of two different wave propagation models has been investigated. One is higher-order Boussinesq equations derived by Zou (1999) and the other is the classic Boussinesq equations, Physical experiments are conducted, three different front slopes (1:10, 1:5 and 1:2) of the shelf are set up in the experiment and their effects on wave propagation are investigated. Comparisons of numerical results with test data are made, the model of higher-order Boussinesq equations agrees much better with the measurements than the model of the classical Boussinesq equations, The results show that the higher-order Boussinesq equations can also be applied to the steeper slope case although the mild slope assumption is employed in the derivation of the higher order terms of higher order Boussinesq equations.

Application of numerical modeling and genetic programming to estimate rock mass modulus of deformation

Estimation of the rock mass modulus of deformation(Em)is one of the most important design parameters in designing many structures in and on rock.This parameter can be obtained by in situ tests,empirical relations between deformation modulus and rock mass classifcation,and estimating from laboratory tests results.In this paper,a back analysis calculation is performed to present an equation for estimation of the rock mass modulus of deformation using genetic programming(GP)and numerical modeling.A database of 40,960 datasets,including vertical stress(rz),horizontal to vertical stresses ratio(k),Poisson’s ratio(m),radius of circular tunnel(r)and wall displacement of circular tunnel on the horizontal diameter(d)for input parameters and modulus of deformation for output,was established.The selected parameters are easy to determine and rock mass modulus of deformation can be obtained from instrumentation data of any size circular galleries.The resulting RMSE of 0.86 and correlation coeffcient of97%of the proposed equation demonstrated the capability of the computer program(CP)generated by GP.

Numerical Model of Internal Solitary Wave Evolution on ImpermeableVariable Seabed in A Stratified Two-Layer Fluid System

In the present study a numerical model developed by Lynett and Liu （2002） is modified to include density difference in a stratified two-layer fluid in a three-dimensional internal wave domain. The internal solitary wave （ISW） in the model is assumed to be weakly nonlinear and weakly dispersive, and the viscosity effects at all boundaries are ignored. The governing equations based on the Navier-Stokes and Euler equations are solved for internal solitary wave propagation over variable seabed topography. Theoretical formulations are established, from which analytical solutions are obtained, in addition to numerical results. Wave profiles from previous experimental studies are compared with the numerical results from the present analytical solutions. Numerical models developed on the basis of the present analytical solutions are better than those developed by Lynett and Liu （2002）. The results of numerical modeling agree well with the experimental data.