TurbModel.H

Go to the documentation of this file.

#ifndef TURBMODEL_H
#define TURBMODEL_H
 
#include "amr-wind/turbulence/TurbulenceModel.H"
#include "amr-wind/CFDSim.H"
 
namespace amr_wind::turbulence {
 
template <typename Transport>
class TurbModel : public TurbulenceModel
{
public:
    using TransportType = Transport;
 
    explicit TurbModel(CFDSim& sim) : m_sim(sim), m_transport(sim) {}
 
    void register_mueff_field(Field& mueff) override { m_mueff = &mueff; }
 
    void register_alphaeff_field(Field& alphaeff) override
    {
        m_alphaeff = &alphaeff;
    }
 
    std::unique_ptr<ScratchField> mu() override { return m_transport.mu(); }
 
    std::unique_ptr<ScratchField> alpha() override
    {
        return m_transport.alpha();
    }
 
    std::unique_ptr<ScratchField>
    scalar_diffusivity(const std::string& name) override
    {
        return m_transport.scalar_diffusivity(name);
    }
 
    Field& mueff() override
    {
        AMREX_ASSERT(m_mueff != nullptr);
        return *m_mueff;
    }
 
    Field& alphaeff() override
    {
        AMREX_ASSERT(m_alphaeff != nullptr);
        return *m_alphaeff;
    }
 
protected:
    CFDSim& m_sim;
 
    Transport m_transport;
 
    Field* m_mueff{nullptr};
 
    Field* m_alphaeff{nullptr};
};
 
} // namespace amr_wind::turbulence
 
#endif /* TURBMODEL_H */