Compressible fluids like air can develop sharp discontinuities in their behavior, similar to shock waves. Mathematicians studying these phenomena work with the 3D compressible Euler equations, which govern how gases flow when you account for changes in pressure, density, and spinning motion. A fundamental question is: given smooth initial conditions, how far can a smooth solution be extended before it breaks down, and what does that breakdown look like geometrically? Previous work by the same authors showed that solutions can develop a specific type of singularity, where the gradient of the fluid variables blows up along a surface, and that this singular surface originates from a special lower-dimensional curve in space called the crease.
This paper focuses on what happens at and near the crease from a different direction. Alongside the singular surface where the solution breaks down, the authors prove that another surface, called a Cauchy horizon, emerges from the crease and travels in a perpendicular direction. On the Cauchy horizon side, the fluid solution remains perfectly smooth, even though it is influenced by the nearby singularity at the crease. Together, the crease, the singular surface, and the Cauchy horizon form the complete boundary of the largest possible region in which the fluid equations can be solved smoothly and uniquely, a region mathematicians call the maximal globally hyperbolic development. Understanding this boundary is critical because without controlling its full structure, you cannot guarantee that the smooth solution is unique.
The broader significance is that this is the first rigorous mathematical description of a Cauchy horizon forming stably in three dimensions without imposing artificial symmetry on the problem. Real fluids have varying entropy and vorticity, quantities that are absent or trivial in simpler model problems, and these features create serious mathematical difficulties. The authors develop new geometric tools, including specially adapted coordinate systems that track the shape of the crease as it evolves, to handle these complications. The results establish precise structural properties along the boundary of the solution domain, which are the key ingredients needed to eventually prove uniqueness of the maximal smooth solution in this physically realistic setting.