# 3.3.2. Particle Transport Model

In this section, the equations describing the evolution of parcels of particles are presented. Models are presented for the particle motion and for heat and mass transfer (ie. evaporation and combustion).

## 3.3.2.1. Particle Acceleration and Trajectories

Particles with densities much greater than that of a the fluid phase (solid or liquid particles in gaseous flows) are primarily affected by drag forces and body forces. In this limit where , the particle acceleration is written [101]:

(3.622)

where and are the component of the particle and gas velocities, respectively, is the vector magnitude of the velocity differences, and are the particle and gas densities, and is the component of the acceleration due to body forces. The particle diameter is and this should be considered to be the equivalent particle diameter corresponding to a spherical particle. The effects of non-sphericity on the acceleration can be accounted for through the drag coefficient, . The gas-velocity to be employed in (3.622) is taken from the Eulerian solution of the continuum field. For turbulent flows, the sum of the mean (resolved) velocity and a perturbation to that mean, accounting for the turbulent fluctuations, both contribute to the gas velocity. The effects of turbulent fluctuations are described in the following section, Particle Dispersion and Turbulence.

For a spherical particle, the drag coefficient is modeled using standard drag coefficient relations:

(3.623)

The particle Reynolds number, , is based on the slip velocity and the particle diameter

(3.624)

Yuen and Chen [116] recommend evaluating the viscosity in (3.624) based on the weighted average of the properties at the gas-phase side of the particle surface (weighted by two-thirds) and the gas-phase properties far from the particle (weighted by one-third), the so called ‘1/3 rule.’ So for properties at the surface and far field identified with a subscripted and , respectively, the viscosity would be

(3.625)

A similar relationship is suggested for other transport coefficients (conductivity, diffusivity, etc.). Note that additional forces are relevant for particles with densities nearer to or less than the continuum phase (bubbles) and for particles with high Reynolds numbers. A comprehensive overview of the force on particles is available from Maxey and Riley [107]. Equation (3.622) can be linearized and written

(3.626)

(3.627)

where is the particle velocity response time. In the small Reynolds number limit where the drag coefficient is inversely proportional to the slip velocity, this linearization is particularly relevant. Given the particle acceleration from Eq. (3.622), the particle trajectories can be determined by integrating the simple ODE

(3.628)

Since all the particles within a parcel are the same size, all parcel trajectories are determined by Eqn. (3.628), subject to turbulence effects described below in Particle Dispersion and Turbulence.

## 3.3.2.2. Particle Dispersion and Turbulence

In CFD modeling of turbulent flows, the full velocity spectrum is generally not resolved. Instead, certain velocity fluctuations are modeled, being represented through the turbulent kinetic energy, , which is half the sum of the squares of the velocity fluctuations. These velocity fluctuations tend to introduce random fluctuations in the particle velocities that result in real particles being dispersed relative to the mean continuum flow [117, 118]. For Lagrangian particle methods, this phenomenon is modeled in two ways: by perturbing the velocity of parcels of particles and by affecting the spatial extent of the parcel itself.

To account for the effects of the velocity fluctuations on the parcels or particles, the random walk model of Gosman and Ioannides [111], as modified by Shuen et al. [112], is employed. In this approach, the gas velocity employed in equations (3.622), (3.624), (3.626), and (3.627) is the sum of the mean gas velocity, , and a fluctuating velocity that is sampled from a normal (Gaussian) velocity distribution with a standard deviation given by . For LES, here would be replaced with the subgrid kinetic energy. Sampling from the inverted cumulative distribution function with a random number uniformly distributed between zero and unity, , gives the appropriate fluctuating velocity. The total and fluctuating gas velocities are then

(3.629)

(3.630)

where is the inverse error function. The time during which a given velocity fluctuation affects a given parcel is determined by the expected time that the particle takes to cross the eddy inducing the given velocity fluctuation. Small particles will tend to stay in an eddy for the duration of the eddy lifetime,

(3.631)

where is the turbulent energy dissipation rate and . Larger particles will have sufficient slip velocity to cross the eddy. The eddy-crossing time is estimated as

(3.632)

where the eddy length scale is . The eddy interaction time, that is the time over which a given velocity perturbation affect in Equations (3.629) and (3.630), is the minimum of and . The particles comprising a parcel are presumed to be distributed about the center of the parcel, tracked by Eq. (3.628), in a normal (Gaussian) manner with the standard deviation in each direction given by , where is either , , or . This distribution is written

(3.633)

where is the total number of particles in the parcel and is the center of the parcel. The spatial extent of the parcel is thus determined by . This term is affected by unresolved turbulent fluctuations that act differently on the particles across the parcel. Following Zhou and Yao [113], the spatial extent of the parcel is the mean square displacement over time

(3.634)

where satisfies the ODE

(3.635)

and is the time over which acts. Equation (3.635) is obtained by subtracting the instantaneous particle equations from the mean particle equations. Note that large particles with large are not dispersed appreciably.

## 3.3.2.3. Mass and Energy Exchange between Particles and the Gas Phase

Particles may exchange mass and energy with the gas phase according to conservation principles across the interface. The physics of mass and energy transfer are described in detail here because the anticipated applications include a wide range of physical phenomena that have not been described together. Included in the phenomena of interest are metal and hydrocarbon particle combustion as well as particle condensation and evaporation. Each of these is anticipated to be strongly energetic in the sense that the product of the evaporation rate with the enthalpy of evaporation and combustion is expected to be substantial. Further, for evaporation and condensation applications, it is expected that vapor pressures will range over a sufficiently wide range that the transition from evaporation to condensation should be correctly described as this often limits evaporation and condensation rates. For metal combustion, the associated temperatures are sufficiently high so that radiative heat transfer must be considered.

The theory for droplet vaporization and combustion has generally been developed based on a large number of simplifications [95, 96] including spherical-symmetry, unity Lewis numbers, droplet surfaces at the boiling temperature, and infinitely fast conduction through the droplet. Recent work provides guidance as to how these assumptions can be relaxed [119].

### 3.3.2.3.1. Theory for spherically symmetric flow

In general, the theory of heat and mass transfer is developed for spherically symmetric systems and then corrected to account for increased transfer associated with advection and asymmetry. In this section, relations are developed based on the assumption of spherical symmetry and corrections is provided in Extension to multiple oxidizers.

(3.636)

can be integrated to give

(3.637)

where is the rate of mass evaporation from the particle, is the Stefan velocity, directed normally away from the particle, and is the local vapor density considering both the particle vapor and the continuum gas concentrations. A coordinate transformation to the variable

(3.638)

or

(3.639)

which represents the ratio of the Stefan velocity to the thermal diffusion velocity, greatly simplifying the species and energy conservation equations. In Equations (3.638) and (3.639), is the thermal conductivity of the vapor, is its specific heat at constant pressure, and is the evaporating species diffusion coefficient. The subscript 0 on the evaporation rate indicates that this evaporation corresponds to that for the spherically-symmetric case; corrections relating the evaporation rate for the spherically-symmetric case to that with finite-slip velocities are provided in Extension to multiple oxidizers. In the spherically-symmetric case, conservation equations for conserved scalars, , can be written

(3.640)

for which an analytic solution

(3.641)

is readily obtained. In non-reacting, non-radiating flows, any mass fraction or temperature can be a conserved scalar. Other conserved scalars are provided below. The appropriate from of (3.638) and (3.639) to be used depends on the variable represented by (3.640). Since the diffusion coefficients appearing in (3.639) are most important near the particle surface, the diffusion coefficient to be employed is that most relevant at the surface. The choice will be clearly identified below. The application of the Dirichlet boundary conditions at both the surface and far from the droplet and the application of a Neumann boundary condition at the surface relate the boundary conditions for the conserved scalar to through

(3.642)

where the subscript indicates quantities evaluated at the droplet surface, the so-called film state. For reacting flows we will consider two such conserved scalars in the present work. The general configuration considered is a fuel droplet reacting in an oxidizing atmosphere; in the event that, for example, an oxidizer droplet is reacting in a fuel atmosphere, then the “oxidizer” and “fuel” described would be switched. Also, if there is no reaction, the results generalize to droplet evaporation where the “fuel” is the evaporating component. Allowing variable properties, the temperature oxidizer coupling function is

(3.643)

where is the molecular weight and is the stoichiometric coefficient of species . The standard enthalpy of combustion for fuel and oxidizer, per unit mass of fuel evaluated at the film temperature, , is

(3.644)

where the last summation is taken over the produces of reaction, . Note that when there is no combustion, the second term in Eqn. (3.644) is ignored. For , the thermal diffusivity is the relevant diffusivity so that Eqn. (3.638) is employed in conjunction with .

The constants and in (3.641) are evaluated using the boundary conditions for the temperature and oxidizer at the droplet surface and in the far field. The boundary conditions employed for temperature are that the heat flux into the particle is balanced by the sum of the enthalpy of vaporization, the heating of the particle and any radiative losses

(3.645)

Here, the enthalpy of vaporization is , the particle specific heat is , and the radiative heat loss over the droplet surface is

(3.646)

where is the particle absorptivity, is the Stefan-Boltzmann constant, and is the incident radiation. The incident radiation is the radiation intensity integrated over all directions, that is, the entire steradian solid angle. If the radiation intensity is , then where is the differential solid angle. If the particle is not opaque, the particle absorptivity will be a function of the particle size [120]. Note that, in Eqn (3.645), the evaporation rate appears without the subscript 0, indicating that this is the evaporation rate corrected for finite-slip velocities as prescribed in (3.661) below. This is appropriate because the heat and mass flux to the surface are both increased by the relative droplet motion, while the other terms in (3.645) are not affected by that. The oxidizer is presumed to not be absorbed by the surface so that a no-flux boundary condition is employed

(3.647)

Because the thermal diffusivity is used in definition for , the ratio of the thermal to mass diffusivity in the form of the oxidizer Lewis number appears in (3.647). Taking the derivative of Eqn. (3.643) evaluated at the droplet surface and using (3.645) and (3.647), we obtain

(3.648)

that will appear in the denominator of Eqn. (3.642). Also the temperature and oxidizer mass fraction must approach their far field values at large radii. Additional assumptions are required to identify the temperature and the oxidizer mass fraction at the surface; these will be discussed later. Applying these boundary conditions to (3.641) for provides an expression for the rate of evaporation in terms of evaluated at the surface of the particle

(3.649)

where the Spalding transfer number associated with is

(3.650)

Here the subscripts and indicate the states at the particle surface (on the gas-phase side of the interface) and the ambient far-field environment, respectively. The denominator in (3.650) represents a variety of potential sinks for the enthalpy at the droplet surface. These sinks include, in the order in which they are written, the enthalpy associated with vaporizing the particle, the radiative losses from the surface, the enthalpy conducted into the particle, the enthalpy flux to the gas-phase due to variable specific heats, and the enthalpy of combustion lost as a consequence of oxidizer leakage. Note that the radiative absorption and emission from the particle surface, , are included here in the surface boundary condition, but the radiative losses from a flame around the droplet must be accounted for by modifying to provide an effective heat release, the heat release decremented by the near flame radiative losses. As far as the droplet is concerned, these radiative losses are far enough away to not affect the film state except to the extent that radiative flux to the surface is affected.

The above expressions comprise a relatively complete definition of the physics of particle evaporation, combustion, and interaction with a radiative field accounting for variable thermophysical properties. These expressions are simplified by making the assumption that no oxidizer penetrates a flame to reach the surface, , and by setting leading to

(3.651)

where the denominator is referred to as the effective enthalpy

(3.652)

Equation (3.651) is employed in the numerical models. A coupling function for fuel and oxidizer is similar written

(3.653)

For a species that is evaporating, the flux boundary condition is

(3.654)

and the fuel-oxidizer Spalding mass transfer number is

(3.655)

where has been assumed as in Eq. (3.651). Note that only one of or will be non-zero based on the current assumption of zero leakage through flames; if combustion is occurring then it will be that is zero. The diffusion coefficient appropriate for the fuel-oxidizer coupling function is that for the fuel so that the second of Eqn (3.638) and (3.639) is used with the diffusion coefficient specifically that of the fuel, and the equivalent of Eqn. (3.649) for the fuel-oxidizer system

(3.656)

### 3.3.2.3.2. Extension to multiple oxidizers

Multiple oxidizers is discussed in detail later, in Multiple Oxidizers.

### 3.3.2.3.3. Correlations for finite slip velocities

The above relationships for the heat and mass transfer are derived for a spherically symmetric field around the droplet and are valid for droplets with zero slip velocity in the absence of buoyancy. Empirical correlations are available in terms of the Nusselt and Sherwood numbers parameterized by Reynolds, Schmidt, and Prandtl numbers to describe the effect of finite slip velocities in modifying the heat and mass transfer by reducing the boundary layer thickness. The Schmidt and Prandtl numbers are defined:

(3.657)

(3.658)

These quantities with the subscript represent appropriate averages for transport properties in the gas phase boundary layer around the particle.These gas-phase quantities are evaluated using an appropriate averaging process that will generally be analogous to Eqn. (3.625). The alternative is to tabulate these quantities using Equations (3.665) and (3.666). The Nusselt number describes a dimensionless heat transfer rate to the droplet for a given difference between the ambient and surface temperature,

(3.659)

The correction to the evaporation rate employed in the present work is based on measurements by Ranz and Marshall (1952)

(3.660)

This can be introduced into (3.650), and the evaporation rate can be written

(3.661)

The Nusselt number for zero slip velocity, , is included in the relations of the previous section, specifically in Eq. (3.650). A similar correlation,

(3.662)

can be used for the Sherwood number,

(3.663)

which is the dimensionless mass transfer coefficient, so that the evaporation rate for finite slip velocities can also be written

(3.664)

where should be taken from Eqn. (3.656).

### 3.3.2.3.4. Evaporation rates and effective diffusivities

The evaporation rate is seen in Eqn (3.649) and (3.656) to depend linearly on an area-weighted gas-phase diffusivity. Since accurate evaluation of this weighting is not feasible for the particle transport model, effect diffusivities,

(3.665)

(3.666)

are defined. The one-third rule defined in Eq. (3.625) provides a rough guideline for evaluating these diffusion coefficients in some cases. The effective thermal diffusion coefficient is combined with Equations (3.649) and (3.661) to give an evaporating rate based on thermal diffusion

(3.667)

The strong dependence of the evaporating rate on the diffusion coefficient, coupled with the fact that the diffusion coefficients depend strongly on variations in the compositions and temperature of the gases around the droplet, mean that the reasonable but judicious choice of diffusion coefficients can often match observed experimental measurements. Similarly, the mass transfer driven evaporation rate can be written with the use of Equations (3.656), (3.662), and (3.664)

(3.668)

It is necessary that the evaporation rate predicted by Eqn. (3.668) be equal to that predicted by (3.667). These equations show that thermal and mass diffusion vary both through their diffusion coefficients and through difference in boundary layer thickness attributable to finite slip velocities. Both these effects are combined in an effective Lewis number, , to give

(3.669)

Equations (3.667), (3.668), and (3.669) are used in the computational model for the evaporation rate. Note that energy and mass conservation must vie the same evaporation rate; this requirement determines the film state as described in the following section, Closure of film state with effective heat transfer coefficient.

### 3.3.2.3.5. Closure of film state with effective heat transfer coefficient

The system described by Equations (3.667) and (3.668) is closed with two additional assumptions. First, the film conditions, and , are related through the Clausius-Clapeyron relationship

(3.670)

where the partial pressure gives the mole fraction through that can subsequently be converted to the mass fraction with the relationship . As provided by Lefebvre [121], such a relationship is

(3.671)

for and zero otherwise. If the critical point temperature is not provided, the code sets to a very large value, essentially making independent of temperature. Second, the droplet heating is related to the difference between the film temperature and the droplet temperature by assuming an effective internal heat transfer coefficient in the form of an internal Nusselt number, , for the particle so that

(3.672)

This internal Nusselt number, which is different from that for the external heat transfer indicated in Eqn. (3.661), can be estimated based on the results of numerical studies where the internal droplet was resolved [119, 122]. There the Nusselt numbers for no circulation and for rapid circulation were identified as 6.58 and 17.9, respectively, and a transition region was also identified based on the liquid Peclet number

(3.673)

where

(3.674)

which is based on the maximum surface velocity, . This transition was empirically fitted [119] to

(3.675)

In [119], both the evaporation rate and the surface temperature were reproduces using Eq. (3.675) with little error compared to simulations incorporating a detailed internal droplet convection model. Naturally, the liquid Peclet number should be zero if the particle is below the freezing temperature of the particle constituent.

The evaporation rates indicated in Equations (3.667) and (3.668) must be equal, subject to the constraints of Equations (3.670) and (3.672) based on the closure approximation employed in this section. Equating Equations (3.667) and (3.668) leads to a nonlinear equation for the surface temperature that is to be solved. This is readily solved using Newton’s method as described here. Newton’s method is an iterative root-finding method written in the form

(3.676)

where is the equation for which the root will be found and the superscript refers to the iteration number. Equating Equations (3.667) and (3.668) and using (3.669) gives

(3.677)

Differentiation gives

(3.678)

(3.679)

(3.680)

(3.681)

(3.682)

In writing these expressions for , it is assumed that and .

Because of the strong nonlinearities in , care must be taken in providing initial conditions to solve these relations. This is conducted through a two-stage procedure. First, the boiling point temperature is identified, then a temperature just under the boiling temperature is used as an initial guess for the iterative solution that determines the film temperature. This procedure arises from a consideration of the shape of the . For realistic temperatures, those for which , both and , are monotonically strongly increasing in magnitude. An initial guess with a temperature that is too low (less than ) results in a prediction, with the Newton’s method, of a very high temperature on the successive iteration due to the small magnitude of the derivative, , for a small T. Typically, this second iteration will result in a temperature for which that leads to negative values of and the iteration diverges into non-physical regimes. However, an initial guess that is within the physically reasonable regime () and yet above the final will reliably converge. Therefore, the initial guess of is used as an initial guess in determining the film temperature. This method has been tested for a wide range of conditions and appears robust except for those scenarios where the denominator of takes on negative values. (Since can be less than and this can result in the denominator of taking on negative values, in which case the iteration may fail.)

Because of the cost and potential stability issues in this method, the fast conduction limit described in Closure for surface state assuming fast conduction is used in the code instead.

### 3.3.2.3.6. Closure for surface state assuming fast conduction

The surface state described in the previous section is the most physically realistic state that can be obtained without solving a differential equation for the heat transfer with the droplet. However, determining this state involves the iterative solution of a system of nonlinear equations. In the previous section a robust method of solving these equations is provided, but a simpler approximation may provide suitable results under certain conditions. This simpler approach is to assume that the heat transfer within the droplet is fast relative to the heating of the droplet. This is equivalent to taking the zero Biot number limit, , in Eqn. (3.672), in which case . With the film state determined by the droplet temperature, the film mass fraction is directly obtained from Equations (3.670) and (3.671), and the mass transfer number is obtained from Eq. (3.655) with the evaporation rate following from Eqn. (3.668). The thermal transfer number is obtained from the mass transfer number by equating Equations (3.668) and (3.667), and the droplet heating is obtained y solving for the enthalpy change of the particle in Eq. (3.650) to obtain

(3.683)

Sirignano [123] and coworkers have demonstrated that for many conditions, this particular limit is a poor approximation for at least some of the particle lifetime. In the present models, this limit is employed in two situations: (a) if the particle temperature is within 1% of the wet bulb temperature, and (b) if the particle temperature exceeds the wet bulb temperature. As the droplet temperature approaches the wet bulb temperature, employing this limit is inconsequential. In the latter case, the rate of droplet cooling is likely to be over predicted, but this is a scenario for which the convergence of the film state otherwise is not guaranteed. In the interest of creating a more robust model, and because this particular situation is not anticipated to be predominant, we employ this simpler limit.

## 3.3.2.4. Conserved scalars and transfer numbers for various applications

In Correlations for finite slip velocities, expressions are provided for the particle evaporation rate, (3.668) and (3.667), using a model system of fuel evaporating from the particle and reactive with an oxidizer that diffuses from the ambient gas. Of significance in those expressions are the transfer numbers defined in Theory for spherically symmetric flow in Equations (3.651) and (3.650), and these in turn are based upon conserved scalars defined in Equations (3.643) and (3.653). In this appendix, expressions for alternate transfer numbers and conserved scalars are provided for two additional systems: a simpler system in which evaporating or condensation occurs without any reaction in the boundary layer (ie. water droplet evaporation or condensation) and a more complicated system in which multiple oxidizers are involved in the oxidation of the evaporated fuel (relevant to metal oxidation).

### 3.3.2.4.1. Simple evaporation and condensation

When species evaporate or condense but do not otherwise react in the boundary layer surrounding the particle, the species mass fraction of the evaporating/condensing species itself is a conserved scalar. In this case, the mass transfer number corresponding to Eq. (3.650) is simply

(3.684)

and the heat transfer number corresponding to (3.651) is

(3.685)

### 3.3.2.4.2. Multiple Oxidizers

For metal oxidation, it is possible to have multiple oxidizers that simultaneously react with the evaporating metal. These can be expressed with a series of parallel single-oxidizer reactions of the form . Useful conserved scalars that can be formed with this set of reactions include

(3.686)

(3.687)

where the enthalpy of reaction of the oxidizer, , with is given by . These are analogous to the conserved scalars defined in Equations (3.643) and (3.653). Using these conserved scalars, a mass transfer number analogous to Eqn. (3.650) is found to be

(3.688)

and the heat transfer number is

(3.689)

These can be employed in Eqn. (3.667) and (3.668) to provide expressions for the particle mass burning rate as a function of the various oxidizer mass fractions far from the particle.

## 3.3.2.5. Energy Exchange between Particles and the Gas Phase Without Mass Transfer

In this section, a simpler scenario is considered where evaporation from and condensation onto particles is presumed to be negligible. For example, metal particles in dry air at ambient temperatures are unlikely to participate in evaporation or condensation. Models are presented in this section to treat these scenarios.

The models described in the previous sections are ill posed to solve these problems because the formulation is based on a balance between diffusion and the Stefan convective velocity, which is proportional to the negligible . The models in this section are used for the heated particles model.

In the event that there is no mass transfer, the particle heating rate is determined based on the balance between conductive and radiative transfer. The conductive heat transfer can be expressed using an effective heat transfer coefficient, thereby taking advantage of available Nusselt number correlations indicated in Simple evaporation and condensation. Then the relationship for the droplet heating is

(3.690)

For the particle with no mass transfer, the closure of the film (or surface) temperature is obtained by equating the external heat flux, described by the right hand side of Eqn. (3.690), with the internal heat flux, as indicated on the right hand side of Eqn. (3.672) to obtain the following quartic constraint, which is solved for

(3.691)

## 3.3.2.6. Further Notes on Radiative Heat Transfer

The expression employed for radiative droplet heating in Eqn. (3.646) is appropriate for relatively large and opaque particles, referred to as the geometric optics limit. For different particles, the expression for the absorptivity will change, but otherwise the expressions remain appropriate. To identify the appropriate absorptivity, it is useful to compare absorption coefficients for particle clouds. The absorption coefficient, , in units of inverse length, is an effective cross-sectional area per volume. For large particles with a condense phase absorptivity give by , the absorption coefficient is the summation of the cross-sectional areas time the absorptivity

(3.692)

For comparison, in the small-particle or Rayleigh limit, the absorption coefficient is proportional to the volume fraction with the proportionality coming from the complex index of refraction. For intermediate particles where the particle optical depth is comparable to the particle radius, intermediate limits are appropriate, and the absorption can range from being proportional to the particle area to being proportional to the particle volume. The appropriate absorption coefficient, and from it the particle emissivity, cases can be determined as described in the available texts [120, 124].

## 3.3.2.7. Input Parameters for Particle Evolution

A large number of parameters are required to specify the evolution of the particles. The parameters provided in the input file are described in table Input parameters related to the particle evolution provided through the input file. In addition, table Variables passed from the gas-phase Eulerian solver required to evolve the particles specifies those variables that must be obtained from the gas-phase continuum flow. Because parcels of particles have finite extent and may span several control volumes, it is sometimes necessary to interpolate values of these gas-phase variables from several control volumes.

Table 3.12 Input parameters related to the particle evolution provided through the input file
Variable | Input name | Units | Description
--- | --- | --- | ---
| INJP_DENp |  | Particle density
| INJP_Tp |  | Particle temperature
| INJP_MWp |  | Molecular weight of fuel or particle component
| INJP_Clp |  | Particle specific heat
| INJP_Emp | none | Particle absorptivity,emissivity
| INJP_Prf | none | Film Prandtl number
| INJP_Scf | none | Film Schmidt number
| INJP_Prefvapp |  | Reference pressure for vaporization for particle
| INJP_Trefvapp |  | Reference temperature for vaporization for particle
| INJP_Hvaprefp |  | Reference enthalpy of vaporization for particle
| INJP_T1cp |  | Critical temperature for particle
| INJP_T1fp |  | Freezing temperature for particle
| INJP_STp |  | Particle surface tension
| INJP_mup |  | Particle viscosity
| INJP_Prl | none | Particle Prandtl number
|  |  | Enthalpy of combustion for vapor species at
|  | mol/ mol Fuel | Stoichiometric coefficients (relative to evaporating species)

Table 3.13 Variables passed from the gas-phase Eulerian solver required to evolve the particles
Variable | Units | Description
--- | --- | ---
|  | Gas-phase density
|  | Pressure
|  | Mean gas velocity
|  | Turbulent kinetic energy
|  | Turbulent kinetic energy dissipation rate
|  | Gas temperature
|  | Gas-phase mean molecular weight
| N/A | Mass fraction of fuel in gas phase
| N/A | Mass fraction of oxidizer in gas phase
|  | Gas-phase specific heat
|  | Gas-phase viscosity
|  | Incident radiation