TerminalVelocity< RT > Struct Template Reference

Provide terminal velocity calculations for atmospheric particles. More...

#include <ERF_TerminalVelocity.H>

Collaboration diagram for TerminalVelocity< RT >:

Public Member Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE auto	viscCoeff (const RT a_T) const noexcept
	Compute viscosity coefficient - Pruppacher & Klett, 1997. More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	RogersYau (const RT a_r) const noexcept
	Terminal velocity - Rogers & Yau, 1989. More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	AtlasUlbrich (const RT a_r) const noexcept
	Terminal velocity - Atlas & Ulbrich, 1977. More...
AMREX_GPU_HOST_DEVICE RT	CloudRainShima (const RT a_r, const RT a_rho, const RT a_P, const RT a_T) const noexcept
	Terminal velocity - cloud/rain droplets Shima et al., 2009, "The super-droplet method for the numerical simulation of clouds and precipitation: A particle-based and probabilistic microphysics model coupled with a non- hydrostatic model." Quart. J. Roy. Meteorol. Soc., 135: 1307-1320 https://github.com/Shima-Lab/SCALE-SDM_BOMEX_Sato2018/blob/master/contrib/SDM/sdm_motion.f90. More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	vz_b_x (const RT a_x) const noexcept
	Compute vz_b(x) More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	vz_c_x (const RT a_x) const noexcept
	Compute vc_b(x) More...
AMREX_GPU_HOST_DEVICE RT	IceBohm (const RT a_m, const RT a_a, const RT a_c, const RT a_rhoi, const RT a_rho, const RT a_T) const noexcept
	Terminal velocity - ice particles (Bohm's theory) Ref: Shima et al., 2020, "Predicting the morphology of ice particles in deep convection using the super-droplet method: development and evaluation of SCALE-SDM 0.2.5-2.2.0, -2.2.1, and -2.2.2.", Geosci. Model Dev., 13, 4107–4157, 2020. More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	MeltBlendWeight (const RT a_ell) const noexcept
	Frick et al. (2013, Geosci. Model Dev. 6, Eq. 15) liquid-water-fraction weight for the fall speed of a melting (mixed) particle. Rises from 0 (dry ice) to 1 (fully melted) as the meltwater mass fraction a_ell goes 0 -> 1. More...
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT	MixedPhaseVterm (const RT a_v_ice, const RT a_v_rain, const RT a_ell) const noexcept
	Terminal velocity of a melting (mixed) ice-water particle, interpolated between the dry-ice (Bohm) fall speed and the equivalent-raindrop fall speed by the meltwater mass fraction (Frick et al. 2013, Eq. 14). More...

Public Attributes
RT	m_rho_w
RT	m_rho_i

Detailed Description

template<typename RT>
struct TerminalVelocity< RT >

Provide terminal velocity calculations for atmospheric particles.

Template Parameters

RT	Real-valued type.

Member Function Documentation

AtlasUlbrich()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::AtlasUlbrich ( const RT  a_r ) const

inlinenoexcept

Terminal velocity - Atlas & Ulbrich, 1977.

    {
        static const RT a = RT(3.778);
        static const RT b = RT(0.67);
        RT D = RT(2)*a_r*RT(1000); // meter -> mm
        return a*(std::exp(b*std::log(D)));
    }

CloudRainShima()

template<typename RT >

AMREX_GPU_HOST_DEVICE RT TerminalVelocity< RT >::CloudRainShima ( const RT  a_r, const RT  a_rho, const RT  a_P, const RT  a_T  ) const

inlinenoexcept

Terminal velocity - cloud/rain droplets Shima et al., 2009, "The super-droplet method for the numerical simulation of clouds and precipitation: A particle-based and probabilistic microphysics model coupled with a non- hydrostatic model." Quart. J. Roy. Meteorol. Soc., 135: 1307-1320 https://github.com/Shima-Lab/SCALE-SDM_BOMEX_Sato2018/blob/master/contrib/SDM/sdm_motion.f90.

Parameters

a_r	radius
a_rho	density
a_P	pressure
a_T	temperature

    {
        RT lb4l = RT(6.62E-6);    // base length[cm] for mean free path  of air molecules
        RT visb4l = RT(1.818E-4); // base dynamic viscosity[g/(cm*s)] for mean free path of air molecules
        RT pb4l = RT(1013.25);    // base pressure[hPa] for mean free path of air molecules
        RT tb4l = RT(293.15);     // base temperature[20C] for mean free path of air molecules
 
        RT P_hPa = a_P / RT(100.0);           // [Pa] -> [hPa]
        RT diameter_cm = RT(2.0)*std::min(a_r,RT(0.004))*RT(100.0);   // [m] -> [cm]
        RT rho_mat_cgs = m_rho_w/RT(1000.0);  // [kg/m^3] -> [g/cm^3]
        RT rho_air_cgs = a_rho/RT(1000.0);    // [kg/m^3] -> [g/cm^3]
        RT grav_cgs = CONST_GRAV * RT(100);   // [m/s^2] -> [cm/s^2]
 
        RT gxdrow = grav_cgs * ( rho_mat_cgs - rho_air_cgs );
        RT viscosity = viscCoeff(a_T);
 
        RT v_terminal = RT(0.0);
        if (diameter_cm < RT(1.9e-3)) {
            // small cloud droplets
            RT sd_l = lb4l * (viscosity/visb4l) * (pb4l/P_hPa) * std::sqrt(a_T/tb4l);
            RT c1  = gxdrow / (RT(18.0)*viscosity);
            RT csc = RT(1.0) + RT(2.510) * (sd_l/diameter_cm);
            v_terminal = c1 * csc * diameter_cm * diameter_cm;
        } else if (diameter_cm < RT(1.07e-1)) {
            // large cloud droplets and small raindrops
            RT sd_l = lb4l * (viscosity/visb4l) * (pb4l/P_hPa) * std::sqrt(a_T/tb4l);
            RT c2 = rho_air_cgs * (RT(4.0)*gxdrow)/(RT(three)*viscosity*viscosity);
            RT nda = c2 * diameter_cm * diameter_cm * diameter_cm;
            RT csc = RT(1.0) + RT(2.510) * ( sd_l / diameter_cm );
            RT nre = csc * std::exp(vz_b_x(std::log(nda)));
            v_terminal = (viscosity*nre)/(rho_air_cgs*diameter_cm);
        } else {
            // large raindrops
            RT tau = RT(1.0) - a_T/RT(647.0960);
            RT sigma = RT(0.2358) * std::exp( RT(1.2560)*std::log(tau) )
                                * ( RT(1.0) - RT(0.6250)*tau );
            if( a_T<RT(267.5) ) {
                tau = std::tanh( (a_T-RT(243.9))/RT(35.35) );
                sigma = sigma - RT(2.854E-3) * tau + RT(1.666E-3);
            }
            sigma = sigma * RT(1.E+3); // N/m => g/s^2
            RT c3 = (RT(4.0)*gxdrow)/(RT(three)*sigma);
            RT bond = c3 * diameter_cm * diameter_cm;
            RT npp =    (rho_air_cgs*rho_air_cgs) * (sigma*sigma*sigma)
                        / (gxdrow*viscosity*viscosity*viscosity*viscosity);
            npp = std::exp( std::log(npp)/RT(6.0) );
            RT nre = npp * std::exp(vz_c_x(std::log(bond*npp)));
            v_terminal = (viscosity*nre)/(rho_air_cgs*diameter_cm);
        }
 
        v_terminal /= RT(100.0); // [cm/s] -> [m/s]
        return v_terminal;
    }

Here is the call graph for this function:

IceBohm()

template<typename RT >

AMREX_GPU_HOST_DEVICE RT TerminalVelocity< RT >::IceBohm ( const RT  a_m, const RT  a_a, const RT  a_c, const RT  a_rhoi, const RT  a_rho, const RT  a_T  ) const

inlinenoexcept

Terminal velocity - ice particles (Bohm's theory) Ref: Shima et al., 2020, "Predicting the morphology of ice particles in deep convection using the super-droplet method: development and evaluation of SCALE-SDM 0.2.5-2.2.0, -2.2.1, and -2.2.2.", Geosci. Model Dev., 13, 4107–4157, 2020.

and

Code: https://doi.org/10.5281/zenodo.3483650 File: [...]/contrib/SDM/sdm_motion.f90 Copyright (c) 2012-2014, Team SCALE, All rights reserved.

Parameters

a_m	mass
a_a	ice particle equatorial radius
a_c	ice particle polar radius
a_rhoi	ice particle density
a_rho	density
a_T	temperature

    {
        RT mu = RT(0.1) * viscCoeff(a_T); // viscCoeff is CGS poise; IceBohm is SI -> Pa.s
 
        auto phi = a_c / a_a; // aspect ratio
        auto d = RT(2.0) * a_a; // characteristic length
        auto q = std::exp( std::exp(-phi) * std::log(a_rhoi/m_rho_i) ); // area ratio
 
        RT X0 = RT(2.8e06);
        auto t3 = RT(PI) * mu*mu * std::max(phi,RT(1.0)) * std::max(std::sqrt(std::sqrt(q)), q);
        auto X = (RT(8.0)*a_m*RT(CONST_GRAV)*a_rho) / t3;
        auto Xp = X * (RT(1.0) + (X/X0)*(X/X0)) / (RT(1.0) + RT(1.6)*(X/X0)*(X/X0));
 
        auto k = std::min( std::max(RT(0.82)+RT(0.18)*phi, RT(0.85)),
                           std::min( RT(0.37) + RT(0.63)/std::sqrt(phi),
                                     RT(1.33)/(std::max(std::log(phi),RT(0.0))+RT(1.19)) ));
 
        auto phi_poly = RT(3.76)-RT(8.41)*phi+RT(9.18)*phi*phi-RT(3.53)*phi*phi*phi;
        auto TVGamma = std::max( RT(1.0), std::min(RT(1.98), phi_poly) );
 
        auto C_DPS = std::max(RT(0.292)*k*TVGamma, RT(0.492) - RT(0.2)/std::sqrt(phi));
        auto C_DP  = std::max(RT(1.0), q*(RT(1.46)*q-RT(0.46))) * C_DPS;
        auto C_DO = RT(4.5) * k*k * std::max(phi, RT(1.0));
 
        auto t2 = RT(1.0) + C_DP/(RT(6.0)*k) * std::sqrt(Xp/C_DP);
        auto beta = std::sqrt(t2) - RT(1.0);
        auto TVgamma = (C_DO-C_DP) / (RT(4.0)*C_DP);
        auto t1 = RT(1.0) + (RT(2.0)*beta*std::exp(-beta*TVgamma)) / ((RT(2.0)+beta)*(RT(1.0)+beta));
        auto N_Re = (RT(6.0)*k/C_DP) * beta*beta * t1;
 
        auto vterm = (mu*N_Re)/(a_rho*d);
        return vterm;
    }

Here is the call graph for this function:

MeltBlendWeight()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::MeltBlendWeight ( const RT  a_ell ) const

inlinenoexcept

Frick et al. (2013, Geosci. Model Dev. 6, Eq. 15) liquid-water-fraction weight for the fall speed of a melting (mixed) particle. Rises from 0 (dry ice) to 1 (fully melted) as the meltwater mass fraction a_ell goes 0 -> 1.

Parameters

a_ell

meltwater (liquid) mass fraction in [0,1]

    {
        return RT(0.246)*a_ell + RT(0.754)*a_ell*a_ell*a_ell*a_ell*a_ell*a_ell*a_ell;
    }

Referenced by TerminalVelocity< RT >::MixedPhaseVterm().

Here is the caller graph for this function:

MixedPhaseVterm()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::MixedPhaseVterm ( const RT  a_v_ice, const RT  a_v_rain, const RT  a_ell  ) const

inlinenoexcept

Terminal velocity of a melting (mixed) ice-water particle, interpolated between the dry-ice (Bohm) fall speed and the equivalent-raindrop fall speed by the meltwater mass fraction (Frick et al. 2013, Eq. 14).

Parameters

a_v_ice	dry-ice (Bohm) fall speed of the ice core
a_v_rain	fall speed of the equivalent raindrop
a_ell	meltwater (liquid) mass fraction in [0,1]

    {
        return a_v_ice + (a_v_rain - a_v_ice) * MeltBlendWeight(a_ell);
    }

Here is the call graph for this function:

RogersYau()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::RogersYau ( const RT  a_r ) const

inlinenoexcept

Terminal velocity - Rogers & Yau, 1989.

    {
        static const RT k1 = RT(1.233e8); // m^{-1} s^{-1}
        return (k1 * a_r * a_r);
    }

viscCoeff()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE auto TerminalVelocity< RT >::viscCoeff ( const RT  a_T ) const

inlinenoexcept

Compute viscosity coefficient - Pruppacher & Klett, 1997.

Parameters

a_T	temperature [K]

    {
        RT T_degC = a_T - RT(273.15); // [K] -> [deg Celsius]
        RT viscosity = RT(0.0);
        if (T_degC >= RT(0)) {
            viscosity = (RT(1.7180) + RT(4.9e-3)*T_degC) * RT(1.0e-4);
        } else {
            viscosity = (  RT(1.7180) + RT(4.9e-3)*T_degC
                         - RT(1.2e-5)*T_degC*T_degC     ) * RT(1.0e-4);
        }
        return viscosity;
    }

Referenced by TerminalVelocity< RT >::CloudRainShima(), and TerminalVelocity< RT >::IceBohm().

Here is the caller graph for this function:

vz_b_x()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::vz_b_x ( const RT  a_x ) const

inlinenoexcept

Compute vz_b(x)

Parameters

a_x

x

    {
        RT x1 = a_x;
        RT x2 = x1 * x1;
        RT x3 = x1 * x2;
        RT y = - RT(0.318657E+1)
               + RT(0.9926960)    * x1
               - RT(0.153193E-2)  * x2
               - RT(0.987059E-3)  * x3
               - RT(0.578878e-3)  * x2*x2
               + RT(0.855176E-4)  * x3*x2
               - RT(0.327815E-5)  * x3*x3;
        return y;
    }

Referenced by TerminalVelocity< RT >::CloudRainShima().

Here is the caller graph for this function:

vz_c_x()

template<typename RT >

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE RT TerminalVelocity< RT >::vz_c_x ( const RT  a_x ) const

inlinenoexcept

Compute vc_b(x)

Parameters

a_x

x

    {
        RT x1 = a_x;
        RT x2 = x1 * x1;
        RT x3 = x1 * x2;
        RT y = - RT(0.500015E+1)
               + RT(0.523778E+1)  * x1
               - RT(0.204914E+1)  * x2
               + RT(0.47529400)   * x3
               - RT(0.542819E-1)  * x2*x2
               + RT(0.238449E-2)  * x3*x2;
        return y;
    }

Referenced by TerminalVelocity< RT >::CloudRainShima().

Here is the caller graph for this function:

Member Data Documentation

m_rho_i

template<typename RT >

RT TerminalVelocity< RT >::m_rho_i

density of ice

Referenced by TerminalVelocity< RT >::IceBohm().

m_rho_w

template<typename RT >

RT TerminalVelocity< RT >::m_rho_w

density of water

Referenced by TerminalVelocity< RT >::CloudRainShima().

---

The documentation for this struct was generated from the following file:

• Source/MaterialProperties/ERF_TerminalVelocity.H