Cutoff Level-Set Mean Curvature G-Equation

• The cutoff level-set mean curvature G-equation is a level-set model for turbulent flame propagation that uses a positive-part operator to prevent unphysical negative laminar speeds.
• It incorporates curvature correction via the Markstein number and utilizes homogenization techniques in cellular flows to derive an effective burning velocity.
• The model poses analytical challenges through non-coercivity and discontinuities at DG = 0, addressed via viscosity solutions and deterministic game representations.

Searching arXiv for recent and foundational papers on the cutoff level-set mean curvature G-equation and related homogenization/level-set curvature analysis. The cutoff level-set mean curvature G-equation is a level-set model for flame front propagation in turbulent combustion in which the front is represented as the zero level set of a scalar function $G(x,t)$, with burnt region $\{G<0\}$ and unburnt region $\{G>0\}$. In its curvature-corrected form, the equation combines advection by an ambient incompressible flow, normal propagation at laminar flame speed, and a mean-curvature correction modulated by the Markstein number $d>0$. The cutoff is imposed through a positive-part operator to prevent a non-physical negative laminar flame speed, yielding the model

$$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$

For two-dimensional periodic incompressible flows, especially cellular flows, the principal analytical question is whether large-time front motion admits a direction-dependent effective burning velocity $\overline H(p)$, equivalently an effective Hamiltonian, so that solutions with planar initial data asymptotically propagate like planar fronts with uniform speed (Gao et al., 2022). The equation also sits within the broader viscosity-solution and level-set framework for motion by mean curvature, where the geometric operator $|DG|\,\mathrm{div}(DG/|DG|)$ is the standard mean-curvature flow operator written in level-set form (Colding et al., 2016).

1. Model, geometric meaning, and cutoff mechanism

The inviscid G-equation is the Hamilton–Jacobi model

$G_t+V(x)\cdot DG+s_L|DG|=0,$

where the front normal speed is the laminar speed plus the normal component of the ambient flow. Incorporating Markstein curvature effect replaces the constant laminar speed by Markstein’s relation

$s_\ell=s_L^0(1-d\kappa)_+,$

with curvature

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$

and, after normalizing $\{G<0\}$0, produces the cutoff mean-curvature G-equation (Gao et al., 2022).

The positive-part operator $\{G<0\}$1 is essential both physically and analytically. Physically, it ensures no “unburning,” because the laminar speed cannot become negative when curvature is large and positive. Analytically, it makes the operator more nonlinear and contributes to discontinuity at $\{G<0\}$2, which is one source of difficulty in homogenization and comparison arguments (Gao et al., 2022). Without cutoff, the term $\{G<0\}$3 can become negative when $\{G<0\}$4 is large and positive, which is regarded as unphysical in the combustion interpretation (Gao et al., 2022, Mitake et al., 2023).

The curvature contribution is the same geometric operator that appears in the level-set formulation of motion by mean curvature. If a level-set function $\{G<0\}$5 represents a moving hypersurface, then

$\{G<0\}$6

and the operator can be written as

$\{G<0\}$7

which makes its degenerate elliptic structure explicit (Colding et al., 2016). In the G-equation setting, the curvature term tends to smooth the front and reduce reaction surface area, generally slowing propagation compared to the inviscid model (Gao et al., 2022).

2. Cellular-flow setting and effective burning velocity

A central setting is the two-dimensional incompressible periodic cellular flow with Hamiltonian

$\{G<0\}$8

and amplitude $\{G<0\}$9, giving the divergence-free velocity field

$\{G>0\}$0

Its streamlines are the level curves of $\{G>0\}$1, and the cell $\{G>0\}$2 contains a single cell (Gao et al., 2022).

For any unit vector $\{G>0\}$3, the main homogenization-type result states that there exists a positive number $\{G>0\}$4 such that the solution of

$\{G>0\}$5

satisfies

$\{G>0\}$6

for a constant $\{G>0\}$7 depending only on the Markstein number $\{G>0\}$8 and the cellular flow amplitude $\{G>0\}$9 (Gao et al., 2022). In the rescaled thin-reaction-zone regime, with $d>0$0 and $d>0$1, the corresponding bound becomes

$d>0$2

The quantity $d>0$3 is the effective burning velocity, also described as the turbulent flame speed or ergodic constant (Gao et al., 2022).

The effective Hamiltonian admits a stationary cell-problem characterization. One seeks a periodic corrector $d>0$4 solving

$d>0$5

in the viscosity sense (Gao et al., 2022). The authors construct correctors by ergodic penalization: for $d>0$6, solve

$d>0$7

and show

$d>0$8

uniformly as $d>0$9 (Gao et al., 2022).

The resulting effective Hamiltonian is strictly positive for $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$0, positively homogeneous of degree one, and continuous in $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$1:

$$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$2

In the inviscid case $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$3, the cellular-flow benchmark due to Xin–Yu gives the sharp law

$$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$4

while with curvature the same $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$5 scaling is conjectured (Gao et al., 2022).

3. Viscosity solutions, homogenization difficulties, and correctors

The cutoff level-set mean curvature G-equation is treated in the viscosity-solution framework. This is consistent with the general theory for level-set mean-curvature equations, where the operator is degenerate parabolic, undefined at $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$6 in classical terms, and nevertheless admits continuous viscosity solutions with uniqueness via comparison principles (Colding et al., 2016). In the cutoff G-equation, the relevant difficulties are sharper because the Hamiltonian is non-coercive in the gradient and the mean-curvature operator is discontinuous at $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$7 (Gao et al., 2022).

The non-coercivity comes from the linear advection term $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$8, while degeneracy arises because the mean-curvature term vanishes at $$G_t+\left(1-d\,\mathrm{div}\!\left(\frac{DG}{|DG|}\right)\right)_+|DG|+V(x)\cdot DG=0.$$9 and is discontinuous as an operator there (Gao et al., 2022). Consequently, the standard perturbed test function method is not directly applicable. The paper therefore establishes existence of ergodic constants and sub- and supersolution correctors through PDE comparison combined with game-theoretic dynamics, rather than through the standard convex/coercive Hamilton–Jacobi homogenization route (Gao et al., 2022).

A key step is a negative linear drift estimate of the form

$\overline H(p)$0

which is then inserted into the Alvarez–Bardi penalization scheme for the stationary $\overline H(p)$1-equation. For small $\overline H(p)$2, one obtains

$\overline H(p)$3

so that $\overline H(p)$4 becomes a viscosity supersolution of a stationary equation with strictly positive right-hand side (Gao et al., 2022). A minimum value principle can then be applied to control minima on suitable regions. Combined with reachability and flow invariance, this yields uniform oscillation bounds on $\overline H(p)$5, which in turn imply uniform convergence of $\overline H(p)$6 to $\overline H(p)$7 (Gao et al., 2022).

This analytical structure places the cutoff curvature G-equation within the modern PDE theory of front propagation, while also showing why it is materially harder than the inviscid G-equation. For incompressible periodic flows, homogenization of the inviscid G-equation is well established, but the cutoff curvature model requires additional structure, especially in exploiting the underlying flow geometry (Gao et al., 2022).

4. Deterministic game representation and streamline geometry

The proof in two-dimensional cellular flow relies on a deterministic game characterization in the sense of Kohn–Serfaty. For time step $\overline H(p)$8 and $\overline H(p)$9 steps, the game trajectory $|DG|\,\mathrm{div}(DG/|DG|)$0 evolves by

$|DG|\,\mathrm{div}(DG/|DG|)$1

where $|DG|\,\mathrm{div}(DG/|DG|)$2 is chosen first by Player I, $|DG|\,\mathrm{div}(DG/|DG|)$3 is then chosen by Player II, and $|DG|\,\mathrm{div}(DG/|DG|)$4 (Gao et al., 2022). The value function satisfies the dynamic programming principle

$|DG|\,\mathrm{div}(DG/|DG|)$5

and as $|DG|\,\mathrm{div}(DG/|DG|)$6 with $|DG|\,\mathrm{div}(DG/|DG|)$7, one has local uniform convergence to the viscosity solution of the cutoff curvature G-equation (Gao et al., 2022).

Allowing the choice $|DG|\,\mathrm{div}(DG/|DG|)$8 encodes the positive-part cutoff. This is the mechanism by which Player I can “not burn” rather than “unburn,” and it is an essential link between the game and the cutoff operator (Gao et al., 2022). The continuum consistency calculation is expressed through the semicontinuous envelopes $|DG|\,\mathrm{div}(DG/|DG|)$9 of the curvature Hamiltonian (Gao et al., 2022).

The geometry of cellular streamlines is then used to design effective reachability strategies. For $G_t+V(x)\cdot DG+s_L|DG|=0,$0, the sets

$G_t+V(x)\cdot DG+s_L|DG|=0,$1

are invariant under the reverse flow $G_t+V(x)\cdot DG+s_L|DG|=0,$2 (Gao et al., 2022). The level sets $G_t+V(x)\cdot DG+s_L|DG|=0,$3 are convex inside cells, boundary strips near cell edges can be driven into interior invariant sets, and transitions between adjacent cells can be performed in uniformly bounded time (Gao et al., 2022). Specifically, the paper proves that from any point $G_t+V(x)\cdot DG+s_L|DG|=0,$4, the curve $G_t+V(x)\cdot DG+s_L|DG|=0,$5 is reachable in time $G_t+V(x)\cdot DG+s_L|DG|=0,$6, and any point on $G_t+V(x)\cdot DG+s_L|DG|=0,$7 can be reached by flowing along the curve in time $G_t+V(x)\cdot DG+s_L|DG|=0,$8 (Gao et al., 2022). It also proves uniform-time passage from boundary strips $G_t+V(x)\cdot DG+s_L|DG|=0,$9 into $s_\ell=s_L^0(1-d\kappa)_+,$0, and uniform-time transitions between neighboring cells (Gao et al., 2022).

These reachability results are combined with a minimum value principle for stationary problems to produce the oscillation bounds needed for the ergodic limit (Gao et al., 2022). An auxiliary construction in the appendix uses an expanding or shrinking ellipse whose boundary speed is strictly less than the PDE’s curvature-plus-advection speed in a strip near a vertical segment with $s_\ell=s_L^0(1-d\kappa)_+,$1, ensuring the front crosses that segment within a fixed time (Gao et al., 2022).

5. Comparisons with related curvature and level-set theories

The cutoff curvature G-equation is closely related to, but distinct from, the pure level-set mean-curvature flow studied in geometric PDE. In the latter, mean-convex arrival-time solutions satisfy the elliptic equation

$s_\ell=s_L^0(1-d\kappa)_+,$2

and there are strong regularity results: in the mean-convex case, the arrival time is twice differentiable everywhere, smooth away from the critical set, and solves the PDE classically everywhere (Colding et al., 2016). Those results are specific to the curvature-only setting and rely on geometric properties of mean-curvature flow (Colding et al., 2016).

For the flame-front equation, the curvature term appears with advection and with cutoff, so the sharp regularity theory of mean-curvature flow does not directly transfer. A plausible implication is that the geometric insight from arrival-time regularity is structurally relevant, but the cutoff G-equation requires separate analysis because advection and the positive-part operator fundamentally alter the stationary and large-time problems (Colding et al., 2016, Gao et al., 2022).

There is also a numerical comparison point. In level-set computation, curvature is typically evaluated from

$s_\ell=s_L^0(1-d\kappa)_+,$3

but standard central differences fail near “kinks,” namely regions equidistant to two or more interfaces where $s_\ell=s_L^0(1-d\kappa)_+,$4 is discontinuous (Lervåg et al., 2014). The paper on curvature calculations shows that such kinks produce spurious oscillations, large curvature spikes of order $s_\ell=s_L^0(1-d\kappa)_+,$5, and incorrect curvature sign, and develops improved discretization strategies including directional differences for normals, Macklin–Lowengrub’s method, Lervåg’s method, and Salac–Lu’s method (Lervåg et al., 2014). This suggests that numerical treatment of the cutoff curvature G-equation is inseparable from robust curvature evaluation when multiple fronts or equal-distance ridges are present.

A later three-dimensional study proposes data-driven mean-curvature correction for level-set methods, using feedforward networks that ingest transformed level-set, gradient, and curvature data on a $s_\ell=s_L^0(1-d\kappa)_+,$6 stencil, together with Gaussian-curvature-based stencil classification (Larios-Cárdenas et al., 2022). That work is formulated for mean-curvature computation rather than homogenization theory, but it demonstrates that under-resolved curvature estimation can be improved substantially relative to baseline finite differences (Larios-Cárdenas et al., 2022). A plausible implication is that any computational study of the cutoff G-equation in three dimensions must account for the numerical fragility of curvature computation independently of the analytical nonexistence issues.

6. Extensions, bifurcation phenomena, and large-time dynamics

The two-dimensional existence theory extends beyond the basic cellular flow. The same approach is described as adaptable to broader two-dimensional periodic incompressible flows

$s_\ell=s_L^0(1-d\kappa)_+,$7

including non-convex cells and “cat’s-eye” dynamics, under reasonable geometric conditions such as nondegenerate critical points (Gao et al., 2022). The extension again depends on geometric reachability across convex or concave boundaries and between islands and unbounded regions (Gao et al., 2022).

In dimensions three and higher, however, the existence of effective burning velocity can fail. For shear flows of the form

$s_\ell=s_L^0(1-d\kappa)_+,$8

the paper on bifurcation proves that in dimension $s_\ell=s_L^0(1-d\kappa)_+,$9 there exists a finite threshold $\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$0 such that the set of intensities admitting effective burning velocity is precisely

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$1

Moreover, there is a saturation threshold $\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$2 such that

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$3

while no effective burning velocity exists for $\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$4 (Mitake et al., 2023). In physical three-dimensional shear flow, corresponding to $\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$5, the result sharpens to

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$6

so the saturation threshold and the breakdown threshold coincide (Mitake et al., 2023).

The cutoff and non-cutoff effective velocities are related by

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$7

whenever the cutoff effective constant exists (Mitake et al., 2023). The mechanism of nonhomogenization is attributed to cutoff degeneracy: the local laminar speed

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$8

can vanish at isolated points where curvature is large, creating localized stagnation sets in the cell problem and obstructing periodic stationary solutions (Mitake et al., 2023). In three-dimensional shear, the characterization of the bifurcation point is linked to regularity theory for two-dimensional minimal-surface type equations, specifically through Simon’s Harnack inequality and a removability analysis for isolated singularities (Mitake et al., 2023).

A different extension introduces a non-negative source term and studies large-time behavior for

$\kappa=\mathrm{div}\!\left(\frac{DG}{|DG|}\right),$9

on the torus, where

$\{G<0\}$00

Under the assumptions

$\{G<0\}$01

the unique viscosity solution converges uniformly as $\{G<0\}$02 to a viscosity solution of the stationary problem (Calderon, 19 Sep 2025). The set $\{G<0\}$03 serves simultaneously as a monotonicity set and a uniqueness set, and the proof uses half-relaxed limits plus a comparison principle on $\{G<0\}$04 (Calderon, 19 Sep 2025). In the radially symmetric case with $\{G<0\}$05, the reduced Hamiltonian

$\{G<0\}$06

is convex in $\{G<0\}$07, allowing an optimal-control representation formula for the solution and its large-time limit (Calderon, 19 Sep 2025). This source-term model is not the same homogenization problem as the cellular-flow setting, but it shows that monotonicity induced by a non-negative source can restore asymptotic convergence in situations where standard Hamilton–Jacobi large-time theory is unavailable (Calderon, 19 Sep 2025).

7. Parameter dependence, misconceptions, and current directions

Several qualitative dependencies are established or conjectured. The effective burning velocity is strictly positive for nonzero propagation direction, positively homogeneous, and continuous in direction $\{G<0\}$08 (Gao et al., 2022). In shear flows, it decreases strictly with the Markstein number $\{G<0\}$09, consistent with the physical smoothing effect of curvature (Gao et al., 2022). If the curvature term is replaced by Laplacian diffusion, then turbulent flame speeds are dramatically slowed down: in two-dimensional cellular flows, $\{G<0\}$10 remains uniformly bounded in $\{G<0\}$11 (Gao et al., 2022). Numerical evidence also indicates that the cutoff is active for moderate intensities, for example $\{G<0\}$12 and $\{G<0\}$13, affecting $\{G<0\}$14 relative to the non-cutoff model (Gao et al., 2022).

One common misconception is that the cutoff is merely a technical regularization. In the cited work it is instead a physical constraint that prevents negative laminar flame speed and thereby excludes “unburning” (Gao et al., 2022, Mitake et al., 2023). Another misconception is that curvature correction is interchangeable with viscous diffusion. The available results distinguish these effects sharply: curvature slows fronts while preserving geometric dependence on level sets, whereas replacing curvature by Laplacian diffusion changes the asymptotic speed law qualitatively, including boundedness in $\{G<0\}$15 for cellular flows (Gao et al., 2022).

A further misconception is that homogenization should persist in higher dimensions whenever it holds in two dimensions. The shear-flow bifurcation result shows otherwise: in dimensions three and higher, effective burning velocity may cease to exist once advection exceeds a threshold, even though the non-cutoff problem homogenizes for all intensities (Mitake et al., 2023). This identifies the cutoff itself, rather than curvature alone, as a mechanism that can stall effective propagation in certain geometries (Mitake et al., 2023).

Current directions follow directly from the cited results. In two-dimensional cellular flow, the precise asymptotic scaling of $\{G<0\}$16 with curvature remains open; the paper conjectures the same $\{G<0\}$17 growth as in the inviscid case and asks whether curvature significantly reduces the prefactor (Gao et al., 2022). In broader two-dimensional flows, the geometry of anisotropy and directional dependence is described as future work (Gao et al., 2022). In three dimensions, an open problem is whether more mixing-dominated flows, such as swirled or chaotic incompressible periodic flows, may still admit effective burning velocities for all intensities despite the shear-flow bifurcation (Mitake et al., 2023). Together, these questions define the present research frontier of the cutoff level-set mean curvature G-equation.