ERF_EBStruct.H

Go to the documentation of this file.

#ifndef ERF_EB_STRUCT_H_
#define ERF_EB_STRUCT_H_
 
#include <string>
#include <iostream>
 
#include <AMReX_ParmParse.H>
#include <AMReX_Print.H>
#include <AMReX_Gpu.H>
 
/**
 * @brief Boundary condition model applied on embedded-boundary surfaces.
 */
enum struct EBBoundaryType {
    SlipWall, NoSlipWall, SurfaceLayer
};
 
/**
 * @brief Embedded-boundary stress component indices.
 */
enum EBTauType {
    tau_eb13, tau_eb23
};
 
/**
 * @brief Grid centering used for embedded-boundary data.
 */
enum EBGridType {
    xface, yface, zface
};
 
/**
 * Container holding EB-related choices
 */
 
struct EBChoice {
  public:
    /**
     * @brief Read embedded-boundary options from the input parameter database.
     * @param pp_prefix ParmParse prefix for the ERF input namespace.
     */
    void init_params(std::string pp_prefix)
    {
        amrex::ParmParse pp(pp_prefix);
 
        static std::string eb_boundary_type_string = "NoSlipWall";
        pp.query("eb_boundary_type",eb_boundary_type_string);
 
        if (!eb_boundary_type_string.compare("SlipWall")) {
            eb_boundary_type = EBBoundaryType::SlipWall;
        } else if (!eb_boundary_type_string.compare("NoSlipWall")) {
            eb_boundary_type = EBBoundaryType::NoSlipWall;
        } else if (!eb_boundary_type_string.compare("SurfaceLayer")) {
            eb_boundary_type = EBBoundaryType::SurfaceLayer;
        } else {
            amrex::Error("Don't know this eb_boundary_type");
        }
    }
 
    /**
     * @brief Print the configured embedded-boundary options.
     */
    void display()
    {
        amrex::Print() << "EB choices: " << std::endl;
        if (eb_boundary_type == EBBoundaryType::SlipWall) {
            amrex::Print() << "    eb_boundary_type        : SlipWall" << std::endl;
        } else if (eb_boundary_type == EBBoundaryType::NoSlipWall) {
            amrex::Print() << "    eb_boundary_type        : NoSlipWall" << std::endl;
        } else if (eb_boundary_type == EBBoundaryType::SurfaceLayer) {
            amrex::Print() << "    eb_boundary_type        : SurfaceLayer" << std::endl;
        }
    }
 
    EBBoundaryType eb_boundary_type = EBBoundaryType::SlipWall;
};
 
/**
 * @var EBChoice::eb_boundary_type
 * @brief Boundary condition model applied on embedded-boundary surfaces.
 */
#endif