---
title: Entropy Correction Artificial Viscosity
url: https://www.emergentmind.com/topics/entropy-correction-artificial-viscosity
type: topic
---

# Entropy Correction Artificial Viscosity

Searching arXiv for recent and foundational papers on entropy correction artificial viscosity.
Entropy correction artificial viscosity (ECAV) is an approach for enforcing a semi-discrete entropy inequality through an entropy dissipative correction term. In high-order discretizations of nonlinear conservation laws, it replaces or supplements more traditional entropy-stable constructions based on algebraic flux differencing, summation-by-parts (SBP) operators, and entropy conservative two-point fluxes. In the ECAV formulation developed for discontinuous Galerkin (DG) methods, the artificial viscosity coefficient is determined from the local violation of a cell entropy inequality and the local entropy dissipation, yielding a parameter-free, cell-local correction that is zero when the entropy inequality is already satisfied [2501.16529]. Subsequent work extended the construction to multiple artificial viscosities, to finite difference discretizations, and to analyses of the viscous DG discretization itself, including contact preservation and an $O(h)$ upper bound on the ECAV coefficient [2604.03158] [2508.21226] [2602.23210].

## 1. Conceptual basis and development

Entropy stable methods have become increasingly popular in the field of computational fluid dynamics because schemes which satisfy some form of a discrete entropy inequality typically behave much more robustly, and do so in a way that is hyperparameter free [2508.21226]. In DG, one established route to entropy stability uses entropy conservative fluxes in an algebraic flux differencing formulation together with SBP discretization matrices. However, explicit expressions for such two-point finite volume fluxes may not be available for all systems, or may be computationally expensive to compute [2501.16529]. ECAV was proposed as an alternative approach to constructing entropy stable DG methods using an entropy correction artificial viscosity, where the artificial viscosity coefficient is determined based on the local violation of a cell entropy inequality and the local entropy dissipation [2501.16529].

The immediate precursor of this line of work is the entropy correction of Abgrall, Offner, and Ranocha (2022), which was reinterpreted and modified into an artificial-viscosity form. The resulting method recovers the same global semi-discrete entropy inequality that is satisfied by entropy stable flux differencing DG methods, while avoiding dependence on bespoke entropy conservative flux functions [2501.16529]. This design choice situates ECAV between classical shock-capturing artificial viscosity methods, which typically use sensors and tunable parameters, and fully entropy-conservative or entropy-stable flux differencing schemes, which derive robustness from discrete thermodynamic structure.

A central organizing principle is that the entropy inequality is enforced by adding only as much dissipation as needed. This minimality is explicit in the coefficient formulas used in both DG and finite difference settings. The method is therefore not a generic viscosity model; it is an entropy correction whose magnitude is computed from an algebraic entropy deficit and an associated dissipative mechanism [2501.16529] [2508.21226].

## 2. Mathematical formulation in high-order DG

The target problem is a system of conservation laws
\[
\frac{\partial \boldsymbol{u}}{\partial t} + \sum_{i = 1}^d \frac{\partial \boldsymbol{f}_i(\boldsymbol{u})}{\partial x_i} = \boldsymbol{0},
\]
equipped with a convex entropy $S(u)$, entropy variables $v(u) = \frac{\partial S}{\partial u}$, and entropy fluxes $F_m(u)$. The entropy inequality reads
\[
\frac{\partial S(u)}{\partial t} + \sum_{m = 1}^d \frac{\partial F_m(u)}{\partial x_m} \leq 0
\]
or, over a DG cell $D$,
\[
\int_D \frac{\partial S(u)}{\partial t} + \int_{\partial D} \sum_{m=1}^d \left( v^T f_m(u) - \psi_m(u) \right) n_m \le 0,
\]
where $\psi_m(u)$ is the entropy potential [2604.03158].

In the monolithic DG construction, artificial viscosity is introduced through a modified PDE of the form
\[
\frac{\partial \boldsymbol{u}}{\partial t} + \sum_{i=1}^d \frac{\partial \bm{f}_i(\boldsymbol{u})}{\partial x_i}
=
\sum_{i,j=1}^d \frac{\partial}{\partial x_i}
\left(
\epsilon_k(\boldsymbol{u}) \bm{K}_{ij} \frac{\partial \bm{v}}{\partial x_j}
\right),
\]
with $\epsilon_k(\boldsymbol{u})$ a per-element, non-negative artificial viscosity coefficient and $\bm{K}$ positive semidefinite. A common choice is Laplacian viscosity with $\bm{K}_{ij} = \delta_{ij} \frac{\partial u}{\partial v}$ [2604.03158].

The DG discretization uses a BR1 treatment of the viscous term. Auxiliary variables $\boldsymbol{\Theta}_i$ approximate entropy-variable gradients, viscous fluxes $\boldsymbol{\sigma}_i$ are formed from $\epsilon_k(u_h)\bm{K}_{ij}\boldsymbol{\Theta}_j$, and the viscous contribution $\boldsymbol{g}_{\mathrm{visc}}$ is inserted into the weak form [2604.03158]. The entropy correction is driven by a local, cell-based entropy residual
\[
\delta_k(u_h) =
\sum_{m=1}^d \int_{D^k}
-
\left(\frac{\partial v_h}{\partial x_m}\right)^T f_m(u_h)
+
\int_{\partial D^k}
\psi_m(\tilde{u}_h) n_m.
\]
In the related DG construction with entropy projection, the residual is written
\[
\delta_k(\bm{u}_h) =
\sum_{m=1}^d
\left(
-\langle \bm{f}_m(\bm{u}_h), \partial_{x_m}\Pi_N \bm{v}(\bm{u}_h) \rangle_{D^k}
+
\langle \psi_m(\tilde{\bm{u}}) n_m, 1 \rangle_{\partial D^k}
\right),
\]
where $\Pi_N$ is the $L^2$ projection onto degree-$N$ polynomials and $\tilde{\bm{u}}$ is the entropy-projected conservative variable [2501.16529].

The minimal elementwise coefficient is then chosen as
\[
\epsilon_k(u_h) \geq
\frac{-\min(0,\delta_k(u_h))}
{\sum_{i,j=1}^d
(
\bm{K}_{ij}\bm{\Theta}_j,\bm{\Theta}_i
)_{D^k}},
\]
with a small regularization of the denominator for numerical stability [2604.03158]. In the DG entropy-projection formulation, the same construction appears as
\[
\epsilon_k(\bm{u}_h) =
\frac{-\min(0,\delta_k(\bm{u}_h))}
{\sum_{i,j=1}^d \langle \bm{K}_{ij}\bm{\Theta}_j,\bm{\Theta}_i\rangle_{D^k}},
\]
and if $\delta_k(\bm{u}_h)\ge 0$, no viscosity is applied [2501.16529].

The entropy dissipation guarantee is explicit: for the BR1 discretization,
\[
\sum_{k} -(\boldsymbol{g}_\text{visc}, v_h)_{D^k}
=
\sum_{k} \sum_{i, j=1}^d
(\epsilon_k(u_h) \bm{K}_{ij} \bm{\Theta}_j, \bm{\Theta}_i)_{D^k}
\geq 0.
\]
Thus, artificial viscosity applied in this way strictly enforces the discrete entropy inequality [2604.03158].

## 3. Viscous discretization, size of the coefficient, and contact preservation

A distinct line of analysis concerns the choice of viscous DG discretization for ECAV. When the artificial viscosity is discretized using a local discontinuous Galerkin (LDG) method, an $O(h)$ upper bound on the ECAV coefficient can be proved [2602.23210]. This upper bound indicates that ECAV does not result in a restrictive time-step condition, so the explicit time-step remains hyperbolic rather than parabolic in scale [2602.23210].

The LDG choice is motivated by structural properties of the discrete gradient operator. In the analysis summarized for [2602.23210], LDG yields a local gradient operator whose norm is bounded below by the true gradient norm, and the nullspace of the LDG gradient for polynomials is only the constant function. By contrast, BR1 can admit non-constant null modes, which can lead to numerical pathologies and spuriously large viscosity coefficients. A plausible implication is that the ECAV mechanism depends not only on the entropy residual, but also on whether the viscous discretization provides a robust denominator in the coefficient formula.

The same analysis shows that ECAV is contact preserving [2602.23210]. With piecewise constant solutions such as stationary contacts, the entropy residual vanishes, so ECAV adds no artificial viscosity. If the underlying DG scheme is contact preserving, ECAV with LDG does not alter that property [2602.23210]. This distinguishes ECAV from many traditional shock-capturing artificial viscosity methods, whose dissipation is selected through problem-dependent indicators and thresholds and can smear contacts.

The contrast with heuristic shock capturing is consequential. Traditional approaches add viscosity based on sensors or smoothness indicators and typically require parameter tuning. ECAV computes its coefficient directly from the algebraic balance of entropy at the element level and injects only as much dissipation as needed to restore the inequality [2602.23210]. In this sense, the artificial viscosity is neither prescribed a priori nor globally calibrated; it is an entropy deficit correction.

## 4. Multiple viscosities and targeted dissipation

The original DG formulation used a single monolithic viscosity coefficient. This was generalized to multiple dissipation models, including viscosity and thermal diffusivity, by writing
\[
\frac{\partial \boldsymbol{u}}{\partial t} + \sum_{i=1}^d \frac{\partial \boldsymbol{f}_i(\boldsymbol{u})}{\partial x_i}
=
\sum_{m=1}^M \boldsymbol{g}_\text{visc}^m,
\]
with each viscous term $\boldsymbol{g}_\text{visc}^m$ assigned its own elementwise coefficient $\epsilon_k^m(u_h)$ [2604.03158]. The coefficients are allocated by solving an elementwise constrained minimization problem,
\[
\min_{\boldsymbol{\epsilon}}
\frac{1}{2}\boldsymbol{\epsilon}^T \bm{W}\boldsymbol{\epsilon}
\quad
\text{subject to}
\quad
\boldsymbol{r}^T \boldsymbol{\epsilon} = \sigma_k(u_h),
\]
where $\sigma_k(u_h) = -\min(0,\delta_k(u_h))$ and $r_m$ denotes the entropy-dissipative contribution of model $m$ [2604.03158]. With equal weighting, the closed-form solution is
\[
\epsilon_i = \frac{\sigma_k(u_h) r_i}{\sum_{j=1}^M r_j^2}.
\]

This construction divides dissipation among models proportional to their relative contribution to entropy dissipation, for minimal overall added dissipation [2604.03158]. The purpose is not merely algebraic convenience. Different artificial viscosity mechanisms can be directed toward different numerical or physical artifacts: Laplacian viscosity for general regularization, thermal diffusivity for false heating or overheating, and spectral vanishing viscosity (SVV) for under-resolved turbulence or shear [2604.03158].

| Model | Role described in the literature | ECAV allocation principle |
|---|---|---|
| Laplacian viscosity | General conservative regularization | Included in the constrained minimization |
| Thermal viscosity | Acts on energy; reduces spurious temperature spikes | Activated through its entropy-dissipative contribution |
| SVV | Targets high modes and under-resolved turbulence/shear | Receives dissipation only where needed |

The reported examples illustrate the selectivity of the multiple-viscosity formulation. In a 1D receding flow, adding thermal artificial viscosity alongside Laplacian artificial viscosity reduces spurious temperature spikes while leaving density and pressure profiles largely unaffected [2604.03158]. In a 2D Kelvin–Helmholtz problem, the SVV coefficient is only activated where needed, in regions of high shear or turbulence, while in a 2D Riemann problem the extra SVV coefficient is essentially zero when the additional mechanism is physically unnecessary [2604.03158]. These results support the characterization of ECAV as minimally dissipative and physically targeted, rather than monolithic.

## 5. Finite-difference ECAV and the relation to knapsack limiting

ECAV was subsequently introduced for finite difference discretizations in an SBP flux-differencing framework [2508.21226]. In that setting, the high-order residual at node $i$ is written
\[
r_i = \sum_j \|n_{ij}\|
\left[
f_{ij}^H + \theta_{ij}(u_i-u_j)
\right],
\]
where $f_{ij}^H$ is the high-order, central, high-accuracy flux, $\theta_{ij}$ is a symmetric non-negative artificial viscosity coefficient, and $\|n_{ij}\|$ are geometric weights [2508.21226]. The cell entropy inequality becomes a nodal constraint,
\[
\sum_{j} \|n_{ij}\|
\left[
(v_j-v_i)^T f_{ij}(\theta_{ij})
-
(\psi(u_j)-\psi(u_i))\cdot \hat{n}_{ij}
\right]
\le 0,
\]
which reduces to a linear inequality $a_i^T\theta_i \ge b_i$ with
\[
a_{ij} = \|n_{ij}\| (v_j-v_i)^T (u_j-u_i)\ge 0.
\]
The minimal local correction solves
\[
\min_{\theta_i \ge 0}\ \|\theta_i\|_2^2
\quad
\text{subject to}
\quad
a_i^T\theta_i \ge b_i,
\]
and has the explicit solution
\[
\theta_i=
\begin{cases}
0, & b_i\le 0,\\[4pt]
\frac{b_i}{a_i^Ta_i}a_i, & b_i>0.
\end{cases}
\]
Conservation is recovered by symmetrization, $\theta_{ij} = \max\{\hat{\theta}_{ij},\hat{\theta}_{ji}\}$ [2508.21226].

This finite-difference version preserves high order accuracy in sufficiently smooth conditions, is entropy stable, and is hyperparameter free [2508.21226]. It also clarifies the relation between ECAV and knapsack limiting (KL). KL blends a high-order scheme with a low-order, positivity preserving, and entropy stable scheme, whereas ECAV adds the minimal local artificial viscosity needed to satisfy the cell entropy inequality [2508.21226]. KL can additionally incorporate positivity constraints for the compressible Euler and Navier–Stokes equations; ECAV, by itself, does not preserve positivity [2508.21226].

The distinction is structural. ECAV is a viscosity-based entropy correction; KL is a constrained flux blending strategy. Yet, for the special case where the low-order flux is local Lax–Friedrichs, the KL flux can be written as a special case of the ECAV form with $\theta_{ij}$ scaled by the maximum wavespeed [2508.21226]. This places both methods within a common framework of local convex optimization driven by entropy admissibility.

## 6. Related interpretations, neighboring methods, and recurring misconceptions

ECAV sits within a broader landscape of entropy-based artificial dissipation, but it is not identical to older entropy viscosity methods. In entropy viscosity, the local viscosity coefficient is typically
\[
\nu(x)=\min(\nu_{\rm max}(x),\nu_{EV}(x)),
\qquad
\nu_{EV}(x)=\alpha_{EV} h^\beta |r_{EV}(x)|,
\]
with $\nu_{\rm max}=\alpha_{\rm max} h \max|u|$ and a scaling exponent $\beta$ that strongly affects performance [1703.06917]. The analysis of shock problems in [1703.06917] shows that $\beta=1$ is preferable to $\beta=2$ for robust shock capturing in high-order methods, especially with moving shocks and multiple shocks of different strengths. ECAV differs in that its coefficient is not set by a residual scaling law with tunable parameters; it is determined by the amount of entropy dissipation required to close a discrete entropy imbalance [2501.16529] [2508.21226].

Another neighboring line concerns entropy correction through filters. In DG discretizations stabilized via SIAC filters, entropy correction is written in a filtered form such as
\[
\mathbf{c}_G = \alpha_{\text{glob}}(\mathbf{I}_G-\mathbf{K}_G)\mathbf{w}_G,
\]
with conservative filter correction and optional dissipative augmentation through an artificial viscosity estimate [2312.05232]. This is an entropy-correction strategy rather than the elementwise coefficient construction that defines ECAV. It nevertheless illustrates a shared theme: entropy correction can be appended to a high-order method while preserving conservation and controlling dissipation [2312.05232].

A further misconception concerns the sign of artificial viscosity. Artificial viscosity has traditionally been interpreted as a positive, spatially acting regularization, and local entropy stability is often linked to pointwise positivity. The data-driven study of linear convection in [2602.07733], however, reports locally negative artificial viscosity near extrema while the numerical solution remains stable and nearly exact. The paper resolves this by reinterpreting artificial viscosity as a space-time closure that compensates unresolved truncation errors while enforcing entropy stability through global dissipation balance rather than pointwise positivity. Within that framework, entropy stability constrains the integrated dissipation budget rather than the pointwise sign of spatial viscosity [2602.07733]. This does not redefine ECAV directly, but it broadens the conceptual context in which entropy-based viscosity is interpreted.

Finally, high-order artificial dissipation in SBP-compatible correction procedure via reconstruction (CPR) schemes provides another adjacent perspective. There, stable artificial dissipation operators are built in the SBP form
\[
-\left(\mathbf{M}^{-1}\mathbf{D}^T\mathbf{M}\mathbf{a}\mathbf{D}\right)^s\mathbf{u},
\]
rather than by naïve discretization, in order to guarantee conservation and discrete $L_2$ stability [1606.00995]. This reinforces a point that also appears in ECAV: the entropy or energy effect of artificial viscosity depends as much on the discretization of the viscous term as on the scalar coefficient itself.

Across these strands, the unifying theme is that entropy correction artificial viscosity is best understood as a minimally dissipative, locally computed mechanism for enforcing a discrete entropy inequality. Its later developments show that the choice of viscous discretization, the possibility of multiple targeted viscosities, and the extension to finite differences are not peripheral refinements; they are integral to how entropy correction is made accurate, contact preserving, and compatible with high-order explicit solvers [2602.23210] [2604.03158] [2508.21226].

Source: https://www.emergentmind.com/topics/entropy-correction-artificial-viscosity