Kinematic waves are solutions to the conservation equation ∂E/∂t + ∇·I = 0, where E is a scalar density field and I is its flux. In one dimension, this becomes ∂E/∂t + ∂I/∂x = 0. For a gas, E = ρ (density) and I = ρv (flux), leading to ∂ρ/∂t + ∂(ρv)/∂x = 0. However, this equation alone is insufficient to determine both ρ and v. Assuming I depends only on E, the equation becomes ∂E/∂t + (dI/dE)∂E/∂x = 0, or ∂u/∂t + c(u)∂u/∂x = 0, where c(u) = dI/dE. This is a first-order quasi-linear PDE. The characteristics of this equation are given by dt/1 = dx/c(u) = du/0. Along characteristics, u is constant, and the solution is given implicitly as u = f[x - c(u)t]. The initial value problem requires specifying u(x, 0) = f(x). The solution is parametric, with x = ξ + c[f(ξ)]t and u = f(ξ). The derivatives of u with respect to t and x are derived, and substituting into the PDE confirms the solution. The validity of the solution depends on whether the solution breaks down. The derivatives become singular when t = -1/(dc/du f'(ξ)). If dc/du f'(ξ) > 0, characteristics do not intersect, but if dc/du f'(ξ) < 0, characteristics intersect, leading to a breakdown at t_b = 1/|dc/du f'(ξ)|_max. For dc/du > 0, the slope of the initial distribution steepens for f'(x) < 0 and flattens for f'(x) > 0. For dc/du < 0, the opposite occurs. An example shows that the solution becomes multivalued after t_b, making it physically invalid. In traffic flow, the density ρ(x,t) and speed v(x,t) are related by q = ρv (flux). The conservation equation is ∂ρ/∂t + ∂q/∂x = 0, with wave velocity c(ρ) = v + ρ ∂v/∂ρ. Since ∂v/∂ρ < 0, c < v, implying the wave propagates backward. For the Lincoln tunnel, data show q = aρ(1 - ρ/ρ_j), with a = 17.2 mph and ρ_j = 250 vpm. The wave velocity c(ρ) = a(1 - 2ρ/Kinematic waves are solutions to the conservation equation ∂E/∂t + ∇·I = 0, where E is a scalar density field and I is its flux. In one dimension, this becomes ∂E/∂t + ∂I/∂x = 0. For a gas, E = ρ (density) and I = ρv (flux), leading to ∂ρ/∂t + ∂(ρv)/∂x = 0. However, this equation alone is insufficient to determine both ρ and v. Assuming I depends only on E, the equation becomes ∂E/∂t + (dI/dE)∂E/∂x = 0, or ∂u/∂t + c(u)∂u/∂x = 0, where c(u) = dI/dE. This is a first-order quasi-linear PDE. The characteristics of this equation are given by dt/1 = dx/c(u) = du/0. Along characteristics, u is constant, and the solution is given implicitly as u = f[x - c(u)t]. The initial value problem requires specifying u(x, 0) = f(x). The solution is parametric, with x = ξ + c[f(ξ)]t and u = f(ξ). The derivatives of u with respect to t and x are derived, and substituting into the PDE confirms the solution. The validity of the solution depends on whether the solution breaks down. The derivatives become singular when t = -1/(dc/du f'(ξ)). If dc/du f'(ξ) > 0, characteristics do not intersect, but if dc/du f'(ξ) < 0, characteristics intersect, leading to a breakdown at t_b = 1/|dc/du f'(ξ)|_max. For dc/du > 0, the slope of the initial distribution steepens for f'(x) < 0 and flattens for f'(x) > 0. For dc/du < 0, the opposite occurs. An example shows that the solution becomes multivalued after t_b, making it physically invalid. In traffic flow, the density ρ(x,t) and speed v(x,t) are related by q = ρv (flux). The conservation equation is ∂ρ/∂t + ∂q/∂x = 0, with wave velocity c(ρ) = v + ρ ∂v/∂ρ. Since ∂v/∂ρ < 0, c < v, implying the wave propagates backward. For the Lincoln tunnel, data show q = aρ(1 - ρ/ρ_j), with a = 17.2 mph and ρ_j = 250 vpm. The wave velocity c(ρ) = a(1 - 2ρ/