MMSForcing.H

Go to the documentation of this file.

#ifndef MMSFORCING_H
#define MMSFORCING_H
 
#include "src/core/FieldRepo.H"
#include "src/equation_systems/icns/MomentumSource.H"
 
namespace kynema_sgf::pde::icns::mms {
 
class MMSForcing : public MomentumSource::Register<MMSForcing>
{
public:
    static std::string identifier() { return "MMSForcing"; }
 
    explicit MMSForcing(const CFDSim& sim);
 
    ~MMSForcing() override = default;
 
    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_mms_vel_source;
};
} // namespace kynema_sgf::pde::icns::mms
 
#endif /* MMSFORCING_H */