Boundary effects on energy dissipation in a cellular automaton model

In this paper, we numerically study energy dissipation caused by traffic in the Nagel-Schreckenberg (NaSch) model with open boundary conditions (OBC). Numerical results show that there is a nonvanishing energy dissipation rate Ed, and no true free-flow phase exists in the deterministic and nondeterministic NaSch models with OBC. In the deterministic case, there is a critical value of the extinction rate $\beta{cd}$ below which Ed increases with increasing $\beta$, but above which Ed abruptly decreases in the case of the speed limit vmax>2. However, when vmax<3, no discontiguous change in Ed occurs. In the nondeterministic case, the dissipated energy has two different contributions: one coming from the randomization, and one from the interactions, which is the only reason for dissipating energy in the deterministic case. The relative contributions of the two dissipation mechanisms are presented in the stochastic NaSch model with OBC. Energy dissipation rate Ed is directly related to traffic phase. Theoretical analyses give an agreement with numerical results in three phases (low-density, high-density and maximum current phase) for the case vmax=1.

💡 Research Summary

The paper investigates energy dissipation in the Nagel‑Schreckenberg (NaSch) cellular automaton traffic model when open boundary conditions (OBC) are imposed. In contrast to the widely studied periodic‑boundary version, OBC introduces an entry probability α at the upstream boundary and an exit (extinction) probability β at the downstream boundary, thereby mimicking realistic road segments with on‑ramps and off‑ramps. The authors perform extensive numerical simulations for both deterministic (p = 0) and stochastic (p > 0) versions of the model, where p denotes the randomization probability that reduces a vehicle’s speed by one cell with probability p each time step.

Key findings for the deterministic case are as follows. Even without randomization, a non‑zero energy dissipation rate Ed persists because vehicles experience abrupt speed changes when they are created at the entrance or removed at the exit. When the maximum allowed speed vmax exceeds 2, Ed exhibits a discontinuous dependence on β: for β below a critical value βcd, Ed grows monotonically with β, reflecting the increasing frequency of vehicles leaving the system; once β surpasses βcd, Ed drops sharply because the downstream bottleneck is relieved, reducing the number of braking events. If vmax ≤ 2, this abrupt transition disappears and Ed varies smoothly with β, indicating that low speed limits suppress the formation of strong boundary‑induced shocks.

In the stochastic case (p > 0) the total dissipation splits into two additive contributions. The first originates from the randomization rule itself: each stochastic deceleration dissipates kinetic energy regardless of interactions. The second stems from vehicle‑vehicle interactions (braking due to insufficient headway), which is the sole source of dissipation in the deterministic limit. By varying p, α, and β the authors quantify the relative weight of these two mechanisms. At low densities (small α) the randomization term dominates, whereas at high densities (large α, small β) interaction‑induced dissipation remains significant even for sizable p.

A special analytical treatment is presented for the simplest case vmax = 1, where each vehicle can only move forward by one cell per time step. In this limit the phase diagram reduces to three well‑known regimes: low‑density (LD), high‑density (HD), and maximum‑current (MC) phases. The authors derive explicit expressions for Ed in each regime. In the LD phase Ed is proportional to pα because vehicles rarely interact and energy loss is almost entirely due to randomization. In the HD phase Ed depends mainly on β, reflecting that the exit bottleneck forces frequent braking near the downstream boundary. In the MC phase both α and β are large enough to sustain the maximal flow; here Ed becomes independent of α and β and is set solely by p, i.e., by the intrinsic stochasticity of driver behavior. The analytical results match the simulation data with high accuracy, confirming the validity of the theoretical approach.

Overall, the study demonstrates that open boundaries eliminate a true free‑flow phase: even in the deterministic limit a finite Ed persists because of boundary‑induced speed adjustments. The magnitude and nature of energy dissipation are controlled by three intertwined factors: (i) the maximum speed limit vmax, which determines whether discontinuous transitions in Ed can occur; (ii) the randomization probability p, which adds a background dissipation independent of traffic density; and (iii) the boundary rates α and β, which set the global density and thus the balance between interaction‑driven and randomization‑driven losses.

The implications are twofold. From a theoretical perspective, the work extends the NaSch model’s phase‑diagram analysis to realistic open systems and clarifies how microscopic braking events translate into macroscopic energy loss. From an applied viewpoint, the results suggest that traffic management strategies aimed at reducing energy consumption should focus not only on smoothing interior flow (e.g., via adaptive cruise control) but also on mitigating boundary bottlenecks—by adjusting ramp metering rates (α) or improving downstream capacity (β). Future research directions proposed include extending the analysis to multi‑lane networks, incorporating heterogeneous driver behavior, and coupling the cellular automaton with real‑time traffic signal control to evaluate how dynamic boundary conditions influence overall energy efficiency.