This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
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.展开更多
The flow of a tetra-hybrid Casson nanofluid(Al_(2)O_(3)-CuO-TiO_(2)-Ag/H_(2)O)over a nonlinear stretching sheet is investigated.The Buongiorno model is used to account for thermophoresis and Brownian motion,while ther...The flow of a tetra-hybrid Casson nanofluid(Al_(2)O_(3)-CuO-TiO_(2)-Ag/H_(2)O)over a nonlinear stretching sheet is investigated.The Buongiorno model is used to account for thermophoresis and Brownian motion,while thermal radiation is incorporated to examine its influence on the thermal boundary layer.The governing partial differential equations(PDEs)are reduced to a system of nonlinear ordinary differential equations(ODEs)with fully non-dimensional similarity transformations involving all independent variables.To solve the obtained highly nonlinear system of differential equations,a novel Clique polynomial collocation method is applied.The analysis focuses on the effects of the Casson parameter,power index,radiation parameter,thermophoresis parameter,Brownian motion parameter,and Lewis number.The key findings show that thermal radiation intensifies the thermal boundary layer,the Casson parameter reduces the velocity,and the Lewis number suppresses the concentration with direct relevance to polymer processing,coating flows,electronic cooling,and biomedical applications.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and ...Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.展开更多
Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
Sheet metal spinning is an incremental forming process for producing axisymmetric thinwalled parts through continuous local deformation under the action of rollers.While studying the spinning process by finite element...Sheet metal spinning is an incremental forming process for producing axisymmetric thinwalled parts through continuous local deformation under the action of rollers.While studying the spinning process by finite element(FE)method,a critical bottleneck is the enormous simulation time.For beating off this challenge,a novel multi-mesh method is developed.The method can dynamically track the movement of rollers and adaptively refine the mesh.Thus,a locally refined quadrilateral computation mesh can be generated in the locally-deforming zone and reduce the unnecessary fine elements outside the locally-deforming zone.In the multi-mesh system,the fine elements and coarse elements are extracted from a storage mesh and a background mesh,respectively.Meanwhile,the hanging nodes in the locally refined mesh are removed by designing 4-refinement templates.Between computation mesh and storage mesh,a bi-cubic parametric surface fitting algorithm and accurate remapping methods are conducted to transmit geometric information and physical fields.The proposed method has been verified by two spinning processes.The results suggest that the method can save time by up to about 67%with satisfactory accuracy,especially for distributions of thickness and strain compared with the fully refined mesh.展开更多
As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-s...As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are reg...In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.展开更多
The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX a...The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.展开更多
At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures...At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.展开更多
The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional ...An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.展开更多
The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improv...The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.展开更多
This study systematically conducted preparation optimization and performance investigations on Co-modified Ce/TiO_(2) catalysts,with a focus on examining how preparation methods and Co loading regulate the catalyst’s...This study systematically conducted preparation optimization and performance investigations on Co-modified Ce/TiO_(2) catalysts,with a focus on examining how preparation methods and Co loading regulate the catalyst’s low-temperature denitrification activity.After identifying optimal preparation parameters via condition screening,multiple characterization techniques-including BET,XRD,XPS,H_(2)-TPR and in situ DRIFTS-were employed to deeply analyze the catalyst’s physicochemical properties and reaction mechanism.Results demonstrated that compared to the impregnation and co-precipitation methods,the Ce-Co_(0.025)/TiO_(2)-SG catalyst(prepared by the sol-gel method with a Co/Ti mass ratio of 0.025)exhibited significantly superior denitrification activity:NO conversion remained stably above 95%in the 225−350℃ temperature range,and it displayed high N_(2) selectivity.Characterization analysis revealed that abundant surface oxygen vacancies,a high proportion of Ce^(3+) species,and prominent acidic sites collectively contributed to enhancing its low-temperature denitrification performance.This work provides reference value for the development of highly efficient low-temperature denitrification catalysts.展开更多
The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase fiel...The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.展开更多
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
文摘This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
基金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 UGC,New Delhi,India for financial assistance via the UGC-Junior Research Fellowship(CSIR-UGC NET JULY 2024)(Student ID:241610090610)。
文摘The flow of a tetra-hybrid Casson nanofluid(Al_(2)O_(3)-CuO-TiO_(2)-Ag/H_(2)O)over a nonlinear stretching sheet is investigated.The Buongiorno model is used to account for thermophoresis and Brownian motion,while thermal radiation is incorporated to examine its influence on the thermal boundary layer.The governing partial differential equations(PDEs)are reduced to a system of nonlinear ordinary differential equations(ODEs)with fully non-dimensional similarity transformations involving all independent variables.To solve the obtained highly nonlinear system of differential equations,a novel Clique polynomial collocation method is applied.The analysis focuses on the effects of the Casson parameter,power index,radiation parameter,thermophoresis parameter,Brownian motion parameter,and Lewis number.The key findings show that thermal radiation intensifies the thermal boundary layer,the Casson parameter reduces the velocity,and the Lewis number suppresses the concentration with direct relevance to polymer processing,coating flows,electronic cooling,and biomedical applications.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金supported by the National Natural Science Foundation of China(No.42372287 and No.U24A20178)the Fundamental Research Funds for the Central Universities CHD(No.2024SHEEAR002)+3 种基金the Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shaanxi Province(No.2020024)the China Postdoctoral Science Foundation(GZC20232955,2024M753472,and 2024MD763937)the Science-Technology Foundation for Young Scientists of Gansu Province,China(No.24JRRA097)the Study of biodiversity survey and limiting factor analysis of Yinkentala(2023ZL01).
文摘Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
基金co-supported by the supports of Guangdong Basic and Applied Basic Research Foundation(No.2019B1515120047)the National Natural Science Foundation of China(No.52130507)。
文摘Sheet metal spinning is an incremental forming process for producing axisymmetric thinwalled parts through continuous local deformation under the action of rollers.While studying the spinning process by finite element(FE)method,a critical bottleneck is the enormous simulation time.For beating off this challenge,a novel multi-mesh method is developed.The method can dynamically track the movement of rollers and adaptively refine the mesh.Thus,a locally refined quadrilateral computation mesh can be generated in the locally-deforming zone and reduce the unnecessary fine elements outside the locally-deforming zone.In the multi-mesh system,the fine elements and coarse elements are extracted from a storage mesh and a background mesh,respectively.Meanwhile,the hanging nodes in the locally refined mesh are removed by designing 4-refinement templates.Between computation mesh and storage mesh,a bi-cubic parametric surface fitting algorithm and accurate remapping methods are conducted to transmit geometric information and physical fields.The proposed method has been verified by two spinning processes.The results suggest that the method can save time by up to about 67%with satisfactory accuracy,especially for distributions of thickness and strain compared with the fully refined mesh.
基金supported by the National Key Research and Development Program of China(2022YFC2904400)Guangxi Science and Technology Major Project(Gui Ke AA23023033)。
文摘As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金funded by National Nature Science Foundation of China(92266203)National Nature Science Foundation of China(52205278)+1 种基金Key Projects of Shijiazhuang Basic Research Program(241791077A)Central Guide Local Science and Technology Development Fund Project of Hebei Province(246Z1022G).
文摘In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.
基金supported by the National Natural Science Foundation of China(Nos.52105384 and U2141215).
文摘The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.
基金This work was supported by the National Key R&D Program of China(Grant No.2022YFE0207000)the National Natural Science Foundation of China(Grant Nos.12372289,11972250,and 12102298)+1 种基金the China Postdoctoral Science Foundation(Grant No.2021M702443)Tianjin Natural Science Foundation(Grant No.22JCZDJC00910).
文摘At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金co-supported by the Foundation of the State Key Laboratory of Aerodynamics(No.RAL202203)the National Key Laboratory Foundation(No.6142202202)the China Postdoctoral Science Foundation(No.2024M754133)。
文摘An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.
基金National Key R&D Program of China under Grant No.2022YFC3003601。
文摘The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.
基金Supported by the National Key Research and Development Program of China (2023YFB4102903)。
文摘This study systematically conducted preparation optimization and performance investigations on Co-modified Ce/TiO_(2) catalysts,with a focus on examining how preparation methods and Co loading regulate the catalyst’s low-temperature denitrification activity.After identifying optimal preparation parameters via condition screening,multiple characterization techniques-including BET,XRD,XPS,H_(2)-TPR and in situ DRIFTS-were employed to deeply analyze the catalyst’s physicochemical properties and reaction mechanism.Results demonstrated that compared to the impregnation and co-precipitation methods,the Ce-Co_(0.025)/TiO_(2)-SG catalyst(prepared by the sol-gel method with a Co/Ti mass ratio of 0.025)exhibited significantly superior denitrification activity:NO conversion remained stably above 95%in the 225−350℃ temperature range,and it displayed high N_(2) selectivity.Characterization analysis revealed that abundant surface oxygen vacancies,a high proportion of Ce^(3+) species,and prominent acidic sites collectively contributed to enhancing its low-temperature denitrification performance.This work provides reference value for the development of highly efficient low-temperature denitrification catalysts.
基金National Science Foundation of China,Grant/Award Number:5227090113Shenzhen Science and Technology Program,Grant/Award Numbers:KCXFZ20230731093901003,KCXFZ20211020163816023The Project of Hetao Shenzhen-Hong Kong Science and Technology Innovation Cooperation Zone,Grant/Award Number:HZQB-KCZYB-2020083。
文摘The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.
