Second order partial differential equations for gradients of orthogonal polynomials in two variables

In this work, we introduce the classical orthogonal polynomials in two variables as the solutions of a matrix second order partial differential equation involving matrix polynomial coefficients, the usual gradient operator, and the divergence operator. Here we show that the successive gradients of these polynomials also satisfy a matrix second order partial differential equation closely related to the first one.     

On the absence of eigenvalues of Maxwell and Lamé systems in exterior domains

The arguments showing non‐existence of eigensolutions to exterior‐boundary value problems associated with systems—such as the Maxwell and Lamé system—rely on showing that such solutions would have to have compact support and therefore—by a unique continuation property—cannot be non‐trivial. Here we will focus on the first part of the argument. For a class of second order elliptic systems it will be shown that L²‐solutions in exterior domains must have compact support. Both the asymptotically isotropic Maxwell system and the Lamé system with asymptotically decaying perturbations can be reduced to this class of elliptic systems. Copyright © 2004 John Wiley & Sons, Ltd.     

Nonelement boundary representation with Bézier surface patches for 3D linear elasticity problems in parametric integral equation system (PIES) and its solving using Lagrange polynomials

In this article, we present a strategy of using rectangular and triangular Bézier surface patches for nonelement representation of 3D boundary geometries for problems of linear elasticity. The boundary generated in this way is directly incorporated in the parametric integral equation system (PIES), which has been developed by the authors. The boundary values on each surface patch are approximated by Lagrange polynomials. Three illustrative examples are presented to confirm the effectiveness of the proposed boundary representation in connection with PIES and to show good accuracy of numerical results.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 34: 51–79, 2018     

Fractional‐order orthogonal Bernstein polynomials for numerical solution of nonlinear fractional partial Volterra integro‐differential equations

In this paper, a new two‐dimensional fractional polynomials based on the orthonormal Bernstein polynomials has been introduced to provide an approximate solution of nonlinear fractional partial Volterra integro‐differential equations. For this aim, the fractional‐order orthogonal Bernstein polynomials (FOBPs) are constructed, and its operational matrices of integration, fractional‐order integration, and derivative in the Caputo sense and product operational matrix are derived. These operational matrices are utilized to reduce the under study problem to a nonlinear system of algebraic equations. Using the approximation of FOBPs, the convergence analysis and error estimate associated to the proposed problem have been investigated. Finally, several examples are included to clarify the validity, efficiency, and applicability of the proposed technique via FOBPs approximation.     

Large time behavior for nonlinear higher order convection–diffusion equations

We study the large time asymptotic behavior, in L^p (1p∞), of higher derivatives D^γu(t) of solutions of the nonlinear equation

(1)

where the integers n and θ are bigger than or equal to 1, a is a constant vector in with . The function ψ is a nonlinearity such that and ψ(0)=0, and is a higher order elliptic operator with nonsmooth bounded measurable coefficients on . We also establish faster decay when .     

Extension theorems for Stokes and Lamé equations for nearly incompressible media and their applications to numerical solution of problems with highly discontinuous coefficients

We prove extension theorems in the norms described by Stokes and Lamé operators for the three‐dimensional case with periodic boundary conditions. For the Lamé equations, we show that the extension theorem holds for nearly incompressible media, but may fail in the opposite limit, i.e. for case of absolutely compressible media. We study carefully the latter case and associate it with the Cosserat problem. Extension theorems serve as an important tool in many applications, e.g. in domain decomposition and fictitious domain methods, and in analysis of finite element methods. We consider an application of established extension theorems to an efficient iterative solution technique for the isotropic linear elasticity equations for nearly incompressible media and for the Stokes equations with highly discontinuous coefficients. The iterative method involves a special choice for an initial guess and a preconditioner based on solving a constant coefficient problem. Such preconditioner allows the use of well‐known fast algorithms for preconditioning. Under some natural assumptions on smoothness and topological properties of subdomains with small coefficients, we prove convergence of the simplest Richardson method uniform in the jump of coefficients. For the Lamé equations, the convergence is also uniform in the incompressible limit. Our preliminary numerical results for two‐dimensional diffusion problems show fast convergence uniform in the jump and in the mesh size parameter. Copyright © 2002 John Wiley & Sons, Ltd.     

Scalar and matrix Riemann–Hilbert approach to the strong asymptotics of Padé approximants and complex orthogonal polynomials with varying weight

We describe methods for the derivation of strong asymptotics for the denominator polynomials and the remainder of Padé approximants for a Markov function with a complex and varying weight. Two approaches, both based on a Riemann–Hilbert problem, are presented. The first method uses a scalar Riemann–Hilbert boundary value problem on a two-sheeted Riemann surface, the second approach uses a matrix Riemann–Hilbert problem. The result for a varying weight is not with the most general conditions possible, but the loss of generality is compensated by an easier and transparent proof.     

A remark on the definitions of viscosity solutions for the integro-differential equations with Lévy operators

We prove the equivalence of the three different definitions of the viscosity solution for the integro-differential equation with the Lévy operator. The two of the definitions are known in the preceding works of the author and the others, and the last one is new. A construction of a sequence of the approximating test functions to the subsolution (or the supersolution) is indispensable for the proof, and it is done explicitly in the paper.     

On comparison of series and numerical solutions for second Painlevé equation

This attempt presents the series solution of second Painlevé equation by homotopy analysis method (HAM). Comparison of HAM solution is provided with that of the Adomian decomposition method (ADM), homotopy perturbation method (HPM), analytic continuation method, and Legendre Tau method. It is revealed that there is very good agreement between the analytic continuation and HAM solutions when compared with ADM, HPM, and Legendre Tau solutions. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Existence of multiplicity harmonic and subharmonic solutions for second‐order quasilinear equation via Poincaré‐Birkhoff twist theorem

In this paper, we investigate the existence and multiplicity of harmonic and subharmonic solutions for second‐order quasilinear equation where , g satisfies the superlinear condition at infinity. We prove that the given equation possesses harmonic and subharmonic solutions by using the phase‐plane analysis methods and a generalized version of the Poincaré‐Birkhoff twist theorem.     

High‐order Padé and singly diagonally Runge‐Kutta schemes for linear ODEs,application to wave propagation problems

In this article, we address the problem of constructing high‐order implicit time schemes for wave equations. We consider two classes of one‐step A‐stable schemes adapted to linear Ordinary Differential Equation (ODE). The first class, which is not dissipative is based on the diagonal Padé approximant of exponential function. For this class, the obtained schemes have the same stability function as Gauss Runge‐Kutta (Gauss RK) schemes. They have the advantage to involve the solution of smaller linear systems at each time step compared to Gauss RK. The second class of schemes are constructed such that they require the inversion of a unique linear system several times at each time step like the Singly Diagonally Runge‐Kutta (SDIRK) schemes. While the first class of schemes is constructed for an arbitrary order of accuracy, the second‐class schemes is given up to order 12. The performance assessment we provide shows a very good level of accuracy for both classes of schemes, and the great interest of considering high‐order time schemes that are faster. The diagonal Padé schemes seem to be more accurate and more robust.     

A Padé family of iterations for the matrix sign function and related problems

In this paper we consider the Pad'e family of iterations for computing the matrix sign function and the Padé family of iterations for computing the matrix p‐sector function. We prove that all the iterations of the Padé family for the matrix sign function have a common convergence region. It completes a similar result of Kenney and Laub for half of the Padé family. We show that the iterations of the Padé family for the matrix p‐sector function are well defined in an analogous common region, depending on p. For this purpose we proved that the Padé approximants to the function (1?z)^?σ, 0<σ<1, are a quotient of hypergeometric functions whose poles we have localized. Furthermore we proved that the coefficients of the power expansion of a certain analytic function form a positive sequence and in a special case this sequence has the log‐concavity property. Copyright © 2011 John Wiley & Sons, Ltd.