Approximate solution of differential equations with the use of asymptotic polynomials

We consider an approximate solution of differential equations with initial and boundary conditions. To find a solution, we use asymptotic polynomials Q _n ^f (x) of the first kind based on Chebyshev polynomials T _n (x) of the first kind and asymptotic polynomials G _n ^f (x) of the second kind based on Chebyshev polynomials U _n (x) of the second kind. We suggest most efficient algorithms for each of these solutions. We find classes of functions for which the approximate solution converges to the exact one. The remainder is represented as an expansion in linear functionals {L _n ^f } in the first case and {M _n ^f } in the second case, whose decay rate depends on the properties of functions describing the differential equation.     

Linear functional equations of the first, second, and third kind in L 2

Under consideration are the functional equations of the first, second, and third kind with operators in wide classes of linear continuous operators in L ₂ containing all integral operators. We propose methods for reducing these equations by linear invertible changes either to linear integral equations of the first kind with nuclear operators or to equivalent linear integral equations of the second kind with quasidegenerate Carleman kernels. Some various approximate methods of solution are applicable to the so-obtained integral equations.     

Legendre spectral Galerkin method for second-kind Volterra integral equations

The Legendre spectral Galerkin method for the Volterra integral equations of the second kind is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the L ₂ norm) will decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical examples are given to illustrate the theoretical results.      

On the integral equation method for the plane mixed boundary value problem of the Laplacian

Although the plane boundary value problem for the Laplacian with given Dirichlet data on one part Γ₂ and given Neumann data on the remaining part Γ₂ of the boundary is the simplest case of mixed boundary value problems, we present several applications in classical mathematical physics. Using Green's formula the problem is converted into a system of Fredholm integral equations for the yet unknown values of the solution u on Γ₂ and the also desired values of the normal derivatie on Γ₁. One of these equations has principal part of the second kind, whereas that one of the other is of the first kind. Since any improvement of constructive methods requires higher regularity of u but, on the other hand, grad u possesses singularities at the collision points Γ₁ ∩ Γ₂ even for C^∞ data, u is decomposed into special singular terms and a regular rest. This is incorporated into the integral equations and the modified system is solved in appropriate Sobolev spaces. The solution of the system requires to solve a Fredholm equation of the first kind on the arc Γ₂ providing an improvement of regularity for the smooth part of u. Since the integral equations form a strongly elliptic system of pseudodifferential operators, the Galerkin procedure converges. Using regular finite element functions on Γ₁ and Γ₂ augmented by the special singular functions we obtain optimal order of asymptotic convergence in the norm corresponding to the energy norm of u and also superconvergence as well as high orders in smoother norms if the given data are smooth (and not the solution).     

On the discretization of differential and Volterra integral equations with variable delay

In this paper we analyze the attainable order ofm-stage implicit (collocation-based) Runge-Kutta methods for differential equations and Volterra integral equations of the second kind with variable delay of the formqt (0<q<1). It will be shown that, in contrast to equations without delay, or equations with constant delay, collocation at the Gauss (-Legendre) points will no longer yield the optimal (local) orderO(h ^2m ). This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC Research Grant OGP0009406).     

Central Units in ℤCp,q

Let C_{p, q} be the semi-direct product of a cyclic group of order q by a cyclic group of order p, and ?C_{p, q} the integral group ring of C_{p, q}. In this article, firstly, we describe the group of normalized central units of ?C_{p, q} as a direct product of two subgroups that we call units of first kind and of second kind. For a class of prime numbers that we call good primes, we construct a multiplicatively independent set which generates the group of units of first kind. Finally, we construct a set of multiplicatively independent units which generates the units of second kind for a larger class of primes.     

Integral equations of the third kind with unbounded operators

We consider linear functional equations of the third kind in L ₂ with arbitrary measurable coefficients and unbounded integral operators with kernels satisfying broad conditions. We propose methods for reducing these equations by linear continuous invertible transformations either to equivalent integral equations of the first kind with nuclear operators or to equivalent integral equations of the second kind with quasidegenerate Carleman kernels. To the integral equations obtained after the reduction, one can apply various exact and approximate methods of solution; in particular, the two approximate methods developed in this article.     

On the limiting behaviour of solutions to boundary integral equations associated with time harmonic wave equations for small frequencies

The treatment of boundary value problems for Helmholtz equation and for the time harmonic Maxwell's equations by boundary integral equations leads to integral equations of the second kind which are uniquely solvable for small positive frequencies λ. However, the integral equations obtained in the limiting case λ = 0 which are related to boundary value problems of potential theory in general are not uniquely solvable since the corresponding boundary value problems are not. By first considering in a general setting of a Banach space X the limiting behaviour of solutions ?_λ to the equation ?_λ – K λ ? λ = fλ as λ → 0 where {K_λ: X → X, λ ∈ (0,α)}, α > 0, denotes a family of compact linear operators such that I - K_λ (I identity) is bijective for λ∈(0,α) whilst I - K₀ is not and ‖ K_λ – K₀‖ →, 0, ‖f_λ – f0‖ → 0, λ → 0, and then applying the results to the boundary integral operators, the limiting behaviour of the integral equations is considered. Thus, the results obtained by Mac Camey for the Helmholtz equation are extended to the case of non-connected boundaries and Werner's results on the integral equations for the Maxwell's equations are extended to the case of multiply connected boundaries.     

On linear Fredholm integro-differential equations of first kind

We study the question whether linear one-dimensional integro-differential equations with constant limits of integration (equations of Fredholm type) containing no free differential expression (equations of first kind) can be reduced to integral equations of first kind and to Fredholm integro-differential equations of second kind.Translated fromVychislitel'naya i Prikladnaya Matematika, Issue 71, 1990, pp. 20–27.     

The exterior dirichlet problem for the Laplace equation

Using a standard application of Green's theorem, the exterior Dirichlet problem for the Laplace equation in three dimensions is reduced to a pair of integral equations. One integral equation is of the second kind and the other is of the first kind. It is known that the integral equation of the second kind is not uniquely solvable, however, it has been demonstrated that the pair of integral equations has a unique solution. The present approach is based on the observation that the known function appearing in the integral equation of the second kind lies in a certain Banach space E which is a proper subspace of the Banach space of continuous complex-valued functions equipped with the maximum norm. Furthermore, it can be shown that the related integral operator when restricted to E has spectral radius less than unity. Consequently, a particular solution to the integral equation of the second kind can be obtained by the method of successive approximations and the unique solution to the problem is then obtained by using the integral equation of the first kind. Comparisons are made between the present algorithm and other known constructive methods. Finally, an example is supplied to illustrate the method of this paper.     

Partition Lattice q-Analogs Related to q-Stirling Numbers

We construct a family of partially ordered sets (posets) that are q-analogs of the set partition lattice. They are different from the q-analogs proposed by Dowling [5]. One of the important features of these posets is that their Whitney numbers of the first and second kind are just the q-Stirling numbers of the first and second kind, respectively. One member of this family [4] can be constructed using an interpretation of Milne [9] for S[n, k] as sequences of lines in a vector space over the Galois field F _q. Another member is constructed so as to mirror the partial order in the subspace lattice.     

Controllability of systems of Stokes equations with one control force: existence of insensitizing controls

In this paper we establish some exact controllability results for systems of two parabolic equations of the Stokes kind. In a first part, we prove the existence of insensitizing controls for the L² norm of the solutions and the curl of solutions of linear Stokes equations. Then, in the limit case where one can expect null controllability to hold for a system of two Stokes equations (namely, when the coupling terms concern first and second order derivatives, respectively), we prove this result for some general couplings.     

A differential game with two pursuers and one evader

This paper is concerned with a coplanar pursuit-evasion problem in which a faster evaderE with constant speedw>1 must pass between two pursuersP ₁,P ₂ having unit speed, the payoff being the distance of closest approach to either one of the pursuers. The control variables are the directions of the velocities ofP ₁,P ₂, andE. The path equations are integrated, and a closed-form solution is obtained in terms of elliptic functions of the first and second kind. A closed-loop solution is given graphically in several diagrams, for different values ofw.     

The ε-Complexity of the Approximate Solution of Integral Equations with Isotropic Kernels

The exact order of ε-complexity is determined in L _p (2 ≤p < ∞) spaces for the second kind of Fredholm integral equations with kernels belonging to an isotropic Sobolev class. Received December 10, 1999, Accepted April 19, 2002     

Lp – Lq Decay Estimates fo the Solutions of Strictly Hyperbolic Equations of Second Order with Increasing in Time Coefficients

One of the most important questions in the theory of nonlinear wave equations is that for global existence of solutions. An essential tool is the Strichartz inequality for special solutions of the wave equation.In the last time different results were proved generalizing the classical one of Strichartz. In the present paper L_p – L_q estimates are proved for the solutions of strictly hyperbolic equations of second order with time dependent coefficients where these are unbounded at infinity. In the first step the WKB method is applied to the construction of a fundamental system of solutions for ordinary differential equations depending on a parameter. In a second step the method of stationary phase yields the asymptotical behaviour of Fourier multipliers with nonstandard phase functions depending on a parameter.     

Fast Construction of Optimal Circulant Preconditioners for Matrices from the Fast Dense Matrix Method

In this paper, we consider solving non-convolution type integral equations by the preconditioned conjugate gradient method. The fast dense matrix method is a fast multiplication scheme that provides a dense discretization matrix A approximating a given integral equation. The dense matrix A can be constructed in O(n) operations and requires only O(n) storage where n is the size of the matrix. Moreover, the matrix-vector multiplication A xcan be done in O(n log n) operations. Thus if the conjugate gradient method is used to solve the discretized system, the cost per iteration is O(n log n) operations. However, for some integral equations, such as the Fredholm integral equations of the first kind, the system will be ill-conditioned and therefore the convergence rate of the method will be slow. In these cases, preconditioning is required to speed up the convergence rate of the method. A good choice of preconditioner is the optimal circulant preconditioner which is the minimizer of C – A _F in Frobenius norm over all circulant matrices C. It can be obtained by taking arithmetic averages of all the entries of A and therefore the cost of constructing the preconditioner is of O(n ²) operations for general dense matrices. In this paper, we develop an O(n log n) method of constructing the preconditioner for dense matrices A obtained from the fast dense matrix method. Application of these ideas to boundary integral equations from potential theory will be given. These equations are ill-conditioned whereas their optimal circulant preconditioned equations will be well-conditioned. The accuracy of the approximation A, the fast construction of the preconditioner and the fast convergence of the preconditioned systems will be illustrated by numerical examples.     

Fast numerical solution of fredholm integral equations with stationary kernels

A fast recursive matrix method for the numerical solution of Fredholm integral equations with stationary kernels is derived. IfN denotes the number of nodal points, the complexity of the algorithm isO(N ²), which should be compared toO(N ³) for conventional algorithms for solving such problems. The method is related to fast algorithms for inverting Toeplitz matrices.Applications to equations of the first and second kind as well as miscellaneous problems are discussed and illustrated with numerical examples. These show that the theoretical improvement in efficiency is indeed obtained, and that no problems with numerical stability or accuracy are encountered.     

NEW STRUCTURES FOR NON-SELFSIMILAR SOLUTIONS OF MULTI-DIMENSIONAL CONSERVATION LAWS

In this article, we get non-selfsimilar elementary waves of the conservation laws in another kind of view, which is different from the usual self-similar transformation. The solution has different global structure. This article is divided into three parts. The first part is introduction. In the second part, we discuss non-selfsimilar elementary waves and their interactions of a class of twodimensional conservation laws. In this case, we consider the case that the initial discontinuity is parabola with u＋ 〉 0, while explicit non-selfsirnilar rarefaction wave can be obtained. In the second part, we consider the solution structure of case u＋ 〈 0. The new solution structures are obtained by the interactions between different elementary waves, and will continue to interact with other states. Global solutions would be very different from the situation of one dimension.     

The initial boundary-value problem with a free surface condition for the penalized equations of aqueous solutions of polymers

In this paper, we study some nonlocal problems for the Kelvin-Voight equations (1) and the penalized Kelvin-Voight equations (2): the first and second initial boundary-value problems and the first and second time periodic boundary problems. We prove that these problems have global smooth solutions of the classW _∞ ¹ (ℝ⁺;W ₂ ^2+k (Ω)),k=1,2,...;Ω⊂ℝ³. Bibliography: 25 titles. Dedicated to N. N. Uraltseva on her jubilee Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 221, 1995, pp. 185–207. Translated by N. A. Karazeeva.     

Structure of the frequency spectrum of nonaxisymmetric vibrations of a transversally isotropic hollow ball

We consider nonaxisymmetric vibrations of a hollow sphere made of a spherically Isotropic material. The original boundary-value problem is reduced to a system of ordinary differential equations by introducing new unknown junctions and expanding the sought solution in a series in spherical functions. The nonaxisymmetric problem of vibrations of a hollow ball decomposes into two independent problems that describe modes of first and second kind. The frequency spectrum for modes of first and second kind is analyzed.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 59, pp. 32–35, 1986.