Nikos Katzourakis

Nikos Katzourakis, Tristan Pryer

Given a map $u : Ω\subseteq \mathbb{R}^n \longrightarrow \mathbb{R}^N$, the $\infty$-Laplacian is the system \[ \label{1} Δ_\infty u \, :=\, \Big(\text{D}u \otimes \text{D}u + |\text{D}u|^2 [\text{D}u]^\bot \! \otimes I \Big) : \text{D}^2 u\, = \, 0. \tag{1} \] \eqref{1} is the model system of vectorial Calculus of Variations in $L^\infty$ and arises as the "Euler-Lagrange" equation in relation to the supremal functional \[ \label{2} E_\infty(u,Ω)\, :=\, \| \text{D}u \|_{L^\infty(Ω)}. \tag{2} \] The scalar case of \eqref{1} has been introduced by Aronsson in the 1960s and by now is relatively classical and well understood. The general system \eqref{1} has been discovered and studied by the first author in a series of recent papers. Supremal functionals are fundamental for applications because they provide more realistic models as opposed to conventional integral models. Herein we provide numerical approximations of solutions to the Dirichlet problem when $n=2$ and $N=2,3$ for certain carefully selected boundary data on the unit square. Our experiments demonstrate interesting and unexpected phenomena occurring in $L^\infty$ and provide insights on the structure of general solutions and the natural separation to phases they present.

Nikos Katzourakis, Tristan Pryer

In [KP16] (arXiv:1605.07880) the authors introduced a second-order variational problem in $L^{\infty}$. The associated equation, coined the $\infty$-Bilaplacian, is a \emph{third order} fully nonlinear PDE given by $Δ^2_\infty u\, := (Δu)^3 | D (Δu) |^2 = 0.$ In this work we build a numerical method aimed at quantifying the nature of solutions to this problem which we call $\infty$-Biharmonic functions. For fixed $p$ we design a mixed finite element scheme for the pre-limiting equation, the $p$-Bilaplacian $Δ^2_p u\, := Δ(| Δu |^{p-2} Δu) = 0.$ We prove convergence of the numerical solution to the weak solution of $Δ^2_p u = 0$ and show that we are able to pass to the limit $p\to\infty$. We perform various tests aimed at understanding the nature of solutions of $Δ^2_\infty u$ and in 1-$d$ we prove convergence of our discretisation to an appropriate weak solution concept of this problem, that of $\mathcal D$-solutions.

Nikos Katzourakis, Tristan Pryer

In this paper we initiate the study of $2$nd order variational problems in $L^\infty$, seeking to minimise the $L^\infty$ norm of a function of the hessian. We also derive and study the respective PDE arising as the analogue of the Euler-Lagrange equation. Given $\mathrm{H}\in C^1(\mathbb{R}^{n\times n}_s)$, for the functional \[ \label{1} \mathrm{E}_\infty(u,\mathcal{O})\, =\, \big\| \mathrm{H}\big(\mathrm{D}^2 u\big) \big\|_{L^\infty(\mathcal{O})}, \ \ \ u\in W^{2,\infty}(Ω),\ \mathcal{O}\subseteq Ω, \tag{1} \] the associated equation is the fully nonlinear 3rd order PDE \[ \label{2} \mathrm{A}^2_\infty u\, :=\,\big(\mathrm{H}_X\big(\mathrm{D}^2u\big)\big)^{\otimes 3}:\big(\mathrm{D}^3u\big)^{\otimes 2}\, =\,0. \tag{2} \] Special cases arise when $\mathrm{H}$ is the Euclidean length of either the full hessian or of the Laplacian, leading to the $\infty$-Polylaplacian and the $\infty$-Bilaplacian respectively. We establish several results for \eqref{1} and \eqref{2}, including existence of minimisers, of absolute minimisers and of "critical point" generalised solutions, proving also variational characterisations and uniqueness. We also construct explicit generalised solutions and perform numerical experiments.

Nikos Katzourakis, Tristan Pryer

Let $Ω$ be an open set. We consider the supremal functional \[ \tag{1} \label{1} \ \ \ \ \ \ \mathrm{E}_\infty (u,\mathcal{O})\, :=\, \| \mathrm D u \|_{L^\infty( \mathcal{O} )}, \ \ \ \mathcal{O} \subseteq Ω\text{ open}, \] applied to locally Lipschitz mappings $u : \mathbb R^n \supseteq Ω\longrightarrow \mathbb R^N$, where $n,N\in \mathbb N$. This is the model functional of Calculus of Variations in $L^\infty$. The area is developing rapidly, but the vectorial case of $N\geq 2$ is still poorly understood. Due to the non-local nature of \eqref{1}, usual minimisers are not truly optimal. The concept of so-called absolute minimisers is the primary contender in the direction of variational concepts. However, these cannot be obtained by direct minimisation and the question of their existence under prescribed boundary data is open when $n,N\geq 2$. Herein we present numerical experiments based on a new method recently proposed by the first author in the papers [33, 35].