Negative Results in Scalar Conservation Laws

• The paper shows that classical L¹ and BV formulations collapse under non-aligned discontinuous and smooth heterogeneous fluxes, leading to singular measures and loss of well-posedness.
• It demonstrates that vanishing-viscosity limits can generate δ-shock type distributions and finite-time L∞ blow-up, undermining standard entropy solution frameworks.
• The research emphasizes the need for broadened solution concepts, such as measure-valued or renormalized solutions, to properly capture and restore uniqueness in these problems.

Negative results for scalar conservation laws characterize precise mechanisms and regimes where the conventional theories of existence, uniqueness, and well-posedness—formulated in terms of $L^1$ or $\mathrm{BV}$ weak solutions or entropy solutions—fail. Recent research demonstrates that ill-posedness can arise from flux discontinuities with non-aligned structure, as well as from smooth heterogeneous fluxes when no $L^{\infty}$ bounds are assumed. These results highlight that the standard frameworks for scalar conservation laws do not universally guarantee well-defined evolution and that pathologies such as finite-time blow-up, measure-valued singularities, and nonuniqueness may generically occur. This article synthesizes the foundational negative results established in (Zajmović, 15 Sep 2025) and (Ghoshal et al., 13 Jan 2026), emphasizing their precise mathematical statements, proofs, and implications.

1. Pathologies for Scalar Conservation Laws with Discontinuous or Heterogeneous Flux

Scalar conservation laws of the form

$$u_t + \nabla_x \cdot f(x,u) = 0,\quad u(x,0) = u_0(x)$$

exhibit dramatically different behaviors depending on the regularity and alignment of the flux function $f(x,u)$. For multidimensional problems with flux discontinuities across interfaces not aligned with the normal (non-aligned), vanishing viscosity limits may not yield $L^1$ or $\mathrm{BV}$ functions, but rather singular measures supported on lower-dimensional sets (Zajmović, 15 Sep 2025). Conversely, in the smooth spatially varying flux case, explicit constructions now show that even for $f(x,u)$ smooth and for bounded initial data, solutions may blow up in $L^{\infty}$ norm and uniqueness of entropy solutions may fail (Ghoshal et al., 13 Jan 2026).

2. Singular Measure Formation and the Breakdown of $L^1$ and $\mathrm{BV}$ Theories

The multidimensional setting with a non-aligned discontinuous flux,

$$f^k(x, u) = \begin{cases} f_L^k(\hat{x}_k, u), & x_1 < \varphi(\hat{x}_1) \ f_R^k(\hat{x}_k, u), & x_1 > \varphi(\hat{x}_1) \end{cases}$$

where $\varphi \in C^2(\mathbb{R}^{d-1})$ is strictly monotone in each variable and the normal jumps do not align, leads to the following negative result (Zajmović, 15 Sep 2025): the vanishing-viscosity approximation

$$u_t^\varepsilon + \nabla_x \cdot f_\varepsilon(x, u^\varepsilon) = \varepsilon \Delta_x u^\varepsilon,\quad u^\varepsilon(x,0)=0$$

generates, in the limit $\varepsilon\to0$, a purely singular measure $\mu$ supported along the hypersurface $\Sigma = \{x_1 = \varphi(\hat{x}_1)\}$, explicitly of the form

$$\mu = M(\hat{x})\,\delta_{x_1 - \varphi(\hat{x})}(dx)\otimes dt$$

with $M(\hat{x}) < 0$ determined by the net flux jump. The mass concentrates as a Dirac layer on the surface $\Sigma$. Therefore, no nonnegative $L^1$ or $\mathrm{BV}$ solution arises in the limit—a phenomenon that structurally prevents the application of entropy or adapted entropy inequalities in the classical sense.

The key point is that, due to the non-alignment, the signed deficit in the fluxes across $\Sigma$ cannot be regularized or eliminated by viscous smoothing. This establishes the necessity to expand the solution concept to Radon measures and $\delta$-shock-type distributions for such conservation laws.

3. Uniqueness Failure and Blow-up in Conservation Laws with Smooth Heterogeneity

For one-dimensional scalar conservation laws with flux $f(x,u)$ smooth in $x$, it was shown in (Ghoshal et al., 13 Jan 2026) that there exist fluxes (e.g., $f(x,u) = x u^2$) and bounded initial data for which solutions develop $L^\infty$ blow-up at finite time. In particular, for $u_0$ monotone with $u_0(0)<0$, the explicit solution along characteristics gives

$p(t) = \frac{1}{t + 1/p(0)}$

so that $p(t)\to -\infty$ as $t \nearrow -1/p(0)$. The solution remains continuous in time in $L^1$, but forms a singularity at a finite blow-up time, precluding any $L^\infty$ bound beyond $t_{\mathrm{blow}}$.

More fundamentally, the Kružkov entropy equalities alone fail to ensure uniqueness in this regime. An explicit example is constructed where infinitely many entropy solutions exist, all continuous in $L^1$ and attainable from stationary profiles via moving interfaces or rarefactions. The classical doubling-of-variables argument for $L^1$-contraction relies crucially on a uniform $L^\infty$ bound; in its absence, the method fails and entropy solutions are not uniquely determined by the initial data.

4. Non-Existence and Global Ill-Posedness with Smooth Flux and Bounded Data

Negative results extend to global existence. For the Cauchy problem

$u_t + (x u^2 + u^4)_x = 0,\quad u(x,0) \equiv -1$

the solution constructed by characteristics exists only up to $t=1$, after which both the solution and its spatial support blow up everywhere and continuation is impossible, even in the weak $L^1_\mathrm{loc}$ sense (Ghoshal et al., 13 Jan 2026). This ill-posedness persists for higher-dimensional analogues, showing that, outside the setting of genuine nonlinearity or multiplicative-type fluxes, scalar conservation laws may fail to admit global weak solutions despite smoothness and bounded initial data.

5. Impact on Solution Concepts and Numerical Methods

The emergence of singular measures or nonuniqueness breaks the foundational framework used for scalar conservation laws:

• No weak-entropy solution exists in $L^1$ or $\mathrm{BV}$ for non-aligned flux discontinuities. Limits of vanishing-viscosity approximate solutions are purely singular (Zajmović, 15 Sep 2025).
• Standard theories based on adapted, Kružkov, or entropy inequalities collapse when $L^\infty$ bounds are lost (Ghoshal et al., 13 Jan 2026).
• Numerically, this implies that $L^1$-convergent schemes will not capture the true measure-valued solution: mass localized in $\delta$-layers on interfaces is smeared, leading to qualitatively incorrect behavior. Specialized discretizations capable of detecting or representing singularities are required—conventional approaches will truncate or lose the singular mass entirely.
• Given these pathologies, solution concepts must be broadened beyond standard entropy solutions, for instance, to measure-valued, $\delta$-shock, or renormalized solutions, with the addition of new admissibility criteria to restore uniqueness when possible.

6. Modified Uniqueness via Structural and Interface Conditions

A remedy for uniqueness failure in the smooth-heterogeneous flux case is to supplement the entropy condition with an interface (Lax-type) condition. In the setting of piecewise-constant fluxes $F^\delta(x,u)$ approximating $F(x,u)$, the front tracking method constructs approximate solutions by solving Riemann problems and accounting for entropic shocks and rarefactions.

The interface condition demands one-sided bounds at points of flux degeneracy,

$g>0 \implies u(a+,t)<+\infty,\;u(b-,t)>-\infty$

and vice versa for $g<0$. This prevents characteristics from emanating from unbounded traces at interfaces, enabling the extension of an $L^1$ contraction estimate and thus restoring uniqueness for entropy solutions under these additional structural conditions (Ghoshal et al., 13 Jan 2026). This approach demonstrates that, while the entropy inequality alone is insufficient, uniqueness can be recovered by careful handling of interface behavior and the introduction of admissibility conditions reflecting the underlying characteristics.

7. Summary Table: Negative Results for Scalar Conservation Laws

Pathology	Example Structure	Consequence for Solutions
Nonaligned flux discontinuity (Zajmović, 15 Sep 2025)	Discontinuous flux in non-normal direction, $f(x,u)$	Vanishing-viscosity limit: singular measure (Dirac layer), no $L^1$ or $\mathrm{BV}$ solution
Smooth flux, $L^\infty$ blow-up (Ghoshal et al., 13 Jan 2026)	$f(x,u) = x u^2$, nonzero initial data	Finite-time blow-up in $L^\infty$; uniqueness failure for entropy solutions
Coercive smooth flux,no existence (Ghoshal et al., 13 Jan 2026)	$f(x,u)=x u^2 + u^4$, $u(x,0)\equiv -1$	No global weak solution; solution unbounded everywhere after finite time

These negative results demonstrate that, for both discontinuous non-aligned fluxes and certain smooth heterogeneous fluxes, the classical $L^1$/$L^\infty$ PDE theory of scalar conservation laws is fundamentally inadequate. The only possibility for a consistent solution theory in these settings is to admit singular or measure-valued solutions and to carefully delineate admissibility and uniqueness via auxiliary structural constraints.