Numerical investigation of friction-heating-pressurization and its control parameters in the shear band of high-speed landslides

High-speed sliding often leads to catastrophic landslides,many of which,in the initial sliding phase before disintegration,experience a friction-induced thermal pressurization effect in the bottom shear band,accelerating the movement of the overlying sliding mass.To quantitatively investigate this complex multiphysical phenomenon,we established a set of equations that describe the variations in temperature and excess pore pressure within the shear band,as well as the conservation of momentum equation for the overlying sliding mass.With a simplified landslide model,we investigated the variations of temperature and excess pore pressure within the shear band and their impacts on the velocity of the overlying sliding mass.On this basis,we studied the impact of seven key parameters on the maximum temperature and excess pore pressure in the shear band,as well as the impact on the velocity of the overlying sliding mass.The simulation results of the standard model show that the temperature and excess pore pressure in the shear band are significantly higher than those in the adjacent areas,and reach the maximum values in the center.Within a few seconds after the start,the maximum excess pore pressure in the shear zone is close to the initial stress,and the shear strength loss rate exceeds 90%.The thermal pressurization mechanism significantly increases the velocity of the overlying sliding mass.The results of parameter sensitivity analysis show that the thermal expansion coefficient has the most significant impact on the temperature and excess pore pressure in the shear band,and the sliding surface dip angle has the most significant impact on the velocity of the overlying sliding mass.The results of this study are of great significance for clarifying the mechanism of thermal pressurization-induced high-speed sliding.

Numerical Investigation on the Harbour Resonance Inside A Harbour with Lateral Cavities

Wave-induced harbour resonance is numerically investigated inside a harbour with lateral cavities.The theoretical solutions for the amplification parameter are compared with the simulated results under varying dimensionless wave numbers in order to verify the simulation model in a rectangular harbour at a constant depth.The results indicate that the numerical model can correctly calculate the natural frequency and the natural wave height.A range of calculations are performed for harbour resonance with one pair of lateral cavities,two pairs of lateral cavities and three pairs of lateral cavities,respectively.The simulated results indicate that the amplitude of the amplification parameter decreases both at the primary natural oscillation and the secondary natural oscillation,as the number of lateral cavities increases.The dimensionless wave number reduces as the number of lateral cavities increases both at the primary natural oscillation and the secondary natural oscillation as well.

Numerical investigation on the aerodynamic drag reduction based on bottom deflectors and streamlined bogies of a high-speed train

The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-ε turbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.

Elastic-viscoplastic behaviors of polymer-blend geocell sheets:Numerical and experimental investigations

Polymer-blend geocell sheets(PBGS)have been developed as substitute materials for manufacturing geocells.Various attempts have been made to test and predict the behaviors of commonly used geogrids,geotextiles,geomembranes,and geocells.However,the elastic-viscoplastic behaviors of novel-developed geocell sheets are still poorly understood.Therefore,this paper investigates the elastic-viscoplastic behaviors of PBGS to gain a comprehensive understanding of their mechanical properties.Furthermore,the tensile load-strain history under various loading conditions is simulated by numerical calculation for widespread utilization.To achieve this goal,monotonic loading tests,short-term creep and stress relaxation tests,and multi-load-path tests(also known as arbitrary loading history tests)are performed using a universal testing machine.The results are simulated using the nonlinear three-component(NLTC)model,which consists of three nonlinear components,i.e.a hypo-elastic component,a nonlinear inviscid component,and a nonlinear viscid component.The experimental and numerical results demonstrate that PBGS exhibit significant elastic-viscoplastic behavior that can be accurately predicted by the NLTC model.Moreover,the tensile strain rates significantly influence the tensile load,with higher strain rates resulting in increased tensile loads and more linear load-strain curves.Also,parametric analysis of the rheological characteristics reveals that the initial tensile strain rates have negligible impact on the results.The rate-sensitivity coefficient of PBGS is approximately 0.163,which falls within the typical range observed in most geosynthetics.

A Numerical Investigation of the Effect of Boundary Conditions on Acoustic Pressure Distribution in a Sonochemical Reactor Chamber

The intensification of physicochemical processes in the sonochemical reactor chamber is widely used in problems of synthesis,extraction and separation.One of the most important mechanisms at play in such processes is the acoustic cavitation due to the non-uniform distribution of acoustic pressure in the chamber.Cavitation has a strong impact on the surface degradation mechanisms.In this work,a numerical calculation of the acoustic pressure distribution inside the reactor chamber was performed using COMSOL Multiphysics.The numerical results have revealed the dependence of the structure of the acoustic pressure field on the boundary conditions for various thicknesses of the piezoelectric transducer.In particular,the amplitude of the acoustic pressure is minimal in the case of absorbing boundaries,and the attenuation becomes more significant as the thickness of the piezoelectric transducer increases.In addition,reflective boundaries play a significant role in the formation and distribution of zones of maximum cavitation activity.

Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Investigations on High-Speed Flash Boiling Atomization of Fuel Based on Numerical Simulations

Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.

Multiparameter Numerical Investigation of Two Types of Moving Interactions Between the Deep-Sea Mining Vehicle Track Plate and Seabed Soil:Digging and Rotating Motions

To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.

Numerical investigation of geostress influence on the grouting reinforcement effectiveness of tunnel surrounding rock mass in fault fracture zones

Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.

Numerical Investigation of Thermal Behavior of CNC Machine Tool and Its Effects on Dimensional Accuracy of Machined Parts

The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.

Experimental and numerical investigation on alternatives to sandy gravel

The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbelly improvised explosive device(IEDs),which must be buried in water-saturated sandy gravel.The use of sandy gravel has some drawbacks,for instance reproducibility,time consumption,and cost.This paper focuses on the investigation of four alternatives to sandy gravel,which could produce similar specific and cumulative impulses:a concrete pot filled with water,a concrete pot filled with quartz sand,a steel pot without filling and a concrete pot filled with glass spheres(diameter 200μm—300μm)and different water contents.The impulses are measured with a ring technology developed at the Fraunhofer EMI.A numerical soil model based on the work of Marrs,2014 and Fi serov a,2006 and considering the soil moisture was used to simulate the experiments with glass spheres at different water contents,showing much better agreement with the experiments than the classical Laine&Sandvik model,even for high saturation levels.These results can be used to create new test conditions at original scale that are more cost-effective,more reproducible and simpler to manage in comparison to the current tests carried out with STANAG sandy gravel.

Numerical investigation of hydro-morphodynamic characteristics of a cascading failure of landslide dams

A cascading failure of landslide dams caused by strong earthquakes or torrential rains in mountainous river valleys can pose great threats to people’s lives,properties,and infrastructures.In this study,based on the three-dimensional Reynoldsaveraged Navier-Stokes equations(RANS),the renormalization group(RNG)k-εturbulence model,suspended and bed load transport equations,and the instability discriminant formula of dam breach side slope,and the explicit finite volume method(FVM),a detailed numerical simulation model for calculating the hydro-morphodynamic characteristics of cascading dam breach process has been developed.The developed numerical model can simulate the breach hydrograph and the dam breach morphology evolution during the cascading failure process of landslide dams.A model test of the breaches of two cascading landslide dams has been used as the validation case.The comparison of the calculated and measured results indicates that the breach hydrograph and the breach morphology evolution process of the upstream and downstream dams are generally consistent with each other,and the relative errors of the key breaching parameters,i.e.,the peak breach flow and the time to peak of each dam,are less than±5%.Further,the comparison of the breach hydrographs of the upstream and downstream dams shows that there is an amplification effect of the breach flood on the cascading landslide dam failures.Three key parameters,i.e.,the distance between the upstream and the downstream dams,the river channel slope,and the downstream dam height,have been used to study the flood amplification effect.The parameter sensitivity analyses show that the peak breach flow at the downstream dam decreases with increasing distance between the upstream and the downstream dams,and the downstream dam height.Further,the peak breach flow at the downstream dam first increases and then decreases with steepening of the river channel slope.When the flood caused by the upstream dam failure flows to the downstream dam,it can produce a surge wave that overtops and erodes the dam crest,resulting in a lowering of the dam crest elevation.This has an impact on the failure occurrence time and the peak breach flow of the downstream dam.The influence of the surge wave on the downstream dam failure process is related to the volume of water that overtops the dam crest and the erosion characteristics of dam material.Moreover,the cascading failure case of the Xiaogangjian and Lower Xiaogangjian landslide dams has also been used as the representative case for validating the model.In comparisons of the calculated and measured breach hydrographs and final breach morphologies,the relative errors of the key dam breaching parameters are all within±10%,which verify the rationality of the model is applicable to real-world cases.Overall,the numerical model developed in this study can provide important technical support for the risk assessment and emergency treatment of failures of cascading landslide dams.

Numerical investigation of the mechanical behavior of the backfill–rock composite structure under triaxial compression

To ensure safe and economical backfill mining,the mechanical response of the backfill–rock interaction system needs to be understood.The numerical investigation of the mechanical behavior of backfill–rock composite structure(BRCS)under triaxial compression,which includes deformation,failure patterns,strength characteristics,and acoustic emission(AE)evolution,was proposed.The models used in the tests have one rough interface,two cement–iron tailings ratios(CTRs),four interface angles(IAs),and three confining pressures(CPs).Results showed that the deformation,strength characteristics,and failure patterns of BRCS under triaxial compression depend on IA,CP,and CTR.The stress–strain curves of BRCS under triaxial compression could be divided into five stages,namely,compaction,elasticity,yield,strain softening,and residual stress.The relevant AE counts have corresponding relationships with different stages.The triaxial compressive strengths of composites increase linearly with the increase of the CP.Furthermore,the CP stress strengthening effect occurs.When the IAs are45°and 60°,the failure areas of composites appear in the interface and backfill.When the IAs are 75°and 90°,the failure areas of composites appear in the backfill,interface,and rock.Moreover,the corresponding failure modes yield the combined shear failure.The research results provide the basis for further understanding of the stability of the BRCS.

Experimental and numerical investigation on failure characteristics and mechanism of coal with different water contents

Water injection, as a widely used technique to prevent coal burst, can restrain the fractured coal seam and released the energy storage. In this study, laboratory tests were frstly carried out on standard coal specimens with fve diferent water contents (i.e., 0%, 0.6%, 1.08%, 1.5%, 2.0%, and 2.3%). The failure mode, fragment size, and energy distribution characteristics of coal specimens were investigated. Experimental results show that strength, elastic strain energy, dissipated energy, brittleness index, as well as impact energy index decrease with increasing water content. Besides, the failure mode transitions gradually from splitting ejection to tensile-shear mixed failure mode as water content increases, and average fragment size shows positively related to water content. Moreover, scanning electron microscope tests results indicate that water in coal sample mainly causes the mineral softening and defects increase. Furthermore, a numerical model containing roadway excavation was established considering the water on coal burst prevention. Modelling results revealed that water injection can reduce degree of coal burst and ejection velocity of coal blocks, while it will raise up the depth of crack zone and surface displacement of roadway. Combined with laboratory tests and numerical results, the micro mechanism, energy mechanism, and engineering signifcance of water injection on coal burst prevention were fnally analyzed.

Numerical and theoretical investigations of heat transfer characteristics in helium-xenon cooled microreactor core

Helium-xenon cooled microreactors are a vital technological solution for portable nuclear reactor power sources.To exam-ine the convective heat transfer behavior of helium-xenon gas mixtures in a core environment,numerical simulations are conducted on a cylindrical coolant channel and its surrounding solid regions.Validated numerical methods are used to determine the effect and mechanisms of power and its distribution,inlet temperature and velocity,and outlet pressure on the distribution and change trend of the axial Nusselt number.Furthermore,a theoretical framework that can describe the effect of power variation on the evolution of the thermal boundary layer is employed to formulate an axial distribution cor-relation for the Nusselt number of the coolant channel,under the assumption of a cosine distribution for the axial power.Based on the simulation results,the correlation coefficients are determined,and a semi-empirical relationship is identified under the corresponding operating conditions.The correlation derived in this study is consistent with the simulations,with an average relative error of 5.3%under the operating conditions.Finally,to improve the accuracy of the predictions near the entrance,a segmented correlation is developed by combining the Kays correlation with the aforementioned correlation.The new correlation reduces the average relative error to 2.9%and maintains satisfactory accuracy throughout the entire axial range of the channel,thereby demonstrating its applicability to turbulent heat transfer calculations for helium-xenon gas mixtures within the core environment.These findings provide valuable insights into the convective heat transfer behavior of a helium-xenon gas mixture in a core environment.

Experimental investigation on coal pore-fracture variation and fractal characteristics synergistically affected by solvents for improving clean gas extraction

Chemical solvents instead of pure water being as hydraulic fracturing fluid could effectively increase permeability and improve clean methane extraction efficiency.However,pore-fracture variation features of lean coal synergistically affected by solvents have not been fully understood.Ultrasonic testing,nuclear magnetic resonance analysis,liquid phase mass spectrometry was adopted to comprehensively analyze pore-fracture change characteristics of lean coal treated by combined solvent(NMP and CS_(2)).Meanwhile,quantitative characterization of above changing properties was conducted using geometric fractal theory.Relationship model between permeability,fractal dimension and porosity were established.Results indicate that the end face fractures of coal are well developed after CS2and combined solvent treatments,of which,end face box-counting fractal dimensions range from 1.1227 to 1.4767.Maximum decreases in ultrasonic longitudinal wave velocity of coal affected by NMP,CS_(2)and combined solvent are 2.700%,20.521%,22.454%,respectively.Solvent treatments could lead to increasing amount of both mesopores and macropores.Decrease ratio of fractal dimension Dsis 0.259%–2.159%,while permeability increases ratio of NMR ranges from 0.1904 to 6.4486.Meanwhile,combined solvent could dissolve coal polar and non-polar small molecules and expand flow space.Results could provide reference for solvent selection and parameter optimization of permeability-enhancement technology.

Recent progresses in the development of tubular segmented-in-series solid oxide fuel cells:Experimental and numerical study

Solid oxide fuel cells(SOFCs)have attracted a great deal of interest because they have the highest efficiency without using any noble metal as catalysts among all the fuel cell technologies.However,traditional SOFCs suffer from having a higher volume,current leakage,complex connections,and difficulty in gas sealing.To solve these problems,Rolls-Royce has fabricated a simple design by stacking cells in series on an insulating porous support,resulting in the tubular segmented-in-series solid oxide fuel cells(SIS-SOFCs),which achieved higher output voltage.This work systematically reviews recent advances in the structures,preparation methods,perform-ances,and stability of tubular SIS-SOFCs in experimental and numerical studies.Finally,the challenges and future development of tubular SIS-SOFCs are also discussed.The findings of this work can help guide the direction and inspire innovation of future development in this field.

Numerical Investigation on Inclined Hydraulic Transport of Large Particles for Deepsea Mining

Hydraulic transport in pipelines is the most promising conveying method for large ore particles in deepsea mining.The dynamic performances of particles during transportation in vertical,inclined and horizontal pipelines are significant for the design of hydraulic transport systems.In the present study,we focus on the statistical characteristics and flow regimes of the mixture composed of ore particles and seawater in the pipelines.Numerical simulations are conducted by using Computational Fluid Dynamics(CFD)and Discrete Element Method(DEM).The influences of inclination angle and particle diameter are evaluated through two sets of numerical tests.The regulation of the inclined transport is totally different from that of the vertical transport,whereas the dynamics of the mixtures in inclined and horizontal pipes are similar.A number of particles accumulate on the pipe wall even with a small inclination angle.Large hydraulic gradient and local concentration would occur when the inclination angle of the pipe is in the range of30°-60°.With the decrease of particle diameter,the particle flow becomes uniform,reflected by the almost uniform particle distribution in the vertical pipe and the clear interface between the suspended load and the bed-load in the inclined pipe.However,small particles will introduce larger local concentrations and hydraulic gradients in the inclined pipe,which is not conducive to particle transport.

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.

Numerical study on local scour characteristics around submarine pipelines in the Yellow River Delta silty sandy soil under waves and currents

Due to their high reliability and cost-efficiency,submarine pipelines are widely used in offshore oil and gas resource engineering.Due to the interaction of waves,currents,seabed,and pipeline structures,the soil around submarine pipelines is prone to local scour,severely affecting their operational safety.With the Yellow River Delta as the research area and based on the renormalized group(RNG)k-εturbulence model and Stokes fifth-order wave theory,this study solves the Navier-Stokes(N-S)equation using the finite difference method.The volume of fluid(VOF)method is used to describe the fluid-free surface,and a threedimensional numerical model of currents and waves-submarine pipeline-silty sandy seabed is established.The rationality of the numerical model is verified using a self-built waveflow flume.On this basis,in this study,the local scour development and characteristics of submarine pipelines in the Yellow River Delta silty sandy seabed in the prototype environment are explored and the influence of the presence of pipelines on hydrodynamic features such as surrounding flow field,shear stress,and turbulence intensity is analyzed.The results indicate that(1)local scour around submarine pipelines can be divided into three stages:rapid scour,slow scour,and stable scour.The maximum scour depth occurs directly below the pipeline,and the shape of the scour pits is asymmetric.(2)As the water depth decreases and the pipeline suspension height increases,the scour becomes more intense.(3)When currents go through a pipeline,a clear stagnation point is formed in front of the pipeline,and the flow velocity is positively correlated with the depth of scour.This study can provide a valuable reference for the protection of submarine pipelines in this area.