therm3d.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_TFem_H
#define PLASK__MODULE_THERMAL_TFem_H
 
#include <plask/plask.hpp>
#include <plask/common/fem.hpp>
 
#include "common.hpp"
 
namespace plask { namespace thermal { namespace tstatic {
 
struct PLASK_SOLVER_API ThermalFem3DSolver: public FemSolverWithMaskedMesh<Geometry3D, RectangularMesh<3>> {
 
  protected:
 
    int loopno;         
    double maxT;        
    double toterr;      
 
    DataVector<double> temperatures;            
 
    DataVector<double> thickness;               
 
    DataVector<Vec<3,double>> fluxes;           
 
    void setMatrix(FemMatrix& A, DataVector<double>& B,
                   const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,double>& btemperature,
                   const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,double>& bheatflux,
                   const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,Convection>& bconvection,
                   const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,Radiation>& bradiation
                  );
 
    template <typename MatrixT>
    void applyBC(MatrixT& A, DataVector<double>& B, const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,double>& btemperature);
 
    template <typename MatrixT>
    MatrixT makeMatrix();
 
    void saveHeatFluxes(); // [W/m^2]
 
    void onInitialize() override;
 
    void onInvalidate() override;
 
  public:
 
    double inittemp;    
    double maxerr;     
 
    // Boundary conditions
    BoundaryConditions<RectangularMesh<3>::Boundary,double> temperature_boundary;      
    BoundaryConditions<RectangularMesh<3>::Boundary,double> heatflux_boundary;         
    BoundaryConditions<RectangularMesh<3>::Boundary,Convection> convection_boundary;   
    BoundaryConditions<RectangularMesh<3>::Boundary,Radiation> radiation_boundary;     
 
    typename ProviderFor<Temperature,Geometry3D>::Delegate outTemperature;
 
    typename ProviderFor<HeatFlux,Geometry3D>::Delegate outHeatFlux;
 
    typename ProviderFor<ThermalConductivity,Geometry3D>::Delegate outThermalConductivity;
 
    ReceiverFor<Heat,Geometry3D> inHeat;
 
    double compute(int loops=1);
 
    void loadConfiguration(XMLReader& source, Manager& manager) override; // for solver configuration (see: *.xpl file with structures)
 
    ThermalFem3DSolver(const std::string& name="");
 
    std::string getClassName() const override { return "thermal.Static3D"; }
 
    ~ThermalFem3DSolver();
 
    double getErr() const { return toterr; }
 
    FemMatrixAlgorithm getAlgorithm() const { return algorithm; }
 
    void setAlgorithm(FemMatrixAlgorithm alg);
 
  protected:
 
    struct ThermalConductivityData: public LazyDataImpl<Tensor2<double>> {
        const ThermalFem3DSolver* solver;
        shared_ptr<const MeshD<3>> dest_mesh;
        InterpolationFlags flags;
        LazyData<double> temps;
        ThermalConductivityData(const ThermalFem3DSolver* solver, const shared_ptr<const MeshD<3>>& 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<3>>& dst_mesh, InterpolationMethod method) const;
 
    const LazyData<Vec<3>> getHeatFluxes(const shared_ptr<const MeshD<3>>& dst_mesh, InterpolationMethod method);
 
    const LazyData<Tensor2<double>> getThermalConductivity(const shared_ptr<const MeshD<3>>& dst_mesh, InterpolationMethod method);
};
 
}}} //namespaces
 
#endif