Physics.H

Go to the documentation of this file.

#ifndef PHYSICS_H
#define PHYSICS_H
 
#include "amr-wind/core/Factory.H"
#include "amr-wind/core/CollMgr.H"
#include "AMReX_MultiFab.H"
#include "AMReX_Geometry.H"
 
namespace amr_wind {
 
class CFDSim;
 
class PhysicsOld
{
public:
    virtual ~PhysicsOld() = default;
 
    virtual std::string physics_name() const = 0;
 
    virtual void
    initialize_fields(int level, const amrex::Geometry& geom) const = 0;
 
    virtual void pre_advance_work() = 0;
 
    virtual void add_momentum_sources(
        const amrex::Geometry& geom,
        const amrex::MultiFab& density,
        const amrex::MultiFab& velocity,
        const amrex::MultiFab& scalars,
        amrex::MultiFab& vel_forces) const = 0;
};
 
class Physics : public Factory<Physics, CFDSim&>
{
public:
    static std::string base_identifier() { return "Physics"; }
 
    ~Physics() override = default;
 
    virtual void initialize_fields(int level, const amrex::Geometry& geom) = 0;
 
    virtual void pre_init_actions() {}
 
    virtual void post_init_actions() = 0;
 
    virtual void post_regrid_actions() = 0;
 
    virtual void pre_advance_work() = 0;
 
    virtual void pre_predictor_work() {}
 
    virtual void post_advance_work() = 0;
 
    virtual void pre_pressure_correction_work() {}
 
    virtual void post_pressure_correction_work() {}
};
 
class PhysicsMgr : public CollMgr<PhysicsMgr, Physics>
{
public:
    template <typename T>
    T& get()
    {
        AMREX_ASSERT(contains(T::identifier()));
        return dynamic_cast<T&>(operator()(T::identifier()));
    }
 
    template <typename T>
    const T& get() const
    {
        AMREX_ASSERT(contains(T::identifier()));
        return dynamic_cast<const T&>(operator()(T::identifier()));
    }
};
 
} // namespace amr_wind
 
#endif /* PHYSICS_H */