Documentation for Modified Weber Number

[RachaelDMueller]

Review paper on modified weber number and/or ask Shihan the significance of a wind speed of 3.7 m/s to better understand this piece of ModuleOil_0D.F90 [4078]:

if (wind .GT. 3.7) then
             WCC=F_WCC(wind)

             !--------------------------
             do i = 1,ClassesNbr
            ! Calculate mean droplet diameter for each oil droplet interval
            DropletDiameterList_(i) = MinDropletDiameter + DeltaDiameter * (i-1) + 0.5 * DeltaDiameter

            call GetP_Star(CharacteristicDiameter,FilmThickness,DropletDiameterList_(i),P_Star)

           if (i==1) then
            Qd(i) = MassOil  * WCC / (1.62*wind*3.14/9.8)
           else
            Qd(i) = (MassOil   * WCC  / (1.62*wind*3.14/9.8) ) - Qd(i-1)
            end if
        end do
            call GetP_Star(CharacteristicDiameter,FilmThickness,MaxDropletDiameter,P_Star)
            QdTotal = MassOil  * WCC / (1.62*wind*3.14/9.8)
  else
      Qd=0.0
      QdTotal=0.0
      DropletDiameterList_(1:ClassesNbr)=0.0
  end if

[RachaelDMueller]

This code doesn’t appear to be called. It is within the subroutine ModifiedWeberDiaList, which is called from GetModifiedWeberDropletDiameter but GetModifiedWeberDropletDiameter isn’t actually called. My understanding is that we are using a modified weber approach to calculate particle droplet size for dispersion, but this is not the location in the code that defines our modified wave number calculation.