Effective Viscous Flux in Compressible Flow

• Effective viscous flux is a scalar quantity combining pressure and variable viscosity effects, enhancing the compactness analysis in compressible flows.
• It leverages the Hodge decomposition to segregate the gradient component of the Newtonian force from divergence-free motions, clarifying its regularization role.
• Its explicit commutator structure, emerging from density mollification, underpins uniform energy estimates and the proof of global weak solutions.

Searching arXiv for the cited paper and closely related work on effective viscous flux in compressible Navier–Stokes. Effective viscous flux is a distinguished scalar combination of the pressure and the Newtonian forcing in the compressible Navier–Stokes system, used to extract compactness and regularity that are not apparent from the primitive variables alone. In the setting studied by Frid, Marroquin, and Nariyoshi, the fluid is barotropic on a smooth bounded domain $\Omega \subset \mathbb{R}^n$, $n \ge 2$, with pressure law $P(\rho)=A\rho^\gamma$, $A>0$, $\gamma>n/2$, and Newtonian stress tensor $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$, where both viscosities depend on a spatial mollification of the density. The effective viscous flux is then defined by

$$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$

and it is characterized through the Hodge decomposition of the Newtonian force. In this formulation, its gradient is exactly the gradient part of the forcing, and its weak continuity becomes the mechanism that restores strong compactness of the density in global weak-solution theory (Frid et al., 2020).

1. Governing equations and formal definition

The underlying system is the compressible Navier–Stokes equations for the density $\rho=\rho(t,x)\ge 0$ and velocity $u=u(t,x)\in\mathbb{R}^n$:

$\partial_t\rho+\operatorname{div}(\rho u)=0,$

$n \ge 2$0

subject to the no-slip boundary condition $n \ge 2$1 and initial data $n \ge 2$2, $n \ge 2$3. The pressure is barotropic, $n \ge 2$4, and the stress is Newtonian with density-dependent viscosities evaluated at the mollified density $n \ge 2$5, where $n \ge 2$6 is a standard mollifier and $n \ge 2$7 is extended by zero outside $n \ge 2$8 (Frid et al., 2020).

The viscosity coefficients satisfy

$n \ge 2$9

with $P(\rho)=A\rho^\gamma$0. Within this framework, the effective viscous flux is introduced by applying the inverse Laplacian to the divergence of the Newtonian stress and subtracting the resulting scalar potential from the pressure:

$P(\rho)=A\rho^\gamma$1

In the constant-viscosity case $P(\rho)=A\rho^\gamma$2, this reduces to the classical Lions–Feireisl flux

$P(\rho)=A\rho^\gamma$3

This definition isolates the scalar part of the viscous forcing that interacts most effectively with the pressure. A plausible implication is that the flux is not merely a bookkeeping device: it reorganizes the momentum equation into a form where the analytically favorable component is explicit.

2. Hodge decomposition and the gradient part of the Newtonian force

Let $P(\rho)=A\rho^\gamma$4. The Helmholtz–Hodge decomposition gives

$P(\rho)=A\rho^\gamma$5

with

$P(\rho)=A\rho^\gamma$6

Equivalently,

$P(\rho)=A\rho^\gamma$7

Substituting this splitting into the momentum equation yields

$P(\rho)=A\rho^\gamma$8

or, in the effective-flux notation,

$P(\rho)=A\rho^\gamma$9

Thus $A>0$0 is exactly the gradient part of the Newtonian forcing $A>0$1 together with the pressure gradient (Frid et al., 2020).

This characterization is central because it separates the compressive component from the divergence-free remainder. The paper explicitly frames the regularizing properties of the effective viscous flux through this decomposition. This suggests that the flux is best understood as a structural object attached to the Helmholtz–Hodge split rather than as an ad hoc scalar combination.

3. Explicit representation and the commutator structure

The density-dependent-viscosity case differs from the constant-viscosity case by a commutator term. Using

$A>0$2

where $A>0$3 are Riesz transforms, the paper obtains the explicit identity

$A>0$4

Here

$A>0$5

is the commutator. When the shear viscosity is constant, $A>0$6, the commutator vanishes and one recovers

$A>0$7

(Frid et al., 2020).

The significance of this formula lies in showing exactly how variable viscosities alter the classical flux. The new term is not an arbitrary correction but a Riesz-transform commutator generated by the spatial dependence of the viscosity through $A>0$8. A common misunderstanding is to treat the effective viscous flux as always equal to $A>0$9; that identity is valid only in the constant-shear-viscosity case. In the present setting, the commutator is part of the definition at the level of explicit representation.

4. A priori estimates and regularizing effect

The analysis begins from the basic energy inequality, which provides the controls

$\gamma>n/2$0

A Bogovskii test function in the momentum equation yields the higher integrability

$\gamma>n/2$1

uniformly in the approximations. The paper then emphasizes uniform bounds on $\gamma>n/2$2 and its gradient. The principal observation is that $\gamma>n/2$3 solves a perturbed transport–elliptic equation whose right-hand side is more regular, or at least weakly continuous in the limit. Formally,

$\gamma>n/2$4

and the relevant terms converge weakly under the natural compactness in

$\gamma>n/2$5

and related spaces (Frid et al., 2020).

The regularizing role of the effective viscous flux is encoded in its weak continuity. If $\gamma>n/2$6 weakly in the natural spaces, then for any $\gamma>n/2$7,

$\gamma>n/2$8

According to the paper, this continuity of $\gamma>n/2$9 is exactly what restores compactness of the density and enforces

$$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$0

in the weak limit.

5. Role in global weak existence

The principal existence statement is Theorem 1.1: if $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$1, $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$2, $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$3, and $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$4 satisfy the mollified-density laws above, then there exists a global renormalized weak energy solution $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$5 of

$$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$6

with $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$7, satisfying the energy inequality and the continuity equation in DiPerna–Lions sense (Frid et al., 2020).

The proof uses a two-level regularization. Artificial viscosity $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$8 is added in the continuity equation and an artificial pressure $$S=\lambda([\rho]_\varepsilon)(\operatorname{div}u)I+2\mu([\rho]_\varepsilon)D(u)$$9 is added in $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$0; the $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$1–$$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$2 regularized system is solved by a standard Galerkin procedure. Uniform estimates are then derived, including energy bounds, higher integrability of $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$3, and bounds for $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$4 uniformly in $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$5. Passing $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$6 at fixed $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$7 relies on compactness and the weak continuity of the effective viscous flux, extended in Lemma 4.1. Passing $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$8 uses truncation and renormalization of the continuity equation together with the weak continuity identity for $$F_{\mathrm{eff}}:=P(\rho)-\operatorname{div}\Delta^{-1}[\operatorname{div}S],$$9, stated in Lemma 5.2, to obtain strong convergence of $\rho=\rho(t,x)\ge 0$0 and remove the artificial pressure.

Within this proof architecture, the effective viscous flux is the central analytical device. The paper states explicitly that without the extra compactness afforded by the weak continuity of $\rho=\rho(t,x)\ge 0$1, one cannot recover strong convergence of $\rho=\rho(t,x)\ge 0$2 and hence cannot conclude $\rho=\rho(t,x)\ge 0$3.

6. Historical placement and analytical significance

The regularizing properties of the effective viscous flux were known from the pioneering works of Serre, Hoff, Vaigant–Kazhikhov, Lions, and Feireisl, among others. In that earlier line of analysis, the flux is characterized as the function whose gradient is the gradient part in the Hodge decomposition of the Newtonian force when the shear viscosity is constant. The present work extends that connection by addressing global weak existence for compressible Navier–Stokes equations in which both viscosities depend on a spatial mollification of the density (Frid et al., 2020).

The analytical novelty is therefore not the abandonment of the classical mechanism but its persistence under a modified constitutive law. In constant-viscosity theories, the effective viscous flux is often presented in the simpler scalar form $\rho=\rho(t,x)\ge 0$4. Here, the same concept survives, but its explicit realization includes the commutator generated by variable viscosity. This suggests continuity between the classical Lions–Feireisl framework and the mollified-density-viscosity setting: the compactness mechanism remains the same in principle, while the representation of the flux becomes more intricate.

7. Conceptual summary and common points of emphasis

The paper organizes the role of the effective viscous flux around a small number of structural facts. First, it is defined through the Hodge decomposition of $\rho=\rho(t,x)\ge 0$5. Second, it admits an explicit commutator formula when viscosities depend on $\rho=\rho(t,x)\ge 0$6. Third, it satisfies a weak continuity property strong enough to recover compactness of the density. Fourth, that compactness is indispensable for identifying the weak limit of the pressure and proving global existence of renormalized weak energy solutions (Frid et al., 2020).

Several points deserve emphasis. The effective viscous flux is not identical to the pressure, nor merely to the divergence of the velocity; it is a specific pressure–viscosity combination adapted to the Newtonian force. Its gradient is the relevant object in the Hodge decomposition, while the scalar flux itself is the quantity used in compactness arguments. The commutator term does not represent a secondary perturbation in the explicit formula; in the variable-viscosity setting it is part of the exact identity. Finally, the weak continuity of $\rho=\rho(t,x)\ge 0$7 is not an auxiliary estimate but the mechanism that restores the strong convergence of $\rho=\rho(t,x)\ge 0$8 required for the passage to the limit.

In this formulation, effective viscous flux functions as the bridge between elliptic structure, transport dynamics, and compactness in compressible-flow existence theory.