vof_bciface.H

Go to the documentation of this file.

#ifndef VOF_BCIFACE_H
#define VOF_BCIFACE_H
 
#include "src/boundary_conditions/BCInterface.H"
 
namespace kynema_sgf {

 
class BCVOF : public BCIface
{
public:
    explicit BCVOF(Field& field) : BCIface(field) {}
 
protected:
    // Only the process of reading BCs needs to be modified
    void set_bcrec() override
    {
        const auto& ibctype = m_field.bc_type();
 
        for (amrex::OrientationIter oit; oit != nullptr; ++oit) {
            auto ori = oit();
            const auto side = ori.faceDir();
            const auto bct = ibctype[ori];
            const int dir = ori.coordDir();
 
            switch (bct) {
            case BC::periodic:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::int_dir);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::int_dir);
                }
                break;
 
            // For the wall cases, the fillpatch BC is effectively a
            // contact-line model
            case BC::pressure_outflow:
            case BC::zero_gradient:
            case BC::no_slip_wall:
            case BC::symmetric_wall:
            case BC::slip_wall:
            case BC::wall_model:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::foextrap);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::foextrap);
                }
                break;
 
                // Never do hoextrap because it is unbounded
 
            case BC::wave_generation:
            case BC::mass_inflow:
            case BC::mass_inflow_outflow:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::ext_dir);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::ext_dir);
                }
                break;
 
            // Fixed gradient BC is not allowed for VOF
            case BC::fixed_gradient:
            default:
                amrex::Abort(
                    "Invalid BC type for VOF encountered: fixed_gradient");
            }
        }
    }
    amrex::Array<const std::string, 3> inflow_udfs() override
    {
        return {"ConstDirichlet", "ConstDirichlet", "ConstDirichlet"};
    }
 
    void read_values() override { read_inflow_values(); }
};
 
class BCDensityMod : public BCIface
{
public:
    explicit BCDensityMod(Field& field) : BCIface(field) {}
 
protected:
    void set_bcrec() override
    {
        const auto& ibctype = m_field.bc_type();
 
        for (amrex::OrientationIter oit; oit != nullptr; ++oit) {
            auto ori = oit();
            const auto side = ori.faceDir();
            const auto bct = ibctype[ori];
            const int dir = ori.coordDir();
 
            switch (bct) {
            case BC::periodic:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::int_dir);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::int_dir);
                }
                break;
 
            // Both extrapolations are consistent between vof and density
            case BC::pressure_outflow:
            case BC::zero_gradient:
            case BC::no_slip_wall:
            case BC::symmetric_wall:
            case BC::slip_wall:
            case BC::wall_model:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::foextrap);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::foextrap);
                }
                break;
 
                // Never do hoextrap because it is unbounded
 
            case BC::wave_generation:
            case BC::mass_inflow:
            case BC::mass_inflow_outflow:
                if (side == amrex::Orientation::low) {
                    set_bcrec_lo(dir, amrex::BCType::ext_dir);
                } else {
                    set_bcrec_hi(dir, amrex::BCType::ext_dir);
                }
                break;
 
            case BC::fixed_gradient:
            default:
                amrex::Abort(
                    "Invalid BC type for density (vof-compatible) encountered: "
                    "fixed_gradient");
            }
        }
    }
 
    amrex::Array<const std::string, 3> inflow_udfs() override
    {
        return {"ConstDirichlet", "ConstDirichlet", "ConstDirichlet"};
    }
 
    void read_values() override { read_inflow_values(); }
};
 
} // namespace kynema_sgf
 
#endif /* VOF_BCIFACE_H */