The aim of this paper is the derivation of a multiphase model of compressible fluids. Each fluid has a different average translational velocity, density, pressure, internal energy as well as the energies related to rotation and vibration. The main difficulty is the description of these various translational, rotational and vibrational motions in the context of a one-dimensional model. The second difficulty is the determination of closure relations for such a system: the ‘drag’ force between inviscid fluids, pressure relaxation rate, vibration and rotation creation rates, etc. The rotation creation rate is particularly important for turbulent flows with shock waves. In order to derive the one-dimensional multiphase model, two different approaches are used. The first one is based on the Hamilton principle. This method gives thermodynamically consistent equations with a clear mathematical structure, coupling the various motions: translation, rotation and vibration. However the relaxation effects have to be added phenomenologically. In order to achieve the closure of the system and its numerical resolution we use the second approach, in which the pure fluid equations are discretized at the microscopic level and then averaged. In this context, the flow is considered to be the annular flow of two turbulent fluids. We also derive the continuous limit of this model which provides explicit formulae for the closure laws. The structure of this system of partial differential equations is the same as the one obtained by the Hamilton principle. The final issue is to determine the rate of energy (or entropy) rotation. We assume that all the entropy creation related to the various relaxation effects after the passage of the shock wave is converted to rotational motion. The one-dimensional model is validated by comparing its predictions with averaged two-dimensional direct numerical results. The problem on which this model is tested is the interaction of a shock wave propagating in a heavy gas with a light gas bubble. The results obtained by the one-dimensional multiphase model are in a very good agreement with the two-dimensional averaged results.

No abstract

No abstract

No abstract

A new numerical method for solving the Serre-Green-Naghdi (SGN) equations describing dispersive waves on shallow water is proposed. From the mathematical point of view, the SGN equations are the Euler-Lagrange equations for a 'master' lagrangian submitted to a differential constraint which is the mass conservation law. One major numerical challenge in solving the SGN equations is the resolution of an elliptic problem at each time instant. This is the most time-consuming part of the numerical method. The idea is to replace the 'master' lagrangian by a one-parameter family of 'augmented' lagrangians, depending on a greater number of variables, for which the corresponding Euler-Lagrange equations are hyperbolic. In such an approach, the 'master' lagrangian is recovered by the augmented lagrangian in some limit (for example, when the corresponding parameter is large). The choice of such a family of augmented lagrangians is proposed and discussed. The corresponding hyperbolic system is numerically solved by a Godunov type method. Numerical solutions are compared with exact solutions to the SGN equations. It appears that the computational time in solving the hyperbolic system is much lower than in the case where the elliptic operator is inverted. The new method is applied, in particular, to the study of 'Favre waves' representing non-stationary undular bores produced after reflection of the fluid flow with a free surface at an immobile wall.

AbstractWe derive a mathematical model of shear flows of shallow water down an inclined plane. The non-dissipative part of the model is obtained by averaging the incompressible Euler equations over the fluid depth. The averaged equations are simplified in the case of weakly sheared flows. They are reminiscent of the compressible non-isentropic Euler equations where the flow enstrophy plays the role of entropy. Two types of enstrophies are distinguished: a small-scale enstrophy generated near the wall, and a large-scale enstrophy corresponding to the flow in the roller region near the free surface. The dissipation is then added in accordance with basic physical principles. The model is hyperbolic, the corresponding ‘sound velocity’ depends on the flow enstrophies. Periodic stationary solutions to this model describing roll waves were obtained. The solutions are in good agreement with the experimental profiles of roll waves measured in Brock’s experiments. In particular, the height of the vertical front of the waves, the shock thickness and the wave amplitude are well captured by the model.

The mathematical model of shear shallow water flows of uniform density is studied. This is a 2D hyperbolic non-conservative system of equations which is reminiscent of a generic Reynoldsaveraged model of barotropic turbulent flows. The model has three families of characteristics corresponding to the propagation of surface waves, shear waves and average flow (contact characteristics). The system is non-conservative : for six unknowns (the fluid depth, two components of the depth averaged horizontal velocity, and three independent components of the symmetric Reynolds stress tensor) one has only five conservation laws (conservation of mass, momentum, energy and mathematical 'entropy'). A splitting procedure for solving such a system is proposed allowing us to define a weak solution. Each split subsystem contains only one family of waves (either surface or shear waves) and contact characteristics. The accuracy of such an approach is tested on exact 2D solutions describing the flows where the velocity is linear with respect to the space variables, and 1D solutions. The capacity of the model to describe the full transition observed in the formation of roll waves : from uniform flow to one-dimensional roll waves, and, finally, to 2D transverse 'fingering' of roll wave profiles is shown.

No abstract