In this paper an approximate equation is derived to describe smooth parts of the stability boundary for linear Hamiltonian systems, depending on arbitrary number of parameters. With this equation, we can obtain parame...In this paper an approximate equation is derived to describe smooth parts of the stability boundary for linear Hamiltonian systems, depending on arbitrary number of parameters. With this equation, we can obtain parameters corresponding to the stability boundary, as well as to the stability and instability domains, provided that one point on the stability boundary is known. Then differential equations describing the evolution of eigenvalues and eigenvectors along a curve on the stability boundary surface are derived. These equations also allow us to obtain curves belonging to the stability boundary. Applications to linear gyroscopic systems are considered and studied with examples.展开更多
The criterion for the points in the parameter space being on the stability boundary of linear Hamiltonian system depending on arbitrary numbers of parameters was given, through the sensitivity analysis of eigenvalues ...The criterion for the points in the parameter space being on the stability boundary of linear Hamiltonian system depending on arbitrary numbers of parameters was given, through the sensitivity analysis of eigenvalues and eigenvectors. The results show that multiple eigenvalues with Jordan chain take a very important role in the stability of Hamiltonian systems.展开更多
A single-axis flux decay model including an excitation control model proposed in [12,14,16] is studied. As the bifurcation parameter P m (input power to the generator) varies, the system exhibits dynamics emerging fro...A single-axis flux decay model including an excitation control model proposed in [12,14,16] is studied. As the bifurcation parameter P m (input power to the generator) varies, the system exhibits dynamics emerging from static and dynamic bifurcations which link with system collapse. We show that the equilibrium point of the system undergoes three bifurcations: one saddle-node bifurcation and two Hopf bifurcations. The state variables dominating system collapse are different for different critical points, and the excitative control may play an important role in delaying system from collapsing. Simulations are presented to illustrate the dynamical behavior associated with the power system stability and collapse. Moreover, by computing the local quadratic approximation of the 5-dimensional stable manifold at an order 5 saddle point, an analytical expression for the approximate stability boundary is worked out.展开更多
The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parab...The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.展开更多
Compressible boundary layers stability on blade cascade suction surface was discussed by wind tunnel experiment and numerical solution. Three dimensional disturbance wave Parabolized Stability Equations (PSE) of ortho...Compressible boundary layers stability on blade cascade suction surface was discussed by wind tunnel experiment and numerical solution. Three dimensional disturbance wave Parabolized Stability Equations (PSE) of orthogonal Curvilinear Coordinates in compressible flow was deducted. The surface pressure of blade in wind tunnel experiment was measured. The Falkner-Skan equation was solved under the boundary conditions of experiment result, and velocity, pressure and temperature of average flow were obtained. Substituted this result for discretization of the PSE Eigenvalue Problem, the stability problem can be solved.展开更多
In this paper. we discuss the existence and stability of solution for two semi-homogeneous boundary value problems. The relative theorems in [1.2] are extended. Meanwhile. we obtain some new results.
Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the MoninObukhov similarity for describing the turbulent fluctu ations of 3D win...Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the MoninObukhov similarity for describing the turbulent fluctu ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * σ v /u * σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day clear sky overcast cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "1/3" slope in unstable stratification and "2/3" slope in stable stratification.展开更多
The elementary task is to calculate the growth rates of disturbances when the e;method in transition prediction is performed. However, there is no unified knowledge to determine the growth rates of disturbances in thr...The elementary task is to calculate the growth rates of disturbances when the e;method in transition prediction is performed. However, there is no unified knowledge to determine the growth rates of disturbances in three-dimensional(3 D) flows. In this paper, we study the relation among the wave parameters of the disturbance in boundary layers in which the imaginary parts of wave parameters are far smaller than the real parts.The generalized growth rate(GGR) in the direction of group velocity is introduced, and the conservation relation of GGR is strictly deduced in theory. This conservation relation manifests that the GGR only depends on the real parts of wave parameters instead of the imaginary parts. Numerical validations for GGR conservation are also provided in the cases of first/second modes and crossflow modes. The application of GGR to the eN method in 3 D flows is discussed, and the puzzle of determining growth rates in 3 D flows is clarified. A convenient method is also proposed to calculate growth rates of disturbances in 3 D flows. Good agreement between this convenient method and existing methods is found except the condition that the angle between the group velocity direction and the x-direction is close to 90?which can be easily avoided in practical application.展开更多
The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, w...The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, which is very effective, was used to study the nonlinear evolution. of T-S disturbance wave at finite amplitudes. Solving nonlinear equations of different modes by using predictor-corrector and iterative approach, which is uncoupled between modes, improving computational accuracy by using high order compact differential scheme, satisfying normalization condition I determining tables of nonlinear terms at different modes, and implementing stably the spatial marching, were included in this method. With different initial amplitudes, the nonlinear evolution of T-S wave was studied. The nonlinear nonparallel results of examples compare with data of direct numerical simulations (DNS) using full Navier-Stokes equations.展开更多
TheMoon is the only celestial body that human beings have visited.The design of the Earth-Moon transfer orbits is a critical issue in lunar exploration missions.In the 21st century,new lunar missions including the con...TheMoon is the only celestial body that human beings have visited.The design of the Earth-Moon transfer orbits is a critical issue in lunar exploration missions.In the 21st century,new lunar missions including the construction of the lunar space station,the permanent lunar base,and the Earth-Moon transportation network have been proposed,requiring low-cost,expansive launch windows and a fixed arrival epoch for any launch date within the launch window.The low-energy and low-thrust transfers are promising strategies to satisfy the demands.This review provides a detailed landscape of Earth-Moon transfer trajectory design processes,from the traditional patched conic to the state-of-the-art low-energy and low-thrust methods.Essential mechanisms of the various utilized dynamic models and the characteristics of the different design methods are discussed in hopes of helping readers grasp thebasic overviewof the current Earth-Moon transfer orbitdesignmethods anda deep academic background is unnecessary for the context understanding.展开更多
This paper proposes a novel Multivariate Quotient-Difference(MQD)method to obtain the approximate analytical solution for AC power flow equations.Therefore,in the online environment,the power flow solutions covering d...This paper proposes a novel Multivariate Quotient-Difference(MQD)method to obtain the approximate analytical solution for AC power flow equations.Therefore,in the online environment,the power flow solutions covering different operating conditions can be directly obtained by plugging values into multiple symbolic variables,such that the power injections and consumptions of selected buses or areas can be independently adjusted.This method first derives a power flow solution through a Multivariate Power Series(MPS).Next,the MQD method is applied to transform the obtained MPS to a Multivariate Pad´e Approximants(MPA)to expand the Radius of Convergence(ROC),so that the accuracy of the derived analytical solution can be significantly increased.In addition,the hypersurface of the voltage stability boundary can be identified by an analytical formula obtained from the coefficients of MPA.This direct method for power flow solutions and voltage stability boundaries is fast for many online applications,since such analytical solutions can be derived offline and evaluated online by only plugging values into the symbolic variables according to the actual operating conditions.The proposed method is validated in detail on New England 39-bus and IEEE 118-bus systems with independent load variations in multi-regions.展开更多
Droop-controlled voltage-source converters(VSCs)can provide frequency and voltage support to power grids.However,during a grid fault,VsCs may experience transient instability,which can be significantly affected by bot...Droop-controlled voltage-source converters(VSCs)can provide frequency and voltage support to power grids.However,during a grid fault,VsCs may experience transient instability,which can be significantly affected by both the control parameters and fault conditions.This mechanism has not been fully elucidated in previous studies.In particular,grid-voltage faults are commonly accompanied by a grid voltage phase-angle jump(VPAJ),which is typically ignored in the evaluation of the transient stability of VSCs.To address this issue,this study comprehensively assesses the impact of the VPAJ and key control parameters on the transient characteristics of VSCs.Furthermore,the critical clearing angle and critical clearing time are quantitatively calculated to define the transient stability boundary.In addition,a transient stability-enhancement control method that considers the transient stability constraints is proposed.Finally,simulations and experimental tests are conducted to validate both the theoretical analysis and proposed method.展开更多
As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely use...As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.展开更多
This article studies the nonlinear evolution of disturbance waves in supersonic nonparallel boundary layer flows by using nonlinear parabolic stability equations (NPSE). An accurate numerical method is developed to ...This article studies the nonlinear evolution of disturbance waves in supersonic nonparallel boundary layer flows by using nonlinear parabolic stability equations (NPSE). An accurate numerical method is developed to solve the equations and march the NPSE in a stable manner. Through computation,are obtained the curves of amplitude and disturbance shape function of harmonic waves. Especially are demonstrated the physical characteristics of nonlinear stability of various harmonic waves,including instantaneous stream wise vortices,spanwise vortices and Λ structure etc,and are used to study and analyze the mechanism of the transition process. The calculated results have evidenced the effectiveness of the proposed NPSE method to research the nonlinear stability of the supersonic boundary layers.展开更多
Subject Code:E01 Conventional polycrystalline materials become harder with decreasing grain size,following the classical Hall—Petch relationship,i.e.,strength increases reversely proportional to the square root of th...Subject Code:E01 Conventional polycrystalline materials become harder with decreasing grain size,following the classical Hall—Petch relationship,i.e.,strength increases reversely proportional to the square root of the grain size.Strengthening occurs due to dislocation pileups at grain boundaries that prevent the dislocations展开更多
The nonparallel effects on the stability of the boundary layer flow was investigated using the Parabolie Stability Equations (PSE). In order to improve the accuracy of the calculation which is very important for the i...The nonparallel effects on the stability of the boundary layer flow was investigated using the Parabolie Stability Equations (PSE). In order to improve the accuracy of the calculation which is very important for the investigation of stability, higher order expansions in orthogonal functions in normal direction and the effective algebraic mapping to deal with the problem of infinite region were used and the way to collocate the boundary point based on its characteristics was adopted. With the effective control of step size in the marching procedure, the special condition was satisfied, and the stability of calculation was assured. From the curves of the neutral stability, the growth rate, the amplitude variation and disturbed velocity profile, the effects of the nonparallelism were given accurately and analyzed detailedly. It is found that the nonparallelism of the flow amplifies the amplitude and growth rate of disturbances, especially for three-dimensional disturbances, even can change the sign of flow stability from stability to instability for some cases. Computed results are in good agreement with the classical experimental results.展开更多
The linear instabilities of incompressible confluent mixing layer and boundary layer were analyzed.The mixing layers include wake,shear layer and their combination.The mean velocity profile of confluent flow is taken ...The linear instabilities of incompressible confluent mixing layer and boundary layer were analyzed.The mixing layers include wake,shear layer and their combination.The mean velocity profile of confluent flow is taken as a superposition of a hyperbolic and exponential function to model a mixing layer and the Blasius similarity solution for a flat plate boundary layer.The stability equation of confluent flow was solved by using the global numerical method.The unstable modes associated with both the mixing and boundary layers were identified.They are the boundary layer mode,mixing layer mode 1(nearly symmetrical mode)and mode 2(nearly anti-symmetrical mode).The interactions between the mixing layer stability and the boundary layer stability were examined.As the mixing layer approaches the boundary layer,the neutral curves of the boundary layer mode move to the upper left,the resulting critical Reynolds number decreases,and the growth rate of the most unstable mode increases.The wall tends to stabilize the mixing layer modes at low frequency.In addition,the mode switching behavior of the relative level of the spatial growth rate between the mixing layer mode 1 and mode 2 with the velocity ratio is found to occur at low frequency.展开更多
We study the stability of unduloids with free boundary in the domain B between two parallel hyperplanes in R^(n+1). If the unduloid has one half of period in B and is sufficiently close to a cylinder, then for 2≤n≤1...We study the stability of unduloids with free boundary in the domain B between two parallel hyperplanes in R^(n+1). If the unduloid has one half of period in B and is sufficiently close to a cylinder, then for 2≤n≤10, it is unstable; while for n≥11, it is stable. If the unduloid has two or more halves of period in B and is sufficiently close to a cylinder, then for all n≥2, it is unstable.展开更多
In this paper, we consider the boundary stabilization of the wave equation with variable coefficients by Riemmannian geometry method subject to a different geometric condition which is motivated by the geometric multi...In this paper, we consider the boundary stabilization of the wave equation with variable coefficients by Riemmannian geometry method subject to a different geometric condition which is motivated by the geometric multiplier identities. Several (multiplier) identities (inequalities) which have been built for constant wave equation by Kormornik and Zuazua are generalized to the variable coefficient case by some computational techniques in Riemmannian geometry, so that the precise estimates on the exponential decay rate are derived from those inequalitities. Also, the exponential decay for the solutions of semilinear wave equation with variable coefficients is obtained under natural growth and sign assumptions on the nonlinearity. Our method is rather general and can be adapted to other evolution systems with variable coefficients (e.g. elasticity plates) as well.展开更多
In this study,the transverse vibration of a traveling beam made of functionally graded material was analyzed.The material gradation was assumed to vary continuously along the thickness direction of the beam in the for...In this study,the transverse vibration of a traveling beam made of functionally graded material was analyzed.The material gradation was assumed to vary continuously along the thickness direction of the beam in the form of power law exponent.The effect of the longitudinally varying tension due to axial acceleration was highlighted,and the dependence of the tension on the finite support rigidity was also considered.A complex governing equation of the functionally graded beam was derived by the Hamilton principle,in which the geometric nonlinearity,material properties and axial load were incorporated.The direct multiscale method was applied to the analysis process of an axially moving functionally graded beam with timedependent velocity,and the natural frequency and solvability conditions were obtained.Based on the conditions,the stability boundaries of subharmonic resonance and combination resonance were obtained.It was found that the dynamic behavior of axial moving beams could be tuned by using the distribution law of the functional gradient parameters.展开更多
基金The project supported by the National Science Foundations of Russia and China (10072012)
文摘In this paper an approximate equation is derived to describe smooth parts of the stability boundary for linear Hamiltonian systems, depending on arbitrary number of parameters. With this equation, we can obtain parameters corresponding to the stability boundary, as well as to the stability and instability domains, provided that one point on the stability boundary is known. Then differential equations describing the evolution of eigenvalues and eigenvectors along a curve on the stability boundary surface are derived. These equations also allow us to obtain curves belonging to the stability boundary. Applications to linear gyroscopic systems are considered and studied with examples.
文摘The criterion for the points in the parameter space being on the stability boundary of linear Hamiltonian system depending on arbitrary numbers of parameters was given, through the sensitivity analysis of eigenvalues and eigenvectors. The results show that multiple eigenvalues with Jordan chain take a very important role in the stability of Hamiltonian systems.
基金Supported by the National Key Basic Research Fund (No.G1998020307)KZCX-2-SW-118 Chinese Academy of Sciences
文摘A single-axis flux decay model including an excitation control model proposed in [12,14,16] is studied. As the bifurcation parameter P m (input power to the generator) varies, the system exhibits dynamics emerging from static and dynamic bifurcations which link with system collapse. We show that the equilibrium point of the system undergoes three bifurcations: one saddle-node bifurcation and two Hopf bifurcations. The state variables dominating system collapse are different for different critical points, and the excitative control may play an important role in delaying system from collapsing. Simulations are presented to illustrate the dynamical behavior associated with the power system stability and collapse. Moreover, by computing the local quadratic approximation of the 5-dimensional stable manifold at an order 5 saddle point, an analytical expression for the approximate stability boundary is worked out.
文摘The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.
文摘Compressible boundary layers stability on blade cascade suction surface was discussed by wind tunnel experiment and numerical solution. Three dimensional disturbance wave Parabolized Stability Equations (PSE) of orthogonal Curvilinear Coordinates in compressible flow was deducted. The surface pressure of blade in wind tunnel experiment was measured. The Falkner-Skan equation was solved under the boundary conditions of experiment result, and velocity, pressure and temperature of average flow were obtained. Substituted this result for discretization of the PSE Eigenvalue Problem, the stability problem can be solved.
文摘In this paper. we discuss the existence and stability of solution for two semi-homogeneous boundary value problems. The relative theorems in [1.2] are extended. Meanwhile. we obtain some new results.
基金provided by the Natural Science Fund for Universities in Jiangsu Province (Grant No 08KJA170002)the Meteorology Fund of the Ministry of Science and Technology [Grant No GYHY(QX) 2007-6-26]+3 种基金the National Natural Science Foundation of China (Grant No 40775012)the National Key Technology R&D Program (Grant No 2008BAC48B01)the Qing-Lan Project for Cloud-Fog-Precipitation-Aerosol Study in Jiangsu Province, the Graduate Student Innovation Plan for the Universities of Jiangsu Province (Grant No CX10B 292Z)a Project Funded by the Priority Academic Development of Jiangsu Higher Education Institutions
文摘Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the MoninObukhov similarity for describing the turbulent fluctu ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * σ v /u * σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day clear sky overcast cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "1/3" slope in unstable stratification and "2/3" slope in stable stratification.
基金Project supported by the National Natural Science Foundation of China(Nos.11672351 and11332007)the National Key R&D Plan(No.2016YFA0401200)the FengLei Youth Innovation Fund of China Aerodynamics Research and Development Center(No.KT-FLJJ-201803)
文摘The elementary task is to calculate the growth rates of disturbances when the e;method in transition prediction is performed. However, there is no unified knowledge to determine the growth rates of disturbances in three-dimensional(3 D) flows. In this paper, we study the relation among the wave parameters of the disturbance in boundary layers in which the imaginary parts of wave parameters are far smaller than the real parts.The generalized growth rate(GGR) in the direction of group velocity is introduced, and the conservation relation of GGR is strictly deduced in theory. This conservation relation manifests that the GGR only depends on the real parts of wave parameters instead of the imaginary parts. Numerical validations for GGR conservation are also provided in the cases of first/second modes and crossflow modes. The application of GGR to the eN method in 3 D flows is discussed, and the puzzle of determining growth rates in 3 D flows is clarified. A convenient method is also proposed to calculate growth rates of disturbances in 3 D flows. Good agreement between this convenient method and existing methods is found except the condition that the angle between the group velocity direction and the x-direction is close to 90?which can be easily avoided in practical application.
文摘The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, which is very effective, was used to study the nonlinear evolution. of T-S disturbance wave at finite amplitudes. Solving nonlinear equations of different modes by using predictor-corrector and iterative approach, which is uncoupled between modes, improving computational accuracy by using high order compact differential scheme, satisfying normalization condition I determining tables of nonlinear terms at different modes, and implementing stably the spatial marching, were included in this method. With different initial amplitudes, the nonlinear evolution of T-S wave was studied. The nonlinear nonparallel results of examples compare with data of direct numerical simulations (DNS) using full Navier-Stokes equations.
基金supported by the National Key Research and Development Program of China(No.2021YFA0717100)the National Natural Science Foundation of China(Nos.12072270 and U2013206).
文摘TheMoon is the only celestial body that human beings have visited.The design of the Earth-Moon transfer orbits is a critical issue in lunar exploration missions.In the 21st century,new lunar missions including the construction of the lunar space station,the permanent lunar base,and the Earth-Moon transportation network have been proposed,requiring low-cost,expansive launch windows and a fixed arrival epoch for any launch date within the launch window.The low-energy and low-thrust transfers are promising strategies to satisfy the demands.This review provides a detailed landscape of Earth-Moon transfer trajectory design processes,from the traditional patched conic to the state-of-the-art low-energy and low-thrust methods.Essential mechanisms of the various utilized dynamic models and the characteristics of the different design methods are discussed in hopes of helping readers grasp thebasic overviewof the current Earth-Moon transfer orbitdesignmethods anda deep academic background is unnecessary for the context understanding.
基金supported by the National Natural Science Foundation of China under Project 52007133 and U22B20100。
文摘This paper proposes a novel Multivariate Quotient-Difference(MQD)method to obtain the approximate analytical solution for AC power flow equations.Therefore,in the online environment,the power flow solutions covering different operating conditions can be directly obtained by plugging values into multiple symbolic variables,such that the power injections and consumptions of selected buses or areas can be independently adjusted.This method first derives a power flow solution through a Multivariate Power Series(MPS).Next,the MQD method is applied to transform the obtained MPS to a Multivariate Pad´e Approximants(MPA)to expand the Radius of Convergence(ROC),so that the accuracy of the derived analytical solution can be significantly increased.In addition,the hypersurface of the voltage stability boundary can be identified by an analytical formula obtained from the coefficients of MPA.This direct method for power flow solutions and voltage stability boundaries is fast for many online applications,since such analytical solutions can be derived offline and evaluated online by only plugging values into the symbolic variables according to the actual operating conditions.The proposed method is validated in detail on New England 39-bus and IEEE 118-bus systems with independent load variations in multi-regions.
基金Supported in part by the National Natural Science Foundation of China under Grant 52307230the Natural Science Foundation of Shandong Province under Grant ZR2023QE156.
文摘Droop-controlled voltage-source converters(VSCs)can provide frequency and voltage support to power grids.However,during a grid fault,VsCs may experience transient instability,which can be significantly affected by both the control parameters and fault conditions.This mechanism has not been fully elucidated in previous studies.In particular,grid-voltage faults are commonly accompanied by a grid voltage phase-angle jump(VPAJ),which is typically ignored in the evaluation of the transient stability of VSCs.To address this issue,this study comprehensively assesses the impact of the VPAJ and key control parameters on the transient characteristics of VSCs.Furthermore,the critical clearing angle and critical clearing time are quantitatively calculated to define the transient stability boundary.In addition,a transient stability-enhancement control method that considers the transient stability constraints is proposed.Finally,simulations and experimental tests are conducted to validate both the theoretical analysis and proposed method.
基金supported by National Key Basic Research Program (973 Program, Grant No. 2011CB706804)National Natural Science Foundation of China (Grant No. 50805093)Science & Technology Commission of Shanghai Municipality, China (Grant No. 09QH1401500)
文摘As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.
基金National Natural Science Foundation of China (10772082)Doctoral Foundation of Ministry of Education of China (20070287005)
文摘This article studies the nonlinear evolution of disturbance waves in supersonic nonparallel boundary layer flows by using nonlinear parabolic stability equations (NPSE). An accurate numerical method is developed to solve the equations and march the NPSE in a stable manner. Through computation,are obtained the curves of amplitude and disturbance shape function of harmonic waves. Especially are demonstrated the physical characteristics of nonlinear stability of various harmonic waves,including instantaneous stream wise vortices,spanwise vortices and Λ structure etc,and are used to study and analyze the mechanism of the transition process. The calculated results have evidenced the effectiveness of the proposed NPSE method to research the nonlinear stability of the supersonic boundary layers.
文摘Subject Code:E01 Conventional polycrystalline materials become harder with decreasing grain size,following the classical Hall—Petch relationship,i.e.,strength increases reversely proportional to the square root of the grain size.Strengthening occurs due to dislocation pileups at grain boundaries that prevent the dislocations
文摘The nonparallel effects on the stability of the boundary layer flow was investigated using the Parabolie Stability Equations (PSE). In order to improve the accuracy of the calculation which is very important for the investigation of stability, higher order expansions in orthogonal functions in normal direction and the effective algebraic mapping to deal with the problem of infinite region were used and the way to collocate the boundary point based on its characteristics was adopted. With the effective control of step size in the marching procedure, the special condition was satisfied, and the stability of calculation was assured. From the curves of the neutral stability, the growth rate, the amplitude variation and disturbed velocity profile, the effects of the nonparallelism were given accurately and analyzed detailedly. It is found that the nonparallelism of the flow amplifies the amplitude and growth rate of disturbances, especially for three-dimensional disturbances, even can change the sign of flow stability from stability to instability for some cases. Computed results are in good agreement with the classical experimental results.
基金supported by the National Natural Science Foundation of China (No. 51476152)
文摘The linear instabilities of incompressible confluent mixing layer and boundary layer were analyzed.The mixing layers include wake,shear layer and their combination.The mean velocity profile of confluent flow is taken as a superposition of a hyperbolic and exponential function to model a mixing layer and the Blasius similarity solution for a flat plate boundary layer.The stability equation of confluent flow was solved by using the global numerical method.The unstable modes associated with both the mixing and boundary layers were identified.They are the boundary layer mode,mixing layer mode 1(nearly symmetrical mode)and mode 2(nearly anti-symmetrical mode).The interactions between the mixing layer stability and the boundary layer stability were examined.As the mixing layer approaches the boundary layer,the neutral curves of the boundary layer mode move to the upper left,the resulting critical Reynolds number decreases,and the growth rate of the most unstable mode increases.The wall tends to stabilize the mixing layer modes at low frequency.In addition,the mode switching behavior of the relative level of the spatial growth rate between the mixing layer mode 1 and mode 2 with the velocity ratio is found to occur at low frequency.
基金supported by National Natural Science Foundation of China (Grant Nos. 11271214 and 11671224)Ben Andrews’ ARC Laureate Fellowship (Grant No. FL150100126)
文摘We study the stability of unduloids with free boundary in the domain B between two parallel hyperplanes in R^(n+1). If the unduloid has one half of period in B and is sufficiently close to a cylinder, then for 2≤n≤10, it is unstable; while for n≥11, it is stable. If the unduloid has two or more halves of period in B and is sufficiently close to a cylinder, then for all n≥2, it is unstable.
基金Grant for Supporting Plan of Tsinghua University.
文摘In this paper, we consider the boundary stabilization of the wave equation with variable coefficients by Riemmannian geometry method subject to a different geometric condition which is motivated by the geometric multiplier identities. Several (multiplier) identities (inequalities) which have been built for constant wave equation by Kormornik and Zuazua are generalized to the variable coefficient case by some computational techniques in Riemmannian geometry, so that the precise estimates on the exponential decay rate are derived from those inequalitities. Also, the exponential decay for the solutions of semilinear wave equation with variable coefficients is obtained under natural growth and sign assumptions on the nonlinearity. Our method is rather general and can be adapted to other evolution systems with variable coefficients (e.g. elasticity plates) as well.
基金The authors acknowledge the support of National Natural Science Foundation of China(Nos.11672187,11572182)Natural Science Foundation of Liaoning Province(201602573)+2 种基金the Key Research Projects of Shanghai Science and Technology Commission(No.18010500100)Innovation Program of Shanghai Education Commission(No.2017-01-07-00-09-E00019)Beiyang Young Scholars of Tian jin University(2019XRX-0027).
文摘In this study,the transverse vibration of a traveling beam made of functionally graded material was analyzed.The material gradation was assumed to vary continuously along the thickness direction of the beam in the form of power law exponent.The effect of the longitudinally varying tension due to axial acceleration was highlighted,and the dependence of the tension on the finite support rigidity was also considered.A complex governing equation of the functionally graded beam was derived by the Hamilton principle,in which the geometric nonlinearity,material properties and axial load were incorporated.The direct multiscale method was applied to the analysis process of an axially moving functionally graded beam with timedependent velocity,and the natural frequency and solvability conditions were obtained.Based on the conditions,the stability boundaries of subharmonic resonance and combination resonance were obtained.It was found that the dynamic behavior of axial moving beams could be tuned by using the distribution law of the functional gradient parameters.