#include <ERF_MOSTUtils.H>

Collaboration diagram for similarity_funs:

Public Member Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_psi_m2 (amrex::Real zeta) const

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_psi_h2 (amrex::Real zeta) const

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_psi_m (amrex::Real zeta) const

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_psi_h (amrex::Real zeta) const

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_phi_m (amrex::Real zeta) const

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_phi_h (amrex::Real zeta) const

Private Attributes
amrex::Real	beta_m {amrex::Real(5.0)}
	Constants from Dyer, BLM, 1974. More...

amrex::Real	beta_h {amrex::Real(5.0)}
	https://doi.org/amrex::Real(10.1007)/BF00240838 More...

amrex::Real	gamma_m {amrex::Real(16.0)}

amrex::Real	gamma_h {amrex::Real(16.0)}

Detailed Description

Structure of similarity functions for Moeng formulation

Member Function Documentation

◆ calc_phi_h()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_phi_h ( amrex::Real zeta ) const

inline

Compute the Businger-Dyer heat gradient function, the exact companion of calc_psi_h: dtheta/dz = theta_* phi_h(zeta) / (kappa z).

See the note in calc_phi_m about the Jimenez pair.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             return one + beta_h * zeta;
         } else {
             return one / std::sqrt(one - gamma_h * zeta);
         }
     }

◆ calc_phi_m()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_phi_m ( amrex::Real zeta ) const

inline

Compute the Businger-Dyer momentum gradient function.

This is the exact companion of calc_psi_m: the two are tied by psi_m(zeta) = int_0^zeta (1 - phi_m(x))/x dx, so differentiating the MOST wind profile u(z) = (u_star/kappa)[ln(z/z0) - psi_m(z/L)] gives du/dz = u_* phi_m(zeta) / (kappa z).

NOTE: the surface_temp* flux calculators use the Jimenez psi pair (calc_psi_m2/calc_psi_h2) rather than this one; their companions are not implemented yet, so a prescribed-surface-temperature run gets the Businger-Dyer gradient here.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             return one + beta_m * zeta;
         } else {
             return one / std::sqrt(std::sqrt(one - gamma_m * zeta));
         }
     }

◆ calc_psi_h()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_psi_h ( amrex::Real zeta ) const

inline

Compute the Businger-Dyer heat stability correction.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             return -beta_h * zeta;
         } else {
             amrex::Real x = std::sqrt(one - gamma_h * zeta);
             return two * std::log(myhalf * (one + x));
         }
     }

Referenced by compute_if_most_target_vel(), near_surface_diagnostics::evaluate_most_factor(), ImmersedForcingBuildings_Scalar(), ImmersedForcingTerrain_Scalar(), ImmersedForcingTerrain_Xmom(), ImmersedForcingTerrain_Ymom(), surface_flux_charnock::iterate_flux(), surface_flux_mod_charnock::iterate_flux(), surface_flux_donelan::iterate_flux(), surface_flux_wave_coupled::iterate_flux(), surface_flux::iterate_flux(), and surface_flux_eb::iterate_flux().

Here is the caller graph for this function:

◆ calc_psi_h2()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_psi_h2 ( amrex::Real zeta ) const

inline

Compute the Jimenez heat stability correction.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             amrex::Real x = std::pow(one + std::pow(zeta, amrex::Real(1.1)), one/amrex::Real(1.1));
             return ( -amrex::Real(5.3)*std::log(zeta + x) );
         } else {
             amrex::Real x = std::sqrt(one - amrex::Real(16.0)*zeta);
             amrex::Real psi_k_h = two * std::log(myhalf * (one + x));
             amrex::Real y = std::pow(one - amrex::Real(34.0)*zeta, one/three);
             amrex::Real psi_c_h = (three/two)*std::log((y*y + y + one)/three)
                                 - std::sqrt(three)*std::atan((two*y + one)/std::sqrt(three))
                                 + PI/std::sqrt(three);
             return ( (psi_k_h + zeta*zeta*psi_c_h) / (one + zeta*zeta) );
         }
     }

Referenced by surface_temp_charnock::iterate_flux(), surface_temp_mod_charnock::iterate_flux(), surface_temp_donelan::iterate_flux(), surface_temp_wave_coupled::iterate_flux(), surface_temp::iterate_flux(), surface_temp_eb::iterate_flux(), erf_cloud_chamber_wall_flux::most_stability_coordinate(), erf_cloud_chamber_wall_flux::most_wall_coefficients(), and erf_surface_moisture::surface_layer_aerodynamic_resistance().

Here is the caller graph for this function:

◆ calc_psi_m()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_psi_m ( amrex::Real zeta ) const

inline

Compute the Businger-Dyer momentum stability correction.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             return -beta_m * zeta;
         } else {
             amrex::Real x = std::sqrt(std::sqrt(one - gamma_m * zeta));
             return two * std::log(myhalf * (one + x)) + std::log(myhalf * (one + x * x)) -
                    two * std::atan(x) + PIoTwo;
         }
     }

Referenced by compute_if_most_target_vel(), ImmersedForcingBuildings_Scalar(), ImmersedForcingTerrain_Scalar(), ImmersedForcingTerrain_Xmom(), ImmersedForcingTerrain_Ymom(), surface_flux_charnock::iterate_flux(), surface_flux_mod_charnock::iterate_flux(), surface_flux_donelan::iterate_flux(), surface_flux_wave_coupled::iterate_flux(), surface_flux::iterate_flux(), and surface_flux_eb::iterate_flux().

Here is the caller graph for this function:

◆ calc_psi_m2()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real similarity_funs::calc_psi_m2 ( amrex::Real zeta ) const

inline

Compute the Jimenez momentum stability correction.

Parameters

[in] zeta stability coordinate z/L

     {
         if (zeta > 0) {
             amrex::Real x = std::pow(one + std::pow(zeta, amrex::Real(2.5)), one/amrex::Real(2.5));
             return ( -amrex::Real(6.1)*std::log(zeta + x) );
         } else {
             amrex::Real x = std::pow(one - amrex::Real(16.0)*zeta, fourth);
             amrex::Real psi_k_m =  two * std::log(myhalf * (one + x)) + std::log(myhalf * (one + x * x)) -
                                    two * std::atan(x) + PIoTwo;
             amrex::Real y = std::pow(one - amrex::Real(10.0)*zeta, one/three);
             amrex::Real psi_c_m = (three/two)*std::log((y*y + y + one)/three)
                                 - std::sqrt(three)*std::atan((two*y + one)/std::sqrt(three))
                                 + PI/std::sqrt(three);
             return ( (psi_k_m + zeta*zeta*psi_c_m) / (one + zeta*zeta) );
         }
     }