ABLMeanBoussinesq.H

Go to the documentation of this file.

#ifndef ABLMEANBOUSSINESQ_H
#define ABLMEANBOUSSINESQ_H
 
#include "src/core/FieldRepo.H"
#include "src/transport_models/TransportModel.H"
#include "src/equation_systems/icns/MomentumSource.H"
#include "src/utilities/FieldPlaneAveraging.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::pde::icns {

class ABLMeanBoussinesq : public MomentumSource::Register<ABLMeanBoussinesq>
{
public:
    static std::string identifier() { return "ABLMeanBoussinesq"; }
 
    explicit ABLMeanBoussinesq(const CFDSim& sim);
 
    ~ABLMeanBoussinesq() override;
 
    void operator()(
        int lev,
        const amrex::MFIter& mfi,
        const amrex::Box& bx,
        FieldState fstate,
        const amrex::Array4<amrex::Real>& src_term) const override;
 
    void mean_temperature_init(const FieldPlaneAveraging& /*tavg*/);
 
    void mean_temperature_update(const FieldPlaneAveraging& /*tavg*/);
 
private:
    const amrex::AmrCore& m_mesh;
 
    amrex::Vector<amrex::Real> m_gravity{0.0_rt, 0.0_rt, -9.81_rt};
 
    amrex::Gpu::DeviceVector<amrex::Real> m_theta_ht;
    amrex::Gpu::DeviceVector<amrex::Real> m_theta_vals;

    const transport::TransportModel& m_transport;
 
    int m_axis{2};
 
    bool m_const_profile{false};

    void read_temperature_profile(const std::string& profile_file_name);
};
} // namespace kynema_sgf::pde::icns
 
#endif /* ABLMEANBOUSSINESQ_H */