Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data imm...Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.展开更多
The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mes...The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.展开更多
The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of...The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.展开更多
This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic...This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.展开更多
The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary....The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.展开更多
Fast and accurate transient stability analysis is crucial to power system operation.With high penetration level of wind power resources,practical dynamic security region(PDSR)with hyper plane expression has outstandin...Fast and accurate transient stability analysis is crucial to power system operation.With high penetration level of wind power resources,practical dynamic security region(PDSR)with hyper plane expression has outstanding advantages in situational awareness and series of optimization problems.The precondition of obtaining accurate PDSR boundary is to locate sufficient points around the boundary(critical points).Therefore,this paper proposes a space division and Wasserstein generative adversarial network with gra-dient penalty(WGAN-GP)based fast generation method of PDSR boundary.First,the typical differential topological characterizations of dynamic security region(DSR)provide strong theoretical foundation that the interior of DSR is hole-free and the boundaries of DSR are tight and knot-free.Then,the space division method is proposed to calculate critical operation area where the PDSR boundary is located,tremen-dously compressing the search space to locate critical points and improving the confidence level of boundary fitting result.Furthermore,the WGAN-GP model is utilized to fast obtain large number of criti-cal points based on learning the data distribution of the small training set aforementioned.Finally,the PDSR boundary with hyperplanes is fitted by the least square method.The case study is tested on the Institute of Electrical and Electronics Engineers(IEEE)39-bus system and the results verify the accuracy and efficiency of the proposed method.展开更多
The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through...The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and ga...To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.展开更多
This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inv...This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inversions(RIs)were the most frequent events(99.8%),followed by surface-based temperature inversions(SBIs,63%),elevated temperature inversions(EIs,50%),and low-level jets(LLJs,41%).(2)Higher near-surface temperature leads to a decrease in EI height,thickness,and intensity.(3)When EIs occurred,RIs were consistently observed beneath the base of the EIs.(4)The intensity of the RIs correlated negatively with surface humidity,whereas the intensity of the LLJs correlated positively with surface winds except in summer.(5)The boundary layer heights(BLHs)are best estimated using the potential temperature gradient method when EIs occur,otherwise,the Richardson number method is the best method.In contrast,the relative humidity and temperature gradient methods exhibit high dispersion or uncertainty.(6)The BLHs showed seasonal variability,peaking in winter(350–1450 m).These findings advance our understanding of ABL dynamics and BLH estimation over tropical islands.展开更多
The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary...The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.展开更多
The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is m...The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between c...THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between crystallographic orientation,grain boundary(GB)proximity,and pore characteristics(size/location).This study compares single-crystal nickel models along[100],[110],and[111]orientations with equiaxed polycrystalline models containing 0,1,and 2.5 nm pores in surface and subsurface configurations.Our results reveal that crystallographic anisotropy manifests as a 24.4%higher elastic modulus and 22.2%greater hardness in[111]-oriented single crystals compared to[100].Pore-GB synergistic effects are found to dominate the deformation behavior:2.5 nm subsurface pores reduce hardness by 25.2%through stress concentration and dislocation annihilation at GBs,whereas surface pores enable mechanical recovery via accelerated dislocation generation post-collapse.Additionally,size-dependent deformation regimes were identified,with 1 nm pores inducing negligible perturbation due to rapid atomic rearrangement,in contrast with persistent softening in 2.5 nm pores.These findings establish atomic-scale design principles for defect engineering in nickel-based aerospace components,demonstrating how crystallographic orientation,pore configuration,and GB interactions collectively govern nanoindentation behavior.展开更多
This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the pred...This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.展开更多
Cobalt-free LiNiO_(2)(LNO)is considered a promising cathode for its high energy density and costeffectiveness.However,its structural instability under deep delithiation severely limits practical application in nextgen...Cobalt-free LiNiO_(2)(LNO)is considered a promising cathode for its high energy density and costeffectiveness.However,its structural instability under deep delithiation severely limits practical application in nextgeneration batteries.Herein,we propose a high-valence Mo6+doping strategy to simultaneously improve mechanical robustness and electrochemical stability.By stabilizing intergranular interfaces,this method effectively suppresses mechanical degradation induced by lattice strain under deep delithiation.The modified cathode exhibits exceptional electrochemical performance,achieving a specific capacity of 234 mAh·g^(-1)at 0.1 C with 83.4% retention over 100 cycles at 45℃ in lithium-ion batteries(LIBs).Notably,it maintains comparable efficacy in all-solid-state batteries(ASSBs),delivering 239 mAh·g^(-1)at 0.05 C and 82.8% retention after 300 cycles.Density functional theory(DFT)calculations demonstrate a pronounced rise in oxygen vacancy formation energy,increasing from 1.42 to 3.27 eV.These findings offer valuable insights into overcoming the kinetic performance limitations of cobalt-free LNO under deep delithiation conditions.展开更多
The operational efficiency of membrane electrode assemblies in direct liquid fuel cells is critically dependent on the fuel purity in the anode compartment.To address the inherent challenge of fuel mixing problem in a...The operational efficiency of membrane electrode assemblies in direct liquid fuel cells is critically dependent on the fuel purity in the anode compartment.To address the inherent challenge of fuel mixing problem in alcohol systems,we propose a rational catalyst design strategy focusing on morphological and compositional optimization.Sodium borohydride-derived PtCuMo alloy aerogels(AA)exhibit abundant grain boundary defects,while solvothermally prepared nanowire arrays(NA)maintain excellent single-crystalline characteristics.Density functional theory calculations demonstrate that engineered grain boundaries can effectively broaden the adsorption energy window for key reaction intermediates,enabling superior adaptability to diverse catalytic pathways.By precisely controlling Cu content,we identified Pt_(3)Cu_(3)Mo_(0.5)AA as the optimal catalyst configuration,demonstrating 150% enhancement in methanol oxidation reaction activity compared to Pt_(3)Cu_(6)Mo_(0.5)NA(1.5 vs.0.6 A·mg_(Pt)^(-1))and 17% improvement in ethanol oxidation reaction performance versus Pt_(3)Cu_(1)Mo_(0.5)NA(0.82 vs.0.70 A·mg_(Pt)^(-1)).Practical application testing using gas diffusion electrodes(anode loading:0.85 mg_(Pt)·cm^(-2))achieved a mass-specific power density of 14.14 W·g_(Pt)^(-1)in 1:1 methanol/ethanol blends,representing a 3.5-fold improvement over commercial Pt/C benchmarks.This work establishes a fundamental framework for developing highperformance,broad-spectrum electrocatalysts in advanced fuel cell systems.展开更多
In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical m...In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.展开更多
基金Supported by the National Key R&D Program of China(2022YFB3303500).
文摘Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.
基金partially supported by the National Natural Science Foundation of China(Nos.92271103,12202191)。
文摘The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.
基金funded by research grant OPUS no.2021/41/B/ST8/02432 entitled Probabilistic entropy in engineering computations sponsored by The National Science Center in Polandthe Institute of Structural Analysis of Poznan University of Technology in the framework of the internal research grant 0411/SBAD/0010.
文摘The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.
基金supported by the Shanxi Scholarship Council of China(Grant No.2023-036)the Natural Science Foundation of Shanxi Province(Grant No.202303021222020).
文摘This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.
基金Supported by the Key Scientific Research Plan of Colleges and Universities in Henan Province(23B140006)。
文摘The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.
基金funded in part by the National Key Research and Development Program of China(2020YFB0905900)in part by Science and Technology Project of State Grid Corporation of China(SGCC)The Key Technologies for Electric Internet of Things(SGTJDK00DWJS2100223).
文摘Fast and accurate transient stability analysis is crucial to power system operation.With high penetration level of wind power resources,practical dynamic security region(PDSR)with hyper plane expression has outstanding advantages in situational awareness and series of optimization problems.The precondition of obtaining accurate PDSR boundary is to locate sufficient points around the boundary(critical points).Therefore,this paper proposes a space division and Wasserstein generative adversarial network with gra-dient penalty(WGAN-GP)based fast generation method of PDSR boundary.First,the typical differential topological characterizations of dynamic security region(DSR)provide strong theoretical foundation that the interior of DSR is hole-free and the boundaries of DSR are tight and knot-free.Then,the space division method is proposed to calculate critical operation area where the PDSR boundary is located,tremen-dously compressing the search space to locate critical points and improving the confidence level of boundary fitting result.Furthermore,the WGAN-GP model is utilized to fast obtain large number of criti-cal points based on learning the data distribution of the small training set aforementioned.Finally,the PDSR boundary with hyperplanes is fitted by the least square method.The case study is tested on the Institute of Electrical and Electronics Engineers(IEEE)39-bus system and the results verify the accuracy and efficiency of the proposed method.
基金supported by the National Natural Science Foundation of China(Grant 12472113).
文摘The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
基金Supported by the China National Natural Science Foundation Project(52274048)Beijing Natural Science Foundation Project(3222037)。
文摘To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.
基金National Key Research and Development Program of China(2023YFC3008002)National Natural Science Foundation of China(U21A6001,42075059)Key Laboratory of Guangdong Province(2020B1212060025)。
文摘This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inversions(RIs)were the most frequent events(99.8%),followed by surface-based temperature inversions(SBIs,63%),elevated temperature inversions(EIs,50%),and low-level jets(LLJs,41%).(2)Higher near-surface temperature leads to a decrease in EI height,thickness,and intensity.(3)When EIs occurred,RIs were consistently observed beneath the base of the EIs.(4)The intensity of the RIs correlated negatively with surface humidity,whereas the intensity of the LLJs correlated positively with surface winds except in summer.(5)The boundary layer heights(BLHs)are best estimated using the potential temperature gradient method when EIs occur,otherwise,the Richardson number method is the best method.In contrast,the relative humidity and temperature gradient methods exhibit high dispersion or uncertainty.(6)The BLHs showed seasonal variability,peaking in winter(350–1450 m).These findings advance our understanding of ABL dynamics and BLH estimation over tropical islands.
基金supported by the National Natural Science Foundation of China(Grant No.52071094).
文摘The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.
基金supported by the Science and Technology Research Project of Henan Province(242102241055)the Industry-University-Research Collaborative Innovation Base on Automobile Lightweight of“Science and Technology Innovation in Central Plains”(2024KCZY315)the Opening Fund of State Key Laboratory of Structural Analysis,Optimization and CAE Software for Industrial Equipment(GZ2024A03-ZZU).
文摘The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金The National Natural Science Foundation of China(Grant No.12462006)Beijing Institute of Structure and Environment Engineering Joint Innovation Fund(No.BQJJ202414).
文摘THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between crystallographic orientation,grain boundary(GB)proximity,and pore characteristics(size/location).This study compares single-crystal nickel models along[100],[110],and[111]orientations with equiaxed polycrystalline models containing 0,1,and 2.5 nm pores in surface and subsurface configurations.Our results reveal that crystallographic anisotropy manifests as a 24.4%higher elastic modulus and 22.2%greater hardness in[111]-oriented single crystals compared to[100].Pore-GB synergistic effects are found to dominate the deformation behavior:2.5 nm subsurface pores reduce hardness by 25.2%through stress concentration and dislocation annihilation at GBs,whereas surface pores enable mechanical recovery via accelerated dislocation generation post-collapse.Additionally,size-dependent deformation regimes were identified,with 1 nm pores inducing negligible perturbation due to rapid atomic rearrangement,in contrast with persistent softening in 2.5 nm pores.These findings establish atomic-scale design principles for defect engineering in nickel-based aerospace components,demonstrating how crystallographic orientation,pore configuration,and GB interactions collectively govern nanoindentation behavior.
基金supported by the National Key R&D Program of China[grant number 2023YFC3008004]。
文摘This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.
基金the Guangxi New Energy Vehicle Laboratory Special Project(No.GuikeAA23062079)the Special Project for Central Guidance of Local Science and Technology Development in Sichuan Province(No.2024ZYD0302).
文摘Cobalt-free LiNiO_(2)(LNO)is considered a promising cathode for its high energy density and costeffectiveness.However,its structural instability under deep delithiation severely limits practical application in nextgeneration batteries.Herein,we propose a high-valence Mo6+doping strategy to simultaneously improve mechanical robustness and electrochemical stability.By stabilizing intergranular interfaces,this method effectively suppresses mechanical degradation induced by lattice strain under deep delithiation.The modified cathode exhibits exceptional electrochemical performance,achieving a specific capacity of 234 mAh·g^(-1)at 0.1 C with 83.4% retention over 100 cycles at 45℃ in lithium-ion batteries(LIBs).Notably,it maintains comparable efficacy in all-solid-state batteries(ASSBs),delivering 239 mAh·g^(-1)at 0.05 C and 82.8% retention after 300 cycles.Density functional theory(DFT)calculations demonstrate a pronounced rise in oxygen vacancy formation energy,increasing from 1.42 to 3.27 eV.These findings offer valuable insights into overcoming the kinetic performance limitations of cobalt-free LNO under deep delithiation conditions.
基金financially supported by the National Natural Science Foundation of China(No.52073214)Guangxi Natural Science Fund for Distinguished Young Scholars(No.2024GXNSFFA010008).
文摘The operational efficiency of membrane electrode assemblies in direct liquid fuel cells is critically dependent on the fuel purity in the anode compartment.To address the inherent challenge of fuel mixing problem in alcohol systems,we propose a rational catalyst design strategy focusing on morphological and compositional optimization.Sodium borohydride-derived PtCuMo alloy aerogels(AA)exhibit abundant grain boundary defects,while solvothermally prepared nanowire arrays(NA)maintain excellent single-crystalline characteristics.Density functional theory calculations demonstrate that engineered grain boundaries can effectively broaden the adsorption energy window for key reaction intermediates,enabling superior adaptability to diverse catalytic pathways.By precisely controlling Cu content,we identified Pt_(3)Cu_(3)Mo_(0.5)AA as the optimal catalyst configuration,demonstrating 150% enhancement in methanol oxidation reaction activity compared to Pt_(3)Cu_(6)Mo_(0.5)NA(1.5 vs.0.6 A·mg_(Pt)^(-1))and 17% improvement in ethanol oxidation reaction performance versus Pt_(3)Cu_(1)Mo_(0.5)NA(0.82 vs.0.70 A·mg_(Pt)^(-1)).Practical application testing using gas diffusion electrodes(anode loading:0.85 mg_(Pt)·cm^(-2))achieved a mass-specific power density of 14.14 W·g_(Pt)^(-1)in 1:1 methanol/ethanol blends,representing a 3.5-fold improvement over commercial Pt/C benchmarks.This work establishes a fundamental framework for developing highperformance,broad-spectrum electrocatalysts in advanced fuel cell systems.
基金supported by Anhui Provincial Natural Science Foundation(2408085QA030)Natural Science Research Project of Anhui Educational Committee,China(2022AH050825)+3 种基金Medical Special Cultivation Project of Anhui University of Science and Technology(YZ2023H2C008)the Excellent Research and Innovation Team of Anhui Province,China(2022AH010052)the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology,China(2021yjrc51)Collaborative Innovation Program of Hefei Science Center,CAS,China(2019HSC-CIP006).
文摘In this paper,a novel method for investigating the particle-crushing behavior of breeding particles in a fusion blanket is proposed.The fractal theory and Weibull distribution are combined to establish a theoretical model,and its validity was verified using a simple impact test.A crushable discrete element method(DEM)framework is built based on the previously established theoretical model.The tensile strength,which considers the fractal theory,size effect,and Weibull variation,was assigned to each generated particle.The assigned strength is then used for crush detection by comparing it with its maximum tensile stress.Mass conservation is ensured by inserting a series of sub-particles whose total mass was equal to the quality loss.Based on the crushable DEM framework,a numerical simulation of the crushing behavior of a pebble bed with hollow cylindrical geometry under a uniaxial compression test was performed.The results of this investigation showed that the particle withstands the external load by contact and sliding at the beginning of the compression process,and the results confirmed that crushing can be considered an important method of resisting the increasing external load.A relatively regular particle arrangement aids in resisting the load and reduces the occurrence of particle crushing.However,a limit exists to the promotion of resistance.When the strain increases beyond this limit,the distribution of the crushing position tends to be isotropic over the entire pebble bed.The theoretical model and crushable DEM framework provide a new method for exploring the pebble bed in a fusion reactor,considering particle crushing.
