Wave splitting of the Timoshenko beam equation in the time domain

In recent years, wave splitting in conjunction with invariant imbedding and Green's function techniques has been applied with great success to a number of interesting inverse and direct scattering problems. The aim of the present paper is to derive a wave splitting for the Timoshenko equation, a fourth order PDE of importance in beam theory. An analysis of the hyperbolicity of the Timoshenko equation and its, in a sense, less physical relatives-the Euler-Bernoulli and the Rayleigh equations-is also provided.     

Invariant imbedding applied to homogeneous two-point boundary-value problems with a singularity of the first kind

It is shown that previous work of Elder can be used to extend the version of invariant imbedding due to M. R. Scott to homogeneous (vector) differential systems having a singularity of the first kind. The boundary conditions considered consist of existence (finite) at the singularity and specified values for some subset of the dependent variables at a second point. The important special case of a second-order equation is discussed in some detail. Computational considerations are discussed and numerical examples are presented.     

Linear self-adjoint multipoint boundary value problems and related approximation schemes

Summary Linear self-adjoint multipoint boundary value problems are investigated. The case of the homogeneous equation is shown to lead to spline solutions, which are then utilized to construct a Green's function for the case of homogeneous boundary conditions. An approximation scheme is described in terms of the eigen-functions of the inverse of the Green's operator and is shown to be optimal in the sense of then-widths of Kolmogorov. Convergence rates are given and generalizations to more general boundary value problems are discussed.     

Invariant imbedding and multiple shooting for the solution of unstable linear boundary value problems

Invariant imbedding has been used to solve unstable linear boundary value problems for a few years. First this method is derived using the theory of characteristics; there the boundary value problem has to be imbedded in a problem of double dimension. If the corresponding Riccati equation has a critical length, one has to repeat the algorithm. A relation between this repeated invariant imbedding and multiple shooting is shown. In examples invariant imbedding, repeated invariant imbedding, multiple shooting and the superposition principle are compared.     

Excitation of a wave field in a triangular domain with impedance boundary conditions

The Helmholtz equation in a closed domain that is an equilateral triangle with inhomogeneous impedance boundary conditions is considered. A functional equation in which the unknown function is the Fourier-image of a wave field on the boundary of the domain is constructed. This functional equation is solved for the case of homogeneous boundary conditions (the problem on eigenvalues), as well as for the case of inhomogeneous boundary conditions in the absence of the resonance. Bibliography: 4 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 250, 1998, pp. 300–318. Translated by A. V. Shanin.     

A one-sweep method for linear elliptic equations over irregular domains

The determination of the configuration of equilibrium in a number of problems in mechanics and structures such as torsion, deflection of elastic membranes,etc., involve the solution of variational problems defined over irregular regions. This problem, in turn, may be reduced to the solution of elliptic differential equations subject to boundary conditions. In this paper, we study a method for the solution of such a problem when the region is of irregular shape. The method consists in solving the problem over a larger, imbedding, rectangular domain subject to appropriate constraints such as to satisfy the conditions of the original problem at the boundary. In this paper, we introduce the constraints by considering appropriate factors on the Green's function of the auxiliary problem. A conveniently discretized version of the problem is then treated by invariant imbedding, yielding some earlier results plus some new ones, namely, a direct one-sweep procedure that minimizes storage requirements. In addition, the present solution appears to be very convenient when the solution is required at a limited number of points. The derivations are specialized to Laplace's equation, but the method can be applied readily to general systems of second-order elliptic equations with no essential modifications. Finally, the existence of the necessary matrices in the imbedding equations is established.     

The method of lines for Poisson's equation with nonlinear or free boundary conditions

Summary The method of lines is used to solve Poisson's equation on an irregular domain with nonlinear or free boundary conditions. The partial differential equation is approximated by a system of second order ordinary differential equations subject to multi-point boundary conditions. The system is solved with an SOR iteration which employs invariant imbedding for each one dimensional problem. An application of the method to a boundary control problem and to a free surface problem arising in electrochemical machining is described. Finally, some theoretical convergence results are presented for a model problem with radiative boundary conditions on fixed boundaries.This work was supported by the U.S. Army Research Office under Grant DA-AG29-76-G-0261     

Application of invariant imbedding to linear multipoint boundary value problems with general boundary conditions

Two extensions of the usual application of invariant imbedding to the solution of linear boundary value problems are presented. The invariant imbedding formulation of a linear two point boundary value problem in which functional relationships are given between the variables at either one or both of the boundary points is presented. Also, extension of invariant imbedding to linear multipoint boundary value problems is given. Using these extensions singly or in combination, a general multipoint boundary value of linear ordinary differential equations can be solved. In addition, the problems of infinite initial conditions and / or indeterminate initial derivatives are resolved. Numerical examples demonstrate the feasibility and accuracy of the method.     

Invariant imbedding and the identification of aquifer parameters

The technique of invariant imbedding is applied to the problem of identifying the parameters in an unconfined aquifer system. This new technique is shown to be a very effective way of converting field observations based upon pumping tests into the desired aquifer parameters. The procedure is straightforward as it requires neither curve plotting nor graphical matching. The parameters to be identified are the hydraulic conductivity and specific storage in an extensive unconfined aquifer system. Results and numerical experiments are presented.Identification is an inverse process whereby the parameters embedded in a differential equation are determined from observations of systems input and output along with appropriate initial and boundary conditions. These parameters are usually noy physically measurable. In general, the governing equation is nonlinear with no closed-form solution.In this paper, this inverse problem is solved by invariant imbedding and quasilinearization. A comparison is made between these methods. The problem of convergence and stability is discussed and demonstrated by numerical experimentation.     

Boundary value problem for a second-order elliptic equation in the exterior of an ellipse

We consider a boundary value problem for a second-order linear elliptic differential equation with constant coefficients in a domain that is the exterior of an ellipse. The boundary conditions of the problem contain the values of the function itself and its normal derivative. We give a constructive solution of the problem and find the number of solvability conditions for the inhomogeneous problem as well as the number of linearly independent solutions of the homogeneous problem. We prove the boundary uniqueness theorem for the solutions of this equation.     

An initial-value method for the determination of Green's functions

The convergence of an initial-value method for computing the Green's function of a class of second-order differential operators is established. The proof relies on an interpolation procedure which is shown to generalize the Nyström method for Fredholm integral equations. The approximate Green's function is related to the solution of a discrete summation equation. The results of Anselone and Moore on collectively compact operators are then applied.This research was partially supported by UNLV Grant No. 001-060-4573.     

Coupling of numerical Green''s matrix and boundary integral equations for the elastodynamic analysis of inhomogeneous bodies on an elastic half-space

A coupled system of integral equations (of the domain and boundary types) is formulated for the elastodynamic response analysis of a locally inhomogeneous body on a homogeneous elastic half-space. The method uses the fundamental solution for homogeneous elastostatics in the inhomogeneous domain owing to the lack of a fundamental solution in inhomogeneous elastodynamics.

The integral representation of displacements in the inhomogeneous domain is formulated with the help of this elastostatic fundamental solution by considering the term induced by the inhomogeneity of materials and the acceleration term as the body force term. Then the Green's matrix is obtained numerically from this integral representation and combined with the ordinary boundary integral equations, which are valid in the exterior homogeneous half-space.

Some numerical examples show the efficiency and the versatility of this coupled method.     

On an integral representation of the solution of the Laplace equation with mixed boundary conditions

We obtain an integral representation of the solution of the Laplace equation with three distinct boundary conditions. Depending on the statement of the problem, the homogeneous boundary value problem may have nontrivial solutions; in other cases, the solution of the homogeneous problem is zero. Note that the inhomogeneous problem is always solvable.     

Invariant imbedding approach for the solution of the minimal surface equation

A new combined technique based on the application of a linearization procedure either (i), the combination of Outer- and Picard-approximation or (ii) the combination of Newton- and Picard-approximation, and invariant imbedding is proposed for obtaining a numerical solution of the minimal surface equation. The existence of inverses of certain matrices appearing in the invariant imbedding equations and the stability of the algorithm are investigated. The minimal surface equation under various boundary conditions and the subsonic fluid flow problem are chosen as test examples for illustrating the method. The numerical results indicate that the proposed method can be used efficiently for solving elliptic problems of a highly nonlinear nature.     

On the solution of the Gellerstedt problem for the Lavrent’ev-Bitsadze equation

We study the solvability of the Gellerstedt problem for the Lavrent??ev-Bitsadze equation under an inhomogeneous boundary condition on the half-circle of the ellipticity domain of the equation, homogeneous boundary conditions on external, internal, and parallel side characteristics of the hyperbolicity domain of the equation, and the transmission conditions on the type change line of the equation.     

Cauchy system for Fredholm resolvents with composite displacement kernels

In recent years, the invariant imbedding approach to initial-value solutions of Fredholm integral equations with degenerate or semidegenerate kernels has been discussed with emphasis. In the present paper, with the aid of invariant imbedding, the Cauchy problem for Fredholm resolvents with composite displacement kernels is reduced to that for generalized Chandrasekhar'sX-functions andY-functions. In other words, it is shown how to convert the initial-value solution of the resolvent into a system of simultaneous nonlinear integrodifferential equations ofX-functions andY-functions of a single argument. From the computational aspect, the result seems to be more tractable than the original ones.Dedicated to R. BellmanThis work was partially supported by the Ministry of Education of Japan, Grant No. 503540. The author is indebted to Professor R. E. Kalaba, University of Southern California, and Dr. H. H. Kagiwada, Hughes Aircraft Company, for their helpful comments and kind interest in the present paper.     

An effective boundary element method for inhomogeneous partial differential equations

A method for removing the domain or volume integral arising in boundary integral formulations for linear inhomogeneous partial differential equations is presented. The technique removes the integral by considering a particular solution to the homogeneous partial differential equation which approximates the inhomogeneity in terms of radial basis functions. The remainder of the solution will then satisfy a homogeneous partial differential equation and hence lead to an integral equation with only boundary contributions. Some results for the inhomogeneous Poisson equation and for linear elastostatics with known body forces are presented.     

Integrable equations for the partition function of the six-vertex model

The partition function of the six-vertex model with the domain-wall boundary condition is considered in the homogeneous and inhomogeneous cases. The determinant representation allows us to show that the partition function is a solution of the Toda equation in the homogeneous case and a solution of the Hirota equation in the inhomogeneous case. Bibliography: 10 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 245, 1997, pp. 207–215. Translated by N. A. Kitanin.     

Multidimensional Kagiwada-Kalaba method of invariant imbedding with general linear implicit boundary conditions

The Kagiwada-Kalaba method of invariant imbedding for multidimensional systems is first derived for the split linear implicit boundary conditions. The justification for the Kagiwada-Kalaba procedure is explained in terms of the special nature of the split linear implicit boundary conditions. Extension of the Kagiwada-Kalaba method from the split linear implicit boundary conditions to general linear implicit boundary conditions is described.     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.