On non-linear evolution equations of higher order – the introduction and application of a novel computational approach

The purpose of the present paper is to introduce a new computational algebraic procedure that can easily be applied to derive class of solutions of non-linear partial differential equations (nPDE) especially of higher order.The crucial step needs an auxiliary variable satisfying some ordinary differential equations (ODE) of first order containing sine, cosine and their hyperbolic varieties introducing to the first time.General transformations are given to determine class of solutions explicitly.The validity and reliability of the method is tested by its application to some important non-linear evolution equations leading to new class of solutions with physical significance.Nevertheless it should be emphasised that this techniques do not need the solution of complicate nODEs as in the case of similarity reduction.Further, the algorithm works efficiently, is clear structured and can be used in any applications independent of the order of the nPDE. For computational purposes the method is appropriate to rewrite it in any computer languages.Therefore, the given novel algebraic approach is suitable for a wider class of nPDE in order to augment the solution manifold by a straightforward alternative approach.     

THE GROWTH OF SOLUTIONS OF SYSTEMS OF COMPLEX NONLINEAR ALGEBRAIC DIFFERENTIAL EQUATIONS

We investigate the problem of growth order of solutions of a type of systems of non-linear algebraic differential equations, and extend some results of the growth order of solutions of algebraic differential equations to systems of algebraic differential equations.     

The Cavalcante-Tenenblat equation - Does the equation admits physical significance?

The main purpose of the given paper is to analyze a less studied third order non-linear partial differential equation, the so-called Cavalcante-Tenenblat equation (CTE) in the following form: .Since general class of solutions are of basic interest a complete characterization of the group properties is given. The traveling wave ‘ansatz’ restricts the solution manifold to special class of solutions and hence, a generalize algorithm is necessary.We determine the Lie point symmetry vector fields and calculate similarity ‘ansätze’. Further, we also derive a few non-linear transformations and some similarity solutions are obtained explicitly. Due to the complexity of some similarity solutions a numerical procedure is of advantage.Moreover, the non-classical case (potential symmetries) is studied to the first time and further, we show how the CTE leads to approximate symmetries and we apply the method to the first time. We call the disturbed equation the CTE-ε equation and we show how to derive new class of solutions.Finally, the equation does not pass the Painlevé-test and is therefore not soluble by the Inverse Scattering Transform Method (IST).Hence, suitable alternative (algebraic) approaches are necessary to derive class of solutions explicitly.     

Complete group classification of systems of two linear second‐order ordinary differential equations: the algebraic approach

We give a complete group classification of the general case of linear systems of two second‐order ordinary differential equations. The algebraic approach is used to solve the group classification problem for this class of equations. This completes the results in the literature on the group classification of two linear second‐order ordinary differential equations including recent results which give a complete group classification treatment of such systems. We show that using the algebraic approach leads to the study of a variety of cases in addition to those already obtained in the literature. We illustrate that this approach can be used as a useful tool in the group classification of this class of equations. A discussion of the subsequent cases and results is given. Copyright © 2014 John Wiley & Sons, Ltd.     

复非线性代数微分方程解的增长级(英文)

本文研究了一类非线性微分方程的解的增长级的问题.利用亚纯函数的Nevanlinna值分布理论和Wiman-Valiron整函数理论的方法,获得了比以往文献更为精确,更为一般的结论,推广了Gol’dberg,Barsegian,Hayman以及Korhonen等作者的一些结果.     

A New Approach for Solving a Class of Delay Fractional Partial Differential Equations

In this article, a new numerical approach has been proposed for solving a class of delay time-fractional partial differential equations. The approximate solutions of these equations are considered as linear combinations of Müntz–Legendre polynomials with unknown coefficients. Operational matrix of fractional differentiation is provided to accelerate computations of the proposed method. Using Padé approximation and two-sided Laplace transformations, the mentioned delay fractional partial differential equations will be transformed to a sequence of fractional partial differential equations without delay. The localization process is based on the space-time collocation in some appropriate points to reduce the fractional partial differential equations into the associated system of algebraic equations which can be solved by some robust iterative solvers. Some numerical examples are also given to confirm the accuracy of the presented numerical scheme. Our results approved decisive preference of the Müntz–Legendre polynomials with respect to the Legendre polynomials.     

Rational ODEs with Movable Algebraic Singularities

A class of second-order rational ordinary differential equations, admitting certain families of formal algebraic series solutions, is considered. For all solutions of these equations, it is shown that any movable singularity that can be reached by analytic continuation along a finite-length curve is an algebraic branch point. The existence of these formal series expansions is straightforward to determine for any given equation in the class considered. We apply the theorem to a family of equations, admitting different kinds of algebraic singularities. As a further application we recover the known fact for generic values of parameters that the only movable singularities of solutions of the Painlevé equations   P_II – P   _VI   are poles.     

Entropy solution theory for fractional degenerate convection–diffusion equations

We study a class of degenerate convection-diffusion equations with a fractional non-linear diffusion term. This class is a new, but natural, generalization of local degenerate convection-diffusion equations, and include anomalous diffusion equations, fractional conservation laws, fractional porous medium equations, and new fractional degenerate equations as special cases. We define weak entropy solutions and prove well-posedness under weak regularity assumptions on the solutions, e.g. uniqueness is obtained in the class of bounded integrable solutions. Then we introduce a new monotone conservative numerical scheme and prove convergence toward the entropy solution in the class of bounded integrable BV functions. The well-posedness results are then extended to non-local terms based on general Lévy operators, connections to some fully non-linear HJB equations are established, and finally, some numerical experiments are included to give the reader an idea about the qualitative behavior of solutions of these new equations.     

The existence of nonoscillatory solutions of higher order nonlinear neutral equations

In this work, we consider the existence of nonoscillatory solutions of variable coefficient higher order nonlinear neutral differential equations. Our results include as special cases some well-known results for linear and nonlinear equations of first, second and higher order. We use the Banach contraction principle to obtain new sufficient conditions for the existence of nonoscillatory solutions.     

Index-aware model order reduction for differential-algebraic equations

We introduce a model order reduction (MOR) procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. The decoupling procedure in differential and algebraic part is based on the projector and matrix chain which leads to the definition of tractability index. The differential part can be reduced by using any MOR method, we use Krylov-based projection methods to illustrate our approach. The reduction on the differential part induces a reduction on the algebraic part. In this paper, we present the method for index-1 differential-algebraic equations. We implement numerically this procedure and show numerical evidence of its validity.     

二阶代数微分方程亚纯解的增长性估计

本文研究了代数微分方程亚纯解的增长级.运用正规族理论,给出了某类二阶代数微分方程亚纯解的增长级的一个估计,该估计依赖于方程的有理函数系数.推广了2001年廖良文与杨重骏的一个结果.     

Algebraic and Differential Nonlinear Superposition Formulas

Knowledge of the Lax pair and the Darboux transformation for a completely integrable system provides an iterative approach for generating exact solutions. This approach involves solving for the eigenfunction of the Lax pair at each step. But this process can be considerably simplified using the Bäcklund transformation and Bianchi's permutability theorem. This allows constructing the so-called nonlinear superposition formula, which provides a new solution of the system in terms of three previous solutions. The advantage of this approach is that the differential order of the nonlinear superposition formulas is lower than that of the Lax pairs, and in some cases, these formulas reduce to algebraic equations. We consider the construction of new nonlinear superposition formulas in the form of both differential equations and algebraic equations.     

非线性波方程准确孤立波解的符号计算

该文将机械化数学方法应用于偏微分方程领域，建立了构造一类非线性发展方程孤立波解的一种统一算法，并在计算机数学系统上加以实现，推导出了一批非线性发展方程的精确孤立波解．算法的基本原理是利用非线性发展方程孤立波解的局部性特点，将孤立波表示为双曲正切函数的多项式．从而将非线性发展方程（组）的求解问题转化为非线性代数方程组的求解问题．利用吴文俊消元法在计算机代数系统上求解非线性代数方程组，最终获得非线性发展方程（组）的准确孤立波解.     

A new Jacobi operational matrix: An application for solving fractional differential equations

In this paper, we derived the shifted Jacobi operational matrix (JOM) of fractional derivatives which is applied together with spectral tau method for numerical solution of general linear multi-term fractional differential equations (FDEs). A new approach implementing shifted Jacobi operational matrix in combination with the shifted Jacobi collocation technique is introduced for the numerical solution of nonlinear multi-term FDEs. The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations which greatly simplifying the problem. The proposed methods are applied for solving linear and nonlinear multi-term FDEs subject to initial or boundary conditions, and the exact solutions are obtained for some tested problems. Special attention is given to the comparison of the numerical results obtained by the new algorithm with those found by other known methods.     

Backward error,sensitivity, and refinement of computed solutions of algebraic Riccati equations

In this paper, a new backward error criterion, together with a sensitivity measure, is presented for assessing solution accuracy of nonsymmetric and symmetric algebraic Riccati equations (AREs). The usual approach to assessing reliability of computed solutions is to employ standard perturbation and sensitivity results for linear systems and to extend them further to AREs. However, such methods are not altogether appropriate since they do not take account of the underlying structure of these matrix equations. The approach considered here is to first compute the backward error of a computed solution X? that measures the amount by which data must be perturbed so that X? is the exact solution of the perturbed original system. Conventional perturbation theory is used to define structured condition numbers that fully respect the special structure of these matrix equations. The new condition number, together with the backward error of computed solutions, provides accurate estimates for the sensitivity of solutions. Optimal perturbations are then used in an iterative refinement procedure to give further more accurate approximations of actual solutions. The results are derived in their most general setting for nonsymmetric and symmetric AREs. This in turn offers a unifying framework through which it is possible to establish similar results for Sylvester equations, Lyapunov equations, linear systems, and matrix inversions.     

Polynomial solutions of certain differential equations arising in physics

Conditions for the existence of polynomial solutions of certain second‐order differential equations have recently been investigated by several authors. In this paper, a new algorithmic procedure is given to determine necessary and sufficient conditions for a differential equation with polynomial coefficients containing parameters to admit polynomial solutions and to compute these solutions. The effectiveness of this approach is illustrated by applying it to determine new solutions of several differential equations of current interest. A comparative analysis is given to demonstrate the advantage of this algorithmic procedure over existing software. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimal Control of One-Dimensional Partial Differential Algebraic Equations with Applications

We present an approach to compute optimal control functions in dynamic models based on one-dimensional partial differential algebraic equations (PDAE). By using the method of lines, the PDAE is transformed into a large system of usually stiff ordinary differential algebraic equations and integrated by standard methods. The resulting nonlinear programming problem is solved by the sequential quadratic programming code NLPQL. Optimal control functions are approximated by piecewise constant, piecewise linear or bang-bang functions. Three different types of cost functions can be formulated. The underlying model structure is quite flexible. We allow break points for model changes, disjoint integration areas with respect to spatial variable, arbitrary boundary and transition conditions, coupled ordinary and algebraic differential equations, algebraic equations in time and space variables, and dynamic constraints for control and state variables. The PDAE is discretized by difference formulae, polynomial approximations with arbitrary degrees, and by special update formulae in case of hyperbolic equations. Two application problems are outlined in detail. We present a model for optimal control of transdermal diffusion of drugs, where the diffusion speed is controlled by an electric field, and a model for the optimal control of the input feed of an acetylene reactor given in form of a distributed parameter system.     

When nonautonomous equations are equivalent to autonomopus ones

We consider nonlinear systems of first order partial differential equations admitting at least two one-parameter Lie groups of transformations with commuting infinitesimal operators. Under suitable conditions it is possible to introduce a variable transformation based on canonical variables which reduces the model in point to autonomous form. Remarkably, the transformed system may admit constant solutions to which there correspond non-constant solutions of the original model. The results are specialized to the case of first order quasilinear systems admitting either dilatation or spiral groups of transformations and a systematic procedure to characterize special exact solutions is given. At the end of the paper the equations of axi-symmetric gas dynamics are considered.     

Integrating factors, adjoint equations and Lagrangians

Integrating factors and adjoint equations are determined for linear and non-linear differential equations of an arbitrary order. The new concept of an adjoint equation is used for construction of a Lagrangian for an arbitrary differential equation and for any system of differential equations where the number of equations is equal to the number of dependent variables. The method is illustrated by considering several equations traditionally regarded as equations without Lagrangians. Noether's theorem is applied to the Maxwell equations.