Two models of Langmuir turbulence, the nonlinear Schro¨dinger equation and the Zakharov equations, are solved numerically for an initial value problem in which the electric field evolves from an almost flat initial condition via the modulational instability and finally saturates into a set of solitons. The two models agree well with each other only when the initial dimensionless electric field has an amplitude less than unity. An analytic soliton gas model consisting of equal‐amplitude, randomly spaced, zero‐speed solitons is remarkably good at reproducing the time‐averaged Fourier spectra in both cases.