Employing the usual multilevel Monte Carlo estimator, we introduce a framework for estimating the solutions of SDEs by an Euler-Maruyama scheme. By considering the expected value of such solutions, we produce simulations using approximately normal random variables, and recover the estimate from the exact normal distribution by use of a multilevel correction, leading to faster simulations without loss of accuracy. We will also highlight this concept in the framework of reduced precision and vectorised computations.

In this talk I will introduce a new mathematical model for the computational modelling of the active contraction of cardiac tissue using stress-assisted conductivity as the main mechanism for mechanoelectrical feedback. The coupling variable is the Kirchhoff stress and so the equations of hyperelasticity are written in mixed form and a suitable finite element formulation is proposed. Next I will introduce a simplified version of the coupled system, focusing on its analysis in terms of solvability and stability of continuous and discrete mixed-primal formulations, and the convergence of these methods will be illustrated through two numerical tests.

In petroleum reservoir simulation, the standard preconditioner is a two-stage process which involves solving a restricted pressure system with AMG. Initially designed for isothermal models, this approach is often used in the thermal case. However, it does not incorporate heat diffusion or the effects of temperature changes on fluid flow through viscosity and density. We seek to develop preconditioners which consider this cross-coupling between pressure and temperature. In order to study the effects of both pressure and temperature on fluid and heat flow, we first consider a model of non-isothermal single phase flow through porous media. For this model, we develop a block preconditioner with an efficient Schur complement approximation. Then, we extend this method for multiphase flow as a two-stage preconditioner.

Nets of lines are line arrangements satisfying very strict intersection conditions. We will see that nets can be defined in a very natural way in algebraic geometry, and, thanks to the strict intersection properties they satisfy, we will see that a lot can be said about classifying them over the complex numbers. Despite this, there are still basic unanswered questions about nets, which we will discuss. 
 

Accurate simulation of electromagnetic scattering by ice crystals in clouds is an important problem in atmospheric physics, with single scattering results feeding directly into the radiative transfer models used to predict long-term climate behaviour. The problem is challenging for numerical simulation methods because the ice crystals in a given cloud can be extremely varied in size and shape, sometimes exhibiting fractal-like geometrical characteristics and sometimes being many hundreds or thousands of wavelengths in diameter. In this talk I will focus on the latter "high-frequency" issue, describing a hybrid numerical-asymptotic boundary element method for the simplified problem of acoustic scattering by penetrable convex polygons, where high frequency asymptotic information is used to build a numerical approximation space capable of achieving fixed accuracy of approximation with frequency-independent computational cost. 

Harmonic map equations are an elliptic PDE system arising from the  
minimisation problem of Dirichlet energies between two manifolds. In  
this talk we present some some recent works concerning the symmetry  
and stability of harmonic maps. We construct a new family of  
''twisting'' examples of harmonic maps and discuss the existence,  
uniqueness and regularity issues. In particular, we characterise of  
singularities of minimising general axially symmetric harmonic maps,  
and construct non-minimising general axially symmetric harmonic maps  
with arbitrary 0- or 1-dimensional singular sets on the symmetry axis.  
Moreover, we prove the stability of harmonic maps from $\mathbb{B}^3$  
to $\mathbb{S}^2$ under $W^{1,p}$-perturbations of boundary data, for  
$p \geq 2$. The stability fails for $p <2$ due to Almgren--Lieb and  
Mazowiecka--Strzelecki.

(Joint work with Prof. Robert M. Hardt.)

This week is a women's only week in which we are joined by the Mirzakhani society, the society for undergraduate women in maths, for lunch.

We are very excited to have another session with invited speakers joining us for the lunch next week. Annika Heckel, Frances Kirwan and Marc Lackenby will all be joining us for a panel discussion on balancing family with academia. All are welcome to join us and to ask questions. 

We hope to see many of you at the lunch - Monday 1-2pm Quillen Room (N3.12).