femT3d.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_TFem3D_H
#define PLASK__MODULE_THERMAL_TFem3D_H
 
#include <plask/plask.hpp>
#include <plask/common/fem.hpp>
 
namespace plask { namespace thermal { namespace dynamic {
 
struct PLASK_SOLVER_API DynamicThermalFem3DSolver: public FemSolverWithMaskedMesh<Geometry3D, RectangularMesh<3>> {
 
  protected:
 
    double maxT;        
 
    DataVector<double> temperatures;           
 
    DataVector<double> thickness;               
 
    DataVector<Vec<3,double>> fluxes;      
 
    void setMatrix(FemMatrix& A, FemMatrix& B, DataVector<double>& F,
                   const BoundaryConditionsWithMesh<RectangularMesh<3>::Boundary,double>& btemperature
                  );
 
    void saveHeatFluxes(); // [W/m^2]
 
    void onInitialize() override;
 
    void onInvalidate() override;
 
  public:
 
    // Boundary conditions
    BoundaryConditions<RectangularMesh<3>::Boundary,double> temperature_boundary;      
 
    typename ProviderFor<Temperature, Geometry3D>::Delegate outTemperature;
 
    typename ProviderFor<HeatFlux, Geometry3D>::Delegate outHeatFlux;
 
    typename ProviderFor<ThermalConductivity, Geometry3D>::Delegate outThermalConductivity;
 
    ReceiverFor<Heat, Geometry3D> inHeat;
 
    double inittemp;       
    double methodparam;   
    double timestep;       
    double elapstime;    
    bool lumping;          
    size_t rebuildfreq;    
    size_t logfreq;        
 
    double compute(double time);
 
    double getElapsTime() const { return elapstime; }
 
    void loadConfiguration(XMLReader& source, Manager& manager) override; // for solver configuration (see: *.xpl file with structures)
 
    DynamicThermalFem3DSolver(const std::string& name="");
 
    std::string getClassName() const override { return "thermal.Dynamic3D"; }
 
    ~DynamicThermalFem3DSolver();
 
  protected:
 
    struct ThermalConductivityData: public LazyDataImpl<Tensor2<double>> {
        const DynamicThermalFem3DSolver* solver;
        shared_ptr<const MeshD<3>> dest_mesh;
        InterpolationFlags flags;
        LazyData<double> temps;
        ThermalConductivityData(const DynamicThermalFem3DSolver* 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
 
} // namespace plask
 
#endif