Two dimensional parabolic stability equations (PSE) are numerically solved using expansions in orthogonal functions in the normal direction.The Chebyshev polynomials approximation,which is a very useful form of ortho...Two dimensional parabolic stability equations (PSE) are numerically solved using expansions in orthogonal functions in the normal direction.The Chebyshev polynomials approximation,which is a very useful form of orthogonal expansions, is applied to solving parabolic stability equations. It is shown that results of great accuracy are effectively obtained.The availability of using Chebyshev approximations in parabolic stability equations is confirmed.展开更多
Parabolized stability equations （PSE） were used to study the evolution of disturbances in compressible boundary layers. The results were compared with those obtained by direct numerical simulations （DNS）, to check...Parabolized stability equations （PSE） were used to study the evolution of disturbances in compressible boundary layers. The results were compared with those obtained by direct numerical simulations （DNS）, to check if the results from PSE method were reliable or not. The results of comparison showed that no matter for subsonic or supersonic boundary layers, results from both the PSE and DNS method agreed with each other reasonably well, and the agreement between temperatures was better than those between velocities. In addition, linear PSE was used to calculate the neutral curve for small amplitude disturbances in a supersonic boundary layer. Compared with those obtained by linear stability theory （LST）, the situation was similar to those for incom- pressible boundary layer.展开更多
The parabolized stability equation （PSE） method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulat...The parabolized stability equation （PSE） method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation （DNS） as applied to a problem of compressible boundary layer with Mach number 6.展开更多
By using characteristic analysis of the linear and nonlinear parabolic stability equations ( PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub- characteristic analysis of PS...By using characteristic analysis of the linear and nonlinear parabolic stability equations ( PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub- characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic., respectively . The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories , the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time , the methods of removing the remained ellipticity are further obtained from the nonlinear PSE .展开更多
An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber ...An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.展开更多
The nonlinear parabolized stability equations(NPSEs)approach is widely used to study the evolution of disturbances in hypersonic boundary layers owing to its high computational efficiency.However,divergence of the NPS...The nonlinear parabolized stability equations(NPSEs)approach is widely used to study the evolution of disturbances in hypersonic boundary layers owing to its high computational efficiency.However,divergence of the NPSEs will occur when disturbances imposed at the inlet no longer play a leading role or when the nonlinear effect becomes very strong.Two major improvements are proposed here to deal with the divergence of the NPSEs.First,all disturbances are divided into two types:dominant waves and non-dominant waves.Disturbances imposed at the inlet or playing a leading role are defined as dominant waves,with all others being defined as non-dominant waves.Second,the streamwise wavenumbers of the non-dominant waves are obtained using the phase-locked method,while those of the dominant waves are obtained using an iterative method.Two reference wavenumbers are introduced in the phase-locked method,and methods for calculating them for different numbers of dominant waves are discussed.Direct numerical simulation(DNS)is performed to verify and validate the predictions of the improved NPSEs in a hypersonic boundary layer on an isothermal swept blunt plate.The results from the improved NPSEs approach are in good agreement with those of DNS,whereas the traditional NPSEs approach is subject to divergence,indicating that the improved NPSEs approach exhibits greater robustness.展开更多
The enhanced mountain-to-plain convective storms in Beijing on 22 May 2021 were simulated using the highresolution Weather Research and Forecasting model,enabling detailed analyses of convective instability characteri...The enhanced mountain-to-plain convective storms in Beijing on 22 May 2021 were simulated using the highresolution Weather Research and Forecasting model,enabling detailed analyses of convective instability characteristics and underlying causes of stability variations.Generalized potential temperature outperformed traditional potential temperature and equivalent potential temperature in capturing instability variations associated with mid-level latent heating and near-surface evaporative cooling.Local instability variance was primarily governed by potential divergence and the advection of potential instability,with these factors exhibiting out-of-phase distributions.Prior to the onset of heavy precipitation,intense downdrafts transported unstable air from higher levels into more stable regions at lower levels,increasing local near-surface instability,which contributed to the formation of heavy precipitation.During the heavy precipitation stage,vertical divergence between slantwise updrafts and downdrafts in the lowmiddle stable layers led to destabilization,supporting sustained convective development within the precipitation area.At the leading edge of the heavy precipitation,instability enhancement was primarily driven by vertical advection,and less stable air in the lower levels was transported upward,enhancing instability at higher levels.展开更多
The paper deals with the criteria for the closed- loop stability of a noise control system in a duct. To study the stability of the system, the model of delay differential equation is derived from the propagation of a...The paper deals with the criteria for the closed- loop stability of a noise control system in a duct. To study the stability of the system, the model of delay differential equation is derived from the propagation of acoustic wave governed by a partial differential equation of hyperbolic type. Then, a simple feedback controller is designed, and its closed- loop stability is analyzed on the basis of the derived model of delay differential equation. The obtained criteria reveal the influence of the controller gain and the positions of a sensor and an actuator on the closed-loop stability. Finally, numerical simulations are presented to support the theoretical results.展开更多
In Ref [1] the asymptotic stability of nonlinear slowly changing system has been discussed .In Ref [2] the instability of solution for the order linear differential equaiton with varied coefficient has been discussed ...In Ref [1] the asymptotic stability of nonlinear slowly changing system has been discussed .In Ref [2] the instability of solution for the order linear differential equaiton with varied coefficient has been discussed .In this paper,we have discussed instability of solution for a class of the third order nonlinear diffeential equation by means of the metod of Refs [1] and [2] .展开更多
We introduce the Dirac equation in four-dimensional gravity which is a generally covariant form. We choose the suitable variable and solve the corresponding equation. To solve such equation and to obtain the correspon...We introduce the Dirac equation in four-dimensional gravity which is a generally covariant form. We choose the suitable variable and solve the corresponding equation. To solve such equation and to obtain the corresponding bispinor, we employ the factorization method which introduces the associated Laguerre polynomial. The asso- ciated Laguerre polynomials help us to write the Dirac equation of four-dimensional gravity in the form of the shape invariance equation. Thus we write the shape invariance condition with respect to the secondary quantum number. Finally, we obtain the spinor wave function and achieve the corresponding stability of condition for the four-dimensional gravity system.展开更多
On the stability analysis of large-scale systems by Lyapunov functions, it is necessary to determine the stability of vector comparison equations. For discrete systems, only the stability of linear autonomous comparis...On the stability analysis of large-scale systems by Lyapunov functions, it is necessary to determine the stability of vector comparison equations. For discrete systems, only the stability of linear autonomous comparison equations was studied in the past. In this paper, various criteria of stability for discrete nonlinear autonomous comparison equations are completely established. Among them, a criterion for asymptotic stability is not only sufficient, but also necessary, from which a criterion on the function class C, is derived. Both of them can be used to determine the unexponential stability, even in the large, for discrete nonlinear (autonomous or nonautonomous) systems. All the criteria are of simple algebraic forms and can be readily used.展开更多
In this paper we are interested in the large time behavior of the nonlinear diffusion equationWe consider functions which allow the equation to possess traveling wave solutions. We first present an existence and uniqu...In this paper we are interested in the large time behavior of the nonlinear diffusion equationWe consider functions which allow the equation to possess traveling wave solutions. We first present an existence and uniqueness as well as some comparison principle result of generalized solutions to the Cauchy problem. Then we give for some threshold results, from which we can see that u=a is stable, while u= 0 or u=1 is unstable under some assumptions, etc.展开更多
The parabolized stability equations （PSEs） for high speed flows, especially supersonic and hypersonic flows, are derived and used to analyze the nonparallel boundary layer stability. The proposed numerical technique...The parabolized stability equations （PSEs） for high speed flows, especially supersonic and hypersonic flows, are derived and used to analyze the nonparallel boundary layer stability. The proposed numerical techniques for solving PSE include the following contents： introducing the efficiently normal transformation of the boundary layer, improving the computational accuracy by using a high-order differential scheme near the wall, employing the predictor-corrector and iterative approach to satisfy the important normalization condition, and implementing the stable spatial marching. Since the second mode dominates the growth of the disturbance in high Mach number flows, it is used in the computation. The evolution and characteristics of the boundary layer stability in the high speed flow are demonstrated in the examples. The effects of the nonparallelizm, the compressibility and the cooling wall on the stability are analyzed. And computational results are in good agreement with the relevant data.展开更多
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.展开更多
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.展开更多
The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory （LST） based on the local parallel hypothesis. Considering t...The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory （LST） based on the local parallel hypothesis. Considering the non-parallelism effect, the parabolized stability equation （PSE） method lacks local characteristic of stability analysis. In this paper, a local stability analysis method considering non-parallelism is proposed, termed as EPSE since it may be considered as an expansion of the PSE method. The EPSE considers variation of the shape function in the streamwise direction. Its local characteristic is convenient for stability analysis. This paper uses the EPSE in a strong non-parallel flow and mode exchange problem. The results agree well with the PSE and the direct numerical simulation （DNS）. In addition, it is found that the growth rate is related to the normalized method in the non-parallel flow. Different results can be obtained using different normalized methods. Therefore, the normalized method must be consistent.展开更多
It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, es...It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation （3D- LPSE） approach, a 3D extension of the two-dimensional LPSE （2D-LPSE）, is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation （DNS） rather than a 2D-eigenvalue problem （EVP） procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.展开更多
This study is to numerically test the interfacial instability of ferrofluid flow under the presence of a vacuum magnetic field.The ferrofluid parabolized stability equations(PSEs)are derived from the ferrofluid stabil...This study is to numerically test the interfacial instability of ferrofluid flow under the presence of a vacuum magnetic field.The ferrofluid parabolized stability equations(PSEs)are derived from the ferrofluid stability equations and the Rosensweig equations,and the characteristic values of the ferrofluid PSEs are given to describe the ellipticity of ferrofluid flow.Three numerical models representing specific cases considering with/without a vacuum magnetic field or viscosity are created to mathematically examine the interfacial instability by the computation of characteristic values.Numerical investigation shows strong dependence of the basic characteristic of ferrofluid Rayleigh-Taylor instability(RTI)on viscosity of ferrofluid and independence of the vacuum magnetic field.For the shock wave striking helium bubble,the magnetic field is not able to trigger the symmetry breaking of bubble but change the speed of the bubble movement.In the process of droplet formation from a submerged orifice,the collision between the droplet and the liquid surface causes symmetry breaking.Both the viscosity and the magnetic field exacerbate symmetry breaking.The computational results agree with the published experimental results.展开更多
The nth-order expansion of the parabolized stability equation （EPSEn） is obtained from the Taylor expansion of the linear parabolized stability equation （LPSE） in the streamwise direction. The EPSE together with t...The nth-order expansion of the parabolized stability equation （EPSEn） is obtained from the Taylor expansion of the linear parabolized stability equation （LPSE） in the streamwise direction. The EPSE together with the homogeneous boundary conditions forms a local eigenvalue problem, in which the streamwise variations of the mean flow and the disturbance shape function are considered. The first-order EPSE （EPSE1） and the second-order EPSE （EPSE2） are used to study the crossflow instability in the swept NLF（2）-0415 wing boundary layer. The non-parallelism degree of the boundary layer is strong. Compared with the growth rates predicted by the linear stability theory （LST）, the results given by the EPSE1 and EPSE2 agree well with those given by the LPSE. In particular, the results given by the EPSE2 are almost the same as those given by the LPSE. The prediction of the EPSE1 is more accurate than the prediction of the LST, and is more efficient than the predictions of the EPSE2 and LPSE. Therefore, the EPSE1 is an efficient ey prediction tool for the crossflow instability in swept-wing boundary-layer flows.展开更多
Parabolized stability equations （PSE） approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether ...Parabolized stability equations （PSE） approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.展开更多
文摘Two dimensional parabolic stability equations (PSE) are numerically solved using expansions in orthogonal functions in the normal direction.The Chebyshev polynomials approximation,which is a very useful form of orthogonal expansions, is applied to solving parabolic stability equations. It is shown that results of great accuracy are effectively obtained.The availability of using Chebyshev approximations in parabolic stability equations is confirmed.
基金Project supported by the National Natural Science Foundation of China (Key Program)(No.10632050)the Science Foundation of Liuhui Center of Applied Mathematics,Nankai University and Tianjin University.
文摘Parabolized stability equations （PSE） were used to study the evolution of disturbances in compressible boundary layers. The results were compared with those obtained by direct numerical simulations （DNS）, to check if the results from PSE method were reliable or not. The results of comparison showed that no matter for subsonic or supersonic boundary layers, results from both the PSE and DNS method agreed with each other reasonably well, and the agreement between temperatures was better than those between velocities. In addition, linear PSE was used to calculate the neutral curve for small amplitude disturbances in a supersonic boundary layer. Compared with those obtained by linear stability theory （LST）, the situation was similar to those for incom- pressible boundary layer.
基金supported by the National Natural Science Foundation of China(Nos.11202147,11472188,11332007,11172203,and 91216111)the Specialized Research Fund(New Teacher Class)for the Doctoral Program of Higher Education(No.20120032120007)
文摘The parabolized stability equation （PSE） method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation （DNS） as applied to a problem of compressible boundary layer with Mach number 6.
基金the National Natural Science Foundation of China (10032050)the National 863 Program Foundation of China (2002AA633100)
文摘By using characteristic analysis of the linear and nonlinear parabolic stability equations ( PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub- characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic., respectively . The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories , the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time , the methods of removing the remained ellipticity are further obtained from the nonlinear PSE .
基金Project supported by the National Natural Science Foundation of China(Nos.11332007,11402167,11672205,and 11732011)the National Key Research and Development Program of China(No.2016YFA0401200)
文摘An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.
基金the National Natural Science Foundation of China(Grant Nos.12072232 and 11672351)the National Key Project of China(Grant No.GJXM92579).
文摘The nonlinear parabolized stability equations(NPSEs)approach is widely used to study the evolution of disturbances in hypersonic boundary layers owing to its high computational efficiency.However,divergence of the NPSEs will occur when disturbances imposed at the inlet no longer play a leading role or when the nonlinear effect becomes very strong.Two major improvements are proposed here to deal with the divergence of the NPSEs.First,all disturbances are divided into two types:dominant waves and non-dominant waves.Disturbances imposed at the inlet or playing a leading role are defined as dominant waves,with all others being defined as non-dominant waves.Second,the streamwise wavenumbers of the non-dominant waves are obtained using the phase-locked method,while those of the dominant waves are obtained using an iterative method.Two reference wavenumbers are introduced in the phase-locked method,and methods for calculating them for different numbers of dominant waves are discussed.Direct numerical simulation(DNS)is performed to verify and validate the predictions of the improved NPSEs in a hypersonic boundary layer on an isothermal swept blunt plate.The results from the improved NPSEs approach are in good agreement with those of DNS,whereas the traditional NPSEs approach is subject to divergence,indicating that the improved NPSEs approach exhibits greater robustness.
基金funded by the Beijing Municipal Science and Technology Commission [grant number Z221100005222012]the Department of Science and Technology of Hebei Province [grant number 22375404D]+2 种基金the Strategic Priority Research Program of the Chinese Academy of Sciences [grant number XDB0760303]the National Natural Science Foundation of China [grant numbers U2233218 and 42275010]the Open Foundation of the Key Open Laboratory of Urban Meteorology,China Meteorological Administration [grant number LUM-2023-06]。
文摘The enhanced mountain-to-plain convective storms in Beijing on 22 May 2021 were simulated using the highresolution Weather Research and Forecasting model,enabling detailed analyses of convective instability characteristics and underlying causes of stability variations.Generalized potential temperature outperformed traditional potential temperature and equivalent potential temperature in capturing instability variations associated with mid-level latent heating and near-surface evaporative cooling.Local instability variance was primarily governed by potential divergence and the advection of potential instability,with these factors exhibiting out-of-phase distributions.Prior to the onset of heavy precipitation,intense downdrafts transported unstable air from higher levels into more stable regions at lower levels,increasing local near-surface instability,which contributed to the formation of heavy precipitation.During the heavy precipitation stage,vertical divergence between slantwise updrafts and downdrafts in the lowmiddle stable layers led to destabilization,supporting sustained convective development within the precipitation area.At the leading edge of the heavy precipitation,instability enhancement was primarily driven by vertical advection,and less stable air in the lower levels was transported upward,enhancing instability at higher levels.
基金the National Natural Science Foundation of China (10532050)
文摘The paper deals with the criteria for the closed- loop stability of a noise control system in a duct. To study the stability of the system, the model of delay differential equation is derived from the propagation of acoustic wave governed by a partial differential equation of hyperbolic type. Then, a simple feedback controller is designed, and its closed- loop stability is analyzed on the basis of the derived model of delay differential equation. The obtained criteria reveal the influence of the controller gain and the positions of a sensor and an actuator on the closed-loop stability. Finally, numerical simulations are presented to support the theoretical results.
文摘In Ref [1] the asymptotic stability of nonlinear slowly changing system has been discussed .In Ref [2] the instability of solution for the order linear differential equaiton with varied coefficient has been discussed .In this paper,we have discussed instability of solution for a class of the third order nonlinear diffeential equation by means of the metod of Refs [1] and [2] .
文摘We introduce the Dirac equation in four-dimensional gravity which is a generally covariant form. We choose the suitable variable and solve the corresponding equation. To solve such equation and to obtain the corresponding bispinor, we employ the factorization method which introduces the associated Laguerre polynomial. The asso- ciated Laguerre polynomials help us to write the Dirac equation of four-dimensional gravity in the form of the shape invariance equation. Thus we write the shape invariance condition with respect to the secondary quantum number. Finally, we obtain the spinor wave function and achieve the corresponding stability of condition for the four-dimensional gravity system.
文摘On the stability analysis of large-scale systems by Lyapunov functions, it is necessary to determine the stability of vector comparison equations. For discrete systems, only the stability of linear autonomous comparison equations was studied in the past. In this paper, various criteria of stability for discrete nonlinear autonomous comparison equations are completely established. Among them, a criterion for asymptotic stability is not only sufficient, but also necessary, from which a criterion on the function class C, is derived. Both of them can be used to determine the unexponential stability, even in the large, for discrete nonlinear (autonomous or nonautonomous) systems. All the criteria are of simple algebraic forms and can be readily used.
文摘In this paper we are interested in the large time behavior of the nonlinear diffusion equationWe consider functions which allow the equation to possess traveling wave solutions. We first present an existence and uniqueness as well as some comparison principle result of generalized solutions to the Cauchy problem. Then we give for some threshold results, from which we can see that u=a is stable, while u= 0 or u=1 is unstable under some assumptions, etc.
文摘The parabolized stability equations （PSEs） for high speed flows, especially supersonic and hypersonic flows, are derived and used to analyze the nonparallel boundary layer stability. The proposed numerical techniques for solving PSE include the following contents： introducing the efficiently normal transformation of the boundary layer, improving the computational accuracy by using a high-order differential scheme near the wall, employing the predictor-corrector and iterative approach to satisfy the important normalization condition, and implementing the stable spatial marching. Since the second mode dominates the growth of the disturbance in high Mach number flows, it is used in the computation. The evolution and characteristics of the boundary layer stability in the high speed flow are demonstrated in the examples. The effects of the nonparallelizm, the compressibility and the cooling wall on the stability are analyzed. And computational results are in good agreement with the relevant data.
基金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.
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.11332007,11172203,and 91216111)
文摘The e-N method is widely used in transition prediction. The amplitude growth rate used in the e-N method is usually provided by the linear stability theory （LST） based on the local parallel hypothesis. Considering the non-parallelism effect, the parabolized stability equation （PSE） method lacks local characteristic of stability analysis. In this paper, a local stability analysis method considering non-parallelism is proposed, termed as EPSE since it may be considered as an expansion of the PSE method. The EPSE considers variation of the shape function in the streamwise direction. Its local characteristic is convenient for stability analysis. This paper uses the EPSE in a strong non-parallel flow and mode exchange problem. The results agree well with the PSE and the direct numerical simulation （DNS）. In addition, it is found that the growth rate is related to the normalized method in the non-parallel flow. Different results can be obtained using different normalized methods. Therefore, the normalized method must be consistent.
基金Project supported by the National Natural Science Foundation of China(Nos.11272183,11572176,11402167,11202147,and 11332007)the National Program on Key Basic Research Project of China(No.2014CB744801)
文摘It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation （3D- LPSE） approach, a 3D extension of the two-dimensional LPSE （2D-LPSE）, is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation （DNS） rather than a 2D-eigenvalue problem （EVP） procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.
基金the National Natural Science Foundation of China(No.11971411)the Research Foundation of Education Bureau of Hunan Province of China(No.18A067)。
文摘This study is to numerically test the interfacial instability of ferrofluid flow under the presence of a vacuum magnetic field.The ferrofluid parabolized stability equations(PSEs)are derived from the ferrofluid stability equations and the Rosensweig equations,and the characteristic values of the ferrofluid PSEs are given to describe the ellipticity of ferrofluid flow.Three numerical models representing specific cases considering with/without a vacuum magnetic field or viscosity are created to mathematically examine the interfacial instability by the computation of characteristic values.Numerical investigation shows strong dependence of the basic characteristic of ferrofluid Rayleigh-Taylor instability(RTI)on viscosity of ferrofluid and independence of the vacuum magnetic field.For the shock wave striking helium bubble,the magnetic field is not able to trigger the symmetry breaking of bubble but change the speed of the bubble movement.In the process of droplet formation from a submerged orifice,the collision between the droplet and the liquid surface causes symmetry breaking.Both the viscosity and the magnetic field exacerbate symmetry breaking.The computational results agree with the published experimental results.
基金supported by the National Natural Science Foundation of China(No.11332007)
文摘The nth-order expansion of the parabolized stability equation （EPSEn） is obtained from the Taylor expansion of the linear parabolized stability equation （LPSE） in the streamwise direction. The EPSE together with the homogeneous boundary conditions forms a local eigenvalue problem, in which the streamwise variations of the mean flow and the disturbance shape function are considered. The first-order EPSE （EPSE1） and the second-order EPSE （EPSE2） are used to study the crossflow instability in the swept NLF（2）-0415 wing boundary layer. The non-parallelism degree of the boundary layer is strong. Compared with the growth rates predicted by the linear stability theory （LST）, the results given by the EPSE1 and EPSE2 agree well with those given by the LPSE. In particular, the results given by the EPSE2 are almost the same as those given by the LPSE. The prediction of the EPSE1 is more accurate than the prediction of the LST, and is more efficient than the predictions of the EPSE2 and LPSE. Therefore, the EPSE1 is an efficient ey prediction tool for the crossflow instability in swept-wing boundary-layer flows.
基金Project supported by the National Natural Science Foundation of China(Nos.10632050,90716007)the Foundation of LIU Hui Center of Applied Mathematics of Nankai University and Tianjin University
文摘Parabolized stability equations （PSE） approach is used to investigate problems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.
