We develop a theory to connect the following three areas: (a) the mean field equations on flat tori, (b) the classical Lame equations and (c) modular forms. A major theme in part I is a classification of developing maps f attached to solutions of the mean field equation according to the types of transformation laws (or monodromy) with respect to the period lattice satisfied by f.

As V. I. Arnold observed in the 1960s, the Euler equations of incompressible fluid flow correspond formally to geodesic equations in a group of volume-preserving diffeomorphisms. Working in an Eulerian framework, we study incompressible flows of shapes as critical paths for action (kinetic energy) along transport paths constrained to be shape densities (characteristic functions). The formal geodesic equations for this problem are Euler equations for incompressible, inviscid potential flow of fluid with zero pressure and surface tension on the free boundary. The problem of minimizing this action exhibits an instability associated with microdroplet formation, with the following outcomes: Any two shapes of equal volume can be approximately connected by an Euler spray---a countable superposition of ellipsoidal geodesics. The infimum of the action is the Wasserstein distance squared, and is almost never attained except in dimension 1. Every Wasserstein geodesic between bounded densities of compact support provides a solution of the (compressible) pressureless Euler system that is a weak limit of (incompressible) Euler sprays. Each such Wasserstein geodesic is also the unique minimizer of a relaxed least-action principle for a two-fluid mixture theory corresponding to incompressible fluid mixed with vacuum.

Using the convex integration technique for the three-dimensional Navier-Stokes equations introduced by T. Buckmaster and V. Vicol, it is shown the existence of non-unique weak solutions for the 3D Navier-Stokes equations with fractional hyperviscosity $(-\Delta)^{\theta}$, whenever the exponent $\theta$ is less than J.-L. Lions' exponent $5/4$, i.e., when $\theta < 5/4$.

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

No abstract uploaded!