---
title: Curvature-Difference Methods Overview
url: https://www.emergentmind.com/topics/curvature-difference-methods
type: topic
---

# Curvature-Difference Methods Overview

Curvature-difference methods denote a heterogeneous set of techniques in which curvature, or a curvature-derived quantity, is estimated, regularized, or explicitly differenced in order to isolate a target geometric signal from confounding effects. In the recent literature summarized here, the term includes coordinate-wise Hessian subtraction for local memorization in diffusion models [2605.26756], kink-avoiding curvature evaluation in level-set methods [1409.6555; 1405.0829], monotone and filtered finite-difference discretizations for affine-curvature flow [1610.08831], semidiscrete curvature quotients for hyperbolic curvature flow [2204.12878], adaptive monitor functions built from curvature and curvature variation for Willmore flow [2601.01433], and curvature estimation from volume fractions via height functions and paraboloid fitting [2304.08643]. A plausible unifying interpretation is that these methods do not treat raw curvature as intrinsically trustworthy; instead, they modify its computation so that the retained signal is less contaminated by overfitting, medial-axis kinks, transport error, or mesh distortion.

## 1. Conceptual scope

Across these works, the operative object is always a curvature-related observable, but the source of error differs sharply by domain. In diffusion models, the confounder is intrinsic data-driven rigidity, which can produce high curvature even without memorization. In interface tracking, the confounder is non-smoothness of the signed-distance field near points equidistant to multiple interfaces. In curvature-driven PDEs, the central issue is nonlinear instability or loss of monotonicity under naive central differencing. In adaptive geometric evolution, curvature and curvature variation are used as monitor variables for mesh redistribution rather than as a force alone. In VOF-based interface reconstruction, the main trade-off is between static accuracy, robustness to advected volume-fraction error, and computational cost.

| Domain | Curvature quantity | Main corrective mechanism |
|---|---|---|
| Diffusion models | Diagonal Hessian of log-density | Subtract unconditional or less-trained baseline |
| Level-set interfaces | $\kappa=\nabla\!\cdot(\nabla\phi/|\nabla\phi|)$ | Avoid stencils crossing kinks; reconstruct locally smooth fields |
| Geometric flows | Discrete curvature quotients or affine-curvature operators | Monotone wide stencils, semidiscrete ODEs, filtered schemes |
| Willmore flow | $\kappa$ and $\partial_s\kappa$ | Monitor-based redistribution or tangential velocity |
| VOF interfaces | Mean curvature from heights or paraboloid fits | Estimator selection by accuracy-cost regime |

This heterogeneity is important. A common misconception is that “curvature-difference” always means subtracting two curvature fields. That interpretation is exact for the diffusion-model construction, but the broader literature also uses the term for finite-difference curvature quotients, wide-stencil reconstructions, or local extraction procedures whose purpose is to prevent curvature evaluation from crossing a singular or non-smooth representation.

## 2. Coordinate-wise curvature differences in diffusion models

The most explicit curvature-difference framework appears in the study of local memorization in text-to-image diffusion models [2605.26756]. At noise level $t$, with state $x_t\in\mathbb R^d$, the conditional and unconditional score networks are
\[
s_\theta(x_t,c)\approx\nabla_{x_t}\log p_\theta(x_t\mid c),
\qquad
s_\theta(x_t)\approx\nabla_{x_t}\log p_\theta(x_t),
\]
and the Hessian of the log-density is
\[
H_\theta(x_t,c)=\nabla_{x_t}s_\theta(x_t,c).
\]
The coordinate-wise curvature is the diagonal entry
\[
c_i(x_t,c)=\bigl[\nabla^2_{x_t}\log p_\theta(x_t\mid c)\bigr]_{ii}
=\bigl[-H_\theta(x_t,c)\bigr]_{ii},
\]
which measures the “sharpness” along coordinate $i$.

The geometric motivation is a curvature-variance correspondence. On a Gaussian, directions of low variance $\sigma^2$ have high curvature $\sim 1/\sigma^2$. The paper formalizes this through the curvature-variance identity
\[
\mathrm{Cov}[x_0\mid x_t]
= \frac{1}{\alpha_t^2}\bigl(\sigma_t^4\,\nabla^2_{x_t}\log p(x_t)+\sigma_t^2I\bigr),
\]
so that $\bigl(-\nabla^2_{x_t}\log p(x_t)\bigr)_{ii}$ is, up to scale, the inverse conditional variance. Local memorization is therefore characterized as a coordinate-wise variance collapse. The complication is that low conditional variance can also arise from intrinsic structure in $p_*(x)$, such as semantically rigid regions or simple textures. The method consequently defines Overfitting-Driven Memorization by subtracting a baseline curvature.

The conditional and baseline curvatures are
\[
h_i^{\text{cond}}(x_t)=\bigl[-H_\theta(x_t,c)\bigr]_{ii},
\qquad
h_i^{\text{base}}(x_t)=
\begin{cases}
\bigl[-H_\theta(x_t)\bigr]_{ii},&\text{unconditional baseline},\\[4pt]
\bigl[-H_{\tilde\theta}(x_t,c)\bigr]_{ii},&\text{less-trained baseline},
\end{cases}
\]
and the curvature-difference score is
\[
\Delta h_i(x_t)=h_i^{\text{cond}}(x_t)-h_i^{\text{base}}(x_t).
\]
The less-trained baseline is an earlier checkpoint, such as SD v1.1 for target SD v1.4. The unconditional baseline reuses the model’s unconditional U-Net used for classifier-free guidance.

A second construction gives a score-difference proxy with a Fisher-information interpretation. If
\[
\delta s(x_t)=s_\theta(x_t,c)-s_\theta(x_t),
\qquad
\Delta s_i(x_t)=\bigl[\delta s(x_t)\bigr]_i^2,
\]
then, at late sampling steps $t\to 0$,
\[
\Delta s_i(x_t)\approx \mathbb E\bigl[\Delta h_i(x_t)\bigr].
\]
For a less-trained baseline, the proxy becomes
\[
\Delta s_i(x_t)=\bigl[s_\theta(x_t,c)-s_{\tilde\theta}(x_t,c)\bigr]_i^2.
\]
This yields a geometric explanation for the widely used score-difference-based detection metric.

Algorithmically, the workflow is concise. A sample is generated with DDIM and classifier-free guidance. At the final step $t^*$, one computes either
\[
\delta s=s_\theta(x_{t^*},c)-s_\theta(x_{t^*},\emptyset)
\]
or
\[
\delta s=s_\theta(x_{t^*},c)-s_{\tilde\theta}(x_{t^*},c).
\]
The score-difference map is $\mathcal M_i=[\delta s]_i^2$. The curvature-difference estimate uses Hutchinson:
\[
\hat{\Delta h}
=
-\,\frac1K\sum_{k=1}^K
v_k\odot
\Bigl(\nabla_{x_{t^*}}\bigl(\delta s\cdot v_k\bigr)\Bigr),
\]
with $K$ random Rademacher vectors, followed by summation across channel indices for a spatial heatmap. The score-difference variant requires two forward/backwards per timestep and has negligible overhead. The curvature-difference variant requires $K$ Hessian-vector products, approximately $K$ additional backward passes. In practice $K=1$–16, and even $K=1$ suffices.

The empirical study uses Stable Diffusion v1.4 with baseline v1.1, v2.1 with baseline v2.0, and Realistic Vision v5.1, together with ground-truth template masks from Webster (2023), global-mem masks, and non-mem masks. Evaluation uses Intersection-over-Union and Pixel-Accuracy over $256\times256$ masks, with global detection by spatial averaging and AUC / TPR@1%FPR.

| Method on SD v1.4 | IoU | ACC |
|---|---:|---:|
| Bright Ending (template-verbatim only) | $\approx 0.75$ | $\approx 0.81$ |
| Raw curvature (template-verbatim only) | $\approx 0.59$ | $\approx 0.60$ |
| $\Delta h$ with unconditional baseline (template-verbatim only) | $\approx 0.90$ | $\approx 0.94$ |
| $\Delta s$ with unconditional baseline (template-verbatim only) | $\approx 0.83$ | $\approx 0.90$ |
| $\Delta h$ on all cases | $\approx 0.95$ | $\approx 0.97$ |
| Bright Ending on all cases | $\approx 0.56$ | not reported |

Aggregating $\mathcal M$ by spatial mean yields AUC $\approx 0.99$ and TPR@1%FPR $\approx 0.98$ for $\Delta h/\Delta s$, far above raw curvature or Bright Ending. Heatmaps from $\Delta h$ are sharply localized on memorized patches, whereas Bright Ending often lights up non-memorized regions. Later sampling steps, with $t$ close to $0$, dramatically improve the signal. Using SD v1.2 instead of v1.1 as baseline yields almost identical performance, and simple smoothing removes rare score-difference outliers.

The principal limitation is scope. The method is tailored to “verbatim” local memorization, specifically situations in which variability collapses coordinate-wise. It does not detect concept-level or style memorization that is globally distributed. It is also sensitive to late-timestep noise bounds and must be evaluated at sufficiently low noise, such as the final DDIM step.

## 3. Kink-robust curvature evaluation in level-set methods

In the level-set literature, curvature-difference methods arise from the failure of the textbook curvature formula when the signed-distance field contains kinks [1409.6555; 1405.0829]. With interface $\Gamma$ represented as the zero contour of $\phi(x,y,t)$, one writes
\[
n=\frac{\nabla\phi}{|\nabla\phi|},
\qquad
\kappa=\nabla\cdot n
=
\nabla\cdot\Bigl(\frac{\nabla\phi}{|\nabla\phi|}\Bigr).
\]
If $\phi$ is reinitialized as a signed-distance function, then it is only $C^0$ across points equidistant to two or more interfaces. In these regions, $|\nabla\phi|$ has a jump. Standard second-order central differences cross the kink, sample non-smooth data, and cease to converge. The reported pathology includes curvature spikes of order $1/\Delta x$, for example approximately $67$ on a mesh with $\Delta x\approx 0.015$, and even spurious sign changes.

Three improved schemes were summarized as avoiding direct differentiation of $\phi$ across a kink. The Macklin–Lowengrub method fits a local quadratic interface parametrization by least squares, computes interface curvature from the smooth fitted polynomial, and extrapolates it back to nearby grid nodes. Lervåg’s method reconstructs a smooth local signed-distance field $\phi_{\mathrm{loc}}$ from a monotone cubic Hermite spline and then applies the standard CD-2 stencils to $\phi_{\mathrm{loc}}$. The Salac–Lu method keeps one level-set field per object, computes curvature on each smooth field separately, and blends them by
\[
\kappa=\frac{\kappa_1|\phi_2|+\kappa_2|\phi_1|}{|\phi_1|+|\phi_2|}.
\]

| Method | Core mechanism | Reported outcome |
|---|---|---|
| MLM | Least-squares quadratic fit near interface | Curvature on contour itself; no differentiation across kink |
| LM | Spline-based local signed-distance reconstruction | Standard stencils applied to smooth $\phi_{\mathrm{loc}}$ |
| SLM | Separate level sets per object and blend curvatures | Avoids medial-axis kink entirely |
| LOLEX | Local extraction, body separation, reinitialization, standard stencil | Robust in 2D and 3D near topological change |

The benchmark behavior is consistent across test cases. For a static disk above a rectangle, CD-2 produces spikes up to $O(1/\Delta x)\approx \pm 67\,\mathrm{m}^{-1}$, while MLM, LM, and SLM all yield nearly uniform $\kappa\approx -4$ with pointwise errors $O(10^{-2})$ or smaller. In a two-drop collision in shear flow, curvature spikes from CD-2 create artificial high pressure in the thin film and prevent coalescence, whereas MLM, LM, and SLM remove the spikes and recover coalescence at the correct time. MLM predicts slightly earlier coalescence, while LM and SLM agree almost exactly.

The LOcaL Level-set Extraction method extends the same logic by isolating a small local subgrid around problematic stencil points. Kink regions are detected by the quality function
\[
Q(x)=|1-|\nabla\phi(x)||,
\]
and whenever a local stencil contains both signs of $\phi$ and $\max Q>\eta$, the algorithm copies $\phi$ onto a small local block, identifies connected bodies by a bodyscan, extracts a separate $\phi_{\mathrm{body}}$ for each, extrapolates ghost cells, reinitializes by solving
\[
\phi_\tau+\operatorname{sgn}(\phi)(|\nabla\phi|-1)=0,
\]
and finally recomputes $\nabla\phi$ and $\kappa$ on the smooth local field.

The quantitative validation is stronger than a mere artifact suppression claim. For a circle above a line with separations from $3.6\Delta x$ down to $0$, the standard method saturates at $1/\Delta x$ once the gap falls below approximately $2.4\Delta x$, whereas LOLEX stays within $O(1\%)$ of the analytical value $\kappa_{\mathrm{circle}}=10$ all the way to contact and converges with grid refinement. For concentric circles of width $1.6\Delta x$, LOLEX matches curve-fit normals to within max-RMS $0.0086$, while central and directional differences are off by $0.92$ and $1.78$ respectively. For a 3D sphere of radius $12.5\Delta x$ hovering $1.2\Delta x$ above a plane, the global $27$-point stencil mis-estimates curvature by $\pm 80\%$ near the kink, whereas LOLEX reduces this to $\pm 2\%$.

The same pattern appears in two-phase flow simulations. Standard curvature produces unphysical pressure oscillations in thinning films, delays rupture, and yields interleaved spurious pressure patches just before film rupture and at neck maximum. LOLEX remains smooth and is reported to incur only a $7\%$ CPU overhead versus the standard curvature routine on a $201^2$ test. The main caveat is that a residual $1$–$2\%$ curvature error remains from local reinitialization, and the method is not conservative, with mass loss remaining $O(10^{-3})$.

## 4. Finite-difference curvature quotients for geometric flows

In numerical geometric evolution, curvature-difference methods appear as stable discretizations of curvature-driven PDEs rather than as local interface diagnostics [1610.08831; 2204.12878]. Two distinct examples are prominent: affine-curvature flow in a level-set formulation and hyperbolic curvature flow for planar curves.

For affine-curvature flow, the governing PDE is
\[
u_t=F[u]=|\nabla u|\,\kappa_a,
\qquad
\kappa_a=k[u]^{1/3},
\]
with Euclidean curvature
\[
k[u]=\nabla\!\cdot\!\Bigl(\frac{\nabla u}{|\nabla u|}\Bigr)
=
\frac{u_{xx}u_y^2-2u_xu_yu_{xy}+u_{yy}u_x^2}{(u_x^2+u_y^2)^{3/2}},
\]
and equivalent form
\[
F[u]=\bigl(|\nabla u|^2\,\Delta_1u\bigr)^{1/3},
\qquad
\Delta_1u=\nabla u\,k[u].
\]
A naive central-difference discretization is observed to blow up near points where $u_{xx}\approx 0$. The remedy is a monotone, viscosity-compatible scheme. One-sided differences define monotone approximations of $|\nabla u|^\pm$, while $\Delta_1u$ is replaced by a wide-stencil median formula
\[
\Delta_1^e u(x,y)=2\,\frac{\mathrm{median}\{u_i-u(x,y)\}}{r^2},
\]
with neighbor values sampled on an approximately circular stencil of radius $r=h\,n_\theta$. Since the cube-root map is not monotone, it is decomposed into nondecreasing pieces and regularized through
\[
A^\delta(p,q)=\operatorname{sgn}(q)\min\{A(p,q),Kp,Lq\},
\]
leading to an elliptic operator $F^{e,\delta}$ and a filtered scheme
\[
F^f[u]=S_\epsilon\bigl(F^a[u],F^e[u]\bigr).
\]
Time stepping uses forward Euler with
\[
u^{n+1}=u^n+\Delta t\,F^{e,\delta}[u^n],
\qquad
\Delta t\le 1/C^h,
\]
where
\[
C^h=\frac{K}{h}+\frac{2L}{(h\,n_\theta)^2}.
\]
The Barles–Souganidis framework then gives convergence to the unique viscosity solution. The filtered scheme attains nearly the accuracy of central differencing in smooth regions while remaining convergent.

For hyperbolic curvature flow in the plane, the continuous normal-flow equation is
\[
x_{tt}+\beta\,x_t
=
\frac{1}{|x_\rho|}\,\partial_\rho\Bigl(\frac{x_\rho}{|x_\rho|}\Bigr)
-
(x_t\cdot\tau_t)\tau.
\]
A semidiscrete curve is represented by nodal values $x_j(t)$ on a periodic grid $\rho_j=jh$, with backward difference
\[
\delta x_j=\frac{x_j-x_{j-1}}{h},
\qquad
q_j=|\delta x_j|,
\qquad
\tau_j=\frac{\delta x_j}{q_j}.
\]
The discrete curvature vector is approximated by
\[
F_j^{(1)}:=\frac{2}{q_j+q_{j+1}}\frac{\tau_{j+1}-\tau_j}{h},
\]
and the scalar curvature can be taken as
\[
\kappa_j^h=\Bigl|\frac{2}{q_j+q_{j+1}}\frac{\tau_{j+1}-\tau_j}{h}\Bigr|.
\]
The resulting nodal ODE system is
\[
\ddot x_j+\beta\,\dot x_j
=
\frac{2}{q_j+q_{j+1}}\frac{\tau_{j+1}-\tau_j}{h}
-
(\dot x_j\cdot \dot\theta_j)\theta_j.
\]

The analysis proves a discrete energy law for $\beta=0$ and second-order convergence in $h$ under regularity assumptions and a uniform lower bound $|x_\rho|\ge c_0>0$:
\[
\max_{0\le t\le T}\Bigl\{\|x(t)-x^h(t)\|_{1,h}+\|\dot x(t)-\dot x^h(t)\|_{0,h}\Bigr\}\le C h^2.
\]
The fully discrete version uses $\Delta t\approx h$ with a semi-implicit update for the curvature term. Numerical experiments starting from a smooth strictly convex ellipse show bounded curvature initially, followed by finite-time blow-up of
\[
K_\infty=\max_j \frac{|\tau_{j+1}-\tau_j|}{q_j},
\]
and development of two corners. Outward normal initial velocity delays blow-up, and damping $\beta>0$ postpones the singularity further.

Taken together, these two lines of work show distinct operational meanings of a curvature-difference method. In affine flow, the emphasis is monotonicity, consistency, and explicit-Euler stability for a degenerate nonlinear PDE. In hyperbolic curvature flow, the emphasis is a direct curvature difference quotient on a moving mesh, coupled to energy control and second-order convergence.

## 5. Adaptive redistribution and tangential-velocity formulations for Willmore flow

For Willmore flow, curvature-difference methods are embedded in adaptivity rather than only in force evaluation [2601.01433]. Let $\Gamma(t)\subset\mathbb R^2$ be a smooth closed planar curve with parametrization $\mathbf X(\rho,t)$ and arc-length derivative $\partial_s=(1/|\partial_\rho \mathbf X|)\partial_\rho$. The scalar curvature and its variation are
\[
\kappa
=
-\,\partial_{ss}\mathbf X\cdot \mathbf n
=
-\,\frac{\partial_{\rho\rho}\mathbf X\cdot \mathbf n}{|\partial_\rho\mathbf X|^2},
\qquad
\partial_s\kappa=\frac{1}{|\partial_\rho\mathbf X|}\partial_\rho\kappa.
\]
The Willmore energy is
\[
W(\Gamma)=\tfrac12\int_\Gamma \kappa^2\,ds,
\]
and the $L^2$-gradient flow satisfies
\[
\partial_t\mathbf X\cdot \mathbf n=V,
\qquad
V=\partial_{ss}\kappa+\tfrac12\,\kappa^3.
\]

Time discretization uses BDF$_k$:
\[
D_t^{(k)}\mathbf X^{n+1}
=
\frac{1}{\Delta t}\sum_{p=0}^k \alpha_p\,\mathbf X^{\,n+1-p},
\]
with BDF1 coefficients $(1,-1)$, BDF2 coefficients $(\tfrac32,-2,\tfrac12)$, and BDF3 coefficients $(\tfrac{11}{6},-3,\tfrac32,-\tfrac13)$. The fully discrete finite-difference method couples the BDF approximation to centered spatial differences for $\mathbf X$, $\kappa$, and $V$. Nonlinearity is treated by Picard or Newton–Picard iteration. No explicit CFL-type constraint is reported beyond choosing $\Delta t$ small enough for convergence of the nonlinear solver.

The first adaptive strategy is weighted-arc-length equidistribution, denoted A-WAR. A positive scalar monitor $m(\rho,t)$ encodes $\kappa$ and $\kappa_s$, and the weighted arc-length density is
\[
m(\rho,t)\,|\partial_\rho\mathbf X|\,d\rho = dw(\rho,t),
\qquad
w(\rho,t)=\int_0^\rho m\,|\partial_q\mathbf X|\,dq.
\]
The equidistribution principle requires
\[
w_i^{\rm tar}=i\,L_w/M,\qquad i=0,\dots,M.
\]
Discrete weights are
\[
\omega_j=\tfrac12(m_j+m_{j+1})\,\Delta s_j,
\qquad
w_j=\sum_{k=0}^{j-1}\omega_k.
\]
Adaptive monitor selection uses
\[
C_0=\max|\kappa|,
\qquad
C_1=\max|\partial_s\kappa|,
\qquad
Q=\frac{\max\Delta s}{\min\Delta s},
\]
with monitor choices switching among
\[
m=1+\alpha|\kappa|,
\qquad
m=1+\alpha\bigl[(1-\beta)|\kappa|+\beta\kappa^2\bigr],
\qquad
m=1+\alpha\kappa^2,
\]
and, when curvature variation is large,
\[
m=1+\alpha|\kappa|+\gamma|\kappa_s|,
\qquad
m=1+\alpha\sqrt{\kappa^2+\gamma\,\kappa_s^2}.
\]

The second strategy, A-BDF$_k$-FDM, removes explicit redistribution by adding a tangential velocity enforcing
\[
M_f(\kappa,\kappa_s)\,\partial_\rho s=\text{constant in }\rho.
\]
Starting from the mesh functional
\[
\mathcal I[\rho]=\frac12\int_\Gamma \frac{1}{M_f}\,(\partial_s\rho)^2\,ds,
\]
the Euler–Lagrange equation gives
\[
\partial_\rho\bigl(M_f\,\partial_\rho s\bigr)=0.
\]
A relaxation dynamics for $\rho$ then yields
\[
\partial_t\mathbf X
=
V\,\mathbf n+V_\tau\,\boldsymbol\tau,
\]
with tangential velocity
\[
V_\tau
=
\frac{P}{\mathcal J\,(M_f|\partial_\rho\mathbf X|)^2}
\,\partial_\rho\Bigl[\frac{1}{M_f|\partial_\rho\mathbf X|}\Bigr]
=
\frac{P}{\mathcal J}\,(M_f|\partial_\rho\mathbf X|)^{-2}
\,\partial_\rho(M_f|\partial_\rho\mathbf X|).
\]

An additional Energy-Stable Correction algorithm introduces
\[
W_c^n=W^n+C,
\qquad
W^n=\tfrac12\sum_i (\kappa_i^n)^2\,|\delta_\rho \mathbf X_i^n|,
\]
and an auxiliary variable $R^n\approx W_c^n$. After a provisional update, one computes
\[
D^{n+1}=\sum_i (\overline V_i^{n+1})^2\,|\delta_\rho \mathbf X_i^n|,
\qquad
R^{n+1}=\frac{R^n}{1+\Delta t\,D^{n+1}/W_c^{n+1}},
\]
so that
\[
R^{n+1}-R^n
=
-\Delta t\,\frac{R^{n+1}}{W_c^{n+1}}\,D^{n+1}\le 0.
\]
This guarantees discrete energy decay at the theoretical level.

The reported numerical example uses the oscillatory initial curve
\[
x(\rho)=\bigl[1+0.3\sin(4\pi\rho)+0.2\cos(8\pi\rho)\bigr]\cos(2\pi\rho),
\]
\[
y(\rho)=\bigl[1+0.3\sin(4\pi\rho)+0.2\cos(8\pi\rho)\bigr]\sin(2\pi\rho).
\]
Standard BDF1-FDM leads to mesh ratio
\[
R_1=\max\Delta s/\min\Delta s
\]
growing to $O(10)$ by $t=10$, with clustering and eventual breakdown. By contrast, A-WAR and A-BDF1-FDM keep both $R_1$ and the weighted ratio
\[
R_2=\max(m_j\Delta s_j)/\min(m_j\Delta s_j)
\]
below $1.5$ for all $t\le 10$, while accurately capturing the evolving curve without mesh tangling.

## 6. Curvature estimation from volume fractions and comparative trade-offs

Curvature estimation from volume fractions provides a useful comparative frame because it isolates the accuracy-cost structure of several interface-curvature estimators on Cartesian meshes [2304.08643]. Four methods are compared: the height function (HF), PLIC-centroidal fitting (PC), PLIC-volumetric fitting (PV), and volumetric fitting (VF).

HF defines a pseudo-normal direction from the largest component of $\nabla\alpha$, forms column sums of the volume fraction $\alpha$, and computes derivatives of the resulting height field by second-order central differences. The mean curvature estimate is
\[
H = - \frac{\,h_{xx} + h_{yy} + h_{xx}h_y^2 + h_{yy}h_x^2 - 2h_{xy}h_xh_y\,}
{2\,(1 + h_x^2 + h_y^2)^{3/2}}.
\]
For exact $\alpha$, HF is $O(\Delta^2)$ in curvature and converges to machine precision in static random paraboloid tests, but it fails more often when $C\Delta>10^{-1}$ because complete columns are unavailable. If $\alpha$ comes from second-order VOF advection and is only $O(\Delta)$ in $L^\infty$, the curvature may degrade to $O(1)$ for small $\Delta$.

PC and PV fit a local paraboloid
\[
z' = a_0+a_1x'+a_2y'+a_3x'^2+a_4x'y'+a_5y'^2
\]
in a rotated frame aligned with the target-cell normal, then evaluate curvature from the fitted coefficients. PC minimizes a weighted least-squares error over neighboring PLIC centroids, using radial and area-projection weights. PV instead matches the actual PLIC planes in a least-squares sense of projected area and assembles a symmetric $6\times 6$ normal system from polygon moments computed by Green’s theorem. VF fits the same paraboloid directly to the volume fraction field by minimizing
\[
J(a)=\sum_{i\in N_c} w_i^R\,[\widetilde\alpha_i(a)-\alpha_i]^2,
\qquad
\widetilde\alpha_i = \frac{1}{\mathrm{Vol}_i}\int_{\Omega_i}\chi_f\,dV,
\]
using exact analytic polyhedron-paraboloid intersections and a Levenberg–Marquardt solve.

| Method | Static behavior | Dynamic / noisy behavior | Cost |
|---|---|---|---:|
| HF | $O(\Delta^2)$ to machine precision | Degrades for advected $\alpha$; column failures for high curvature | $1\times$ |
| PC | $O(\Delta^2)\to O(\Delta^1)$ at low $C\Delta$ | Highest spurious currents of the three dynamic methods | $\sim 10\times$ |
| PV | $O(\Delta^2)\to O(\Delta^1)$ only at very low $C\Delta$ | Best balance; lowest overall spurious current in 3D translating droplet | $\sim 10\times$ |
| VF | $O(\Delta^2)$ to machine precision | Dynamic tests not performed | $\sim 100\times$ |

The detailed results sharpen these trade-offs. In static random paraboloids, HF and VF both reach second-order convergence down to machine precision. PC and PV are second-order down to approximately $C\Delta\sim 10^{-2}$, after which they transition to $O(\Delta^1)$ behavior. Under perturbations of $\alpha$ with $k=0.1$, PV maintains lower error than HF and PC across all $C\Delta$, while VF errors are less than or approximately equal to PV errors. Increasing the PV stencil from $S=3$ to $S=5$ extends convergence down to $C\Delta\sim 10^{-4}$, with only minor further gain at $S=7$.

In dynamic tests, HF gives the lowest spurious velocities for a static droplet at high Laplace number, but in a 3D translating droplet it produces higher capillary numbers when $\alpha$ has large transport errors. PV gives the lowest $\mathrm{Ca}_{\mathrm{rms}}/\mathrm{Ca}_{\max}$ among the three tested dynamic methods and is recommended as the default curvature estimator when coupling to a two-phase Navier–Stokes solver. HF remains attractive as a first pass, with PV as a backup where height columns fail. VF is principally reserved for static or near-static configurations where extremely high-fidelity curvature is required and its two-orders-of-magnitude cost penalty is acceptable.

A broader implication, suggested by the combined literature, is that curvature-difference methods are best understood not as a single numerical recipe but as a design principle. They intervene precisely where naive curvature evaluation is confounded: by intrinsic data rigidity in diffusion models, by non-differentiability in signed-distance fields, by loss of monotonicity in nonlinear PDE discretizations, by mesh distortion in fourth-order geometric flows, or by transport noise in VOF reconstructions. The recurring result is not that curvature becomes universally easy to compute, but that reliable curvature information depends on explicitly controlling the mechanism that corrupts it.

Source: https://www.emergentmind.com/topics/curvature-difference-methods