ABLModulatedPowerLaw.H

Go to the documentation of this file.

#ifndef ABLMODULATEDPOWERLAW_H
#define ABLMODULATEDPOWERLAW_H
 
#include "amr-wind/core/Field.H"
#include "amr-wind/CFDSim.H"
#include "AMReX_Gpu.H"
#include "amr-wind/utilities/trig_ops.H"
#include "amr-wind/core/vs/vector_space.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace amr_wind {

class ABLModulatedPowerLaw
{
 
public:
    explicit ABLModulatedPowerLaw(CFDSim& /*sim*/);

    void post_init_actions();
 
    void pre_advance_work();
 
    void post_advance_work();
 
    void set_velocity(
        const int lev,
        const amrex::Real time,
        const Field& fld,
        amrex::MultiFab& mfab,
        const int dcomp = 0,
        const int orig_comp = 0) const;
 
    void set_temperature(
        const int lev,
        const amrex::Real time,
        const Field& fld,
        amrex::MultiFab& mfab) const;
 
    bool is_active() const { return m_activate_mpl; }
 
private:
    const CFDSim& m_sim;
    const amr_wind::SimTime& m_time;
    const FieldRepo& m_repo;
    const amrex::AmrCore& m_mesh;
    Field& m_velocity;
    Field& m_temperature;
 
    amrex::Real m_zoffset{0.0_rt};
    amrex::Real m_zref{90};
    amrex::Real m_shear_exp{0.1_rt};
    amrex::Real m_umax_factor{1.2_rt};
    amrex::Real m_bulk_velocity{15.0_rt};
    amrex::Real m_shearlayer_height{600.0_rt};
    amrex::Real m_shearlayer_smear_thickness{30.0_rt};
    amrex::Real m_wind_speed{8.0_rt};
    amrex::Real m_wind_direction{270.0_rt};
    amrex::Vector<amrex::Real> m_uvec{8.0_rt, 0.0_rt, 0.0_rt};
 
    amrex::Real m_start_time{0.0_rt};
    amrex::Real m_stop_time{std::numeric_limits<amrex::Real>::max()};
    amrex::Real m_degrees_per_sec{0.02_rt};

    amrex::Vector<amrex::Real> m_theta_heights;
    amrex::Vector<amrex::Real> m_theta_values;
 
    // Device copies of the above arrays
    amrex::Gpu::DeviceVector<amrex::Real> m_thht_d;
    amrex::Gpu::DeviceVector<amrex::Real> m_thvv_d;
 
    amrex::Real m_delta_t{0.8_rt};
    amrex::Real m_theta_cutoff_height{250.0_rt};
    amrex::Real m_theta_gauss_mean{0.0_rt};
    amrex::Real m_theta_gauss_var{1.0_rt};
 
    bool m_activate_mpl{false};
};
 
} // namespace amr_wind
 
#endif /* ABLMODULATEDPOWERLAW_H */