pub fn bgk_equilibrium(
rho: f64,
ux: f64,
uy: f64,
w: f64,
cx: f64,
cy: f64,
) -> f64Expand description
Compute the BGK equilibrium distribution for one velocity direction.
Formula: f_eq = w · ρ · (1 + (c·u)/cs² + (c·u)²/(2cs⁴) − u²/(2cs²))
rho– macroscopic densityux,uy– macroscopic velocity componentsw– lattice weight for this directioncx,cy– discrete velocity components for this direction