This component makes it possible to model the dynamic behavior of the droplets of a fog. I.e. the size of the droplets and their speed of fall.
The fog appears when the amount of vapour per volume (the vapour density) present in a gas node exceeds its maximum admissible value (the saturation density of the vapour). Once this limit is reached, any additional gain in moisture (or reduction in vapour pressure) would imply a liquid phase in the form of spherical drops of liquid suspended in the air node. The main objective of this model is not to finely represents the behavior of the real fog but rather to aviod sensible fraction of liquid in the gas node.
It is assumed that:
First, the amount of liquid (droplets) derives from the difference between density of the water species and the saturation density of water in gas at current condition. Then the diameter of the droplet is computed from the equality between the amount of liquid and the product of the number of droplets and the volume of one droplet and the density of liquid (water) of the droplet.
Where:
d_sat
is the saturation density of the condensable species in the gas mixture p_sat
is the saturation pressure of the condensable species in the gas mixture R
is the perfect gas law constant MM
is the molar mass of the condensable species T
is the temperature of the gas mixture V
is the volume of the gas node d_drop
is the diameter of a droplet The dynamic behavior of the droplet derives from the fundammental dynamic principe where the balance of forces is the weight, the drag and the buoyancy. As droplets are non inertial, the sum of these 3 forces is zero. Therefore, the velocity of droplets derives:
Where:
Fdrag
is the visquous force for a sphere in a laminar dfluid
is the density of the gas mixture of the gas node. As the current model is modelling only moist, the ddroplet
is fixed at 1.2 kg.m-3μfluid
is the dynamic viscosity of the gas mixture of the gas nodeVel
is the fall velocity of dropletsP
is the weight of a dropletddroplet
is the density of the liquid that the droplets is made of. As the current model is modelling only droplet of water, the ddroplet
is fixed at 1000 kg.m-3The drops falling 'on the floor of the gas node' are removed from the system gas node (the control volume indeed). The flow lost is therefore:
The liquid surface is a difficult parameter to calculate in the sense that the vertical distribution of the drops is unknown. The mass flow rate is therefore computed from the volume ratio. First, imagine that there is no slipping between the droplet and the carrier phase (the surrounding gas). Their velocities are thus equal. The number density is known as it is an input parameter. Knowing the diameter of the droplets and their number densities, the volume of liquid per total volume is determined. Finally, the volume flow rate leaving the node is computed with the assumption of no slip condition.