International Conference on Free Boundary Problems: Theory and Applications
Monday 23rd June 2014 to Friday 27th June 2014
09:00 to 09:55  Registration  INI 1  
09:55 to 10:00  Welcome from John Toland (INI Director)  INI 1  
10:00 to 10:45 
R Goldstein (University of Cambridge) Plenary Lecture 1: How a Volvox embryo turns itself inside out Chair: GQ Chen
During the growth of daughter colonies of the multicellular alga Volvox the spherical embryos must turn themselves inside out to complete their development. This process of 'inversion' has many features in common with gastrulation, the process by which an initially convex spherical shell of animal cells develops an invagination, leading to the formation of a gastric system. In both cases it is understood that cell shape changes play a major role in guiding the process, but quantification of the dynamics, and formulation of a mathematical description of the process, have been lacking. In this talk I will describe advances my group has made recently on both fronts. Using the technique of SPIM (selective plane illumination microscopy) we have obtained the first realtime threedimensional timelapse movies of inversion in Volvox, using several species displaying distinct morphological events. The beginnings of an elastic theory of these processes will also be descri bed.

INI 1  
11:00 to 11:30  Morning Coffee  
11:30 to 12:15 
B Wagner (Technische Universität Berlin) Plenary Lecture 2: Unstaedy nonuniform base states and their stability Chair: GQ Chen
In this talk we consider several pattern forming systems, ranging from phase separation of polymer blends, selfassembly of crystalline films to dewetting of polymer films. These systems all have unstaedy nonuniform base states. We develop asymptotic techniques to analyse their associated linear stability problems and derive expressions for predicting the dominant wavelength of the pattern.

INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:15 
Plenary Lecture 3: Free Boundaries and Fluid Mixing at the Micro Level Chair: JL Vázquez
Turbulent mixing often occurs with immiscible fluids or with miscible fluids over rapid time scales, so that the flow is locally inhomogeneous at a micro level for periods of interest. We start with a review of problems in which such flows arise.
The flow regions in which the mixing occurs can generally be identified reliably; examples will be given. The challenge for current research is to describe the microscopic and inhomogeneous mixture in a statistical sense. Full resolution of the flows is generally out of the question and will remain so for decades. Thus we are interested in statistical properties of the flow that are stable and appear to converge under mesh refinement, with sufficient detail in the statistical description (for example a pdf or cdf (cumulative distribution function) to support reaction processes in flows of engineering interest. A first step, generally insufficient, is to compute means and variances of fluctuating processes. This goal is still in the future. Partial results leading in this direction will be presented. We formulate a notion of stochastic convergence and present numerical algorithms which appear to be convergent in this metric. We introduce theoretical ideas related to convergence based on the renormalization group, including the important notion that the solution, at the LES level of resolution of necessity considered here, is not unique. In other words, the usual standard of convergence under mesh refinement is not sufficient to guarantee a simulation in agreement with experimental data. We discuss methods to mitigate this serious obstacle to scientific progress. Basically, experiments are essential to select the correct nonunique solution and the algorithm and its adjustable parameters to reach this goal. While flows of interest are commonly at high Reynolds numbers outside the regime of relevant experiments, the expansion parameter is 1/Reynolds number. In terms of this parameter, the perturbation from experiment to applications is small, within normally accepted ranges for perturbative extensions of validation regimes. "cambridge.14.abs" 37L, 2201C 
INI 1  
14:20 to 14:50  Afternoon Tea  
14:50 to 15:20 
Some Free Boundary Problems for Flows at high Re Session: Focus Session Fluids (parallel session) Chair: A Lacey
This talk will describe some models and open problems for (i) disturbance waves in annular 2phase flow and (ii) bubble collapse.

INI 1  
14:50 to 15:20 
Dynamics of an Euler beam with unilateral constraints Session: Focus Session Recent Developments 1 (parallel session) Chair: D Bucur
We study the vibrations of an elastic beam between rigid obstacles. The non penetrability condition leads to a description of the dynamics as a hyperbolic fourth order variational inequality. For this free boundary problem we construct a sequence of approximate solutions by combining some classical space discretizations with timestepping schemes especially suited to unilateral contact problems for discrete mechanical systems. We prove the stability and the convergence of these numerical methods and we obtain an existence result for our original problem under very general assumptions on the geometry of the obstacles.

INI 2  
15:20 to 15:50 
Nondegeneracy in the Obstacle Problem with a Degenerate Force Term Session: Focus Session Recent Developments 1 (parallel session) Chair: D Bucur
In this talk I present the proof of the optimal nondegeneracy of the solution $u$ of the obstacle problem $\triangle u=f\chi_{\{u>0\}}$ in a bounded domain $D\subset\mathbb{R}^{n}$, where we only require $f$ to have a nondegeneracy of the type $f(x)\geq\lambda\vert (x_1,\cdots,x_p)\vert^{\alpha}$ for some $\lambda>0$, $1\leq p\leq n$ (an integer) and $\alpha>0$. We prove optimal uniform $(2+\alpha)$th order and nonuniform quadratic nondegeneracy, more precisely we prove that there exists $C>0$ (depending only on $n$, $p$ and $\alpha$) such that for $x$ a free boundary point and $r>0$ small enough we have $\sup_{\partial B_r(x)}u\geq C\lambda (r^{2+\alpha}+\vert(x_1,\cdots,x_p)\vert^{\alpha}r^{2})$. I also present the proof of the optimal growth with the assumption $\vert f(x)\vert\leq\Lambda\vert (x_1,\cdots,x_p)\vert^{\alpha}$ for some $\Lambda\geq 0$ and the porosity of the free boundary.
Preprint: http://www.newton.ac.uk/preprints/NI14045.pdf 
INI 2  
15:20 to 15:50 
M Fontelos (ICMAT) Shapes of charged drops in an electric field and Rayleigh jets Session: Focus Session Fluids (parallel session) Chair: A Lacey
When a viscous drop is electrically charged or subject to an external electric field, it may undergo deformations and, occasionally produce thin jets (called Rayleigh jets) that are of practical interest. We study the drop's shapes and related bifurcation problems, as well as the possible mechanisms for jet formation.

INI 1  
16:00 to 16:30 
J Oliver (University of Oxford) On contactline dynamics with mass transfer Session: Focus Session Fluids (parallel session) Chair: A Lacey
We investigate the effect of mass transfer on the evolution of a thin (twodimensional) partially wetting drop. While the effects of viscous dissipation, capillarity, slip and uniform mass transfer are taken into account, other effects, such as gravity, surface tension gradients, vapour transport and heat transport, are neglected in favour of mathematical tractability. Our matchedasymptotic analysis reveals that the leadingorder outer formulation and contactline law that is selected in the smallslip limit depends delicately on both the sign and size of the mass transfer flux, leading in particular to novel generalisations of Tanner's law. We analyse the resulting evolution of the drop and report good agreement with numerical simulations.
CoAuthors:
Matthew Saxton (University of Oxford), Jonathan Whiteley (University of Oxford), Dominic Vella (University of Oxford); John King (University of Nottingham).

INI 1  
16:00 to 16:30 
Spreading profile and nonlinear Stefan problems Session: Focus Session Recent Developments 1 (parallel session) Chair: D Bucur
I will report some recent progresses (in joint works with Z. Guo, B. Lou, H. Matsuzawa, H. Matano, K. Wang, M. Zhou etc.) on the study of a general nonlinear Stefan problem, used as a model for the understanding of a variety of spreading phenomena, where the unknown function u(t,x) represents the density of concentration of a certain (chemical or biological) species at time t and space location x, with the free boundary standing for the spreading front. Such spreading problems are usually modeled by the corresponding Cauchy problem, which has attracted extensive research starting from the wellknown 1937 paper of KolmogorovPetrovskiPiskunov. We will discuss the similarity and differences of the longtime behavior of these two types of mathematical models by closely examining their spreading profiles.

INI 2  
16:30 to 17:00 
Surfactants in twophase flow Session: Focus Session Fluids (parallel session) Chair: A Lacey
Coauthors: Harald Garcke (Universitat Regensburg), Bjorn Stinner (University of Warwick)
Emulsification is an important industrial process that involves mixing two or more fluids that normally are unmixable. More precisely, in the process of emulsification, it is desirable to have stable dispersions of one fluid in another. Common examples of emulsions are milk, fire extinguishers and hand cream. The mixture is thermodynamically unstable and will progressively revert back to their unmixed states over time. Surface active agents (or surfactants) are often added to increase the stability of the mixture and hence there is great interest, especially in industrial applications, to understand the influence of surfactants on the dynamics of multifluid systems. In this talk, I will outline the derivation of two new models on soluble surfactants in two phase flow. The first is a sharp interface model that describes the interfaces with moving hypersurfaces, and the second is a phase field model. Both models are thermodynamically consistent and generalise existing models in the literature. I will then discuss the relationship between two models, which is supported with some numerical simulations. 
INI 1  
16:30 to 17:00 
A model for the quasistatic crack growth in hydraulic fracture Session: Focus Session Recent Developments 1 (parallel session) Chair: D Bucur
We present a variational model for the quasistatic crack growth in hydraulic fracture in the framework of the energy formulation of rateindependent processes. The cracks are assumed to lie on a prescribed plane and to satisfy a very weak regularity assumption.

INI 2  
17:00 to 18:00  Welcome Wine Reception and Poster Session 
09:00 to 09:45 
Plenary Lecture 4 Chair: J King 
INI 1  
09:50 to 10:20  Morning Coffee  
10:20 to 10:50 
Multiscale analyses of tissue growth and front propagation Session: Focus Session FBPs in Biology (parallel session) Chair: R Kornhuber
The derivation of continuum models which represent underlying discrete or microscale phenomena is emerging as an important part of mathematical biology: integration between subcellular, cellular and tissuelevel behaviour is crucial to understanding tissue growth and mechanics. I will consider the application of a multiscale method to two problems on this theme.
Firstly a new macroscale description of nutrientlimited tissue growth, which is formulated as a microscale movingboundary problem within a porous medium, is introduced. A multiscale homogenisation method is employed to enable explicit accommodation of the influence of the underlying microscale tissue structure, and its evolution, on the macroscale dynamics. A challenging consideration in continuum models of tissue is the accommodation of (spatiallydiscrete) cellsignalling events, a feature of which being the progression of moving fronts of cellsignalling activity across a lattice. New (continuum) analyses of monotone waves in a discrete diffusion equation are presented, and extended to modulated fronts exhibited in cell signalling models. 
INI 1  
10:20 to 10:50 
Weak shock diffraction Session: Focus Session Mixed FBPs (parallel session) Chair: M Feldman
Coauthor: Allen Tesdall (CUNY)
We study the diffraction of a weak, selfsimilar shock in two space dimensions near a point where its shock strength approaches zero and the shock turns continuously into an expansion wavefront. For example, this happened when a weak shock hits a semiinfinite screen. The local asymptotic solution satisfies the unsteady transonic small disturbance equation. We also consider a related halfspace problem where a shock whose strength approaches zero reflects off a ``soft'' boundary. Numerical solutions show a complex reflection pattern similar to one that occurs in the Guderley Mach reflection of weak shocks. 
INI 2  
10:50 to 11:20 
Dielectric Boundary in Biomolecular Solvation Session: Focus Session FBPs in Biology (parallel session) Chair: R Kornhuber
A dielectric boundary in a biomolecular system is a solutesolvent (e.g., proteinwater) interface that defines the dielectric coefficient to be one value in the solute region and another in solvent. The inhomogeneous dielectric medium gives rise to an effective dielectric boundary force that is crucial to the biomolecular conformation and dynamics. This lecture presents a precise definition and explicit formula of the dielectric boundary force based on the classical PoissonBoltzmann theory of electrostatic interactions. These results are used to study the motion of a cylindrical dielectric boundary driven by the competition between the surface tension, dielectric boundary force, and solvent viscous force. Implications of the mathematical findings to biomolecular conformational stabilities are discussed.

INI 1  
10:50 to 11:20 
Shock diffraction by convex cornered wedges Session: Focus Session Mixed FBPs (parallel session) Chair: M Feldman
Coauthor: Prof. GuiQiang G. Chen (University of Oxford)
In this talk, I would like to present one of our current research projects, that is on the mathematical analysis of shock diffraction by convex cornered wedges. The existence of the regular configuration is established up to the critical wedge angle, which should be the criterion of the transition between the regular configuration and the Mach configuration. This work is joint with Prof. G.Q. Chen. 
INI 2  
11:30 to 12:00 
The Signorini problem, fractional Laplacians and the language of semigroups Session: Focus Session Mixed FBPs (parallel session) Chair: M Feldman
The Signorini problem can be equivalently formulated as a thin obstacle problem for an elastic membrane. The resulting free boundary problem turns out to be equivalent to the obstacle problem for the fractional Laplacian on the whole space. We will show how to understand this problem under the light of the language of semigroups that I developed in my PhD thesis (2010). In particular, we are able to consider different kinds of Signorini problems that are equivalent to obstacle problems for fractional powers of operators different than the Laplacian on the whole space. Boundary conditions of different kinds (Dirichlet, Neumann, periodic) and radial solutions can also be treated with this unified language. Another advantage is that this language avoids the use of the Fourier transform. The basic regularity results (Harnack inequalities, Schauder estimates) for these fractional nonlocal operators can be studied by means of the generalization of the CaffarelliSilvestre extensio n problem that I proved in my PhD thesis. It turns out that the solution for the extension problem can be written in terms of the heat semigroup.

INI 2  
11:30 to 12:00 
Numerical investigations of long range segregation systems Session: Focus Session FBPs in Biology (parallel session) Chair: R Kornhuber
We investigate numerical approximation for a class of elliptic and parabolic competitiondiffusion systems of long range segregation models for two and more competing species.
Also we study the long term behavior for parabolic system. We prove that as the competition rate goes to infinity the solution converges, along with suitable sequences, to a spatially long range segregated state satisfying some differential inequalities.
Moreover, we provide numerical simulations for parabolic and elliptic cases.

INI 1  
12:00 to 12:30 
B Fang (Shanghai Jiao Tong University) On Stability of Steady Transonic Shocks in Supersonic Flow around a Wedge Session: Focus Session Mixed FBPs (parallel session) Chair: M Feldman
In this talk we are concerned with the stability of steady transonic shocks in supersonic flow around a wedge. When a uniform supersonic flow comes against a straight wedge with small vertex angle, a shockfront attached to the wedge appears. It is indicated in the book "Supersonic Flow and Shock Waves" by Courant and Friedrichs that there are two admissible shock solutions satisfying both RankineHugoniot conditions and the entropy condition. The weaker shock solution may be transonic, while the stronger one must be transonic. In this talk, we shall present stability results for 2D and MD flows.
This talk is based on the joint works with Prof. G.Q. Chen, and Prof. S.X. Chen. 
INI 2  
12:00 to 12:30 
C Venkataraman (University of Sussex) Moving boundary problems in cell biology Session: Focus Session FBPs in Biology (parallel session) Chair: R Kornhuber
We discuss the mathematical formulation and numerical solution of free and moving boundary problems that arise in the study of cell motility. We consider the dynamics of the cell membrane and species that reside on the cell membrane. We also present novel finite element methods for the simulation of the model equations and report on numerical results.

INI 1  
12:30 to 13:30  Sandwich lunch at INI  
13:30 to 13:35 
J King & C Elliot Welcome and Introduction Chair: C Elliot 
INI 1  
13:35 to 13:45 
Free Boundary Problems in Biology and Medicine Chair: C Elliot 
INI 1  
13:45 to 14:10 
How the Geometry of the Cell Boundary Couples Cellular Blebbing and Actin Based Protrusions Chair: C Elliot 
INI 1  
14:10 to 14:35 
Mathematical Modelling of Tissue Growth in a Perfusion Bioreactor Chair: C Elliot 
INI 1  
14:35 to 15:00 
J Lowengrub (University of California, Irvine) Feedback, Lineages and Cancer Therapy Chair: C Elliot 
INI 1  
15:00 to 15:20  Tea and Coffee Break  
15:20 to 15:45 
Free Boundary Problems in Models for Bacterial Biofilms Chair: C Elliot 
INI 1  
15:45 to 16:10 
Modelling the Spatial Temporal Dynamics of Small Monomeric G Proteins Chair: C Elliot 
INI 1  
16:10 to 17:00 
Open Discussion Chair: C Elliot 
INI 1  
17:00 to 18:00  Drinks Reception and Networking 
09:00 to 09:45 
Plenary Lecture 5: Obstacle problems in geometric evolutions Chair: H Shahgholian
I will consider geometric evolutions of fronts in the presence of obstacles. I will discuss in particular the mean curvature flow and the Willmore flow, presenting some new results and some open questions.

INI 1  
09:50 to 10:20  Morning Coffee  
10:20 to 11:05 
P Markowich (KAUST) Plenary Lecture 6: On Boltzmanntype and Free Boundary Models for Price Formation Chair: H Shahgholian
We present a theory for a price formation free boundary model formulated by Lasry and Lions. Also we introduce a mesoscopic Boltzmanntype model for price formation in economic markets of single goods, which gives the LasryLions model in the limit of large trading frequency.

INI 1  
11:15 to 11:45 
Symmetries of Boundary Value Problems: Definitions, Algorithms and Applications to Physically Motivated Problems Session: Focus Session Recent Developments 2A (parallel session) Chair: N Wolanski
One may note that the symmetrybased methods were not widely used for solving boundaryvalue problems (BVPs). To the best of our knowledge, the first rigorous definition of Lie's invariance for BVPs was formulated by George Bluman in early 1970s and applied to some classical BVPs.
However, Bluman's definition cannot be directly applied to BVPs of more general form, for example, to those involving boundary conditions on the moving surfaces, which are described by unknown functions. In our recent papers, a new definition of Lie's invariance of BVP with a wide range of boundary conditions (including those at infinity and moving surfaces) was formulated. Moreover, an algorithm of the group classification for the given class of BVPs was worked out. The definition and algorithm were applied to some classes of nonlinear twodimensional and multidimensional BVPs of Stefan type with the aim to show their efficiency. In particular, the group classification problem for these classes of BVPs was solved, reductions to BVPs of lower dimensionality were constructed and examples of exact solutions with physical meaning were found.
Very recently, the definition and algorithms were extended on the case of conditional invariance for BVPs and applied to some nonlinear BVPs.
This research was supported by a Marie Curie International Incoming Fellowship within the 7th European Community Framework Programme.

INI 1  
11:15 to 11:45 
Generalized Neumann solutions for the twophase fractional Lam\'eClapeyronStefan problems Session: Focus Session Recent Developments 2B (parallel session) Chair: K Yeressian Negarchi
A generalized Neumann solution for the twophase fractional Lam\'eClapeyronStefan problem for a semiinfinite material with constant boundary and initial conditions is obtained. In this problem, the two governing diffusion equations and a governing condition for the free boundary include a fractional time derivative in the Caputo sense of order 0 Moreover, we also obtain a new generalized Neumann solutions for the twophase fractional Lam\'eClapeyronStefan problems when a heat flux or a convective boundary condition is imposed on the fixed face x=0. These solutions are obtained when an inequality for the coefficient which characterizes the heat flux or the convective boundary condition is satisfied. In these cases, when a? 1 we also recover the explicit solution given in Tarzia, Quart. Appl. Math., 39 (1981), 491499 and in MATSerie A, 8 (2004), 2127. The first part of this presentation is a joint work with S.D. Roscani (Univ. Rosario, Argentina). 
INI 2  
11:45 to 12:15 
J Gwinner (Universität der Bundeswehr München) Free Boundary Value Problems in Contact Mechanics Session: Focus Session Recent Developments 2A (parallel session) Chair: N Wolanski
In this talk we report on recent progress in the solution of various nonsmooth free boundary value problems that arise in contact mechanics, including unilateral contact with Tresca friction and nonmonotone adhesion/delamination problems. We analyse these problems and employ for their numerical treatment optimization and regularization techniques combined with fi nite element discretisation.
This presentation is based on the recent papers [1, 2] and the recent PhD Thesis [3] of N. Ovcharova.
[1] J.Gwinner, hpFEM convergence for unilateral contact problems with Tresca friction in plane linear elastostatics, J. Comput. Appl. Math. 254 (2013) 175 184 [2] N. Ovcharova and J.Gwinner, A study of regularization techniques of nondi erentiable optimization in view of application to hemivariational inequalities, J. Optim Theory Appl., DOI 10.1007/s109570140521y (2014) [3] N. Ovcharova, Regularization Methods and Finite Element Approximation of Some Hemivariational Inequalities with Applications to Nonmonotone Contact Problems, Ph.D. Thesis, Universitaet der Bundeswehr Muenchen, 2012.

INI 1  
11:45 to 12:15 
Convergence of optimal control problems governed by second kind parabolic variational inequalities Session: Focus Session Recent Developments 2B (parallel session) Chair: K Yeressian Negarchi
We consider a family of optimal control problems where the control variable is given by a free boundary condition of Neumann type. This family is governed by parabolic variational inequalities of the second kind. We prove the strong convergence of the optimal controls and state systems associated to this family to a similar optimal control problem.

INI 2  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:15 
K Deckelnick (OttovonGuerickeUniversität Magdeburg) Plenary Lecture 7: Double obstacle phase field approach for an elliptic inverse problem with discontinuous coefficients Chair: J Lowengrub
We consider the inverse problem of recovering interfaces where the diffusion coefficient in an elliptic PDE has jump discontinuities. We employ a least squares approach together with a perimeter regularization. A suitable relaxation of the perimeter leads to a sequence of CahnHilliard type functionals for which we obtain a $\Gamma$convergence result. Using a finite element discretization of the elliptic PDE and a suitable
adjoint problem we derive an iterative method in order to approximate discrete critical points. We prove convergence of the iteration and present results of numerical tests.
This is joint work with C.M. Elliott (Warwick) and V. Styles (Sussex).

INI 1  
14:20 to 14:50  Afternoon Tea  
14:50 to 15:20 
T Ranner (University of Leeds) Computational methods for an optimal partition problem on surfaces Session: Focus Session Numerics for FBPS (parallel session) Chair: C Elliott
In this talk I will explore computational techniques for solving a free boundary problem posed on a surface. The problem is to divide a surface into regions such that the sum of the first eigenvalue of the Dirichlet LaplaceBeltrami operator over each region is minimized. Different relaxations of this problem will be explored. Each takes the form of large system of partial differential equations which will be solved using algorithms designed for high performance computing techniques.

INI 1  
14:50 to 15:20 
Optimal Regularity and the Free Boundary in the Parabolic Signorini Problem Session: Focus Session Analysis 1 (parallel session) Chair: H Shahgholian
Coauthors: Nicola Garofalo (University of Padova), Arshak Petrosyan (Purdue University), Tung To (Purdue University)
In this talk we will give an overview of our comprehensive treatment of the parabolic Signorini problem, based on a generalization of Almgren's monotonicity of the frequency. In particular, we will discuss the optimal regularity of solutions, the classi fication of free boundary points, the regularity of the regular set, and the structure of the singular set. 
INI 2  
15:20 to 15:50 
Multimaterial structural topology optimization using phase field methods Session: Focus Session Numerics for FBPS (parallel session) Chair: C Elliott
Coauthors: Luise Blank (University of Regensburg), Hassan FarshbafShaker (WIAS Berlin), Harald Garcke (University of Regensburg)
A phase field approach for structural topology optimization which allows for topology changes and multiple materials is analyzed. First order optimality conditions are rigorously derived and several numerical results for mean compliance problems are presented. 
INI 1  
15:20 to 15:50 
D Bucur (Université de Savoie) Isoperimetric inequalities and free discontinuity problems Session: Focus Session Analysis 1 (parallel session) Chair: H Shahgholian
The FaberKrahn inequality for the first eigenvalue of the Laplace operator with Robin boundary conditions was recently proved by D. Daners in the context of Lipschitz sets. This talk introduces a new approach to deal with this type of isoperimetric inequalities and is based on recent advances on the regularity theory for general free discontinuity problems, in particular on a new monotonicity formula jointly obtained with S. Luckhaus. As a consequence, we not only extend isoperimetric inequalities to arbitrary sets, but also obtain a range of new ones (joint work with A. Giacomini).

INI 2  
16:00 to 16:30 
A stable discretization for the simulation of twophase flow Session: Focus Session Numerics for FBPS (parallel session) Chair: C Elliott
Coauthors: Harald Garcke (University of Regensburg), Michael Hinze (University of Hamburg)
In [Abels, Garcke, Gruen, M3AS, 22, 2012] a model for the simulation of variable density twophase flow using the diffuse interface approach is presented. Due to an additional transport term in the NavierStokes system this model enjoys consistency with thermodynamics. To preserve this important property in the numerical simulation we propose a time discretization that uses a carefully designed weak formulation and delivers a time discrete variant of the continuous energy equality. The spatial discretization is enriched by a suitable post processing step to ensure the energy inequality to hold in the fully discrete setting. An aposteriori adaptive procedure is proposed that deals with the errors both in the NavierStokes and the CahnHilliard system. 
INI 1  
16:00 to 16:30 
KA Lee (Seoul National University) Curvature Flows with Free Boundaries Session: Focus Session Analysis 1 (parallel session) Chair: H Shahgholian
In this talk, we are going to discuss alpha Gauss Curvature Flows with free boundaries. For a=1n, a flat spot on convex manifold will persists for some time, which is similar as the waiting time in Porous Medium Equations. Naturally, the interface between strictly convex region and flat spot will form a free boundary that will moves along the time. We will show the relation between optimal regularity and bounded nondegenerate speed of free boundaries. The optimal regularity and the smoothness of the free boundary will be discussed. We will also consider the dynamics and regularity of free boundary of the manifold with a concave region.

INI 2  
16:30 to 17:00 
L Banas ([University of Bielefeld]) Numerical approximation of phasefield models for multiphase flow Session: Focus Session Numerics for FBPS (parallel session) Chair: C Elliott
Flow of mixtures of incompressible fluids with complex fluid interactions (surface tensions, contact angles) can be described by a system of (multicomponent) CahnHilliardNavierStokes equations. We propose finite element based numerical approximations of some nonsmooth phasefield models for mixtures of incompressible fluids with variable densities and viscosities. We discuss theoretical and practical issues related to the proposed numerical approximations and present computational experiments to demonstrate main features of the considered models.

INI 1  
16:30 to 17:00 
A Karakhanyan (University of Edinburgh) The Proximity of Free Boundary to Wedges and Planes Session: Focus Session Analysis 1 (parallel session) Chair: H Shahgholian
The talk is about the geometry of free boundary striking the centre of the half ball B_1^+ where a wedge like boundary data is prescribed near the origin. It will be seen that the effect of the wedgelike data dooms the asymptotic shape of the free boundary and the homogeneous global solutions turn out to have wedgelike free boundary too. A delicate analysis involving monotonicity formula techniques is employed to furnish the main results, which are from a joint paper with H. Shahgholian (to appear in the Transactions of the AMS).

INI 2  
17:20 to 18:00  FRBW04 organiser's meeting  
19:30 to 22:00  Conference Dinner at Christ's College 
09:00 to 09:45 
J Ball (University of Oxford) Plenary Lecture 8: Interfaces arising from solid phase transformations Chair: J Toland
Solid phase transformations give rise to a variety of unusual kinds of interfaces between different phases, some only observed in recent experiments. The lecture will discuss ways of describing and predicting these, and related questions of nonlinear analysis.

INI 1  
09:50 to 10:20  Morning Coffee  
10:20 to 10:50 
Free BoundaryValue Problems for Liquid Crystals Session: Focus Session Material Science (parallel session) Chair: B Wagner
We discuss modelling strategies for free boundaryvalue problems arising in liquid crystal applications. We have recently modelled stable liquid crystal patterns in shallow cuboidshaped nematic wells and have computed phase diagrams as a function of geometrical aspect ratio and anchoring strength on well boundaries. Recent work by Telo da Gama group (CFTC, Lisbon) shows that our computed patterns are also relevant to wedgeshaped nematicfilled wells at the isotropicnematic transition temperature. At the isotropicnematic transition temperature, the isotropic and nematic phases are separated by a free boundary, depending on the geometry and the surface effects at play. We discuss modelling approaches to this freeboundary problem and the role of such free boundaries in applications. This is joint work with Chong Luo, Samo Kralj and Samo Kralj.

INI 1  
10:20 to 10:50 
Free boundary problems in glacial hydrology Session: Focus Session Recent Developments 3 (parallel session) Chair: D Danielli
Many aspects of a glacier’s dynamics depend on the evolution of an interface between slowly creeping ice and fast flowing water. As ice can be modelled as a very viscous fluid, variations on free boundary Stokes flow arise in the modelling of such interfaces.
In this talk I will focus on the shape of a channel that is cut through or under the ice by viscous dissipation of meltwater. Networks of these channels form a vital part of a glacier’s hydrological system, which carries meltwater from the surface to the glacier margin, where mixing with dense ocean water has a strong effect on melting at the ice face. The evolving cross section of the channel is related to the problem of a contracting or expanding bubble in two dimensional Stokes flow, which allows us to derive analytic results. I will also discuss the impact of a meltwater source, such as a channel, on the spatial distribution of melting at the ice face.

INI 2  
10:50 to 11:20 
On the production of highporosity metallic solids Session: Focus Session Material Science (parallel session) Chair: B Wagner
Coauthor: Stephen Davis (Northwestern University)
Highporosity metallic solids are useful as lightweight materials in many engineering disciplines. However, batch processing techniques for producing such materials by solidification of a molten metallic foam are problematic, and it is not currently possible to control the porosity distribution of the final product {\it a priori}. The molten metallic foam is inherently unstable; the thin liquid bridges between bubbles drain rapidly and rupture due to intermolecular forces, leading to bubble coalescence and largescale topological rearrangement. In this talk we will consider the competition between coalesecence and freezing in these dynamically evolving foams, examining in particular the coupling between the microscale hydrodynamics in the molten liquid films and the progression of a solidification front. The foam is modelled as a coupled system of free boundary problems involving both liquid/gas and liquid/solid interfaces, which we solve using a combination of asymptotic and numerical techniques. An understanding of these microscale processes motivates other protocols for solid foam production where the porosity distribution of the final product can be more effectively controlled. 
INI 1  
10:50 to 11:20 
A classical Perron method for existence of smooth solutions to boundary value and obstacle problems for boundarydegenerate elliptic operators via holomorphic m Session: Focus Session Recent Developments 3 (parallel session) Chair: D Danielli
We prove existence of solutions to boundary value problems and obstacle problems for degenerateelliptic, linear, secondorder partial differential operators with partial Dirichlet boundary conditions using a new version of the Perron method. The elliptic operators considered have a degeneracy along a portion of the domain boundary which is similar to the degeneracy of a model linear operator identified by Daskalopoulos and Hamilton (1998) in their study of the porous medium equation or the degeneracy of the Heston operator (Heston, 1993) in mathematical fi nance. Our Perron method relies on weak and strong maximum principles for degenerateelliptic operators, concepts of continuous subsolutions and supersolutions for boundary value and obstacle problems for degenerateelliptic operators, and maximum and comparison principle estimates previously developed by us.

INI 2  
11:30 to 12:00 
A phase field approach for shape optimization in fluids Session: Focus Session Material Science (parallel session) Chair: B Wagner
Coauthor: Harald Garcke (University of Regensburg)
We consider the problem of shape optimization with a general objective functional using the Stokes equations as state constraints. By combining a porous medium and a phase field approach we obtain a wellposed problem in a diffuse interface setting. We discuss wellposedness and optimality conditions and relate the phase field problem to a perimeter penalized sharp interface problem. 
INI 1  
11:30 to 12:00 
Existence and Qualitative Properties of Grounds States to the ChoquardType Equations Session: Focus Session Recent Developments 3 (parallel session) Chair: D Danielli
Coauthor: Jean Van Schaftingen (LouvainlaNeuve, Belgium)
The Choquard equation, also known as the Hartree equation or nonlinear SchrodingerNewton equation is a stationary nonlinear Schrodinger type equation where the nonlinearity is coupled with a nonlocal convolution term given by an attractive gravitational potential. We present sharp Liouvilletype theorems on nonexistence of positive supersolutions of such equations in exterior domains and consider existence, positivity, symmetry and optimal decay properties of ground state solutions under various assumptions on the decay of the external potential and the shape of the nonlinearity. We also discuss the existence of semiclassical solutions to the equation. 
INI 2  
12:00 to 12:30 
M Glicksman (University of Florida) Deterministic Pattern Formation in DiffusionLimited Systems Session: Focus Session Material Science (parallel session) Chair: B Wagner
An interface evolving under local equilibrium develops gradients in the GibbsThomsonHerring interface temperature distribution that provide tangential energy fluxes. The LeibnizReynolds transport theorem exposes a 4thorder, netzero energy field (the 'Bias' field) that autonomously deposits and removes capillarymediated thermal energy. Where energy is released locally, the freezing rate is persistently retarded, and where energy is removed, the rate is enhanced. These contravening dynamic field responses balance at points (roots) where the surface Laplacian of the chemical potential vanishes, inducing an inflection, or ‘curling’, of the interface. Interfacial inflection couples to the main transport fields producing pattern branching, folding, and complexity.
Precision noisefree numerical schemes, including integral equation sharpinterface solvers (J. Lowengrub, S. Li) and, recently, three noisefree phasefield simulations (A. Mullis, M. Zaeem, K. Reuther) independently confirm that pattern branching initiates at locations predicted using analytical methods for smooth, noisefree starting shapes in 2D. A limit cycle may develop as the interface and its energy field coevolve, synchronizing the inflection points to produce classical dendritic structures. Noise and stochastics play no direct role in the proposed deterministic mechanism of branching and pattern morphogenesis induced by persistent 'perturbations'. 
INI 1  
12:00 to 12:30 
CH Cheng (National Central University) Some moving boundary value problems consisting of viscous incompressible fluids moving and interacting with nonlinear elastic shells Session: Focus Session Recent Developments 3 (parallel session) Chair: D Danielli
In this talk, I will briefly talk about my past works on the interaction of incompressible NavierStokes with either the biofluid shells or the Koiter shells. The biofluid shells can be used to model the cell membranes, while the Koiter shells are used to model general elastic thin shells that have equilibrium state. The dynamics of both kinds of shells are described by nonlinear PDEs derived from energy functionals depending on geometric structures of the shells. The wellposedness of the problems concerning the coupling of the fluids and the shells are nontrivial, and we established the existence and uniqueness of solutions in Sobolev spaces.

INI 2  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:15 
AL Bertozzi (University of California, Los Angeles) Plenary Lecture 9: Curvature flow on graphs for large data classification Chair: JF Rodrigues
In the continuum, close connections exist between mean curvature ow, the AllenCahn (AC) partial di erential equation, and the MerrimanBenceOsher (MBO) threshold dynamics scheme. Graph analogues of these processes have recently seen a rise in popularity as relaxations of NPcomplete combinatorial problems, which demands deeper theoretical underpinnings of the graph processes. We discuss several applications including supervised and unsupervised machine learning and community detection in social networks. We discuss connections to spectral graph theory and fast algorithms and some recent results for curvature flow on graphs and open problems.

INI 1  
14:20 to 14:50  Afternoon Tea  
14:50 to 15:20 
M Chipot ([University of Zürich]) Obstacle problems in unbounded domains and related issues Session: Focus Session Analysis 2 (parallel session) Chair: JL Vázquez
We will present a formulation of obstacle problems in unbounded domains when the energy method does not work i.e. when the force does not belong to H1(O). We will report also on some new results for the pure Neumann problem in unbounded cylinders.

INI 1  
14:50 to 15:20 
Oscillating facets Session: Focus Session Recent Developments 4 (parallel session) Chair: A Karakhanyan
Coauthor: Piotr Rybka (The University of Warsaw)
We study a very singular onedimensional parabolic problem with initial data in the energy space. We study facet creation and extinction of solutions caused by the evolution of facets. 
INI 2  
15:20 to 15:50 
C Heinemann (WeierstraßInstitut für Angewandte Analysis und Stochastik) Analysis of coupled phasefield models describing damage phenomena Session: Focus Session Recent Developments 4 (parallel session) Chair: A Karakhanyan
Coauthors: Christiane Kraus (WIAS Berlin), Elisabetta Rocca (WIAS Berlin / University of Pavia)
In this talk we are going to investigate evolution inclusions describing damage processes in elastic media which are coupled with further physical phenomena such as heat conduction or phase separation. The main difficulty in the analysis of these PDE systems is the irreversibility condition for the damage progression. To handle this highly nonlinear behavior we develop an appropriate notion of weak solution consisting of a variational inequality and a full energy inequality. Results ensuring the existence of weak solutions will be presented. The talk is based on joint works with C. Kraus (WIAS) and E. Rocca (WIAS/University of Pavia). 
INI 2  
15:20 to 15:50 
Pointwise estimates and regularity for the parabolic obstacle problem Session: Focus Session Analysis 2 (parallel session) Chair: JL Vázquez
Coauthor: Régis Monneau (ENPC)
I will discuss pointwise regularity properties for parabolic equations. In particular, I will talk about a method that can be used to obtain a second order Taylor expansion. I will also discuss how the same method can be applied to the parabolic obstacle problem in order to study the regularity of the free boundary. 
INI 1  
16:00 to 16:30 
Free boundaries in fractional filtration equations Session: Focus Session Analysis 2 (parallel session) Chair: JL Vázquez
Coauthors: Arturo de Pablo (U. Carlos III de Madrid), Ana Rodríguez (U. Politécnica de Madrid), Juan Luis Vázquez (U. Autónoma de Madrid)
In this talk we will review some recent results on the Cauchy problem for the fractional filtration equation $\partial_t u+(\Delta)^{\sigma/2}\varphi(u)$, $\sigma\in(0,2)$, where the nonlinearity $\varphi$ satisfies some not very restrictive conditions. Solutions to these problems become instantaneously positive if the initial data are nonnegative. However, a free boundary emerges in certain cases in which there is a distinguished value, as in the fractional Stefan problem, $\varphi(u)=(u1)_+$, or in the ‘mesa’ limit, $m\to\infty$, for the fractional porous media equation, $\varphi(u)=u^m$. 
INI 1  
16:00 to 16:30 
S Mitchell (University of Limerick) Numerical and asymptotic solutions of vertical continuous casting with and without superheat Session: Focus Session Recent Developments 4 (parallel session) Chair: A Karakhanyan
Coauthor: Michael Vynnycky (Royal Institute of Technology (Sweden))
In a continuous casting process, such as the strip casting of copper, molten metal first passes through a watercooled mould region, before being subjected to a high cooling rate further downstream. Consequently, the molten metal solidifies and the solidified metal is withdrawn at a uniform casting speed. Industrialists need to understand the factors influencing product quality and process productivity. Of key significance is the heat transfer that occurs during solidification, particularly the location of the interface between molten metal and solid. The modelling of the continuous casting of metals is known to involve the complex interaction of nonisothermal fluid and solid mechanics. Typically, the flow in the molten metal is turbulent, and it is generally believed that a computational fluid dynamics (CFD) approach is necessary in order to correctly capture the heat transfer characteristics. However, we can show that an asymptotically reduced version of the CFDbased model, which neglects this turbulence, gives predictions for the pool depth, local temperature profiles and mould wall heat flux that agree very well with results of the original CFD model. This reduced model can be described as a steady state 2D heat flow Stefan problem, with a degenerate initial condition and nonstandard Neumanntype boundary condition. If we assume the incoming metal is at the melt temperature then we obtain a onephase model but with potentially two stages, depending whether the metal is fully solidified before leaving the mould. However, in reality, the incoming temperature is greater than the melt temperature, termed as including superheat, and this leads to a twophase model with a presolidification stage, where the second phase only first appears after a finite delay. In this work we highlight some numerical challenges in solving the systems with and without superheat. The Keller box finitedifference scheme is used, along with a boundary immobilisation method. 
INI 2  
16:30 to 17:00 
On a thermodynamically consistent Stefan problem with variable surface energy Session: Focus Session Analysis 2 (parallel session) Chair: JL Vázquez
A thermodynamically consistent twophase Stefan problem with temperature dependent surface tension is studied. It is shown that this problem generates a local semiflow on a welldefined state manifold. Moreover, stability and instability results of equilibrium configurations will be presented. It will be pointed out that surface heat capacity has a striking effect on the stability behavior of multiple equilibria.

INI 1  
16:30 to 17:00 
Higher order convergent trial methods for Bernoulli's free boundary problem Session: Focus Session Recent Developments 4 (parallel session) Chair: A Karakhanyan
Coauthor: Helmut Harbrecht (University of Basel)
Free boundary problem is a partial differential equation to be solved in a domain, a part of whose boundary is unknown – the socalled free boundary. Beside the standard boundary conditions that are needed in order to solve the partial differential equation, an additional boundary condition is imposed at the free boundary. One aims thus to determine both, the free boundary and the solution of the partial differential equation. This work is dedicated to the solution of the generalized exterior Bernoulli free boundary problem which is an important model problem for developing algorithms in a broad band of applications such as optimal design, fluid dynamics, electromagnentic shaping etc. For its solution the trial method, which is a fixedpoint type iteration method, has been chosen. The iterative scheme starts with an initial guess of the free boundary. Given one boundary condition at the free boundary, the boundary element method is applied to compute an approximation of the violated boundary data. The free boundary is then updated such that the violated boundary condition is satisfied at the new boundary. Taylor's expansion of the violated boundary data around the actual boundary yields the underlying equation, which is formulated as an optimization problem for the sought update function. When a target tolerance is achieved, the iterative procedure stops and the approximate solution of the free boundary problem is detected. The efficiency of the trial method as well as its speed of convergence depends significantly on the update rule for the free boundary, and thus on the violated boundary condition. This talk focuses on the trial method with violated Dirichlet boundary data and on the development of higher order convergent versions of the trial method with the help of shape sensitivity analysis. 
INI 2 
09:00 to 09:45 
Plenary Lecture 10: Absence of the interface splash singularity for the twofluid Euler equations Chair: M Feldman
An interface splash singularity occurs when a locally smooth fluid interface selfintersects. Such interface singularities occur for onefluid interfaces in the Euler equations and other fluids models.
By means of elementary arguments in Lagrangian coordinates, we prove that such a singularity cannot occur in finitetime for a twofluid interface evolved by either the incompressible Euler equations (with surface tension) or the Muskat equations. By assuming that such a singularity can occur, we find a sharp blowup rate for the vorticity, and characterize the geometry of the evolving interface. This leads to a contradiction, showing that such a singularity can occur. This is joint work with D. Coutand. 
INI 1  
10:00 to 10:45 
J Lowengrub (University of California, Irvine) Plenary Lecture 11: Numerical simulation of endocytosis Chair: M Feldman
Coauthors: Sebastian Aland (UC Irvine/TU Dresden), Jun Allard (UC Irvine)
Many cell processes involve the formation of membrane vesicles from a larger membrane, including endocytosis, interorganelle transport and virus entry. These events are typically orchestrated by curvatureinducing molecules attached to the membrane, such as clathrin and bardomain proteins. Recent reports demonstrate that in some circumstances vesicles can form de novo in a few milliseconds, e.g., ultrafast endocytosis at the neurological synapse. Membrane dynamics at these scales (millisecond, nanometer) are dominated by hydrodynamic interactions, as the membrane pushes the intracellular and extracellular fluids around to accommodate curvature. To study this problem, we develop new diffuse interface models for the dynamics of inextensible vesicles in a viscous fluid with stiff, curvatureinducing molecules. A new feature of this work is the implementation of the local inextensibility condition by using a local Lagrange multiplier harmonically extended off the interface. To make the method even more robust, we develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. This is critical to accurately capturing hydrodynamic effects during endocytosis. By varying the membrane coverage of curvatureinducing molecules, we find that there is a cri tical (smallest) neck radius and a critical (fastest) budding time. 
INI 1  
11:00 to 11:30  Morning Coffee  
11:30 to 12:15 
O Jensen (University of Manchester) Plenary Lecture 12: Slamming in flexiblechannel flows Chair: J King
Coauthor: Feng Xu (University of Nottingham)
Largeamplitude selfexcited oscillations of highReynoldsnumber flow in a long flexiblewalled channel can exhibit vigorous slamming motion, whereby the channel is almost completely occluded over a very short interval in space and time. Treating the flexible channel wall as an inertialess elastic membrane, this nearsingular behaviour is exhibited in twodimensional NavierStokes simulations and can be captured in a reduced onedimensional PDE model (Stewart et al. J. Fluid Mech. 662:288, 2010). The properties of the rigid parts of the system, upstream and downstream of the membrane, play a major role in determining the onset of oscillations. In order to investigate the extreme flow structure that arises during a brief slamming event, we systematically reduce the PDE model to a thirdorder nonlinear algebraicdifferential system, which identifies the likely dominant physical balances. 
INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:15 
M Feldman (University of WisconsinMadison) Plenary Lecture 13: Shock Reflection, von Neumann conjectures, and free boundary problems Chair: C Elliott
Coauthor: GuiQiang Chen (Oxford)
We discuss shock reflection problem for compressible gas dynamics, and von Neumann conjectures on transition between regular and Mach reflections. Then we will talk about recent results on existence of regular reflection solutions for potential flow equation up to the detachment angle, and discuss some techniques. The approach is to reduce the shock reflection problem to a free boundary problem for a nonlinear equation of mixed elliptichyperbolic type. Open problems will also be discussed. 
INI 1  
14:30 to 15:00  Afternoon Tea  
15:00 to 15:45 
Plenary Lecture 14: Free boundary problems for mechanical models of tumor growth Chair: C Elliott
Mathematical models of tumor growth, now commonly used, present several
levels of complexity, both in terms of the biomedical ingredients and the
mathematical description. The simplest ones contain competition for space
using purely fluid mechanical concepts. Another possible ingredient is the
supply of nutrients. The models can describe the tissue either at the level
of cell densities, or at the scale of the solid tumor, in this latter case
by means of a free boundary problem.
We first formulate a free boundary model of HeleShaw type, a variant
including growth terms, starting from the description at the cell level and
passing to a certain singular limit which leads to a HeleShaw type problem.
A detailed mathematical analysis of this purely mechanical model is
performed. Indeed, we are able to prove strong convergence in passing to the
limit, with various uniform gradient estimates; we also prove uniqueness for
the limit problem. At variance with the classical HeleShaw problem, here
the geometric motion governed by the pressure is not sufficient to
completely describe the dynamics.
Using this theory as a basis, we go on to consider a more complex model
including nutrients. Here, technical difficulties appear, that reduce the
generality and detail of the description.
We prove uniqueness for the system, a main mathematical difficulty.
Joint work with Benoit Perthame, Paris, and Fernando Quiros, Madrid

INI 1  
15:45 to 16:00  Closing Remarks  INI 1 