diff --git a/amr-wind/equation_systems/icns/source_terms/DragForcing.cpp b/amr-wind/equation_systems/icns/source_terms/DragForcing.cpp index e54295ec4d..473ce6d78a 100644 --- a/amr-wind/equation_systems/icns/source_terms/DragForcing.cpp +++ b/amr-wind/equation_systems/icns/source_terms/DragForcing.cpp @@ -27,10 +27,10 @@ DragForcing::DragForcing(const CFDSim& sim) pp.query("sponge_south", m_sponge_south); pp.query("sponge_north", m_sponge_north); pp.query("is_laminar", m_is_laminar); - amrex::ParmParse pp_abl("ABL"); + amrex::ParmParse pp_abl("ABL"); pp_abl.query("mol_length", m_mol_length); pp_abl.query("mo_gamma_m", m_gamma_m); - pp_abl.query("mo_beta_m", m_beta_m); + pp_abl.query("mo_beta_m", m_beta_m); const auto& phy_mgr = m_sim.physics_manager(); if (phy_mgr.contains("ABL")) { const auto& abl = m_sim.physics_manager().get(); @@ -91,7 +91,7 @@ void DragForcing::operator()( } else { const amrex::Real x = std::sqrt(std::sqrt(1 - m_beta_m * zeta)); psi_m_cell = 2.0 * std::log(0.5 * (1.0 + x)) + log(0.5 * (1 + x * x)) - - 2.0 * std::atan(x) + utils::half_pi(); + 2.0 * std::atan(x) + utils::half_pi(); } //! Cell Value const auto& prob_lo = geom.ProbLoArray(); diff --git a/amr-wind/equation_systems/tke/source_terms/KransAxell.H b/amr-wind/equation_systems/tke/source_terms/KransAxell.H index 037e391fa2..18b2c83be1 100644 --- a/amr-wind/equation_systems/tke/source_terms/KransAxell.H +++ b/amr-wind/equation_systems/tke/source_terms/KransAxell.H @@ -40,7 +40,7 @@ private: amrex::Real m_kappa{0.41}; amrex::Real m_mol_length{10000}; amrex::Real m_gamma_m{5.0}; - amrex::Real m_beta_m{16.0}; + amrex::Real m_beta_m{16.0}; amrex::Real m_sponge_start{600}; amrex::Real m_ref_tke{1e-10}; amrex::Vector m_gravity{0.0, 0.0, -9.81}; diff --git a/amr-wind/equation_systems/tke/source_terms/KransAxell.cpp b/amr-wind/equation_systems/tke/source_terms/KransAxell.cpp index 590ccdf7f3..545b72c60b 100644 --- a/amr-wind/equation_systems/tke/source_terms/KransAxell.cpp +++ b/amr-wind/equation_systems/tke/source_terms/KransAxell.cpp @@ -26,7 +26,7 @@ KransAxell::KransAxell(const CFDSim& sim) pp.query("surface_temp_flux", m_heat_flux); pp.query("mol_length", m_mol_length); pp.query("mo_gamma_m", m_gamma_m); - pp.query("mo_beta_m", m_beta_m); + pp.query("mo_beta_m", m_beta_m); pp.query("meso_sponge_start", m_sponge_start); { amrex::ParmParse pp_incflow("incflo"); @@ -71,7 +71,7 @@ void KransAxell::operator()( } else { const amrex::Real x = std::sqrt(std::sqrt(1 - m_beta_m * zeta)); psi_m = 2.0 * std::log(0.5 * (1.0 + x)) + log(0.5 * (1 + x * x)) - - 2.0 * std::atan(x) + 2 * std::atan(1.0); + 2.0 * std::atan(x) + 2 * std::atan(1.0); } amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { amrex::Real bcforcing = 0; @@ -80,7 +80,7 @@ void KransAxell::operator()( const amrex::Real z = problo[2] + (k + 0.5) * dx[2]; if (k == 0) { const amrex::Real m = std::sqrt(ux * ux + uy * uy); - const amrex::Real ustar = m * kappa / ( std::log(z / z0) + psi_m); + const amrex::Real ustar = m * kappa / (std::log(z / z0) + psi_m); const amrex::Real rans_b = std::pow( std::max(heat_flux, 0.0) * kappa * z / std::pow(Cmu, 3), (2.0 / 3.0)); @@ -115,7 +115,8 @@ void KransAxell::operator()( const amrex::Real uy = vel(i, j, k, 1); const amrex::Real z = 0.5 * dx[2]; amrex::Real m = std::sqrt(ux * ux + uy * uy); - const amrex::Real ustar = m * kappa / ( std::log(z / z0) + psi_m); + const amrex::Real ustar = + m * kappa / (std::log(z / z0) + psi_m); const amrex::Real rans_b = std::pow( std::max(heat_flux, 0.0) * kappa * z / std::pow(Cmu, 3), (2.0 / 3.0)); diff --git a/amr-wind/wind_energy/ABLWallFunction.H b/amr-wind/wind_energy/ABLWallFunction.H index 5a13c8f3ef..40b0e807be 100644 --- a/amr-wind/wind_energy/ABLWallFunction.H +++ b/amr-wind/wind_energy/ABLWallFunction.H @@ -105,7 +105,7 @@ private: std::string m_wall_shear_stress_type{"moeng"}; std::string m_surface_flux_method{"uniform"}; amrex::Real m_mol_length{10000}; - amrex::Real m_surface_roughness{0.1}; + amrex::Real m_surface_roughness{0.1}; }; } // namespace amr_wind