BCConstantHeatFlux.h

Go to the documentation of this file.

#ifndef BOUNDARY_CONDITION_CONSTANT_HEAT_FLUX_H
#define BOUNDARY_CONDITION_CONSTANT_HEAT_FLUX_H
 
#include <cmath>
#include <memory>
#include <numbers>
#include <string>
#include <valarray>
#include <vector>
 
#include "GenericBoundaryCondition.h"
#include "IniParser.h"
#include "SimulationClassBase.h"
 
template <typename T>
class BCConstantHeatFlux : public GenericBoundaryCondition<T>
{
    private:
        T m_heat_flux{};
        T m_area{};
 
    public:
        void writeBoundaryCondition() override;
 
        BCConstantHeatFlux(SimulationClassBase<T>* sim, std::vector<int> boundary_points,
                           std::vector<std::vector<int>> connected_points, int facet_nr, T area,
                           std::string face_prefix);
        ~BCConstantHeatFlux() override = default;
};
 
template <typename T>
BCConstantHeatFlux<T>::BCConstantHeatFlux(SimulationClassBase<T>* sim,
                                          std::vector<int> boundary_points,
                                          std::vector<std::vector<int>> connected_points,
                                          int facet_nr, T area, std::string face_prefix)
    : GenericBoundaryCondition<T>{sim, boundary_points, connected_points, facet_nr, face_prefix},
      m_area{area}
{
    // Only const heat flux case right now
    try
    {
        m_heat_flux = this->m_sim->m_simulation_config.getDoubleParameters(
            this->m_face_indicator_prefix + "boundary_heat_flux");
    }
    catch (const BadInput& e)
    {
        std::cerr
            << "[BCConstantHeatFlux]: Couldn't find heat flux value in simulation config under: "
            << this->m_face_indicator_prefix << "boundary_heat_flux.\n";
        throw BadInput(
            static_cast<std::string>("[BCConstantHeatFlux]: Couldn't find heat flux value "
                                     "in simulation config under: ")
            + this->m_face_indicator_prefix + static_cast<std::string>("boundary_heat_flux."));
    }
}
 
template <typename T>
void BCConstantHeatFlux<T>::writeBoundaryCondition()
{
    for (int i{0}; i < this->m_boundary_points.size(); i++)
    {
        this->m_forcing_vector_elements[i] += m_heat_flux * m_area;
    }
}
 
#endif