By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as w...By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as well as duration of spontaneous combustion from large-scale spontaneous combustion experiment, 'three zones' of spontaneous combustion were partitioned and mining conditions to avoid spontaneous combustion were obtained. The above method was employed to partition 'three zones' in gob of fully mechanized top-coal caving long wall face and got fairly good result. Calculation of the above method is much smaller than simulating the whole process of coal spontaneous combustion, but the prediction precision can satisfy the demand of predicting and extinguishing spontaneous combustion in mining.展开更多
It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period (two dates), the attendance of budge...It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period (two dates), the attendance of budgetary restrictions and the conditions of general equilibrium (without production). It uses the technique of selection of candidate points in a simulation process with a portfolio of efficient assets and a hypothesis for the process of determining the returns of the securities: A GARCH process. By this technique, it compares the stochastic volatility patterns between the artificial series, obtained in the simulation, and the real series of household aggregate consumption in the US and Brazil.展开更多
Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the...Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the form of exponent function with radial distance, and considering the nonlinear constitutive relation of stress and strain, the stiffness matrix is established. The mechanics behavior of the pile—soil interface is simulated with a nonlinear interface element. This method can truly express the behavior of the pile-soil system. The load-settlement relation Q-S curves of two big diameter prototype piles on bearing test are analyzed, and satisfying results are obtained. This method is reasonable in theory and feasible in engineering.展开更多
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.展开更多
With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity ...With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity achievement during drilling have become inevitable.These challenges necessitate a technology that can accurately detect adjacent wells in real time during drilling operations.As current borehole acoustic reflection imaging technology heavily relies on cable-based logging,it cannot perform real-time detection of adjacent wells during drilling,thereby limiting the drilling efficiency.This study proposes a new adjacent-well-acoustic-detection-while-drilling method that integrates wireline borehole acoustic reflection imaging with drilling technology,along with an adjacent-well imaging method based on compressed sensing(CS).Together,these methods enable high-resolution,real-time detection of the adjacent target wells during drilling,ensuring safe and efficient underground drilling operations.The finite-difference method was used to simulate three-dimensional numerical models under drilling conditions for two scenarios—with and without target wells adjacent to the drilling well.Experimental validation was conducted in a water tank using an adjacent-well-acoustic-detection-while-drilling tool.The simulated target well was imaged using the CS method,and the imaging results were compared with those obtained from numerical and physical simulations,thereby validating the feasibility of the proposed acoustic detection and imaging methods.The results demonstrate that as the radial distance from the target well increases,the PP echo exhibits delayed arrival times and approaches a plane wave while exhibiting amplitude attenuation.Conversely,a linear increase in the target well diameter advances the PP echo arrival time and enhances its amplitude proportionally.When the target and drilling wells are approximately parallel with a small intersection angle,PP echoes yield better detection results than SS echoes;when the wells are coplanar with a large intersection angle,SS echoes provide better detection results.The receiver element aligned with the target well's azimuth detects all echo modes with the earliest arrival times and highest amplitudes.The adjacent-well imaging method based on CS offers very high spatial resolution,with target wells appearing as local amplitude maxima.This feature enables the precise determination of their azimuth and inclination relative to the drilling wells.The findings offer a solid physical and methodological foundation for real-time detection of adjacent wells during drilling operations and demonstrate enormous theoretical and engineering application potential.展开更多
Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patt...Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patterns around bridge-road transitions under different wind directions remain limited.In this study,numerical simulations are employed to study the influence of wind angle on the flow field structure,wind profiles,horizontal wind speed distribution,and sand deposition-erosion behavior around the bridge-road transition.The results indicate that as the wind angle increases,the leeward vortex recirculation zone expands,the peak of the vertical wind speed profile increases,and the horizontal wind speed decreases.With respect to deposition-erosion patterns,the sand deposition length on the windward sides of both the bridge pier and the bridge-road transition decreases with increasing wind angle.In contrast,the deposition length on the leeward side increases.Specifically,when the wind angle increases from 60°to 90°,the windward deposition length at the bridge-road transition decreases from 1.99H to 1.49H(H,the clear height of bridge),while the leeward deposition expands from 0.70H to 13.28H.In addition,the bridge-road transition section serves as a channel for sand transport,promoting the migration of sand toward the subgrade.A smaller wind angle enhances the lateral guidance effect on the wind-sand flow,resulting in a greater sand deposition length on the leeward side of the subgrade.Therefore,the prevention and control of sand hazards in the bridge-road transition section should not be neglected.It is necessary to regularly remove accumulated sand to prevent its spread from the transition section to the subgrade,which could threaten traffic safety.展开更多
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.展开更多
Geological CO2 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 CO2 injection and stora...Geological CO2 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 CO2 injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO2 trapping mechanisms.In the base case,a constant injection rate of 3500 m3/d over fifteen years resulted in a cumulative injection of 19.2×106 m3.The CO2 plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO2 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 CO2 storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO2 capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.展开更多
The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact ...The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact of kinematic parameters remains relatively unexplored.This study focuses on the swimming posture and movement characteristics of octopuses,formulating a synchronized flexible undulation equation for their eight arms.Based on this,computational fluid dynamics combined with dynamic mesh technology is employed to numerically simulate their propulsion mechanism and swimming characteristics.Furthermore,by varying the hover coast time ratio and duty cycle,the influence of different kinematic parameters on propulsion efficiency is explored(close phase after hover coast time ratio is defined as DSC;open phase after hover coast time ratio is defined as DSO,and the duty cycle is defined as DC).The results indicate that the periodic opening and closing motion of the octopus leads to the alternating generation and dissipation of reverse vortices in the wake field,with the jet effect between them being the primary source of propulsion.When DSC=0.30,the forward speed of the octopus after 1 s is only 6%lower than that without hover coast behavior.When DSO=0.30,the peak thrust coefficient of the octopus reaches 0.82,which is 15%higher than that without hover coast behavior,indicating the highest burst acceleration.On the other hand,when DC=0.67,under asymmetric periodic swimming with fast-closing-slow-opening,the average thrust coefficient reaches 0.26,an increase of 23.8%compared to uniform opening-closing periodic swimming,while the average lateral force coefficient decreases by 66.7%.In this scenario,the octopus can achieve better forward propulsion and maintain stable movement.The findings of this study provide valuable insights for research on the intermittent swimming behavior of octopuses and underwater vehicles.展开更多
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.展开更多
This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Adve...This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Advector interface reconstruction technique,atomization characteristics are simulated and analyzed for different values of We.The results indicate that the geometry induces turbulent jets,which drive turbulent atomization through the shear interactions at the gas–liquid interface.The key observed phenomena include the interaction of impact waves with liquid sheet perforation and the breakup of web of ligaments,both of which are prominent under high backpressure conditions.A novel method,based on the threshold velocity of spray droplet groups,is employed to quantitatively measure the spreading angle,showing that the angle increases with We in both front and side views.Additionally,the Sauter mean diameter of droplets follows power-law scaling with exponents of-1/3 in the upstream region and-1/2 in the downstream region,while the droplet size distribution conforms to a log-normal profile.This research provides valuable insights into interface evolution and droplet characteristics during impingingjet atomization under high backpressure,offering essential guidance for optimizing industrial atomization processes.展开更多
With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving larg...With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.展开更多
The pot cover effect can induce various forms of distress in cover layer engineering,such as salt heave,cracking,and differential settlement,with water vapor migration being the primary cause.However,current research ...The pot cover effect can induce various forms of distress in cover layer engineering,such as salt heave,cracking,and differential settlement,with water vapor migration being the primary cause.However,current research on the pot cover effect in saline soils rarely takes into account the water vapor transport process.Therefore,elucidating the coupled transport mechanisms of water,vapor,heat,and salt in saline soils under this effect is crucial for the prevention and control of related engineering hazards.This study developed a numerical model describing the coupled water-vapor-heat-salt transport in unsaturated saline sulfate soil and validated its reliability through laboratory unidirectional freezing column tests.Based on this model,a numerical analysis was conducted to investigate the formation mechanism of the pot cover effect during the unidirectional freezing of the saline soil.The results indicate that the moisture and salt fields exhibit a typical bimodal distribution pattern,with peaks located at the soil surface and the freezing front,respectively.Compared with the initial water content of 19% and initial salt content of 1%,the total water content at the surface and freezing front increased by 21% and 13%,respectively,while the total salt content rose by 1.25% and 0.5%,respectively.Liquid water flux upward in the unfrozen zone,while it approaches zero within the frozen zone.In contrast,both vapor flux and solute flux migrate upward throughout the entire soil column,reaching their maximum values at the freezing front.Compared to models that neglect vapor transport,the simulated total moisture content at the surface was 12% higher in the model accounting for vapor movement,indicating that water vapor migration is a key factor contributing to moisture accumulation at the surface.The findings of this study can provide a theoretical basis for preventing engineering hazards associated with the pot cover effect in saline soils.展开更多
Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrus...Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrust,in addition to precisely controlling the separation attitude.To address these issues,this paper introduces a multi-gasbag propulsion system with a high aspect ratio that coordinates multiple gasbags to generate sufficient thrust.By adjusting the inlet size of the gasbag,the separation behavior of the release unit can be accurately controlled.A multidimensional two-phase flow model is established,accompanied by both combustion and flow experiments and a double-gasbag propulsion experiment.The results demonstrate that the proposed mathematical model is accurate,effectively captures the pressure fluctuations and spatiotemporal distribution of flow field parameters,and determines the separation attitude of the release unit.For the cases studied in this paper,the pressure at the gasbag inlet(z=650 mm)is the dominant factor during the gasbag propulsion response,causing the release unit to rotate counterclockwise when the gasbag inlet sizes are identical.Increasing the inlet size at z=50 mm compensates for the adverse effects of uneven axial pressure distribution,thereby achieving a neutral separation for the release unit.When the radii r1 and r2 vary between 2 and 12 mm,the angular velocity and attitude angle of the release unit are found to range from-15.50 to 15.20 rad/s and from-0.109 to 0.106 rad,respectively.展开更多
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.展开更多
Direct numerical simulation of spatially developing turbulent boundary layers with periodic blowing or suction through a series of inclined slots in a control region was conducted.The wall-generated Reynolds shear str...Direct numerical simulation of spatially developing turbulent boundary layers with periodic blowing or suction through a series of inclined slots in a control region was conducted.The wall-generated Reynolds shear stress(RSS),i.e.,the RSS on the wall in the control region,is generated through the control scheme.The effect of the wall-generated RSS on friction drag was examined via a series of simulations.The Reynolds numbers of the turbulent boundary layers investigated vary from 300-860 on the basis of the external flow velocity and the momentum thickness.The proposed control scheme was used to verify the relationship between the wall-generated RSS and skin friction drag,and it was found that a wall-generated negative RSS(net positive)increases the skin friction drag,whereas a wall-generated positive RSS(net negative)reduces it.The proposed control method can provide a high drag reduction rate,and even negative resistance and backflow can be observed.展开更多
Gas condensate reservoirs constitute important natural gas resources;however,their development is frequently hindered by condensate banking and complex multiphase flow behavior.Naturally fractured gas condensate reser...Gas condensate reservoirs constitute important natural gas resources;however,their development is frequently hindered by condensate banking and complex multiphase flow behavior.Naturally fractured gas condensate reservoirs present additional challenges because their dualporosity and dual-permeability structure induces strong phase redistribution and nonuniform flow between matrix and fracture systems,thereby complicating reservoir characterization and compositional simulation.In this study,integrated laboratory experiments and numerical simulations were performed for a deep,rich,naturally fractured gas condensate reservoir.Depletion,diffusion,and core flooding experiments involving CO2,N2,and dry gas injection were conducted using fractured core samples.A dual-porosity and dual-permeability compositional model incorporating a five-spot well pattern was established to evaluate condensate liquid recovery and to quantify mass transfer between matrix and fracture networks.The effect of matrix-fracture permeability contrast on production performance was systematically analyzed.The results indicate that matrix permeability is a primary parameter controlling recovery in gas condensate reservoirs.The ratio of matrix-fracture permeability contrasts exerts a stronger influence on condensate liquid recovery than on natural gas recovery.Pressure maintenance through gas injection is critical for improving recovery performance.When reservoir pressure declines below the dew-point pressure,early gas injection is recommended to mitigate condensate accumulation in the near-well region.Among the injected gases evaluated,CO2 demonstrated superior pressure maintenance performance compared with N2 and dry gas.展开更多
Soft soil is widely distributed and has complex origins.Shield tunnels are inevitably constructed within soft soil interlayers,and under seismic action,tunnels may be subject to severe damages.On the basis of actual e...Soft soil is widely distributed and has complex origins.Shield tunnels are inevitably constructed within soft soil interlayers,and under seismic action,tunnels may be subject to severe damages.On the basis of actual engineering,this study utilized dynamic triaxial testing to investigate the dynamic properties of soft soil.Using the PIMY constitutive model,the seismic subsidence characteristics of the soft soil were characterized,and a refined finite element model was established to study the mutual influence mechanism between the soft soil layer and shield tunnels via the open-source software framework OpenSees.The results demonstrate that soil exhibits a softening effect under dynamic loading;soft soil with better structural integrity is less prone to seismic subsidence;and the greater the inertial force acting on the soft soil,the greater the likelihood of settlement.Under seismic action,the presence of the shield tunnel exacerbates the settlement of the soft soil,as the surrounding soil experiences significant inertial forces from the tunnel structure,hindering drainage and accelerating the accumulation of pore water pressure;The soft soil itself has large deformation and displacement under the action of earthquake,which leads to the great stress,deformation,and displacement of the structure.The arch foot position of the tunnel is identified as the most vulnerable to damage.展开更多
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 predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simu...The predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simulated by a modified Leslie-Gower predator-prey model.Using the theory of reaction-diffusion equations,a priori estimate,existence,uniqueness and stability conditions of the positive steady state solution are established.Furthermore,numerical simulations are conducted to quantitatively analyze the dynamical behavior.The research shows that as long as the Allee effect constant satisfies the appropriate relationship and the growth rates of predator and prey are appropriately large,the predator and prey can not only coexist,but also the coexistence mode is unique and stable under low predation-rate.In addition,the numerical simulations show that the coexistence may be stable under high predation-rate.Meanwhile,with the increase of predation rate,the population density of predators will decrease.展开更多
基金Supported by Natural Science Program of Shaanxi Province Education Department (05JK261)
摘要By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as well as duration of spontaneous combustion from large-scale spontaneous combustion experiment, 'three zones' of spontaneous combustion were partitioned and mining conditions to avoid spontaneous combustion were obtained. The above method was employed to partition 'three zones' in gob of fully mechanized top-coal caving long wall face and got fairly good result. Calculation of the above method is much smaller than simulating the whole process of coal spontaneous combustion, but the prediction precision can satisfy the demand of predicting and extinguishing spontaneous combustion in mining.
摘要It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period (two dates), the attendance of budgetary restrictions and the conditions of general equilibrium (without production). It uses the technique of selection of candidate points in a simulation process with a portfolio of efficient assets and a hypothesis for the process of determining the returns of the securities: A GARCH process. By this technique, it compares the stochastic volatility patterns between the artificial series, obtained in the simulation, and the real series of household aggregate consumption in the US and Brazil.
摘要Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the form of exponent function with radial distance, and considering the nonlinear constitutive relation of stress and strain, the stiffness matrix is established. The mechanics behavior of the pile—soil interface is simulated with a nonlinear interface element. This method can truly express the behavior of the pile-soil system. The load-settlement relation Q-S curves of two big diameter prototype piles on bearing test are analyzed, and satisfying results are obtained. This method is reasonable in theory and feasible in engineering.
基金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.
基金supported in part by the National Natural Science Foundation of China under Grant Nos.12334019,12274465,12504558 and U25B20244in part by the China Postdoctoral Science Foundation under Grant No.2025M770469+1 种基金in part by the Postdoctoral Fellowship Program of CPSF under Grant No.GZC20251952in part by the Science Foundation of China University of Petroleum,Beijing under Grant No.2462025XKBH014.
摘要With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity achievement during drilling have become inevitable.These challenges necessitate a technology that can accurately detect adjacent wells in real time during drilling operations.As current borehole acoustic reflection imaging technology heavily relies on cable-based logging,it cannot perform real-time detection of adjacent wells during drilling,thereby limiting the drilling efficiency.This study proposes a new adjacent-well-acoustic-detection-while-drilling method that integrates wireline borehole acoustic reflection imaging with drilling technology,along with an adjacent-well imaging method based on compressed sensing(CS).Together,these methods enable high-resolution,real-time detection of the adjacent target wells during drilling,ensuring safe and efficient underground drilling operations.The finite-difference method was used to simulate three-dimensional numerical models under drilling conditions for two scenarios—with and without target wells adjacent to the drilling well.Experimental validation was conducted in a water tank using an adjacent-well-acoustic-detection-while-drilling tool.The simulated target well was imaged using the CS method,and the imaging results were compared with those obtained from numerical and physical simulations,thereby validating the feasibility of the proposed acoustic detection and imaging methods.The results demonstrate that as the radial distance from the target well increases,the PP echo exhibits delayed arrival times and approaches a plane wave while exhibiting amplitude attenuation.Conversely,a linear increase in the target well diameter advances the PP echo arrival time and enhances its amplitude proportionally.When the target and drilling wells are approximately parallel with a small intersection angle,PP echoes yield better detection results than SS echoes;when the wells are coplanar with a large intersection angle,SS echoes provide better detection results.The receiver element aligned with the target well's azimuth detects all echo modes with the earliest arrival times and highest amplitudes.The adjacent-well imaging method based on CS offers very high spatial resolution,with target wells appearing as local amplitude maxima.This feature enables the precise determination of their azimuth and inclination relative to the drilling wells.The findings offer a solid physical and methodological foundation for real-time detection of adjacent wells during drilling operations and demonstrate enormous theoretical and engineering application potential.
基金supported by the Major Science and Technology Projects of Inner Mongolia Autonomous Region‘Open Bidding for Selecting the Best Candidates’(2024JBGS0009-01)the National Natural Science Foundation of China(42461011)the Key Research and Development Program of Gansu Province(25YFGA040).
摘要Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patterns around bridge-road transitions under different wind directions remain limited.In this study,numerical simulations are employed to study the influence of wind angle on the flow field structure,wind profiles,horizontal wind speed distribution,and sand deposition-erosion behavior around the bridge-road transition.The results indicate that as the wind angle increases,the leeward vortex recirculation zone expands,the peak of the vertical wind speed profile increases,and the horizontal wind speed decreases.With respect to deposition-erosion patterns,the sand deposition length on the windward sides of both the bridge pier and the bridge-road transition decreases with increasing wind angle.In contrast,the deposition length on the leeward side increases.Specifically,when the wind angle increases from 60°to 90°,the windward deposition length at the bridge-road transition decreases from 1.99H to 1.49H(H,the clear height of bridge),while the leeward deposition expands from 0.70H to 13.28H.In addition,the bridge-road transition section serves as a channel for sand transport,promoting the migration of sand toward the subgrade.A smaller wind angle enhances the lateral guidance effect on the wind-sand flow,resulting in a greater sand deposition length on the leeward side of the subgrade.Therefore,the prevention and control of sand hazards in the bridge-road transition section should not be neglected.It is necessary to regularly remove accumulated sand to prevent its spread from the transition section to the subgrade,which could threaten traffic safety.
摘要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.
基金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 CO2 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 CO2 injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO2 trapping mechanisms.In the base case,a constant injection rate of 3500 m3/d over fifteen years resulted in a cumulative injection of 19.2×106 m3.The CO2 plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO2 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 CO2 storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO2 capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.
基金supported by the National Key R&D Program of China(Grant No.2024YFD2400200).
摘要The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact of kinematic parameters remains relatively unexplored.This study focuses on the swimming posture and movement characteristics of octopuses,formulating a synchronized flexible undulation equation for their eight arms.Based on this,computational fluid dynamics combined with dynamic mesh technology is employed to numerically simulate their propulsion mechanism and swimming characteristics.Furthermore,by varying the hover coast time ratio and duty cycle,the influence of different kinematic parameters on propulsion efficiency is explored(close phase after hover coast time ratio is defined as DSC;open phase after hover coast time ratio is defined as DSO,and the duty cycle is defined as DC).The results indicate that the periodic opening and closing motion of the octopus leads to the alternating generation and dissipation of reverse vortices in the wake field,with the jet effect between them being the primary source of propulsion.When DSC=0.30,the forward speed of the octopus after 1 s is only 6%lower than that without hover coast behavior.When DSO=0.30,the peak thrust coefficient of the octopus reaches 0.82,which is 15%higher than that without hover coast behavior,indicating the highest burst acceleration.On the other hand,when DC=0.67,under asymmetric periodic swimming with fast-closing-slow-opening,the average thrust coefficient reaches 0.26,an increase of 23.8%compared to uniform opening-closing periodic swimming,while the average lateral force coefficient decreases by 66.7%.In this scenario,the octopus can achieve better forward propulsion and maintain stable movement.The findings of this study provide valuable insights for research on the intermittent swimming behavior of octopuses and underwater vehicles.
基金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.
基金partly supported by the National Natural Science Foundation of China(Nos.U23B6009 and 12272050)。
摘要This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Advector interface reconstruction technique,atomization characteristics are simulated and analyzed for different values of We.The results indicate that the geometry induces turbulent jets,which drive turbulent atomization through the shear interactions at the gas–liquid interface.The key observed phenomena include the interaction of impact waves with liquid sheet perforation and the breakup of web of ligaments,both of which are prominent under high backpressure conditions.A novel method,based on the threshold velocity of spray droplet groups,is employed to quantitatively measure the spreading angle,showing that the angle increases with We in both front and side views.Additionally,the Sauter mean diameter of droplets follows power-law scaling with exponents of-1/3 in the upstream region and-1/2 in the downstream region,while the droplet size distribution conforms to a log-normal profile.This research provides valuable insights into interface evolution and droplet characteristics during impingingjet atomization under high backpressure,offering essential guidance for optimizing industrial atomization processes.
基金Project(2025ZD1009704)supported by the National Science and Technology Major Project of ChinaProjects(2023JJ30659,2022JJ30706)supported by Hunan Provincial Natural Science Foundation,China。
摘要With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.
基金supported by the National Natural Science Foundation of China(Grant Nos.52168052,12362032,52168052)the Gansu Provincial University Industry Support Program(Grant No.2025CYZC-033)+1 种基金the Science Foundation for Distinguished Young Scholars of Gansu Province(Grant No.24JRRA167)the Gansu Province Longyuan Young Talents Project(Zhang Mingli,2025)。
摘要The pot cover effect can induce various forms of distress in cover layer engineering,such as salt heave,cracking,and differential settlement,with water vapor migration being the primary cause.However,current research on the pot cover effect in saline soils rarely takes into account the water vapor transport process.Therefore,elucidating the coupled transport mechanisms of water,vapor,heat,and salt in saline soils under this effect is crucial for the prevention and control of related engineering hazards.This study developed a numerical model describing the coupled water-vapor-heat-salt transport in unsaturated saline sulfate soil and validated its reliability through laboratory unidirectional freezing column tests.Based on this model,a numerical analysis was conducted to investigate the formation mechanism of the pot cover effect during the unidirectional freezing of the saline soil.The results indicate that the moisture and salt fields exhibit a typical bimodal distribution pattern,with peaks located at the soil surface and the freezing front,respectively.Compared with the initial water content of 19% and initial salt content of 1%,the total water content at the surface and freezing front increased by 21% and 13%,respectively,while the total salt content rose by 1.25% and 0.5%,respectively.Liquid water flux upward in the unfrozen zone,while it approaches zero within the frozen zone.In contrast,both vapor flux and solute flux migrate upward throughout the entire soil column,reaching their maximum values at the freezing front.Compared to models that neglect vapor transport,the simulated total moisture content at the surface was 12% higher in the model accounting for vapor movement,indicating that water vapor migration is a key factor contributing to moisture accumulation at the surface.The findings of this study can provide a theoretical basis for preventing engineering hazards associated with the pot cover effect in saline soils.
基金supported by the National Natural Science Foundation of China(Grant No.52406186)。
摘要Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrust,in addition to precisely controlling the separation attitude.To address these issues,this paper introduces a multi-gasbag propulsion system with a high aspect ratio that coordinates multiple gasbags to generate sufficient thrust.By adjusting the inlet size of the gasbag,the separation behavior of the release unit can be accurately controlled.A multidimensional two-phase flow model is established,accompanied by both combustion and flow experiments and a double-gasbag propulsion experiment.The results demonstrate that the proposed mathematical model is accurate,effectively captures the pressure fluctuations and spatiotemporal distribution of flow field parameters,and determines the separation attitude of the release unit.For the cases studied in this paper,the pressure at the gasbag inlet(z=650 mm)is the dominant factor during the gasbag propulsion response,causing the release unit to rotate counterclockwise when the gasbag inlet sizes are identical.Increasing the inlet size at z=50 mm compensates for the adverse effects of uneven axial pressure distribution,thereby achieving a neutral separation for the release unit.When the radii r1 and r2 vary between 2 and 12 mm,the angular velocity and attitude angle of the release unit are found to range from-15.50 to 15.20 rad/s and from-0.109 to 0.106 rad,respectively.
摘要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(Grant Nos.12362022 and 11602123)the Natural Science Foundation of Inner Mongolia(Grant Nos.2020MS01012 and 2022QN05016)+2 种基金the Scientific Research Projects of Colleges and Universities in Inner Mongolia(Grant No.NJZY21158)the Project of Ordos City Key Innovation Team(R&D Team for Drum-type Small and Micro Gas Turbines)the Key Special Fund of Inner Mongolia Autonomous Region's Science and Technology Revitalization Initiative(Grant No.2021EEDSCXSFQZD009)。
摘要Direct numerical simulation of spatially developing turbulent boundary layers with periodic blowing or suction through a series of inclined slots in a control region was conducted.The wall-generated Reynolds shear stress(RSS),i.e.,the RSS on the wall in the control region,is generated through the control scheme.The effect of the wall-generated RSS on friction drag was examined via a series of simulations.The Reynolds numbers of the turbulent boundary layers investigated vary from 300-860 on the basis of the external flow velocity and the momentum thickness.The proposed control scheme was used to verify the relationship between the wall-generated RSS and skin friction drag,and it was found that a wall-generated negative RSS(net positive)increases the skin friction drag,whereas a wall-generated positive RSS(net negative)reduces it.The proposed control method can provide a high drag reduction rate,and even negative resistance and backflow can be observed.
基金the National Natural Science Foundation of China(NO.52004032).
摘要Gas condensate reservoirs constitute important natural gas resources;however,their development is frequently hindered by condensate banking and complex multiphase flow behavior.Naturally fractured gas condensate reservoirs present additional challenges because their dualporosity and dual-permeability structure induces strong phase redistribution and nonuniform flow between matrix and fracture systems,thereby complicating reservoir characterization and compositional simulation.In this study,integrated laboratory experiments and numerical simulations were performed for a deep,rich,naturally fractured gas condensate reservoir.Depletion,diffusion,and core flooding experiments involving CO2,N2,and dry gas injection were conducted using fractured core samples.A dual-porosity and dual-permeability compositional model incorporating a five-spot well pattern was established to evaluate condensate liquid recovery and to quantify mass transfer between matrix and fracture networks.The effect of matrix-fracture permeability contrast on production performance was systematically analyzed.The results indicate that matrix permeability is a primary parameter controlling recovery in gas condensate reservoirs.The ratio of matrix-fracture permeability contrasts exerts a stronger influence on condensate liquid recovery than on natural gas recovery.Pressure maintenance through gas injection is critical for improving recovery performance.When reservoir pressure declines below the dew-point pressure,early gas injection is recommended to mitigate condensate accumulation in the near-well region.Among the injected gases evaluated,CO2 demonstrated superior pressure maintenance performance compared with N2 and dry gas.
基金funded by the Key Research and Development Projects of Guangdong Province(2019B111108001)China and the Guangdong Basic and Applied Basic Research Foundation(2024A1515010045)+3 种基金China and the Basic Research Program of Jiangsu(BK20231217 and BK20220265)the independent research project of the State Key Laboratory of Subtropical Building and Urban Science(2023ZB15)China the Key Laboratory of Geomechanics and Geotechnical Engineering Safety,the Chinese Academy of Sciences(SKLGME023001)China and the Open Research Fund Program of the State Key Laboratory of Hydroscience and Engineering(SKLGME-KF-2025-D-02),China。
摘要Soft soil is widely distributed and has complex origins.Shield tunnels are inevitably constructed within soft soil interlayers,and under seismic action,tunnels may be subject to severe damages.On the basis of actual engineering,this study utilized dynamic triaxial testing to investigate the dynamic properties of soft soil.Using the PIMY constitutive model,the seismic subsidence characteristics of the soft soil were characterized,and a refined finite element model was established to study the mutual influence mechanism between the soft soil layer and shield tunnels via the open-source software framework OpenSees.The results demonstrate that soil exhibits a softening effect under dynamic loading;soft soil with better structural integrity is less prone to seismic subsidence;and the greater the inertial force acting on the soft soil,the greater the likelihood of settlement.Under seismic action,the presence of the shield tunnel exacerbates the settlement of the soft soil,as the surrounding soil experiences significant inertial forces from the tunnel structure,hindering drainage and accelerating the accumulation of pore water pressure;The soft soil itself has large deformation and displacement under the action of earthquake,which leads to the great stress,deformation,and displacement of the structure.The arch foot position of the tunnel is identified as the most vulnerable to damage.
基金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 the National Natural Science Foundation of China(11961030),the Natural Science Foundation of Shaanxi Province(2022JM-034)。
摘要The predation mechanism of invertebrates(e.g.,Tortanus dextrilobatus)on plankton in aquatic population ecosystem is a significant research topic.In this paper,the interaction between invertebrates and plankton is simulated by a modified Leslie-Gower predator-prey model.Using the theory of reaction-diffusion equations,a priori estimate,existence,uniqueness and stability conditions of the positive steady state solution are established.Furthermore,numerical simulations are conducted to quantitatively analyze the dynamical behavior.The research shows that as long as the Allee effect constant satisfies the appropriate relationship and the growth rates of predator and prey are appropriately large,the predator and prey can not only coexist,but also the coexistence mode is unique and stable under low predation-rate.In addition,the numerical simulations show that the coexistence may be stable under high predation-rate.Meanwhile,with the increase of predation rate,the population density of predators will decrease.