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.展开更多
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 conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete bounda...In conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete boundary unit cells(i.e.,boundary defects)even though the bulk polarization is zero,which challenges the conventional understanding of HOTIs.Here,based on a Kekul´e-distorted honeycomb lattice with incomplete unit cells,we reveal that incomplete unit cells exhibit fractional charges through the analysis of Wannier centers by developing a compensation method and creating the concept of Wannier center domain(WCD)which is the smallest region that one Wannier center occupies.This method compensates for the missing parts of these boundary incomplete unit cells with additional WCDs to make them complete.The compensated WCDs automatically carry the corresponding charge,and this charge together with that of the incomplete unit cell constitutes the total charge of the complete unit cell after compensation.We conclude that the emergence of corner states is attributed to the filling anomaly,which is a fundamental mechanism.Our results refresh the understanding of HOTIs,especially those with structural discontinuities,and provide a novel design for topological states which have application value in producing optical functional devices.展开更多
We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary,subject to the kinematic and vorticity boundary conditio...We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary,subject to the kinematic and vorticity boundary conditions on the non-flat boundary.We observe that,under the nonhomogeneous boundary conditions,the pressure p can be still recovered by solving the Neumann problem for the Poisson equation.Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions.Based on this,we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition(without the incompressibility condition),which establishes a velocity mapping.Then we develop apriori estimates for the velocity mapping,especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions,which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem.Finally,we establish that,when the viscosity coefficient tends zero,the strong solutions of the initial-boundary value problem in R^n(n≥3)with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.展开更多
Due to the coexistence of compressibility,viscosity,and threedimensional effects,laminar flow is difficult to maintain for high-speed boundary layer on complex geometries.The unstable disturbance waves in the boundary...Due to the coexistence of compressibility,viscosity,and threedimensional effects,laminar flow is difficult to maintain for high-speed boundary layer on complex geometries.The unstable disturbance waves in the boundary layer are excited and rapidly increase during the receptivity process,so sufficiently large Reynolds stress causes the basic flow velocity profile to change,and the formation of turbulence is inevitable.展开更多
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 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 this paper,we consider incompressible Navier-Stokes/Cahn-Hilliard system with the generalized Navier boundary condition and the dynamic boundary condition in a channel,which can describe the interaction between a b...In this paper,we consider incompressible Navier-Stokes/Cahn-Hilliard system with the generalized Navier boundary condition and the dynamic boundary condition in a channel,which can describe the interaction between a binary material and the walls of the physical domain.We prove the global-in-time existence and uniqueness of strong solutions to this initial boundary value problem in a 2D channel domain.展开更多
Vorticity is locally generated on a boundary at a rate measured by the boundary vorticity flux(BVF),which can be further decomposed into the sum of the orbital rotation and the generalized spin(specifically,the sum of...Vorticity is locally generated on a boundary at a rate measured by the boundary vorticity flux(BVF),which can be further decomposed into the sum of the orbital rotation and the generalized spin(specifically,the sum of shear and streaming vorticity)under the field description.For an incompressible viscous flow interacting with a stationary wall,the full expressions of the boundary fluxes of the orbital rotation and the spin are derived,for the first time,to elucidate their boundary creation mechanisms.Then,these new findings are successfully extended to the study of the boundary enstrophy dynamics,as well as the Lyman vorticity dynamics as an alternative interpretation to the boundary vorticity dynamics.Interestingly,it is found that the boundary coupling of the longitudinal and transverse processes is only embodied in the boundary spin flux,which is definitely not responsible for the generation of the boundary orbital-rotation flux.In addition,the boundary fluxes of enstrophy are directly associated with the boundary source of the second principal invariant of the velocity gradient tensor(VGT)and the two quadratic forms representing the spin-geometry interaction.The present exposition provides a new perspective and an additional dimension for understanding the vorticity dynamics on boundaries,which could be valuable in clarifying the formation mechanisms of near-wall coherent structures and flow noise at the fundamental level.展开更多
Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction len...Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction lengths have not been sufficiently investigated. First, this study presents a theoretical scaling analysis and validates it through wind tunnel experiments. It conducts detailed control volume analysis of mass conservation, considering the differences between inviscid and viscous cases. Then, three models for analysing interaction length under gradual expansion waves are derived. Related experiments using schlieren photography are conducted to validate the models in a Mach 2.73 flow. The interaction scales are captured at various relative distances between the shock impingement location and the expansion regions with wedge angles ranging from 12° to 15° and expansion angles of 9°, 12°, and 15°. Three trend lines are plotted based on different expansion angles to depict the relationship between normalised interaction length and normalised interaction strength metric. In addition, the relationship between the coefficients of the trend line and the expansion angles is introduced to predict the interaction length influenced by gradual expansion waves. Finally, the estimation of normalised interaction length is derived for various coefficients within a unified form.展开更多
The Beijing 325 m meteorological tower stands as a pivotal research platform for exploring atmospheric boundary layer physics and atmospheric chemistry.With a legacy spanning 45 years,the tower has played a crucial ro...The Beijing 325 m meteorological tower stands as a pivotal research platform for exploring atmospheric boundary layer physics and atmospheric chemistry.With a legacy spanning 45 years,the tower has played a crucial role in unraveling the complexities of urban air pollution,atmospheric processes,and climate change in Beijing,China.This review paper provides a comprehensive overview of the measurements on the tower over the past two decades.Through long-term comprehensive observations,researchers have elucidated the intricate relationships between anthropogenic emissions,meteorological dynamics,and atmospheric composition,shedding light on the drivers of air pollution and its impacts on public health.The vertical measurements on the tower also enable detailed investigations into boundary layer dynamics,turbulent mixing,and pollutant dispersion,providing invaluable data for validating chemical transport models.Key findings from the tower’s research include the identification of positive feedback mechanisms between aerosols and the boundary layer,the characterization of pollutant sources and transport pathways,the determination of fluxes of gaseous and particulate species,and the assessment of the effectiveness of pollution control measures.Additionally,isotopic measurements have provided new insights into the sources and formation processes of particulate matter and reactive nitrogen species.Finally,the paper outlines future directions for tower-based research,emphasizing the need for long-term comprehensive measurements,the development of innovative tower platforms,and integration of emerging technologies.展开更多
In order to better describe the phenomenon of biological invasion,this paper introduces a free boundary model of biological invasion.Firstly,the right free boundary is added to the equation with logistic terms.Secondl...In order to better describe the phenomenon of biological invasion,this paper introduces a free boundary model of biological invasion.Firstly,the right free boundary is added to the equation with logistic terms.Secondly,the existence and uniqueness of local solutions are proved by the Sobolev embedding theorem and the comparison principle.Finally,according to the relevant research data and contents of red fire ants,the diffusion area and nest number of red fire ants were simulated without external disturbance.This paper mainly simulates the early diffusion process of red fire ants.In the early diffusion stage,red fire ants grow slowly and then spread over a large area after reaching a certain number.展开更多
The combination of dual-main-phase(DMP)(Nd,Ce)-Fe-B magnets and grain boundary diffusion process(GBDP)is currently a research topic for obtaining high-cost performance materials in rare earth permanent magnet fields.T...The combination of dual-main-phase(DMP)(Nd,Ce)-Fe-B magnets and grain boundary diffusion process(GBDP)is currently a research topic for obtaining high-cost performance materials in rare earth permanent magnet fields.The novel structural features of GBDP(Nd,Ce)-Fe-B magnets give a version of different domain reversal processes from those of non-diffused magnets.In this work,the in-situ magnetic domain evolution of the DMP magnets was observed at elevated temperatures,and the temperature demagnetization and coercivity mechanism of the GBDP dual-main-phase(Nd,Ce)-Fe-B magnets are discussed.The results show that the shell composition of different types of grains in DMP magnets is similar,while the magnetic microstructure results indicate the Ce-rich grains tend to demagnetize first.Dy-rich shell with a high anisotropic field caused by GBDP leads to an increase in the nucleation field,which enhances the coercivity.It is found that much more grains exhibit single domain characteristics in the remanent state for GBDP dual-main-phase(Nd,Ce)-Fe-B magnets.In addition,the grains that undergo demagnetization first are Ce-rich or Nd-rich grains,which is different from that of non-diffused magnets.These results were not found in previous studies but can be intuitively characterized from the perspective of magnetic domains in this work,providing a new perspective and understanding of the performance improvement of magnetic materials.展开更多
Grain boundary diffusion technology is pivotal in the preparation of high-performance NdFeB magnets.This study investigates the factors that affect the efficiency of grain boundary diffusion,starting from the properti...Grain boundary diffusion technology is pivotal in the preparation of high-performance NdFeB magnets.This study investigates the factors that affect the efficiency of grain boundary diffusion,starting from the properties of the diffusion matrix.Through the adjustment of the sintering process,we effectively prepared magnets with varied densities that serve as the matrix for grain boundary diffusion with TbH,diffusion.The mobility characteristics of the Nd-rich phase during the densification stage are leveraged to ensure a more extensive distribution of heavy rare earth elements within the magnets.According to the experimental results,the increase in coercivity of low-density magnets after diffusion is significantly greater than that of relatively high-density magnets.The coercivity values measured are 805.32 kA/m for low-density magnets and 470.3 kA/m for high-density magnets.Additionally,grain boundary diffusion notably enhances the density of initial low-density magnets,addressing the issue of low density during the sintering stage.Before the diffusion treatment,the Nd-rich phases primarily concentrate at the triangular grain boundaries,resulting in an increased number of cavity defects in the magnets.These cavity defects contain atoms in a higher energy state,making them more prone to transition.Consequently,the diffusion activation energy at the void defects is lower than the intracrystalline diffusion activation energy,accelerating atom diffusion.The presence of larger cavities also provides more space for atom migration,thereby promoting the diffusion process.After the diffusion treatment,the proportion of bulk Nd-rich phases significantly decreases,and they infiltrate between the grains to fill the cavity defects,forming continuous fine grain boundaries.Based on these observations,the study aims to explore how to utilize this information to develop an efficient technique for grain boundary diffusion.展开更多
Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymme...Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymmetric tilt boundaries are hard to be accounted for based on traditional theoretical models,and the corresponding solute segregation is complex.Herein,atomic structures of a specific asymmetric boundary on{1012}TBs were reveled using aberration-corrected high-angle annular dark-field scanning transmission electron microscopy(HAADF-STEM),molecular dynamics(MD)and density functional theory(DFT)simulations.Reaction between<a60>M dislocations and the{1012}TB can generate a~61°/25°asymmetric tilt boundary.The segregation of Gd and Zn atoms is closely related to the aggregateddislocations and the interfacial interstices of the asymmetric tilt boundary,which is energetically favorable in reducing the total system energy.展开更多
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.展开更多
To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced v...To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.展开更多
The microstructure of(Nd,Ce)-Fe-B sintered magnets with different diffusion depths was characterized by a magnetic force microscope,and the relationship between the magnetic properties and the local structure of grain...The microstructure of(Nd,Ce)-Fe-B sintered magnets with different diffusion depths was characterized by a magnetic force microscope,and the relationship between the magnetic properties and the local structure of grain boundary diffused magnets is discussed.The domains perpendicular to the c-axis(easy magnetization direction)show a typical maze-like pattern,while those parallel to the c-axis show the characte ristics of plate domains.The significant gradient change is shown in the concentration of Dy with the direction of diffusion from the surface to the interior.Dy diffuses along grain boundaries and(Dy,Nd)_(2)Fe_(14)B layer with a high anisotropy field formed around the grains.Through in-situ electron probe micro-analysis/magnetic force microscopy(EPMA/MFM),it is found that the average domain width decreases,and the proportion of single domain grains increases as diffusion depth increases.This is caused by both the change of concentration and distribution of Dy.The grain boundary diffusion process changes the microstructure and microchemistry inside the magnet,and these local magnetism differences can be reflected by the configuration of the magnetic domain structure.展开更多
The grain boundary diffusion process(GBDP)has proven to be an effective method for enhancing the coercivity of sintered Nd-Fe-B magnets.However,the limited diffusion depth and thicker shell struc-ture have impeded the...The grain boundary diffusion process(GBDP)has proven to be an effective method for enhancing the coercivity of sintered Nd-Fe-B magnets.However,the limited diffusion depth and thicker shell struc-ture have impeded the further development of magnetic properties.Currently,the primary debates re-garding the mechanism of GBDP with Tb revolve around the dissolution-solidification mechanism and the atomic substitution mechanism.To clarify this mechanism,the microstructure evolution of sintered Nd-Fe-B magnets during the heating process of GBDP has been systematically studied by quenching at different tem peratures.In this study,it was found that the formation of TbFe_(2) phase is related to the dis-solution of _(2)Fe_(14)B grains during GBDP with Tb.The theory of mixing heat and phase separation further confirms that the Nd_(2)Fe_(14)B phase dissolves to form a mixed phase of Nd and TbFe_(2),which then solidifies into the(Nd,Tb)_(2)Fe_(14)B phase.Based on the discovery of the TbFe_(2) phase,the dissolution-solidification mechanism is considered the primary mechanism for GBDP.This is supported by the elemental content of the two typical core-shell structures observed.展开更多
基金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.
基金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 by the Natural Science Basic Research Program of Shaanxi Province (Grant Nos.2024JC-JCQN-06 and2025JC-QYCX-006)the National Natural Science Foundation of China (Grant No.12474337)Chinese Academy of Sciences Project (Grant Nos.E4BA270100,E4Z127010F,E4Z6270100,and E53327020D)。
文摘In conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete boundary unit cells(i.e.,boundary defects)even though the bulk polarization is zero,which challenges the conventional understanding of HOTIs.Here,based on a Kekul´e-distorted honeycomb lattice with incomplete unit cells,we reveal that incomplete unit cells exhibit fractional charges through the analysis of Wannier centers by developing a compensation method and creating the concept of Wannier center domain(WCD)which is the smallest region that one Wannier center occupies.This method compensates for the missing parts of these boundary incomplete unit cells with additional WCDs to make them complete.The compensated WCDs automatically carry the corresponding charge,and this charge together with that of the incomplete unit cell constitutes the total charge of the complete unit cell after compensation.We conclude that the emergence of corner states is attributed to the filling anomaly,which is a fundamental mechanism.Our results refresh the understanding of HOTIs,especially those with structural discontinuities,and provide a novel design for topological states which have application value in producing optical functional devices.
基金supported in part by the National Science Foundation under Grants DMS-0807551,DMS-0720925,and DMS-0505473the Natural Science Foundationof China(10728101)supported in part by EPSRC grant EP/F029578/1
文摘We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary,subject to the kinematic and vorticity boundary conditions on the non-flat boundary.We observe that,under the nonhomogeneous boundary conditions,the pressure p can be still recovered by solving the Neumann problem for the Poisson equation.Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions.Based on this,we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition(without the incompressibility condition),which establishes a velocity mapping.Then we develop apriori estimates for the velocity mapping,especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions,which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem.Finally,we establish that,when the viscosity coefficient tends zero,the strong solutions of the initial-boundary value problem in R^n(n≥3)with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.
文摘Due to the coexistence of compressibility,viscosity,and threedimensional effects,laminar flow is difficult to maintain for high-speed boundary layer on complex geometries.The unstable disturbance waves in the boundary layer are excited and rapidly increase during the receptivity process,so sufficiently large Reynolds stress causes the basic flow velocity profile to change,and the formation of turbulence is inevitable.
基金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.
基金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.
基金supported by the Key Project of the NSFC(12131010)the NSFC(12271032)+2 种基金supported by the NSFC(12371205)the NSF of Guangdong Province(2025A1515012026)supported by the NSF of Guangdong Province(2024A1515013238)。
文摘In this paper,we consider incompressible Navier-Stokes/Cahn-Hilliard system with the generalized Navier boundary condition and the dynamic boundary condition in a channel,which can describe the interaction between a binary material and the walls of the physical domain.We prove the global-in-time existence and uniqueness of strong solutions to this initial boundary value problem in a 2D channel domain.
基金Project supported by the National Natural Science Foundation of China(No.12402262)。
文摘Vorticity is locally generated on a boundary at a rate measured by the boundary vorticity flux(BVF),which can be further decomposed into the sum of the orbital rotation and the generalized spin(specifically,the sum of shear and streaming vorticity)under the field description.For an incompressible viscous flow interacting with a stationary wall,the full expressions of the boundary fluxes of the orbital rotation and the spin are derived,for the first time,to elucidate their boundary creation mechanisms.Then,these new findings are successfully extended to the study of the boundary enstrophy dynamics,as well as the Lyman vorticity dynamics as an alternative interpretation to the boundary vorticity dynamics.Interestingly,it is found that the boundary coupling of the longitudinal and transverse processes is only embodied in the boundary spin flux,which is definitely not responsible for the generation of the boundary orbital-rotation flux.In addition,the boundary fluxes of enstrophy are directly associated with the boundary source of the second principal invariant of the velocity gradient tensor(VGT)and the two quadratic forms representing the spin-geometry interaction.The present exposition provides a new perspective and an additional dimension for understanding the vorticity dynamics on boundaries,which could be valuable in clarifying the formation mechanisms of near-wall coherent structures and flow noise at the fundamental level.
基金co-supported by the National Natural Science Foundation of China (No. 12172175)the National Science and Technology Major Project, China (No. J2019-II0014-0035)the Science Center for Gas Turbine Project, China (Nos. P2022-C-II-002-001, P2022-A-II-002-001)
文摘Cowl-induced incident Shock Wave/Boundary Layer Interactions (SWBLI) under the influence of gradual expansion waves are frequently observed in supersonic inlets. However, the analysis and prediction of interaction lengths have not been sufficiently investigated. First, this study presents a theoretical scaling analysis and validates it through wind tunnel experiments. It conducts detailed control volume analysis of mass conservation, considering the differences between inviscid and viscous cases. Then, three models for analysing interaction length under gradual expansion waves are derived. Related experiments using schlieren photography are conducted to validate the models in a Mach 2.73 flow. The interaction scales are captured at various relative distances between the shock impingement location and the expansion regions with wedge angles ranging from 12° to 15° and expansion angles of 9°, 12°, and 15°. Three trend lines are plotted based on different expansion angles to depict the relationship between normalised interaction length and normalised interaction strength metric. In addition, the relationship between the coefficients of the trend line and the expansion angles is introduced to predict the interaction length influenced by gradual expansion waves. Finally, the estimation of normalised interaction length is derived for various coefficients within a unified form.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDB0760200)the National Natural Science Foundation of China(Grant Nos.42330605 and 42377101).
文摘The Beijing 325 m meteorological tower stands as a pivotal research platform for exploring atmospheric boundary layer physics and atmospheric chemistry.With a legacy spanning 45 years,the tower has played a crucial role in unraveling the complexities of urban air pollution,atmospheric processes,and climate change in Beijing,China.This review paper provides a comprehensive overview of the measurements on the tower over the past two decades.Through long-term comprehensive observations,researchers have elucidated the intricate relationships between anthropogenic emissions,meteorological dynamics,and atmospheric composition,shedding light on the drivers of air pollution and its impacts on public health.The vertical measurements on the tower also enable detailed investigations into boundary layer dynamics,turbulent mixing,and pollutant dispersion,providing invaluable data for validating chemical transport models.Key findings from the tower’s research include the identification of positive feedback mechanisms between aerosols and the boundary layer,the characterization of pollutant sources and transport pathways,the determination of fluxes of gaseous and particulate species,and the assessment of the effectiveness of pollution control measures.Additionally,isotopic measurements have provided new insights into the sources and formation processes of particulate matter and reactive nitrogen species.Finally,the paper outlines future directions for tower-based research,emphasizing the need for long-term comprehensive measurements,the development of innovative tower platforms,and integration of emerging technologies.
基金Supported by National Natural Science Foundation of China(12101482)Postdoctoral Science Foundation of China(2022M722604)+2 种基金General Project of Science and Technology of Shaanxi Province(2023-YBSF-372)The Natural Science Foundation of Shaan Xi Province(2023-JCQN-0016)Shannxi Mathmatical Basic Science Research Project(23JSQ042)。
文摘In order to better describe the phenomenon of biological invasion,this paper introduces a free boundary model of biological invasion.Firstly,the right free boundary is added to the equation with logistic terms.Secondly,the existence and uniqueness of local solutions are proved by the Sobolev embedding theorem and the comparison principle.Finally,according to the relevant research data and contents of red fire ants,the diffusion area and nest number of red fire ants were simulated without external disturbance.This paper mainly simulates the early diffusion process of red fire ants.In the early diffusion stage,red fire ants grow slowly and then spread over a large area after reaching a certain number.
基金supported by the National Key Research and Development Program of China(Nos.2021YFB3503003,2021YFB3503100,and 2022YFB3505401).
文摘The combination of dual-main-phase(DMP)(Nd,Ce)-Fe-B magnets and grain boundary diffusion process(GBDP)is currently a research topic for obtaining high-cost performance materials in rare earth permanent magnet fields.The novel structural features of GBDP(Nd,Ce)-Fe-B magnets give a version of different domain reversal processes from those of non-diffused magnets.In this work,the in-situ magnetic domain evolution of the DMP magnets was observed at elevated temperatures,and the temperature demagnetization and coercivity mechanism of the GBDP dual-main-phase(Nd,Ce)-Fe-B magnets are discussed.The results show that the shell composition of different types of grains in DMP magnets is similar,while the magnetic microstructure results indicate the Ce-rich grains tend to demagnetize first.Dy-rich shell with a high anisotropic field caused by GBDP leads to an increase in the nucleation field,which enhances the coercivity.It is found that much more grains exhibit single domain characteristics in the remanent state for GBDP dual-main-phase(Nd,Ce)-Fe-B magnets.In addition,the grains that undergo demagnetization first are Ce-rich or Nd-rich grains,which is different from that of non-diffused magnets.These results were not found in previous studies but can be intuitively characterized from the perspective of magnetic domains in this work,providing a new perspective and understanding of the performance improvement of magnetic materials.
基金Project supported by the National Natural Science Foundation of China(52361033)National Key Research and Development Program(2022YFB3505400)+3 种基金Ministry of Industry and Information Technology Heavy Rare Earth Special Use of Sintered NdFeB Project(TC220H06J)Academic and Technical Leaders in Major Disciplines in Jiangxi Province(2022BCJ23007)Jiangxi Province Science and Technology Cooperation Key Project(20212BDH80007)Jiangxi Graduate Student Innovation Special Fund Project(YC2023-B213)。
文摘Grain boundary diffusion technology is pivotal in the preparation of high-performance NdFeB magnets.This study investigates the factors that affect the efficiency of grain boundary diffusion,starting from the properties of the diffusion matrix.Through the adjustment of the sintering process,we effectively prepared magnets with varied densities that serve as the matrix for grain boundary diffusion with TbH,diffusion.The mobility characteristics of the Nd-rich phase during the densification stage are leveraged to ensure a more extensive distribution of heavy rare earth elements within the magnets.According to the experimental results,the increase in coercivity of low-density magnets after diffusion is significantly greater than that of relatively high-density magnets.The coercivity values measured are 805.32 kA/m for low-density magnets and 470.3 kA/m for high-density magnets.Additionally,grain boundary diffusion notably enhances the density of initial low-density magnets,addressing the issue of low density during the sintering stage.Before the diffusion treatment,the Nd-rich phases primarily concentrate at the triangular grain boundaries,resulting in an increased number of cavity defects in the magnets.These cavity defects contain atoms in a higher energy state,making them more prone to transition.Consequently,the diffusion activation energy at the void defects is lower than the intracrystalline diffusion activation energy,accelerating atom diffusion.The presence of larger cavities also provides more space for atom migration,thereby promoting the diffusion process.After the diffusion treatment,the proportion of bulk Nd-rich phases significantly decreases,and they infiltrate between the grains to fill the cavity defects,forming continuous fine grain boundaries.Based on these observations,the study aims to explore how to utilize this information to develop an efficient technique for grain boundary diffusion.
基金supported by the Scientific and Technological Developing Scheme of Jilin Province under grants no.YDZJ202301ZYTS538the Chinese Academy of Sciences Youth Innovation Promotion Association under grants number 2023234+3 种基金the National Natural Science Foundation of China under grants number U21A20323the Scientific and Technological Developing Scheme of Jilin Province under grants no.SKL202302038the Major Scientific and Technological Projects of Hebei Province under grants No.23291001Zthe Scientific and Technology Project of Hanjiang District.
文摘Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymmetric tilt boundaries are hard to be accounted for based on traditional theoretical models,and the corresponding solute segregation is complex.Herein,atomic structures of a specific asymmetric boundary on{1012}TBs were reveled using aberration-corrected high-angle annular dark-field scanning transmission electron microscopy(HAADF-STEM),molecular dynamics(MD)and density functional theory(DFT)simulations.Reaction between<a60>M dislocations and the{1012}TB can generate a~61°/25°asymmetric tilt boundary.The segregation of Gd and Zn atoms is closely related to the aggregateddislocations and the interfacial interstices of the asymmetric tilt boundary,which is energetically favorable in reducing the total system energy.
基金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(52201350,52201394,and 52271301)the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.
基金Project supported by the National Key Research and Development Program of China(2021YFB3503003,2021YFB3503100,2022YFB3505401)。
文摘The microstructure of(Nd,Ce)-Fe-B sintered magnets with different diffusion depths was characterized by a magnetic force microscope,and the relationship between the magnetic properties and the local structure of grain boundary diffused magnets is discussed.The domains perpendicular to the c-axis(easy magnetization direction)show a typical maze-like pattern,while those parallel to the c-axis show the characte ristics of plate domains.The significant gradient change is shown in the concentration of Dy with the direction of diffusion from the surface to the interior.Dy diffuses along grain boundaries and(Dy,Nd)_(2)Fe_(14)B layer with a high anisotropy field formed around the grains.Through in-situ electron probe micro-analysis/magnetic force microscopy(EPMA/MFM),it is found that the average domain width decreases,and the proportion of single domain grains increases as diffusion depth increases.This is caused by both the change of concentration and distribution of Dy.The grain boundary diffusion process changes the microstructure and microchemistry inside the magnet,and these local magnetism differences can be reflected by the configuration of the magnetic domain structure.
基金supported by the National Key Research and Development Program of China(2022YFB3505503)the National Natural Science Foundation of China(52201230)+2 种基金the Key R&D Program of Shandong Province(2022CXGC020307)the China Postdoctoral Science Foundation(2022M71204)the Beijing NOVA Program(Z211100002121092).
文摘The grain boundary diffusion process(GBDP)has proven to be an effective method for enhancing the coercivity of sintered Nd-Fe-B magnets.However,the limited diffusion depth and thicker shell struc-ture have impeded the further development of magnetic properties.Currently,the primary debates re-garding the mechanism of GBDP with Tb revolve around the dissolution-solidification mechanism and the atomic substitution mechanism.To clarify this mechanism,the microstructure evolution of sintered Nd-Fe-B magnets during the heating process of GBDP has been systematically studied by quenching at different tem peratures.In this study,it was found that the formation of TbFe_(2) phase is related to the dis-solution of _(2)Fe_(14)B grains during GBDP with Tb.The theory of mixing heat and phase separation further confirms that the Nd_(2)Fe_(14)B phase dissolves to form a mixed phase of Nd and TbFe_(2),which then solidifies into the(Nd,Tb)_(2)Fe_(14)B phase.Based on the discovery of the TbFe_(2) phase,the dissolution-solidification mechanism is considered the primary mechanism for GBDP.This is supported by the elemental content of the two typical core-shell structures observed.
