Efficient hybrid method for solving special type of nonlinear partial differential equations

In this article, an efficient hybrid method has been developed for solving some special type of nonlinear partial differential equations. Hybrid method is based on tanh–coth method, quasilinearization technique and Haar wavelet method. Nonlinear partial differential equations have been converted into a nonlinear ordinary differential equation by choosing some suitable variable transformations. Quasilinearization technique is used to linearize the nonlinear ordinary differential equation and then the Haar wavelet method is applied to linearized ordinary differential equation. A tanh–coth method has been used to obtain the exact solutions of nonlinear ordinary differential equations. It is easier to handle nonlinear ordinary differential equations in comparison to nonlinear partial differential equations. A distinct feature of the proposed method is their simple applicability in a variety of two‐ and three‐dimensional nonlinear partial differential equations. Numerical examples show better accuracy of the proposed method as compared with the methods described in past. Error analysis and stability of the proposed method have been discussed.     

One method for solving systems of nonlinear partial differential equations

A method for reducing systems of partial differential equations to corresponding systems of ordinary differential equations is proposed. A system of equations describing two-dimensional, cylindrical, and spherical flows of a polytropic gas; a system of dimensionless Stokes equations for the dynamics of a viscous incompressible fluid; a system of Maxwell’s equations for vacuum; and a system of gas dynamics equations in cylindrical coordinates are studied. It is shown how this approach can be used for solving certain problems (shockless compression, turbulence, etc.).     

Multidimensional Inverse Scattering for First-Order Systems

A method for solving the inverse problem for a class of multidimensional first-order systems is given. The analysis yields equations which the scattering data must satisfy; these equations are natural candidates for characterizing admissible scattering data. The results are used to solve the multidimensional N-wave resonant interaction equations.     

Runge–Kutta Based Procedure for the Optimal Control of Differential-Algebraic Equations

A new approach for optimization of control problems defined by fully implicit differential-algebraic equations is described in the paper. The main feature of the approach is that system equations are substituted by discrete-time implicit equations resulting from the integration of the system equations by an implicit Runge–Kutta method. The optimization variables are parameters of piecewise constant approximations to control functions; thus, the control problem is reduced to the control space only. The method copes efficiently with problems defined by large-scale differential-algebraic equations.     

A computer science approach for solving elliptic differential equations

A method that combines computer science and numerical analysis techniques for the solution of elliptic partial differential equations in domains consisting of the union of rectangles is presented. It is applied to Laplace's equation and to the Navier–Stokes equations in several domains. The method is very flexible and can be easily modified to solve other equations. © 1995 John Wiley & Sons, Inc.     

Polygons of differential equations for finding exact solutions

A method for finding exact solutions of nonlinear differential equations is presented. Our method is based on the application of polygons corresponding to nonlinear differential equations. It allows one to express exact solutions of the equation studied through solutions of another equation using properties of the basic equation itself. The ideas of power geometry are used and developed. Our approach has a pictorial interpretation, which is illustrative and effective. The method can be also applied for finding transformations between solutions of differential equations. To demonstrate the method application exact solutions of several equations are found. These equations are: the Korteveg–de Vries–Burgers equation, the generalized Kuramoto–Sivashinsky equation, the fourth-order nonlinear evolution equation, the fifth-order Korteveg–de Vries equation, the fifth-order modified Korteveg–de Vries equation and the sixth-order nonlinear evolution equation describing turbulent processes. Some new exact solutions of nonlinear evolution equations are given.     

Boundary-Value Problems for Second-Order Equations and for Systems of First-Order Equations

The oblique derivative problem for harmonic functions under violation of the Shapiro–Lopatinsky condition is considered as well as some multi-dimensional analogues of the Cauchy–Riemann system. These problems are reduced to nonelliptic pseudo-differential equations. A method generalizing the regularization of singular integral equations is also presented.     

Convergence analysis of a family of Steffensen-type methods for generalized equations

A class of Steffensen-type algorithms for solving generalized equations on Banach spaces is proposed. Using well-known fixed point theorem for set-valued maps [A.L. Dontchev, W.W. Hager, An inverse function theorem for set-valued maps, Proc. Amer. Math. Soc. 121 (1994) 481-489] and some conditions on the first-order divided difference, we provide a local convergence analysis. We also study the perturbed problem and we present a new regula-falsi-type method for set-valued mapping. This study follows the works on the Secant-type method presented in [S. Hilout, A uniparametric Secant-type methods for nonsmooth generalized equations, Positivity (2007), submitted for publication; S. Hilout, A. Piétrus, A semilocal convergence of a Secant-type method for solving generalized equations, Positivity 10 (2006) 673-700] and extends the results related to the resolution of nonlinear equations [M.A. Hernández, M.J. Rubio, The Secant method and divided differences Hölder continuous, Appl. Math. Comput. 124 (2001) 139-149; M.A. Hernández, M.J. Rubio, Semilocal convergence of the Secant method under mild convergence conditions of differentiability, Comput. Math. Appl. 44 (2002) 277-285; M.A. Hernández, M.J. Rubio, ω-Conditioned divided differences to solve nonlinear equations, in: Monogr. Semin. Mat. García Galdeano, vol. 27, 2003, pp. 323-330; M.A. Hernández, M.J. Rubio, A modification of Newton's method for nondifferentiable equations, J. Comput. Appl. Math. 164/165 (2004) 323-330].     

Factorization method for linear and quasilinear singularly perturbed boundary value problems for ordinary differential equations

For linear singularly perturbed boundary value problems, we come up with a method that reduces solving a differential problem to a discrete (difference) problem. Difference equations, which are an exact analog of differential equations, are constructed by the factorization method. Coefficients of difference equations are calculated by solving Cauchy problems for first-order differential equations. In this case nonlinear Ricatti equations with a small parameter are solved by asymptotic methods, and solving linear equations reduces to computing quadratures. A solution for quasilinear singularly perturbed equations is obtained by means of an implicit relaxation method. A solution to a linearized problem is calculated by analogy with a linear problem at each iterative step. The method is tested against solutions to the known Lagerstrom-Cole problem.     

Block methods for nonlinear Volterra integral equations

A method used to solve linear Fredholm equations is adapted to obtain a generalization of block methods for solving nonlinear Volterra integral equations. Conditions for such methods to be convergent andA-stable are given along with some numerical examples.     

Mixed problem for pseudoparabolic integro-differential equation with degenerate kernel

A mixed problem for a certain nonlinear third-order intregro-differential equation of the pseudoparabolic type with a degenerate kernel is considered. The method of degenerate kernel is essentially used and developed and the Fourier method of variable separation is employed for this equation. A system of countable systems of algebraic equations is first obtained; after it is solved, a countable system of nonlinear integral equations is derived. The method of sequential approximations is used to prove the theorem on the unique solvability of the mixed problem.     

Domain-imbedding alternating direction method for linear elliptic equations on irregular regions using collocation

A new method is presented for solving elliptic partial differential equations over two-dimensional irregular regions. The scheme imbeds the irregular region in a rectangle, and then uses an alternating direction iteration to solve the resulting system of linear equations. Collocation with cubic Hermite splines is used for discretization. The method is shown to be equivalent to a multiboundary alternating direction method. A theory of convergence for a simplified case is given, details of implementation are discussed, and two numerical illustrations are presented. © 1993 John Wiley & Sons, Inc.     

Ritz method for equations with small parameters for higher derivatives

The problem of convergence of the Ritz method is considered for positive definite operational equations of the form $$A_\varepsilon u \equiv (\varepsilon {\rm A}_1 + A_0 )u = f$$ containing small parameters ? for the principal part. For specific natural conditions it is proved that the Ritz method, used for an approximate solution to such equations, converges to an exact solution in a metric with quadratic form uniformly with respect to the parameter ?.     

Quadratically Convergent Method for Simultaneously Approaching the Roots of Polynomial Solutions of a Class of Differential Equations: Application to Orthogonal Polynomials

This paper investigates the application of the method introduced by L. Pasquini (1989) for simultaneously approaching the zeros of polynomial solutions to a class of second-order linear homogeneous ordinary differential equations with polynomial coefficients to a particular case in which these polynomial solutions have zeros symmetrically arranged with respect to the origin. The method is based on a family of nonlinear equations which is associated with a given class of differential equations. The roots of the nonlinear equations are related to the roots of the polynomial solutions of differential equations considered. Newton's method is applied to find the roots of these nonlinear equations. In (Pasquini, 1994) the nonsingularity of the roots of these nonlinear equations is studied. In this paper, following the lines in (Pasquini, 1994), the nonsingularity of the roots of these nonlinear equations is studied. More favourable results than the ones in (Pasquini, 1994) are proven in the particular case of polynomial solutions with symmetrical zeros. The method is applied to approximate the roots of Hermite–Sobolev type polynomials and Freud polynomials. A lower bound for the smallest positive root of Hermite–Sobolev type polynomials is given via the nonlinear equation. The quadratic convergence of the method is proven. A comparison with a classical method that uses the Jacobi matrices is carried out. We show that the algorithm derived by the proposed method is sometimes preferable to the classical QR type algorithms for computing the eigenvalues of the Jacobi matrices even if these matrices are real and symmetric.     

A series method for solving nonlinear two-point boundary value problems

Summary A new method for solving nonlinear boundary value problems based on Taylor-type expansions generated by the use of Lie series is derived and applied to a set of test examples. A detailed discussion is given of the comparative performance of this method under various conditions. The method is of theoretical interest but is not applicable, in its present form, to real life problems; in particular, because of the algebraic complexity of the expressions involved, only scalar second order equations have been discussed, though in principle systems of equations could be similarly treated. A continuation procedure based on this method is suggested for future investigation.     

A Sinc–Collocation method for the linear Fredholm integro-differential equations

A Sinc–Collocation method for solving linear integro-differential equations of the Fredholm type is discussed. The integro-differential equations are reduced to a system of algebraic equations and Q-R method is used to establish numerical procedures. The convergence rate of the method is O( e ^{- k?N} )O{\left( {e^{{ - k{\sqrt N }}} } \right)} . Numerical results are included to confirm the efficiency and accuracy of the method even in the presence of singularities and a comparison with the rationalized Haar wavelet method is made.     

Unconditionally stable explicit methods for parabolic equations

Summary This paper discussesrational Runge-Kutta methods for stiff differential equations of high dimensions. These methods are explicit and in addition do not require the computation or storage of the Jacobian. A stability analysis (based onn-dimensional linear equations) is given. A second orderA ₀-stable method with embedded error control is constructed and numerical results of stiff problems originating from linear and nonlinear parabolic equations are presented.     

Stability of implicit-explicit linear multistep methods for ordinary and delay differential equations

Stability properties of implicit-explicit (IMEX) linear multistep methods for ordinary and delay differential equations are analyzed on the basis of stability regions defined by using scalar test equations. The analysis is closely related to the stability analysis of the standard linear multistep methods for delay differential equations. A new second-order IMEX method which has approximately the same stability region as that of the IMEX Euler method, the simplest IMEX method of order 1, is proposed. Some numerical results are also presented which show superiority of the new method.      

Adomian decomposition and Padé approximate for solving differential-difference equation

In this paper, we apply the Adomian decomposition method and Padé-approximate to solving the differential-difference equations (DDEs) for the first time. A simple but typical example is used to illustrate the validity and the great potential of the Adomian decomposition method (ADM) in solving DDEs. Comparisons are made between the results of the proposed method and exact solutions. The results show that ADM is an attractive method in solving the differential-difference equations.     

A simple moving mesh method for one- and two-dimensional phase-field equations

A simple moving mesh method is proposed for solving phase-field equations. The numerical strategy is based on the approach proposed in Li et al. [J. Comput. Phys. 170 (2001) 562–588] to separate the mesh-moving and PDE evolution. The phase-field equations are discretized by a finite-volume method, and the mesh-moving part is realized by solving the conventional Euler–Lagrange equations with the standard gradient-based monitors. Numerical results demonstrate the accuracy and effectiveness of the proposed algorithm.