Past OxPDE Lunchtime Seminar

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.

I will discuss the investigatation of highly nonlinear solitary waves in heterogeneous one-dimensional granular crystals using numerical computations, asymptotics, and experiments. I will focus primarily on periodic arrangements of particles in experiments in which stiffer/heavier stainless stee are alternated with softer/lighter ones.

The governing model, which is reminiscent of the Fermi-Pasta-Ulam lattice, consists of a set of coupled ordinary differential equations that incorporate Hertzian interactions between adjacent particles. My collaborators and I find good agreement between experiments and numerics and gain additional insight by constructing an exact compaction solution to a nonlinear partial differential equation derived using long-wavelength asymptotics. This research encompasses previously-studied examples as special cases and provides key insights into the influence of heterogeneous, periodic lattice on the properties of the solitary waves.

I will briefly discuss more recent work on lattices consisting of randomized arrangements of particles, optical versus acoustic modes, and the incorporation of dissipation.

I will introduce the moving frame approach to the analysis of invariant variational problems and the evolution of differential invariants under invariant submanifold flows. Applications will include differential geometric flows, integrable systems, and image processing.

I will describe a topological approach to the motion planning problem of

robotics which leads to a new homotopy invariant of topological spaces

reflecting their "navigational complexity". Technically, this invariant is

defined as the genus (in the sense of A. Schwartz) of a specific fibration.