The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we pro...The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.展开更多
Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy...Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.展开更多
The mechanical behavior of nonwoven fabrics as reinforcement in cementitious composites remains insufficiently explored,particularly from a numerical modeling perspective,despite their growing interest as sustainable ...The mechanical behavior of nonwoven fabrics as reinforcement in cementitious composites remains insufficiently explored,particularly from a numerical modeling perspective,despite their growing interest as sustainable alternatives to conventional textiles.This study presents a simplified,engineering-oriented numerical modeling framework for reproducing the flexural mechanical response of cementitious composites reinforced with flax nonwoven fabric.Four-point bending(flexural)behavior of nonwoven fabric–reinforced cementitious composites was numerically simulated using ANSYS software.The model is developed using Finite Element Analysis(FEA)and incorporates a Representative Volume Element(RVE)approach to account for the heterogeneous fiber–matrix interaction.The required material properties were iteratively calibrated using existing experimental data for three composite configurations comprising 4,5,and 6 layers.The proposed model demonstrated agreement with experimental results within the investigated configurations,achieving normalized root mean square errors(nRMSE)of 5.2%,4.6%,and 2.07%for the respective configurations.Furthermore,correlations between material parameters and geometric factors were identified,providing preliminary insights for estimating model input properties from easily measurable variables.Finally,sensitivity analyses were performed to evaluate the influence of key geometric and material parameters on the structural response,offering valuable insights for the optimized design of nonwoven fabric-reinforced cementitious composites.展开更多
Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-al...Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-altitude regions.However,owing to the complex effects of low ambient pressure on strongly coupled atomized field sources,it is difficult to fully describe the comprehensive behaviour of such sources theoretically,which limits the further development of random splashing numerical models.In this paper,a refined random splashing numerical model characterized by low ambient pressure is developed based on experimental results and applied to high-altitude earth‒rockfill dam projects.Compared with the reference ambient pressure condition(P0=101.457 kPa),which corresponds to the same flood discharge flow,a decrease in ambient pressure by 0.1 P0leads to a maximum change rate not exceeding 10 m for the characteristic boundary of the 10 mm/h atomized rain intensity line at the QX Hydropower Station.This observation also applies to both the 40 mm/h and 10 mm/h atomized rain intensity lines at the RM Hydropower Station.For the two groups of flip bucket types designed for the RM Hydropower station,the loads associated with atomized protection are predominantly concentrated on the left bank.The maximum height of the 10 mm/h atomized rain intensity line affected by atomized rain ranges from 0.83 to 0.85 times the maximum dam height of 315 m.The distance between the farthest downstream boundary and the Spillway No.3 outlet is between 666.80 and 692.40 metres.Since the flip bucket shape variations only slightly affect the atomization zone extent,further optimization is needed.The study can provide valuable methodological and decision-making support for safeguarding against existing and potential impacts within areas affected by flood discharge atomization from high-altitude hydropower stations.展开更多
The predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simu...The predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simulated by a modified Leslie-Gower predator-prey model.Using the theory of reaction-diffusion equations,a priori estimate,existence,uniqueness and stability conditions of the positive steady state solution are established.Furthermore,numerical simulations are conducted to quantitatively analyze the dynamical behavior.The research shows that as long as the Allee effect constant satisfies the appropriate relationship and the growth rates of predator and prey are appropriately large,the predator and prey can not only coexist,but also the coexistence mode is unique and stable under low predation-rate.In addition,the numerical simulations show that the coexistence may be stable under high predation-rate.Meanwhile,with the increase of predation rate,the population density of predators will decrease.展开更多
Mathematical model-based accurate evaluation of the remediation process at organic pollution sites serves as an efficient approach to the management and remediation of contaminant source zones.Numerical and upscaled a...Mathematical model-based accurate evaluation of the remediation process at organic pollution sites serves as an efficient approach to the management and remediation of contaminant source zones.Numerical and upscaled analytical solution models are effective mathematical methods for reproducing the Dense Nonaqueous Phase Liquid(DNAPL)remediation process.However,in the current design of pollutant removal schemes,effective mass transfer models for characterizing the elution behaviors of contaminants remain lacking.In this study,two mathematical methods integrated with improved mass transfer models were employed to simulate the multi-stage contaminant elution behaviors under two distinct scenarios:A mixed-source region subjected to continuous water flushing and a residual DNAPL source treated with shorter-duration pulse flushing of the ethanol solution.Both the improved numerical model and upscaled analytical solution model demonstrated enhanced accuracy,which was attributed to the incorporation of solubilization mechanisms into mass transfer processes and the adoption of a multi-source region division method.The Mean Absolute Errors(MAE)of the numerical simulation for the two scenarios were 20.68 mg/L and 6.93 mg/L,respectively,whereas those of the upscaled model were 33.29 mg/L and 8.60 mg/L,respectively.Comparing the two improved models,the numerical model exhibited higher accuracy,while the upscaled model was characterized by faster computation speed and fewer input parameters.展开更多
The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combi...The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.展开更多
The Sichuan Basin,adjacent to the Himalayas in China,exhibits complex structural deformation and stress regimes,necessitating systematic structural and stress analysis.The Luzhou region in the southern Sichuan Basin i...The Sichuan Basin,adjacent to the Himalayas in China,exhibits complex structural deformation and stress regimes,necessitating systematic structural and stress analysis.The Luzhou region in the southern Sichuan Basin is cut by three sets of faults,which trend NNE(NE)-SSW(SW),E(NEE)-W(SWW),and N(NNW)-S(SSE).Late Ordovician to early Silurian strata lie between the basement and cover sequences,and the faults are mainly extrusion thrust faults.Integrated geophysical data and finite element numerical modelling reveal that the maximum horizontal principal stress(SH) is 76-131 MPa,the minimum horizontal principal stress(Sh) is 64-113 MPa,and the vertical stress(Sv) is 69-125 MPa,consistent with a strike-slip regime(SH>Sv> Sh).Consequently,the present-day stress regime differs from that of existing thrust-fault models.The NNE(NE)-SSW(SW) trending faults have the greatest impact on the magnitude and orientation of the stress,and the E(NEE)-W(SWW) faults have the greatest impact on the horizontal stress differences(ΔS).ΔS is negatively correlated with the total gas content and negatively correlated with borehole deformation in seven wells in the Luzhou region;therefore,E(NEE)-W(SWW) trending faults should be avoided during petroleum exploration in the Luzhou region.展开更多
Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and ...Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.展开更多
This work generalizes the subdiffusive Black-Scholes model by introducing the variable exponent in order to provide adequate descriptions for the option pricing,where the variable exponent may account for the variatio...This work generalizes the subdiffusive Black-Scholes model by introducing the variable exponent in order to provide adequate descriptions for the option pricing,where the variable exponent may account for the variation of the memory property.In addition to standard nonlinear-to-linear transformation,we apply a further spatial-temporal transformation to convert the model to a more tractable form in order to circumvent the difficulties caused by the"non-positive,non-monotonic"variable-exponent memory kernel.An interesting phenomenon is that the spatial transformation not only eliminates the advection term but naturally turns the original noncoercive spatial operator into a coercive one due to the specific structure of the Black-Scholes model,which thus avoids imposing constraints on coefficients.Then we perform numerical analysis for both the semi-discrete and fully discrete schemes to support numerical simulation.Numerical experiments are carried out to substantiate the theoretical results.展开更多
Metallic protective structures(e.g.,beams and plates)are widely used against impact and blast loadings.Their precise dynamic responses are critical for design and service,especially when unintended preloading or prest...Metallic protective structures(e.g.,beams and plates)are widely used against impact and blast loadings.Their precise dynamic responses are critical for design and service,especially when unintended preloading or prestress caused by accidental deformation is present.In this study,the effects of prestress on the structural deformation and springback behaviors of a fully clamped metallic beam subjected to a subsequent impact load are systematically investigated.A combined research approach,consisting of analytical modeling constructed by a simplified stringhinge model(SSHM)that accounts for the roles of structural hinge and string components as well as double-solver coupling numerical simulation incorporating implicit and explicit solvers simultaneously,is employed,which is validated against existing experimental results.The influence of material strain hardening is considered.The presence of prestress can improve the impact resistance of a beam by reducing its peak deflection and increasing the structural springback of the beam,owing mainly to altered beam geometries and,initially,the stress state as the beam is deformed.Using the analytical modeling of the SSHM,the roles of the components of the hinge and string under various loading scenarios are subsequently delineated,particularly in terms of the development process of the mechanical performance of each component within the out-of-plane deformation and springback stages.During the deformation process of the impacted target,the roles of the bending moment and membrane force vary with increasing midspan deflection.These roles also change when the pretension intensity is increased.展开更多
This study developed a two-dimensional storm surge model for hydrodynamic simulations in port engineering,utilizing an improved Local Time-Stepping(LTS)scheme.The model implements unstructured triangular grids with lo...This study developed a two-dimensional storm surge model for hydrodynamic simulations in port engineering,utilizing an improved Local Time-Stepping(LTS)scheme.The model implements unstructured triangular grids with localized refinement within the engineering area,enhancing computational efficiency through the improved LTS algorithm.Implementation in a Qingdao port demonstrated that,compared with the conventional Global Time-Stepping(GTS)scheme,the LTS approach enhanced computational efficiency by 5.08 times and 3.30 times before and after construction,respectively,reducing computation time by 30−40 hours.Validation results confirm the model’s high accuracy under both astronomical tide and storm surge conditions.Simulations of storm surges during Severe Typhoon Muifa(1109)and Super Typhoon Lekima(2019)further validated the model’s effectiveness.The results indicated that the overall simulation trend of storm surge aligned closely with the observed patterns of water level fluctuations,yielding highly satisfactory results.The Root Mean Square Error(RMSE)for Typhoon Lekima was 0.0083 m,while for Typhoon Muifa it was 0.0059 m,further demonstrating the model’s accuracy and applicability.The study analyzed the storm surge flow field post-construction using Typhoon Lekima as a case study.This research demonstrates the LTS model’s significant potential and promising applications in storm surge simulations and marine port construction.展开更多
Pipelines,as critical infrastructure for oiland gas transportation,require precise evaluation of peak loads in displacement-prone zones to ensure operational safety.The current design guidelines for lateral peak soil ...Pipelines,as critical infrastructure for oiland gas transportation,require precise evaluation of peak loads in displacement-prone zones to ensure operational safety.The current design guidelines for lateral peak soil resistance(ALA-2oO1 and PRCI-2oo9)are based on earlyanalytical studies with limited simulations and physical test data.These guidelines fail to adequately account for the coupled effects of soil friction and cohesion while also overlooking asymmetric soil constraints.These limitations raise significant concerns regarding their applicability in practical engineering scenarios,necessitatingthe development of more accurate analytical methods.The present study combines full-scale lateral pipesoil interaction tests with finite element modeling via the coupled Eulerian-Lagrangian approach in ABAQUS/Explicit.After validation,parametric studies were conducted toestablish a comprehensive database of lateral peak soil resistances.Based on the observed resistance development patterns,the lateral peak resistance calculation equation in ALA-2001 was modified to provide a more accurate analytical model capable of better reflecting real-world pipe-soil interaction behavior.The reliability of the proposed model was confirmed through independentphysical tests,demonstrating its significant value for pipeline engineering design and safety assessment.展开更多
Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at...Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at home and abroad into this formula,minimum self-supporting arch formulas under different surrounding rock classes were obtained.Based on the actual engineering case of a dual-mode shield tunnel,a numerical model for the tunnel boring machine excavation mode was established to verify the theoretical formulas.Next,three surrounding rock classes,four soil layer thickness gradients,and twelve overburden thickness gradients were designed,resulting in 144 models formed by the combination of the three factors.Uniform tests were conducted,and the pressure arch heights under different surrounding rock classes were obtained.The results show that in the theoretical formulas,the tunnel radius has a linear positive correlation with the pressure arch height,while the tunnel depth has a linear positive correlation with the square of the pressure arch height.According to numerical simulation results,the pressure arch height increases with the increase of the overburden thickness and then tends toward a critical value of twice the tunnel diameter.Finally,the results of the numerical model are in good agreement with those calculated using the theoretical formulas,verifying the rationality of the established theoretical formulas.展开更多
Marine-continental transitional shale(McTS)gas holds excellent gas-generating hydrocarbon basis and exploration potential.Conducting quantitative analysis on the evolution of shale gas content and the coupled relation...Marine-continental transitional shale(McTS)gas holds excellent gas-generating hydrocarbon basis and exploration potential.Conducting quantitative analysis on the evolution of shale gas content and the coupled relationship between hydrocarbon generation and storage during geological history is essential for a profound understanding of shale gas enrichment mechanisms.This studyestablishes integrated models for hydrocarbon generation evolution,porosity evolution,and gas occurrence in Type Ill organicrich MCTS through a synergistic experimental approach combining multi-temperature methane isothermal adsorption experiments andgold-tube pyrolysis experiments on low-maturity shale samples.Simulating a variety of real and virtual burial histories and thermal histories,the evolution process ofgas content in McTSwas reconstructed and the influence of various geological conditions during burial on gas content evolution was clarified.The results indicate that a seven-stage evolution(AG)of gas content in McTS from the Shanxi Formation,Southern North China Basin.Critical thresholds include:(1)dissolution-enhancedreservoir modification at vitrinite reflectance(EasyRo)=1.0%,(2)adsorbed gas saturation at EasyRo=1.3%,(3)dual saturation of free andadsorbedgas at EasyRo=2.0%,(4)15%30%gas loss through expulsion during overmature stages(EasyRo>2.0%),and(5)partial freeto-adsorbed gas conversion triggered by tectonic uplift.Total organic carbon(ToC)content and overpressure exhibit positive correlations with gas content,while tectonic uplift magnitude shows a negative impact.The influence of maximum burial depth,paleo-heat flow,andgeothermalgradient demonstrate complex nonlinear relationships on gas content.展开更多
Marine forecasting is critical for navigation safety and disaster prevention.However,traditional ocean numerical forecasting models are often limited by substantial errors and inadequate capture of temporal-spatial fe...Marine forecasting is critical for navigation safety and disaster prevention.However,traditional ocean numerical forecasting models are often limited by substantial errors and inadequate capture of temporal-spatial features.To address the limitations,the paper proposes a TimeXer-based numerical forecast correction model optimized by an exogenous-variable attention mechanism.The model treats target forecast values as internal variables,and incorporates historical temporal-spatial data and seven-day numerical forecast results from traditional models as external variables based on the embedding strategy of TimeXer.Using a self-attention structure,the model captures correlations between exogenous variables and target sequences,explores intrinsic multi-dimensional relationships,and subsequently corrects endogenous variables with the mined exogenous features.The model’s performance is evaluated using metrics including MSE(Mean Squared Error),MAE(Mean Absolute Error),RMSE(Root Mean Square Error),MAPE(Mean Absolute Percentage Error),MSPE(Mean Square Percentage Error),and computational time,with TimeXer and PatchTST models serving as benchmarks.Experiment results show that the proposed model achieves lower errors and higher correction accuracy for both one-day and seven-day forecasts.展开更多
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.展开更多
Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models...Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models for high-volatility(Na and Se),medium-volatility(Fe and Cu),and low-volatility(Ni and Cr)impurities,revealing dual threshold effects on impurity removal:low-to-medium temperatures(£550℃)effectively suppress volatilization,while elevated temperatures promote co-evaporation.At the optimal 550℃,tellurium purity reaches 5 N 8 with>90%yield.Spatial fractionation analysis demonstrates high-volatility impurities enriching in the upper condensation zone(X/L<0.25),whereas medium/low-volatility impurities accumulate in the lower zone and residues.Remarkably,residual impurities show significant enrichment versus raw materials-Na 4.04 times,Fe 13.3 times,Cu 53.6 times,Ni 7.17 times,and Cr 15.12 times,through formation of non-volatile compounds/solid solutions.Temperature-space coupling effects drive distinct deviation patterns:fluctuations in high-volatility impurities vs.temperature-progressive deviations in medium/low-volatility species.The impurity concentration of model-experiment deviations(mean±SD)are quantified as:Se 0.265±0.12,Na 0.224±0.02,Cu 0.146±0.06,Fe 0.133±0.13,Cr 0.101±0.07,and Ni 0.122±0.08.展开更多
This study proposed a numerical model based on a hybrid scheme using finite-difference and finite-volume methods to simulate nonlinear wave propagation from deep to shallow water.Six shock-capturing-type reconstructio...This study proposed a numerical model based on a hybrid scheme using finite-difference and finite-volume methods to simulate nonlinear wave propagation from deep to shallow water.Six shock-capturing-type reconstruction schemes,namely the second-order Monotonic Upwind Scheme for Conservation Laws(MUSCL),the second-order MUSCL-Total Variation Diminishing(TVD)scheme with the van Leer limiter,the standard fourth-order MUSCL-TVD scheme with the Minmod limiter,the improved fourth-order MUSCL-TVD scheme,the fifth-order Weighted Essentially Non-oscillatory(WENO)scheme,and the fifth-order TVD scheme with the Superbee slope limiter,were used to discretize the flux terms in the governing equations.Regular wave propagation in water of uniform depth,and over a submerged bar or a slope,was simulated using these schemes.Comparison of numerical results with theoretical solutions or experimental data illustrates the simulation performance of the six reconstruction schemes.展开更多
High-voltage electric pulse(HVEP)rock fragmentation has demonstrated substantial potential for sustainable fracturing of hard rocks owing to its energy efficiency.The transient nature and highly disruptive characteris...High-voltage electric pulse(HVEP)rock fragmentation has demonstrated substantial potential for sustainable fracturing of hard rocks owing to its energy efficiency.The transient nature and highly disruptive characteristics of its physical fracturing process render experimental investigation of the underlying rock-breaking mechanisms challenging.However,existing numerical studies lack comprehensive models that precisely link electrical breakdown phenomena with mechanical disintegration processes.This study combines COMSOL electrical breakdown simulations with four-dimension lattice spring model(4D-LSM)mechanical analysis to establish a coupled HVEP rock fragmentation model.The core concept of the model construction is to import the temperature field of the plasma channel obtained from the electrical breakdown into the mechanical solver to realize the precise connection between the two stages.The validated numerical model elucidates the full process of HVEP-induced fragmentation under varying electrical parameters.Furthermore,the effects of confining pressure and mineral grain size on fragmentation behavior have been investigated.Finally,parametric simulations across 25 electrical parameter combinations demonstrate the critical role of electrode spacing optimization in achieving energy-efficient rock fragmentation.These findings provide a predictive tool for designing efficient HVEP systems in deep resource extraction and mineral processing engineering.展开更多
基金supported by Excellent Research Group Project for Multiphase Evolution in Hyper-Gravity of the National Natural Science Foundation of China(No.52588202)。
摘要The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.
基金financially supported by the National Key Research and Development Program of China (No. 2023YFB3812601)the National Natural Science Foundation of China (No. 51925401)the Young Elite Scientists Sponsorship Program by CAST, China (No. 2022QNRC001)。
摘要Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.
基金funded by the research project PID2022-137156OB-I00financed by MCIN/AEI/10.13039/501100011033/FEDER,EU.
摘要The mechanical behavior of nonwoven fabrics as reinforcement in cementitious composites remains insufficiently explored,particularly from a numerical modeling perspective,despite their growing interest as sustainable alternatives to conventional textiles.This study presents a simplified,engineering-oriented numerical modeling framework for reproducing the flexural mechanical response of cementitious composites reinforced with flax nonwoven fabric.Four-point bending(flexural)behavior of nonwoven fabric–reinforced cementitious composites was numerically simulated using ANSYS software.The model is developed using Finite Element Analysis(FEA)and incorporates a Representative Volume Element(RVE)approach to account for the heterogeneous fiber–matrix interaction.The required material properties were iteratively calibrated using existing experimental data for three composite configurations comprising 4,5,and 6 layers.The proposed model demonstrated agreement with experimental results within the investigated configurations,achieving normalized root mean square errors(nRMSE)of 5.2%,4.6%,and 2.07%for the respective configurations.Furthermore,correlations between material parameters and geometric factors were identified,providing preliminary insights for estimating model input properties from easily measurable variables.Finally,sensitivity analyses were performed to evaluate the influence of key geometric and material parameters on the structural response,offering valuable insights for the optimized design of nonwoven fabric-reinforced cementitious composites.
基金funded by the State Key Laboratory of Hydraulic Engineering Intelligent Construction and Operation,Tianjin University(Grants No.HESS-2216)Nanxun scholars program of ZJWEU(Grant No.RC2024011034)+9 种基金Joint Funds of the Zhejiang Provincial Natural Science Foundation of China(Grant No.LZJWY22E090003)Project funded by China Postdoctoral Science Foundation(Grant No.2023M732603)Science and Technology Projects of Xizang Autonomous Region,China(Grant No.XZ202501ZY0109)National Natural Science Foundation of China(Grant No.U1765202)Natural Science Foundation of Tianjin(Grant No.22JCYBJC01180)Major Science and Technology Program of Zhejiang Province(Grant No.2021C03019)Zhejiang Key Laboratory of River-Lake Water Network Health Restoration(Grant No.HHSWKF202501)Central Guidance Funds for Science and Technology Local Development Project(Grant No.2025ZY01091)Scientific Research Fund Key Project of Zhejiang Institute of Hydraulics&Estuary(Zhejiang Institute of Marine Planning&Design)(Grant No.ZIHE25Z002)University-Level Key Course of Zhejiang University of Water Resources and Electric Power(Grant No.ZDKC202319).
摘要Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-altitude regions.However,owing to the complex effects of low ambient pressure on strongly coupled atomized field sources,it is difficult to fully describe the comprehensive behaviour of such sources theoretically,which limits the further development of random splashing numerical models.In this paper,a refined random splashing numerical model characterized by low ambient pressure is developed based on experimental results and applied to high-altitude earth‒rockfill dam projects.Compared with the reference ambient pressure condition(P0=101.457 kPa),which corresponds to the same flood discharge flow,a decrease in ambient pressure by 0.1 P0leads to a maximum change rate not exceeding 10 m for the characteristic boundary of the 10 mm/h atomized rain intensity line at the QX Hydropower Station.This observation also applies to both the 40 mm/h and 10 mm/h atomized rain intensity lines at the RM Hydropower Station.For the two groups of flip bucket types designed for the RM Hydropower station,the loads associated with atomized protection are predominantly concentrated on the left bank.The maximum height of the 10 mm/h atomized rain intensity line affected by atomized rain ranges from 0.83 to 0.85 times the maximum dam height of 315 m.The distance between the farthest downstream boundary and the Spillway No.3 outlet is between 666.80 and 692.40 metres.Since the flip bucket shape variations only slightly affect the atomization zone extent,further optimization is needed.The study can provide valuable methodological and decision-making support for safeguarding against existing and potential impacts within areas affected by flood discharge atomization from high-altitude hydropower stations.
基金Supported by the National Natural Science Foundation of China(11961030),the Natural Science Foundation of Shaanxi Province(2022JM-034)。
摘要The predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simulated by a modified Leslie-Gower predator-prey model.Using the theory of reaction-diffusion equations,a priori estimate,existence,uniqueness and stability conditions of the positive steady state solution are established.Furthermore,numerical simulations are conducted to quantitatively analyze the dynamical behavior.The research shows that as long as the Allee effect constant satisfies the appropriate relationship and the growth rates of predator and prey are appropriately large,the predator and prey can not only coexist,but also the coexistence mode is unique and stable under low predation-rate.In addition,the numerical simulations show that the coexistence may be stable under high predation-rate.Meanwhile,with the increase of predation rate,the population density of predators will decrease.
基金supported by the Qinghai Institute of Technology"Kunlun Talents"Talent Introduction Research Project(2023-QLGKLYCZX-002)the Qinghai Institute of Technology New Teacher Thesis Extension Research and Cultivation Project(2023011wys005)the National Natural Science Foundation of China(42177077).
摘要Mathematical model-based accurate evaluation of the remediation process at organic pollution sites serves as an efficient approach to the management and remediation of contaminant source zones.Numerical and upscaled analytical solution models are effective mathematical methods for reproducing the Dense Nonaqueous Phase Liquid(DNAPL)remediation process.However,in the current design of pollutant removal schemes,effective mass transfer models for characterizing the elution behaviors of contaminants remain lacking.In this study,two mathematical methods integrated with improved mass transfer models were employed to simulate the multi-stage contaminant elution behaviors under two distinct scenarios:A mixed-source region subjected to continuous water flushing and a residual DNAPL source treated with shorter-duration pulse flushing of the ethanol solution.Both the improved numerical model and upscaled analytical solution model demonstrated enhanced accuracy,which was attributed to the incorporation of solubilization mechanisms into mass transfer processes and the adoption of a multi-source region division method.The Mean Absolute Errors(MAE)of the numerical simulation for the two scenarios were 20.68 mg/L and 6.93 mg/L,respectively,whereas those of the upscaled model were 33.29 mg/L and 8.60 mg/L,respectively.Comparing the two improved models,the numerical model exhibited higher accuracy,while the upscaled model was characterized by faster computation speed and fewer input parameters.
基金Guangzhou Metro Scientific Research Project(No.JT204-100111-23001)Chongqing Municipal Special Project for Technological Innovation and Application Development(No.CSTB2022TIAD-KPX0101)Science and Technology Research and Development Program of China State Railway Group Co.,Ltd.(No.N2023G045)。
摘要The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.
基金supported by the National Natural Science Foundation of China (42372179)the Geological Survey Project of China (DD20240200605)the China National Petroleum Corporation Technology Project (2023ZZ02)。
摘要The Sichuan Basin,adjacent to the Himalayas in China,exhibits complex structural deformation and stress regimes,necessitating systematic structural and stress analysis.The Luzhou region in the southern Sichuan Basin is cut by three sets of faults,which trend NNE(NE)-SSW(SW),E(NEE)-W(SWW),and N(NNW)-S(SSE).Late Ordovician to early Silurian strata lie between the basement and cover sequences,and the faults are mainly extrusion thrust faults.Integrated geophysical data and finite element numerical modelling reveal that the maximum horizontal principal stress(SH) is 76-131 MPa,the minimum horizontal principal stress(Sh) is 64-113 MPa,and the vertical stress(Sv) is 69-125 MPa,consistent with a strike-slip regime(SH>Sv> Sh).Consequently,the present-day stress regime differs from that of existing thrust-fault models.The NNE(NE)-SSW(SW) trending faults have the greatest impact on the magnitude and orientation of the stress,and the E(NEE)-W(SWW) faults have the greatest impact on the horizontal stress differences(ΔS).ΔS is negatively correlated with the total gas content and negatively correlated with borehole deformation in seven wells in the Luzhou region;therefore,E(NEE)-W(SWW) trending faults should be avoided during petroleum exploration in the Luzhou region.
基金supported by the China Postdoctoral Science Foundation(CPSF)(Grant No.2024M762769)the Natural Science Basic Research Program of Shaanxi(Grant No.2024JC-YBQN-0333)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20232230).
摘要Slopes are likely to fail in areas with frequent rainfall and earthquakes.The deformation characteristics of unsaturated slopes subjected to post-rainfall earthquakes are investigated using centrifuge model tests and finite element analyses.Three tests of the slope deformation under earthquake and post-rainfall earthquakes are first studied using image analysis techniques.Then,based on an elastoplastic constitutive model,numerical simulations are carried out using the finite element method and compared with the centrifuge test results.Finally,a parametric study is performed to clarify the effects of antecedent rainfall on earthquake-induced slope deformation.The results show that slope deformation caused by post-rainfall earthquakes differs from that caused by earthquakes without antecedent rainfall.The seepage flow and soil strength of the slope are affected by previous rainfall conditions,such as intensity and duration,which directly influence the slope deformation caused by the subsequent earthquake.Soil displacement and strain become greater and the slip surface is more noticeable during the post-rainfall earthquake of higher intensity.In addition,the time interval between the rainfall and the earthquake has a considerable impact on the detailed characteristics of the slope deformation,and the significant deformation occurs at the time of lowest soil strength when seepage flow reaches the lower part of the slope.Moreover,the repeated intermittent rainfall greatly affects the subsequent earthquake-induced slope deformation,the main characteristics of which are closely related to the changes in saturation and strength of the slope.However,with the prolonged time gap between each round of rainfall,the earthquake-induced slope deformation becomes insignificant.
基金supported by the National Social Science Foundation of China(24BTJ006)the Taishan Scholars Program of Shandong Province(tsqn202306250).
摘要This work generalizes the subdiffusive Black-Scholes model by introducing the variable exponent in order to provide adequate descriptions for the option pricing,where the variable exponent may account for the variation of the memory property.In addition to standard nonlinear-to-linear transformation,we apply a further spatial-temporal transformation to convert the model to a more tractable form in order to circumvent the difficulties caused by the"non-positive,non-monotonic"variable-exponent memory kernel.An interesting phenomenon is that the spatial transformation not only eliminates the advection term but naturally turns the original noncoercive spatial operator into a coercive one due to the specific structure of the Black-Scholes model,which thus avoids imposing constraints on coefficients.Then we perform numerical analysis for both the semi-discrete and fully discrete schemes to support numerical simulation.Numerical experiments are carried out to substantiate the theoretical results.
基金supported by the National Natural Science Foundation of China(Grant No.12502178)the Youth Foundation of Rocket Force University of Engineering(Grant No.2024QN-B018)the Innovation Project of Fundamental Technology Frontier(Grant No.2025QYCX-MS-03-024)。
摘要Metallic protective structures(e.g.,beams and plates)are widely used against impact and blast loadings.Their precise dynamic responses are critical for design and service,especially when unintended preloading or prestress caused by accidental deformation is present.In this study,the effects of prestress on the structural deformation and springback behaviors of a fully clamped metallic beam subjected to a subsequent impact load are systematically investigated.A combined research approach,consisting of analytical modeling constructed by a simplified stringhinge model(SSHM)that accounts for the roles of structural hinge and string components as well as double-solver coupling numerical simulation incorporating implicit and explicit solvers simultaneously,is employed,which is validated against existing experimental results.The influence of material strain hardening is considered.The presence of prestress can improve the impact resistance of a beam by reducing its peak deflection and increasing the structural springback of the beam,owing mainly to altered beam geometries and,initially,the stress state as the beam is deformed.Using the analytical modeling of the SSHM,the roles of the components of the hinge and string under various loading scenarios are subsequently delineated,particularly in terms of the development process of the mechanical performance of each component within the out-of-plane deformation and springback stages.During the deformation process of the impacted target,the roles of the bending moment and membrane force vary with increasing midspan deflection.These roles also change when the pretension intensity is increased.
基金financially supported by the National Natural Science Foundation of China(Grant No.52071306)the Natural Science Foundation of Zhejiang Province(Grant No.LZ22E090003).
摘要This study developed a two-dimensional storm surge model for hydrodynamic simulations in port engineering,utilizing an improved Local Time-Stepping(LTS)scheme.The model implements unstructured triangular grids with localized refinement within the engineering area,enhancing computational efficiency through the improved LTS algorithm.Implementation in a Qingdao port demonstrated that,compared with the conventional Global Time-Stepping(GTS)scheme,the LTS approach enhanced computational efficiency by 5.08 times and 3.30 times before and after construction,respectively,reducing computation time by 30−40 hours.Validation results confirm the model’s high accuracy under both astronomical tide and storm surge conditions.Simulations of storm surges during Severe Typhoon Muifa(1109)and Super Typhoon Lekima(2019)further validated the model’s effectiveness.The results indicated that the overall simulation trend of storm surge aligned closely with the observed patterns of water level fluctuations,yielding highly satisfactory results.The Root Mean Square Error(RMSE)for Typhoon Lekima was 0.0083 m,while for Typhoon Muifa it was 0.0059 m,further demonstrating the model’s accuracy and applicability.The study analyzed the storm surge flow field post-construction using Typhoon Lekima as a case study.This research demonstrates the LTS model’s significant potential and promising applications in storm surge simulations and marine port construction.
基金co-financed by the Key Science and Technology Project of Ministry of Emergency Management of the People's Republic of China(Grant No.2024EMST090903)the National KeyR&D Program ofChina(Grant No.2022YFC3070100)the Young Elite Scientists Sponsorship Program by Beijing Association for Science and Technology,China(Grant No.BYESS2023261).
摘要Pipelines,as critical infrastructure for oiland gas transportation,require precise evaluation of peak loads in displacement-prone zones to ensure operational safety.The current design guidelines for lateral peak soil resistance(ALA-2oO1 and PRCI-2oo9)are based on earlyanalytical studies with limited simulations and physical test data.These guidelines fail to adequately account for the coupled effects of soil friction and cohesion while also overlooking asymmetric soil constraints.These limitations raise significant concerns regarding their applicability in practical engineering scenarios,necessitatingthe development of more accurate analytical methods.The present study combines full-scale lateral pipesoil interaction tests with finite element modeling via the coupled Eulerian-Lagrangian approach in ABAQUS/Explicit.After validation,parametric studies were conducted toestablish a comprehensive database of lateral peak soil resistances.Based on the observed resistance development patterns,the lateral peak resistance calculation equation in ALA-2001 was modified to provide a more accurate analytical model capable of better reflecting real-world pipe-soil interaction behavior.The reliability of the proposed model was confirmed through independentphysical tests,demonstrating its significant value for pipeline engineering design and safety assessment.
基金The National Natural Science Foundation of China(No.52478426)the Natural Science Foundation of Hunan Province(No.2024JJ5428).
摘要Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at home and abroad into this formula,minimum self-supporting arch formulas under different surrounding rock classes were obtained.Based on the actual engineering case of a dual-mode shield tunnel,a numerical model for the tunnel boring machine excavation mode was established to verify the theoretical formulas.Next,three surrounding rock classes,four soil layer thickness gradients,and twelve overburden thickness gradients were designed,resulting in 144 models formed by the combination of the three factors.Uniform tests were conducted,and the pressure arch heights under different surrounding rock classes were obtained.The results show that in the theoretical formulas,the tunnel radius has a linear positive correlation with the pressure arch height,while the tunnel depth has a linear positive correlation with the square of the pressure arch height.According to numerical simulation results,the pressure arch height increases with the increase of the overburden thickness and then tends toward a critical value of twice the tunnel diameter.Finally,the results of the numerical model are in good agreement with those calculated using the theoretical formulas,verifying the rationality of the established theoretical formulas.
基金supported by theNational Natural Science Foundation of China(No.42472210 and U25D9024)geological survey project of the China Geological Survey Oil and Gas Survey(No.[2024]02-07-03)the Grants-in-Aid of American Association of Petroleum Geologists(AAPG).
摘要Marine-continental transitional shale(McTS)gas holds excellent gas-generating hydrocarbon basis and exploration potential.Conducting quantitative analysis on the evolution of shale gas content and the coupled relationship between hydrocarbon generation and storage during geological history is essential for a profound understanding of shale gas enrichment mechanisms.This studyestablishes integrated models for hydrocarbon generation evolution,porosity evolution,and gas occurrence in Type Ill organicrich MCTS through a synergistic experimental approach combining multi-temperature methane isothermal adsorption experiments andgold-tube pyrolysis experiments on low-maturity shale samples.Simulating a variety of real and virtual burial histories and thermal histories,the evolution process ofgas content in McTSwas reconstructed and the influence of various geological conditions during burial on gas content evolution was clarified.The results indicate that a seven-stage evolution(AG)of gas content in McTS from the Shanxi Formation,Southern North China Basin.Critical thresholds include:(1)dissolution-enhancedreservoir modification at vitrinite reflectance(EasyRo)=1.0%,(2)adsorbed gas saturation at EasyRo=1.3%,(3)dual saturation of free andadsorbedgas at EasyRo=2.0%,(4)15%30%gas loss through expulsion during overmature stages(EasyRo>2.0%),and(5)partial freeto-adsorbed gas conversion triggered by tectonic uplift.Total organic carbon(ToC)content and overpressure exhibit positive correlations with gas content,while tectonic uplift magnitude shows a negative impact.The influence of maximum burial depth,paleo-heat flow,andgeothermalgradient demonstrate complex nonlinear relationships on gas content.
基金supported by the National Key Research and Development Program Project(2023YFC3107804)Planning Fund Project of Humanities and Social Sciences Research of the Ministry of Education(24YJA880097)the Graduate Education Reform Project in North China University of Technology(217051360025XN095-17)。
摘要Marine forecasting is critical for navigation safety and disaster prevention.However,traditional ocean numerical forecasting models are often limited by substantial errors and inadequate capture of temporal-spatial features.To address the limitations,the paper proposes a TimeXer-based numerical forecast correction model optimized by an exogenous-variable attention mechanism.The model treats target forecast values as internal variables,and incorporates historical temporal-spatial data and seven-day numerical forecast results from traditional models as external variables based on the embedding strategy of TimeXer.Using a self-attention structure,the model captures correlations between exogenous variables and target sequences,explores intrinsic multi-dimensional relationships,and subsequently corrects endogenous variables with the mined exogenous features.The model’s performance is evaluated using metrics including MSE(Mean Squared Error),MAE(Mean Absolute Error),RMSE(Root Mean Square Error),MAPE(Mean Absolute Percentage Error),MSPE(Mean Square Percentage Error),and computational time,with TimeXer and PatchTST models serving as benchmarks.Experiment results show that the proposed model achieves lower errors and higher correction accuracy for both one-day and seven-day forecasts.
基金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.
基金Project(2023YFC2907904)supported by the National Key Research and Development Program of ChinaProject(2024RC3032)supported by the Science and Technology Innovation Program of Hunan Province,China+2 种基金Project(2024JJ4057)supported by the Natural Science Foundation of Hunan Province,ChinaProject(YDZX2024091)supported by the Shandong Provincial Technology Innovation Guidance Plan,ChinaProject(YPML-20240502090)supported by the Key R&D Program of Yunnan Precious Metals Laboratory,China。
摘要Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models for high-volatility(Na and Se),medium-volatility(Fe and Cu),and low-volatility(Ni and Cr)impurities,revealing dual threshold effects on impurity removal:low-to-medium temperatures(£550℃)effectively suppress volatilization,while elevated temperatures promote co-evaporation.At the optimal 550℃,tellurium purity reaches 5 N 8 with>90%yield.Spatial fractionation analysis demonstrates high-volatility impurities enriching in the upper condensation zone(X/L<0.25),whereas medium/low-volatility impurities accumulate in the lower zone and residues.Remarkably,residual impurities show significant enrichment versus raw materials-Na 4.04 times,Fe 13.3 times,Cu 53.6 times,Ni 7.17 times,and Cr 15.12 times,through formation of non-volatile compounds/solid solutions.Temperature-space coupling effects drive distinct deviation patterns:fluctuations in high-volatility impurities vs.temperature-progressive deviations in medium/low-volatility species.The impurity concentration of model-experiment deviations(mean±SD)are quantified as:Se 0.265±0.12,Na 0.224±0.02,Cu 0.146±0.06,Fe 0.133±0.13,Cr 0.101±0.07,and Ni 0.122±0.08.
基金supported by the National Natural Science Foundation of China(Nos.51679132,11602136 and U22A20216)the Science and Technology Commission of Shanghai Municipality(No.21ZR1427000)the Shanghai Frontiers Science Center of‘Full Penetration’Far-Reaching Offshore Ocean Energy and Power.
摘要This study proposed a numerical model based on a hybrid scheme using finite-difference and finite-volume methods to simulate nonlinear wave propagation from deep to shallow water.Six shock-capturing-type reconstruction schemes,namely the second-order Monotonic Upwind Scheme for Conservation Laws(MUSCL),the second-order MUSCL-Total Variation Diminishing(TVD)scheme with the van Leer limiter,the standard fourth-order MUSCL-TVD scheme with the Minmod limiter,the improved fourth-order MUSCL-TVD scheme,the fifth-order Weighted Essentially Non-oscillatory(WENO)scheme,and the fifth-order TVD scheme with the Superbee slope limiter,were used to discretize the flux terms in the governing equations.Regular wave propagation in water of uniform depth,and over a submerged bar or a slope,was simulated using these schemes.Comparison of numerical results with theoretical solutions or experimental data illustrates the simulation performance of the six reconstruction schemes.
基金financial support from the National Natural Science Foundation of China(Nos.52209144 and 12472405).
摘要High-voltage electric pulse(HVEP)rock fragmentation has demonstrated substantial potential for sustainable fracturing of hard rocks owing to its energy efficiency.The transient nature and highly disruptive characteristics of its physical fracturing process render experimental investigation of the underlying rock-breaking mechanisms challenging.However,existing numerical studies lack comprehensive models that precisely link electrical breakdown phenomena with mechanical disintegration processes.This study combines COMSOL electrical breakdown simulations with four-dimension lattice spring model(4D-LSM)mechanical analysis to establish a coupled HVEP rock fragmentation model.The core concept of the model construction is to import the temperature field of the plasma channel obtained from the electrical breakdown into the mechanical solver to realize the precise connection between the two stages.The validated numerical model elucidates the full process of HVEP-induced fragmentation under varying electrical parameters.Furthermore,the effects of confining pressure and mineral grain size on fragmentation behavior have been investigated.Finally,parametric simulations across 25 electrical parameter combinations demonstrate the critical role of electrode spacing optimization in achieving energy-efficient rock fragmentation.These findings provide a predictive tool for designing efficient HVEP systems in deep resource extraction and mineral processing engineering.