In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary c...In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary conditions on the flat bottom.Assuming“Reynolds roughness regime”,in which the thickness of the domain is very small compared to the wavelength of the roughness(i.e.a very slight roughness),we rigorously derive a generalized Reynolds equation for pressure,clearly showing the roughness-induced effects.Moreover,we give expressions for the average velocity and microrotation.展开更多
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.展开更多
Plastic strain in polycrystalline metals is highly localized in grain boundaries(GBs),slip bands(SBs)and twins.While extensive research has focused on intra-granular deformation mechanisms such as slip and twinning,st...Plastic strain in polycrystalline metals is highly localized in grain boundaries(GBs),slip bands(SBs)and twins.While extensive research has focused on intra-granular deformation mechanisms such as slip and twinning,strain localization at GBs has been largely overlooked.In this study,high-resolution digital image correlation(HRDIC)was employed to capture the strain distribution and its evolution during tension in an extruded pure Mg sheet.Particular attention was paid to strain localization at GBs and its governing factors.Results reveal that,at 3%applied strain,approximately 10%of GBs were categorized as extremely-high-strain GBs(defined as the GB where at least 20 data points have an effective shear strain(ε_(eff))value exceeding the 99th percentile of the overallεeff distribution),and the majority(84%)of them were observed to deform at even 0.5%applied strain.This suggests that early-stage deformation plays a critical role in subsequent GB strain localization.The mean strain value and grain boundary sliding(GBS)displacement of GBs increased significantly with applied strain,with progressively accelerating increasing rates observed in most instances.Most(~62%)GBs exhibiting slip transfer showed low strain,while a small fraction(~8%)of them exhibited extremely high strain.This indicates that slip transfer can mitigate GB strain localization in most cases.However,complex local conditions are also critical,and case-by-case analysis is essential.Moreover,GBs with misorientation angles ranging from 50°to 80°were found to be more likely to exhibit extremely high strain.This work provides valuable insights into GB strain localization,which is critical for further understanding the plastic deformation of polycrystalline Mg.展开更多
In this paper,we study the issue of controlling a rotating flexible body-beam system(RFBBS)which consists of a tip mass attached to the free-end and a rigid disk attached to the clamped-end of an Euler-Bernoulli beam....In this paper,we study the issue of controlling a rotating flexible body-beam system(RFBBS)which consists of a tip mass attached to the free-end and a rigid disk attached to the clamped-end of an Euler-Bernoulli beam.The boundary control input is affected by both unknown disturbance and nonlinear input backlash.First,the input backlash is considered as desired control input combined with a nonlinear input error,converting it to an external disturbance,and then,the control signal is designed through the energy-based control method.Next,the closed-loop system’s stability is analysed through Lyapunov direct method.Finally,the efficacy of the proposed control scheme is tested through numerical simulations utilizing the finite difference method.展开更多
In this paper,we study the nonlinear Riemann boundary value problem with square roots that is represented by a Cauchy-type integral with kernel density in variable exponent Lebesgue spaces.We discuss the odd-order zer...In this paper,we study the nonlinear Riemann boundary value problem with square roots that is represented by a Cauchy-type integral with kernel density in variable exponent Lebesgue spaces.We discuss the odd-order zero-points distribution of the solutions and separate the single valued analytic branch of the solutions with square roots,then convert the problem to a Riemann boundary value problem in variable exponent Lebesgue spaces and discuss the singularity of solutions at individual zeros belonging to curve.We consider two types of cases those where the coefficient is Hölder and those where it is piecewise Hölder.Then we solve the Hilbert boundary value problem with square roots in variable exponent Lebesgue spaces.By discussing the distribution of the odd-order zero-points for solutions and the method of symmetric extension,we convert the Hilbert problem to a Riemann boundary value problem.The equivalence of the transformation is discussed.Finally,we get the solvable conditions and the direct expressions of the solutions in variable exponent Lebesgue spaces.展开更多
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.展开更多
As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises s...As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.展开更多
Physics-informed neural networks(PINNs)have been shown as powerful tools for solving partial differential equations(PDEs)by embedding physical laws into the network training.Despite their remarkable results,complicate...Physics-informed neural networks(PINNs)have been shown as powerful tools for solving partial differential equations(PDEs)by embedding physical laws into the network training.Despite their remarkable results,complicated problems such as irregular boundary conditions(BCs)and discontinuous or high-frequency behaviors remain persistent challenges for PINNs.For these reasons,we propose a novel two-phase framework,where a neural network is first trained to represent shape functions that can capture the irregularity of BCs in the first phase,and then these neural network-based shape functions are used to construct boundary shape functions(BSFs)that exactly satisfy both essential and natural BCs in PINNs in the second phase.This scheme is integrated into both the strong-form and energy PINN approaches,thereby improving the quality of solution prediction in the cases of irregular BCs.In addition,this study examines the benefits and limitations of these approaches in handling discontinuous and high-frequency problems.Overall,our method offers a unified and flexible solution framework that addresses key limitations of existing PINN methods with higher accuracy and stability for general PDE problems in solid mechanics.展开更多
In this article,we show the existence,uniqueness and stability of bounded solutions to the following quasilinear problems with mean curvature operator(φ'(x′(t)))′=f(t,x),t≥t_(0),lim_(t→∞)x(t)=ψ_(0),lim_(t→...In this article,we show the existence,uniqueness and stability of bounded solutions to the following quasilinear problems with mean curvature operator(φ'(x′(t)))′=f(t,x),t≥t_(0),lim_(t→∞)x(t)=ψ_(0),lim_(t→∞)x′(t)e^(t)=0,where t_(0) and ψ_(0) are real constants,φ(s)=s/√1−s^(2),s∈R with s∈(−1,1),f:[t_(0),∞)×R→R satisfies the Lipschitz or Osgood-type conditions.展开更多
This study investigates the influence of hydrogen concentration at grain boundaries on the sensitivity of polycrystalline iron to hydrogen embrittlement using molecular dynamics simulations.These simulations reveal th...This study investigates the influence of hydrogen concentration at grain boundaries on the sensitivity of polycrystalline iron to hydrogen embrittlement using molecular dynamics simulations.These simulations reveal the diffusion behavior of hydrogen atoms at grain boundaries and their consequential impact on the hydrogen embrittlement sensitivity of iron alloys.The findings indicate that as the hydrogen concentration increases,both the yield strength and ultimate tensile strength of Fe-H alloys exhibit a declining trend.Moreover,the capture of hydrogen atoms at the grain boundaries significantly influences the fracture toughness of the material and promotes the formation and propagation of cracks.This study provides a novel theoretical basis for understanding and predicting the hydrogen embrittlement behavior of iron-based materials in hydrogen-rich environments,offering valuable insights for the design and development of Fe alloys with enhanced resistance to hydrogen embrittlement.展开更多
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.展开更多
Unmanned aerial vehicles(UAVs)are widely utilized in area coverage tasks due to their flexibility and efficiency in geo-graphic information acquisition.However,complex boundary conditions in actual water area maps oft...Unmanned aerial vehicles(UAVs)are widely utilized in area coverage tasks due to their flexibility and efficiency in geo-graphic information acquisition.However,complex boundary conditions in actual water area maps often reduce coverage efficiency.To address this issue,this paper proposes a map preprocessing algorithm that linearizes boundary lines and processes concave areas into concave polygons,followed by gridding the map.Additionally,a collaborative area coverage method for UAV swarms is introduced based on region partitioning,which considers the comprehensive cost of energy consumption and time.An improved Hungarian algorithm is utilized for region partitioning,and a Dubins-A*-based plow-ing area full coverage path planning method is proposed to achieve path smoothing and collaborative coverage of each partition.Two sets of simulation experiments are conducted.The first experiment verifies the effectiveness of the map preprocessing algorithm,and the second compares the proposed collaborative area coverage algorithm with other methods,demonstrating its performance advantages.展开更多
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.展开更多
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.展开更多
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.展开更多
文摘In this paper,we study the asymptotic behavior of the micropolar fluid flow through a thin domain,assuming zero Dirichlet boundary condition on the top boundary,which is rapidly oscillating,and non-standard boundary conditions on the flat bottom.Assuming“Reynolds roughness regime”,in which the thickness of the domain is very small compared to the wavelength of the roughness(i.e.a very slight roughness),we rigorously derive a generalized Reynolds equation for pressure,clearly showing the roughness-induced effects.Moreover,we give expressions for the average velocity and microrotation.
基金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 Natural Science Foundation of China(Nos.52571157,52171125 and 52301152)Sichuan Science and Technology Program(No.2024NSFSC0193).
文摘Plastic strain in polycrystalline metals is highly localized in grain boundaries(GBs),slip bands(SBs)and twins.While extensive research has focused on intra-granular deformation mechanisms such as slip and twinning,strain localization at GBs has been largely overlooked.In this study,high-resolution digital image correlation(HRDIC)was employed to capture the strain distribution and its evolution during tension in an extruded pure Mg sheet.Particular attention was paid to strain localization at GBs and its governing factors.Results reveal that,at 3%applied strain,approximately 10%of GBs were categorized as extremely-high-strain GBs(defined as the GB where at least 20 data points have an effective shear strain(ε_(eff))value exceeding the 99th percentile of the overallεeff distribution),and the majority(84%)of them were observed to deform at even 0.5%applied strain.This suggests that early-stage deformation plays a critical role in subsequent GB strain localization.The mean strain value and grain boundary sliding(GBS)displacement of GBs increased significantly with applied strain,with progressively accelerating increasing rates observed in most instances.Most(~62%)GBs exhibiting slip transfer showed low strain,while a small fraction(~8%)of them exhibited extremely high strain.This indicates that slip transfer can mitigate GB strain localization in most cases.However,complex local conditions are also critical,and case-by-case analysis is essential.Moreover,GBs with misorientation angles ranging from 50°to 80°were found to be more likely to exhibit extremely high strain.This work provides valuable insights into GB strain localization,which is critical for further understanding the plastic deformation of polycrystalline Mg.
基金supported in part by the National Natural Science Fundation of China under Grant Nos.62403263 and 62373207in part by the Natural Science Fundation of Qingdao,China under Grant No.24-4-4-zrjj-88-jch+1 种基金in part by the Team Plan for Youth Innovation of Universities in Shandong Province under Grant No.2024KJH148in part by the Foundation of Key Laboratory of Autonomous Systems and Networked Control(South China University of Technology),Ministry of Education under Grant No.2024A01.
文摘In this paper,we study the issue of controlling a rotating flexible body-beam system(RFBBS)which consists of a tip mass attached to the free-end and a rigid disk attached to the clamped-end of an Euler-Bernoulli beam.The boundary control input is affected by both unknown disturbance and nonlinear input backlash.First,the input backlash is considered as desired control input combined with a nonlinear input error,converting it to an external disturbance,and then,the control signal is designed through the energy-based control method.Next,the closed-loop system’s stability is analysed through Lyapunov direct method.Finally,the efficacy of the proposed control scheme is tested through numerical simulations utilizing the finite difference method.
基金supported by the National Natural Science Foundation of China(11601525)the Natural Science Foundation of Hunan Province(2024JJ5412),the Changsha Municipal Natural Science Foundation(kq2402193).
文摘In this paper,we study the nonlinear Riemann boundary value problem with square roots that is represented by a Cauchy-type integral with kernel density in variable exponent Lebesgue spaces.We discuss the odd-order zero-points distribution of the solutions and separate the single valued analytic branch of the solutions with square roots,then convert the problem to a Riemann boundary value problem in variable exponent Lebesgue spaces and discuss the singularity of solutions at individual zeros belonging to curve.We consider two types of cases those where the coefficient is Hölder and those where it is piecewise Hölder.Then we solve the Hilbert boundary value problem with square roots in variable exponent Lebesgue spaces.By discussing the distribution of the odd-order zero-points for solutions and the method of symmetric extension,we convert the Hilbert problem to a Riemann boundary value problem.The equivalence of the transformation is discussed.Finally,we get the solvable conditions and the direct expressions of the solutions in variable exponent Lebesgue spaces.
基金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 National Natural Science Foundation of China(No.12372233)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University,China(No.25GH01020005)the“111 Project”of China(No.B17037)。
文摘As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.
基金Project supported by the Basic Science Research Program through the National Research Foundation(NRF)of Korea funded by the Ministry of Science and ICT(No.RS-2024-00337001)。
文摘Physics-informed neural networks(PINNs)have been shown as powerful tools for solving partial differential equations(PDEs)by embedding physical laws into the network training.Despite their remarkable results,complicated problems such as irregular boundary conditions(BCs)and discontinuous or high-frequency behaviors remain persistent challenges for PINNs.For these reasons,we propose a novel two-phase framework,where a neural network is first trained to represent shape functions that can capture the irregularity of BCs in the first phase,and then these neural network-based shape functions are used to construct boundary shape functions(BSFs)that exactly satisfy both essential and natural BCs in PINNs in the second phase.This scheme is integrated into both the strong-form and energy PINN approaches,thereby improving the quality of solution prediction in the cases of irregular BCs.In addition,this study examines the benefits and limitations of these approaches in handling discontinuous and high-frequency problems.Overall,our method offers a unified and flexible solution framework that addresses key limitations of existing PINN methods with higher accuracy and stability for general PDE problems in solid mechanics.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12361040,12061064)the National Science Foundation of Gansu Province(Grant No.22JR5RA264)State Scholarship Fund(Grant No.20230862021).
文摘In this article,we show the existence,uniqueness and stability of bounded solutions to the following quasilinear problems with mean curvature operator(φ'(x′(t)))′=f(t,x),t≥t_(0),lim_(t→∞)x(t)=ψ_(0),lim_(t→∞)x′(t)e^(t)=0,where t_(0) and ψ_(0) are real constants,φ(s)=s/√1−s^(2),s∈R with s∈(−1,1),f:[t_(0),∞)×R→R satisfies the Lipschitz or Osgood-type conditions.
基金supported by the National Science Fund for Distinguished Young Scholars(No.52425404).
文摘This study investigates the influence of hydrogen concentration at grain boundaries on the sensitivity of polycrystalline iron to hydrogen embrittlement using molecular dynamics simulations.These simulations reveal the diffusion behavior of hydrogen atoms at grain boundaries and their consequential impact on the hydrogen embrittlement sensitivity of iron alloys.The findings indicate that as the hydrogen concentration increases,both the yield strength and ultimate tensile strength of Fe-H alloys exhibit a declining trend.Moreover,the capture of hydrogen atoms at the grain boundaries significantly influences the fracture toughness of the material and promotes the formation and propagation of cracks.This study provides a novel theoretical basis for understanding and predicting the hydrogen embrittlement behavior of iron-based materials in hydrogen-rich environments,offering valuable insights for the design and development of Fe alloys with enhanced resistance to hydrogen embrittlement.
基金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.
基金National Natural Science Foundation of China(62402020,62303022)Beijing Nova Program(20240484720)+1 种基金Project of Cultivation for Young Top-Notch Talents of Beijing Municipal Institutions(BPHR202203043)BTBU Digital Business Platform Project byBMEC.
文摘Unmanned aerial vehicles(UAVs)are widely utilized in area coverage tasks due to their flexibility and efficiency in geo-graphic information acquisition.However,complex boundary conditions in actual water area maps often reduce coverage efficiency.To address this issue,this paper proposes a map preprocessing algorithm that linearizes boundary lines and processes concave areas into concave polygons,followed by gridding the map.Additionally,a collaborative area coverage method for UAV swarms is introduced based on region partitioning,which considers the comprehensive cost of energy consumption and time.An improved Hungarian algorithm is utilized for region partitioning,and a Dubins-A*-based plow-ing area full coverage path planning method is proposed to achieve path smoothing and collaborative coverage of each partition.Two sets of simulation experiments are conducted.The first experiment verifies the effectiveness of the map preprocessing algorithm,and the second compares the proposed collaborative area coverage algorithm with other methods,demonstrating its performance advantages.
文摘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.
基金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.
基金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.
