Past OxPDE Lunchtime Seminar

The aim of this talk is to explain how to construct solutions to a

relativistic transport equation via a time discrete scheme based on an

optimal transportation problem.

First of all, I will present a joint work with J. Bertrand, where we prove the existence of an optimal map

for the Monge-Kantorovich problem associated to relativistic cost functions.

Then, I will explain a joint work with Robert McCann, where

we study the limiting process between the discrete and the continuous

equation.

I will talk about $W^{2,1}$ regularity for strictly convex Aleksandrov solutions to the Monge Amp\`ere equation

\[

\det D^2 u =f

\]

where $f$ satisfies $\log f\in L^{\infty} $. Under the previous assumptions in the 90's Caffarelli was able to prove that $u \in C^{1,\alpha}$ and that $u\in W^{2,p}$ if $|f-1|\leq \varepsilon(p)$. His results however left open the question of Sobolev regularity of $u$ in the general case in which $f$ is just bounded away from $0$ and infinity. In a joint work with Alessio Figalli we finally show that actually $|D^2u| \log^k |D^2 u| \in L^1$ for every positive $k$.

\\

If time will permit I will also discuss some question related to the $W^{2,1}$ stability of solutions of Monge-Amp\`ere equation and optimal transport maps and some applications of the regularity to the study of the semi-geostrophic system, a simple model of large scale atmosphere/ocean flows (joint works with Luigi Ambrosio, Maria Colombo and Alessio Figalli).

In the first part of the talk I will introduce a notion of convexity ($\mathcal{X}$-convexity) which applies to any given family of vector fields: the main model which we have in mind is the case of vector fields satisfying the H\"ormander condition.

Then I will give a PDE-characterization for $\mathcal{X}$-convex functions using a viscosity inequality for the intrinsic Hessian and I will derive bounds for the intrinsic gradient and intrinsic local Lipschitz-continuity for this class of functions.\\

In the second part of the talk I will introduce a notion of subdifferential for any given family of vector fields (namely $\mathcal{X}$-subdifferential) and show that a non empty $\mathcal{X}$-subdifferential at any point characterizes the class of $\mathcal{X}$-convex functions.

As application I will prove a Jensen-type inequality for $\mathcal{X}$-convex functions in the case of Carnot-type vector fields. {\em (Joint work with Martino Bardi)}.

An overview is given of some key issues and definitions in the Calculus of Variations, with a focus on lower semicontinuity and quasiconvexity. Some well known results and instructive counterexamples are also discussed. We then move to consider variational problems in the BV setting, and present a new lower semicontinuity result for quasiconvex integrals of subquadratic growth. The proof of this requires some interesting techniques, such as obtaining boundedness properties for an extension operator, and exploiting fine properties of Sobolev maps.

We will give a survey on some recent results on travel tomography which consists in determining the index of refraction of a medium by measuring the travel times of sound waves going through the medium. In differential geometry this is known as the boundary rigidity problem. We will also consider the related problem of tensor tomography which consists in determining a function, a vector field or tensors of higher rank from their integrals along geodesics.

We consider an active scalar equation with singular drift velocity that is motivated by a model for the geodynamo. We show that the non-diffusive equation is ill-posed in the sense of Hadamard in Sobolev spaces. In contrast, the critically diffusive equation is globally well-posed. This work is joint with Vlad Vicol.

Discrete problems have a habit of being beautiful but difficult. This can be true even of discrete problems whose continuous analogues are easy. For example: computing the surface area of a sphere of radius N^{1/2} in k-dimensional Euclidean space (easy). Counting the number of representations of an integer N as a sum of k squares (historically hard). In this talk we'll survey a menagerie of discrete analogues of operators arising in harmonic analysis, including singular integral operators (such as the Hilbert transform), maximal functions, and fractional integral operators. In certain cases we can learn everything we want to know about the discrete operator immediately, from its continuous analogue. In other cases the discrete operator requires a completely new approach. We'll see what makes a discrete operator easy/hard to treat, and outline some of the methods that are breaking new ground, key aspects of which come from number theory. In particular, we will highlight the roles played by theta functions, exponential sums, Waring's problem, and the circle method of Hardy and Littlewood. No previous knowledge of singular integral operators or the circle method will be assumed.

Solutions which are time-bounded in L^3 up to time T can be continued

past this time, by a landmark result of Escauriaza-Seregin-Sverak,

extending Serrin's criterion. On the other hand, the local Cauchy

theory holds up to solutions in BMO^-1; we aim at describing how one

can obtain intermediate regularity results, assuming a priori bounds

in negative regularity Besov spaces.

This is joint work with J.-Y. Chemin, Isabelle Gallagher and Gabriel

Koch.

We prove the existence of weak solutions to steady Navier Stokes equations

$$\text{div}\, \sigma+f=\nabla\pi+(\nabla u)u.$$

Here $u:\mathbb{R}^2\supset \Omega\rightarrow \mathbb{R}^2$ denotes

the velocity field satisfying $\text{div}\, u=0$,

$f:\Omega\rightarrow\mathbb{R}^2$ and

$\pi:\Omega\rightarrow\mathbb{R}$ are external volume force and

pressure, respectively. In order to model the behavior of

Prandtl-Eyring fluids we assume

$$\sigma= DW(\varepsilon (u)),\quad W(\varepsilon)=|\varepsilon|\log

(1+|\varepsilon|).$$

A crucial tool in our approach is a modified Lipschitz truncation

preserving the divergence of a given function.

I will describe recent work with Charles Fefferman on a

construction of families of analytic almost-sharp fronts for SQG. These

are special solutions of SQG which have a very sharp transition in a

very thin layer. One of the main difficulties of the construction is the

fact that there is no formal limit for the family of equations. I will

show how to overcome this difficulty, linking the result to joint work

with C. Fefferman and Kevin Luli on the existence of a "spine" for

almost-sharp fronts. This is a curve, defined for every time slice by a

measure-theoretic construction, that describes the evolution of the

almost-sharp front.

Modern mathematical approaches to plasticity result in non-convex variational problems for which the standard methods of the calculus of variations are not applicable. In this contribution we consider geometrically nonlinear crystal elasto-plasticity in two dimensions with one active slip system. In order to derive information about macroscopic material behavior the relaxation of the corresponding incremental problems is studied. We focus on the question if realistic systems with an elastic energy leading to large penalization of small elastic strains can be well-approximated by models based on the assumption of rigid elasticity. The interesting finding is that there are qualitatively different answers depending on whether hardening is included or not. In presence of hardening we obtain a positive result, which is mathematically backed up by Γ-convergence, while the material shows very soft macroscopic behavior in case of no hardening. The latter is due to the vanishing relaxation for a large class of applied loads.

This is joint work with Sergio Conti and Georg Dolzmann.

We propose a general framework for the study of $L^1$ contractive semigroups of solutions to conservation laws with discontinuous flux. Developing the ideas of a number of preceding works we claim that the whole admissibility issue is reduced to the selection of a family of "elementary solutions", which are certain piecewise constant stationary weak solutions. We refer to such a family as a "germ". It is well known that (CL) admits many different $L^1$ contractive semigroups, some of which reflects different physical applications. We revisit a number of the existing admissibility (or entropy) conditions and identify the germs that underly these conditions. We devote specific attention to the anishing viscosity" germ, which is a way to express the "$\Gamma$-condition" of Diehl. For any given germ, we formulate "germ-based" admissibility conditions in the form of a trace condition on the flux discontinuity line $x=0$ (in the spirit of Vol'pert) and in the form of a family of global entropy inequalities (following Kruzhkov and Carrillo). We characterize those germs that lead to the $L^1$-contraction property for the associated admissible solutions. Our approach offers a streamlined and unifying perspective on many of the known entropy conditions, making it possible to recover earlier uniqueness results under weaker conditions than before, and to provide new results for other less studied problems. Several strategies for proving the existence of admissible solutions are discussed, and existence results are given for fluxes satisfying some additional conditions. These are based on convergence results either for the vanishing viscosity method (with standard viscosity or with specific viscosities "adapted" to the choice of a germ), or for specific germ-adapted finite volume schemes.

This is joint work with Boris Andreianov and Nils Henrik Risebro.

The reconstruction of the early universe amounts to recovering the tiny density fluctuations of the early universe (shortly after the "big bang") from the current observation of the matter distribution in the universe. Following Zeldovich, Peebles and, more recently Frisch and collaboratoirs, we use a newtonian gravitational model with time dependent coefficients taking into accont general relativity effects. Due to the (remarkable) convexity of the corresponding action, the reconstruction problem apparently reduces to a straightforward convex minimization problem. Unfortunately, this approach completely ignores the mass concentration effects due to gravitational instabilities.

In this lecture, we show a way of modifying the action in order to take concentrations into account. This is obtained through a (questionable) modification of the gravitation model,

by substituting the fully nonlinear Monge-Amp`ere equation for the linear Poisson equation. (This is a reasonable approximation in the sense that it makes exact some approximate solutions advocated by Zeldovich for the original gravitational model.) Then the action can be written as a perfect square in which we can input mass concentration effects in a canonical way, based on the theory of gradient flows with convex potentials and somewhat related to the concept of self-dual Lagrangians developped by Ghoussoub. A fully discrete algorithm is introduced for the EUR problem in one space dimension.

Considering kinetic equations (Boltzmann, BGK, say...) in the small mean free path regime lead to conservation laws (the Euler system, typically) When the problem is set in a domain, boundary layers might occur due to the fact that incoming fluxes could be far from equilibrium states. We consider the problem from a numerical perspective and we propose a definition of numerical fluxes for the Euler system which is intended to account for the formation of these boundary layers.

Free-discontinuity problems describe situations where the solution of

interest is defined by a function and a lower dimensional set consisting

of the discontinuities of the function. Hence, the derivative of the

solution is assumed to be a "small function" almost everywhere except on

sets where it concentrates as a singular measure.

This is the case, for instance, in certain digital image segmentation

problems and brittle fracture models.

In the first part of this talk we show new preliminary results on

the existence of minimizers for inverse free-discontinuity problems, by

restricting the solutions to a class of functions with piecewise Lipschitz

discontinuity set.

If we discretize such situations for numerical purposes, the inverse

free-discontinuity problem in the discrete setting can be re-formulated as

that of finding a derivative vector with small components at all but a few

entries that exceed a certain threshold. This problem is similar to those

encountered in the field of "sparse recovery", where vectors

with a small number of dominating components in absolute value are

recovered from a few given linear measurements via the minimization of

related energy functionals.

As a second result, we show that the computation of global minimizers in

the discrete setting is an NP-hard problem.

With the aim of formulating efficient computational approaches in such

a complicated situation, we address iterative thresholding algorithms that

intertwine gradient-type iterations with thresholding steps which were

designed to recover sparse solutions.

It is natural to wonder how such algorithms can be used towards solving

discrete free-discontinuity problems. This talk explores also this

connection, and, by establishing an iterative thresholding algorithm for

discrete inverse free-discontinuity problems, provides new insights on

properties of minimizing solutions thereof.

Tonelli gave the first rigorous treatment of one-dimensional variational problems, providing conditions for existence and regularity of minimizers over the space of absolutely continuous functions.  He also proved a partial regularity theorem, asserting that a minimizer is everywhere differentiable, possible with infinite derivative, and that this derivative is continuous as a map into the extended real line.  Some recent work has lowered the smoothness assumptions on the Lagrangian for this result to various Lispschitz and H\"older conditions.  In this talk we will discuss the partial regularity result, construct examples showing that mere continuity of the Lagrangian is an insufficient condition.

It has long been known that many materials are crystalline when in their energy-minimizing states. Two of the most common crystalline structures are the face-centred cubic (fcc) and hexagonal close-packed (hcp) crystal lattices. Here we introduce the problem of crystallization from a mathematical viewpoint and present an outline of a proof that the ground state of a large system of identical particles, interacting under a suitable potential, behaves asymptotically like fcc or hcp, as the number of particles tends to infinity. An interesting feature of this result is that it holds under no initial assumption on the particle positions. The talk is based upon a joint work in progress with Florian Theil.

It is well-known that Morrey's quasiconvexity is closely related to gradient Young measures,

i.e., Young measures generated by sequences of gradients in

$L^p(\Omega;\mathbb{R}^{m\times n})$. Concentration effects,

however, cannot be treated by Young measures. One way how to describe both oscillation and

concentration effects in a fair generality are the so-called DiPerna-Majda measures.

DiPerna and Majda showed that having a sequence $\{y_k\}$ bounded in $L^p(\Omega;\mathbb{R}^{m\times n})$,$1\le p$ $0$.

We consider thin films of a cholesteric liquid-crystal material subject to an applied electric field.  In such materials, the liquid-crystal "director" (local average orientation of the long axis of the molecules) has an intrinsic tendency to rotate in space; while the substrates that confine the film tend to coerce a uniform orientation.

The electric field encourages certain preferred orientations of the director as well, and these competing influences give rise to several different stable equilibrium states of the director field, including spatially uniform, translation invariant (functions only of position across the cell gap) and periodic (with 1-D or 2-D periodicity in the plane of the film).  These structures depend on two principal control parameters: the ratio of the cell gap to the intrinsic "pitch" (spatial period of rotation) of the cholesteric and the magnitude of the applied voltage.

We report on numerical work (not complete) on the bifurcation and phase behavior of this system.  The study was motivated by potential applications involving switchable gratings and eyewear with tunable transparency. We compare our results with experiments conducted in the Liquid Crystal Institute at Kent State University.

Numerical homogenization/upscaling for problems with multiple scales have attracted increasing attention in recent years. In particular, problems with non-separable scales pose a great challenge to mathematical analysis and simulation.

In this talk, we present some rigorous results on homogenization of divergence form scalar and vectorial elliptic equations with $L^\infty$ rough coefficients which allow for a continuum of scales. The first approach is based on a new type of compensation phenomena for scalar elliptic equations using the so-called ``harmonic coordinates''. The second approach, the so-called ``flux norm approach'' can be applied to finite dimensional homogenization approximations of both scalar and vectorial problems with non-separated scales. It can be shown that in the flux norm, the error associated with approximating the set of solutions of the PDEs with rough coefficients, in a properly defined finite-dimensional basis, is equal to the error associated with approximating the set of solutions of the same type of PDEs with smooth coefficients in a standard finite element space. We will also talk about the ongoing work on the localization of the basis in the flux norm approach.

{\bf Keble Workshop on Partial Differential Equations

in Science and Engineering}

\\

\\Place: Roy Griffiths Room in the ARCO Building, Keble College

\\Time: 1:00pm-5:10pm, Thursday, June 10.

\\

\\

Program:\\

\\ 1:00-1:20pm: Coffee and Tea

\\

\\ 1:20-2:10pm: Prof. Walter Craig (Joint with OxPDE Lunchtime Seminar)

\\

\\ 2:20-2:40pm Prof. Mikhail Feldman

\\

\\ 2:50-3:10pm Prof. Paul Taylor

\\

\\ 3:20-3:40pm Coffee and Biscuits

\\

\\ 3:40-4:00pm: Prof. Sir John Ball

\\

\\ 4:10-4:30pm: Dr. Apala Majumdar

\\

\\ 4:40-5:00pm: Prof. Robert Pego

\\

\\ 5:10-6:00pm: Free Discussion

\\

\\{\bf Titles and Abstracts:}

\\

1.{\bf Title: On the singular set of the Navier-Stokes equations

\\ Speaker: Prof. Walter Craig, McMaster University, Canada}

\\ Abstract:\\

The Navier-Stokes equations are important in

fluid dynamics, and a famous mathematics problem is the

question as to whether solutions can form singularities.

I will describe these equations and this problem, present

three inequalities that have some implications as to the

question of singularity formation, and finally, give a

new result which is a lower bound on the size of the

singular set, if indeed singularities exist.

\\

\\{\bf 2. Title: Shock Analysis and Nonlinear Partial Differential Equations of Mixed Type.

\\ Speaker: Prof. Mikhail Feldman, University of Wisconsin-Madison, USA}

\\

\\ Abstract:\\ Shocks in gas or compressible fluid arise in various physical

situations, and often exhibit complex structures. One example is reflection

of shock by a wedge. The complexity of reflection-diffraction configurations

was first described by Ernst Mach in 1878. In later works, experimental and

computational studies and asymptotic analysis have shown that various patterns

of reflected shocks may occur, including regular and Mach reflection. However,

many fundamental issues related to shock reflection are not understood,

including transition between different reflection patterns. For this reason

it is important to establish mathematical theory of shock reflection,

in particular existence and stability of regular reflection solutions for PDEs

of gas dynamics. Some results in this direction were obtained recently.

\\

In this talk we start by discussing examples of shocks in supersonic and

transonic flows of gas. Then we introduce the basic equations of gas dynamics:

steady and self-similar compressible Euler system and potential flow equation.

These equations are of mixed elliptic-hyperbolic type. Subsonic and supersonic

regions in the flow correspond to elliptic and hyperbolic regions of solutions.

Shocks correspond to certain discontinuities in the solutions. We discuss some

results on existence and stability of steady and self-similar shock solutions,

in particular the recent work (joint with G.-Q. Chen) on global existence of

regular reflection solutions for potential flow. We also discuss open problems

in the area.

\\

\\{\bf 3. Title: Shallow water waves - a rich source of interesting solitary wave

solutions to PDEs

\\ Speaker: Prof. Paul H. Taylor, Keble College and Department of Engineering Science, Oxford}

\\

\\Abstract:\\ In shallow water, solitary waves are ubiquitous: even the wave crests

in a train of regular waves can be modelled as a succession of solitary waves.

When successive crests are of different size, they interact when the large ones

catch up with the smaller. Then what happens? John Scott Russell knew by experiment

in 1844, but answering this question mathematically took 120 years!

This talk will examine solitary wave interactions in a range of PDEs, starting

with the earliest from Korteweg and De Vries, then moving onto Peregrine's

regularized long wave equation and finally the recently introduced Camassa-Holm

equation, where solitary waves can be cartoon-like with sharp corners at the crests.

For each case the interactions can be described using the conserved quantities,

in two cases remarkably accurately and in the third exactly, without actually

solving any of the PDEs.

The methodology can be extended to other equations such as the various versions

of the Boussinesq equations popular with coastal engineers, and perhaps even

the full Euler equations.

\\

{\bf 4. Title: Austenite-Martensite interfaces

\\ Speaker: Prof. Sir John Ball, Queen's College and Mathematical Institute, Oxford}

\\

\\Abstract:\\ Many alloys undergo martensitic phase transformations

in which the underlying crystal lattice undergoes a change of shape

at a critical temperature. Usually the high temperature phase (austenite)

has higher symmetry than the low temperature phase (martensite).

In order to nucleate the martensite it has to somehow fit geometrically

to the austenite. The talk will describe different ways in which this

occurs and how they may be studied using nonlinear elasticity and

Young measures.

\\

\\{\bf 5. Title: Partial Differential Equations in Liquid Crystal Science and

Industrial Applications

\\ Speaker: Dr. Apala Majumdar, Keble College and Mathematical Institute, Oxford}

\\

\\Abstract:\\

Recent years have seen a growing demand for liquid crystals in modern

science, industry and nanotechnology. Liquid crystals are mesophases or

intermediate phases of matter between the solid and liquid phases of

matter, with very interesting physical and optical properties.

We briefly review the main mathematical theories for liquid crystals and

discuss their analogies with mathematical theories for other soft-matter

phases such as the Ginzburg-Landau theory for superconductors. The

governing equations for the static and dynamic behaviour are typically

given by systems of coupled elliptic and parabolic partial differential

equations. We then use this mathematical framework to model liquid crystal

devices and demonstrate how mathematical modelling can be used to make

qualitative and quantitative predictions for practical applications in

industry.

\\

\\{\bf 6. Title: Bubble bath, shock waves, and random walks --- Mathematical

models of clustering

\\Speaker: Prof. Robert Pego, Carnegie Mellon University, USA}

\\Abstract:\\ Mathematics is often about abstracting complicated phenomena into

simple models. This talk is about equations that model aggregation

or clustering phenomena --- think of how aerosols form soot particles

in the atmosphere, or how interplanetary dust forms comets, planets

and stars. Often in such complex systems one observes universal trend

toward self-similar growth. I'll describe an explanation for this

phenomenon in two simple models describing: (a) ``one-dimensional

bubble bath,'' and (b) the clustering of random shock waves.

In this talk, we describe new profile decompositions for bounded sequences in Banach spaces of functions defined on $\mathbb{R}^d$. In particular, for "critical spaces" of initial data for the Navier-Stokes equations, we show how these can give rise to new proofs of recent regularity theorems such as those found in the works of Escauriaza-Seregin-Sverak and Rusin-Sverak. We give an update on the state of the former and a new proof plus new results in the spirit of the latter. The new profile decompositions are constructed using wavelet theory following a method of Jaffard.

In this talk we discuss the solution of the elastodynamic

equations in a bounded domain with hereditary-type linear

viscoelastic constitutive relation. Existence, uniqueness, and

regularity of solutions to this problem is demonstrated

for those viscoelastic relaxation tensors satisfying the condition

of being completely monotone. We then consider the non-self-adjoint

and non-linear eigenvalue problem associated with the

frequency-domain form of the elastodynamic equations, and show how

the time-domain solution of the equations can be expressed in

terms of an eigenfunction expansion.

We consider a 3-dimensional elastic continuum whose material points

can experience no displacements, only rotations. This framework is a

special case of the Cosserat theory of elasticity. Rotations of

material points of the continuum are described mathematically by

attaching to each geometric point an orthonormal basis which gives a

field of orthonormal bases called the coframe. As the dynamical

variables (unknowns) of our theory we choose the coframe and a

density.

In the first part of the talk we write down the general dynamic

variational functional of our problem. In doing this we follow the

logic of classical linear elasticity with displacements replaced by

rotations and strain replaced by torsion. The corresponding

Euler-Lagrange equations turn out to be nonlinear, with the source

of this nonlinearity being purely geometric: unlike displacements,

rotations in 3D do not commute.

In the second part of the talk we present a class of explicit

solutions of our Euler-Lagrange equations. We call these solutions

plane waves. We identify two types of plane waves and calculate

their velocities.

In the third part of the talk we consider a particular case of our

theory when only one of the three rotational elastic moduli, that

corresponding to axial torsion, is nonzero. We examine this case in

detail and seek solutions which oscillate harmonically in time but

depend on the space coordinates in an arbitrary manner (this is a

far more general setting than with plane waves). We show [1] that

our second order nonlinear Euler-Lagrange equations are equivalent

to a pair of linear first order massless Dirac equations. The

crucial element of the proof is the observation that our Lagrangian

admits a factorisation.

[1] Olga Chervova and Dmitri Vassiliev, "The stationary Weyl

equation and Cosserat elasticity", preprint http://arxiv.org/abs/1001.4726

In this talk we consider the Boltzmann equation arising in gas dynamics with long-range interactions. Mathematically, it involves bilinear singular integral operators known as collision operators with non-cutoff collision kernels. As for the associated Cauchy problem, we develop a theory of weak solutions and present some of its a priori estimates related with physical quantities including the energy and moments.

The aim of this talk is to analyze energy functionals concentrated on the jump set of 2D vector fields of unit length and of vanishing divergence.

The motivation of this study comes from thin-film micromagnetics where these functionals correspond to limiting wall-energies. The main issue consists in characterizing the wall-energy density (the cost function) so that the energy functional is lower semicontinuous (l.s.c.). The key point resides in the concept of entropies due to the scalar conservation law implied by our vector fields. Our main result identifies appropriate cost functions

associated to certain sets of entropies. In particular, certain power cost functions lead to l.s.c. energy functionals.

A second issue concerns the existence of minimizers of such energy functionals that we prove via a compactness result. A natural question is whether the viscosity solution is a minimizing configuration. We show that in general it is not the case for nonconvex domains.

However, the case of convex domains is still open. It is a joint work with Benoit Merlet, Ecole Polytechnique (Paris).

One of important subjects in the study of transonic flow is to understand a global structure of flow through a convergent-divergent nozzle so called a de Laval nozzle. Depending on the pressure at the exit of the de Laval nozzle, various patterns of flow may occur. As an attempt to understand such a phenomenon, we introduce a new potential flow model called 'non-isentropic potential flow system' which allows a jump of the entropy across a shock, and use this model to rigorously prove the unique existence and the stability of transonic shocks for a fixed exit pressure. This is joint work with Mikhail Feldman.

I will discuss recent results on dispersive estimates for linear PDE's with time dependent coefficients. Then I will discuss how such

estimates can be used to study stability of nonlinear solitary waves and resonance phenomena.

Given a bounded Lipschitz domain, we consider the Dirichlet problem with boundary data in Besov spaces

for divergence form strongly elliptic systems of arbitrary order with bounded complex-valued coefficients.

The main result gives a sharp condition on the local mean oscillation of the coefficients of the differential operator

and the unit normal to the boundary (automatically satisfied if these functions belong to the space VMO)

which guarantee that the solution operator associated with this problem is an isomorphism.

According to the Ernst-Geroch reduction, in an axially symmetric stationary electrovac spacetime, the Einstein-Maxwell equations reduce to a harmonic map problem with singular boundary data. I will discuss the “regularity” of the reduced harmonic maps near the boundary and its implication on the regularity of the corresponding spacetimes.

We consider the three dimensional Navier-Stokes equations with a large initial data and

we prove the existence of a global smooth solution. The main feature of the initial data

is that it varies slowly in the vertical direction and has a norm which blows up as the

small parameter goes to zero. In the language of geometrical optics, this type of

initial data can be seen as the ``ill prepared" case. Using analytical-type estimates

and the special structure of the nonlinear term of the equation we obtain the existence

of a global smooth solution generated by this large initial data. This talk is based on a

work in collaboration with J.-Y. Chemin and I. Gallagher and on a joint work with Z.

Zhang.

We address the effect of extreme geometry on a non-convex variational problem motivated by recent investigations of magnetic domain walls trapped by sharp thin necks. We prove the existence of local minimizers representing geometrically constrained walls under suitable symmetry assumptions on the domains and provide an asymptotic characterization of the wall profile. The asymptotic behavior, which depends critically on the scaling of length and width of the neck, turns out to be qualitatively different from the higher-dimensional case and a richer variety of regimes is shown to exist.

We present an alternative viewpoint on recent studies of regularity of solutions to the Navier-Stokes equations in critical spaces. In particular, we prove that mild solutions which remain bounded in the

space $\dot H^{1/2}$ do not become singular in finite time, a result which was proved in a more general setting by L. Escauriaza, G. Seregin and V. Sverak using a different approach. We use the method of "concentration-compactness" + "rigidity theorem" which was recently developed by C. Kenig and F. Merle to treat critical dispersive equations. To the authors' knowledge, this is the first instance in which this method has been applied to a parabolic equation. This is joint work with Carlos Kenig.

The orientational order of a nematic liquid crystal in a spatially inhomogeneous flow situation is best described by a Q-tensor field because of the defects that inevitably occur. The evolution is determined by two equations. The flow is governed by a generalised Stokes equation in which the divergence of the stress tensor also depends on Q and its time derivative. The evolution of Q is governed by a convection-diffusion type equation that contains terms nonlinear in Q that stem from a Landau-de Gennes potential.

In this talk, I will show how the most general evolution equations can be derived from a dissipation principle. Based on this, I will identify a specific model with three viscosity coefficients that allows the contribution of the orientation to the viscous stress to be cast in the form of a Q-dependent body force. This leads to a convenient time-discretised strategy for solving the flow-orientation problem using two alternating steps. First, the flow field of the Stokes flow is computed for a given orientation field. Second, with the given flow field, one time step of the orientation equation is carried out. The new orientation field is then used to compute a new body force which is again used in the Stokes equation and so forth.

For some simple test applications at low Reynolds numbers, it is found that the non-homogeneous orientation of the nematic liquid crystal leads to non-linear flow effects similar to those known from Newtonian flow at high Reynolds numbers.

I will address the celebrated and long standing “No-Hair” conjecture that aims for

the classification of stationary, regular, electro-vacuum black hole space-times.

Besides reviewing some of the necessary concepts from general relativity I will

focus on the analysis of the singular harmonic map to which the source free Einstein-Maxwell

equations reduce in the stationary and axisymmetric case.

The accuracy of the Nos-Hoover thermostat to sample the Gibbs measure depends on the

dynamics being ergodic. It has been observed for a long time that this dynamics is

actually not ergodic for some simple systems, such as the harmonic oscillator.

In this talk, we rigourously prove the non-ergodicity of the Nos-Hoover thermostat, for

the one-dimensional harmonic oscillator.

We will also show that, for some multidimensional systems, the averaged dynamics for the limit

of infinite thermostat "mass" has many invariants, thus giving

theoretical support for either non-ergodicity or slow ergodization.

Our numerical experiments for a two-dimensional central force problem

and the one-dimensional pendulum problem give evidence for

non-ergodicity.

We also present numerical experiments for the Nose-Hoover chain with

two thermostats applied to the one-dimensional harmonic

oscillator. These experiments seem to support the non-ergodicity of the

dynamics if the masses of the reservoirs are large enough and are

consistent with ergodicity for more moderate masses.

Joint work with Frederic Legoll and Richard Moeckel

Scalar balance laws with monostable reaction, possibly non-convex flux, and

viscosity $\varepsilon$ are known to admit so-called entropy travelling fronts for all velocities greater than or equal to an $\varepsilon$-dependent minimal value, both when $\varepsilon$ is positive, when all fronts are smooth, and for $\varepsilon =0$, when the possibly non-convex flux results in fronts of speed close to the minimal value typically having discontinuities where jump conditions hold.

I will discuss the vanishing-viscosity limit of these fronts.

The Cosserat brothers’ ingenuous and powerful idea (presented in several papers in the Comptes Rendus at the turn of the 20th century) of solving elasticity problems by considering the homogeneous Navier equations as an eigenvalue problem is presented. The theory was taken up by Mikhlin in the 1970’s who rigorously studied it in the context of spectral analysis of pde’s, and who also presented a representation theorem for the solution of the boundary-value problems of displacement and traction in elasticity as a convergent series of the ( orthogonal and complete in the Sobolev space H1) Cosserat eigenfunctions. The feature of this representation is that the dependence of the solution on geometry, material constants and loading is provided explicitly. Recent work by the author and co-workers obtains the set of eigenfunctions for the spherical shell and compares them with the Cosserat expressions, and further explores applications and a new variational principle. In cases that the loading is orthogonal to some of the eigenfunctions, the form of the solution can be greatly simplified. Applications, in addition to elasticity theory, include thermoelasticity, poroelesticity, thermo-viscoelasticity, and incompressible Stokes flow; several examples are presented, with comparisons to known solutions, or new solutions.

The aim of my talk is to present the work of my PhD Thesis and my current research. It is concerned with local/global bifurcation of standing wave solutions to some nonlinear Schr\"odinger equations in $\mathbb{R}^N \ (N\geq1)$ and with stability properties of these solutions. The equations considered have a nonlinearity of the form $V(x)|\psi|^{p-1}\psi$, where $V:\mathbb{R}^N\to\mathbb{R}$ decays at infinity and is subject to various assumptions. In particular, $V$ could be singular at the origin.

Local/global smooth branches of solutions are obtained for the stationary equation by combining variational techniques and the implicit function theorem. The orbital stability of the corresponding standing waves is studied by means of the abstract theory of Grillakis, Shatah and Strauss.

We are interested in optimizing the compliance of an elastic structure when the applied forces are partially unknown or submitted to perturbations, the so-called "robust compliance".

For linear elasticity,the compliance is a solution to a minimizing problem of the energy. The robust compliance is then a min-max, the minimum beeing taken amongst the possible displacements and the maximum amongst the perturbations. We show that this problem is well-posed and easy to compute.

We then show that the problem is relatively easy to differentiate with respect to the domain and to compute the steepest direction of descent.

The levelset algorithm is then applied and many examples will explain the different mathematical and technical difficulties one faces when one

tries to tackle this problem.

Motivated by the tensile experiments on titanium alloys of Petrinic et al

(2006), which show the formation of cracks through the formation and

coalescence of voids in ductile fracture, we consider the problem of

formulating a variational model in nonlinear elasticity compatible both

with cavitation and with the appearance of discontinuities across

two-dimensional surfaces. As in the model for cavitation of Müller and

Spector (1995) we address this problem, which is connected to the

sequential weak continuity of the determinant of the deformation gradient

in spaces of functions having low regularity, by means of adding an

appropriate surface energy term to the elastic energy. Based upon

considerations of invertibility we are led to an expression for the

surface energy that admits a physical and a geometrical interpretation,

and that allows for the formulation of a model with better analytical

properties. We obtain, in particular, important regularity properites of

the inverses of deformations, as well as the weak continuity of the

determinants and the existence of minimizers. We show further that the

creation of surface can be modelled by carefully analyzing the jump set of

the inverses, and we point out some connections between the analysis of

cavitation and fracture, the theory of SBV functions, and the theory of

cartesian currents of Giaquinta, Modica and Soucek. (Joint work with

Carlos Mora-Corral, Basque Center for Applied Mathematics).

Modelling the behaviour of light in photonic crystal fibres requires

solving 2nd-order elliptic eigenvalue problems with discontinuous

coefficients. The eigenfunctions of these problems have limited

regularity. Therefore, the planewave expansion method would appear to

be an unusual choice of method for such problems. In this talk I

examine the convergence properties of the planewave expansion method as

well as demonstrate that smoothing the coefficients in the problem (to

get more regularity) introduces another error and this cancels any

benefit that smoothing may have.

The relaxation of a free-energy functional which describes the

order-strain interaction in nematic liquid crystal elastomers is obtained

explicitly. We work in the regime of small strains (linearized

kinematics). Adopting the uniaxial order tensor theory or Frank

model to describe the liquid crystal order, we prove that the

minima of the relaxed functional exhibit an effective biaxial

microstructure, as in de Gennes tensor model. In particular, this

implies that the response of the material is soft even if the

order of the system is assumed to be fixed. The relaxed energy

density satisfies a solenoidal quasiconvexification formula.