Hydrodynamic Gradient Expansion Diverges beyond Bjorken Flow

The gradient expansion is the fundamental organizing principle underlying relativistic hydrodynamics, yet understanding its convergence properties for general nonlinear flows has posed a major challenge. We introduce a simple method to address this question in a class of fluids modeled by Israel-Stewart-type relaxation equations. We apply it to (1+1)-dimensional flows and provide numerical evidence for factorially divergent gradient expansions. This generalizes results previously only obtained for (0+1)-dimensional comoving flows, notably Bjorken flow. We also demonstrate that the only known nontrivial case of a convergent hydrodynamic gradient expansion at the nonlinear level relies on Bjorken flow symmetries and becomes factorially divergent as soon as these are relaxed. Finally, we show that factorial divergence can be removed using a momentum space cutoff, which generalizes a result obtained earlier in the context of linear response.