To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRE...The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.展开更多
Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuse...Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.展开更多
Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie conditio...Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.展开更多
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 glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction...The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.展开更多
The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under resho...The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under reshock impact is investigated using a two-component discrete Boltzmann method.The hydrodynamic and thermodynamic characteristics of the fluid system are comprehensively analyzed from the perspectives of the density gradient,vorticity,kinetic energy,mixing degree,mixing width and non-equilibrium intensity.Simulation results indicate that for longer relaxation time,the diffusion and dissipation are enhanced,the physical gradients decrease,and the growth of the interface is suppressed.Furthermore,the non-equilibrium manifestations show complex patterns,driven by the competitive physical mechanisms of the diffusion,dissipation,shock wave,rarefaction wave,transverse wave and fluid instabilities.These findings provide valuable insight into the fundamental mechanism of compressible fluid flows.展开更多
Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres wi...Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
Gas storage in abandoned mines is one way to reuse waste space resources.The surrounding rock of gas storage reservoirs in underground roadways undergoes damage and deformation under the cyclic loading of gas charging...Gas storage in abandoned mines is one way to reuse waste space resources.The surrounding rock of gas storage reservoirs in underground roadways undergoes damage and deformation under the cyclic loading of gas charging and discharging,which can pose a risk to the safety of the reservoirs.This study establishes a true triaxial numerical model of rock mass with the discrete element method(DEM)and explores the crack evolution of surrounding rock of underground gas storage during cyclic loading and unloading.Also,a damage evolution model in numerical analysis considering residual deformation is developed to explain the experimental results.As was revealed,cyclic loading and unloading resulted in fatigue damage in the specimen and caused strength deterioration of the specimen.During the loading process,the uniformly distributed force chains of the rock mass redistributed,evolving gradually to mostly transverse force chains.This contributed to the appearance of blank areas in the force chains when through cracks appear.The ratio of tensile cracks to shear cracks gradually decreases and finally stabilizes at 7:1.The damage evolution model considering residual strain can be mutually verified with the numerical simulation results.Based on the DEM model,it was found that there was a certain threshold of confining pressure.When the confining pressure exceeded 30 MPa,the deformation to ductility of sandstone samples began to accelerate,with a greater residual strength.This study provides a theoretical basis for analyzing the long-term mechanical behavior of surrounding rock of gas storage in abandoned mines.展开更多
Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is sti...Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.展开更多
To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hyd...Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.展开更多
Based on mathematical orthogonality and mechanical equilibrium,a deformation energy decomposition method for classical isotropic square and cube elements is proposed by considering the physical parameters of materials...Based on mathematical orthogonality and mechanical equilibrium,a deformation energy decomposition method for classical isotropic square and cube elements is proposed by considering the physical parameters of materials.By aid of this method,the comprehensive deformation energy of planar discrete elastomers can be decomposed into five basic deformation energies,and the comprehensive deformation energy of spatial discrete elastomers can be decomposed into eighteen basic deformation energies.The quantification and visualization of structural deformation performance can be realized.According to the magnitude of different deformation energy in the same element,the decomposition diagram is drawn,which can visually display the area dominated by each basic deformation energy.The cloud diagram is drawn based on the distribution of specific deformation energy in different elements,which can be used to analyze the gradient change of deformation energy in the structure.Finally,the deformation properties of cantilever beam and four-sided consolidation plate are analyzed by deformation energy decomposition method.The correctness and superiority of this method are verified by comparing with the results of strain energy decomposition.展开更多
The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is ess...The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is essential for enhancing the performance and durability of SOFCs.Accurate prediction of mechanical properties of porous electrode can be achieved by microscale finite element modeling based on three-dimensional(3D)microstructures,which requires expensive 3D tomography techniques and massive computational resources.In this study,we proposed a cost-effective alternative approach to access the mechanical properties of porous electrodes,with the elastic properties of La_(0.6)Sr_(0.4)Co_(0.2)Fe_(0.8)O_(3-δc)athode serving as a case study.Firstly,a stochastic modeling was used to reconstruct 3D microstructures from two-dimensional(2D)cross-sections as an alternative to expensive tomography.Then,the discrete element method(DEM)was used to predict the elastic properties of porous ceramics based on the discretized 3D microstructures reconstructed by stochastic modeling.Based on 2D microstructure and the elastic properties calculated by the DEM modeling of the 3D reconstructed porous microstructures,a convolutional neural network(CNN)based deep learning model was built to predict the elastic properties rapidly from 2D microstructures.The proposed combined framework can be implemented with limited computational resources and provide a basis for rapid prediction of mechanical properties and parameter estimation for multiscale modeling of SOFCs.展开更多
Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sen...Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sensing technology was used to perform strain monitoring,but there is an urgent need to determine how to transform strain into displacement.In this study,optical frequency domain reflectometry was applied in laboratory tests.Aiming at the shear process and center settlement process of shield tunnel segments,two kinds of quantitative calculation methods were put forward to carry out a quantitative analysis.Meanwhile,the laboratory test process was simulated numerically utilizing the discrete element numerical analysis method.Optical fiber,an atypical geotechnical material,was innovatively applied for discrete element modeling and numerical simulation.The results show that the measured displacement of the dial gauge,the calculated results of the numerical model,and the displacement quantitatively calculated from the optical fiber data agree with each other in general.The latter two methods can potentially be utilized in engineering application of deformation monitoring at shield tunnel joints,but need to be further calibrated and adjusted in detail.展开更多
Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between ...Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between coarser and finer particles,this research differentiates the stiffness distribution of gap-graded soils from the perspective of contact and particle types.Results indicate that particle property disparity significantly influence the small-strain stiffness characteristics,consequently altering the overall stiffness distribution in gap-graded soil specimens.Specifically,with the equivalent coarser particle property,an increase in particle Young's modulus of finer particles results in an augmentation of small-strain stiffness values,alongside an increased stiffness distribution contribution from finer particles.Nevertheless,this study reveals that even with a higher particle Young's modulus of finer particles,the proportion of small-strain stiffness transferred by finer particles remains consistently lower than their volume fraction.Furthermore,the proportion of stiffness transferred by finer particles may fall below their contribution to stress transmission.This investigation accentuates the subtle yet significant effects of particle property variations on small strain stiffness and its subsequent distribution,providing a foundation for advancing the significance of particle property disparities in evaluating soil responses.展开更多
It is a major challenge for the airframe-inlet design of modern combat aircrafts,as the flow and electromagnetic wave propagation in the inlet of stealth aircraft are very complex.In this study,an aerodynamic/stealth ...It is a major challenge for the airframe-inlet design of modern combat aircrafts,as the flow and electromagnetic wave propagation in the inlet of stealth aircraft are very complex.In this study,an aerodynamic/stealth optimization design method for an S-duct inlet is proposed.The upwind scheme is introduced to the aerodynamic adjoint equation to resolve the shock wave and flow separation.The multilevel fast multipole algorithm(MLFMA)is utilized for the stealth adjoint equation.A dorsal S-duct inlet of flying wing layout is optimized to improve the aerodynamic and stealth characteristics.Both the aerodynamic and stealth characteristics of the inlet are effectively improved.Finally,the optimization results are analyzed,and it shows that the main contradiction between aerodynamic characteristics and stealth characteristics is the centerline and crosssectional area.The S-duct is smoothed,and the cross-sectional area is increased to improve the aerodynamic characteristics,while it is completely opposite for the stealth design.The radar cross section(RCS)is reduced by phase cancelation for low frequency conditions.The method is suitable for the aerodynamic/stealth design of the aircraft airframe-inlet system.展开更多
To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 ...To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.展开更多
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金funded by the National Natural Science Foundation of China(No.51978048).
文摘The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.
基金the project of the State Key Laboratory of Petroleum Resources and Engineering(No.PRE/open-2307)the CNOOC Research Institute(No.2020PFS-03).
文摘Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.
基金a U.S. Federal Railroad Administration (FRA)BAA project,titled “Mitigation of Differential Movement at Railway Transitions for High-Speed Passenger Rail and Joint Passenger/Freight Corridors”the financial support provided by the China Scholarship Council (CSC),which funded Zhongyi Liu’s and Wenjing Li’s time and research efforts for this study
文摘Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.
基金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 National Natural Science Foundation of China(Grant No.52304003)the Natural Science Foundation of Sichuan Province(Grant No.2024NSFSC0961)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230090).
文摘The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.
基金supported by the National Natural Science Foundation of China(Grant No.U2242214)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2024A1515010927)+3 种基金the Humanities and Social Science Foundation of the Ministry of Education in China(Grant No.24YJCZH163)the Fujian Provincial Units Special Funds for Education and Research(Grant No.2022639)Fundamental Research Funds for the Central Universities,Sun Yatsen University(Grant No.24qnpy044)partly supported by the Open Research Fund of Key Laboratory of Analytical Mathematics and Applications(Fujian Normal University),Ministry of Education,P.R.China(Grant No.JAM2405)。
文摘The Richtmyer–Meshkov(RM)instability plays an important role in various natural and engineering fields such as inertial confinement fusion.In this study,the effect of relaxation time on the RM instability under reshock impact is investigated using a two-component discrete Boltzmann method.The hydrodynamic and thermodynamic characteristics of the fluid system are comprehensively analyzed from the perspectives of the density gradient,vorticity,kinetic energy,mixing degree,mixing width and non-equilibrium intensity.Simulation results indicate that for longer relaxation time,the diffusion and dissipation are enhanced,the physical gradients decrease,and the growth of the interface is suppressed.Furthermore,the non-equilibrium manifestations show complex patterns,driven by the competitive physical mechanisms of the diffusion,dissipation,shock wave,rarefaction wave,transverse wave and fluid instabilities.These findings provide valuable insight into the fundamental mechanism of compressible fluid flows.
基金the financial support provided by Tianfu Yongxing Laboratory Organized Research Project Funding(No.2023CXXM01)the ARC linkage program(No.LP200100420).
文摘Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金National Natural Science Foundation of China,Grant/Award Numbers:U22A20598,52104107National Key Research and Development Program of China,Grant/Award Numbers:2023YFC2907300,2019YFE0118500,2019YFC1904304Natural Science Foundation of Jiangsu Province,Grant/Award Number:BK20200634。
文摘Gas storage in abandoned mines is one way to reuse waste space resources.The surrounding rock of gas storage reservoirs in underground roadways undergoes damage and deformation under the cyclic loading of gas charging and discharging,which can pose a risk to the safety of the reservoirs.This study establishes a true triaxial numerical model of rock mass with the discrete element method(DEM)and explores the crack evolution of surrounding rock of underground gas storage during cyclic loading and unloading.Also,a damage evolution model in numerical analysis considering residual deformation is developed to explain the experimental results.As was revealed,cyclic loading and unloading resulted in fatigue damage in the specimen and caused strength deterioration of the specimen.During the loading process,the uniformly distributed force chains of the rock mass redistributed,evolving gradually to mostly transverse force chains.This contributed to the appearance of blank areas in the force chains when through cracks appear.The ratio of tensile cracks to shear cracks gradually decreases and finally stabilizes at 7:1.The damage evolution model considering residual strain can be mutually verified with the numerical simulation results.Based on the DEM model,it was found that there was a certain threshold of confining pressure.When the confining pressure exceeded 30 MPa,the deformation to ductility of sandstone samples began to accelerate,with a greater residual strength.This study provides a theoretical basis for analyzing the long-term mechanical behavior of surrounding rock of gas storage in abandoned mines.
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
基金Australian Research Council Linkage Program(LP200301404)for sponsoring this researchthe financial support provided by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology,SKLGP2021K002)National Natural Science Foundation of China(52374101,32111530138).
文摘Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.
基金supported by the National Natural Science Foundation of China(Grant Nos.52208322 and 51878621)the Postdoctoral Science Foundation of China(Grant No.2022M712905)the Key Research Projects of Henan Higher Education Institutions(Grant No.22A560005).
文摘Based on mathematical orthogonality and mechanical equilibrium,a deformation energy decomposition method for classical isotropic square and cube elements is proposed by considering the physical parameters of materials.By aid of this method,the comprehensive deformation energy of planar discrete elastomers can be decomposed into five basic deformation energies,and the comprehensive deformation energy of spatial discrete elastomers can be decomposed into eighteen basic deformation energies.The quantification and visualization of structural deformation performance can be realized.According to the magnitude of different deformation energy in the same element,the decomposition diagram is drawn,which can visually display the area dominated by each basic deformation energy.The cloud diagram is drawn based on the distribution of specific deformation energy in different elements,which can be used to analyze the gradient change of deformation energy in the structure.Finally,the deformation properties of cantilever beam and four-sided consolidation plate are analyzed by deformation energy decomposition method.The correctness and superiority of this method are verified by comparing with the results of strain energy decomposition.
基金supported by the National Natural Science Foundation of China(Grant Nos.12172104 and 11932005)the Talent Recruitment Project of Guangdong(2021QN02L892)+3 种基金the Stable Supporting Fund of Shenzhen(GXWD20231130153335002)the Shccig-Qinling Program(SMYJY202300140C)the program of Innovation Team in Universities and Colleges in Guangdong(2021KCXTD006)Development and Reform Commission of Shenzhen(XMHT20220103004).
文摘The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is essential for enhancing the performance and durability of SOFCs.Accurate prediction of mechanical properties of porous electrode can be achieved by microscale finite element modeling based on three-dimensional(3D)microstructures,which requires expensive 3D tomography techniques and massive computational resources.In this study,we proposed a cost-effective alternative approach to access the mechanical properties of porous electrodes,with the elastic properties of La_(0.6)Sr_(0.4)Co_(0.2)Fe_(0.8)O_(3-δc)athode serving as a case study.Firstly,a stochastic modeling was used to reconstruct 3D microstructures from two-dimensional(2D)cross-sections as an alternative to expensive tomography.Then,the discrete element method(DEM)was used to predict the elastic properties of porous ceramics based on the discretized 3D microstructures reconstructed by stochastic modeling.Based on 2D microstructure and the elastic properties calculated by the DEM modeling of the 3D reconstructed porous microstructures,a convolutional neural network(CNN)based deep learning model was built to predict the elastic properties rapidly from 2D microstructures.The proposed combined framework can be implemented with limited computational resources and provide a basis for rapid prediction of mechanical properties and parameter estimation for multiscale modeling of SOFCs.
基金National Natural Science Foundation of China,Grant/Award Numbers:41977218,42222707State Key Laboratory for GeoMechanics and Deep Underground Engineering,Grant/Award Number:SKLGDUEK2117。
文摘Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sensing technology was used to perform strain monitoring,but there is an urgent need to determine how to transform strain into displacement.In this study,optical frequency domain reflectometry was applied in laboratory tests.Aiming at the shear process and center settlement process of shield tunnel segments,two kinds of quantitative calculation methods were put forward to carry out a quantitative analysis.Meanwhile,the laboratory test process was simulated numerically utilizing the discrete element numerical analysis method.Optical fiber,an atypical geotechnical material,was innovatively applied for discrete element modeling and numerical simulation.The results show that the measured displacement of the dial gauge,the calculated results of the numerical model,and the displacement quantitatively calculated from the optical fiber data agree with each other in general.The latter two methods can potentially be utilized in engineering application of deformation monitoring at shield tunnel joints,but need to be further calibrated and adjusted in detail.
基金Financial supports from the PolyU Distinguished Postdoctoral Fellowship Scheme are highly appreciatedsupported by the National Natural Science Foundation of China (Grant No.52201008)the Fundamental Research Funds for the Central Universities,the State Key Laboratory of Particle Detection and Electronics (Grant No.SKLPDE-KF-202311).
文摘Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between coarser and finer particles,this research differentiates the stiffness distribution of gap-graded soils from the perspective of contact and particle types.Results indicate that particle property disparity significantly influence the small-strain stiffness characteristics,consequently altering the overall stiffness distribution in gap-graded soil specimens.Specifically,with the equivalent coarser particle property,an increase in particle Young's modulus of finer particles results in an augmentation of small-strain stiffness values,alongside an increased stiffness distribution contribution from finer particles.Nevertheless,this study reveals that even with a higher particle Young's modulus of finer particles,the proportion of small-strain stiffness transferred by finer particles remains consistently lower than their volume fraction.Furthermore,the proportion of stiffness transferred by finer particles may fall below their contribution to stress transmission.This investigation accentuates the subtle yet significant effects of particle property variations on small strain stiffness and its subsequent distribution,providing a foundation for advancing the significance of particle property disparities in evaluating soil responses.
文摘It is a major challenge for the airframe-inlet design of modern combat aircrafts,as the flow and electromagnetic wave propagation in the inlet of stealth aircraft are very complex.In this study,an aerodynamic/stealth optimization design method for an S-duct inlet is proposed.The upwind scheme is introduced to the aerodynamic adjoint equation to resolve the shock wave and flow separation.The multilevel fast multipole algorithm(MLFMA)is utilized for the stealth adjoint equation.A dorsal S-duct inlet of flying wing layout is optimized to improve the aerodynamic and stealth characteristics.Both the aerodynamic and stealth characteristics of the inlet are effectively improved.Finally,the optimization results are analyzed,and it shows that the main contradiction between aerodynamic characteristics and stealth characteristics is the centerline and crosssectional area.The S-duct is smoothed,and the cross-sectional area is increased to improve the aerodynamic characteristics,while it is completely opposite for the stealth design.The radar cross section(RCS)is reduced by phase cancelation for low frequency conditions.The method is suitable for the aerodynamic/stealth design of the aircraft airframe-inlet system.
基金supported by the National Key R&D Program of China (Grant No.2022YFC3003401)the National Natural Science Foundation of China (Grant Nos.42041006 and 42377137).
文摘To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.
