14.2.3. Time Schemes in Multiphase Flow

In many multiphase applications, the process can vary spatially as well as temporally. In order to accurately model multiphase flow, both higher-order spatial and time discretization schemes are necessary. In addition to the first-order time scheme in Ansys Fluent, the second-order time scheme is available in the Mixture and Eulerian multiphase models, and with the VOF Implicit Scheme.


Important:  The second-order time scheme cannot be used with the VOF Explicit Schemes.


The second-order time scheme has been adapted to all the transport equations, including mixture phase momentum equations, energy equations, species transport equations, turbulence models, phase volume fraction equations, the pressure correction equation, and the granular flow model. In multiphase flow, a general transport equation (similar to that of Equation 23–16) may be written as

(14–5)

Where is either a mixture (for the mixture model) or a phase variable, is the phase volume fraction (unity for the mixture equation), is the mixture phase density, is the mixture or phase velocity (depending on the equations), is the diffusion term, and is the source term.

As a fully implicit scheme, this second-order time-accurate scheme achieves its accuracy by using an Euler backward approximation in time (see Equation 23–24). The general transport equation, Equation 14–5 is discretized as

(14–6)

Equation 14–6 can be written in simpler form:

(14–7)

where

This scheme is easily implemented based on Ansys Fluent’s existing first-order Euler scheme. It is unconditionally stable, however, the negative coefficient at the time level , of the three-time level method, may produce oscillatory solutions if the time steps are large.

This problem can be eliminated if a bounded second-order scheme is introduced. See Bounded Second-Order Implicit Time Integration for more information.