BoussinesqBuoyancy.H

Go to the documentation of this file.

#ifndef BOUSSINESQBUOYANCY_H
#define BOUSSINESQBUOYANCY_H
 
#include "amr-wind/core/FieldRepo.H"
#include "amr-wind/equation_systems/icns/MomentumSource.H"
 
namespace amr_wind::pde::icns {
 
class BoussinesqBuoyancy : public MomentumSource::Register<BoussinesqBuoyancy>
{
public:
    static std::string identifier() { return "BoussinesqBuoyancy"; }
 
    explicit BoussinesqBuoyancy(const CFDSim& sim);
 
    ~BoussinesqBuoyancy() override;
 
    void operator()(
        const int lev,
        const amrex::MFIter& mfi,
        const amrex::Box& bx,
        const FieldState fstate,
        const amrex::Array4<amrex::Real>& src_term) const override;
 
private:
    const Field& m_temperature;
    const Field* m_vof;
 
    amrex::Vector<amrex::Real> m_gravity{0.0, 0.0, -9.81};
 
    amrex::Real m_ref_theta{300.0};
 
    amrex::Real m_beta{0.0};
 
    bool m_is_vof{false};
};
} // namespace amr_wind::pde::icns
 
#endif /* BOUSSINESQBUOYANCY_H */