IB.H

Go to the documentation of this file.

#ifndef IB_H
#define IB_H
 
#include "src/core/Physics.H"
#include "src/core/FieldRepo.H"

 
namespace kynema_sgf {
 
class Field;
 
namespace ib {
 
class ImmersedBoundaryModel;

class IB : public Physics::Register<IB>
{
public:
    static std::string identifier() { return "IB"; }
 
    explicit IB(CFDSim& /*sim*/);
 
    ~IB() override;
 
    void
    initialize_fields(int /*level*/, const amrex::Geometry& /*geom*/) override
    {}
 
    void pre_init_actions() override;
 
    void post_init_actions() override;
 
    void post_regrid_actions() override;
 
    void pre_advance_work() override;
 
    void post_advance_work() override;
 
    void pre_pressure_correction_work() override;
 
    void post_pressure_correction_work() override;
 
protected:
 
    [[nodiscard]] int num_ibs() const { return static_cast<int>(m_ibs.size()); }
 
    virtual void prepare_outputs();
 
private:
    static void update_positions();
 
    void update_velocities();
 
    void compute_forces();
 
    CFDSim& m_sim;
 
    std::vector<std::unique_ptr<ImmersedBoundaryModel>> m_ibs;

    Field& m_ib_levelset;

    Field& m_ib_normal;
};
 
} // namespace ib
} // namespace kynema_sgf
 
#endif /* IB_H */