Defogging

The defogging model is based on solving an additional scalar transport equation that represents the mass fraction of water vapor. A source/sink term for the scalar is considered for condensation/evaporation of the fog layer as well as the latent heat required for transition.

When there is a difference between the water vapor content at the fog layer and the cell next to this surface, the model calculates a rate of evaporation or condensation depending on the conditions. The assumptions are:

  • The vapor content in the air does not affect the thermal properties of the vapor-air mixture.
  • The water vapor mass is neglected with respect to the total mass in a cell.

The rate of mass transfer per unit surface [kg/s-m2] is:

Figure 1. EQUATION_DISPLAY
m˙=ρgβgCG  lnCG(1-Cs)CS(1-Cg)
(207)

where ρ g is density [kg/m3] .

Figure 2. EQUATION_DISPLAY
β g = C emp D V / L R e 0.8 S c 0.43 [ m/s ]
(208)

where C emp is the empirical constant useful for calibration (0.05-0.9) for windshield applications, L is a characteristic length (cubic root of cell volume next to fog layer boundary) [m], D V is diffusion of vapor in air [m2/s] , R e is the Reynolds number (which is based on the gas state and the characteristic length L ), and S c is the Schmidt number, and:

Figure 3. EQUATION_DISPLAY
C G = 0.622 + 0.378 C g
(209)
Figure 4. EQUATION_DISPLAY
C S = 0.622 + 0.378 C s
(210)

where C g is the actual concentration of vapor in air and C s is the saturation concentration of vapor in air.

For a given temperature, the saturation pressure is calculated from the following expression:

Figure 5. EQUATION_DISPLAY
ps(T)=611.85e17.502(T273.15)T32.25
(211)

The amount of evaporated/condensed vapor in one time-step is:

Figure 6. EQUATION_DISPLAY
dm˙=m˙dt
(212)

The thickness of the liquid film is updated as:

Figure 7. EQUATION_DISPLAY
s = s t - 1 + d m / ρ l
(213)