In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over ...In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.展开更多
A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulat...A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulate the plasma arc. The argon-nitrogen plasma arc is simulated for different arc currents and gas flow rates of argon. Various combinations of arc core radius and arc length, which correspond to a given torch power, are predicted. A most feasible combination of the same, which corresponds to an actual physical situation of the arc inside the torch, is identified using the thermodynamic principle of minimum entropy production for a particular torch power. The effect of the arc current and gas flow rate on the plasma arc characteristics and torch efficiency is explained. The effect of the nitrogen content in the plasma gas on the torch power and efficiency is clearly detected. Predicted torch efficiencies are comparable to the measured ones and the effect of the arc current and gas flow rate on predicted and measured efficiencies is almost similar. The efficiency of the torch, cathode and anode losses and core temperature and velocity at the nozzle exit are reported for five different cases.展开更多
Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock st...Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation.This may lead to problems in the evaluation of rock structure stability and safe life.Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross-scale view.This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures,that is,the homogenization theory,the hierarchical approach,and the concurrent approach.For these approaches,their benefits,drawbacks,and application scope are underlined.Despite the considerable attempts that have been made,some key issues still result in multiple challenges.Therefore,this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future.The review results show that,in addition to numerical techniques,for example,high-performance computing,more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.展开更多
Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implemen...Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.展开更多
To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets ...To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.展开更多
Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to en...Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.展开更多
Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ra...Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ratios of merging water bodies.This study investigated the mixing structure at open channel confluences using three-dimensional numerical modeling.A comprehensive three-dimensional numerical model was developed and validated against a dataset obtained from controlled laboratory experiments.This dataset incorporated three-dimensional time-averaged velocity measurements.The skew-induced and stress-induced equation systems were adopted as the core governing equations,providing a framework for simulating various scenarios.A total of ten different cases were analyzed.The results highlighted the effect of discharge ratios on turbulence,lateral and vertical vorticities,and the distribution of mixing,which intensified with higher magnitudes of discharge ratios.The mixing structure,driven by velocity gradients and vorticity,revealed the significant role of lateral and vertical vorticities in determining hydrodynamic behaviors and mixing distributions at confluences.Specifically,the momentum ratio of incoming flows governed the spatial evolution of mixing processes.This study revealed that the distribution of mixing served as a key indicator for identifying the formation of mid-channel scours.High normalized velocities induced toward the left bank led to the superelevation of the water surface,enhancing the potential for bed material and the formation of significant scour holes beneath the elevated water surface.This novel approach provides a deeper understanding of the mixing patterns at confluences,particularly in scenarios with equilibrated discharge ratios but in different magnitudes.展开更多
Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade compone...Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.展开更多
A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of ...A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.展开更多
Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels ha...Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels has an important bearing on sediment transport, riverbed evolution, and pollutant transport study. In this article, different turbulence closure schemes i.e., the mixing-length model and the k-ε model with different pressure solution techniques i. e., hydrostatic assumptions and dynamic pressure treatments are applied to study the helical secondary flows in an experiment curved channel. The agreements of vertically-averaged velocities between the simulated results obtained by using different turbulence models with different pressure solution techniques and the measured data are satisfactory. Their discrepancies with respect to surface elevations, superelevations and secondary flow patterns are discussed.展开更多
This paper summarizes the development of a three-dimensional numerical model for analyzing single geocell-reinforced soil.In this model,the infill soil was modeled using the Duncan-Chang model,which can simulate non-l...This paper summarizes the development of a three-dimensional numerical model for analyzing single geocell-reinforced soil.In this model,the infill soil was modeled using the Duncan-Chang model,which can simulate non-linearity and stress-dependency of soil.Geocell was modeled using linearly elastic plate elements,which can carry both bending and membrane stresses.A linear interface stress-strain relationship with a MohrCoulomb yield criterion was adopted to model the interface friction between the geocell wall and the soil.By modeling the geocell and the soil separately,the interaction between the soil and the geocell can be accurately simulated.To verify this model,a plate load test was conducted in the laboratory,in which a 12-cmthick sand layer reinforced by a single geocell was subjected to a vertical load from a circular steel plate.The load-displacement curves and the horizontal tensile strain of the geocell were recorded during the test.A numerical model was created according to the setup of the load test.The numerical results compared reasonably well with the test data.展开更多
With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rock...With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rockbursts.To investigate the influence of tunnel diameter on the deformation and failure characteristics of surrounding rock,large-sized rocklike gypsum specimens were tested using a selfdeveloped true triaxial rockburst loading system containing circular tunnels with three different diameters(D=0.07 m,0.11 m,and 0.15 m).Acoustic emission monitoring,together with a miniature intelligent camera,was employed to analyze the entire process,focusing on macroscopic failure patterns,fragment characteristics,and underlying failure mechanisms.In addition,theoretical analyses were carried out and combined with numerical simulations to investigate the differences in energy evolution associated with rockburst physical models.The results indicate that:(1)The rockburst process with different tunnel diameters consistently evolved through three distinct stages—initial particle ejection,crack propagation accompanied by flake spalling,and,finally,fragment ejection leading to the formation of a‘V'-shaped notch.(2)Increasing tunnel diameter reduces rockburst failure load while increasing surrounding rock damage extent,total mass and average size of ejected fragments.Additionally,shear failure proportion decreases with tensile failure becoming increasingly dominant.(3)Larger tunnel diameters reduce the attenuation rate of elastic strain energy,thereby expanding the zone of elastic strain energy accumulation and disturbance and creating conditions for larger volume rockburst.(4)Larger tunnel diameters result in a smaller principal stress ratio at equivalent distances in the surrounding rock,indicating a higher likelihood of tensile failure.(5)Numerical analyses further reveal that larger tunnel diameters reduce the maximum elastic strain energy density around the tunnel,lowering the energy released per unit volume of rockburst fragments and their ejection velocities.However,both the total failure volume and overall energy release from rockburst increase.Model experiments with different tunnel diameters are of great significance for optimizing engineering design and parameter selection,as well as guiding tunnel construction under complex geological conditions.展开更多
The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbule...The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbulence models struggle to make accurate predictions for subsonic and supersonic flows in nozzles.In this study,we explored a novel model,the algebraic stress model k-kL-ARSM+J,to enhance the accuracy of turbulence numerical simulations.This new model was used to conduct numerical simulations of the design and off-design performance of a 3D supersonic asymmetric truncated nozzle designed in our laboratory,with the aim of providing a realistic pattern of changes.The research indicates that,compared to linear eddy viscosity turbulence models such as k-kL and shear stress transport(SST),the k-kL-ARSM+J algebraic stress model shows better accuracy in predicting the performance of supersonic nozzles.Its predictions were identical to the experimental values,enabling precise calculations of the nozzle.The performance trends of the nozzle are as follows:as the inlet Mach number increases,both thrust and pitching moment increase,but the rate of increase slows down.Lift peaks near the design Mach number and then rapidly decreases.With increasing inlet pressure,the nozzle thrust,lift,and pitching moment all show linear growth.As the flight altitude rises,the internal flow field within the nozzle remains relatively consistent due to the same supersonic nozzle inlet flow conditions.However,external to the nozzle,the change in external flow pressure results in the nozzle exit transitioning from over-expanded to under-expanded,leading to a shear layer behind the nozzle that initially converges towards the nozzle center and then diverges.展开更多
Systematic bias is a type of model error that can affect the accuracy of data assimilation and forecasting that must be addressed.An online bias correction scheme called the sequential bias correction scheme(SBCS),was...Systematic bias is a type of model error that can affect the accuracy of data assimilation and forecasting that must be addressed.An online bias correction scheme called the sequential bias correction scheme(SBCS),was developed using the6 h average bias to correct the systematic bias during model integration.The primary purpose of this study is to investigate the impact of the SBCS in the high-resolution China Meteorological Administration Meso-scale(CMA-MESO)numerical weather prediction(NWP)model to reduce the systematic bias and to improve the data assimilation and forecast results through this method.The SBCS is improved upon and applied to the CMA-MESO 3-km model in this study.Four-week sequential data assimilation and forecast experiments,driven by rapid update and cycling(RUC),were conducted for the period from 2–29 May 2022.In terms of the characteristics of systematic bias,both the background and analysis show diurnal bias,and these large biases are affected by complex underlying surfaces(e.g.,oceans,coasts,and mountains).After the application of the SBCS,the results of the data assimilation show that the SBCS can reduce the systematic bias of the background and yield a neutral to slightly positive result for the analysis fields.In addition,the SBCS can reduce forecast errors and improve forecast results,especially for surface variables.The above results indicate that this scheme has good prospects for high-resolution regional NWP models.展开更多
Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains in...Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains insufficiently optimized,demanding more profound mechanistic studies and engineering refinements.To address the difficulties in resolving the flow inhomogeneity at the stack scale,this study establishes a multi-physics field coupling model and analyzes the pressure distributions,flow rate differences,active substance concentration,and electrochemical characteristics.The results show that the uneven cell pressure distribution is a key factor affecting the consistency of the system performance,and the increase in the flow rate improves the reactant homogeneity,with both the average concentration and the uniformity factor increasing with the flow rate.In contrast,high current densities lead to an increased imbalance between electrochemical depletion and reactant replenishment,resulting in a significant decrease in reactant concentration in the under-ribs region.In addition,a higher flow rate can expand the high-current-density region where the stack operates efficiently.This study provides a theoretical basis for optimizing the design of the stack components.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques desig...Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques designed for tertiary oil recovery have garnered significant attention,with microgel flooding emerging as a particularly prominent area of research.Despite its promise,the complex mechanisms underlying microgel flooding have been rarely investigated numerically.This study aims to address these gaps by characterizing the distribution of microgel concentration and viscosity within different pore structures.To enhance the accuracy of these characterizations,the viscosity of microgels is adjusted to account for the shear effects induced by flow rate and the swelling effects driven by salinity variations.The absolute permeability of the rock and the relative permeability of both oil and microgel are also analyzed to elucidate the mechanisms of microgel flooding.Additionally,a connectivity model is employed to achieve a quantitative representation of fluid flow capacity.The proposed model is validated through conceptual examples and applied to real oilfield blocks,demonstrating its accuracy and practical applicability.展开更多
The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope ...The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.展开更多
The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying ge...The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.展开更多
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...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.展开更多
基金supported by the Natural Science Foundation of China(No.41574127)the China Postdoctoral Science Foundation(No.2017M622608)the project for the independent exploration of graduate students at Central South University(No.2017zzts008)
文摘In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.
文摘A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulate the plasma arc. The argon-nitrogen plasma arc is simulated for different arc currents and gas flow rates of argon. Various combinations of arc core radius and arc length, which correspond to a given torch power, are predicted. A most feasible combination of the same, which corresponds to an actual physical situation of the arc inside the torch, is identified using the thermodynamic principle of minimum entropy production for a particular torch power. The effect of the arc current and gas flow rate on the plasma arc characteristics and torch efficiency is explained. The effect of the nitrogen content in the plasma gas on the torch power and efficiency is clearly detected. Predicted torch efficiencies are comparable to the measured ones and the effect of the arc current and gas flow rate on predicted and measured efficiencies is almost similar. The efficiency of the torch, cathode and anode losses and core temperature and velocity at the nozzle exit are reported for five different cases.
基金National Natural Science Foundation of China,Grant/Award Numbers:52192691,52192690。
文摘Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation.This may lead to problems in the evaluation of rock structure stability and safe life.Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross-scale view.This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures,that is,the homogenization theory,the hierarchical approach,and the concurrent approach.For these approaches,their benefits,drawbacks,and application scope are underlined.Despite the considerable attempts that have been made,some key issues still result in multiple challenges.Therefore,this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future.The review results show that,in addition to numerical techniques,for example,high-performance computing,more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.
基金supported by grants received by the first author and third author from the Institute of Eminence,Delhi University,Delhi,India,as part of the Faculty Research Program via Ref.No./IoE/2024-25/12/FRP.
文摘Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.
基金financially supported by the Ministry of Science and Technology of China(Nos.2022YFF0801201,2021YFC2900300)the National Natural Science Foundation of China(Nos.41872245,U1911202)the Guangdong Basic and Applied Basic Research Foundation(No.2020A1515010666)。
文摘To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.
基金supported by the National Natural Science Foundation of China(No.92371206)the Postgraduate Scientific Research Innovation Project of Hunan Province,China(No.CX2023063).
文摘Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.
文摘Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ratios of merging water bodies.This study investigated the mixing structure at open channel confluences using three-dimensional numerical modeling.A comprehensive three-dimensional numerical model was developed and validated against a dataset obtained from controlled laboratory experiments.This dataset incorporated three-dimensional time-averaged velocity measurements.The skew-induced and stress-induced equation systems were adopted as the core governing equations,providing a framework for simulating various scenarios.A total of ten different cases were analyzed.The results highlighted the effect of discharge ratios on turbulence,lateral and vertical vorticities,and the distribution of mixing,which intensified with higher magnitudes of discharge ratios.The mixing structure,driven by velocity gradients and vorticity,revealed the significant role of lateral and vertical vorticities in determining hydrodynamic behaviors and mixing distributions at confluences.Specifically,the momentum ratio of incoming flows governed the spatial evolution of mixing processes.This study revealed that the distribution of mixing served as a key indicator for identifying the formation of mid-channel scours.High normalized velocities induced toward the left bank led to the superelevation of the water surface,enhancing the potential for bed material and the formation of significant scour holes beneath the elevated water surface.This novel approach provides a deeper understanding of the mixing patterns at confluences,particularly in scenarios with equilibrated discharge ratios but in different magnitudes.
基金Project supported by the National Natural Science Foundation of China(Nos.12372071 and 12372070)the Aeronautical Science Fund of China(No.2022Z055052001)the Foundation of China Scholarship Council(No.202306830079)。
文摘Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.
基金supported by the National Natural Science Foundation of China(Grant No.52375340,51975263,52405366).
文摘A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.
基金supported by the National Natural Science Foundation of China (Grant No.50479034)the Natural Science Foundation of Tianjin (Grant No.09YFSZSF02100)+1 种基金the financial support of the USDA Agriculture Research Service under Specific Research Agreement (Grant No. 58-6408-2-0062)(monitored by the USDA-ARS National Sedimentation Laboratory)the US State Department Agency for International Development under Agreement (Grant No.EE-G-00-02-00015-00) and the University of Mississippi
文摘Most natural rivers are curved channels, where the turbulent flows have a complex helical pattern, as has been extensively studied both numerically and experimentally. The helical flow structure in curved channels has an important bearing on sediment transport, riverbed evolution, and pollutant transport study. In this article, different turbulence closure schemes i.e., the mixing-length model and the k-ε model with different pressure solution techniques i. e., hydrostatic assumptions and dynamic pressure treatments are applied to study the helical secondary flows in an experiment curved channel. The agreements of vertically-averaged velocities between the simulated results obtained by using different turbulence models with different pressure solution techniques and the measured data are satisfactory. Their discrepancies with respect to surface elevations, superelevations and secondary flow patterns are discussed.
基金This research was funded jointly by the University of Kansas,Transportation Research Institute from Grant#DT0S59-06-G-00047,provided by the US Department of Transportation-Research and Innovative Technology Administration and PRS Mediterranean,Inc.in Israel.
文摘This paper summarizes the development of a three-dimensional numerical model for analyzing single geocell-reinforced soil.In this model,the infill soil was modeled using the Duncan-Chang model,which can simulate non-linearity and stress-dependency of soil.Geocell was modeled using linearly elastic plate elements,which can carry both bending and membrane stresses.A linear interface stress-strain relationship with a MohrCoulomb yield criterion was adopted to model the interface friction between the geocell wall and the soil.By modeling the geocell and the soil separately,the interaction between the soil and the geocell can be accurately simulated.To verify this model,a plate load test was conducted in the laboratory,in which a 12-cmthick sand layer reinforced by a single geocell was subjected to a vertical load from a circular steel plate.The load-displacement curves and the horizontal tensile strain of the geocell were recorded during the test.A numerical model was created according to the setup of the load test.The numerical results compared reasonably well with the test data.
基金funded by the National Natural Science Foundation of China(Nos.42077228,52174085)。
文摘With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rockbursts.To investigate the influence of tunnel diameter on the deformation and failure characteristics of surrounding rock,large-sized rocklike gypsum specimens were tested using a selfdeveloped true triaxial rockburst loading system containing circular tunnels with three different diameters(D=0.07 m,0.11 m,and 0.15 m).Acoustic emission monitoring,together with a miniature intelligent camera,was employed to analyze the entire process,focusing on macroscopic failure patterns,fragment characteristics,and underlying failure mechanisms.In addition,theoretical analyses were carried out and combined with numerical simulations to investigate the differences in energy evolution associated with rockburst physical models.The results indicate that:(1)The rockburst process with different tunnel diameters consistently evolved through three distinct stages—initial particle ejection,crack propagation accompanied by flake spalling,and,finally,fragment ejection leading to the formation of a‘V'-shaped notch.(2)Increasing tunnel diameter reduces rockburst failure load while increasing surrounding rock damage extent,total mass and average size of ejected fragments.Additionally,shear failure proportion decreases with tensile failure becoming increasingly dominant.(3)Larger tunnel diameters reduce the attenuation rate of elastic strain energy,thereby expanding the zone of elastic strain energy accumulation and disturbance and creating conditions for larger volume rockburst.(4)Larger tunnel diameters result in a smaller principal stress ratio at equivalent distances in the surrounding rock,indicating a higher likelihood of tensile failure.(5)Numerical analyses further reveal that larger tunnel diameters reduce the maximum elastic strain energy density around the tunnel,lowering the energy released per unit volume of rockburst fragments and their ejection velocities.However,both the total failure volume and overall energy release from rockburst increase.Model experiments with different tunnel diameters are of great significance for optimizing engineering design and parameter selection,as well as guiding tunnel construction under complex geological conditions.
基金supported by the Zhejiang Provincial Key Research and Development Program of China(No.2020C01020).
文摘The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbulence models struggle to make accurate predictions for subsonic and supersonic flows in nozzles.In this study,we explored a novel model,the algebraic stress model k-kL-ARSM+J,to enhance the accuracy of turbulence numerical simulations.This new model was used to conduct numerical simulations of the design and off-design performance of a 3D supersonic asymmetric truncated nozzle designed in our laboratory,with the aim of providing a realistic pattern of changes.The research indicates that,compared to linear eddy viscosity turbulence models such as k-kL and shear stress transport(SST),the k-kL-ARSM+J algebraic stress model shows better accuracy in predicting the performance of supersonic nozzles.Its predictions were identical to the experimental values,enabling precise calculations of the nozzle.The performance trends of the nozzle are as follows:as the inlet Mach number increases,both thrust and pitching moment increase,but the rate of increase slows down.Lift peaks near the design Mach number and then rapidly decreases.With increasing inlet pressure,the nozzle thrust,lift,and pitching moment all show linear growth.As the flight altitude rises,the internal flow field within the nozzle remains relatively consistent due to the same supersonic nozzle inlet flow conditions.However,external to the nozzle,the change in external flow pressure results in the nozzle exit transitioning from over-expanded to under-expanded,leading to a shear layer behind the nozzle that initially converges towards the nozzle center and then diverges.
基金supported by the National Natural Science Foundation of China(Grant Nos.U2242213,U2142213,42305167,42175105)。
文摘Systematic bias is a type of model error that can affect the accuracy of data assimilation and forecasting that must be addressed.An online bias correction scheme called the sequential bias correction scheme(SBCS),was developed using the6 h average bias to correct the systematic bias during model integration.The primary purpose of this study is to investigate the impact of the SBCS in the high-resolution China Meteorological Administration Meso-scale(CMA-MESO)numerical weather prediction(NWP)model to reduce the systematic bias and to improve the data assimilation and forecast results through this method.The SBCS is improved upon and applied to the CMA-MESO 3-km model in this study.Four-week sequential data assimilation and forecast experiments,driven by rapid update and cycling(RUC),were conducted for the period from 2–29 May 2022.In terms of the characteristics of systematic bias,both the background and analysis show diurnal bias,and these large biases are affected by complex underlying surfaces(e.g.,oceans,coasts,and mountains).After the application of the SBCS,the results of the data assimilation show that the SBCS can reduce the systematic bias of the background and yield a neutral to slightly positive result for the analysis fields.In addition,the SBCS can reduce forecast errors and improve forecast results,especially for surface variables.The above results indicate that this scheme has good prospects for high-resolution regional NWP models.
基金supported by National Natural Science Foundation of China(No.524B2078,12426307,51906203)Guangdong Major Project of Basic and Applied Basic Research(2023B0303000002)+6 种基金Guangdong Basic and Applied Basic Research Foundation(2023B1515120005)Natural Science Foundation of Shenzhen(JCYJ20241202125327036,JCYJ20240813100103005)Shenzhen Engineering Research Center of Redox Flow Battery for Energy Storage(XMHT20230208003)Research Project on Medium-and Long-Duration Flow Battery Energy Storage Technology(2024KJTW0015)China Association for Science and Technology(OR2308010)High level of special funds(G03034K001)supported by the Center for Computational Science and Engineering at the Southern University of Science and Technology.
文摘Long-duration energy storage has become critical for renewable energy integration.While redox flow batteries,especially vanadium-based systems,are scaling up in capacity,their performance at the stack level remains insufficiently optimized,demanding more profound mechanistic studies and engineering refinements.To address the difficulties in resolving the flow inhomogeneity at the stack scale,this study establishes a multi-physics field coupling model and analyzes the pressure distributions,flow rate differences,active substance concentration,and electrochemical characteristics.The results show that the uneven cell pressure distribution is a key factor affecting the consistency of the system performance,and the increase in the flow rate improves the reactant homogeneity,with both the average concentration and the uniformity factor increasing with the flow rate.In contrast,high current densities lead to an increased imbalance between electrochemical depletion and reactant replenishment,resulting in a significant decrease in reactant concentration in the under-ribs region.In addition,a higher flow rate can expand the high-current-density region where the stack operates efficiently.This study provides a theoretical basis for optimizing the design of the stack components.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金supported by the National Natural Science Foundation project“Micro-Scale Effect of Oil-Gas Flow and the Mechanism of Enhancing Shale Oil Recovery by Natural Gas Injection”(No.52074317)。
文摘Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques designed for tertiary oil recovery have garnered significant attention,with microgel flooding emerging as a particularly prominent area of research.Despite its promise,the complex mechanisms underlying microgel flooding have been rarely investigated numerically.This study aims to address these gaps by characterizing the distribution of microgel concentration and viscosity within different pore structures.To enhance the accuracy of these characterizations,the viscosity of microgels is adjusted to account for the shear effects induced by flow rate and the swelling effects driven by salinity variations.The absolute permeability of the rock and the relative permeability of both oil and microgel are also analyzed to elucidate the mechanisms of microgel flooding.Additionally,a connectivity model is employed to achieve a quantitative representation of fluid flow capacity.The proposed model is validated through conceptual examples and applied to real oilfield blocks,demonstrating its accuracy and practical applicability.
基金Funded by the National Natural Science Foundation of China Academy of Engineering Physics and Jointly Setup"NSAF"Joint Fund(No.U1430119)。
文摘The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.
基金supported by a grant from the Research Grant Council of Hong Kong Special Administrative Region(Project No.11207724).
文摘The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.
基金support from the National Key R&D plan(Grant No.2022YFC3004303)the National Natural Science Foundation of China(Grant No.42107161)+3 种基金the State Key Laboratory of Hydroscience and Hydraulic Engineering(Grant No.2021-KY-04)the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering(sklhse-2023-C-01)the Open Research Fund Program of Key Laboratory of the Hydrosphere of the Ministry of Water Resources(mklhs-2023-04)the China Three Gorges Corporation(XLD/2117).
文摘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.
