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.展开更多
Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at vary...In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.展开更多
Burden is one of the main parameters in blast design.However,field tests,either single-or multi-hole blasts,used to determine an appropriate burden,are difficult to capture crack propagation,evolution of breakage angl...Burden is one of the main parameters in blast design.However,field tests,either single-or multi-hole blasts,used to determine an appropriate burden,are difficult to capture crack propagation,evolution of breakage angle,and the mechanism governing these processes in the rock.In this study,a single-hole bench blasting model is developed using LS-DYNA software to comprehensively investigate the relationship between burden and rock breakage.The simulation results show that the breakage angle decreases with the increase in burden,and the blasted volume reaches a peak value with a burden of 4 m.Meanwhile,backbreak distance increases with increasing burden.The optimum burden in this simulation is found to be 4.0 m,as the ratio of burden to blasthole diameter is equal to 20.62 and the ratio of burden to bench height is 0.44,based on a comprehensive analysis of the blasted volume,average damage,and total damage.Under the optimum burden condition,tensile stress wave regions are simultaneously generated at the free surfaces of both the bench top and bench slope,allowing more effective utilization of the two free surfaces and resulting in a more uniform damage distribution within the burden region.展开更多
Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On ...Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.展开更多
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°.展开更多
Adhesively bonded joints are widely used in modern lightweight structures due to their high strengthto-weight ratio and design flexibility.However,the reliable non-destructive evaluation of bond integrity remains a si...Adhesively bonded joints are widely used in modern lightweight structures due to their high strengthto-weight ratio and design flexibility.However,the reliable non-destructive evaluation of bond integrity remains a significant challenge.This study presents a numerical investigation of adhesively bonded joints with different adhesive properties using ultrasonic guided waves.The main focus of the investigation is to evaluate the feasibility of using guided waves to assess bond integrity,particularly for detecting challenging weak bonds.For this purpose,a theoretical analysis of dispersion curves was conducted,revealing that the S0 Lamb wave mode is significantly sensitive to variations in adhesive properties in the 300-700 kHz frequency range.Finite element modelling was used to analyse the propagation of guided waves in two scenarios:an adhesively bonded aluminum structure and a more complex configuration-adhesively bonded lap joints.The Short-Time Fourier Transform(STFT)was used to process the obtained results and determine the group velocities of guided waves.By analysing the group velocity characteristics,their dependence on the adhesive properties was identified.In the first scenario,a clear separation of S0 modes from A0 modes was observed in the STFT analysis,with a decrease in group velocity as adhesive stiffness increased.For the more complex lap joint scenario,the separation between A0 and S0 modes was less distinct.However,the analysis of the average group velocity shows a dependence of average group velocity on adhesive properties.This is similar to the first scenario.There is a decrease in average group velocity as adhesive stiffness increases.The results obtained demonstrate that guided wavebased methods have a high potential for non-destructive evaluation of adhesively bonded structures,including the detection of weak bonds.展开更多
Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numer...Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numerically investigate the reaction process of hydrocarbon-containing VOCs in RCO using computational fluid dynamics(CFD)simulation.To obtain the conversion characteristics of multi-component hydrocarbons,the effects of intake load,equivalence ratio,and the composition of multi-component hydrocarbons on the flow,heat transfer,and conversion rate of the reactor were analyzed.A feasibility study plan targeting the hard-to-convert components was also proposed.The results indicated that as the load increases,the conversion rates of the various components decrease,while the reaction rates increase.Moreover,increasing the flow velocity intensifies turbulence and enhances the collision frequency between the gas and the wall surfaces.This,in turn,amplifies the resistance effect of the porous medium.As the equivalence ratio of VOCs to oxygen increases,the oxygen-deficient condition leads to a decrease in the molecular weight of the hydrocarbons involved in the reaction.The reaction temperature also shows a downward trend.A comparative analysis of the catalytic combustion characteristics of multi-component VOCs and single-component gases reveals that adding ethane and propane can facilitate methane oxidation.展开更多
Prenatal exposure to bisphenols and metals has raised significant concerns regarding their potential impact on fetal development,particularly the risk of fetal chromosome numerical abnormalities(CNA).In this case-cont...Prenatal exposure to bisphenols and metals has raised significant concerns regarding their potential impact on fetal development,particularly the risk of fetal chromosome numerical abnormalities(CNA).In this case-control study,we analyzed bisphenol and metal concentrations in amniotic fluid of high-risk pregnant women undergoing amniocentesis.Concentrations of bisphenols and metals were measured using ultra-performance liquid chromatography-tandem mass spectrometry and inductively coupled plasma-mass spectrometry,respectively.Logistic regression and quantile-based g-computation were applied to evaluate individual and combined effects,while dose-response relationships were assessed using restricted cubic splines.Our findings indicated that bisphenol S(BPS),bisphenol Z(BPZ),bisphenol AF(BPAF),antimony(Sb),and vanadium(V)were significantly associated with an increased risk of CNA when analyzed individually,whereas manganese,iron,copper(Cu),nickel(Ni),and zinc(Zn)were significantly and inversely associated with CNA risk.Combined exposure to bisphenol and metal mixtures was associated with an increased risk of CNA in multi-pollutant models.Cu and Ni exhibited a positive additive interaction.Furthermore,BPS,BPZ,and BPAF were individually associated with an increased risk of Down syndrome,while Zn was associated with a decreased risk of Down syndrome.BPS,Sb,V,and Zn were individually associated with an increased risk of Klinefelter syndrome.These findings underscore the potential role of prenatal bisphenol and metal exposure in the pathogenesis of fetal CNA,highlighting both additive and synergistic effects.展开更多
This study presents a systematic numerical analysis of wind loads on offshore photovoltaic(PV)panels.A computational fluid dynamics(CFD)model,incorporating a free-surface wave boundary condition,is developed and valid...This study presents a systematic numerical analysis of wind loads on offshore photovoltaic(PV)panels.A computational fluid dynamics(CFD)model,incorporating a free-surface wave boundary condition,is developed and validated against experimental data.Parametric investigations quantify the effects of wind speed,panel tilt angle,clearance,and wave characteristics on the aerodynamic coefficients(drag,lift,and moment).Results indicate that all force coefficients increase with wind speed,with the lift coefficient being most sensitive to wave action.While a larger tilt angle intensifies airflow disturbance and amplifies the coefficients,this effect is more pronounced over flat ground than above a wavy surface.As clearance increases,the drag coefficient fluctuates before rising,the lift coefficient exhibits a trough-shaped response,and the moment coefficient increases monotonically,with values consistently higher over waves.Furthermore,the aerodynamic coefficients generally decrease with greater wave height.The maximum wind load occurs directly above the wave trough,and the aerodynamic force coefficient varies non-monotonically with wave position,first decreasing and then increasing.These findings offer practical guidance for the structural design and safety assurance of offshore PV systems.展开更多
Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, w...Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.展开更多
Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagati...Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagation behavior.To address the unclear mechanisms governing fracture penetration across coal-gangue interfaces,this study employs the Continuum-Discontinuum Element Method(CDEM)to simulate and analyze the vertical propagation of hydraulic fractures initiating within coal seams,based on geomechanical parameters derived from the deep Benxi Formation coal seams in the southeastern Ordos Basin.The investigation systematically examines the influence of geological and operational parameters on cross-interfacial fracture growth.Results demonstrate that vertical stress difference,elastic modulus contrast between coal and gangue layers,interfacial stress differential,and interfacial cohesion at coal-gangue interfaces are critical factors governing hydraulic fracture penetration through these interfaces.High vertical stress differences(>3 MPa)inhibit interfacial dilation,promoting predominant crosslayer fracture propagation.Reduced interfacial stress contrasts and enhanced interfacial cohesion facilitate fracture penetration across interfaces.Furthermore,smaller elastic modulus contrasts between coal and gangue correlate with increased interfacial aperture.Finally,lower injection rates effectively suppress vertical fracture propagation in deep coal reservoirs.This study elucidates the characteristics and mechanisms governing cross-layer fracture propagation in coal–rock composites with interbedded partings,and delineates the dynamic evolution laws and dominant controlling factors involved.Thefindings provide critical theoretical insights for the optimization of fracture design and the efficient development of deep coalbed methane reservoirs.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
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.展开更多
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.展开更多
As binary geological media,soil-rock mixtures(SRMs)exhibit a distinct gradational composition,leading to their unique mechanical behaviors.To appraise the stability of SRM slopes,it is essential to determine equivalen...As binary geological media,soil-rock mixtures(SRMs)exhibit a distinct gradational composition,leading to their unique mechanical behaviors.To appraise the stability of SRM slopes,it is essential to determine equivalent parameters of SRMs,which are typically obtained through experimental and numerical methods.In contrasted to other numerical methods,the numerical manifold method(NMM)is more effective in addressing SRM problems.This is because the high-precision regular mathematical meshes in NMM can be used without aligning with the soil-rock interfaces and boundaries of SRMs.In the current research,the equivalent strength parameters of SRMs,i.e.the equivalent cohesion ce and internal friction angleϕ_(e),are determined using NMM.Initially,an NMM triaxial numerical model is established and validated based on triaxial experiments.Subsequently,the soil and rock parameters are derived through parameter inversion.Moreover,the impacts of rock content,size,shape and rock blocks'major-axis orientation on ce andϕ_(e) of SRMs are thoroughly examined using the NMM triaxial numerical model.Additionally,a fitting function is proposed to linkϕ_(e) to the rock content and size of SRMs.When other influencing factors are fixed,the above fitting model leads to the following conclusions:(1)the predictedϕ_(e) of SRMs increase with the increase of rock content;and(2)SRM samples with smaller rocks display a higher predictedϕ_(e).展开更多
Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers of...Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.展开更多
Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and stora...Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO_(2) trapping mechanisms.In the base case,a constant injection rate of 3500 m^(3)/d over fifteen years resulted in a cumulative injection of 19.2×10^(6) m^(3).The CO_(2) plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO_(2) in the reservoir was 10.6 million tonnes(Mt),representing less than 10% of its theoretical capacity.The plume was projected to reach the entrapment crest and the top of the reservoir within a century,indicating secure long-term containment.Structural,stratigraphic,and residual trapping dominate in Reservoir A(approximately 90%).Anticlinal closures with thick overlying mudstones in the Zhujiang Formation provide effective seals,further enhancing storage security.Reservoir properties and heterogeneity play a crucial role in controlling CO_(2) storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO_(2) capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.展开更多
Functionally graded material(FGM)plates are widely used in various engineering structures owing to their tailor-made mechanical properties,whereas cracked homogeneous plates constitute a canonical setting in fracture ...Functionally graded material(FGM)plates are widely used in various engineering structures owing to their tailor-made mechanical properties,whereas cracked homogeneous plates constitute a canonical setting in fracture mechanics analysis.These two classes of problems respectively embody material non-uniformity and geometric discontinuity,thereby imposing more stringent requirements on numerical methods in terms of high-order field continuity and accurate defect representation.Based on the classical Kirchhoff-Love plate theory,a numerical manifold method(MLS-NMM)incorporating moving least squares(MLS)interpolation is developed for bending analysis of FGM plates and fracture simulation of homogeneous plates with defects.The method constructs an H^(2)-regular approximation with high-order continuous weighting functions and,combined with the separation of mathematical and physical covers,establishes a unified framework that accurately handles material gradients and cracks without mesh reconstruction.For the crack tip,a singular physical cover incorporating the Williams asymptotic field is introduced to achieve local enrichment,enabling the natural capture of displacement discontinuity and stress singularity.Stress intensity factors are extracted using the interaction integral method,and the dimensionless J-integral shows a maximum relative error below 1.2%compared with the reference solution.Numerical results indicate that MLS-NMM exhibits excellent convergence performance:using 676 mathematical nodes,the nondimensional central deflection of both FGM and homogeneous plates agrees with reference solutions with a maximum relative error below 0.81%,and no shear locking occurs.A systematic analysis reveals that for a simply supported on all four edges(SSSS)FGM square plate with a/h=10,the nondimensional central deflection increases by 212%as the gradient index nrises from 0 to 5.For a homogeneous plate containing a central crack with c/a=0.6,the nondimensional central deflection increases by approximately 46%compared with the intact plate.Under weak boundary constraints(e.g.,SFSF),the deformation is markedly amplified,with the deflection reaching more than three times that under strong constraints(SCSC).The proposed method provides an efficient,reconstruction-free numerical tool for high-accuracy bending and fracture analyses of FGM and cracked thin-plate structures.展开更多
Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteri...Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteristics and flow field evolution laws of overlapping rotor configurations in hovering conditions through numerical simulation methods.The research method involves constructing a computational model for rotor flow fields and aerodynamic characteristics based on the Reynolds-averaged Navier-Stokes(RANS)equations and the Spalart-Allmaras(S-A)turbulence model.The dynamic simulation of rotor rotational motion was achieved by using the moving nested grid technology.The reliability of the computational method was ensured through the grid independence verification and the comparison with experimental data.The research results indicate that in overlapping rotor systems,rotorⅡexperiences a decrease in thrust,significant power fluctuations,and reduced hovering efficiency due to continuous interference from the adjacent rotor’s wake and blade-vortex interactions.Blade-tip vortices undergo breakage,fusion,and secondary rolling in the overlapping region,forming large-scale turbulent structures that lead to attenuation of the induced velocity field and aerodynamic efficiency losses.Additionally,the interaction between the rotor downwash and the fuselage triggers a“fountain effect”and a sudden increase in surface pressure on the fuselage,exacerbating flow field distortion.Based on the aforementioned mechanisms,the safe flight of overlapping rotor configurations can be achieved by optimizing the configuration strategy of the rotational speed phase difference between adjacent blades.This study provides a theoretical basis for the rotor layout design and the aerodynamic performance enhancement of heavy-load eVTOL aircraft.展开更多
基金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.
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
基金The National Natural Science Foundation of China(Grant No.52201376)the Natural Science Foundation of Hubei Province,China(Grant No.2023AFB683).
文摘In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.
基金supported by the European Union in the frame of Horizon Europe AVANTIS project,Grant Agreement No.101137552.
文摘Burden is one of the main parameters in blast design.However,field tests,either single-or multi-hole blasts,used to determine an appropriate burden,are difficult to capture crack propagation,evolution of breakage angle,and the mechanism governing these processes in the rock.In this study,a single-hole bench blasting model is developed using LS-DYNA software to comprehensively investigate the relationship between burden and rock breakage.The simulation results show that the breakage angle decreases with the increase in burden,and the blasted volume reaches a peak value with a burden of 4 m.Meanwhile,backbreak distance increases with increasing burden.The optimum burden in this simulation is found to be 4.0 m,as the ratio of burden to blasthole diameter is equal to 20.62 and the ratio of burden to bench height is 0.44,based on a comprehensive analysis of the blasted volume,average damage,and total damage.Under the optimum burden condition,tensile stress wave regions are simultaneously generated at the free surfaces of both the bench top and bench slope,allowing more effective utilization of the two free surfaces and resulting in a more uniform damage distribution within the burden region.
文摘Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.
基金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°.
基金supported by the Research Council of Lithuania(LMTLT),agreement no.S-MIP-22-5.
文摘Adhesively bonded joints are widely used in modern lightweight structures due to their high strengthto-weight ratio and design flexibility.However,the reliable non-destructive evaluation of bond integrity remains a significant challenge.This study presents a numerical investigation of adhesively bonded joints with different adhesive properties using ultrasonic guided waves.The main focus of the investigation is to evaluate the feasibility of using guided waves to assess bond integrity,particularly for detecting challenging weak bonds.For this purpose,a theoretical analysis of dispersion curves was conducted,revealing that the S0 Lamb wave mode is significantly sensitive to variations in adhesive properties in the 300-700 kHz frequency range.Finite element modelling was used to analyse the propagation of guided waves in two scenarios:an adhesively bonded aluminum structure and a more complex configuration-adhesively bonded lap joints.The Short-Time Fourier Transform(STFT)was used to process the obtained results and determine the group velocities of guided waves.By analysing the group velocity characteristics,their dependence on the adhesive properties was identified.In the first scenario,a clear separation of S0 modes from A0 modes was observed in the STFT analysis,with a decrease in group velocity as adhesive stiffness increased.For the more complex lap joint scenario,the separation between A0 and S0 modes was less distinct.However,the analysis of the average group velocity shows a dependence of average group velocity on adhesive properties.This is similar to the first scenario.There is a decrease in average group velocity as adhesive stiffness increases.The results obtained demonstrate that guided wavebased methods have a high potential for non-destructive evaluation of adhesively bonded structures,including the detection of weak bonds.
基金supported by National Key Research&Development Program of China(2022YFB4101500).
文摘Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numerically investigate the reaction process of hydrocarbon-containing VOCs in RCO using computational fluid dynamics(CFD)simulation.To obtain the conversion characteristics of multi-component hydrocarbons,the effects of intake load,equivalence ratio,and the composition of multi-component hydrocarbons on the flow,heat transfer,and conversion rate of the reactor were analyzed.A feasibility study plan targeting the hard-to-convert components was also proposed.The results indicated that as the load increases,the conversion rates of the various components decrease,while the reaction rates increase.Moreover,increasing the flow velocity intensifies turbulence and enhances the collision frequency between the gas and the wall surfaces.This,in turn,amplifies the resistance effect of the porous medium.As the equivalence ratio of VOCs to oxygen increases,the oxygen-deficient condition leads to a decrease in the molecular weight of the hydrocarbons involved in the reaction.The reaction temperature also shows a downward trend.A comparative analysis of the catalytic combustion characteristics of multi-component VOCs and single-component gases reveals that adding ethane and propane can facilitate methane oxidation.
基金supported by the Key Project of Natural Science Foundation of Tianjin(No.23JCZDJC00330)Tianjin Municipal Education Commission Scientific Research Program(No.2022ZD056).
文摘Prenatal exposure to bisphenols and metals has raised significant concerns regarding their potential impact on fetal development,particularly the risk of fetal chromosome numerical abnormalities(CNA).In this case-control study,we analyzed bisphenol and metal concentrations in amniotic fluid of high-risk pregnant women undergoing amniocentesis.Concentrations of bisphenols and metals were measured using ultra-performance liquid chromatography-tandem mass spectrometry and inductively coupled plasma-mass spectrometry,respectively.Logistic regression and quantile-based g-computation were applied to evaluate individual and combined effects,while dose-response relationships were assessed using restricted cubic splines.Our findings indicated that bisphenol S(BPS),bisphenol Z(BPZ),bisphenol AF(BPAF),antimony(Sb),and vanadium(V)were significantly associated with an increased risk of CNA when analyzed individually,whereas manganese,iron,copper(Cu),nickel(Ni),and zinc(Zn)were significantly and inversely associated with CNA risk.Combined exposure to bisphenol and metal mixtures was associated with an increased risk of CNA in multi-pollutant models.Cu and Ni exhibited a positive additive interaction.Furthermore,BPS,BPZ,and BPAF were individually associated with an increased risk of Down syndrome,while Zn was associated with a decreased risk of Down syndrome.BPS,Sb,V,and Zn were individually associated with an increased risk of Klinefelter syndrome.These findings underscore the potential role of prenatal bisphenol and metal exposure in the pathogenesis of fetal CNA,highlighting both additive and synergistic effects.
基金supported by China Postdoctoral Science Foundation(Grant No.2024M752865)Postdoctoral Fellowship Program of CPSF(Grant No.GZC20241531)+2 种基金Shandong Provincial Higher Education Institutions Youth Plan Team(2022KJ081)the Double First-Class Discipline Construction Fund Project of Harbin Institute of Technology at Weihai(2023SYLCB04)the Open Funding of the Research Center of Civil,Hydraulic and Power Engineering of Xizang(XZA202405CHP2002B).
文摘This study presents a systematic numerical analysis of wind loads on offshore photovoltaic(PV)panels.A computational fluid dynamics(CFD)model,incorporating a free-surface wave boundary condition,is developed and validated against experimental data.Parametric investigations quantify the effects of wind speed,panel tilt angle,clearance,and wave characteristics on the aerodynamic coefficients(drag,lift,and moment).Results indicate that all force coefficients increase with wind speed,with the lift coefficient being most sensitive to wave action.While a larger tilt angle intensifies airflow disturbance and amplifies the coefficients,this effect is more pronounced over flat ground than above a wavy surface.As clearance increases,the drag coefficient fluctuates before rising,the lift coefficient exhibits a trough-shaped response,and the moment coefficient increases monotonically,with values consistently higher over waves.Furthermore,the aerodynamic coefficients generally decrease with greater wave height.The maximum wind load occurs directly above the wave trough,and the aerodynamic force coefficient varies non-monotonically with wave position,first decreasing and then increasing.These findings offer practical guidance for the structural design and safety assurance of offshore PV systems.
文摘Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.
文摘Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagation behavior.To address the unclear mechanisms governing fracture penetration across coal-gangue interfaces,this study employs the Continuum-Discontinuum Element Method(CDEM)to simulate and analyze the vertical propagation of hydraulic fractures initiating within coal seams,based on geomechanical parameters derived from the deep Benxi Formation coal seams in the southeastern Ordos Basin.The investigation systematically examines the influence of geological and operational parameters on cross-interfacial fracture growth.Results demonstrate that vertical stress difference,elastic modulus contrast between coal and gangue layers,interfacial stress differential,and interfacial cohesion at coal-gangue interfaces are critical factors governing hydraulic fracture penetration through these interfaces.High vertical stress differences(>3 MPa)inhibit interfacial dilation,promoting predominant crosslayer fracture propagation.Reduced interfacial stress contrasts and enhanced interfacial cohesion facilitate fracture penetration across interfaces.Furthermore,smaller elastic modulus contrasts between coal and gangue correlate with increased interfacial aperture.Finally,lower injection rates effectively suppress vertical fracture propagation in deep coal reservoirs.This study elucidates the characteristics and mechanisms governing cross-layer fracture propagation in coal–rock composites with interbedded partings,and delineates the dynamic evolution laws and dominant controlling factors involved.Thefindings provide critical theoretical insights for the optimization of fracture design and the efficient development of deep coalbed methane reservoirs.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金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 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.
基金supported by the National Natural Science Foundation of China(Grant Nos.12272393 and 52130905).
文摘As binary geological media,soil-rock mixtures(SRMs)exhibit a distinct gradational composition,leading to their unique mechanical behaviors.To appraise the stability of SRM slopes,it is essential to determine equivalent parameters of SRMs,which are typically obtained through experimental and numerical methods.In contrasted to other numerical methods,the numerical manifold method(NMM)is more effective in addressing SRM problems.This is because the high-precision regular mathematical meshes in NMM can be used without aligning with the soil-rock interfaces and boundaries of SRMs.In the current research,the equivalent strength parameters of SRMs,i.e.the equivalent cohesion ce and internal friction angleϕ_(e),are determined using NMM.Initially,an NMM triaxial numerical model is established and validated based on triaxial experiments.Subsequently,the soil and rock parameters are derived through parameter inversion.Moreover,the impacts of rock content,size,shape and rock blocks'major-axis orientation on ce andϕ_(e) of SRMs are thoroughly examined using the NMM triaxial numerical model.Additionally,a fitting function is proposed to linkϕ_(e) to the rock content and size of SRMs.When other influencing factors are fixed,the above fitting model leads to the following conclusions:(1)the predictedϕ_(e) of SRMs increase with the increase of rock content;and(2)SRM samples with smaller rocks display a higher predictedϕ_(e).
基金supported by Supported by the Scientific Research Foundation for High-Level Talents of Zhoukou Normal University(ZKNUC2024018).
文摘Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.
基金funded by the National Natural Science Foundation of China(Grant No.42102169)the Shaanxi Province key Research and Development Project(Grant No.2023-ZDLSF-64)+1 种基金the Scientific Research and Technological Development Project of China National Logging Corporation(Grant No.25ZYCJSG013-2504)the Youth Science and Technology Special Fund of PetroChina(Grant No.2024DQ03172)。
文摘Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO_(2) trapping mechanisms.In the base case,a constant injection rate of 3500 m^(3)/d over fifteen years resulted in a cumulative injection of 19.2×10^(6) m^(3).The CO_(2) plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO_(2) in the reservoir was 10.6 million tonnes(Mt),representing less than 10% of its theoretical capacity.The plume was projected to reach the entrapment crest and the top of the reservoir within a century,indicating secure long-term containment.Structural,stratigraphic,and residual trapping dominate in Reservoir A(approximately 90%).Anticlinal closures with thick overlying mudstones in the Zhujiang Formation provide effective seals,further enhancing storage security.Reservoir properties and heterogeneity play a crucial role in controlling CO_(2) storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO_(2) capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.
基金supported by Beijing Natural Science Foundation(L233025)。
文摘Functionally graded material(FGM)plates are widely used in various engineering structures owing to their tailor-made mechanical properties,whereas cracked homogeneous plates constitute a canonical setting in fracture mechanics analysis.These two classes of problems respectively embody material non-uniformity and geometric discontinuity,thereby imposing more stringent requirements on numerical methods in terms of high-order field continuity and accurate defect representation.Based on the classical Kirchhoff-Love plate theory,a numerical manifold method(MLS-NMM)incorporating moving least squares(MLS)interpolation is developed for bending analysis of FGM plates and fracture simulation of homogeneous plates with defects.The method constructs an H^(2)-regular approximation with high-order continuous weighting functions and,combined with the separation of mathematical and physical covers,establishes a unified framework that accurately handles material gradients and cracks without mesh reconstruction.For the crack tip,a singular physical cover incorporating the Williams asymptotic field is introduced to achieve local enrichment,enabling the natural capture of displacement discontinuity and stress singularity.Stress intensity factors are extracted using the interaction integral method,and the dimensionless J-integral shows a maximum relative error below 1.2%compared with the reference solution.Numerical results indicate that MLS-NMM exhibits excellent convergence performance:using 676 mathematical nodes,the nondimensional central deflection of both FGM and homogeneous plates agrees with reference solutions with a maximum relative error below 0.81%,and no shear locking occurs.A systematic analysis reveals that for a simply supported on all four edges(SSSS)FGM square plate with a/h=10,the nondimensional central deflection increases by 212%as the gradient index nrises from 0 to 5.For a homogeneous plate containing a central crack with c/a=0.6,the nondimensional central deflection increases by approximately 46%compared with the intact plate.Under weak boundary constraints(e.g.,SFSF),the deformation is markedly amplified,with the deflection reaching more than three times that under strong constraints(SCSC).The proposed method provides an efficient,reconstruction-free numerical tool for high-accuracy bending and fracture analyses of FGM and cracked thin-plate structures.
基金supported by the National Natural Science Foundation of China(No.11872211)。
文摘Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteristics and flow field evolution laws of overlapping rotor configurations in hovering conditions through numerical simulation methods.The research method involves constructing a computational model for rotor flow fields and aerodynamic characteristics based on the Reynolds-averaged Navier-Stokes(RANS)equations and the Spalart-Allmaras(S-A)turbulence model.The dynamic simulation of rotor rotational motion was achieved by using the moving nested grid technology.The reliability of the computational method was ensured through the grid independence verification and the comparison with experimental data.The research results indicate that in overlapping rotor systems,rotorⅡexperiences a decrease in thrust,significant power fluctuations,and reduced hovering efficiency due to continuous interference from the adjacent rotor’s wake and blade-vortex interactions.Blade-tip vortices undergo breakage,fusion,and secondary rolling in the overlapping region,forming large-scale turbulent structures that lead to attenuation of the induced velocity field and aerodynamic efficiency losses.Additionally,the interaction between the rotor downwash and the fuselage triggers a“fountain effect”and a sudden increase in surface pressure on the fuselage,exacerbating flow field distortion.Based on the aforementioned mechanisms,the safe flight of overlapping rotor configurations can be achieved by optimizing the configuration strategy of the rotational speed phase difference between adjacent blades.This study provides a theoretical basis for the rotor layout design and the aerodynamic performance enhancement of heavy-load eVTOL aircraft.
