Past OxPDE Lunchtime Seminar

{\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.