amr-wind/equation_systems/icns/source_terms/BoussinesqBuoyancy.cpp
/** Boussinesq buoyancy source term for ABL simulations
 *
 *  Reads in the following parameters from `BoussinesqBuoyancy` namespace:
 *
 *  - `reference_temperature` (Mandatory) temperature (`T0`) in Kelvin
 *  - `thermal_expansion_coeff` Optional, default = `1.0 / T0`
 *  - `gravity` acceleration due to gravity (m/s)
 */
Boussinesq buoyancy source term for ABL simulations
Reads in the following parameters from `BoussinesqBuoyancy` namespace.


Boussinesq buoyancy source term for ABL simulations Reads in the following parameters from ` BoussinesqBuoyancy ` namespace .
BoussinesqBuoyancy::BoussinesqBuoyancy(const CFDSim& sim)
    : m_temperature(sim.repo().get_field("temperature"))
{
    amrex::ParmParse pp_boussinesq_buoyancy(this->identifier());
    pp_boussinesq_buoyancy.get("reference_temperature", m_ref_theta);
    if (pp_boussinesq_buoyancy.contains("thermal_expansion_coeff")) {
        pp_boussinesq_buoyancy.get("thermal_expansion_coeff", m_beta);
    } else {
        m_beta = 1.0 / m_ref_theta;
    }
    amrex::ParmParse pp_incflo("incflo");
    pp_incflo.queryarr("gravity", m_gravity);
}
BoussinesqBuoyancy :: BoussinesqBuoyancy ( const CFDSim & sim ) : m_temperature ( sim . repo ( ) . get_field ( " " ) ) { amrex :: ParmParse pp_boussinesq_buoyancy ( this -> identifier ( ) ) ; pp_boussinesq_buoyancy . get ( " " , m_ref_theta ) ; if ( pp_boussinesq_buoyancy . contains ( " " ) ) { pp_boussinesq_buoyancy . get ( " " , m_beta ) ; } else { m_beta = 1.0 / m_ref_theta ; } amrex :: ParmParse pp_incflo ( " " ) ; pp_incflo . queryarr ( " " , m_gravity ) ; }
Boussinesq buoyancy source term for ABL simulations
  Reads in the following parameters from `BoussinesqBuoyancy` namespace:
Boussinesq buoyancy source term for ABL simulations Reads in the following parameters from ` BoussinesqBuoyancy ` namespace :
// FIXME: gravity in `incflo` namespace