A boundary expansion of solutions to nonlinear singular elliptic equations

In this paper we establish an asymptotic expansion near the boundary for solutions to the Dirichlet problem of elliptic equations with singularities near the boundary.This expansion formula shows the singularity profile of solutions at the boundary.We deal with both linear and nonlinear elliptic equations,including fully nonlinear elliptic equations and equations of Monge-Ampère type.     

Singularities of divisors of low degree on abelian varieties

Building on previous work of Kollár, Ein, Lazarsfeld, and Hacon, we show that ample divisors of low degree on an abelian variety have mild singularities in case the abelian variety is simple or the degree of the polarization is two. Olivier Debarre was visiting the University of Michigan when part of this work was done, with support from William Fulton and Robert Lazarsfeld. Christopher Hacon was partially supported by NSA research grant no: MDA904-03-1-0101 and by a grant from the Sloan Foundation.     

ASYMPTOTIC SOLUTIONS OF THE ELASTIC EQUATION IN THE CASE OF TOTAL REFLECTION

In this paper we construct outgoing asymptotic solutions of the anisotropic elastic wave equation coinciding with given data on the boundary when total reflection happens. In this case the solutions consist of the usual hyperbolic parts and the elliptic parts exponentially decreasing in the normal direction to the boundary. The main assertion is that sum of those parts can cover any boundary data because of the elasticity. The proof is reduced to the complex analysis of the symbol of the elastic operator. By means of the asymptotic solutions we can make parametrices and study propagation of the singularities in the case of total reflection.

     

Steady-state vibrations of an unbounded linear piezoelectric medium

We consider fundamental (Dirichlet and Neumann-type) boundary value problems in a theory of generalized plane strain for the steady-state vibrations of an infinite piezoelectric medium with transversely isotropic symmetry (6 mm). Using integral equation methods with the appropriate Sommerfeld-type radiation conditions, we prove existence and uniqueness results for the corresponding exterior boundary value problems. Exact solutions are obtained in the form of integral potentials. (Received: September 27, 2005)     

Multiscale asymptotic expansion and finite element methods for the mixed boundary value problems of second order elliptic equation in perforated domains

In this paper, we will discuss the mixed boundary value problems for the second order elliptic equation with rapidly oscillating coefficients in perforated domains, and will present the higher-order multiscale asymptotic expansion of the solution for the problem, which will play an important role in the numerical computation . The convergence theorems and their rigorous proofs will be given. Finally a multiscale finite element method and some numerical results will be presented. This work is Supported by National Natural Science Foundation of China (grant # 10372108, # 90405016), and Special Funds for Major State Basic Research Projects( grant # TG2000067102)     

Existence and Exponential Decay of Solutions to a Quasilinear Thermoelastic Plate System

We consider a quasilinear PDE system which models nonlinear vibrations of a thermoelastic plate defined on a bounded domain in , n ≤ 3. Existence of finite energy solutions describing the dynamics of a nonlinear thermoelastic plate is established. In addition asymptotic long time behavior of weak solutions is discussed. It is shown that finite energy solutions decay exponentially to zero with the rate depending only on the (finite energy) size of initial conditions. The proofs are based on methods of weak compactness along with nonlocal partial differential operator multipliers which supply the sought after “recovery” inequalities. Regularity of solutions is also discussed by exploiting the underlying analyticity of the linearized semigroup along with a related maximal parabolic regularity [1, 16, 44]. The research of I. Lasiecka has been partially supported by DMS-NSF Grant Nr 0606882. S. Maad was supported by the Swedish Research Council and by the European Union under the Marie Curie Fellowship MEIF-CT-2005-024191.     

Solutions for the prescribing mean curvature equation

By variational methods, for a kind of Yamabe problem whose scalar curvature vanishes in the unit ball BN and on the boundary S^N-1 the mean curvature is prescribed, we construct multi-peak solutions whose maxima are located on the boundary as the parameter tends to 0^＋ under certain assumptions. We also obtain the asymptotic behaviors of the solutions.     

Determination of jumps of distributions by differentiated means

Differentiated means are defined in order to find formulas for jumps of distributions. We analyze two types of jumps occurring in the notions of distributional jump behavior and symmetric jump behavior. We start by defining what we call Riesz differentiated means for numerical series, then the differentiated means are extended to distributional evaluations for the Schwartz class of tempered distributions. The jumps of tempered distributions are completely determined by the differentiated means of the Fourier transform. We also find formulas for the jumps in terms of the asymptotic behavior of partial derivatives of harmonic representations and harmonic conjugate functions. Applications to Fourier series are given. The second author gratefully acknowledges support by the Louisiana State Board of Regents grant LEQSF(2005-2007)-ENH-TR-21.     

Thermal effects in orthotropic porous elastic beams

This paper is concerned with the linear theory of anisotropic porous elastic bodies. The extension and bending of orthotropic porous elastic cylinders subjected to a plane temperature field is investigated. The work is motivated by the recent interest in the using of the orthotropic porous elastic solid as model for bones and various engineering materials. First, the thermoelastic deformation of inhomogeneous beams whose constitutive coefficients are independent of the axial coordinate is studied. Then, the extension and bending effects in orthotropic cylinders reinforced by longitudinal rods are investigated. The three-dimensional problem is reduced to the study of two-dimensional problems. The method is used to solve the problem of an orthotropic porous circular cylinder with a special kind of inhomogeneity.      

Iterative solution of a singular convection-diffusion perturbation problem

An iterative domain decomposition method is developed to solve a singular perturbation problem. The problem consists of a convection-diffusion equation with a discontinuous (piecewise-constant) diffusion coefficient, and the problem domain is decomposed into two subdomains, on each of which the coefficient is constant. After showing that the boundary value problem is well posed, we indicate a specific numerical implementation of the iterative technique that combines the finite element method on one subdomain with the method of matched asymptotic expansions on the other subdomain. This procedure extends work by Carlenzoli and Quarteroni, which was originally intended for a boundary layer problem with an outer region and an inner region. Our extension carries over to a problem where the domain consists of the outer and inner boundary layer regions plus a region in which the diffusion coefficient is constant and significant in magnitude. An unexpected benefit of our new implementation is its efficiency, which is due to the fact that at each iteration the problem needs to be solved explicitly only on one subdomain. It is only when the final approximation on the entire domain is desired that the matched asymptotic expansions approximation need be computed on the second subdomain. Two-dimensional convergence results and numerical results illustrating the method for a two-dimensional test problem are given.Received: February 12, 2004     

Performance Study and Analysis of Parallel Multilevel Preconditioners

The significant gap between peak and realized performance of parallel systems motivates the need for performance analysis. In order to predict the performance of a class of parallel multilevel ILU preconditioner (PBILUM), we build two performance prediction models for both the preconditioner construction phase and the solution phase. These models combine theoretical features of the preconditioners with estimates on computation cost, communications overhead, etc. Experimental simulations show that our model predication based on certain reasonable assumptions is close to the simulation results. The models may be used to predict the performance of this class of parallel preconditioners.*The research work of the authors was supported in part by the U.S. National Science Foundation under grants CCR-9988165, CCR-0092532, ACR-0202934, and ACR-234270, by the U.S. Department of Energy Office of Science under grant DE-FG02-02ER45961, by the Kentucky Science & Engineering Foundation under grant KSEF-02-264-RED-002.     

Global stability of the Armstrong-Frederick model with periodic biaxial inputs

The paper is concerned with the study of plasticity models described by differential equations with stop and play operators. We suggest sufficient conditions for the global stability of a unique periodic solution for the scalar models and for the vector models with biaxial inputs of a particular form, namely the sum of a uniaxial function and a constant term. For another class of simple biaxial inputs, we present an example of the existence of unstable periodic solutions. The paper was written during the research stay of D. Rachinskii at the Technical University Munich supported by the research fellowship from the Alexander von Humboldt Foundation. His work was partially supported by the Russian Science Support Foundation, Russian Foundation for Basic Research (Grant No. 01-01-00146, 03-01-00258), and the Grants of the President of Russia (Grant No. MD-87.2003.01, NS-1532.2003.1). The support is gratefully acknowledged.     

Law of large numbers and central limit theorem for Donsker's delta function of diffusions I

Limit theorems in the space of Hida distributions, similar to the law of large numbers and the central limit theorem, are shown for composites of the Dirac distribution with solutions of one-dimensional, non-degenerate Itô equations.Supported by National Science Foundation under grant DMS-9001859.Supported by the Louisiana Education Quality Support Fund under grant (91–93) RD-A-08.Supported by the Council on Research of Louisiana State University.     

Decay estimates for a thermoelastic system in waveguides

We investigate the linear system of thermoelasticity, consisting of an elasticity equation and a heat conduction equation, in a waveguide Ω=(0,1)×Rn−1, with certain boundary conditions. We consider the cases of homogeneous and inhomogeneous systems and prove decay estimates of the solutions, which are a key ingredient to showing the global existence of solutions to non-linear thermoelasticity, after having decomposed the solutions into various parts. We also give a simplified proof to the representation of the solutions to the Cauchy problem of thermoelasticity.     

Three-scale asymptotics for a diffusion problem coupled with the wave equation

In this article, we present an asymptotic analysis of waves of elastic stress in an infinite solid whose boundary is subject to a rapid thermal load. The problem under consideration couples the wave equation and the heat equation, and the asymptotic approximation of the solution requires three-scaled variables. The asymptotic approximation is supplied with a rigorous remainder estimate and is illustrated numerically.     

Estimating the largest singular values of large sparse matrices via modified moments

We describe a procedure for determining a few of the largest singular values of a large sparse matrix. The method by Golub and Kent which uses the method of modified moments for estimating the eigenvalues of operators used in iterative methods for the solution of linear systems of equations is appropriately modified in order to generate a sequence of bidiagonal matrices whose singular values approximate those of the original sparse matrix. A simple Lanczos recursion is proposed for determining the corresponding left and right singular vectors. The potential asynchronous computation of the bidiagonal matrices using modified moments with the iterations of an adapted Chebyshev semi-iterative (CSI) method is an attractive feature for parallel computers. Comparisons in efficiency and accuracy with an appropriate Lanczos algorithm (with selective re-orthogonalization) are presented on large sparse (rectangular) matrices arising from applications such as information retrieval and seismic reflection tomography. This procedure is essentially motivated by the theory of moments and Gauss quadrature.This author's work was supported by the National Science Foundation under grants NSF CCR-8717492 and CCR-910000N (NCSA), the U.S. Department of Energy under grant DOE DE-FG02-85ER25001, and the Air Force Office of Scientific Research under grant AFOSR-90-0044 while at the University of Illinois at Urbana-Champaign Center for Supercomputing Research and Development.This author's work was supported by the U.S. Army Research Office under grant DAAL03-90-G-0105, and the National Science Foundation under grant NSF DCR-8412314.     

A bifurcation analysis of the nonlinear parametric programming problem

The structure of solutions to the nonlinear parametric programming problem with a one dimensional parameter is analyzed in terms of the bifurcation behavior of the curves of critical points and the persistence of minima along these curves. Changes in the structure of the solution occur at singularities of a nonlinear system of equations motivated by the Fritz John first-order necessary conditions. It has been shown that these singularities may be completely partitioned into seven distinct classes based upon the violation of one or more of the following: a complementarity condition, a constraint qualification, and the nonsingularity of the Hessian of the Lagrangian on a tangent space. To apply classical bifurcation techniques to these singularities, a further subdivision of each case is necessary. The structure of curves of critical points near singularities of lowest (zero) codimension within each case is analyzed, as well as the persistence of minima along curves emanating from these singularities. Bifurcation behavior is also investigated or discussed for many of the subcases giving rise to a codimension one singularity.This work was supported by the National Science Foundation through NSF Grants DMS-85-10201 and DMS-87-04679 and by the Air Force Office of Scientific Research through grant number AFOSR-88-0059.     

A generalized Green's formula for elliptic problems in domains with edges

The usual Green's formula connected with the operator of a boundary-value problem fails when both of the solutions u and v that occur in it have singularities that are too strong at a conic point or at an edge on the boundary of the domain. We deduce a generalized Green's formula that acquires an additional bilinear form in u and v and is determined by the coefficients in the expansion of solutions near singularities of the boundary. We obtain improved asymptotic representations of solutions in a neighborhood of an edge of positive dimension, which together with the generalized Green's formula makes it possible, for example, to describe the infinite-dimensional kernel of the operator of an elliptic problem in a domain with edge. Bibliography: 14 titles. Translated fromProblemy Matematicheskogo Analiza, No. 13, 1992, pp. 106–147.     

Perturbation results for nearly uncoupled Markov chains with applications to iterative methods

Summary The standard perturbation theory for linear equations states that nearly uncoupled Markov chains (NUMCs) are very sensitive to small changes in the elements. Indeed, some algorithms, such as standard Gaussian elimination, will obtain poor results for such problems. A structured perturbation theory is given that shows that NUMCs usually lead to well conditioned problems. It is shown that with appropriate stopping, criteria, iterative aggregation/disaggregation algorithms will achieve these structured error bounds. A variant of Gaussian elimination due to Grassman, Taksar and Heyman was recently shown by O'Cinneide to achieve such bounds.Supported by the National Science Foundation under grant CCR-9000526 and its renewal, grant CCR-9201692. This research was done in part, during the author's visit to the Institute for Mathematics and its Applications, 514 Vincent Hall, 206 Church St. S.E., University of Minnesota, Minneapolis, MN 55455, USA     

The central limit theorem for random perturbations of rotations

Summary.   We prove a functional central limit theorem for stationary random sequences given by the transformations