As the trajectory of a particle is computed, FLUENT keeps track of the heat, mass, and momentum gained or lost by the particle stream that follows that trajectory and these quantities can be incorporated in the subsequent continuous phase calculations. Thus, while the continuous phase always impacts the discrete phase, you can also incorporate the effect of the discrete phase trajectories on the continuum. This two-way coupling is accomplished by alternately solving the discrete and continuous phase equations until the solutions in both phases have stopped changing. This interphase exchange of heat, mass, and momentum from the particle to the continuous phase is depicted qualitatively in Figure 22.9.1.
The momentum transfer from the continuous phase to the discrete phase is computed in FLUENT by examining the change in momentum of a particle as it passes through each control volume in the FLUENT model. This momentum change is computed as
|=||viscosity of the fluid|
|=||density of the particle|
|=||diameter of the particle|
|Re||=||relative Reynolds number|
|=||velocity of the particle|
|=||velocity of the fluid|
|=||mass flow rate of the particles|
|=||other interaction forces|
This momentum exchange appears as a momentum sink in the continuous phase momentum balance in any subsequent calculations of the continuous phase flow field and can be reported by FLUENT as described in Section 22.16.
The heat transfer from the continuous phase to the discrete phase is computed in FLUENT by examining the change in thermal energy of a particle as it passes through each control volume in the FLUENT model. In the absence of a chemical reaction (i.e., for all particle laws except Law 5) the heat exchange is computed as
|=||mass of the particle on cell entry (kg)|
|=||mass of the particle on cell exit (kg)|
|=||heat capacity of the particle (J/kg-K)|
|=||heat of pyrolysis as volatiles are evolved (J/kg)|
|=||temperature of the particle on cell entry (K)|
|=||temperature of the particle on cell exit (K)|
|=||reference temperature for enthalpy (K)|
|=||latent heat at reference conditions (J/kg)|
The latent heat at the reference conditions for droplet particles is computed as the difference of the liquid and gas standard formation enthalpies, and can be related to the latent heat at the boiling point as follows:
|=||heat capacity of gas product species (J/kg-K)|
|=||boiling point temperature (K)|
|=||latent heat at the boiling point temperature (J/kg)|
For the volatile part of the combusting particles, some constraints are applied to ensure that the enthalpy source terms do not depend on the particle history. The formulation should be consistent with the mixing of two gas streams, one consisting of the fluid and the other consisting of the volatiles. Hence is derived by applying a correction to , which accounts for different heat capacities in the particle and gaseous phase:
|=||particle initial temperature (K)|
The mass transfer from the discrete phase to the continuous phase is computed in FLUENT by examining the change in mass of a particle as it passes through each control volume in the FLUENT model. The mass change is computed simply as
This mass exchange appears as a source of mass in the continuous phase continuity equation and as a source of a chemical species defined by you. The mass sources are included in any subsequent calculations of the continuous phase flow field and are reported by FLUENT as described in Section 22.16.
Under-Relaxation of the Interphase Exchange Terms
Note that the interphase exchange of momentum, heat, and mass is under-relaxed during the calculation, so that
where is the under-relaxation factor for particles/droplets that you can set in the Solution Controls panel. The default value for is 0.5. This value may be reduced to improve the stability of coupled calculations. Note that the value of does not influence the predictions obtained in the final converged solution.
Two options exist when updating the new particle source terms , and . The first option is to compute the new source terms and the particle source terms, , and , at the same time. The second option is to update the new source terms, , and , every flow iteration, while the particle source terms, , and , are calculated every Discrete Phase Model iteration. The latter option is recommended for transient flows, where the particles are updated once per flow time step.
Interphase Exchange During Stochastic Tracking
When stochastic tracking is performed, the interphase exchange terms, computed via Equations 22.9-1 to 22.9-8, are computed for each stochastic trajectory with the particle mass flow rate, , divided by the number of stochastic tracks computed. This implies that an equal mass flow of particles follows each stochastic trajectory.
Interphase Exchange During Cloud Tracking
When the particle cloud model is used, the interphase exchange terms are computed via Equations 22.9-1 to 22.9-8 based on ensemble-averaged flow properties in the particle cloud. The exchange terms are then distributed to all the cells in the cloud based on the weighting factor defined in Equation 22.2-55.