In the latter half of the 20th century, physicists Belinski, Khalatnikov and Lifshitz (BKL) proposed a general ansatz for solutions to the Einstein equations possessing a (spacelike) singularity. They suggest that, near the singularity, the evolution of the spacetime geometry at different spatial points decouples and is well-approximated by a system of autonomous nonlinear ODEs, and further that general orbits of these ODEs resemble a (chaotic) cascade of heteroclinic orbits called “BKL bounces”. In this talk, we present recent work verifying the validity of BKL’s heuristics in a large class of symmetric, but spatially inhomogeneous, spacetimes which exhibit (up to one) BKL bounce on causal curves reaching the singularity. In particular, we prove AVTD behavior (i.e. decoupling) even in the presence of inhomogeneous BKL bounces. The proof uses nonlinear ODE analysis coupled to hyperbolic energy estimates, and one hopes our methods may be applied more generally.