Flux-Balance Laws: Principles & Applications

• Flux-Balance Laws are mathematical relations that equate the time change of a density with the divergence of a flux and accompanying source terms, fundamental in fields like continuum mechanics and asymptotic gravity.
• They underpin analysis in PDE theory by linking entropy production, weak formulations, and numerical stability, as seen in studies of hyperbolic balance and conservation laws.
• Applied in areas from self-force calculations in Kerr spacetimes to constraint-based metabolic models, these laws translate physical symmetries and resource allocations into actionable balance equations.

Searching arXiv for recent and foundational papers on flux-balance laws across self-force, asymptotic gravity, PDEs, and constraint-based modeling. Searching "flux-balance laws self-force Kerr Carter constant" Flux-balance laws are relations that connect the change of a density, charge, or orbital invariant to flux through the relevant boundary and, when present, to source or production terms. In one common local form they appear as $\partial_t u+\nabla\cdot F(u)=S(u,x,t)$, while in asymptotic gravity they relate variations of charges on cuts of $\mathscr I^+$ to radiative fluxes, and in self-force theory they relate the secular drift of Kerr invariants to radiation through future null infinity and the future event horizon (Abgrall et al., 2021, Fiorucci et al., 30 Apr 2025, Grant et al., 2022). A recurrent theme is that the same system may admit a true conservation law only in special symmetry or source-free regimes; otherwise one obtains a balance law.

1. Definition and structural principles

In continuum formulations, a balance law expresses the temporal change of a density plus the divergence of a flux as a source or production term. In Rational Extended Thermodynamics, a balance system is written as

$(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$

or equivalently

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$

with $F_i^0$ the densities, $F_i^A$ the flux components, and $\Pi_i$ the production terms (Preston, 2010). Within that framework, “supplementary balance laws” are additional laws $(K^\mu)_{,\mu}=Q$ satisfied by all solutions; they include the original balances, the entropy balance, and symmetry-induced balances, and they are characterized by the Lagrange–Liu relations

$d_y K^\mu=\lambda^i\,d_y F_i^\mu.$

The same formalism links entropy density, main fields, and production through

$$\frac{\partial h^0}{\partial w^j}=\lambda^j,\qquad E=\lambda^i\Pi_i\ge 0$$

(Preston, 2010).

This broad viewpoint already shows why the term “flux-balance law” is not restricted to one discipline. In some settings the law is local and differential, in others it is global and boundary-based, and in others it is an emergent law for averaged quantities. A common misconception is to treat “balance law” and “conservation law” as interchangeable. The papers considered here consistently distinguish them: sources, explicit coordinate dependence, radiation leakage through asymptotic boundaries, hyperstress terms, or broken scaling symmetry convert conservation into balance (Lazar et al., 2013, Preston, 2010).

2. Hyperbolic balance laws, entropy, and weak formulations

In PDE theory, balance laws are commonly studied in the form

$\mathscr I^+$0

or, in one space dimension,

$\mathscr I^+$1

on the quarter-plane (Colombo et al., 26 Mar 2025, Sahoo et al., 2022). For the latter problem, introducing

$\mathscr I^+$2

transforms the equation into a conservation law for $\mathscr I^+$3, and the solution admits a Lax–Oleinik type representation

$\mathscr I^+$4

with admissible sets $\mathscr I^+$5 and $\mathscr I^+$6 controlling whether minimizing curves reach the initial line or the boundary $\mathscr I^+$7 (Sahoo et al., 2022). The boundary analysis is nontrivial because minimizing curves may spend time on the boundary, and the paper classifies boundary points into Types I–III in order to verify the Bardos–le Roux–Nédélec condition.

For general scalar balance laws, the entropy-production viewpoint packages all entropy inequalities into a single operator. Given a bounded function $\mathscr I^+$8, the entropy production $\mathscr I^+$9 acts on an entropy $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$0; for distributional solutions and $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$1,

$(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$2

where $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$3 are the Kružkov entropy productions (Colombo et al., 26 Mar 2025). The same work proves that positivity of all $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$4 is equivalent to positivity of $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$5 for every convex entropy, and extends the framework to space- and time-dependent or complex-valued entropies. This places entropy admissibility, weak solvability, and Fourier representations of entropy production in a single operator-theoretic framework.

The integral balance interpretation of weak solutions is also sharpened by flux regularity results. For nonlinear hyperbolic conservation laws, the time-integrated boundary flux is shown to be Lipschitz continuous with respect to suitable perturbations of the boundary, so that a weak solution indeed satisfies the integral balance law on smooth bounded domains (Ben-Artzi et al., 2020). In one-dimensional strictly hyperbolic systems with source, vanishing-viscosity solutions are SBV for all but at most countably many times under genuine nonlinearity, and this extends to the broader non-degeneracy condition; when linear degeneracy is present, full SBV regularity fails but an SBV-like statement survives at the level of the eigenvalue functions (Ancona et al., 2024). Stability theory also extends beyond local fluxes: for systems with spatially nonlocal interactions and temporal memory, entropy solutions satisfy quantitative $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$6 stability estimates with respect to perturbations in fluxes, kernels, sources, and initial data (Aggarwal et al., 3 May 2026).

3. Self-force theory and Kerr flux-balance laws

In scalar self-force theory on Kerr, a radiating point particle of scalar charge $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$7 and dynamical mass $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$8 moves on a worldline $(F_i^\mu(s(x)))_{,\mu}=\Pi_i(s(x)),$9, and the field satisfies

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$0

Using the Detweiler–Whiting decomposition $F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$1, the first-order scalar self-force is

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$2

and the instantaneous drift of the symmetry-generated integrals obeys

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$3

while for the Carter constant

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$4

(Grant et al., 2022). The central difficulty is that the Carter constant is not directly associated with a Noether current built from a spacetime isometry.

The resolution proceeds through symplectic currents, symmetry operators, and action-angle variables. The scalar-field symplectic current is

$F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$5

and bound Kerr geodesics admit action-angle coordinates $F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$6 with $F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$7 and $F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$8 (Grant et al., 2022). This yields exact integrated balance laws for the actions and, after orbit averaging and replacing $F_i^0{}_{,t}+F_i^A{}_{,x_A}=\Pi_i,$9 by the globally defined radiative field $F_i^0$0, explicit asymptotic mode sums at $F_i^0$1 and $F_i^0$2. With $F_i^0$3 and

$F_i^0$4

the averaged law for the Carter constant becomes

$F_i^0$5

which supplies a flux-balance law for $F_i^0$6 even though no direct stress-energy Noether current reproduces $F_i^0$7 (Grant et al., 2022).

The gravitational spinning-body extension uses the regular metric perturbation and the gravitational symplectic current instead of a scalar field. In the pole–dipole approximation with the Tulczyjew–Dixon spin supplementary condition, the paper derives a quasi-local averaged law

$F_i^0$8

and, off resonance, an invariant flux law

$F_i^0$9

for the four geodesic-like invariants $F_i^A$0 (Grant, 2024). Here $F_i^A$1 is the generalized Carter constant including the Killing–Yano correction. The fifth invariant, the Rüdiger constant, remains coupled to constraint variables rather than yielding a practical flux law. This suggests that, at linear order in spin, the generalized Carter constant is accessible through asymptotic radiation, whereas the full five-invariant system is not yet closed in a practical boundary-flux form.

4. Asymptotic gravitational charges, memory, and extensions beyond flat space

At future null infinity of asymptotically flat spacetimes, flux-balance laws govern the Bondi mass aspect $F_i^A$2, the angular momentum aspect $F_i^A$3, and the higher BMS charges built from the shear $F_i^A$4 and news $F_i^A$5. In a Carrollian-holographic derivation, local Carroll, Weyl, and diffeomorphism invariance at $F_i^A$6 imply the Bondi evolution equations

$F_i^A$7

together with the corresponding angular-momentum balance equation, and the cut charge

$F_i^A$8

obeys a flux law sourced by the news (Fiorucci et al., 30 Apr 2025). This is a boundary-intrinsic derivation of the same balance structure usually obtained from bulk Einstein equations.

Within the standard Bondi framework, the Poincaré and proper BMS flux-balance laws distinguish supermomentum, super-angular momentum, and super-center-of-mass charges, and they encode the displacement, spin, and center-of-mass memory effects defined from the shear (Compère et al., 2019). Their multipolar realization is exact in terms of radiative symmetric tracefree moments. In that formulation, fluxes of energy, angular momentum, and octupole super-angular momentum arise at $F_i^A$9, fluxes of quadrupole supermomentum arise at $\Pi_i$0, and fluxes of momentum, center-of-mass, and octupole super-center-of-mass arise at $\Pi_i$1 (Compère et al., 2019). The same work also scrutinizes the prescriptions for angular momentum and center-of-mass and singles out the covariant choice that aligns with the symplectic framework.

The covariant phase-space and Wald–Zoupas formalism extends these laws beyond general relativity. In diffeomorphism-invariant extensions such as luminal Horndeski theory, charges on cuts $\Pi_i$2 satisfy

$\Pi_i$3

and the flux receives both the general-relativistic news contribution and a scalar contribution proportional to $\Pi_i$4 (Maibach et al., 11 Jan 2026). In Brans–Dicke theory this reduces to an explicit scalar-flux correction to the supertranslation law. A related cosmological deformation occurs in de Sitter spacetime: at quadrupolar order the $\Pi_i$5 flux-balance laws govern dilatations, rotations, spatial translations, and cosmological boosts, and the flat limit recovers the standard Poincaré laws at future null infinity (Compère et al., 2024). The de Sitter dilatation flux can be written in two distinct negative-definite forms, making precise how curvature deforms the notion of radiative energy.

5. Constraint-based metabolic formulations

In systems biology, the phrase “flux-balance” refers to stoichiometric mass balance under metabolic constraints. Standard flux balance analysis imposes

$\Pi_i$6

and typically maximizes the biomass-reaction flux $\Pi_i$7 (Mori et al., 2016). Constrained Allocation Flux Balance Analysis adds a proteome-wide constraint derived from growth-law data,

$\Pi_i$8

where the coefficients encode catabolic, enzymatic, and ribosomal proteome costs (Mori et al., 2016). In this usage, flux-balance laws are not asymptotic charge laws but steady-state mass-balance constraints augmented by resource-allocation laws. The model predicts a respiratory-to-fermentative crossover and quantitatively relates growth rate, uptake, secretion, and yield.

Dynamic enzyme-cost Flux Balance Analysis extends this setting by adding slow-time balances for external species and macromolecules while maintaining quasi-steady internal metabolism:

$\Pi_i$9

Enzyme capacities appear as linear constraints

$(K^\mu)_{,\mu}=Q$0

and robust deFBA handles uncertainty through a scenario tree with non-anticipativity in the first interval and a robust horizon $(K^\mu)_{,\mu}=Q$1 (Lindhorst et al., 2016). This suggests a distinct but mathematically recognizable interpretation of flux-balance laws: balance of stoichiometric fluxes, enzyme capacities, and biomass synthesis under uncertainty rather than boundary-flux accounting.

6. Continuum-mechanical and numerical realizations

In micromorphic elasticity, translational, rotational, and dilatational symmetries generate fluxes whose divergences are the configurational balance laws of the medium. The translational flux is the Eshelby stress tensor

$(K^\mu)_{,\mu}=Q$2

the rotational flux is the angular momentum tensor $(K^\mu)_{,\mu}=Q$3, and the scaling flux is the dilatational vector $(K^\mu)_{,\mu}=Q$4 (Lazar et al., 2013). Their surface integrals define the $(K^\mu)_{,\mu}=Q$5-, $(K^\mu)_{,\mu}=Q$6-, and $(K^\mu)_{,\mu}=Q$7-integrals. In homogeneous isotropic media without body forces or couples, the $(K^\mu)_{,\mu}=Q$8- and $(K^\mu)_{,\mu}=Q$9-integrals become conservation laws, whereas the dilatational law is generally broken by the hyperstress term $d_y K^\mu=\lambda^i\,d_y F_i^\mu.$0 (Lazar et al., 2013). This is an explicit example of the general principle that internal length scales or constitutive structure convert a conservation law into a balance law.

Numerical analysis recasts the same principle at the discrete level. For multidimensional linear hyperbolic balance laws, the Global Flux formulation rewrites the spatial operator through one-dimensional and mixed derivatives of global primitives, for example

$d_y K^\mu=\lambda^i\,d_y F_i^\mu.$1

and the resulting finite-element schemes are stationarity preserving because the stabilizations act on the same global-flux combinations that define the steady kernel (Barsukow et al., 3 Oct 2025). With Gauss–Lobatto nodes, the line and row integrators are Lobatto IIIA tables and yield nodal steady-state errors $d_y K^\mu=\lambda^i\,d_y F_i^\mu.$2 for $d_y K^\mu=\lambda^i\,d_y F_i^\mu.$3 (Barsukow et al., 3 Oct 2025). Residual distribution schemes implement the same balance at the element level through

$d_y K^\mu=\lambda^i\,d_y F_i^\mu.$4

or through boundary numerical fluxes, and the framework encompasses SUPG, DG, finite-volume, and well-balanced shallow-water discretizations (Abgrall et al., 2021). A plausible implication is that, in numerical settings as in analytic ones, a flux-balance law is most effective when the discrete operator, the source treatment, and the steady-state kernel are made exactly compatible.