From structure to dynamics in chemical reaction networks
Abstract
Chemical reaction network (CRN) theory focusses on making claims about dynamical behaviours of reaction networks which are, as far as possible, dependent on the network structure but independent of model details such as functions chosen and parameter values. The claims are generally about the existence, nature and stability of limit sets, and the possibility of bifurcations, in models of CRNs with particular structural features. The methodologies developed can often be applied to large classes of models occurring in biology and engineering, including models whose origins are not chemical in nature. Many results have a natural algorithmic formulation. Apart from the potential for application, the results are often pleasing mathematically for their power and generality.
This talk will concern some recent themes in CRN theory, particularly focussed on how the presence or absence of particular subnetworks ("motifs") influences allowed dynamical behaviours in ODE models of a CRN. A number of recent results take the form: "a CRN containing no subnetworks satisfying condition X cannot display behaviour of type Y"; but also, in the opposite direction, "if a CRN contains a subnetwork satisfying condition X, then some model of this CRN from class C admits behaviour of type Y". The proofs of such results draw on a variety of techniques from analysis, algebra, combinatorics, and convex geometry. I'll describe some of these results, outline their proofs, and sketch some current challenges in this area.
Cavitation and fracture in soft solids
Abstract
Cracks in many soft solids behave very differently to the classical picture of fracture, where cracks are long and thin, with damage localised to a crack tip. In particular, small cracks in soft solids become highly rounded — almost circular — before they start to extend. However, despite being commonplace, this is still not well understood. We use a phase-separation technique in soft, stretched solids to controllably nucleate and grow small, nascent cracks. These give insight into the soft failure process. In particular, our results suggest fracture occurs in two regimes. When a crack is large, it obeys classical linear-elastic fracture mechanics, but when it is small it grows in a new, scale-free way at a constant driving stress.
OCIAM Dinner at Christ Church, Oxford
Here's a quick note about the location and dress code for Wednesday's OCIAM event at Christ Church.
The Lecture will take place in the Michael Dummett Lecture Theatre, which is in Blue Boar Quad at 6pm. Please enter via the lodge and ask for directions.
Pre-dinner drinks at 7:15pm and dinner at 7:45pm itself will take place in the Lee Building (in the Freind room = SCR dining room. Yes, e before i.)
Given that we will be in the SCR dining room, please dress smartly (Jacket and tie for the gents, please. No jeans.)
[[{"fid":"54605","view_mode":"preview","fields":{"format":"preview"},"link_text":"Christ Church map.pdf","type":"media","field_deltas":{"1":{"format":"preview"}},"attributes":{"class":"media-element file-preview","data-delta":"1"}}]]
Poroelastic propagation and pancakes: understanding why supraglacial lakes spread but Venutian lava domes stop
Abstract
Many fluid flows in natural systems are highly complex, with an often beguilingly intricate and confusing detailed structure. Yet, as with many systems, a good deal of insight can be gained by testing the consequences of simple mathematical models that capture the essential physics. We’ll tour two such problems. In the summer melt seasons in Greenland, lakes form on the surface of the ice which have been observed to rapidly drain. The propagation of the meltwater in the subsurface couples the elastic deformation of the ice and, crucially, the flow of water within the deformable subglacial till. In this case the poroelastic deformation of the till plays a subtle, but crucial, role in routing the surface meltwater which spreads indefinitely, and has implications for how we think about large-scale motion in groundwater aquifers or geological carbon storage. In contrast, when magma erupts onto the Earth’s surface it flows before rapidly cooling and crystallising. Using analogies from the kitchen we construct, and experimentally test, a simple model of what sets the ultimate extent of magmatic intrusions on Earth and, as it turns out, on Venus. The results are delicious! In both these cases, we see how a simplified mathematical analysis provides insight into large scale phenomena.
Concentrations of solutions to compressible Navier-Stokes equations
Abstract
This work is devoted to the study of the following boundary value problem for compressible Navier-Stokes equations\begin{align*}&\begin{aligned}[b] \partial_t(\varrho \mathbf{u})+\text{div}(\varrho \mathbf u\otimes\mathbf u) &+\nabla p(\rho)\\&= \text{div} \mathbb S(\mathbf u)+\varrho\, \mathbf f\quad\text{ in }\Omega\times (0,T),\end{aligned} \\[6pt]&\partial_t\varrho+\text{div}(\varrho \mathbf u)=0\quad\text{ in }\Omega\times (0,T), \\[6pt]&\begin{aligned}[c] &\mathbf u=0\quad\text{ on }\partial \Omega\times( 0,T), \\ &\mathbf u(x,0)=\mathbf u_0(x)\quad\text{ in } \Omega,\\&\varrho(x,0)= \varrho_0(x) \quad\text{ in } \Omega, \end{aligned}\end{align*} where $\Omega$ is a bounded domain in $\mathbb R^d$, $d=2,3$, $\varrho_0>0$, $\mathbf u_0$, $\mathbf f$ are given functions, $p(\varrho)=\varrho^\gamma$, $\mathbb S(\mathbf u)=\mu(\nabla\mathbf{u}+\nabla\mathbf{u}^\top)+\lambda \text{div } \mathbf{u}$, $\mu, \lambda$ are positive constants. We consider the endpoint cases $\gamma=3/2$, $d=3$ and $\gamma=1$, $d=2$, when the energy estimate does not guarantee the integrability of the kinetic energy density with an exponent greater than 1, which leads to the so-called concentration problem. In order to cope with this difficulty we develop new approach to the problem. Our method is based on the estimates of the Newton potential of $p(\varrho)$. We prove that the kinetic energy density in 3-dimensional problem with $\gamma=3/2$ is bounded in $L\log L^\alpha$ Orlitz space and obtain new estimates for the pressure function. In the case $d=2$ and $\gamma=1$ we prove the existence of the weak solution to the problem. We also discuss the structure of concentrations for rotationally-symmetric and stationary solutions.
Diffeomorphic Approximation of W^{1,1} Planar Sobolev Homeomorphisms
Abstract
Let $\Omega\subseteq\mathbb{R}^2$ be a domain and let $f\in W^{1,1}(\Omega,\mathbb{R}^2)$ be a homeomorphism (between $\Omega$ and $f(\Omega)$). Then there exists a sequence of smooth diffeomorphisms $f_k$ converging to $f$ in $W^{1,1}(\Omega,\mathbb{R}^2)$ and uniformly. This is a joint result with A. Pratelli.
Hypoelliptic Laplacian, Brownian motion and the trace formula
Abstract
The hypoelliptic Laplacian is a family of operators indexed by $b \in \mathbf{R}^*_+$, acting on the total space of the tangent bundle of a Riemannian manifold, that interpolates between the ordinary Laplacian as $b \to 0$ and the generator of the geodesic flow as $b \to +\infty$. These operators are not elliptic, they are not self-adjoint, they are hypoelliptic. One can think of the total space of the tangent bundle as the phase space of classical mechanics; so that the hypoelliptic Laplacian produces an interpolation between the geodesic flow and its quantisation. There is a dynamical counterpart, which is a natural interpolation between classical Brownian motion and the geodesic flow.
The hypoelliptic deformation preserves subtle invariants of the Laplacian. In the case of locally symmetric spaces (which are defined via Lie groups), the deformation is essentially isospectral, and leads to geometric formulas for orbital integrals, a key ingredient in Selberg's trace formula.
In a first part of the talk, I will describe the geometric construction of the hypoelliptic Laplacian in the context of de Rham theory. In a second part, I will explain applications to the trace formula.
Ginzburg–Landau functionals with a general compact vacuum manifold on planar domains
Abstract
Ginzburg–Landau type functionals provide a relaxation scheme to construct harmonic maps in the presence of topological obstructions. They arise in superconductivity models, in liquid crystal models (Landau–de Gennes functional) and in the generation of cross-fields in meshing. For a general compact manifold target space we describe the asymptotic number, type and location of singularities that arise in minimizers. We cover in particular the case where the fundamental group of the vacuum manifold in nonabelian and hence the singularities cannot be characterized univocally as elements of the fundamental group. The results unify the existing theory and cover new situations and problems.
This is a joint work with Antonin Monteil and Rémy Rodiac (UCLouvain, Louvain- la-Neuve, Belgium)
On boundary value problem for steady Navier-Stokes system in 2D exterior domains
Abstract
We study solutions to stationary Navier-Stokes system in two dimensional exterior domains, namely, existence of these solutions and their asymptotical behavior. The talk is based on the recent joint papers with K. Pileckas and R. Russo where the uniform boundedness and uniform convergence at infinity for arbitrary solution with finite Dirichlet integral were established. Here no restrictions on smallness of fluxes are assumed, etc. In the proofs we develop the ideas of the classical papers of Gilbarg & H.F. Weinberger (Ann. Scuola Norm.Pisa 1978) and Amick (Acta Math. 1988).