A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modelling has been applied to propose a model for dilute 2PHF turbulent flows.Phase interaction terms with a clear physical meaning enter the equations and the formalism provides some guidelinesfor the avoidance of closure assumptions or the rational approximation of these terms.Continuous phase averagedcontinuity,momentum,turbulent kinetic energy and turbulence dissipation rate equations have been rigorously andsystematically obtained in a single step.These equations display a structure similar to that for single-Phase flows.It is also assumed that dispersed phase dynamics is well described by a probability density function (pdf) equationand Eulerian continuity,momentum and fluctuating kinetic energy equations for the dispersed Phase are deduced.An extension of the standard κ-εturbulence model for the continuous phase is used.A gradient transport modelis adopted for the A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modeling has been applied to propose a model for dilute 2PHF turbulent flows.Phase interaction terms with a clear physical meaning enter the equations and the formalism provide some guidelinesfor the avoidance of closure assumptions or the rational approximation of these terms.Continuous phase averagedcontinuity, momentum, turbulent kinetic energy and turbulence dissipation rate rate equations have been rigorously and systematically obtained in a single step. These equations display a structure similar to that for single-phase flows.It is also described that dispersed phase dynamics is well described by a probability density function (pdf) equation and Eulerian continuity, momentum and fluctuating kinetic energy equations for the dispersed Phase are deduced. Extension of the standard κ-εturbulence model for the continuous phase is used. A gradient transport model is adopted for the