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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is int...In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.展开更多
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.展开更多
The Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions is considered.The existence of the global attractor is proved and the long time behavior of the trajectories,namely,the convergence ...The Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions is considered.The existence of the global attractor is proved and the long time behavior of the trajectories,namely,the convergence to steady states,is studied.展开更多
In this paper,we continue the earlier work[Lu,L,&Wang,D.L.(2017).Dynamic boundary feed-back of a pendulum coupled with a viscous damped wave equation.In Proceedings of the 36th Chinese Control Conference(CCQ)(pp.1...In this paper,we continue the earlier work[Lu,L,&Wang,D.L.(2017).Dynamic boundary feed-back of a pendulum coupled with a viscous damped wave equation.In Proceedings of the 36th Chinese Control Conference(CCQ)(pp.1676-1680)]on study the stability of a pendulum coupled with a viscous damped wave equation model.This time we get the exponential stability result which is much better than the previous strong stability.By a detailed spectral analysis and opera-tor separation,we establish the Riesz basis property as well as the spectrum determined growth condition for the system.Finally,the exponential stability of the system is achieved.展开更多
In this paper,the Peng-Robinson equation of state with dynamic boundary conditions is discussed,which considers the interactions with solid walls.At first,the model is introduced and the regularization method on the n...In this paper,the Peng-Robinson equation of state with dynamic boundary conditions is discussed,which considers the interactions with solid walls.At first,the model is introduced and the regularization method on the nonlinear term is adopted.Next,The scalar auxiliary variable(SAV)method in temporal and finite element method in spatial are used to handle the Peng-Robinson equation of state.Then,the energy dissipation law of the numerical method is obtained.Also,we acquire the convergence of the discrete SAV finite element method(FEM).Finally,a numerical example is provided to confirm the theoretical result.展开更多
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,展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
This paper presents a design of boundary controllers implemented at the top end for global stabilization of a marine riser in a three dimensional space under environmental loadings. Based on the energy approach, nonli...This paper presents a design of boundary controllers implemented at the top end for global stabilization of a marine riser in a three dimensional space under environmental loadings. Based on the energy approach, nonlinear partial differential equations of motion, including bending-bending and longitudinal-bending couplings for the risers are derived. The couplings cause mutual effects between the three independent directions in the riser＇s motions, and make it difficult to minimize its vibrations. The Lyapunov direct method is employed to design the boundary controller. It is shown that the proposed boundary controllers can effectively reduce the riser＇s vibration. Stability analysis of the closed-loop system is performed using the Lyapunov direct method. Numerical simulations illustrate the results.展开更多
基金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.
基金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.
基金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 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.
基金supported by NSFC(Grant No.11871441)supported by NSF(Grant No.DMS-2012669).
文摘In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.
基金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.
文摘The Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions is considered.The existence of the global attractor is proved and the long time behavior of the trajectories,namely,the convergence to steady states,is studied.
基金supported by Beijing Excellent Talents Train-ing Project Foundation and School Key Projects for Science and Technology[2017000020124G053 and 2020Z170-KXZ].
文摘In this paper,we continue the earlier work[Lu,L,&Wang,D.L.(2017).Dynamic boundary feed-back of a pendulum coupled with a viscous damped wave equation.In Proceedings of the 36th Chinese Control Conference(CCQ)(pp.1676-1680)]on study the stability of a pendulum coupled with a viscous damped wave equation model.This time we get the exponential stability result which is much better than the previous strong stability.By a detailed spectral analysis and opera-tor separation,we establish the Riesz basis property as well as the spectrum determined growth condition for the system.Finally,the exponential stability of the system is achieved.
基金the National Natural Science Foundation of China(No.11871441).
文摘In this paper,the Peng-Robinson equation of state with dynamic boundary conditions is discussed,which considers the interactions with solid walls.At first,the model is introduced and the regularization method on the nonlinear term is adopted.Next,The scalar auxiliary variable(SAV)method in temporal and finite element method in spatial are used to handle the Peng-Robinson equation of state.Then,the energy dissipation law of the numerical method is obtained.Also,we acquire the convergence of the discrete SAV finite element method(FEM).Finally,a numerical example is provided to confirm the theoretical result.
基金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,
基金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.
基金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.
文摘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.
基金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.
基金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.
基金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.
文摘This paper presents a design of boundary controllers implemented at the top end for global stabilization of a marine riser in a three dimensional space under environmental loadings. Based on the energy approach, nonlinear partial differential equations of motion, including bending-bending and longitudinal-bending couplings for the risers are derived. The couplings cause mutual effects between the three independent directions in the riser＇s motions, and make it difficult to minimize its vibrations. The Lyapunov direct method is employed to design the boundary controller. It is shown that the proposed boundary controllers can effectively reduce the riser＇s vibration. Stability analysis of the closed-loop system is performed using the Lyapunov direct method. Numerical simulations illustrate the results.
