In this paper, we first present constructing a Lyapunov function for (1. 1) and then we show the asymptotic stability in the large of the trivial solution x=0 for case p≡ 0,and the boundedness result of the sol...In this paper, we first present constructing a Lyapunov function for (1. 1) and then we show the asymptotic stability in the large of the trivial solution x=0 for case p≡ 0,and the boundedness result of the solutions of (1 .1 ) for case p≠0. These results improve sveral well-known results.展开更多
An alternating direction implicit (ADI) Galerkin method with moving finite element spaces is formulated for a class of second order hyperbolic equations in two space variables. A priori H 1 error estimate is derived.
A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illus...A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
In this work, we apply a hyperbola function method to solve the nonlinear family of third order Korteweg-de Vries equations. Exact travelling wave solutions are obtained and expressed in terms of hyperbolic functions ...In this work, we apply a hyperbola function method to solve the nonlinear family of third order Korteweg-de Vries equations. Exact travelling wave solutions are obtained and expressed in terms of hyperbolic functions and trigonometric functions. The method used is a promising method to solve other nonlinear evaluation equations.展开更多
In this paper, by defining an appropriate Lyapunov functional, we obtain sufficient conditions for which all solutions of certain real non-autonomous third order nonlinear differential equations are asymptotically sta...In this paper, by defining an appropriate Lyapunov functional, we obtain sufficient conditions for which all solutions of certain real non-autonomous third order nonlinear differential equations are asymptotically stable and bounded. The results obtained improve and extend some known results in the literature.展开更多
We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, ...We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.展开更多
The dynamic characteristics and the efficiency of the Ostrowski’s method allow it to be crowned as an excellent tool for solving nonlinear problems.This article shows different versions of the classic method that all...The dynamic characteristics and the efficiency of the Ostrowski’s method allow it to be crowned as an excellent tool for solving nonlinear problems.This article shows different versions of the classic method that allow it to be applied to a wide range of engineering problems.Among them stands out the derivative-free definition applying divided differences,the introduction of memory and its extension to the resolution of nonlinear systems of equations.All of these versions are compared in a numerical simulations section where the results obtained are compared with other classic methods.展开更多
In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations （SDEs）. Six fully explicit methods （TSM 1 -- TSM 6） are given in this paper. Their order of strong convergence ...In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations （SDEs）. Six fully explicit methods （TSM 1 -- TSM 6） are given in this paper. Their order of strong convergence is proved. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of Ito SDEs.展开更多
A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of ...A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.展开更多
In this article, Crank-Nicolson method is used to study the variable order fractional cable equation. The variable order fractional derivatives are described in the Riemann- Liouville and the Griinwald-Letnikov sense....In this article, Crank-Nicolson method is used to study the variable order fractional cable equation. The variable order fractional derivatives are described in the Riemann- Liouville and the Griinwald-Letnikov sense. The stability analysis of the proposed technique is discussed. Numerical results are provided and compared with exact solutions to show the accuracy of the proposed technique.展开更多
In this paper,we introduce the dissipative spectral methods(DSM)for the first order linear hyperbolic equations in one dimension.Specifically,we consider the Fourier DSM for periodic problems and the Legendre DSM for ...In this paper,we introduce the dissipative spectral methods(DSM)for the first order linear hyperbolic equations in one dimension.Specifically,we consider the Fourier DSM for periodic problems and the Legendre DSM for equations with the Dirichlet boundary condition.The error estimates of the methods are shown to be quasioptimal for variable-coefficients equations.Numerical results are given to verify high accuracy of the DSM and to compare the proposed schemes with some high performance methods,showing some superiority in long-term integration for the periodic case and in dealing with limited smoothness near or at the boundary for the Dirichlet case.展开更多
Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bou...Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bounded or unbounded domain in R~3 and obtain the error estimates of the approximate solution in energy norm and local maximum norm.展开更多
A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordin...A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.展开更多
Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a cl...Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.展开更多
The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework i...The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework is the automatic way of getting arbitrarily high order methods,which can be put in the Runge-Kutta(RK)form.The drawback is the larger computational cost with respect to the most used RK methods.To reduce such cost,in an explicit setting,we propose an efcient modifcation:we introduce interpolation processes between the DeC iterations,decreasing the computational cost associated to the low order ones.We provide the Butcher tableaux of the new modifed methods and we study their stability,showing that in some cases the computational advantage does not afect the stability.The fexibility of the novel modifcation allows nontrivial applications to PDEs and construction of adaptive methods.The good performances of the introduced methods are broadly tested on several benchmarks both in ODE and PDE contexts.展开更多
In this paper, Goursat’s problems for: linear and nonlinear hyperbolic equations of second-order, systems of nonlinear hyperbolic equations and fourth-order linear hyperbolic equations in which the attached condition...In this paper, Goursat’s problems for: linear and nonlinear hyperbolic equations of second-order, systems of nonlinear hyperbolic equations and fourth-order linear hyperbolic equations in which the attached conditions are given on the characteristics curves are transformed in such a manner that the Adomian decomposition method (ADM) can be applied. Some examples with closed-form solutions are studied in detail to further illustrate the proposed technique, and the results obtained indicate this approach is indeed practical and efficient.展开更多
文摘In this paper, we first present constructing a Lyapunov function for (1. 1) and then we show the asymptotic stability in the large of the trivial solution x=0 for case p≡ 0,and the boundedness result of the solutions of (1 .1 ) for case p≠0. These results improve sveral well-known results.
基金the National Natural Sciences Foundation of China
文摘An alternating direction implicit (ADI) Galerkin method with moving finite element spaces is formulated for a class of second order hyperbolic equations in two space variables. A priori H 1 error estimate is derived.
文摘A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
文摘In this work, we apply a hyperbola function method to solve the nonlinear family of third order Korteweg-de Vries equations. Exact travelling wave solutions are obtained and expressed in terms of hyperbolic functions and trigonometric functions. The method used is a promising method to solve other nonlinear evaluation equations.
文摘In this paper, by defining an appropriate Lyapunov functional, we obtain sufficient conditions for which all solutions of certain real non-autonomous third order nonlinear differential equations are asymptotically stable and bounded. The results obtained improve and extend some known results in the literature.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)the National Basic Research Program of China(Grant No.2012CB025904)the State Key Laboratory of Science and Engineering Computing and the Center for High Performance Computing of Northwestern Polytechnical University,China
文摘We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.
文摘The dynamic characteristics and the efficiency of the Ostrowski’s method allow it to be crowned as an excellent tool for solving nonlinear problems.This article shows different versions of the classic method that allow it to be applied to a wide range of engineering problems.Among them stands out the derivative-free definition applying divided differences,the introduction of memory and its extension to the resolution of nonlinear systems of equations.All of these versions are compared in a numerical simulations section where the results obtained are compared with other classic methods.
文摘In this paper we discuss two-stage Miistein methods for solving Ito stochastic differential equations （SDEs）. Six fully explicit methods （TSM 1 -- TSM 6） are given in this paper. Their order of strong convergence is proved. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of Ito SDEs.
文摘A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.
文摘In this article, Crank-Nicolson method is used to study the variable order fractional cable equation. The variable order fractional derivatives are described in the Riemann- Liouville and the Griinwald-Letnikov sense. The stability analysis of the proposed technique is discussed. Numerical results are provided and compared with exact solutions to show the accuracy of the proposed technique.
基金supported by National Natural Science Foundation of China(11171209)Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50101)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20060280010)Graduate Innovative Foundation of Shanghai University(SHUCX091048).
文摘In this paper,we introduce the dissipative spectral methods(DSM)for the first order linear hyperbolic equations in one dimension.Specifically,we consider the Fourier DSM for periodic problems and the Legendre DSM for equations with the Dirichlet boundary condition.The error estimates of the methods are shown to be quasioptimal for variable-coefficients equations.Numerical results are given to verify high accuracy of the DSM and to compare the proposed schemes with some high performance methods,showing some superiority in long-term integration for the periodic case and in dealing with limited smoothness near or at the boundary for the Dirichlet case.
基金China State Major Key Project for Basic Researches
文摘Combining difference method and boundary integral equation method,we propose a new numerical method for solving initial-boundary value problem of second order hyperbolic partial differential equations defined on a bounded or unbounded domain in R~3 and obtain the error estimates of the approximate solution in energy norm and local maximum norm.
基金The National Natural Science Foundation of China (No10471023)
文摘A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.
基金supported by ARO grant W911NF-04-1-0291,NSF grant DMS-0510345 and AFOSR grant FA9550-05-1-0123.
文摘Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.
文摘The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework is the automatic way of getting arbitrarily high order methods,which can be put in the Runge-Kutta(RK)form.The drawback is the larger computational cost with respect to the most used RK methods.To reduce such cost,in an explicit setting,we propose an efcient modifcation:we introduce interpolation processes between the DeC iterations,decreasing the computational cost associated to the low order ones.We provide the Butcher tableaux of the new modifed methods and we study their stability,showing that in some cases the computational advantage does not afect the stability.The fexibility of the novel modifcation allows nontrivial applications to PDEs and construction of adaptive methods.The good performances of the introduced methods are broadly tested on several benchmarks both in ODE and PDE contexts.
文摘In this paper, Goursat’s problems for: linear and nonlinear hyperbolic equations of second-order, systems of nonlinear hyperbolic equations and fourth-order linear hyperbolic equations in which the attached conditions are given on the characteristics curves are transformed in such a manner that the Adomian decomposition method (ADM) can be applied. Some examples with closed-form solutions are studied in detail to further illustrate the proposed technique, and the results obtained indicate this approach is indeed practical and efficient.
