similarity_funs Struct Reference

#include <ERF_MOSTStress.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

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_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(), surface_flux_charnock::iterate_flux(), surface_flux_mod_charnock::iterate_flux(), surface_flux_donelan::iterate_flux(), surface_temp_donelan::iterate_flux(), surface_flux_wave_coupled::iterate_flux(), surface_temp_wave_coupled::iterate_flux(), surface_flux::iterate_flux(), surface_flux_eb::iterate_flux(), make_mom_sources(), and make_sources().

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::iterate_flux(), and surface_temp_eb::iterate_flux().

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(), surface_flux_charnock::iterate_flux(), surface_flux_mod_charnock::iterate_flux(), surface_flux_donelan::iterate_flux(), surface_temp_donelan::iterate_flux(), surface_flux_wave_coupled::iterate_flux(), surface_temp_wave_coupled::iterate_flux(), surface_flux::iterate_flux(), surface_flux_eb::iterate_flux(), make_mom_sources(), and make_sources().

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) );
        }
    }

Referenced by surface_temp_charnock::iterate_flux(), surface_temp_mod_charnock::iterate_flux(), surface_temp::iterate_flux(), and surface_temp_eb::iterate_flux().

Here is the caller graph for this function:

Member Data Documentation

beta_h

amrex::Real similarity_funs::beta_h {amrex::Real(5.0)}

private

https://doi.org/amrex::Real(10.1007)/BF00240838

Referenced by calc_psi_h().

beta_m

amrex::Real similarity_funs::beta_m {amrex::Real(5.0)}

private

Constants from Dyer, BLM, 1974.

Referenced by calc_psi_m().

gamma_h

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

private

Referenced by calc_psi_h().

gamma_m

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

private

Referenced by calc_psi_m().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H