A shape-optimization technique for the capillary surface problem

Summary The problem of the construction of an equilibrium surface taking the surface tension into account leads to Laplace-Young equation which is a nonlinear elliptic free-boundary problem. In contrast to Orr et al. where an iterative technique is used for direct solution of the equation for problems with simple geometry, we propose here an alternative approach based on shape optimization techniques. The shape of the domain of the liquid is varied to attain the optimality condition. Using optimal control theory to derive expressions for the gradient, a numerical scheme is proposed and simple model problems are solved to validate the scheme.     

Tracking poles and representing Hankel operators directly from data

Summary We propose and analyse a method of estimating the poles near the unit circleT of a functionG whose values are given at a grid of points onT: we give an algorithm for performing this estimation and prove a convergence theorem. The method is to identify the phase for an estimate by considering the peaks of the absolute value ofG onT, and then to estimate the modulus by seeking a bestL ² fit toG over a small arc by a first order rational function. These pole estimates lead to the construction of a basis ofL ² which is well suited to the numerical representation of the Hankel operator with symbolG and thereby to the numerical solution of the Nehari problem (computing the bestH , analytic, approximation toG relative to theL norm), as analysed in [HY]. We present the results of numerical tests of these algorithms.Partially supported by grants from the AFOSR and NSF     

Heat “polynomials” for a singular differential operator on (0, ∞)

We generalize the theory of the heat polynomials introduced by P. V. Rosenbloom and D. V. Widder for a more general class of singular differential operator on (0, ). The heat polynomials associated with the Bessel operator and studied by D. T. Haimo appear as a particular case in this paper. In the special cases of second derivative and Bessel operators the heat polynomials are in fact polynomials inx andt, however, this property does not hold in general.Communicated by Tom. H. Koornwinder.     

Truncation error bounds for modified continued fractions with applications to special functions

Much recent work has been done to investigate convergence of modified continued fractions (MCF's), following the proof by Thron and Waadeland [35] in 1980 that a limit-periodic MCFK(a _n , 1;x ₁), with andnth approximant

Continued fractions associated with the Newton-Padé table

Summary We discuss first the block structure of the Newton-Padé table (or, rational interpolation table) corresponding to the double sequence of rational interpolants for the data{(z _k, h(z_k)} _k =0. (The (m, n)-entry of this table is the rational function of type (m,n) solving the linearized rational interpolation problem on the firstm+n+1 data.) We then construct continued fractions that are associated with either a diagonal or two adjacent diagonals of this Newton-Padé table in such a way that the convergents of the continued fractions are equal to the distinct entries on this diagonal or this pair of diagonals, respectively. The resulting continued fractions are generalizations of Thiele fractions and of Magnus'sP-fractions. A discussion of an some new results on related algorithms of Werner and Graves-Morris and Hopkins are also given.Dedicated to the memory of Helmut Werner (1931–1985)     

Building of mathematical models of continuum media on the basis of the invariance principle

The system of differential equations which describes the motion of continuum media of gas, liquid, Reiner-Rievling-type liquid, etc., is considered.% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacqaHbp% GCdaWgaaWcbaGaamiDaaqabaGccqGHRaWkcaqGKbGaaeyAaiaabAha% caqGOaGaeqyWdiNaaeyDaiaabMcacaqG9aGaaeimaiaabUdaaeaacq% aHbpGCcaGGBbGaaeyDamaaBaaaleaacaqG0baabeaakiabgUcaRiaa% cIcacaqG1bGaeyyXICTaey4bIeTaaiykaiaabwhacaGGDbGaeyOeI0% IaamizaiaadMgacaWG2bGaey4dIuTaaiikaiabgEGirlaabwhacaGG% PaGaey4kaSIaey4bIeTaamiCaiaacUdaaeaacaWGWbWaaSbaaSqaai% aadshaaeqaaOGaey4kaSIaaeyDaiabgwSixlabgEGirlaadchacqGH% RaWkcaWGhbGaaeizaiaabMgacaqG2bGaaeiiaiaabwhacqGHRaWkca% WGibGaeqOXdyMaeyypa0JaaGimaiaac6caaaaa!7268!\[\begin{gathered} \rho _t + {\text{div(}}\rho {\text{u) = 0;}} \hfill \\ \rho [{\text{u}}_{\text{t}} + ({\text{u}} \cdot \nabla ){\text{u}}] - div\prod (\nabla {\text{u}}) + \nabla p; \hfill \\ p_t + {\text{u}} \cdot \nabla p + G{\text{div u}} + H\phi = 0. \hfill \\ \end{gathered} \]Solving the problem of its group classification, we obtained all the state equations which lead to the expansion of the main group assumed by the initial equations under the arbitrary elements , G, H.     

Adaptive boundary element methods for strongly elliptic integral equations

Summary Integral operators are nonlocal operators. The operators defined in boundary integral equations to elliptic boundary value problems, however, are pseudo-differential operators on the boundary and, therefore, provide additional pseudolocal properties. These allow the successful application of adaptive procedures to some boundary element methods. In this paper we analyze these methods for general strongly elliptic integral equations and obtain a-posteriori error estimates for boundary element solutions. We also apply these methods to nodal collocation with odd degree splines. Some numerical examples show that these adaptive procedures are reliable and effective.This work was carried out while Dr. De-hao Yu was an Alexander-von-Humboldt-Stiftung research fellow at the University of Stuttgart in 1987, 1988     

On an application of Hall's representatives theorem to a finite geometry problem

It is well-known thatn points not belonging to a hyperplane determine at leastn hyperplanes. The possible configurations of hyperplanes in the case when the number of hyperplanes is equal ton are known, too. In this paper we obtain these results by means of Hall's representatives theorem. The setting is that of a finite geometry.     

Differential form methods in the theory of variational systems and Lagrangian field theories

This paper will attempt to unify diverse material from physics and engineering in terms of differential forms on manifolds. A variational system will be defined by means of a scalar-valued differential form on a manifold and an ideal in the Grassmann algebra of differential forms on that manifold to serve as constraints. Two types of extremal submanifolds will be defined. The first-called the Euler-Lagrange extremals-will be defined by a method that is the generalization of the classical methods in the calculus of variations. The second—a generalization of a method used by Cartan in his treatise Leçons sur les invariants intégraux-will define extremals as integral submanifolds of an exterior differential system invariently attached to the variational system. As examples, the variational systems attached to string theories in Riemannian manifolds and Yang-Mills fields will be discussed from this differential form point of view. Finally, as application, the differential geometric properties and definition of energy will be presented from the differential form point of view.This work was supported by a grant from the Applied Mathematics program of the National Science Foundation.