Highlights • DSMC simulations of the Enceladus South polar plumes. • Parametrization of the outgassing flow parameters (number density, speed components, kinetic and rotational temperature) at 10 km above ground, where the flow becomes collisionless. • Parametric study of the outgassing flow parameters variations as a function of the vent exit conditions: vent radius, water mixture mass flow, water vapor/water ice mass ratio, water ice grain radius, water vapor and water ice exit speed, vent exit angle and flow temperature. • The DSMC results that have been run for this paper are available at an online repository. Abstract The field parameters – number density, velocity components and temperature – for the Enceladus geysers, and possibly similar jets from other bodies, such as Europa, Ceres or comets, are expensive to obtain using physically correct simulation methods, such as the Direct Simulation Monte Carlo (DSMC) approach. It would be very useful to be able to correctly reproduce all the different states that the flow is undergoing while expanding into vacuum, from a high density and collisional state near the surface to free-molecular and collisionless at high altitude without resorting to expensive DSMC simulations in every case. In this work, we consider a two-phase water vapor/grain mixture exiting a circular vent, assuming a uniform radius of the water ice grains, and study how the field parameters can be fitted at an altitude of 10 km, where the flow is collisionless. To do so, we define simple functional forms for each of these fitted parameters, and we study how their coefficients vary as a function of the vent exit parameters, i.e. the vent radius, the water mixture mass flow, the water vapor/water ice mass ratio, the water ice grain radius, the water vapor and water ice exit speed, the vent exit angle and flow temperature. We define polynomial approximations to model these variations. We show that all the vent parameters have nearly-independent influences on the radial profiles at 10 km, except for the water vapor and water ice exit speed, for which we considered cross-correlations. We finally show that the geyser field parameters can be reconstructed using our parametric study for variations of the vent parameters within the range considered here, and in a time frame of a few milliseconds. The results of the parametrizations presented in this study can now be used to propagate the geyser field parameters using computationally inexpensive free molecular/ballistic codes up to higher altitudes. The DSMC results that have been run for this paper are available at an online repository. [ABSTRACT FROM AUTHOR]