  The prediction of the drop size distribution (DSD) resulting from liquid
atomization is key to the optimization of multi-phase flows, from gas-turbine
propulsion, through agriculture, to healthcare. Obtaining high-fidelity data of
liquid atomization, either experimentally or numerically, is expensive, which
makes the exploration of the design space difficult. First, to tackle these
challenges, we propose a framework to predict the DSD of a liquid spray based
on data as a function of the spray angle, the Reynolds number, and the Weber
number. Second, to guide the design of liquid atomizers, the model accurately
predicts the volume of fluid contained in drops of specific sizes whilst
providing uncertainty estimation. To do so, we propose a Gaussian process
regression (GPR) model, which infers the DSD and its uncertainty form the
knowledge of its integrals, and of its first moment, i.e., the mean drop
diameter. Third, we deploy multiple GPR models to estimate these quantities at
arbitrary points of the design space from data obtained from a large number of
numerical simulations of a flat fan spray. The kernel used for reconstructing
the DSD incorporates prior physical knowledge, which enables the prediction of
sharply peaked and heavy-tailed distributions. Fourth, we compare our method
with a benchmark approach, which estimates the DSD by interpolating the
frequency polygon of the binned drops with a GPR. We show that our integral
approach is significantly more accurate, especially in the tail of the
distribution (i.e., large, rare drops), and it reduces the bias of the density
estimator by up to ten times. Finally, we discuss physical aspects of the
model's predictions and interpret them against experimental results from the
literature. This work opens opportunities for modelling drop size distribution
in multiphase flows from data.
