High fidelity modeling and simulation of moderately dense sprays at relatively low cost is still a major challenge for many applications. For that purpose, we introduce a new multi-fluid model based on a two-size moment formalism in sections, which are size intervals of discretization. It is derived from a Boltzmann type equation taking into account drag, evaporation and coalescence, which are representative of the complex terms that arise in multi-physics environments. The closure of the model comes from a reconstruction of the distribution. A piecewise affine reconstruction in size is thoroughly analyzed in terms of stability and accuracy, a key point for a high-fidelity and reliable description of the spray. Robust and accurate numerical methods are then developed, ensuring the realizability of the moments. The model and method are proven to describe the spray with a high accuracy in size and size-conditioned variables, resorting to a lower number of sections compared to one size moment methods. Moreover, robustness is ensured with efficient and tractable algorithms despite the numerous couplings and various algebra thanks to a tailored overall strategy. This strategy is successfully tested on a difficult 2D unsteady case, which proves the efficiency of the modeling and numerical choices.

We describe possible solutions for a stationary flow of two superposed fluids between two close surfaces in relative motion. Physically, this study is within the lubrication framework, in which it is of interest to predict the relative positions of the lubricant and the air in the device. Mathematically, we observe that this problem corresponds to finding the interface between the two fluids, and we prove that this interface can be viewed as a square root of a polynomial of degree at most 6. We solve this equation using an original method. First, we check that our results are consistent with previous work. Next, we use this solution to answer some physically relevant questions related to the lubrication setting. For instance, we obtain theoretical and numerical results, which can predict the occurrence of a full film with respect to physical parameters (fluxes, shear velocity, viscosities). In particular, we present a figure giving the number of stationary solutions depending on the physical parameters. Moreover, we give some indications for a better understanding of the multi-fluid case.