Comlab

Partial differential equations with a nonlinear pointwise constraint defined through a manifold occur in a variety of applications: The magnetization of a ferromagnet can be described by a unit length vector field and the orientation of the rod-like molecules that constitute a liquid crystal is often modeled by a vector field that attains its values in the real projective plane thus respecting the head-to-tail symmetry of the molecules. Other applications arise in geometric

modeling, quantum mechanics, and general relativity. Simple examples reveal that it is impossible to satisfy pointwise constraints exactly by lowest order finite elements. For two model problems we discuss the practical realization of the constraint, the efficient solution of the resulting nonlinear systems of equations, and weak accumulation of approximations at exact solutions.

We focus on domain decomposition methods for systems of PDEs (versus scalar PDEs). The Smith factorization (a "pure" algebra tool) is used systematically to derive new domain decompositions methods for symmetric and unsymmetric systems of PDEs: the compressible Euler equations, the Stokes and Oseen (linearized Navier-Stokes) problem. We will focus on the Stokes system. In two dimensions the key idea is the transformation of the Stokes problem into a scalar bi-harmonic problem. We show, how a proposed domain decomposition method for the bi-harmonic problem leads to a domain decomposition method for the Stokes equations which inherits the convergence behavior of the scalar problem. Thus, it is sufficient to study the convergence of the scalar algorithm. The same procedure can also be applied to the three-dimensional Stokes problem.

The butterfly-shaped Lorenz attractor is a fractal set made up of infinitely many periodic orbits. Ever since Lorenz (1963) introduced a system of three simple ordinary differential equations, much of the discussion of his system and its strange attractor has adopted a dynamical point of view. In contrast, we allow time to be a complex variable and look upon such solutions of the Lorenz system as analytic functions. Formal analysis gives the form and coefficients of the complex singularities of the Lorenz system. Very precise (> 500 digits) numerical computations show that the periodic orbits of the Lorenz system have singularities which obey that form exactly or very nearly so. Both formal analysis and numerical computation suggest that the mathematical analysis of the Lorenz system is a problem in analytic function theory. (Joint work with S. Sahutoglu).

This is not intended to be a systematic History, but a selection of highlights, with some digressions, including:

The early days of the Computing Lab;

How the coming of the Computer changed some of the ways we do Computation;

A problem from the Study Groups;

Influence of the computing environment (hardware and software);

Convergence analysis for the heat equation, then and now.

Dirichlet to Neumann maps and their generalizations are exceptionally useful tools in the study of eigenvalue problems for ODEs and PDEs. They also have real physical significance through their occurrence in electrical impedance tomography, with applications to medical imagine, landmine detection and non-destructive testing. This talk will review some of the basic properties of Dirichlet to Neumann maps, some new abstract results which make it easier to use them for a wide variety of models, and some analytical/numerical results which depend on them, including detection and elimination of spectral pollution.