Time integration and steady-state continuation for 2d lubrication equations

Lubrication equations allow to describe many structurin processes of thin liquid films. We develop and apply numerical tools suitable for their analysis employing a dynamical systems approach. In particular, we present a time integration algorithm based on exponential propagation and an algorithm for steady-state continuation. In both algorithms a Cayley transform is employed to overcome numerical problems resulting from scale separation in space and time. An adaptive time-step allows to study the dynamics close to hetero- or homoclinic connections. The developed framework is employed on the one hand to analyse different phases of the dewetting of a liquid film on a horizontal homogeneous substrate. On the other hand, we consider the depinning of drops pinned by a wettability defect. Time-stepping and path-following are used in both cases to analyse steady-state solutions and their bifurcations as well as dynamic processes on short and long time-scales. Both examples are treated for two- and three-dimensional physical settings and prove that the developed algorithms are reliable and efficient for 1d and 2d lubrication equations, respectively.Comment: 33 pages, 16 figure

On the numerical approximation of p-biharmonic and ∞-biharmonic functions

In [KP16] (arXiv:1605.07880) the authors introduced a second-order variational problem in L∞. The associated equation, coined the ∞-Bilaplacian, is a \emph{third order} fully nonlinear PDE given by Δ2∞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 ∞-Biharmonic functions. For fixed p we design a mixed finite element scheme for the pre-limiting equation, the p-Bilaplacian Δ2pu:=Δ(|Δu|p−2Δu)=0. We prove convergence of the numerical solution to the weak solution of Δ2pu=0 and show that we are able to pass to the limit p→∞. We perform various tests aimed at understanding the nature of solutions of Δ2∞u and in 1-d we prove convergence of our discretisation to an appropriate weak solution concept of this problem, that of -solutions

The eigenvalue problem for the ∞-Bilaplacian

We consider the problem of finding and describing minimisers of the Rayleigh quotient Λ∞:=infu∈W2,∞(Ω)∖{0}∥Δu∥L∞(Ω)∥u∥L∞(Ω), Λ∞:=infu∈W2,∞(Ω)∖{0}‖Δu‖L∞(Ω)‖u‖L∞(Ω), where Ω⊆RnΩ⊆Rn is a bounded C1,1C1,1 domain and W2,∞(Ω)W2,∞(Ω) is a class of weakly twice differentiable functions satisfying either u=0u=0 on ∂Ω∂Ω or u=|Du|=0u=|Du|=0 on ∂Ω∂Ω . Our first main result, obtained through approximation by LpLp -problems as p→∞p→∞ , is the existence of a minimiser u∞∈W2,∞(Ω)u∞∈W2,∞(Ω) satisfying {Δu∞∈Λ∞Sgn(f∞)Δf∞=μ∞ a.e. in Ω, in D′(Ω), {Δu∞∈Λ∞Sgn(f∞) a.e. in Ω,Δf∞=μ∞ in D′(Ω), for some f∞∈L1(Ω)∩BVloc(Ω)f∞∈L1(Ω)∩BVloc(Ω) and a measure μ∞∈M(Ω)μ∞∈M(Ω) , for either choice of boundary conditions. Here Sgn is the multi-valued sign function. We also study the dependence of the eigenvalue Λ∞Λ∞ on the domain, establishing the validity of a Faber–Krahn type inequality: among all C1,1C1,1 domains with fixed measure, the ball is a strict minimiser of Ω↦Λ∞(Ω)Ω↦Λ∞(Ω) . This result is shown to hold true for either choice of boundary conditions and in every dimension

Nonexistence results for nonlocal equations with critical and supercritical nonlinearities

We prove nonexistence of nontrivial bounded solutions to some nonlinear problems involving nonlocal operators of the form $$Lu(x)=\sum a_{ij}\partial_{ij}u+{\rm PV}\int_{\R^n}(u(x)-u(x+y))K(y)dy.$$ These operators are infinitesimal generators of symmetric L\'evy processes. Our results apply to even kernels $K$ satisfying that $K(y)|y|^{n+\sigma}$ is nondecreasing along rays from the origin, for some $\sigma\in(0,2)$ in case $a_{ij}\equiv0$ and for $\sigma=2$ in case that $(a_{ij})$ is a positive definite symmetric matrix. Our nonexistence results concern Dirichlet problems for $L$ in star-shaped domains with critical and supercritical nonlinearities (where the criticality condition is in relation to $n$ and $\sigma$). We also establish nonexistence of bounded solutions to semilinear equations involving other nonlocal operators such as the higher order fractional Laplacian $(-\Delta)^s$ (here $s>1$) or the fractional $p$-Laplacian. All these nonexistence results follow from a general variational inequality in the spirit of a classical identity by Pucci and Serrin