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.展开更多
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.展开更多
Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbank...Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbanks were improved. Then a two-dimensional numerical model of the channel deformation was developed, consisting of a 2D flow and sediment transport submodel and bank-erosion submodels of different soil riverbanks. In the model, a new technique for updating the bank geometry during the bed evolution was presented, which combines closely two kinds of submodels. The proposed model is capable of not only predicting the processes of flood routing and longitudinal channel deformation in natural rivers, but also simulating the processes of lateral channel deformation, especially the processes of lateral erosion and failure of cohesive, non-cohesive and composite riverbanks.展开更多
Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subs...Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and effic...Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.展开更多
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 intensiti...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.展开更多
The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research ...The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.展开更多
In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two th...In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.展开更多
This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare ...This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.展开更多
To optimize the excavation of rock using underground blasting techniques,a reliable and simplified approach for modeling rock fragmentation is desired.This paper presents a multistep experimentalnumerical methodology ...To optimize the excavation of rock using underground blasting techniques,a reliable and simplified approach for modeling rock fragmentation is desired.This paper presents a multistep experimentalnumerical methodology for simplifying the three-dimensional(3D)to two-dimensional(2D)quasiplane-strain problem and reducing computational costs by more than 100-fold.First,in situ tests were conducted involving single-hole and free-face blasting of a dolomite rock mass in a 1050-m-deep mine.The results were validated by laser scanning.The craters were then compared with four analytical models to calculate the radius of the crushing zone.Next,a full 3D model for single-hole blasting was prepared and validated by simulating the crack length and the radius of the crushing zone.Based on the stable crack propagation zones observed in the 3D model and experiments,a 2D model was prepared.The properties of the high explosive(HE)were slightly reduced to match the shape and number of radial cracks and crushing zone radius between the 3D and 2D models.The final methodology was used to reproduce various cut-hole blasting scenarios and observe the effects of residual cracks in the rock mass on further fragmentation.The presence of preexisting cracks was found to be crucial for fragmentation,particularly when the borehole was situated near a free rock face.Finally,an optimization study was performed to determine the possibility of losing rock continuity at different positions within the well in relation to the free rock face.展开更多
Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-...Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-forming characteristics require further study.Therefore,taking the No.8 deep coal seam in the central-eastern region of Ordos Basin as an example,this study investigated the geologic characteristics of CBM accumulations to establish a numerical model.The evolution of the burial and accumulation of CBM in the area was reconstructed.The modeling results suggest that the No.8 coal seam experienced continuous subsidence from the Late Cretaceous to the Triassic,alternating subsidence and uplift during the Jurassic,rapid burial throughout the Early Cretaceous,and continuous uplift since the Late Cretaceous.The coal reached its maximum maturity at the end of the Early Cretaceous.Furthermore,CBM generation in the region was divided into four stages of thermal eventsdbiogenic and early thermogenic gas,cracking of light oil into gas,cracking of the remaining kerogen into gas,and hydrocarbon generation ceasingdwhich accelerated coal maturity and generation.The adsorption capacity presented an overall declining trend prior to the end of the Cretaceous,followed by a rapid increase since the Late Cretaceous.As for adsorption mass evolution,the CBM successively un-derwent unsaturated minor adsorption,unsaturated rapid-rising adsorption,saturated decreasing adsorption,and saturated rising adsorption.The in-situ gas mass was found to be controlled by a combination of generation,adsorption,and expulsion of hydrocarbons,with its present-day value being 9-29×10^(4)t/km^(2)and the corresponding gas volume per ton of coal being 12-28 m^(3)/t.Moreover,free gas evolution initially showed an increasing trend,followed by a decline,ultimately accounting for 11%-28%of the total gas content.展开更多
Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great ...Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.展开更多
基金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.
基金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.
基金the National Key Research Program of Basic Science Foundation (Grant No. G19990436) partially by the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant No. 50221903).
文摘Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbanks were improved. Then a two-dimensional numerical model of the channel deformation was developed, consisting of a 2D flow and sediment transport submodel and bank-erosion submodels of different soil riverbanks. In the model, a new technique for updating the bank geometry during the bed evolution was presented, which combines closely two kinds of submodels. The proposed model is capable of not only predicting the processes of flood routing and longitudinal channel deformation in natural rivers, but also simulating the processes of lateral channel deformation, especially the processes of lateral erosion and failure of cohesive, non-cohesive and composite riverbanks.
基金supported by National Natural Science Foundation of China(72288101,72361137002,and 72101018)the Dutch Research Council(NWO Grant 482.22.01).
文摘Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.
基金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 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.
基金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.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China (grant numbers: 41874160 and 92055213)。
文摘Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.
基金funding received by a grant from the Natural Sciences and Engineering Research Council of Canada(NSERC)(Grant No.CRDPJ 469057e14).
文摘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.
基金supported by the Research Grants Council of Hong Kong Special Administrative Region Government of China(Grant Nos.15210322 and R5037-18)the financial support(Grant No.86902-00000240)from Shenzhen University.
文摘The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.
基金supported by project XJZ2023050044,A2309002 and XJZ2023070052.
文摘In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.
基金This research has been conducted under SEGTRANS project,funded by the Centre for Industrial Technological Development(CDTI,Government of Spain).
文摘This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.
文摘To optimize the excavation of rock using underground blasting techniques,a reliable and simplified approach for modeling rock fragmentation is desired.This paper presents a multistep experimentalnumerical methodology for simplifying the three-dimensional(3D)to two-dimensional(2D)quasiplane-strain problem and reducing computational costs by more than 100-fold.First,in situ tests were conducted involving single-hole and free-face blasting of a dolomite rock mass in a 1050-m-deep mine.The results were validated by laser scanning.The craters were then compared with four analytical models to calculate the radius of the crushing zone.Next,a full 3D model for single-hole blasting was prepared and validated by simulating the crack length and the radius of the crushing zone.Based on the stable crack propagation zones observed in the 3D model and experiments,a 2D model was prepared.The properties of the high explosive(HE)were slightly reduced to match the shape and number of radial cracks and crushing zone radius between the 3D and 2D models.The final methodology was used to reproduce various cut-hole blasting scenarios and observe the effects of residual cracks in the rock mass on further fragmentation.The presence of preexisting cracks was found to be crucial for fragmentation,particularly when the borehole was situated near a free rock face.Finally,an optimization study was performed to determine the possibility of losing rock continuity at different positions within the well in relation to the free rock face.
基金funded by the applied research and technology project titled“Research on Key Technologies for Optimization Design of Deep Coalbed Methane Development”conducted by PetroChina Company Limited(Grant No.2023ZZ1804).
文摘Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-forming characteristics require further study.Therefore,taking the No.8 deep coal seam in the central-eastern region of Ordos Basin as an example,this study investigated the geologic characteristics of CBM accumulations to establish a numerical model.The evolution of the burial and accumulation of CBM in the area was reconstructed.The modeling results suggest that the No.8 coal seam experienced continuous subsidence from the Late Cretaceous to the Triassic,alternating subsidence and uplift during the Jurassic,rapid burial throughout the Early Cretaceous,and continuous uplift since the Late Cretaceous.The coal reached its maximum maturity at the end of the Early Cretaceous.Furthermore,CBM generation in the region was divided into four stages of thermal eventsdbiogenic and early thermogenic gas,cracking of light oil into gas,cracking of the remaining kerogen into gas,and hydrocarbon generation ceasingdwhich accelerated coal maturity and generation.The adsorption capacity presented an overall declining trend prior to the end of the Cretaceous,followed by a rapid increase since the Late Cretaceous.As for adsorption mass evolution,the CBM successively un-derwent unsaturated minor adsorption,unsaturated rapid-rising adsorption,saturated decreasing adsorption,and saturated rising adsorption.The in-situ gas mass was found to be controlled by a combination of generation,adsorption,and expulsion of hydrocarbons,with its present-day value being 9-29×10^(4)t/km^(2)and the corresponding gas volume per ton of coal being 12-28 m^(3)/t.Moreover,free gas evolution initially showed an increasing trend,followed by a decline,ultimately accounting for 11%-28%of the total gas content.
基金This research was funded by the National Natural Science Foundation of China(52079106,52009104,51609199)the National Key Research and Development Program of China(2016YFC0402704).
文摘Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.
