IBModel.H

Go to the documentation of this file.

#ifndef IBMODEL_H
#define IBMODEL_H
 
#include "amr-wind/core/Factory.H"
#include "amr-wind/immersed_boundary/IBTypes.H"
#include "amr-wind/immersed_boundary/IBOps.H"
#include "amr-wind/core/MultiParser.H"
 
namespace amr_wind {
 
class CFDSim;
 
namespace ib {

class ImmersedBoundaryModel
    : public Factory<
          ImmersedBoundaryModel,
          CFDSim&,
          const std::string&,
          const int>
{
public:
    static std::string base_identifier() { return "ImmersedBoundaryModel"; }
 
    ~ImmersedBoundaryModel() override = default;
 
    [[nodiscard]] virtual const std::string& label() const = 0;
 
    [[nodiscard]] virtual int id() const = 0;
 
    virtual void read_inputs(const ::amr_wind::utils::MultiParser&) = 0;
 
    virtual void init_ib() = 0;
 
    virtual void update_positions() = 0;
 
    virtual void update_velocities() = 0;
 
    virtual void compute_forces() = 0;
 
    virtual void prepare_outputs(const std::string&) = 0;
 
    virtual void write_outputs() = 0;
};

template <typename GeomTrait>
class IBModel : public ImmersedBoundaryModel::Register<IBModel<GeomTrait>>
{
private:
    typename GeomTrait::DataType m_data;

    ops::ProcessOutputsOp<GeomTrait> m_out_op;
 
public:
    static std::string identifier() { return GeomTrait::identifier(); }

    IBModel(CFDSim& sim, const std::string& label, const int id)
        : m_data(sim, label, id), m_out_op(m_data)
    {}

    [[nodiscard]] const std::string& label() const override
    {
        return m_data.info().label;
    }

    [[nodiscard]] int id() const override { return m_data.info().id; }
 
    void read_inputs(const ::amr_wind::utils::MultiParser& pp) override
    {
        ops::ReadInputsOp<GeomTrait>()(m_data, pp);
        m_out_op.read_io_options(pp);
    }
 
    void update_positions() override {}
 
    void update_velocities() override { ops::UpdateVelOp<GeomTrait>()(m_data); }
 
    void compute_forces() override { ops::ComputeForceOp<GeomTrait>()(m_data); }
 
    void prepare_outputs(const std::string& out_dir) override
    {
        m_out_op.prepare_outputs(out_dir);
    }
 
    void write_outputs() override { m_out_op.write_outputs(); }
 
    void init_ib() override { ops::InitDataOp<GeomTrait>()(m_data); }
};
 
} // namespace ib
} // namespace amr_wind
 
#endif /* IBMODEL_H */