The microstructures of titanium（Ti）, an attractive tritium（T） storage material, will affect the evolution process of the retained helium（He）. Understanding the diffusion behavior of He at the atomic scale is cru...The microstructures of titanium（Ti）, an attractive tritium（T） storage material, will affect the evolution process of the retained helium（He）. Understanding the diffusion behavior of He at the atomic scale is crucial for the mechanism of material degradation. The novel diffusion behavior of He has been reported by molecular dynamics（MD） simulation for the bulk hcp-Ti system and the system with grain boundary（GB）. It is observed that the diffusion of He in the bulk hcp-Ti is significantly anisotropic（the diffusion coefficient of the [0001] direction is higher than that of the basal plane）,as represented by the different migration energies. Different from convention, the GB accelerates the diffusion of He in one direction but not in the other. It is observed that a twin boundary（TB） can serve as an effective trapped region for He.The TB accelerates diffusion of He in the direction perpendicular to the twinning direction（TD）, while it decelerates the diffusion in the TD. This finding is attributable to the change of diffusion path caused by the distortion of the local favorable site for He and the change of its number in the TB region.展开更多
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.展开更多
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.展开更多
Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to th...Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to the concealment nature of interfacial interactions.This study establishes an equivalent shear model for a bolt-resin-rock anchoring system and conducts direct shear tests under dynamic normal load(DNL)boundary from both laboratory experiments and discrete element method(DEM)simulations.The research investigates the influence of normal dynamic load amplitude(An)and rock type on shear strength parameters,elucidating the evolutionary characteristics and underlying mechanisms of shear load and normal displacement fluctuations induced by cyclic normal loading,with maximum shear load decreasing by 36.81%to 46.94%as An increases from 10%to 70%when rock type varies from coal to limestone.Through analysis of strain field evolution,the critical impact of rock type on localization of shear failure surface is revealed,with systematic summarization of differentiated wear characteristics,failure modes,and key controlling factors associated with shear failure surface.Mesoscopic investigations enabled by DEM simulations uncover the nonuniform distribution of contact force chains within the material matrix and across the anisotropic interfaces under various DNL boundaries,clarify rock type dependent crack propagation pathways,and quantitatively assess the damage extent of shear failure surface,with the anisotropic interface damage factor increasing from 34.9%to 56.6%as An rises from 10%to 70%,and decreasing from 49.6%to 23.4%as rock type varies from coal to limestone.展开更多
The new AUV driven by multi-vectored thrusters not only has unique kinematic characteristics during the actual cruise but also exists uncertain factors such as hydrodynamic coefficients perturbation and unknown interf...The new AUV driven by multi-vectored thrusters not only has unique kinematic characteristics during the actual cruise but also exists uncertain factors such as hydrodynamic coefficients perturbation and unknown interference of tail fluid, which bring difficult to the stability of the AUV＇s control system. In order to solve the nonlinear term and unmodeled dynamics existing in the new AUV＇s attitude control and the disturbances caused by the external marine environment, a second-order sliding mode controller with double-loop structure that considering the dynamic characteristics of the rudder actuators is designed, which improves the robustness of the system and avoids the control failure caused by the problem that the design theory of the sliding mode controller does not match with the actual application conditions. In order to avoid the loss of the sliding mode caused by the amplitude and rate constraints of the rudder actuator in the new AUV＇s attitude control, the dynamic boundary layer method is used to adjust the sliding boundary layer thickness so as to obtain the best anti-chattering effects. Then the impacts of system parameters, rudder actuator＇s constraints and boundary layer on the sliding mode controller are computed and analyzed to verify the effectiveness and robustness of the sliding mode controller based on dynamic boundary layer. The computational results show that the original divergent second-order sliding mode controller can still effectively implement the AUV＇s attitude control through dynamically adjusting the sliding boundary layer thickness. The dynamic boundary layer method ensures the stability of the system and does not exceed the amplitude constraint of the rudder actuator, which provides a theoretical guidance and technical support for the control system design of the new AUV in real complex sea conditions.展开更多
Objective In geo-marine science,the generalized bottom boundary layer（BBL）represents a layer between sediments and seawater.The BBL plays an important role in geological,geobiochemical,geophysical and geotechnical r...Objective In geo-marine science,the generalized bottom boundary layer（BBL）represents a layer between sediments and seawater.The BBL plays an important role in geological,geobiochemical,geophysical and geotechnical research because it is the connection region of hydrosphere,展开更多
To study the influence of construction interfaces on dynamic characteristics of roller compacted concrete dams(RCCDs),mechanical properties of construction interfaces are firstly analyzed. Then, the viscous-spring art...To study the influence of construction interfaces on dynamic characteristics of roller compacted concrete dams(RCCDs),mechanical properties of construction interfaces are firstly analyzed. Then, the viscous-spring artificial boundary(VSAB) is adopted to simulate the radiation damping of their infinite foundations, and based on the Marc software, a simplified seismic motion input method is presented by the equivalent nodal loads. Finally, based on the practical engineering of a RCC gravity dam, effects of radiation damping and construction interfaces on the dynamic characteristics of dams are investigated in detail. Analysis results show that dynamic response of the RCC gravity dam significantly reduces about 25% when the radiation damping of infinite foundation is considered. Hot interfaces and the normal cold interfaces have little influence on the dynamic response of the RCC gravity dam.However, nonlinear fracture along the cold interfaces at the dam heel will occur under the designed earthquake if the cold interfaces are combined poorly. Therefore, to avoid the fractures along the construction interfaces under the potential super earthquakes,combination quality of the RCC layers should be significantly ensured.展开更多
The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special fe...The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.展开更多
Dynamic strain-induced transformation of the low carbon steel Q（235） at 770℃ and 850℃ leads to fine ferrite grains. The microstructure characterization and mechanism of the fine ferrite grain were studied by scann...Dynamic strain-induced transformation of the low carbon steel Q（235） at 770℃ and 850℃ leads to fine ferrite grains. The microstructure characterization and mechanism of the fine ferrite grain were studied by scanning electron microscopy （SEM）, transmission electron microscopy （TEM） and electron backscattered diffraction （EBSD） technique. The results show that strain-induced microstructure is the mixed microstructure of ferrite and pearlite, with cementite randomly distributed on ferrite grain boundaries and the grains interiors. EBSD images of grain boundaries demonstrate that high angle grain boundaries （HAGBs） are dominant in both of the deformation induced microstructures occurring below and above A（e3） , with only a few low angle grain boundaries （LAGBs） existing in the grain interiors. It implies that the dynamic strain-induced transformation （DSIT） happens above and below A（e3） temperature and has the same phase transition mechanisms. The refinement of ferrite is the cooperative effect of DSIT and continuous dynamic recrystallization （CDRX） of ferrite. Besides, DSIT is deemed as an incomplete carbon diffusion phase transition through the analysis of microstructure and the previous simulated results. The strengths of the Q（235） steel with refined ferrite and pearlite structure get doubled than the initial state without treated by DSIT and the residual stress in the refined structure is partly responsible for the ductility loss.展开更多
The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to pre...The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.展开更多
The decay estimations of the solution to an elliptic equation with dynamical boundary condition is considered.We proved that,for suitable initial datum,the energy of the solution decays "in time" exponential...The decay estimations of the solution to an elliptic equation with dynamical boundary condition is considered.We proved that,for suitable initial datum,the energy of the solution decays "in time" exponentially if p=0,whereas the decay is polynomial order if p>0.展开更多
A one dimensional model is developed for defective gap mode（DGM）with two types of boundary conditions：conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of...A one dimensional model is developed for defective gap mode（DGM）with two types of boundary conditions：conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of spectral gaps equals to the modulation factor,which is consistent with previous studies.For a periodic system with local defects introduced by the boundary conditions,it shows that the conducting-mesh-induced DGM is always well confined by spectral gaps while the conducting-sleeve-induced DGM is not.The defect location can be a useful tool to dynamically control the frequency and spatial periodicity of DGM inside spectral gaps.This controllability can be potentially applied to the interaction between gap eigenmodes and energetic particles in fusion plasmas,and optical microcavities and waveguides in photonic crystals.展开更多
The stress wave propagation law and dynamic buckling critical velocity are formulated and solved by considering a general axial connecting boundary for a slender elastic straight rod impacted by a rigid body. The infl...The stress wave propagation law and dynamic buckling critical velocity are formulated and solved by considering a general axial connecting boundary for a slender elastic straight rod impacted by a rigid body. The influence of connecting stiffness on the critical velocity is investigated with varied impactor mass and buckling time. The influences of rod length and rod mass on the critical velocity are also discussed. It is found that greater connecting stiffness leads to larger stress amplitude, and further results in lower critical velocity. It is particularly noteworthy that when the connecting stiffness is less than a certain value, dynamic buckling only occurs before stress wave reflects off the connecting end. It is also shown that longer rod with larger slenderness ratio is easier to buckle, and the critical velocity for a larger-mass rod is higher than that for a lighter rod with the same geometry.展开更多
The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary...The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.展开更多
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.展开更多
In this paper,we address exponential stabilization of transmission problem of the wave equation with dynamical boundary conditions.We consider waves passing from a medium in which the speed is a_(1)into a medium in wh...In this paper,we address exponential stabilization of transmission problem of the wave equation with dynamical boundary conditions.We consider waves passing from a medium in which the speed is a_(1)into a medium in which the speed a_(2)in the case of total internal reflection and show that such a system can be controlled by introducing both dynamical boundary control along the exterior boundary and distributed control near the transmission boundary.We obtain that the system is exponential stabilization without any restriction on a_(1),a_(2)and the transmission boundary.展开更多
A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in t...A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.展开更多
In this paper,by applying theoretical method to the governing equations of compressible viscous flow,we derive the theoretical formula of the boundary dilatation flux(BDF)on a flexible wall,which generalizes the most ...In this paper,by applying theoretical method to the governing equations of compressible viscous flow,we derive the theoretical formula of the boundary dilatation flux(BDF)on a flexible wall,which generalizes the most recent work of Mao et al.(Acta Mechanica Sinica 38(2022)321583)for a stationary wall.Different boundary sources of dilatation are explicitly identified,revealing not only the boundary generation mechanisms of vortex sound and entropy sound,but also some additional sources due to the surface vorticity,surface angular velocity,surface acceleration and surface curvature.In particular,the generation mechanism of dilatation at boundary due to the coupled divergence terms is highlighted,namely,the product of the surface velocity divergence(▽_(■B)·U)and the vorticity-induced skin friction divergence(V_(■B)·τ_(ω)).The former is attributed to the surface flexibility while the latter characterizes the footprints of near-wall coherent structures.Therefore,by properly designing the surface velocity distribution,the dilatation generation at the boundary could be controlled for practical purpose in near-wall compressible viscous flows.展开更多
Open source feld operation and manipulation（OpenFOAM）is one of the most prevalent open source computational fluid dynamics（CFD）software.It is very convenient for researchers to develop their own codes based on the...Open source feld operation and manipulation（OpenFOAM）is one of the most prevalent open source computational fluid dynamics（CFD）software.It is very convenient for researchers to develop their own codes based on the class library toolbox within OpenFOAM.In recent years,several density-based solvers within OpenFOAM for supersonic/hypersonic compressible flow are coming up.Although the capabilities of these solvers to capture shock wave have already been verifed by some researchers,these solvers still need to be validated comprehensively as commercial CFD software.In boundary layer where diffusion is the dominant transportation manner,the convective discrete schemes＇capability to capture aerothermal variables,such as temperature and heat flux,is different from each other due to their own numerical dissipative characteristics and from viewpoint of this capability,these compressible solvers within OpenFOAM can be validated further.In this paper,frstly,the organizational architecture of density-based solvers within OpenFOAM is analyzed.Then,from the viewpoint of the capability to capture aerothermal variables,the numerical results of several typical geometrical felds predicted by these solvers are compared with both the outcome obtained from the commercial software Fastran and the experimental data.During the computing process,the Roe,AUSM＋（Advection Upstream Splitting Method）,and HLLC（Harten-Lax-van Leer-Contact）convective discrete schemes of which the spatial accuracy is 1st and 2nd order are utilized,respectively.The compared results show that the aerothermal variables are in agreement with results generated by Fastran and the experimental data even if the1st order spatial precision is implemented.Overall,the accuracy of these density-based solvers can meet the requirement of engineering and scientifc problems to capture aerothermal variables in diffusion boundary layer.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant No.51501119)the Scientific Research Starting Foundation for Younger Teachers of Sichuan University,China(Grant No.2015SCU11058)+1 种基金the National Magnetic Confinement Fusion Science Program of China(Grant No.2013GB109002)the Cooperative Research Project "Research of Diffusion Behaviour of He in Grain Boundary of HCP-Titanium",China
文摘The microstructures of titanium（Ti）, an attractive tritium（T） storage material, will affect the evolution process of the retained helium（He）. Understanding the diffusion behavior of He at the atomic scale is crucial for the mechanism of material degradation. The novel diffusion behavior of He has been reported by molecular dynamics（MD） simulation for the bulk hcp-Ti system and the system with grain boundary（GB）. It is observed that the diffusion of He in the bulk hcp-Ti is significantly anisotropic（the diffusion coefficient of the [0001] direction is higher than that of the basal plane）,as represented by the different migration energies. Different from convention, the GB accelerates the diffusion of He in one direction but not in the other. It is observed that a twin boundary（TB） can serve as an effective trapped region for He.The TB accelerates diffusion of He in the direction perpendicular to the twinning direction（TD）, while it decelerates the diffusion in the TD. This finding is attributable to the change of diffusion path caused by the distortion of the local favorable site for He and the change of its number in the TB region.
基金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 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.
基金support from the National Natural Science Foundation of China(Nos.51504247,52174092,51904290,and 52074259)the Natural Science Foundation of Jiangsu Province,China(No.BK20220157)+1 种基金the Fundamental Research Funds for the Central Universities,China(No.2022YCPY0202)the China University of Mining and Technology(CUMT)Open Sharing Fund for Large-scale Instruments and Equipment(No.DYGX-2025-47)is gratefully acknowledged.
文摘Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to the concealment nature of interfacial interactions.This study establishes an equivalent shear model for a bolt-resin-rock anchoring system and conducts direct shear tests under dynamic normal load(DNL)boundary from both laboratory experiments and discrete element method(DEM)simulations.The research investigates the influence of normal dynamic load amplitude(An)and rock type on shear strength parameters,elucidating the evolutionary characteristics and underlying mechanisms of shear load and normal displacement fluctuations induced by cyclic normal loading,with maximum shear load decreasing by 36.81%to 46.94%as An increases from 10%to 70%when rock type varies from coal to limestone.Through analysis of strain field evolution,the critical impact of rock type on localization of shear failure surface is revealed,with systematic summarization of differentiated wear characteristics,failure modes,and key controlling factors associated with shear failure surface.Mesoscopic investigations enabled by DEM simulations uncover the nonuniform distribution of contact force chains within the material matrix and across the anisotropic interfaces under various DNL boundaries,clarify rock type dependent crack propagation pathways,and quantitatively assess the damage extent of shear failure surface,with the anisotropic interface damage factor increasing from 34.9%to 56.6%as An rises from 10%to 70%,and decreasing from 49.6%to 23.4%as rock type varies from coal to limestone.
基金supported by National Hi-tech Research and Development Program of China (863 Program,Grant No.2006AA09Z235)Hunan Provincial Innovation Foundation for Postgraduate of China (Grant No. CX2009B003)
文摘The new AUV driven by multi-vectored thrusters not only has unique kinematic characteristics during the actual cruise but also exists uncertain factors such as hydrodynamic coefficients perturbation and unknown interference of tail fluid, which bring difficult to the stability of the AUV＇s control system. In order to solve the nonlinear term and unmodeled dynamics existing in the new AUV＇s attitude control and the disturbances caused by the external marine environment, a second-order sliding mode controller with double-loop structure that considering the dynamic characteristics of the rudder actuators is designed, which improves the robustness of the system and avoids the control failure caused by the problem that the design theory of the sliding mode controller does not match with the actual application conditions. In order to avoid the loss of the sliding mode caused by the amplitude and rate constraints of the rudder actuator in the new AUV＇s attitude control, the dynamic boundary layer method is used to adjust the sliding boundary layer thickness so as to obtain the best anti-chattering effects. Then the impacts of system parameters, rudder actuator＇s constraints and boundary layer on the sliding mode controller are computed and analyzed to verify the effectiveness and robustness of the sliding mode controller based on dynamic boundary layer. The computational results show that the original divergent second-order sliding mode controller can still effectively implement the AUV＇s attitude control through dynamically adjusting the sliding boundary layer thickness. The dynamic boundary layer method ensures the stability of the system and does not exceed the amplitude constraint of the rudder actuator, which provides a theoretical guidance and technical support for the control system design of the new AUV in real complex sea conditions.
基金supported by the National Natural Science Foundation of China (grants No.41427803 amd 41272316)
文摘Objective In geo-marine science,the generalized bottom boundary layer（BBL）represents a layer between sediments and seawater.The BBL plays an important role in geological,geobiochemical,geophysical and geotechnical research because it is the connection region of hydrosphere,
基金Projects(20120094110005,20120094130003)supported by the Research Fund for the Doctoral Program of Higher Education of ChinaProjects(51379068,51139001,51279052,51209077,51179066)supported by the National Natural Science Foundation of China+1 种基金Project(NCET-11-0628)supported by the Program for New Century Excellent Talents in University,ChinaProjects(201201038,201101013)supported by the Public Welfare Industry Research Special Fund Project of Ministry of Water Resources of China
文摘To study the influence of construction interfaces on dynamic characteristics of roller compacted concrete dams(RCCDs),mechanical properties of construction interfaces are firstly analyzed. Then, the viscous-spring artificial boundary(VSAB) is adopted to simulate the radiation damping of their infinite foundations, and based on the Marc software, a simplified seismic motion input method is presented by the equivalent nodal loads. Finally, based on the practical engineering of a RCC gravity dam, effects of radiation damping and construction interfaces on the dynamic characteristics of dams are investigated in detail. Analysis results show that dynamic response of the RCC gravity dam significantly reduces about 25% when the radiation damping of infinite foundation is considered. Hot interfaces and the normal cold interfaces have little influence on the dynamic response of the RCC gravity dam.However, nonlinear fracture along the cold interfaces at the dam heel will occur under the designed earthquake if the cold interfaces are combined poorly. Therefore, to avoid the fractures along the construction interfaces under the potential super earthquakes,combination quality of the RCC layers should be significantly ensured.
基金The project supported by the National Natural Science Foundation of China (50579081)the Australian Research Council (DP0452681)The English text was polished by Keren Wang
文摘The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.
基金support from the National Natural Science Foundation of China (NSFC) under Grantb No. 50871109
文摘Dynamic strain-induced transformation of the low carbon steel Q（235） at 770℃ and 850℃ leads to fine ferrite grains. The microstructure characterization and mechanism of the fine ferrite grain were studied by scanning electron microscopy （SEM）, transmission electron microscopy （TEM） and electron backscattered diffraction （EBSD） technique. The results show that strain-induced microstructure is the mixed microstructure of ferrite and pearlite, with cementite randomly distributed on ferrite grain boundaries and the grains interiors. EBSD images of grain boundaries demonstrate that high angle grain boundaries （HAGBs） are dominant in both of the deformation induced microstructures occurring below and above A（e3） , with only a few low angle grain boundaries （LAGBs） existing in the grain interiors. It implies that the dynamic strain-induced transformation （DSIT） happens above and below A（e3） temperature and has the same phase transition mechanisms. The refinement of ferrite is the cooperative effect of DSIT and continuous dynamic recrystallization （CDRX） of ferrite. Besides, DSIT is deemed as an incomplete carbon diffusion phase transition through the analysis of microstructure and the previous simulated results. The strengths of the Q（235） steel with refined ferrite and pearlite structure get doubled than the initial state without treated by DSIT and the residual stress in the refined structure is partly responsible for the ductility loss.
基金Supported by the Key Program of National Natural Science Foundation of China(No.51138001)the Science Fund for Creative Research Groups of National Natural Science Foundation of China(No.51121005)+2 种基金the Fundamental Research Funds for the Central Universities(DUT13LK16)the Young Scientists Fund of National Natural Science Foundation of China(No.51109134)China Postdoctoral Science Foundation(No.2011M500814)
文摘The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.
基金Supported by the National Natural Science Foundation of China(10671182)Supported by the Natural Science Foundation of Henan Province(0611053300+1 种基金200510463024)Supported by the Young Skeleton Teacher Project of the Higher School of Henan Province
文摘The decay estimations of the solution to an elliptic equation with dynamical boundary condition is considered.We proved that,for suitable initial datum,the energy of the solution decays "in time" exponentially if p=0,whereas the decay is polynomial order if p>0.
基金supported by National Natural Science Foundation of China(No.11405271)
文摘A one dimensional model is developed for defective gap mode（DGM）with two types of boundary conditions：conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of spectral gaps equals to the modulation factor,which is consistent with previous studies.For a periodic system with local defects introduced by the boundary conditions,it shows that the conducting-mesh-induced DGM is always well confined by spectral gaps while the conducting-sleeve-induced DGM is not.The defect location can be a useful tool to dynamically control the frequency and spatial periodicity of DGM inside spectral gaps.This controllability can be potentially applied to the interaction between gap eigenmodes and energetic particles in fusion plasmas,and optical microcavities and waveguides in photonic crystals.
文摘The stress wave propagation law and dynamic buckling critical velocity are formulated and solved by considering a general axial connecting boundary for a slender elastic straight rod impacted by a rigid body. The influence of connecting stiffness on the critical velocity is investigated with varied impactor mass and buckling time. The influences of rod length and rod mass on the critical velocity are also discussed. It is found that greater connecting stiffness leads to larger stress amplitude, and further results in lower critical velocity. It is particularly noteworthy that when the connecting stiffness is less than a certain value, dynamic buckling only occurs before stress wave reflects off the connecting end. It is also shown that longer rod with larger slenderness ratio is easier to buckle, and the critical velocity for a larger-mass rod is higher than that for a lighter rod with the same geometry.
基金supported by the National Natural Science Foundation of China(Grant No.52071094).
文摘The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.
基金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.
基金supported by the National Natural Science Foundation of China(No.11671240)Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi(No.2024L300)。
文摘In this paper,we address exponential stabilization of transmission problem of the wave equation with dynamical boundary conditions.We consider waves passing from a medium in which the speed is a_(1)into a medium in which the speed a_(2)in the case of total internal reflection and show that such a system can be controlled by introducing both dynamical boundary control along the exterior boundary and distributed control near the transmission boundary.We obtain that the system is exponential stabilization without any restriction on a_(1),a_(2)and the transmission boundary.
文摘A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.
文摘In this paper,by applying theoretical method to the governing equations of compressible viscous flow,we derive the theoretical formula of the boundary dilatation flux(BDF)on a flexible wall,which generalizes the most recent work of Mao et al.(Acta Mechanica Sinica 38(2022)321583)for a stationary wall.Different boundary sources of dilatation are explicitly identified,revealing not only the boundary generation mechanisms of vortex sound and entropy sound,but also some additional sources due to the surface vorticity,surface angular velocity,surface acceleration and surface curvature.In particular,the generation mechanism of dilatation at boundary due to the coupled divergence terms is highlighted,namely,the product of the surface velocity divergence(▽_(■B)·U)and the vorticity-induced skin friction divergence(V_(■B)·τ_(ω)).The former is attributed to the surface flexibility while the latter characterizes the footprints of near-wall coherent structures.Therefore,by properly designing the surface velocity distribution,the dilatation generation at the boundary could be controlled for practical purpose in near-wall compressible viscous flows.
基金supported by the National Natural Science Foundation of China(Grant Nos.51176038,51121004)
文摘Open source feld operation and manipulation（OpenFOAM）is one of the most prevalent open source computational fluid dynamics（CFD）software.It is very convenient for researchers to develop their own codes based on the class library toolbox within OpenFOAM.In recent years,several density-based solvers within OpenFOAM for supersonic/hypersonic compressible flow are coming up.Although the capabilities of these solvers to capture shock wave have already been verifed by some researchers,these solvers still need to be validated comprehensively as commercial CFD software.In boundary layer where diffusion is the dominant transportation manner,the convective discrete schemes＇capability to capture aerothermal variables,such as temperature and heat flux,is different from each other due to their own numerical dissipative characteristics and from viewpoint of this capability,these compressible solvers within OpenFOAM can be validated further.In this paper,frstly,the organizational architecture of density-based solvers within OpenFOAM is analyzed.Then,from the viewpoint of the capability to capture aerothermal variables,the numerical results of several typical geometrical felds predicted by these solvers are compared with both the outcome obtained from the commercial software Fastran and the experimental data.During the computing process,the Roe,AUSM＋（Advection Upstream Splitting Method）,and HLLC（Harten-Lax-van Leer-Contact）convective discrete schemes of which the spatial accuracy is 1st and 2nd order are utilized,respectively.The compared results show that the aerothermal variables are in agreement with results generated by Fastran and the experimental data even if the1st order spatial precision is implemented.Overall,the accuracy of these density-based solvers can meet the requirement of engineering and scientifc problems to capture aerothermal variables in diffusion boundary layer.
