---
title: Well-Balanced Gas-Kinetic Scheme
url: https://www.emergentmind.com/topics/well-balanced-gas-kinetic-scheme-gks
type: topic
---

# Well-Balanced Gas-Kinetic Scheme

Well-balanced gas-kinetic scheme (GKS) denotes a family of kinetic-flux numerical formulations in which BGK-type interface evolution, source discretization, and, in recent variants, mesh motion are coupled so that prescribed steady states are preserved exactly or to roundoff. In the current literature, the equilibria of interest include hydrostatic isothermal states for gravitational gas dynamics, lake-at-rest states for shallow water equations (SWEs), and layered rest states for two-layer SWEs. The topic spans symplecticity-preserving BGK methods for gravitational hydrodynamics, unified gas-kinetic schemes (UGKS) for multiscale flow transport under gravity, well-balanced Navier–Stokes GKS with global reconstruction, compact high-order shallow-water GKS on unstructured meshes, and recent arbitrary Lagrangian–Eulerian (ALE) and adaptive-mesh extensions [1012.3369], [1608.08730], [1904.00178], [2306.04467], [2508.14216], [2510.14673].

## 1. Definition, equilibrium targets, and problem setting

The phrase *well-balanced* is not tied to a single PDE or a single equilibrium. In gravitational gas dynamics, Xiao, Cai, and Xu define a well-balanced scheme as one that evolves an isolated gravitational system under any initial condition to an isothermal hydrostatic equilibrium state and keeps that solution thereafter [1608.08730]. In shallow-water formulations, the target is typically the lake-at-rest solution, characterized by zero velocity and constant free-surface elevation. In the two-layer shallow-water setting, the preserved equilibrium is layered rest, with the upper-layer thickness constant and the lower-layer free surface balanced against bottom topography [2306.04467].

| System | Equilibrium target | Representative condition |
|---|---|---|
| Gravitational gas dynamics | Isothermal hydrostatic equilibrium | $\nabla p = -\rho \nabla \Phi$, $\mathbf{U}=0$ |
| Single-layer SWE | Lake-at-rest | $u=v=0$, $\eta=h+b=\mathrm{const}$ |
| Two-layer SWE | Layered rest | $U_1=V_1=U_2=V_2=0$, $h_2=\mathrm{const}$, $h_1+B=\mathrm{const}$ |
| ALE shallow water on moving mesh | Lake-at-rest on arbitrarily moving meshes | Discrete flux, source, mesh-motion, and bottom-advection terms cancel while preserving GCL |

For the single-layer inviscid SWE with bathymetry $b(x,y)$, the governing equations are
$$
\partial_t h + \partial_x (h u) + \partial_y (h v) = 0,
$$
$$
\partial_t (h u) + \partial_x \left(h u^2 + \frac{1}{2} g h^2\right) + \partial_y (h u v) = - g h \partial_x b,
$$
$$
\partial_t (h v) + \partial_x (h u v) + \partial_y \left(h v^2 + \frac{1}{2} g h^2\right) = - g h \partial_y b,
$$
with lake-at-rest given by $u=v=0$ and $\eta=h+b=\mathrm{const}$ [2510.14673]. For gas dynamics under gravity, the hydrostatic balance couples pressure to gravitational potential $\Phi$ through $\nabla p=-\rho\nabla\Phi$, and isothermal equilibria yield exponential density profiles [1608.08730], [1904.00178].

The central numerical difficulty is that naive discretizations typically do not preserve these equilibria because flux gradients, source terms, and, on moving meshes, geometry-induced terms do not cancel exactly at the discrete level. Well-balanced GKS addresses this by embedding the equilibrium structure directly into reconstruction variables, kinetic interface solutions, or both [1904.00178], [2508.14216].

## 2. Kinetic formulation and interface evolution

A defining feature of GKS is that fluxes are not obtained from an approximate Riemann solver alone, but from moments of a time-dependent distribution function. For gravitational gas dynamics, the BGK-type kinetic equation with external force is
$$
\frac{\partial f}{\partial t} + \boldsymbol{u}\cdot\nabla_{\boldsymbol{x}} f + \boldsymbol{g}\cdot\nabla_{\boldsymbol{u}} f = \frac{f^{+} - f}{\tau},
$$
where $\boldsymbol{g}=-\nabla\Phi$ and $\tau$ is the collision time [1608.08730]. In the one-dimensional notation used in that work,
$$
f_t + u f_x + \phi_x f_u = \frac{f^{+}-f}{\tau}.
$$
The corresponding time-evolving interface distribution is obtained from the analytic integral solution along phase-space characteristics:
$$
f(\boldsymbol{x},\boldsymbol{u},t) = e^{-t/\tau} f_0\!\left(\boldsymbol{x} - \boldsymbol{u} t - \frac{1}{2}\boldsymbol{g} t^2,\; \boldsymbol{u} - \boldsymbol{g} t \right)
+ \int_0^t e^{-(t-s)/\tau} f^{eq}\!\left(\boldsymbol{x} - \boldsymbol{u}(t-s) - \frac{1}{2}\boldsymbol{g}(t-s)^2,\; \boldsymbol{u} - \boldsymbol{g}(t-s), s \right) ds.
$$
This scale-dependent interface solution is the basis of UGKS multiscale behavior, since the ratio $\Delta t/\tau$ continuously connects free transport and hydrodynamic relaxation [1608.08730].

For SWE, the kinetic model is adapted to hydrostatic pressure through a Maxwellian-like equilibrium. In single-layer formulations,
$$
g = h \left(\frac{\lambda}{\pi}\right)\exp\left[-\lambda (\mathbf{u}-\mathbf{U})^2\right], \qquad \lambda=\frac{1}{G h},
$$
with collision invariants $\boldsymbol{\psi}=(1,u,v)^T$ and moments
$$
W=\int f\,\boldsymbol{\psi}\,d\Xi, \qquad (F,G)^T=\int f\,\boldsymbol{\psi}\,\mathbf{u}\,d\Xi
$$
[2508.14216]. In the ALE compact shallow-water scheme, the same kinetic logic is used with a Maxwellian-like equilibrium
$$
g = h (\lambda/\pi)\exp(-\lambda(\xi-U)^2), \qquad \lambda = 1/(g h),
$$
and interface moments yield both fluxes and interface states together with their time derivatives [2510.14673].

Several well-balanced GKS variants exploit the fact that kinetic evolution supplies more than a flux. The compact GKS for two-layer SWEs and the ALE compact GKS both use the interface distribution to obtain time-dependent interface macroscopic states, gradients, and temporal derivatives; these are then recycled into source-term discretization, gradient updates, and one-stage or multistage time integration [2306.04467], [2510.14673]. In the ALE formulation, face motion enters directly through the relative particle speed, and the face-normal kinetic flux becomes
$$
F_n^{ALE} = \int (\xi\cdot n - w_n)\,\psi\, f(x_f,t,\xi)\, d\xi,
$$
which is consistent with the hydrodynamic ALE flux $F(U)\cdot n - U w_n$ [2510.14673].

## 3. Mechanisms used to achieve well-balancedness

The literature contains several distinct, technically precise mechanisms for enforcing discrete balance.

The symplecticity-preserving BGK formulation models the gravitational potential as a piecewise step function with a jump at each cell interface. Particle transmission and reflection across the potential barrier are governed by single-particle energy conservation,
$$
\frac{1}{2}u^2+\Phi = \frac{1}{2}u'^2+\Phi',
$$
together with Liouville’s theorem and the symplecticity-preserving relation $u' du' = u du$. This yields exact high-order moment mappings across the potential jump. That construction is used to prove the necessity of an exact Maxwellian for preserving hydrostatic isothermal states, as well as conservation of total mass and total energy including gravitational potential energy [1012.3369].

A different route is taken by the well-balanced Navier–Stokes GKS with gravitational potential. There, an auxiliary variable $\alpha$ is introduced through
$$
\rho = \rho_{ref}\exp\!\left(-\frac{\phi}{\alpha}\right),
$$
so that $\alpha=RT$ becomes constant at isothermal hydrostatic equilibrium. The momentum equation is reformulated so that, at equilibrium, the numerical momentum flux and the numerical source term each vanish separately rather than merely balancing one another. The resulting well-balanced property therefore relies on global reconstruction in the variables $\{\alpha,U,T,\phi\}$ and on a flux-source splitting that enforces zero numerical flux and zero numerical source at the target state [1904.00178].

In shallow-water schemes, the dominant mechanism is hydrostatic reconstruction at the level of free-surface variables and source discretization. On arbitrary quadrilateral meshes with hanging nodes, the well-balanced GKS with space-time adaptive mesh refinement reconstructs the free-surface elevation $\widetilde{h}=h+B$ rather than $h$,
$$
h(\mathbf{x})=\widetilde{h}(\mathbf{x})-B(\mathbf{x}), \qquad
\nabla h(\mathbf{x})=\nabla \widetilde{h}(\mathbf{x})-\nabla B(\mathbf{x}),
$$
and uses a source discretization on two triangular subcells inside each quadrilateral. On hanging edges, exact balancing under lake-at-rest requires two Gauss points rather than a single midpoint, because otherwise coarse-fine pressure-flux contributions become mismatched [2508.14216].

The two-layer shallow-water compact GKS extends this philosophy to coupled hydrostatic sources. Its target equilibrium satisfies
$$
U_1=V_1=U_2=V_2=0,\qquad h_2=\mathrm{const},\qquad h_1+B=\mathrm{const},
$$
and the source terms contain both bed-slope and interfacial pressure coupling. The method uses water-level reconstruction so that $\nabla(h_1+B)=\mathbf{0}$ and $\nabla h_2=\mathbf{0}$ at edge quadrature points, causing the discrete flux-gradient and source integrals to cancel within each control volume [2306.04467].

On moving meshes, well-balancedness requires one further compatibility relation: bathymetry is stationary in the physical frame but advects relative to the moving mesh according to
$$
\partial_t b + w\cdot\nabla b = 0.
$$
The ALE compact GKS therefore advances bottom topography with the same space-time discretization as the hydrodynamic variables. For lake-at-rest with piecewise linear $b$, the paper proves
$$
\mathrm{Tot}(L_{1,j}^n)+\mathrm{Tot}(L_{3,j}^n)=0
$$
for hydrodynamic and source terms, and
$$
\mathrm{Tot}(L_{2,j}^{h})+\mathrm{Tot}(L_{2,j}^{B})=\eta\,\Delta|\Omega_j|,
$$
so that $(h_j^{n+1}+B_j^{n+1})|\Omega_j^{n+1}|=\eta |\Omega_j^{n+1}|$ and the lake-at-rest state is preserved to roundoff on arbitrarily moving meshes [2510.14673].

## 4. High-order formulations, compact reconstruction, ALE geometry, and adaptivity

Well-balancedness in GKS is now tightly coupled to high-order space-time discretization. The two-layer compact GKS on triangular meshes uses a compact Hermite-type reconstruction based on cell averages and cell-averaged gradients. For fourth-order spatial accuracy it reconstructs a complete degree-3 polynomial
$$
P^3(\boldsymbol{x})=\sum_{k=0}^{9} a_k \varphi_k(\boldsymbol{x}),
$$
and combines it with lower-order polynomials through a compact WENO-like nonlinear blending. The same kinetic evolution supplies the time-dependent interface states needed to update cell-averaged gradients by Gauss’s theorem [2306.04467].

The recent ALE shallow-water scheme adopts a nonlinear fourth-order compact reconstruction on unstructured moving meshes. Its cubic polynomial
$$
P^3(x)=\sum_{k=1}^{10} a_k \phi_k(x)
$$
is obtained from constrained least squares on a compact stencil consisting of face- and vertex-neighbor cells. Because the interface kinetic model provides $f^n$, $f_t^n$, $W^n$, $W_t^n$, $F^n$, and $F_t^n$ in one stage, the method attains second-order temporal accuracy with the one-stage update
$$
U_j^{n+1} |\Omega_j^{n+1}| = U_j^n |\Omega_j^n| + L_j^n \Delta t + \frac{dL_j^n}{dt}\frac{\Delta t^2}{2}.
$$
This is paired with a direct physical-mesh update, so the method avoids remapping [2510.14673].

ALE consistency is not only a matter of replacing $F\cdot n$ by $F\cdot n-Uw_n$. The moving-mesh formulation also requires exact accounting of the swept face area during each time step. In the cited ALE compact GKS, straight faces with linearly moving endpoints yield an explicit geometric contribution involving
$$
(V_1 \times V_2)\cdot k
$$
in the time derivative of the mesh-motion operator. For a constant state, the total mesh-motion-induced update on a moving edge is exactly equal to the swept area $\Delta S_l$, which is the discrete GCL identity [2510.14673].

Adaptive mesh refinement introduces a different geometric complication. The well-balanced shallow-water GKS with space-time adaptive mesh refinement is built on arbitrary quadrilateral meshes with hanging nodes, managed by a quadtree structure via p4est with Morton indexing and strict 2:1 balance. Temporal adaptivity is achieved by level-dependent local time stepping:
$$
\Delta t = 2^{l_{\max}} \Delta t_{\min},\qquad
\Delta t_l = 2^{l_{\max}-l}\Delta t_{\min},\qquad
N_l = 2^l.
$$
Coarse-fine interface fluxes are stored during subcycling and used for conservative synchronization when the coarse cell reaches the corresponding intermediate time [2508.14216].

Across these formulations, compactness is repeatedly emphasized. In the SWE literature, compact GKS means that only immediate neighbors are needed for reconstruction, while discontinuous left/right interface states and their spatiotemporal derivatives are produced by the same kinetic evolution. This reduces stencil width on unstructured meshes and directly couples reconstruction, flux evaluation, and source treatment [2306.04467], [2510.14673].

## 5. Numerical behavior and validated application domains

The empirical record reported in the cited papers covers hydrostatic preservation, small-perturbation propagation, strong discontinuities, wet/dry fronts, moving meshes, and multiscale non-equilibrium transport.

For gravitational gas dynamics, UGKS preserves hydrostatic backgrounds while resolving small perturbations and regime transitions. In the one-dimensional pressure-perturbation test with $\phi_x=-1$, $\rho_0=p_0=\exp(-x)$, and $\mathrm{Kn}=10^{-4}$, the method maintains the hydrostatic background and resolves the perturbation. In the gravity-modified Sod problem, simulations at $\mathrm{Kn}=10^{-4}$, $10^{-2}$, and $1$ display the continuous transition from Euler to collisionless Boltzmann behavior. In the lid-driven cavity under gravity, the method reports a non-equilibrium phenomenon in the transition regime: heat flux correlated with gravity, flowing from the colder upper region to the hotter lower region, contrary to Fourier’s law [1608.08730].

For Navier–Stokes under gravity, the auxiliary-variable well-balanced GKS preserves hydrostatic solutions exactly in one-dimensional tests with the potentials $\phi_1(x)=x+1$, $\phi_2(x)=x^2+1$, and $\phi_3(x)=\sin(2\pi x)+2$ on 100 uniform cells up to $t=2.0$. In long-time viscous runs from a highly non-balanced initial condition, the system converges to zero velocity and constant temperature, with final velocities below $10^{-14}$ at $t=1000$. The same formulation accurately resolves small perturbations riding on equilibrium; normalized solutions for perturbation amplitudes $\eta=10^{-5}$ and $10^{-13}$ collapse, with deviations at $\eta=10^{-14}$ occurring at approximately $10^{-16}$, consistent with double-precision round-off [1904.00178].

In shallow-water applications, machine-precision lake-at-rest preservation is reported on both adaptive and moving meshes. With dynamic AMR on arbitrary quadrilateral meshes with hanging nodes, the well-balanced GKS produces $L_1(h)\approx 10^{-15}$–$10^{-16}$ and $L_1(hU)$, $L_1(hV)\approx 10^{-15}$–$10^{-14}$ at $t=0.1$, $1.0$, and $10.0$. In a two-dimensional dam-break efficiency study, the uniform mesh uses 97,792 cells and 240 s CPU time, whereas AMR uses 18,682 cells and 100 s, and space-time AMR uses 18,682 cells and 60 s, corresponding to speedups of approximately $2.4\times$ and $4.0\times$ [2508.14216].

On moving unstructured meshes, the ALE compact GKS verifies both discrete GCL and well-balancedness. For lake-at-rest over the linear bottom
$$
b(x,y)=0.05+0.075(x+y),
$$
errors for $\eta$ and momentum stabilize below $10^{-9}$ on a moving mesh. For a nonlinear Gaussian hump, the $\eta$ and momentum errors stabilize below $10^{-3}$, with the slight discrepancy localized near the crest and attributed to the reconstruction of nonlinear bathymetry on moving cells. In a small-perturbation test over non-flat bottom, comparisons with a fixed-mesh reference show discrepancies $\le 6\times 10^{-4}$ in $\eta$ away from the bottom crest. In an irregular-domain dam-break, the adaptive moving mesh attains accuracy comparable to a refined fixed mesh near the dam region while using about one quarter of the degrees of freedom [2510.14673].

The two-layer compact GKS further demonstrates that well-balancedness is compatible with strongly coupled hydrostatic sources, internal-wave dynamics, and wet/dry fronts on triangular meshes. The paper reports well-balanced lake-at-rest errors below $10^{-8}$ up to $t=100$, accurate resolution of weak and strong discontinuity Riemann problems, good agreement with references in equal-density and light-over-dense dam-breaks, and robust bore capture in irregular wet and dry bed configurations [2306.04467].

## 6. Limitations, common misconceptions, and extensions

A common misconception is to treat well-balancedness as synonymous with general robustness, positivity preservation, or exact shock resolution. The cited literature does not support that equivalence. Several shallow-water papers emphasize robustness near bores and wet/dry fronts but state that positivity-preserving details are not the focus or that no formal positivity-preserving proof is given [2306.04467], [2508.14216], [2510.14673]. Conversely, the Navier–Stokes gravity formulation explicitly notes that its momentum discretization is non-conservative because part of the gravitational source is merged into the flux; small defects appear near rarefactions in shock-tube comparisons, and the recommended application range is low-speed continuous flows [1904.00178].

Another misconception is that well-balancedness is only a stationary-mesh issue. The ALE shallow-water formulation shows that equilibrium preservation on moving meshes additionally requires compatibility with the geometric conservation law and with bottom-topography advection relative to the mesh. Uniform-state preservation is tied to exact swept-area accounting, and lake-at-rest preservation on moving cells requires the discrete mesh-motion and bathymetry updates to preserve free-surface elevation [2510.14673].

The present literature also makes clear that the target equilibrium matters. The UGKS construction under gravity targets isothermal hydrostatic equilibrium; extensions to non-isothermal steady states would require modified reconstruction that embeds the appropriate steady relations, such as simultaneous gravitational balance and conductive equilibrium [1608.08730]. Similarly, the auxiliary-variable Navier–Stokes GKS is designed around the fact that $\alpha$ becomes constant at isothermal hydrostatic equilibrium, not at arbitrary stratified steady states [1904.00178].

Methodological limitations are formulation-specific. In the shallow-water AMR scheme, the implementation is second-order, and higher-order extensions would require more complex reconstructions; the dissipation parameters in
$$
\tau=C_1\Delta t+C_2\left|\frac{h_l^2-h_r^2}{h_l^2+h_r^2}\right|\Delta t
$$
may require tuning, and extreme wet/dry topographies may remain challenging [2508.14216]. In the ALE compact GKS, residual terms associated with spatial nonuniformity and face-length or normal changes are neglected in $dL_{1,j}^n/dt$ for simplicity, and for nonlinear bathymetry exact piecewise-linear reconstruction is not maintained after moving-mesh updates, which leads to bounded but nonzero free-surface errors of at most $O(10^{-3})$ in the reported tests [2510.14673]. In the two-layer compact GKS, the source discretization is second-order in space for efficiency, which limits global convergence order despite fourth-order reconstruction and a two-stage fourth-order temporal update [2306.04467].

At the same time, the extension paths are explicit in the cited works. The ALE shallow-water paper states that the methodology generalizes to three-dimensional SWEs and to compressible Euler/Navier–Stokes equations on moving meshes [2510.14673]. The AMR paper notes that high-order compact GKS ideas can be adapted to SWE and that additional physics such as friction, wind stress, and Coriolis terms follow the same source-in-flux framework [2508.14216]. The UGKS under gravity points toward non-isothermal equilibria, radiative coupling, and multi-species kinetics, while the two-layer compact GKS remarks that multilayer SWEs, friction, Coriolis forcing, adaptive mesh refinement, and implicit time stepping are natural continuations [1608.08730], [2306.04467].

Taken together, the well-balanced GKS literature shows a consistent pattern: equilibrium preservation is achieved not by a single universal trick, but by matching the kinetic interface model, reconstruction variables, source quadrature, and geometric treatment to the target steady state. The resulting schemes retain the characteristic GKS advantage of time-dependent interface evolution while extending that structure to gravitational forcing, non-flat bathymetry, unstructured high-order reconstruction, moving meshes, and space-time adaptivity [1012.3369], [1608.08730], [1904.00178], [2306.04467], [2508.14216], [2510.14673].

Source: https://www.emergentmind.com/topics/well-balanced-gas-kinetic-scheme-gks