therm2d.hpp

Go to the documentation of this file.

/*
 * This file is part of PLaSK (https://plask.app) by Photonics Group at TUL
 * Copyright (c) 2022 Lodz University of Technology
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, version 3.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 */
#ifndef PLASK__MODULE_THERMAL_THERM2D_H
#define PLASK__MODULE_THERMAL_THERM2D_H
 
#include <plask/plask.hpp>
#include <plask/common/fem.hpp>
 
#include "common.hpp"
 
namespace plask { namespace thermal { namespace tstatic {
 
template <typename Geometry2DType>
struct PLASK_SOLVER_API ThermalFem2DSolver : public FemSolverWithMaskedMesh<Geometry2DType, RectangularMesh<2>> {
  protected:
    int loopno;     
    double maxT;    
    double toterr;  
 
    DataVector<double> temperatures;  
 
    DataVector<double> thickness;  
 
    DataVector<Vec<2, double>> fluxes;  
 
    void setMatrix(FemMatrix& A,
                   DataVector<double>& B,
                   const BoundaryConditionsWithMesh<RectangularMesh<2>::Boundary, double>& btemperature,
                   const BoundaryConditionsWithMesh<RectangularMesh<2>::Boundary, double>& bheatflux,
                   const BoundaryConditionsWithMesh<RectangularMesh<2>::Boundary, Convection>& bconvection,
                   const BoundaryConditionsWithMesh<RectangularMesh<2>::Boundary, Radiation>& bradiation);
 
    template <typename MatrixT> MatrixT makeMatrix();
 
    void saveHeatFluxes();  // [W/m^2]
 
    void onInitialize() override;
 
    void onInvalidate() override;
 
  public:
    // Boundary conditions
    BoundaryConditions<RectangularMesh<2>::Boundary, double>
        temperature_boundary;  
    BoundaryConditions<RectangularMesh<2>::Boundary, double>
        heatflux_boundary;  
    BoundaryConditions<RectangularMesh<2>::Boundary, Convection> convection_boundary;  
    BoundaryConditions<RectangularMesh<2>::Boundary, Radiation> radiation_boundary;    
 
    typename ProviderFor<Temperature, Geometry2DType>::Delegate outTemperature;
 
    typename ProviderFor<HeatFlux, Geometry2DType>::Delegate outHeatFlux;
 
    typename ProviderFor<ThermalConductivity, Geometry2DType>::Delegate outThermalConductivity;
 
    ReceiverFor<Heat, Geometry2DType> inHeat;
 
    double maxerr;    
    double inittemp;  
 
    double compute(int loops = 1);
 
    double getErr() const { return toterr; }
 
    void loadConfiguration(XMLReader& source,
                           Manager& manager) override;  // for solver configuration (see: *.xpl file with structures)
 
    ThermalFem2DSolver(const std::string& name = "");
 
    std::string getClassName() const override;
 
    ~ThermalFem2DSolver();
 
  public:
    struct ThermalConductivityData : public LazyDataImpl<Tensor2<double>> {
        const ThermalFem2DSolver* solver;
        shared_ptr<const MeshD<2>> dest_mesh;
        InterpolationFlags flags;
        LazyData<double> temps;
        ThermalConductivityData(const ThermalFem2DSolver* solver, const shared_ptr<const MeshD<2>>& dst_mesh);
        Tensor2<double> at(std::size_t i) const override;
        std::size_t size() const override;
    };
 
    const LazyData<double> getTemperatures(const shared_ptr<const MeshD<2>>& dest_mesh, InterpolationMethod method) const;
 
    const LazyData<Vec<2>> getHeatFluxes(const shared_ptr<const MeshD<2>>& dest_mesh, InterpolationMethod method);
 
    const LazyData<Tensor2<double>> getThermalConductivity(const shared_ptr<const MeshD<2>>& dest_mesh, InterpolationMethod method);
};
 
}}}  // namespace plask::thermal::tstatic
 
#endif