Past OxPDE Lunchtime Seminar

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.

The talk will survey old and recent applications of variational techniques in studying the existence, stability and bifurcations of time harmonic, localized in space solutions of the nonlinear Schroedinger equation (NLS). Such solutions are called solitons, when the equation is space invariant, and bound-states, when it is not. Due to the Hamiltonian structure of NLS, solitons/bound-states can be characterized as critical points of the energy functional restricted to sets of functions with fixed $L^2$ norm.

In general, the energy functional is not convex, nor is the set of functions with fixed $L^2$ norm closed under weak convergence. Hence the standard variational arguments fail to imply existence of global minimizers. In addition for ``critical" and ``supercritical" nonlinearities the restricted energy functional is not bounded from below. I will first review the techniques used to overcome these drawbacks.

Then I will discuss recent results in which the characterizations of bound-states as critical points (not necessarily global minima) of the restricted energy functional is used to show their orbital stability/instability with respect to the nonlinear dynamics and symmetry breaking phenomena as the $L^2$ norm of the bound-state is varied.

We consider a class of energy functionals containing a small parameter ε and a long-range interaction. Such functionals arise from models for phase separation in diblock copolymers and from stationary solutions of FitzHugh–Nagumo type systems.

On an interval of arbitrary length, we show that every global minimizer is periodic, and provide asymptotic expansions for the periods.

In 2D, periodic hexagonal structures are observed in experiments in certain di-block

copolymer melts. Using the modular function and an heuristic reduction of a mathematical model, we present a mathematical account of a hexagonal pattern selection observed in di-block copolymer melts.

We also consider the sharp interface problem arising in the singular limit,

and prove the existence and the nondegeneracy of solutions whose interface is a distorted circle in a two-dimensional bounded domain without any assumption on the symmetry of the domain.

We study the Neumann and regularity boundary value problems for a divergence form elliptic equation in the plane. We assume the gradient

of the coefficient matrix satisfies a Carleson measure condition and consider data in L^p, 1

The Navier-Stokes equation with a non-linear viscous term will be considered, p is the exponent of non-linearity.

An existence theorem is proved for the case when the convection term is not subordinate to the viscous

term, in particular for the previously open case p

There has been much recent progress in extending Griffith's criterion for

crack growth into mathematical models for quasi-static crack evolution

that are well-posed, in the sense that there exist solutions that can be

numerically approximated. However, mathematical progress in dynamic

fracture (crack growth consistent with Griffith's criterion, together with

elastodynamics) has been meager. We describe some recent results on a

phase-field model of dynamic fracture, as well as some models based on a

"sharp interface" instead of a phase-field.

Some possible strategies for showing existence for these last models will

also be described.

We discuss how various types of orientational and

translational ordering in different liquid crystal phases are

described by macroscopic tensor order parameters. In

particular, we consider a mean-field molecular-statistical

theory of the transition from the orthogonal uniaxial smectic

phase and the tilted biaxial phase composed of biaxial

molecules. The relationship between macroscopic order

parameters, molecular invariant tensors and the symmetry of

biaxial molecules is discussed in detail. Finally we use

microscopic and macroscopic symmetry arguments to consider the

mechanisms of the ferroelectric ordering in tilted smectic

phases determined by molecular chirality.

We study the asymptotics of the Stokes problem in cylinders becoming unbounded in the direction of their axis. We consider

especially the case where the forces are independent of the axis coordinate and the case where they are periodic along the axis, but the same

techniques also work in a more general framework.

We present in detail the case of constant forces (in the axial direction) since it is probably the most interesting for applications and also

because it allows to present the main ideas in the simplest way. Then we briefly present the case of periodic forces on general periodic domains. Finally, we give a result under much more general assumptions on the applied forces.