add isotopic ventilation ratio to physics (f_heavy to f_light) (#1532)

Co-authored-by: Sylwester Arabas <[email protected]>

Author: AgnieszkaZaba

PySDM/formulae.py

@@ -23,7 +23,7 @@
 from PySDM.dynamics.terminal_velocity.gunn_and_kinzer import TpDependent
 
 
-class Formulae:  # pylint: disable=too-few-public-methods,too-many-instance-attributes
+class Formulae:  # pylint: disable=too-few-public-methods,too-many-instance-attributes,too-many-statements
     def __init__(  # pylint: disable=too-many-locals
         self,
         *,
@@ -55,6 +55,7 @@ def __init__(  # pylint: disable=too-many-locals
         isotope_diffusivity_ratios: str = "Null",
         isotope_relaxation_timescale: str = "Null",
         isotope_temperature_inference: str = "Null",
+        isotope_ventilation_ratio: str = "Null",
         optical_albedo: str = "Null",
         optical_depth: str = "Null",
         particle_shape_and_density: str = "LiquidSpheres",
@@ -97,6 +98,7 @@ def __init__(  # pylint: disable=too-many-locals
         self.isotope_diffusivity_ratios = isotope_diffusivity_ratios
         self.isotope_relaxation_timescale = isotope_relaxation_timescale
         self.isotope_temperature_inference = isotope_temperature_inference
+        self.isotope_ventilation_ratio = isotope_ventilation_ratio
         self.particle_shape_and_density = particle_shape_and_density
         self.air_dynamic_viscosity = air_dynamic_viscosity
         self.terminal_velocity = terminal_velocity

PySDM/physics/__init__.py

@@ -37,6 +37,7 @@
     isotope_diffusivity_ratios,
     isotope_relaxation_timescale,
     isotope_temperature_inference,
+    isotope_ventilation_ratio,
     latent_heat_vapourisation,
     latent_heat_sublimation,
     optical_albedo,

PySDM/physics/constants_defaults.py

@@ -519,6 +519,8 @@
 PRUPPACHER_RASMUSSEN_1979_XTHRES = 1.4 * si.dimensionless
 """
 [Pruppacher & Rasmussen 1979](https://doi.org/10.1175/1520-0469%281979%29036%3C1255:AWTIOT%3E2.0.CO;2)
+also in
+[Beard & Pruppacher 1971](https://doi.org/10.1175/1520-0469%281971%29028%3C1455:AWTIOT%3E2.0.CO;2)
 """  # pylint: disable=line-too-long
 PRUPPACHER_RASMUSSEN_1979_CONSTSMALL = 1.0 * si.dimensionless
 """ 〃 """

PySDM/physics/isotope_diffusivity_ratios/grahams_law.py

@@ -8,7 +8,7 @@ def __init__(self, _):
         pass
 
     @staticmethod
-    def ratio_2H(const, temperature):  # pylint: disable=unused-argument
+    def ratio_2H_heavy_to_light(const, temperature):  # pylint: disable=unused-argument
         return (
             (2 * const.M_1H + const.M_16O) / (const.M_2H + const.M_1H + const.M_16O)
         ) ** const.ONE_HALF

PySDM/physics/isotope_diffusivity_ratios/hellmann_and_harvey_2020.py

@@ -9,7 +9,7 @@ def __init__(self, _):
         pass
 
     @staticmethod
-    def ratio_2H(const, temperature):
+    def ratio_2H_heavy_to_light(const, temperature):
         return (
             const.HELLMANN_HARVEY_EQ6_COEFF0
             + const.HELLMANN_HARVEY_EQ6_COEFF1
@@ -19,7 +19,7 @@ def ratio_2H(const, temperature):
         )
 
     @staticmethod
-    def ratio_17O(const, temperature):
+    def ratio_17O_heavy_to_light(const, temperature):
         return (
             const.HELLMANN_HARVEY_EQ7_COEFF0
             + const.HELLMANN_HARVEY_EQ7_COEFF1
@@ -29,7 +29,7 @@ def ratio_17O(const, temperature):
         )
 
     @staticmethod
-    def ratio_18O(const, temperature):
+    def ratio_18O_heavy_to_light(const, temperature):
         return (
             const.HELLMANN_HARVEY_EQ8_COEFF0
             + const.HELLMANN_HARVEY_EQ8_COEFF1

PySDM/physics/isotope_diffusivity_ratios/stewart_1975.py

@@ -17,7 +17,7 @@ def __init__(self, _):
         pass
 
     @staticmethod
-    def ratio_2H(const, temperature):  # pylint: disable=unused-argument
+    def ratio_2H_heavy_to_light(const, temperature):  # pylint: disable=unused-argument
         return (
             (
                 (2 * const.M_1H + const.M_16O)
@@ -30,7 +30,7 @@ def ratio_2H(const, temperature):  # pylint: disable=unused-argument
         ) ** const.ONE_HALF
 
     @staticmethod
-    def ratio_18O(const, temperature):  # pylint: disable=unused-argument
+    def ratio_18O_heavy_to_light(const, temperature):  # pylint: disable=unused-argument
         return (
             ((2 * const.M_1H + const.M_16O) * (const.Md + 2 * const.M_1H + const.M_18O))
             / (

PySDM/physics/isotope_kinetic_fractionation_factors/craig_gordon.py

@@ -13,5 +13,5 @@ def __init__(self, _):
 
     @staticmethod
     def alpha_kinetic(*, relative_humidity, turbulence_parameter_n, delta_diff, theta):
-        """delta_diff = 1 - heavy_to_light_diffusivity_ratio"""
+        """delta_diff = 1 - diffusivity_ratio_heavy_to_light"""
         return 1 + theta * turbulence_parameter_n * delta_diff * (1 - relative_humidity)

PySDM/physics/isotope_kinetic_fractionation_factors/jouzel_and_merlivat_1984.py

@@ -10,12 +10,12 @@ def __init__(self, _):
 
     @staticmethod
     def alpha_kinetic(
-        alpha_equilibrium, saturation_over_ice, heavy_to_light_diffusivity_ratio
+        alpha_equilibrium, saturation_over_ice, diffusivity_ratio_heavy_to_light
     ):
         """eq. (11)"""
         return saturation_over_ice / (
             alpha_equilibrium
-            / heavy_to_light_diffusivity_ratio
+            / diffusivity_ratio_heavy_to_light
             * (saturation_over_ice - 1)
             + 1
         )

PySDM/physics/isotope_ventilation_ratio/__init__.py

@@ -0,0 +1,7 @@
+"""
+Isotope ventilation coefficient formulae
+"""
+
+from PySDM.impl.null_physics_class import Null
+from .neglect import Neglect
+from .brutsaert_1982 import Brutsaert1982

PySDM/physics/isotope_ventilation_ratio/brutsaert_1982.py

@@ -0,0 +1,16 @@
+"""
+based on [Brutsaert 1982](https://doi.org/10.1007/978-94-017-1497-6) Springer Netherlands
+statement about ventilation coefficient for heavy isotopes on pp. 92-93.
+"""
+
+
+class Brutsaert1982:  # pylint: disable=too-few-public-methods
+    def __init__(self, _):
+        pass
+
+    @staticmethod
+    def ratio_heavy_to_light(ventilation_coefficient, diffusivity_ratio_heavy_to_light):
+        """heavy to light isotope ventilation ratio"""
+        return (
+            1 - diffusivity_ratio_heavy_to_light ** (1 / 3)
+        ) / ventilation_coefficient + diffusivity_ratio_heavy_to_light ** (1 / 3)