An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approxi...An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.展开更多
Based on the work of paper [1], we propose a modified Levenberg-Marquardt algoithm for solving singular system of nonlinear equations F(x) = 0, where F(x) : Rn - Rn is continuously differentiable and F'(x) is Lips...Based on the work of paper [1], we propose a modified Levenberg-Marquardt algoithm for solving singular system of nonlinear equations F(x) = 0, where F(x) : Rn - Rn is continuously differentiable and F'(x) is Lipschitz continuous. The algorithm is equivalent to a trust region algorithm in some sense, and the global convergence result is given. The sequence generated by the algorithm converges to the solution quadratically, if ||F(x)||2 provides a local error bound for the system of nonlinear equations. Numerical results show that the algorithm performs well.展开更多
The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equatio...The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical results confirm the theoretical prediction of the exponential rate of convergence. The result in this work seems to be the first successful spectral approach (with theoretical justification) for the Volterra type equations.展开更多
This paper deals with nonlinear expectations. The author obtains a nonlinear gen- eralization of the well-known Kolmogorov’s consistent theorem and then use it to con- struct ?ltration-consistent nonlinear expectatio...This paper deals with nonlinear expectations. The author obtains a nonlinear gen- eralization of the well-known Kolmogorov’s consistent theorem and then use it to con- struct ?ltration-consistent nonlinear expectations via nonlinear Markov chains. Com- pared to the author’s previous results, i.e., the theory of g-expectations introduced via BSDE on a probability space, the present framework is not based on a given probabil- ity measure. Many fully nonlinear and singular situations are covered. The induced topology is a natural generalization of Lp-norms and L∞-norm in linear situations. The author also obtains the existence and uniqueness result of BSDE under this new framework and develops a nonlinear type of von Neumann-Morgenstern representation theorem to utilities and present dynamic risk measures.展开更多
High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of ...High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of mu(= h/lambda, depth to deep-water wave length ratio) and epsilon(= a/h, wave amplitude to depth ratio) for velocity potential, particle velocity vector, pressure and the Boussinesq-type equations for surface elevation eta and horizontal velocity vector (U) over right arrow at any given level in water are given. Then, the exact explicit expressions to the fourth order of mu are derived. Finally, the linear solutions of eta, (U) over right arrow, C (phase-celerity) and C-g (group velocity) for a constant water depth are obtained. Compared with the Airy theory, excellent results can be found even for a water depth as large as the wave legnth. The present high-order models are applicable to nonlinear regular and irregular waves in water of any varying depth (from shallow to deep) and bottom slope (from mild to steep).展开更多
It is a comparatively convenient technique to investigate the motion of a particle with the help of the differential geometry the-ory,rather than directly decomposing the motion in the Cartesian coordinates.The new mo...It is a comparatively convenient technique to investigate the motion of a particle with the help of the differential geometry the-ory,rather than directly decomposing the motion in the Cartesian coordinates.The new model of three-dimensional (3D) guidance problem for interceptors is presented in this paper,based on the classical differential geometry curve theory.Firstly,the kinematical equations of the line of sight (LOS) are gained by carefully investigating the rotation principle of LOS,the kinematic equations of LOS are established,and the concepts of curvature and torsion of LOS are proposed.Simultaneously,the new relative dynamic equations between interceptor and target are constructed.Secondly,it is found that there is an instan-taneous rotation plane of LOS (IRPL) in the space,in which two-dimensional (2D) guidance laws could be constructed to solve 3D interception guidance problems.The spatial 3D true proportional navigation (TPN) guidance law could be directly introduced in IRPL without approximation and linearization for dimension-reduced 2D TPN.In addition,the new series of augmented TPN (APN) and LOS angular acceleration guidance laws (AAG) could also be gained in IRPL.After that,the dif-ferential geometric guidance commands (DGGC) of guidance laws in IRPL are advanced,and we prove that the guidance commands in arc-length system proposed by Chiou and Kuo are just a special case of DGGC.Moreover,the performance of the original guidance laws will be reduced after the differential geometric transformation.At last,an exoatmospheric intercep-tion is taken for simulation to demonstrate the differential geometric modeling proposed in this paper.展开更多
A nonconforming H^1-Calerkin mixed finite element method is analyzed for Sobolev equations on anisotropic meshes. The error estimates are obtained without using Ritz-Volterra projection.
The equations of motion of an insect with flapping wings are derived and then simplified to that of a flying body using the "rigid body" assumption. On the basis of the simplified equations of motion, the longitudin...The equations of motion of an insect with flapping wings are derived and then simplified to that of a flying body using the "rigid body" assumption. On the basis of the simplified equations of motion, the longitudinal dynamic flight stability of four insects (hoverfly, cranefly, dronefly and hawkmoth) in hovering flight is studied (the mass of the insects ranging from 11 to 1,648 mg and wingbeat frequency from 26 to 157Hz). The method of computational fluid dynamics is used to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis are used to solve the equations of motion. The validity of the "rigid body" assumption is tested and how differences in size and wing kinematics influence the applicability of the "rigid body" assumption is investigated. The primary findings are: (1) For insects considered in the present study and those with relatively high wingbeat frequency (hoverfly, drone fly and bumblebee), the "rigid body" assumption is reasonable, and for those with relatively low wingbeat frequency (cranefly and howkmoth), the applicability of the "rigid body" assumption is questionable. (2) The same three natural modes of motion as those reported recently for a bumblebee are identified, i.e., one unstable oscillatory mode, one stable fast subsidence mode and one stable slow subsidence mode. (3) Approximate analytical expressions of the eigenvalues, which give physical insight into the genesis of the natural modes of motion, are derived. The expressions identify the speed derivative Mu (pitching moment produced by unit horizontal speed) as the primary source of the unstable oscillatory mode and the stable fast subsidence mode and Zw (vertical force produced by unit vertical speed) as the primary source of the stable slow subsidence mode.展开更多
We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement ...We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement learning with the gradient of the solution playing the role of the policy function,and the loss function given by the error between the prescribed terminal condition and the solution of the BSDE.The policy function is then approximated by a neural network,as is done in deep reinforcement learning.Numerical results using TensorFlow illustrate the efficiency and accuracy of the studied algorithm for several 100-dimensional nonlinear PDEs from physics and finance such as the Allen–Cahn equation,the Hamilton–Jacobi–Bellman equation,and a nonlinear pricing model for financial derivatives.展开更多
This paper gives several fundamental theorems for the stability, uniform stability, asymptotic stability and uniform asymptotic stability. Those theorems allow the derivative of Lyapunov functions to be positive on ce...This paper gives several fundamental theorems for the stability, uniform stability, asymptotic stability and uniform asymptotic stability. Those theorems allow the derivative of Lyapunov functions to be positive on certain sets,relax the restriction about the rate of change of state variable in a system to be bounded in Marachkov's theorem and extend the related results in [4—7].展开更多
A new analytical method for springback of small curvature plane bending is addressed with unloading rule of classical elastic-plastic theory and principle of strain superposition.We start from strain analysis of plane...A new analytical method for springback of small curvature plane bending is addressed with unloading rule of classical elastic-plastic theory and principle of strain superposition.We start from strain analysis of plane bending which has initial curvature,and the theoretic derivation is on the widely applicable basic hypotheses.The results are unified to geometry constraint equations and springback equation of plane bending,which can be evolved to straight beam plane bending and pure bending.The expanding and setting round process is one of the situations of plane bending,which is a bend-stretching process of plane curved beam.In the present study,springback equation of plane bending is used to analyze the expanding and setting round process,and the results agree with the experimental data.With a reasonable prediction accuracy,this new analytical method for springback of plane bending can meet the needs of applications in engineering.展开更多
The optimal control problem of fully coupled forward-backward stochastic systems is put forward. A necessary condition, called maximum principle, for an optimal control of the problem with the control domain being con...The optimal control problem of fully coupled forward-backward stochastic systems is put forward. A necessary condition, called maximum principle, for an optimal control of the problem with the control domain being convex is proved.展开更多
In order to establish a well-balanced scheme, 2D shallow water equations were transformed and solved by using the Finite Volume Method (FVM) with unstructured mesh. The numerical flux from the interface between cell...In order to establish a well-balanced scheme, 2D shallow water equations were transformed and solved by using the Finite Volume Method (FVM) with unstructured mesh. The numerical flux from the interface between cells was computed with an exact Riemann solver, and the improved dry Riemann solver was applied to deal with the wet/dry problems. The model was verified through computing some typical examples and the tidal bore on the Qiantang River. The results show that the scheme is robust and accurate, and could be applied extensively to engineering problems.展开更多
The Mei symmetry of Tzénoff equations under the infinitesimal transformations of groups is studied in this paper. The definition and the criterion equations of the symmetry are given. If the symmetry is a Noether...The Mei symmetry of Tzénoff equations under the infinitesimal transformations of groups is studied in this paper. The definition and the criterion equations of the symmetry are given. If the symmetry is a Noether symmetry, then the Noether conserved quantity of the Tzénoff equations can be obtained by the Mei symmetry.展开更多
Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for...Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for the same class of nonlinear NDDE; Review of the related concepts and results on RK methods; Information on the asymptotic stability and global stability of the induced NRK method.展开更多
This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equatio...This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equation, generalized omega-equation, and departure from fields obtained by potential vorticity (PV) inversion. The basic thoery, assumptions as well as implementation and limitations for each of the tools are all discussed. These tools are applied to high—resolution mesoscale model data to assess the role of unbalanced dynamics in the generation of a mesoscale gravity wave event over the East Coast of the United States. Comparison of these tools in this case study shows that these various methods agree to a large extent with each other though they differ in details. Key words Unbalanced flow - Geostrophic adjustment - Gravity waves - Nonlinear balance equation - Potential vorticity inversion - Omega equations - Rossby number This research was conducted under support from NSF grant ATM-9700626 of the United States. The numerical computations described herein were performed on the Cray T90 at the North Carolina Supercomputing Center and the Cray supercomputer at the NCAR Scientific Computing Division, which also provided the initialization fields for the MM5. Thanks are extended to Mark Stoelinga at University of Washington for the RIP post-processing package.展开更多
Solving the nonlinear model of an aeroengine is converted to an optimization problem, and thus some optimization search methods can be used. An approach to solving the nonlinear model of an aeroengine by use of the g...Solving the nonlinear model of an aeroengine is converted to an optimization problem, and thus some optimization search methods can be used. An approach to solving the nonlinear model of an aeroengine by use of the genetic algorithm (GA) is developed. By comparison with N R algorithm, the accuracy of the values of initial guesses is not required for GA. Especially, the approach developed can be used when no priori knowledges of the values of initial guesses are availabe, and the convergence is improved significantly. GA properly combined with N R algorithm can increase the convergence speed.展开更多
B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-...B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differentialequations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs ofother type which appear in practice.展开更多
A series of stability, contractivity and asymptotic stability results of the solutions to nonlinear stiff Volterra functional differential equations (VFDEs) in Banach spaces is obtained, which provides the unified the...A series of stability, contractivity and asymptotic stability results of the solutions to nonlinear stiff Volterra functional differential equations (VFDEs) in Banach spaces is obtained, which provides the unified theoretical foundation for the stability analysis of solutions to nonlinear stiff problems in ordinary differential equations(ODEs), delay differential equations(DDEs), integro-differential equations(IDEs) and VFDEs of other type which appear in practice.展开更多
基金supported by the National Natural Science Foundation of China No.10671184
文摘An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.
文摘Based on the work of paper [1], we propose a modified Levenberg-Marquardt algoithm for solving singular system of nonlinear equations F(x) = 0, where F(x) : Rn - Rn is continuously differentiable and F'(x) is Lipschitz continuous. The algorithm is equivalent to a trust region algorithm in some sense, and the global convergence result is given. The sequence generated by the algorithm converges to the solution quadratically, if ||F(x)||2 provides a local error bound for the system of nonlinear equations. Numerical results show that the algorithm performs well.
基金supported by CERG Grants of Hong Kong Research Grant CouncilFRG grants of Hong Kong Baptist University
文摘The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical results confirm the theoretical prediction of the exponential rate of convergence. The result in this work seems to be the first successful spectral approach (with theoretical justification) for the Volterra type equations.
基金Project supported by the National Natural Science Foundation of China(No.10131040).
文摘This paper deals with nonlinear expectations. The author obtains a nonlinear gen- eralization of the well-known Kolmogorov’s consistent theorem and then use it to con- struct ?ltration-consistent nonlinear expectations via nonlinear Markov chains. Com- pared to the author’s previous results, i.e., the theory of g-expectations introduced via BSDE on a probability space, the present framework is not based on a given probabil- ity measure. Many fully nonlinear and singular situations are covered. The induced topology is a natural generalization of Lp-norms and L∞-norm in linear situations. The author also obtains the existence and uniqueness result of BSDE under this new framework and develops a nonlinear type of von Neumann-Morgenstern representation theorem to utilities and present dynamic risk measures.
文摘High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of mu(= h/lambda, depth to deep-water wave length ratio) and epsilon(= a/h, wave amplitude to depth ratio) for velocity potential, particle velocity vector, pressure and the Boussinesq-type equations for surface elevation eta and horizontal velocity vector (U) over right arrow at any given level in water are given. Then, the exact explicit expressions to the fourth order of mu are derived. Finally, the linear solutions of eta, (U) over right arrow, C (phase-celerity) and C-g (group velocity) for a constant water depth are obtained. Compared with the Airy theory, excellent results can be found even for a water depth as large as the wave legnth. The present high-order models are applicable to nonlinear regular and irregular waves in water of any varying depth (from shallow to deep) and bottom slope (from mild to steep).
文摘It is a comparatively convenient technique to investigate the motion of a particle with the help of the differential geometry the-ory,rather than directly decomposing the motion in the Cartesian coordinates.The new model of three-dimensional (3D) guidance problem for interceptors is presented in this paper,based on the classical differential geometry curve theory.Firstly,the kinematical equations of the line of sight (LOS) are gained by carefully investigating the rotation principle of LOS,the kinematic equations of LOS are established,and the concepts of curvature and torsion of LOS are proposed.Simultaneously,the new relative dynamic equations between interceptor and target are constructed.Secondly,it is found that there is an instan-taneous rotation plane of LOS (IRPL) in the space,in which two-dimensional (2D) guidance laws could be constructed to solve 3D interception guidance problems.The spatial 3D true proportional navigation (TPN) guidance law could be directly introduced in IRPL without approximation and linearization for dimension-reduced 2D TPN.In addition,the new series of augmented TPN (APN) and LOS angular acceleration guidance laws (AAG) could also be gained in IRPL.After that,the dif-ferential geometric guidance commands (DGGC) of guidance laws in IRPL are advanced,and we prove that the guidance commands in arc-length system proposed by Chiou and Kuo are just a special case of DGGC.Moreover,the performance of the original guidance laws will be reduced after the differential geometric transformation.At last,an exoatmospheric intercep-tion is taken for simulation to demonstrate the differential geometric modeling proposed in this paper.
基金Supported by the National Natural Science Foundation of China(No.10671184).
文摘A nonconforming H^1-Calerkin mixed finite element method is analyzed for Sobolev equations on anisotropic meshes. The error estimates are obtained without using Ritz-Volterra projection.
基金The project supported by the National Natural Science Foundation of China(10232010 and 10472008)
文摘The equations of motion of an insect with flapping wings are derived and then simplified to that of a flying body using the "rigid body" assumption. On the basis of the simplified equations of motion, the longitudinal dynamic flight stability of four insects (hoverfly, cranefly, dronefly and hawkmoth) in hovering flight is studied (the mass of the insects ranging from 11 to 1,648 mg and wingbeat frequency from 26 to 157Hz). The method of computational fluid dynamics is used to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis are used to solve the equations of motion. The validity of the "rigid body" assumption is tested and how differences in size and wing kinematics influence the applicability of the "rigid body" assumption is investigated. The primary findings are: (1) For insects considered in the present study and those with relatively high wingbeat frequency (hoverfly, drone fly and bumblebee), the "rigid body" assumption is reasonable, and for those with relatively low wingbeat frequency (cranefly and howkmoth), the applicability of the "rigid body" assumption is questionable. (2) The same three natural modes of motion as those reported recently for a bumblebee are identified, i.e., one unstable oscillatory mode, one stable fast subsidence mode and one stable slow subsidence mode. (3) Approximate analytical expressions of the eigenvalues, which give physical insight into the genesis of the natural modes of motion, are derived. The expressions identify the speed derivative Mu (pitching moment produced by unit horizontal speed) as the primary source of the unstable oscillatory mode and the stable fast subsidence mode and Zw (vertical force produced by unit vertical speed) as the primary source of the stable slow subsidence mode.
文摘We study a new algorithm for solvingparabolic partial differential equations(PDEs)and backward stochastic differential equations(BSDEs)in high dimension,which is based on an analogy between the BSDE and reinforcement learning with the gradient of the solution playing the role of the policy function,and the loss function given by the error between the prescribed terminal condition and the solution of the BSDE.The policy function is then approximated by a neural network,as is done in deep reinforcement learning.Numerical results using TensorFlow illustrate the efficiency and accuracy of the studied algorithm for several 100-dimensional nonlinear PDEs from physics and finance such as the Allen–Cahn equation,the Hamilton–Jacobi–Bellman equation,and a nonlinear pricing model for financial derivatives.
文摘This paper gives several fundamental theorems for the stability, uniform stability, asymptotic stability and uniform asymptotic stability. Those theorems allow the derivative of Lyapunov functions to be positive on certain sets,relax the restriction about the rate of change of state variable in a system to be bounded in Marachkov's theorem and extend the related results in [4—7].
基金supported by the National Natural Science Foundation of China(Grant No.50805126)the Natural Science Foundation of Hebei Province(Grant No.E2009000389)
文摘A new analytical method for springback of small curvature plane bending is addressed with unloading rule of classical elastic-plastic theory and principle of strain superposition.We start from strain analysis of plane bending which has initial curvature,and the theoretic derivation is on the widely applicable basic hypotheses.The results are unified to geometry constraint equations and springback equation of plane bending,which can be evolved to straight beam plane bending and pure bending.The expanding and setting round process is one of the situations of plane bending,which is a bend-stretching process of plane curved beam.In the present study,springback equation of plane bending is used to analyze the expanding and setting round process,and the results agree with the experimental data.With a reasonable prediction accuracy,this new analytical method for springback of plane bending can meet the needs of applications in engineering.
文摘The optimal control problem of fully coupled forward-backward stochastic systems is put forward. A necessary condition, called maximum principle, for an optimal control of the problem with the control domain being convex is proved.
基金Project supported by the Natural Science Foundation of Zhejiang Province (Grant No: M403054).
文摘In order to establish a well-balanced scheme, 2D shallow water equations were transformed and solved by using the Finite Volume Method (FVM) with unstructured mesh. The numerical flux from the interface between cells was computed with an exact Riemann solver, and the improved dry Riemann solver was applied to deal with the wet/dry problems. The model was verified through computing some typical examples and the tidal bore on the Qiantang River. The results show that the scheme is robust and accurate, and could be applied extensively to engineering problems.
基金Project supported by the National Natural Science Foundation of China (Grant No 10372053) and the Natural Science Foundation of Henan Province, China (Grant No 0311011400). We are grateful for the instruction and help of Professor Mei F X, in Beijing Institute of Technology.
文摘The Mei symmetry of Tzénoff equations under the infinitesimal transformations of groups is studied in this paper. The definition and the criterion equations of the symmetry are given. If the symmetry is a Noether symmetry, then the Noether conserved quantity of the Tzénoff equations can be obtained by the Mei symmetry.
文摘Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for the same class of nonlinear NDDE; Review of the related concepts and results on RK methods; Information on the asymptotic stability and global stability of the induced NRK method.
文摘This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equation, generalized omega-equation, and departure from fields obtained by potential vorticity (PV) inversion. The basic thoery, assumptions as well as implementation and limitations for each of the tools are all discussed. These tools are applied to high—resolution mesoscale model data to assess the role of unbalanced dynamics in the generation of a mesoscale gravity wave event over the East Coast of the United States. Comparison of these tools in this case study shows that these various methods agree to a large extent with each other though they differ in details. Key words Unbalanced flow - Geostrophic adjustment - Gravity waves - Nonlinear balance equation - Potential vorticity inversion - Omega equations - Rossby number This research was conducted under support from NSF grant ATM-9700626 of the United States. The numerical computations described herein were performed on the Cray T90 at the North Carolina Supercomputing Center and the Cray supercomputer at the NCAR Scientific Computing Division, which also provided the initialization fields for the MM5. Thanks are extended to Mark Stoelinga at University of Washington for the RIP post-processing package.
基金Aeronautic Science Foundation of China ( 0 0 C5 2 0 3 0 ) and National Doctoral Education Foundation ( 2 0 0 0 0 2 870 1)
文摘Solving the nonlinear model of an aeroengine is converted to an optimization problem, and thus some optimization search methods can be used. An approach to solving the nonlinear model of an aeroengine by use of the genetic algorithm (GA) is developed. By comparison with N R algorithm, the accuracy of the values of initial guesses is not required for GA. Especially, the approach developed can be used when no priori knowledges of the values of initial guesses are availabe, and the convergence is improved significantly. GA properly combined with N R algorithm can increase the convergence speed.
文摘B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differentialequations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs ofother type which appear in practice.
基金This work was supported by the National High-Tech ICF Committee in Chinathe National Natural Science Foundation of China(Grant No.10271100).
文摘A series of stability, contractivity and asymptotic stability results of the solutions to nonlinear stiff Volterra functional differential equations (VFDEs) in Banach spaces is obtained, which provides the unified theoretical foundation for the stability analysis of solutions to nonlinear stiff problems in ordinary differential equations(ODEs), delay differential equations(DDEs), integro-differential equations(IDEs) and VFDEs of other type which appear in practice.