A mixed hyperbolic-parabolic, non conservative, Reynolds Stress Model (RSM), is studied.
It is based on an underlying set of Langevin equations, and allows to describe turbulent mixing, including transient demixing effects as well as incomplete mixing.
Its mathematical structure is analysed, and specific regimes, related to acoustic-like, Riemann-type, or self-similar solutions, are identified. A second-order accurate numerical scheme is proposed in arbitrary curvilinear geometry.
Its accuracy and convergence behaviour are tested by comparison with analytical solutions in the different regimes.
The numerical scheme can be generalized to multi-dimensional configurations, with potentially cylindrical symmetry, on unstructured meshes.