A new type of Galerkin finite element for first-order initial-value problems(IVPs)is proposed.Both the trial and test functions employ the same m-degreed polynomials.The adjoint equation is used to eliminate one degre...A new type of Galerkin finite element for first-order initial-value problems(IVPs)is proposed.Both the trial and test functions employ the same m-degreed polynomials.The adjoint equation is used to eliminate one degree of freedom(DOF)from the test function,and then the so-called condensed test function and its consequent condensed Galerkin element are constructed.It is mathematically proved and numerically verified that the condensed element produces the super-convergent nodal solutions of O(h^(2m+2)),which is equivalent to the order of accuracy by the conventional element of degree m+1.Some related properties are addressed,and typical numerical examples of both linear and nonlinear IVPs of both a single equation and a system of equations are presented to show the validity and effectiveness of the proposed element.展开更多
A fractal approximation algorithm is developed to obtain approximate solutions to an inverse initial-value problem IVP(inverse IVP) for the differential equation. Numerical computational results are presented to demon...A fractal approximation algorithm is developed to obtain approximate solutions to an inverse initial-value problem IVP(inverse IVP) for the differential equation. Numerical computational results are presented to demonstrate the effectiveness of this algorithm for solving inverse IVP for a class of specific differential equations.展开更多
Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to...Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to result in uncomfortably many numerical operations and high memory requirements. It is shown in this work that performance is substantially enhanced by the introduction of algorithms for temporal and spatial subdomains in combination with sparse matrix methods. The accuracy and efficiency of the recently developed time spectral, generalized weighted residual method (GWRM) are compared to that of the explicit Lax-Wendroff and implicit Crank-Nicolson methods. Three initial-value PDEs are employed as model problems;the 1D Burger equation, a forced 1D wave equation and a coupled system of 14 linearized ideal magnetohydrodynamic (MHD) equations. It is found that the GWRM is more efficient than the time-stepping methods at high accuracies. The advantageous scalings Nt<sup style="margin-left:-6px;">1.0Ns<sup style="margin-left:-6px;">1.43 and Nt<sup style="margin-left:-6px;">0.0Ns<sup style="margin-left:-6px;">1.08 were obtained for CPU time and memory requirements, respectively, with Nt and Ns denoting the number of temporal and spatial subdomains. For time-averaged solution of the two-time-scales forced wave equation, GWRM performance exceeds that of the finite difference methods by an order of magnitude both in terms of CPU time and memory requirement. Favorable subdomain scaling is demonstrated for the MHD equations, indicating a potential for efficient solution of advanced initial-value problems in, for example, fluid mechanics and MHD.展开更多
Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisatio...Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisation of weighted residual methods to the time and parameter domains [1]. A semi-analytical Chebyshev polynomial ansatz is employed, and the problem reduces to determine the coefficients of the ansatz from linear or nonlinear algebraic systems of equations. In order to avoid large memory storage and computational cost, it is preferable to subdivide the temporal and spatial domains into subdomains. Methods and examples of this article demonstrate how this can be achieved.展开更多
A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this general...A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.展开更多
The concept of approximate generalized conditional symmetry （A GCS） as a generalization to both approximate Lie point symmetry and generalized conditional symmetry is introduced, and it is applied to study the pertu...The concept of approximate generalized conditional symmetry （A GCS） as a generalization to both approximate Lie point symmetry and generalized conditional symmetry is introduced, and it is applied to study the perturbed nonlinear diffusion-convection equations. Complete classification of those perturbed equations which admit cerrain types of AGCSs is derived. Some approximate solutions to the resulting equations can be obtained via the AGCS and the corresponding unperturbed equations.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
In this paper, an almost P-stable two-step sixth-order Hybrid method with phase-lag of order infinity and a class explicit eighth-order Obreckoff methods with phase-lag of order 10-24 are developed for the numerical i...In this paper, an almost P-stable two-step sixth-order Hybrid method with phase-lag of order infinity and a class explicit eighth-order Obreckoff methods with phase-lag of order 10-24 are developed for the numerical integration of the special second-order periodic initial-value problems. These methods have the advantage of higher algebraic order and considerably smaller phase-tag compared with some methods in [1-6]. Numerical examples indicate that these new methods are more accurate than methods developed by [1-6].展开更多
The approximate generalized conditional symmetries of nonlinear filtration equation with a small parameter are studied. The initial-value problem of the equation can be transformed to perturbed Cauchy problem of pertu...The approximate generalized conditional symmetries of nonlinear filtration equation with a small parameter are studied. The initial-value problem of the equation can be transformed to perturbed Cauchy problem of perturbed first-order ordinary dif- ferential equations systems and approximate solutions are obtained by using these symmetries.展开更多
Direct and inverse scattering problems connected with the wave equation in non-homogeneous bounded domains constitute challenging actual subjects for both mathematicians and engineers. Among them one can mention, for ...Direct and inverse scattering problems connected with the wave equation in non-homogeneous bounded domains constitute challenging actual subjects for both mathematicians and engineers. Among them one can mention, for example, inverse source problems in seismology, nondestructive archeological probing, mine prospecting, inverse initial-value problems in acoustic tomography, etc. In spite of its crucial importance, almost all of the available rigorous investigations concern the case of unbounded simple domains such as layered planar or cylindrical or spherical structures. The main reason for the lack of the works related to non-homogeneous bounded structures is the extreme complexity of the explicit expressions of the Green’s functions. The aim of the present work consists in discovering some universal properties of the Green’s functions in question, which reduce enormously the difficulties arising in various applications. The universality mentioned here means that the properties are not depend on the geometrical and physical properties of the configuration. To this end one considers first the case when the domain is partially-homogeneous. Then the results are generalized to the most general case. To show the importance of the universal properties in question, they are applied to an inverse initial-value problem connected with photo-acoustic tomography.展开更多
We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the...We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the method.展开更多
Many problems in applied mathematics lead to ordinary differential equation. In this paper, a considerable refinement and improvement of the Euler＇s method obtained using PSO （particle swarm optimization） was prese...Many problems in applied mathematics lead to ordinary differential equation. In this paper, a considerable refinement and improvement of the Euler＇s method obtained using PSO （particle swarm optimization） was presented. PSO is a technique based on the cooperation between particles. The exchange of information between these particles allows to resolve difficult problems. This approach is carefully handled and tested with an illustrated example.展开更多
The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study a...The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.展开更多
According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified new way proposed by Luo, the un-conventional Hamilton-type variational principles of holonomic...According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified new way proposed by Luo, the un-conventional Hamilton-type variational principles of holonomic conservative system in analytical mechanics can be established systematically. This unconventional Hamilton-type variational principle can fully characterize the initial-value problem of analytical mechanics, so that it is an important innovation for the Hamilton-type variational principle. In this paper, an important integral relation is given, which can be considered as the expression of the generalized principle of virtual work for analytical mechanics in mechanics. Based on this relation, it is possible not only to obtain the principle of virtual work of holonomic conservative system in analytical mechanics, but also to derive systematically the complementary functionals for three-field and two-field unconventional variational principles, and the functional for the one-field one by the generalized Legendre transformation given in this paper. Further, with this new approach, the intrinsic relationship among various principles can be explained clearly. Meanwhile, the unconventional Hamilton-type variational principles of nonholonomic conservative system in analytical mechanics can also be established systematically in this paper.展开更多
In this paper, the authors consider the Harry-Dym equation on the line with decaying initial value. They construct the solution of the Harry-Dym equation via the solution of a 2 × 2 matrix Riemann-Hilbert problem...In this paper, the authors consider the Harry-Dym equation on the line with decaying initial value. They construct the solution of the Harry-Dym equation via the solution of a 2 × 2 matrix Riemann-Hilbert problem in the complex plane. Further, onecusp soliton solution is expressed in terms of the Riemann-Hilbert problem.展开更多
In this paper, the Fokas unified method is used to analyze the initial-boundary value for the ChenLee-Liu equation i？tu ＋ ？xxu-i｜u2｜？xu = 0 on the half line（-∞, 0] with decaying initial value. Assuming that th...In this paper, the Fokas unified method is used to analyze the initial-boundary value for the ChenLee-Liu equation i？tu ＋ ？xxu-i｜u2｜？xu = 0 on the half line（-∞, 0] with decaying initial value. Assuming that the solution u（x, t） exists, we show that it can be represented in terms of the solution of a matrix Riemann-Hilbert problem formulated in the plane of the complex spectral parameter λ. The jump matrix has explicit（x, t） dependence and is given in terms of the spectral functions{a（λ）, b（λ）}and{A（λ）, B（λ）}, which are obtained from the initial data u0（x） = u（x, 0） and the boundary data g0（t） = u（0, t）, g1（t） = ux（0, t）, respectively. The spectral functions are not independent,but satisfy a so-called global relation.展开更多
In this paper, two families of high accuracy explicit two-step methods with minimal phase-lag are developed for the numerical integration of special secondorder periodic initial-value problems. In comparison with some...In this paper, two families of high accuracy explicit two-step methods with minimal phase-lag are developed for the numerical integration of special secondorder periodic initial-value problems. In comparison with some methods in [1, 4,6], the advantage of these methods has a higher accuracy and minimal phaselag. The methods proposed in this paper can be considered as a generalization of some methods in [1,3,4]. Numerical examples indicate that these new methods are generally more accurate than the methods used in [3,6].展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.51878383 and51378293)。
文摘A new type of Galerkin finite element for first-order initial-value problems(IVPs)is proposed.Both the trial and test functions employ the same m-degreed polynomials.The adjoint equation is used to eliminate one degree of freedom(DOF)from the test function,and then the so-called condensed test function and its consequent condensed Galerkin element are constructed.It is mathematically proved and numerically verified that the condensed element produces the super-convergent nodal solutions of O(h^(2m+2)),which is equivalent to the order of accuracy by the conventional element of degree m+1.Some related properties are addressed,and typical numerical examples of both linear and nonlinear IVPs of both a single equation and a system of equations are presented to show the validity and effectiveness of the proposed element.
文摘A fractal approximation algorithm is developed to obtain approximate solutions to an inverse initial-value problem IVP(inverse IVP) for the differential equation. Numerical computational results are presented to demonstrate the effectiveness of this algorithm for solving inverse IVP for a class of specific differential equations.
文摘Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to result in uncomfortably many numerical operations and high memory requirements. It is shown in this work that performance is substantially enhanced by the introduction of algorithms for temporal and spatial subdomains in combination with sparse matrix methods. The accuracy and efficiency of the recently developed time spectral, generalized weighted residual method (GWRM) are compared to that of the explicit Lax-Wendroff and implicit Crank-Nicolson methods. Three initial-value PDEs are employed as model problems;the 1D Burger equation, a forced 1D wave equation and a coupled system of 14 linearized ideal magnetohydrodynamic (MHD) equations. It is found that the GWRM is more efficient than the time-stepping methods at high accuracies. The advantageous scalings Nt<sup style="margin-left:-6px;">1.0Ns<sup style="margin-left:-6px;">1.43 and Nt<sup style="margin-left:-6px;">0.0Ns<sup style="margin-left:-6px;">1.08 were obtained for CPU time and memory requirements, respectively, with Nt and Ns denoting the number of temporal and spatial subdomains. For time-averaged solution of the two-time-scales forced wave equation, GWRM performance exceeds that of the finite difference methods by an order of magnitude both in terms of CPU time and memory requirement. Favorable subdomain scaling is demonstrated for the MHD equations, indicating a potential for efficient solution of advanced initial-value problems in, for example, fluid mechanics and MHD.
文摘Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisation of weighted residual methods to the time and parameter domains [1]. A semi-analytical Chebyshev polynomial ansatz is employed, and the problem reduces to determine the coefficients of the ansatz from linear or nonlinear algebraic systems of equations. In order to avoid large memory storage and computational cost, it is preferable to subdivide the temporal and spatial domains into subdomains. Methods and examples of this article demonstrate how this can be achieved.
文摘A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10371098 and 10447007, the Natural Science Foundation of Shaanxi Province （No 2005A13）, and the Special Research Project of Educational Department of Shaanxi Province （No 03JK060）.
文摘The concept of approximate generalized conditional symmetry （A GCS） as a generalization to both approximate Lie point symmetry and generalized conditional symmetry is introduced, and it is applied to study the perturbed nonlinear diffusion-convection equations. Complete classification of those perturbed equations which admit cerrain types of AGCSs is derived. Some approximate solutions to the resulting equations can be obtained via the AGCS and the corresponding unperturbed equations.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
文摘In this paper, an almost P-stable two-step sixth-order Hybrid method with phase-lag of order infinity and a class explicit eighth-order Obreckoff methods with phase-lag of order 10-24 are developed for the numerical integration of the special second-order periodic initial-value problems. These methods have the advantage of higher algebraic order and considerably smaller phase-tag compared with some methods in [1-6]. Numerical examples indicate that these new methods are more accurate than methods developed by [1-6].
基金Supported by the National Natural Science Foundation of China(11326167,11226195)the Natural Science Foundation of Henan Province(142300410354)Foundation of Henan Educational Committee(13A110119,15B110012)
文摘The approximate generalized conditional symmetries of nonlinear filtration equation with a small parameter are studied. The initial-value problem of the equation can be transformed to perturbed Cauchy problem of perturbed first-order ordinary dif- ferential equations systems and approximate solutions are obtained by using these symmetries.
文摘Direct and inverse scattering problems connected with the wave equation in non-homogeneous bounded domains constitute challenging actual subjects for both mathematicians and engineers. Among them one can mention, for example, inverse source problems in seismology, nondestructive archeological probing, mine prospecting, inverse initial-value problems in acoustic tomography, etc. In spite of its crucial importance, almost all of the available rigorous investigations concern the case of unbounded simple domains such as layered planar or cylindrical or spherical structures. The main reason for the lack of the works related to non-homogeneous bounded structures is the extreme complexity of the explicit expressions of the Green’s functions. The aim of the present work consists in discovering some universal properties of the Green’s functions in question, which reduce enormously the difficulties arising in various applications. The universality mentioned here means that the properties are not depend on the geometrical and physical properties of the configuration. To this end one considers first the case when the domain is partially-homogeneous. Then the results are generalized to the most general case. To show the importance of the universal properties in question, they are applied to an inverse initial-value problem connected with photo-acoustic tomography.
文摘We present the optimal homotopy asymptotic method (OHAM) to find the numerical solution of the second order initial value problems of Bratu-type. We solve some examples to illustrate the validity and efficiency of the method.
文摘Many problems in applied mathematics lead to ordinary differential equation. In this paper, a considerable refinement and improvement of the Euler＇s method obtained using PSO （particle swarm optimization） was presented. PSO is a technique based on the cooperation between particles. The exchange of information between these particles allows to resolve difficult problems. This approach is carefully handled and tested with an illustrated example.
文摘The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.
基金the National Natural Science Foundation of China(Grant Nos. 10172097 & 10272034)the Science Foundation for Doctoral Program of Ministry of Education of China (Grant No. 20030558025)
文摘According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified new way proposed by Luo, the un-conventional Hamilton-type variational principles of holonomic conservative system in analytical mechanics can be established systematically. This unconventional Hamilton-type variational principle can fully characterize the initial-value problem of analytical mechanics, so that it is an important innovation for the Hamilton-type variational principle. In this paper, an important integral relation is given, which can be considered as the expression of the generalized principle of virtual work for analytical mechanics in mechanics. Based on this relation, it is possible not only to obtain the principle of virtual work of holonomic conservative system in analytical mechanics, but also to derive systematically the complementary functionals for three-field and two-field unconventional variational principles, and the functional for the one-field one by the generalized Legendre transformation given in this paper. Further, with this new approach, the intrinsic relationship among various principles can be explained clearly. Meanwhile, the unconventional Hamilton-type variational principles of nonholonomic conservative system in analytical mechanics can also be established systematically in this paper.
基金supported by the National Natural Science Foundation of China(No.11271079)
文摘In this paper, the authors consider the Harry-Dym equation on the line with decaying initial value. They construct the solution of the Harry-Dym equation via the solution of a 2 × 2 matrix Riemann-Hilbert problem in the complex plane. Further, onecusp soliton solution is expressed in terms of the Riemann-Hilbert problem.
基金Supported by the National Natural Science Foundation of China(No.11271008,61072147,11671095)SDUST Research Fund(No.2018TDJH101)
文摘In this paper, the Fokas unified method is used to analyze the initial-boundary value for the ChenLee-Liu equation i？tu ＋ ？xxu-i｜u2｜？xu = 0 on the half line（-∞, 0] with decaying initial value. Assuming that the solution u（x, t） exists, we show that it can be represented in terms of the solution of a matrix Riemann-Hilbert problem formulated in the plane of the complex spectral parameter λ. The jump matrix has explicit（x, t） dependence and is given in terms of the spectral functions{a（λ）, b（λ）}and{A（λ）, B（λ）}, which are obtained from the initial data u0（x） = u（x, 0） and the boundary data g0（t） = u（0, t）, g1（t） = ux（0, t）, respectively. The spectral functions are not independent,but satisfy a so-called global relation.
文摘In this paper, two families of high accuracy explicit two-step methods with minimal phase-lag are developed for the numerical integration of special secondorder periodic initial-value problems. In comparison with some methods in [1, 4,6], the advantage of these methods has a higher accuracy and minimal phaselag. The methods proposed in this paper can be considered as a generalization of some methods in [1,3,4]. Numerical examples indicate that these new methods are generally more accurate than the methods used in [3,6].
