---
title: Willmore Equation for Graphs
url: https://www.emergentmind.com/topics/willmore-equation-for-graphs
type: topic
---

# Willmore Equation for Graphs

Searching arXiv for recent and foundational papers on the Willmore equation for graphs.
The Willmore equation for graphs is the Euler–Lagrange equation associated with the Willmore functional when a surface in $\mathbb{R}^3$ is represented nonparametrically as the graph of a height function. In the graphical setting, one studies surfaces of the form $\Sigma=\{(x,y,u(x,y)):(x,y)\in\Omega\}$, with $\Omega\subset\mathbb{R}^2$ or $\Omega=\mathbb{R}^2$, and rewrites the geometric equation in terms of $u$, its gradient, and its Hessian. This produces a quasilinear fourth-order PDE whose structure can be expressed either intrinsically through the induced metric and the Laplace–Beltrami operator or explicitly in divergence form. The subject connects geometric analysis, elliptic and parabolic fourth-order PDE, symmetry reduction, Bernstein-type rigidity, boundary value problems, variational relaxation, and numerical approximation [1110.3221], [2509.21018], [2603.27848], [1410.5547], [1503.01275], [1111.3043].

## 1. Geometric formulation for graphical surfaces

Let $\Sigma\subset\mathbb{R}^3$ be represented as a graph
\[
\Sigma=\{(x,y,u(x,y)):(x,y)\in\Omega\},
\]
with $u$ smooth. In the notation used across the cited works, one writes $Du=(u_x,u_y)$ and
\[
v=\sqrt{1+|Du|^2},\qquad Q=\sqrt{1+|Du|^2},
\]
so that $v=Q$ [1110.3221], [2509.21018], [2603.27848], [1503.01275]. The induced metric is
\[
g_{ij}=\delta_{ij}+u_i u_j,
\]
with inverse
\[
g^{ij}=\delta^{ij}-\frac{u_i u_j}{v^2}
\]
or equivalently
\[
g^{ij}=\delta_{ij}-\frac{u_i u_j}{Q^2},
\]
and the area element is
\[
d\mu_g=v\,dx\,dy,\qquad dS=Q\,dx,
\]
depending on the notation of the source [1110.3221], [2509.21018], [2603.27848], [1503.01275]. With the upward orientation, the unit normal is
\[
\nu=n=\frac{1}{Q}(-u_1,-u_2,1)=\frac{1}{v}(-Du,1),
\]
and for an entire graph the Gauss map image lies in the upper hemisphere $S^2_+$ [1110.3221], [2509.21018], [2603.27848], [1503.01275].

The second fundamental form coefficients are
\[
h_{ij}=A_{ij}=\frac{u_{ij}}{Q}=\frac{u_{ij}}{v},
\]
and the mean curvature is taken with the convention
\[
H=k_1+k_2,
\]
so minimal graphs satisfy $H=0$ [1110.3221]. In graph variables,
\[
H=\operatorname{div}\!\left(\frac{\nabla u}{Q}\right)
=\frac{\Delta u}{Q}-\frac{\nabla u\cdot(D^2u\,\nabla u)}{Q^3},
\]
and the Gauss curvature is
\[
K=\frac{\det D^2u}{Q^4}=\frac{\det\nabla^2u}{v^4}
\]
[1110.3221], [2509.21018], [2603.27848], [1503.01275]. The squared norm of the second fundamental form satisfies the Gauss equation
\[
|A|^2=H^2-2K
\]
in the convention of Luo and Sun [1110.3221]. The trace-free second fundamental form $A^o$ is not used directly in that note, but the standard decomposition recorded there is
\[
|A|^2=|A^o|^2+\frac12 H^2,
\]
hence
\[
|A^o|^2=|A|^2-\frac12 H^2=\frac{H^2}{2}-2K
\]
[1110.3221].

For a scalar function $f$ on the graph, the Laplace–Beltrami operator is
\[
\Delta_g f=\frac{1}{Q}\,\operatorname{div}\!\left(Q\left(I-\frac{\nabla u\otimes\nabla u}{Q^2}\right)\nabla f\right),
\]
and in particular Luo and Sun record for $H$ the formula
\[
\Delta_g H=\frac{1}{v}\operatorname{div}\!\left(\left(vI-\frac{\nabla u\otimes\nabla u}{v}\right)\nabla H\right)
\]
[1110.3221], [2509.21018].

## 2. Willmore functional and the graphical Euler–Lagrange equation

For a two-dimensional closed surface $\Sigma$ with immersion $f:\Sigma\to\mathbb{R}^3$, the Willmore functional is written in Luo–Sun as
\[
W(f)=\frac14\int_\Sigma H^2\,d\mu_g,
\]
while other cited graph papers write the graph energy as
\[
W[u]=\frac14\int_\Omega H(u)^2\,Q\,dx
\]
or
\[
W_0(u)=\int_\Omega H^2 v\,dx,
\]
reflecting different normalization conventions [1110.3221], [2603.27848], [1503.01275]. The convention $H=k_1+k_2$ explains why the factor $\frac14$ appears in some formulations and why the graph formula for $H$ carries no additional $\frac12$ [1110.3221].

The Willmore equation in $\mathbb{R}^3$ is
\[
\Delta_g H+\frac12 H^3-2HK=0,
\]
equivalently
\[
\Delta_g H+2H\left(\frac14 H^2-K\right)=0,
\]
and in the convention of [1110.3221] it is also equivalent to
\[
\Delta_g H+H|A^o|^2=0
\]
[1110.3221], [2509.21018], [2603.27848]. For graphs, this is a quasilinear fourth-order PDE for $u$. A standard divergence formulation, derived in Deckelnick–Dziuk and used in several later works, is
\[
\operatorname{div}\!\left\{
\frac{1}{Q}
\Big(I-\frac{Du\otimes Du}{Q^2}\Big)\nabla(QH)
-\frac{H^2}{2Q}\,Du
\right\}=0,
\]
with
\[
H=\operatorname{div}\!\left(\frac{Du}{Q}\right),\qquad K=\frac{\det D^2u}{Q^4}
\]
[1110.3221], [2509.21018], [2603.27848], [1410.5547].

This divergence identity is central because it exposes the PDE as a conservation law for a vector field built from $u$. In the low-regularity boundary-value theory of [2509.21018], isolating the biharmonic term yields
\[
\Delta^2u=D_i\,b_1^{\,i}[u]+D_{ij}^2\,b_2^{\,ij}[u],
\]
where
\[
b_1[u]=D^2u\star D^2u\star \sum_{k=1}^3 Q^{-2k-1}P_{2k-1}(\nabla u),
\]
\[
b_2[u]=D^2u\star \sum_{k=1}^2 Q^{-2k-1}P_{2k}(\nabla u)
+ D^2u\star P_2(\nabla u)\star \big(Q(1+Q)\big)^{-1},
\]
and the pointwise bounds
\[
|b_1[u]|\le C\,|\nabla u|\,|D^2u|^2,\qquad
|b_2[u]|\le C\,|\nabla u|^2\,|D^2u|
\]
hold [2509.21018]. The same paper emphasizes that all monomials are at least cubic and that the highest derivatives enter at most quadratically in $b_1$ and linearly in $b_2$.

A complementary nondivergence expansion is used in the flow paper [2603.27848]:
\[
\Delta_{\Gamma(u)} H + 2 H \left(\frac14 H^2 - K\right)
= L(\nabla u)D^4u + R(\nabla u, D^2u, D^3u),
\]
where the principal fourth-order part is
\[
L(\nabla u)D^4u = \sum_{k+\ell=4} L_{k\ell}(\nabla u)\,\partial_{x_1}^k \partial_{x_2}^{\ell} u.
\]
Its coefficients depend only on $Du$ via $Q$, and satisfy the ellipticity bound
\[
\frac{1}{(1+\|\nabla u\|_{C^0}^2)^2}|\xi|^4
\le \sum_{k+\ell=4} L_{k\ell}(\nabla u)\,\xi_1^k\xi_2^\ell
\le 4 |\xi|^4.
\]
Thus $L(\nabla u)$ is uniformly elliptic while $|Du|$ is bounded [2603.27848].

## 3. Entire graphs and Bernstein-type rigidity

A central rigidity result for the Willmore equation on graphs is the Bernstein-type theorem of Luo and Sun. They consider a smooth entire graph
\[
\Sigma=\{(x,y,u(x,y)):(x,y)\in\mathbb{R}^2\}
\]
satisfying the Willmore equation
\[
\Delta_g H+\frac12 H^3-2HK=0
\]
and the integrability assumption
\[
\int_\Sigma H^2\,d\mu_g<\infty.
\]
Their theorem states: “Every smooth, entire graphical solution of (2.1) with finite $L^2$ norm of the mean curvature is a plane” [1110.3221].

The argument reduces the result to the earlier Chen–Lamm theorem that “Every smooth, entire graphical solution of (2.1) with finite $L^2$ norm of the second fundamental form is a plane” [1110.3221]. Luo and Sun prove that for any smooth entire graph in $\mathbb{R}^3$ with square integrable mean curvature,
\[
\int_\Sigma K\,d\mu_g=0.
\]
Combined with
\[
|A|^2=H^2-2K,
\]
this yields
\[
\int_\Sigma |A|^2\,d\mu_g=\int_\Sigma H^2\,d\mu_g,
\]
so $A\in L^2(\Sigma)$ if and only if $H\in L^2(\Sigma)$, and Chen–Lamm then applies [1110.3221].

The proof proceeds through two geometric ingredients. First, using the calibration argument of Colding–Minicozzi as cited in [1110.3221], one obtains quadratic area growth:
\[
|\Sigma\cap B(R)|\le C_2R^2.
\]
Second, because the Gauss map of a global graph takes values in the contractible upper hemisphere $S^2_+$, the area form $\beta$ on $S^2$ satisfies $\beta=d\alpha$ on $S^2_+$, and
\[
K\,d\mu_g=n^*\beta=d(n^*\alpha).
\]
With a compactly supported cutoff $\eta$, the estimate
\[
\left|\int_\Sigma \eta^2 K\sqrt{g}\,dx\wedge dy\right|
\le 4C_\alpha
\left(\int_\Sigma \eta^2|A|^2\,d\mu_g\right)^{1/2}
\left(\int_\Sigma |\nabla\eta|^2\,d\mu_g\right)^{1/2}
\]
is derived, and a logarithmic cutoff $\eta_\sigma$ satisfying
\[
\int_\Sigma |\nabla\eta_\sigma|^2\,d\mu_g\le \frac{C_5}{\log\sigma}
\]
implies
\[
\left|\int_\Sigma K\,d\mu_g\right|
=\lim_{\sigma\to\infty}\left|\int_\Sigma \eta_\sigma^2K\,d\mu_g\right|=0
\]
[1110.3221].

This theorem is explicitly described in [1110.3221] as a Bernstein-type theorem for Willmore graphs, replacing the minimal surface condition by the Willmore Euler–Lagrange equation and assuming $H\in L^2(\Sigma)$. Minimal graphs are included as the case $H\equiv 0$.

## 4. Radial symmetry, ODE reduction, and inverted catenoids

For radially symmetric graphs $u=u(r)$, $r=\sqrt{x^2+y^2}$, Chen and Li derive a sharp reduction of the graphical Willmore equation to an ODE [1410.5547]. Writing
\[
w(r)=u'(r),\qquad v(r)=\sqrt{1+w(r)^2},
\]
the mean curvature becomes
\[
H(r)=\left(\frac{w}{v}\right)'+\frac{1}{r}\frac{w}{v}
=\frac{w'}{v^3}+\frac{w}{rv},
\]
and smoothness at the origin implies
\[
H(0)=2w'(0)
\]
[1410.5547].

The divergence-form graph equation becomes a radial divergence condition. If the associated radial scalar is denoted by $f(r)$, then
\[
(rf(r))'=0.
\]
For smooth graphs on a full disk this forces $f\equiv 0$, whereas on a punctured disk one obtains
\[
f(r)=\frac{\lambda}{r}
\]
for a constant $\lambda$ [1410.5547]. In the smooth case, Chen and Li obtain the radial Willmore ODE
\[
w''+\frac1r w'
=
\frac{5w}{2(1+w^2)}(w')^2
+\frac{w(1+w^2+\tfrac12 w^4)}{r^2},
\]
and in the punctured-disk case the inhomogeneous version
\[
w''+\frac1r w'
=
\frac{5w}{2(1+w^2)}(w')^2
+\frac{w(1+w^2+\tfrac12 w^4)}{r^2}
+\frac{\lambda}{r}
\]
[1410.5547].

Their classification theorem states that if $u(r)$ is a smooth solution to the graphic Willmore equation on a disk centered at the origin, with $u(0)=c$ and mean curvature at $(0,0,c)$ equal to $2a$, then either $a=0$ and $u$ is constant, or $a\neq 0$ and the graph is a spherical cap contained in the half-sphere of radius
\[
R=\frac{1}{|a|},
\]
given explicitly by
\[
u(r)=c-\operatorname{sign}(a)\sqrt{R^2-r^2},\qquad 0\le r<R.
\]
In particular, smooth radially symmetric entire Willmore graphs in $\mathbb{R}^3$ must be flat [1410.5547].

For smooth radial solutions on a punctured disk with $u\in C^1(D(\rho))$, there exist a constant $\lambda$ and a function $\beta\in C^0([0,\rho))$ such that
\[
u''(r)=\frac{\lambda}{2}\log r+\beta(r),
\]
and the graph of $u$ is contained in a translated graphical piece of an inverted catenoid uniquely determined by $\lambda$ and $\beta(0)$ [1410.5547]. The inverted catenoid enters because minimal surfaces are mapped to Willmore surfaces under inversion. Chen and Li parametrize the catenoid by
\[
F_c(t,\theta)=\big(|c|\cosh t\cos\theta,\;|c|\cosh t\sin\theta,\;ct\big),
\]
apply inversion
\[
I_a(x,y,z)=\frac{(x,y,z-a)}{|(x,y,z-a)|^2},
\]
and identify the punctured-disk radial Willmore graph with the outermost graphical piece $E_{c,a}$ of the inverted catenoid, with
\[
c=-\frac{\lambda}{2}
\]
and $a$ determined by matching the asymptotics [1410.5547].

The same paper shows that if
\[
\int_{D(\rho)} H^2\,dx\,dy<\infty,
\]
then a radial solution on the punctured disk extends across the puncture as a $C^{1,\alpha}$ function for any $\alpha\in(0,1)$ and fits the inverted-catenoid model [1410.5547].

## 5. Boundary value problems and low-regularity existence theory

For graphs over bounded domains, the Willmore equation is naturally coupled to clamped boundary data. In the modern boundary-value literature, “clamped” means
\[
u|_{\partial\Omega}=g_0,\qquad \partial_\nu u|_{\partial\Omega}=g_1,
\]
so the boundary position and tangent half-planes are fixed along $\partial\Omega$ [2603.27848], [2509.21018]. In the notation of [2509.21018], the quantity
\[
f_1:=\nu g_1+\nabla_{\mathrm{tan}}g_0
\]
represents the trace of $\nabla u$ on $\partial\Omega$ and enters the estimates.

A major recent development is the low-regularity theory of graphical Willmore boundary problems. In [2509.21018], the Willmore equation is rewritten in divergence form so that it can be treated in weighted second-order Sobolev spaces. The resulting framework weakens the regularity assumptions on both the boundary and the Dirichlet data to the $C^{1+\alpha}$-class while the solution remains smooth in the interior, and extends existence theory to domains with merely Lipschitz boundaries within a purely weighted Sobolev framework [2509.21018].

The same paper treats two regimes. For $\partial\Omega\in C^{1+\alpha}$, $g_0\in C^{1+\alpha}(\partial\Omega)$, and $g_1\in C^\alpha(\partial\Omega)$, smallness is imposed on the slope:
\[
\|g_0'\|_{C^0(\partial\Omega)}+\|g_1\|_{C^0(\partial\Omega)}<\delta,
\]
with an a priori bound
\[
\|g_0'\|_{C^\alpha}+\|g_1\|_{C^\alpha}\le K.
\]
Then there exists a solution
\[
u\in C^{1+\beta}(\overline\Omega)\cap C^\infty(\Omega)
\]
of the clamped problem for any $\beta\in(0,\alpha)$ [2509.21018]. For Lipschitz $\partial\Omega$, the analysis is carried out in weighted Sobolev spaces
\[
W^{2,a}_p(\Omega),\qquad p>2,\quad 0<a<1-\frac2p,
\]
with clamped data in the trace space
\[
\dot W^{1+s}_p(\partial\Omega),\qquad s=1-a-\frac1p\in(0,1),
\]
and smallness required on
\[
\|\nu g_1+\nabla_{\mathrm{tan}}g_0\|_{L^\infty(\partial\Omega)}<\delta.
\]
Under either a BMO-modulo-VMO smallness assumption on the outward normal field or a restricted $(p,a)$ range depending on the Lipschitz constant, one obtains a variational solution
\[
u\in W^{2,a}_p(\Omega)\cap C^\infty(\Omega)
\]
[2509.21018].

A distinct but related advance is the low-regularity flow theory for the Willmore flow of graphs with clamped boundary data [2603.27848]. For a graph over a bounded domain, the $L^2$-gradient flow of the Willmore energy is
\[
\partial_t u
=
-\,Q\left\{\Delta_{\Gamma(u)} H+2H\left(\tfrac14H^2-K\right)\right\}
=
- Q\, \operatorname{div}\!\left\{
\frac{1}{Q}
\Big(I-\frac{Du\otimes Du}{Q^2}\Big)\nabla(QH)
-\frac{H^2}{2Q}\, Du
\right\},
\]
with
\[
\langle \partial_t f,N\rangle
=
-\left(\Delta_g H + 2 H \left(\tfrac14 H^2 - K\right)\right)
\]
in immersion form [2603.27848]. The paper develops time-weighted parabolic Hölder spaces that allow derivatives above parabolic order $s$ to blow up as $t\to0$ at controlled rates, thereby avoiding the classical fourth-order compatibility condition at $t=0$ [2603.27848].

Its short-time existence theorem states that if $\partial\Omega\in C^{4+\alpha}$, $g_0\in C^{4+\alpha}(\partial\Omega)$, and $g_1\in C^{3+\alpha}(\partial\Omega)$, then:

- if $u_0\in C^{1+\alpha}(\overline\Omega)$, there is a unique solution
  \[
  u\in C^{4+\alpha,1+\alpha/4}_{1+\alpha}(Q_T);
  \]
- if $u_0\in C^{0,1}(\overline\Omega)$ with $\|\nabla u_0\|_{L^\infty}\le M$, there exists a solution
  \[
  u\in C^{4+\alpha,1+\alpha/4}_{1}(Q_T)
  \]
  for some $T>0$ [2603.27848].

In the small-data Lipschitz regime, the same work proves global existence, uniform gradient bounds, and exponential convergence to a stationary solution. Specifically, if
\[
\|\nabla u_0\|_{L^\infty(\Omega)}
+\|g_0\|_{C^{4+\alpha}(\partial\Omega)}
+\|g_1\|_{C^{3+\alpha}(\partial\Omega)}<C,
\]
then the flow admits a global solution and
\[
\|\nabla u(\cdot,t)\|_{L^\infty(\Omega)}\le \overline C
\quad\text{for all } t\ge 0.
\]
Moreover, there exists a unique stationary solution $u_\infty\in C^{4+\alpha}(\overline\Omega)$ satisfying the elliptic Willmore equation with the prescribed boundary data, and for every $\beta\in(0,\alpha)$,
\[
\|u(\cdot,t)-u_\infty\|_{C^{4+\beta}(\overline\Omega)}
\le \widetilde C\,e^{-\widetilde C (t-1)}
\quad\text{for all } t\ge 1
\]
[2603.27848].

## 6. Variational relaxation, boundary conditions, and generalized graph classes

The Willmore equation for graphs also appears as the Euler–Lagrange equation in variational problems with prescribed boundary data. In [1503.01275], one considers a bounded $C^2$ domain $\Omega\subset\mathbb{R}^2$ with outward unit normal $\nu$, and the graph functional
\[
W_y(u)=\int_\Omega H^2\,v\,dx-y\int_\Omega K\,v\,dx,
\]
where
\[
H=\operatorname{div}\!\Big(\frac{\nabla u}{v}\Big),\qquad
K=\frac{\det D^2u}{v^4},
\]
and the pure Willmore functional is
\[
W_0(u)=\int_\Omega H^2v\,dx
\]
[1503.01275]. For smooth critical graphs, the interior Euler–Lagrange equation is given there as
\[
\Delta_g H+2H(H^2-K)=0.
\]
This is the form adopted in that paper and differs from the normalization used in [1110.3221], a convention issue that those sources explicitly note.

Two boundary conditions are distinguished in [1503.01275]. The Dirichlet, or clamped, condition prescribes
\[
u|_{\partial\Omega}=\varphi,\qquad \nabla u\cdot \nu|_{\partial\Omega}=\nabla\varphi\cdot \nu.
\]
Under these constraints, the first variation yields the interior Willmore equation and no additional natural boundary condition. The total Gaussian curvature is then determined by the boundary data and topology through Gauss–Bonnet:
\[
\int_\Omega K\,v\,dx + 2\int_{\partial\Gamma(u)} \kappa_g\,ds = 2\pi\,\chi(\Omega),
\]
where the geodesic curvature $\kappa_g$ of the boundary curve on the graph is explicitly determined by $(\Omega,\varphi)$ [1503.01275].

The Navier, or hinged, condition prescribes only
\[
u|_{\partial\Omega}=\varphi,
\]
and the natural boundary condition is
\[
H|_{\partial\Omega}=2yK_N,
\]
with $K_N$ the normal curvature of the boundary curve in the graph surface [1503.01275].

A principal result of [1503.01275] is that for $H^2$-regular graphs, bounds for the Willmore energy imply area and diameter bounds. In particular, if $u\in H^2(\Omega)$ with $(u-\varphi)\in H_0^2(\Omega)$, then there exists $C$ depending on $\Omega$ and $\|\varphi\|_{W^{2,1}(\partial\Omega)}$ such that
\[
\sup_{x\in\Omega}|u(x)|+\int_\Omega v\,dx\le C\big(W_0(u)^2+1\big)
\]
[1503.01275]. The paper then studies the $L^1$-lower semicontinuous relaxation
\[
W(u):=\inf\left\{\liminf_{k\to\infty}W_0(u_k):u_k\in H^2(\Omega),\ (u_k-\varphi)\in H_0^2(\Omega),\ u_k\to u\text{ in }L^1(\Omega)\right\}.
\]
If a sequence has bounded Willmore energy, then after passing to a subsequence it converges in $L^1(\Omega)$ to some
\[
u\in BV(\Omega)\cap L^\infty(\Omega),
\]
and the absolutely continuous contribution
\[
W_a(u)=\int_\Omega H_a^2\,Q^a\,dx
\]
provides a lower bound:
\[
W_a(u)\le \liminf_{k\to\infty}W_0(u_k)
\]
[1503.01275].

The same paper proves the existence of a minimizer
\[
u\in BV(\Omega)\cap L^\infty(\Omega)
\]
for the relaxed energy and shows that finite relaxed Willmore energy implies the attainment of the Dirichlet boundary data in an appropriate trace sense on the non-vertical part of the boundary [1503.01275]. A plausible implication is that the graphical Willmore theory naturally extends beyond the classical $H^2$ graph class to limit configurations with vertical parts, although the paper formulates the structural statements carefully in BV language rather than as a generalized PDE theory.

## 7. Flow, anisotropy, and numerical formulations

The Willmore equation for graphs is also studied through its gradient flow and through anisotropic generalizations. In the anisotropic graph framework of [1111.3043], one considers a rectangle $\Omega=(0,L_1)\times(0,L_2)$ and a convex, positive, $1$-homogeneous anisotropy $\gamma$. The anisotropic mean curvature is defined by
\[
H_\gamma(u)=\operatorname{div}\!\big(\nabla_p\gamma(\nabla u,-1)\big),
\]
and the anisotropic Willmore functional is
\[
\mathcal W_\gamma[u]
=
\frac12\int_\Omega \big(H_\gamma(u)\big)^2\,Q(u)\,dx.
\]
Introducing
\[
w_\gamma(u)=Q(u)\,H_\gamma(u),
\]
the stationary Euler–Lagrange equation is written in vector-conservative form as
\[
\operatorname{div}\!\Big(\mathbb E_\gamma(\nabla u,-1)\,\nabla w_\gamma(u)\Big)
-\frac12\operatorname{div}\!\Big(\frac{w_\gamma(u)^2}{Q(u)^3}\,\nabla u\Big)=0
\quad\text{in }\Omega,
\]
where
\[
\mathbb E_\gamma(\nabla u,-1)=D_p^2\gamma(\nabla u,-1)
\]
[1111.3043].

For the isotropic density
\[
\gamma_{\mathrm{iso}}(p,-1)=\sqrt{1+|p|^2},
\]
this reduces to
\[
\operatorname{div}\!\left(
\frac{1}{Q}\left(I-\frac{\nabla u}{Q}\otimes\frac{\nabla u}{Q}\right)\nabla w_{\mathrm{iso}}
\right)
-\frac12\operatorname{div}\!\left(\frac{w_{\mathrm{iso}}^2}{Q^3}\nabla u\right)=0,
\]
with
\[
w_{\mathrm{iso}}=Q\operatorname{div}\!\left(\frac{\nabla u}{Q}\right)
\]
[1111.3043]. This is the same structural identity as the isotropic graphical Willmore equation recorded in [1110.3221], [2509.21018], and [2603.27848], written in the notation of the anisotropic paper.

The anisotropic Willmore flow is formulated as
\[
\partial_t u
=
-\,Q(u)\,\Bigg[
\operatorname{div}\!\Big(\mathbb E_\gamma(\nabla u,-1)\,\nabla w_\gamma(u)\Big)
-\frac12\operatorname{div}\!\Big(\frac{w_\gamma(u)^2}{Q(u)^3}\,\nabla u\Big)
\Bigg],
\]
coupled to
\[
w_\gamma=Q\,\operatorname{div}\!\big(\nabla_p\gamma(\nabla u,-1)\big)
\]
[1111.3043]. For zero Dirichlet boundary data, the paper proves the continuous energy equality
\[
\int_\Omega \frac{(u_t)^2}{Q}\,dx
+\frac12\frac{d}{dt}\int_\Omega H_\gamma(u)^2\,Q(u)\,dx
=0,
\]
so the energy decays monotonically along the flow [1111.3043]. The same work defines weak solutions, derives a complementary finite volume discretization, reformulates it in finite-difference form, and proves the discrete energy equality
\[
\big((u_t^h)^2,1/Q^h\big)_h + \frac{d}{dt}\big((H_\gamma^h)^2,Q^h\big)_h = 0
\]
for zero Dirichlet data [1111.3043].

The bounded-domain isotropic flow theory of [2603.27848] supplies the analytic counterpart to this numerical and anisotropic perspective. Along the flow with time-independent Dirichlet data, it proves
\[
\frac{d}{dt}W(u(t))
=
-\frac12\int_\Omega |\partial_t u|^2 Q\,dx\le 0,
\]
hence
\[
\int_0^\infty\int_\Omega |\partial_t u|^2Q\,dx\,dt\le 2W(u(0))
\]
[2603.27848]. That identity is the precise gradient-flow energy law underlying the global convergence result.

Taken together, these works show that the Willmore equation for graphs admits several mutually reinforcing formulations: intrinsic geometric form via $\Delta_gH$ and $K$, divergence form suited to weak and low-regularity analysis, biharmonic-relative forms suited to elliptic estimates and fixed-point arguments, ODE reductions under radial symmetry, relaxed variational formulations in $BV$, and conservative parabolic forms suited to flow theory and structure-preserving numerics [1110.3221], [1410.5547], [1503.01275], [1111.3043], [2509.21018], [2603.27848].

Source: https://www.emergentmind.com/topics/willmore-equation-for-graphs