Calculus of Variations and its Applications
Org:
Andrew Colinet (University of Toronto) and
Dominik Stantejsky (McMaster)
[
PDF]
 BLAISE BOURDIN, McMaster University
Phasefield approximation of diffusiondriven fracture [PDF]

This talk focusses on a class of problems where crack propagation is driven by a diffusion process. This general framework encompasses a broad range of phenomena including thermal and dessication cracks, fracture in fluidsaturated porous media, or fracture of materials undergoing phasechange.
The main difficulty in building a rigorous phasefield model of such problems is the different time scales involved in the fracture and diffusion processes. The former is often assumed to remain in an equilibrium state at all time, whereas the later is inherently driven by its nonequilibrium nature.
$\Gamma$convergence, commonly used to derive the convergence of phasefield models to their sharp interface counterparts do not provide any insight on outofequilibrium evolutions. Instead, we propose to reformulate the diffusion problem in terms of the minimizing motion of an energy, and the coupled problem as a PDEconstrained optimization problem. We then propose compatible phasefield approximations of the fracture and diffusion process can be derived, and the convergence of the constrained minimization problem can be proved.
 DENIS BRAZKE, University of Heidelberg
$\Gamma$–limit for a sharp interface model related to pattern formation on biomembranes [PDF]

We derive a macroscopic limit for a sharp interface version of a model proposed by Komura, Shimokawa and Andelman to investigate pattern formation in biomembranes due to competition of chemical and mechanical forces. We identify sub and supercrital parameter regimes and show with the introduction of the autocorrelation function that the ground state energy leads to the isoperimetric problem in the subcritical regime, which is interpreted to not form fine scale patterns
This is joint work with Hans Knüpfer and Anna MarciniakCzochra.
 ALMUT BURCHARD, University of Toronto
Symmetrybreaking in isodiametric capacitor problems [PDF]

A classical theorem of Szegö says that balls maximize electrostatic capacity among sets of given diameter (while minimizing capacity among sets of given volume). On the other hand,
balls do not maximize Rieszcapacity (with a Rieszpotential $x^{\lambda}$,
for fixed $\lambda>0$) among sets of given diameter in high dimensions. Thus symmetrybreaking
occurs as the Riesz kernel transitions from the Newtonian case $(\lambda=n2$) to the logarithmic case (corresponding to $\lambda=0$), once $\lambda$ is small enough relative to $n$. (Joint work with R. Choksi, E. HessChilds, and A. Martìnez.)
 CARRIE CLARK, University of Illinois UrbanaChampaign
Droplet formation in a nonlocal aggregation model [PDF]

The study of aggregation in the physical sciences has produced a rich class of nonlocal shape optimization problems. In this talk, we will discuss droplet formation in energy minimizing configurations for a family of interaction kernels which have a "wellbarrier" type shape. Short distance attraction, combined with mid distance repulsion, and long distance neutrality drives the separation into droplets.
 ANDREW COLINET, McMaster University
Zeroth Order Limiting Behaviour of the GinzburgLandau Functional [PDF]

For $\Omega\subseteq\mathbb{R}^{2}$ there is an extensive literature concerning the limiting
behaviour of the GinzburgLandau energy,
\begin{equation*}
E_{\varepsilon}(u)=\int_{\Omega}\!{}\Bigl\{\frac{1}{2}\nabla{}u^{2}
+\frac{1}{4\varepsilon^{2}}(u^{2}1)^{2}\Bigr\},
\end{equation*}
as $\varepsilon\to0^{+}$.
In such works, it is shown that if a sequence of functions has vorticity concentrate, as
$\varepsilon\to0^{+}$, about a finite collection of interior points of $\Omega$ then the
GinzburgLandau energy converges, after renormalizing, to the total variation of a measure
supported over the same interior points.
However, much less is known when the vorticity of solutions is permitted to concentrate about
points along the boundary.
We consider this question for a connected open set $\Omega\subseteq\mathbb{R}^{2}$ with
$C^{2,1}$ boundary and we prove that similar conclusions to the interior case remain true
up to the boundary provided the functions we consider satisfy suitable boundary restrictions.
In addition, we also show that there are necessary topological restrictions on the vorticity.
 ROBERT HASLHOFER, University of Toronto
Classification of compact ancient noncollapsed flows in $R^4$ [PDF]

To capture singularities under mean curvature flow one wants to understand all ancient solutions. In addition to shrinkers and translators one also encounters ancient ovals, namely compact noncollapsed solutions that are not selfsimilar. In this talk, I will explain that any bubblesheet oval for the mean curvature flow in $R^4$, up to scaling and rigid motion, either is the $O(2)\times O(2)$symmetric ancient oval constructed by White, or belongs to the oneparameter family of $Z_2^2 \times O(2)$symmetric ancient ovals constructed by Du and myself. In particular, this seems to be the first instance of a classification result for geometric flows that are neither cohomogeneityone nor selfsimilar. This is joint work with Beomjun Choi, Toti Daskalopoulos, Wenkui Du and Natasa Sesum. I will also briefly mention the noncompact case, which is joint work with Kyeongsu Choi and Or Hershkovits.
 NICHOLAS KEVLAHAN, McMaster University
Data assimilation for bathymetry in the nonlinear shallow water equations [PDF]

The shallow water equations (SWE) are a widely used model for the propagation of surface waves in oceans, lakes and rivers. Common applications include modelling the propagation of tsunami waves, storm surges and flooding. We consider the problem of determining under which conditions a spacetime variational data assimilation approach based on observations of the free surface is able to reconstruct the bathymetry to a given accuracy (e.g. sufficient for modelling wave propagation). We then use densitybased global sensitivity analysis (GSA) to assess the sensitivity of the surface wave and reconstruction error to model parameters and second order adjoint analysis (SOA) to analyze the sensitivity of the surface wave error, given the reconstructed bathymetry, to perturbations in the observations.
This is joint work with Bartek Protas (McMaster University) and Ramsha Khan (University of Stockholm)
 SULLIVAN MACDONALD, McMaster University
Degenerate Ellipticity and Hypoellipticity for Divergence Operators [PDF]

We discuss differential operators of the form \(L=\mathrm{div}(Q\nabla \cdot)\), where \(Q\) is a nonnegative definite symmetric matrixvalued function. Even though \(L\) fails to be elliptic at points where \(Q\) is singular, in some settings we can recover regularity results for weak solutions to the Dirichlet problem
\begin{align*}
Lu=f & \quad x\in\Omega, \ u=0 & \quad x\in\partial\Omega.
\end{align*}
In this talk, I will present some recent joint work with S. Rodney concerning sufficient conditions for a priori boundedness of weak solutions to these problems. We show that even if a very weak Sobolev inequality with gain on Orlicz scale holds, one can still recover boundedness under appropriate hypotheses on the data \(f\).
On the other hand, if rapid degeneracy of \(Q\) forbids a Sobolev inequality with gain on Orlicz scale from holding, it would still be useful to understand qualitative behaviours of the operator \(L\). To this end, I also discuss recent work at finding sufficient conditions for hypoellipticity of secondorder divergence operators by using decompositions of Holder continuous functions into sums of squares. This work complements some recent results by Sawyer and Korobenko, which show that \(L=\mathrm{div}(Q\nabla\cdot)\) is hypoelliptic when the matrix \(Q\) admits a suitable decomposition.
 DMITRY PELINOVKY, McMaster University
 BARTEK PROTAS, McMaster University
Searching for Singularities in NavierStokes Flows Using Variational Optimization Methods [PDF]

This investigation concerns a systematic computational search for potentially singular behavior in 3D NavierStokes flows. Enstrophy $\mathcal{E}(t)$ serves as a convenient indicator of the regularity of solutions  as long as this quantity remains finite, the solutions are guaranteed to be smooth and satisfy the equations in the classical sense. Another wellknown conditional regularity result are the LadyzhenskayaProdiSerrin conditions asserting that the quantity $\mathcal{L}_{q,p}:=\int_0^T \ \mathbf{u}(t) \_{L^q(\Omega)}^p \, dt$, where $2/p+3/q \le 1$, $q > 3$, must remain bounded if the solution is smooth on the interval $[0,T]$. However, there are no finite a priori bounds available for these quantities and hence the regularity problem for the 3D NavierStokes system remains open. To quantify the maximum possible growth of $\mathcal{E}(T)$ and $\mathcal{L}_{q,p}$, we consider families of variational PDE optimization problems in which initial conditions are sought subject to certain constraints so that these quantities in the resulting NavierStokes flows are maximized. These problems are solved computationally using a largescale adjointbased gradient approach. By solving these problems for a broad range of parameter values we demonstrate that the maximum growth of $\mathcal{E}(T)$ and $\mathcal{L}_{q,p}$ appears finite. Thus, in the worstcase scenarios the two quantities remain bounded for all times and there is no evidence for singularity formation in finite time.
[Joint work with Dongfang Yun and Di Kang]
 IVAN SALGADO, University of Toronto
Approximate Solutions to the Superconducting Interface Model [PDF]

The superconducting interface model is a semilinear hyperbolic system of PDEs introduced in 2016 by Kyle Thompson. It proposes a more tractable, yet closely related alternative to a 1984 model of Edward Witten for cosmic strings carrying a superconducting current. In the superconducting interface model, we consider the system
\begin{equation*}
\begin{cases}
\epsilon^2(\partial_t^2\varphi  \Delta_x \varphi) + \lambda_\varphi(\varphi^21)\varphi + \beta \varphi\sigma^2=0
\\epsilon^2(\partial_t^2\sigma \Delta_x \sigma) + \lambda_\sigma(\sigma^2m_\sigma^2)\sigma + \beta\varphi^2\sigma=0
\end{cases},
\end{equation*}
where $0 < \epsilon\ll 1$, $(\lambda_\varphi, \lambda_\sigma,m_\sigma,\beta)\in (0,\infty)^4$ are parameters, and
\begin{align*}
\varphi &:[0,T] \times \mathbb{R}^n\to \mathbb{R} \text{ and } \sigma :[0,T]\times \mathbb{R}^n\to \mathbb{C},
\end{align*}
for some $T>0$ and $n\in \mathbb{N}$ with $n\geq2$. We are interested in solutions $(\varphi,\sigma)$ such that
\begin{align*}
\varphi \approx
\begin{cases}
+1 & \text{ in } \mathcal{O}^+
\1 & \text{ in } \mathcal{O}^
\end{cases},
\end{align*}
where $\mathcal{O}^+$ and $\mathcal{O}^$ are disjoint open subsets of $[0,T]\times \mathbb{R}^n$ separated by an ``interface'' with thickness of order $\epsilon$, and the currentcarrying field $\sigma$ decays exponentially away from the interface. The problem is of particular interest when $\sigma$ showcases a clear interaction with the geometry of the interface.
The purpose of this talk is to present a methodology for finding solutions to the superconducting interface model by first constructing approximate solutions, and then linearizing the system of PDEs around these approximations. We will talk about some of the results regarding the construction of the approximate solutions and the laws of motion which represent the coupling between the current supported around the interface and the geometry of the interface.
 DOMINIK STANTEJSKY, McMaster University
A finite element approach for minimizing line and surface energies arising in the study of singularities in liquid crystals [PDF]

I will present an algorithm designed to calculate minimizers $T$ of a geometric energy arising in the theory of liquid crystal colloids. The energy involves the two dimensional area of $T$ outside an obstacle, a contribution from $T$ on the obstacle surface, and the length of the boundary $\partial T$ reduced by a prescribed curve to make the problem nontrivial. It can be seen as a generalization of both the obstacle and Plateau problem. We discretize the energy by a finite element method and apply an ADMM scheme to carry out the minimization. We validate our algorithm in the case of a spherical obstacle and give examples of minimizing configurations in the case of a peanut and croissantshaped obstacle.
 ZHICHAO WANG, The University of British Columbia
Minmax minimal hypersurfaces with higher multiplicity [PDF]

Recently, X. Zhou proved that the AlmgrenPitts minmax solution has multiplicity one for bumpy metrics (Multiplicity One Theorem). In this talk, we exhibit the first set of examples of nonbumpy metrics on the $(n+1)$sphere ($2\leq n\leq 6$) in which the varifold associated with the twoparameter minmax construction must be a multiplicitytwo minimal $n$sphere. This is proved by a new areaandseparation estimate for certain minimal hypersurfaces with Morse index two inspired by an early work of ColdingMinicozzi. This is a joint work with X. Zhou.
© Canadian Mathematical Society