temperature.H

Go to the documentation of this file.

#ifndef TEMPERATURE_H
#define TEMPERATURE_H
 
#include "src/equation_systems/PDETraits.H"
#include "src/equation_systems/SchemeTraits.H"
#include "src/equation_systems/PDEHelpers.H"
#include "src/equation_systems/PDE.H"
#include "src/equation_systems/temperature/TemperatureSource.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::pde {



struct Temperature : ScalarTransport
{
    using MLDiffOp = amrex::MLABecLaplacian;
    using SrcTerm = TemperatureSource;
 
    static std::string pde_name() { return "Temperature"; }
    static std::string var_name() { return "temperature"; }
 
    static constexpr amrex::Real default_bc_value = 0.0_rt;
 
    static constexpr int ndim = 1;
    static constexpr bool multiply_rho = true;
    static constexpr bool has_diffusion = true;
    static constexpr bool need_nph_state = true;
};

template <>
struct TurbulenceOp<Temperature>
{
    TurbulenceOp(turbulence::TurbulenceModel& tmodel, PDEFields& fields)
        : m_tmodel(tmodel), m_fields(fields)
    {}
 
    void operator()() { m_tmodel.update_alphaeff(m_fields.mueff); }
 
    turbulence::TurbulenceModel& m_tmodel;
    PDEFields& m_fields;
};
 
} // namespace kynema_sgf::pde
 
#endif /* TEMPERATURE_H */