Solver performing calculations in 2D Cartesian or Cylindrical space using finite element method. More...

#include <solvers/thermal/dynamic/femT2d.hpp>

Inheritance diagram for plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >:

Collaboration diagram for plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >:

Classes
struct	ThermalConductivityData

Public Member Functions
double	compute (double time)
	Run temperature calculations.

double	getElapsTime () const
	Get calculations elapsed time.

void	loadConfiguration (XMLReader &source, Manager &manager) override
	Load configuration from given `source`.

	DynamicThermalFem2DSolver (const std::string &name="")

std::string	getClassName () const override
	Get name of solver.

	~DynamicThermalFem2DSolver ()

std::string	getClassName () const
	Get name of solver.

std::string	getClassName () const
	Get name of solver.

Public Member Functions inherited from plask::FemSolverWithMaskedMesh< Geometry2DType, RectangularMesh< 2 > >
	FemSolverWithMaskedMesh (const std::string &name="")

EmptyElementsHandling	getEmptyElements () const
	Are we using full mesh?

void	setEmptyElements (EmptyElementsHandling val)
	Set whether we should use full mesh.

bool	parseFemConfiguration (XMLReader &reader, Manager &manager)

void	setupMaskedMesh ()

void	onInitialize ()
	Initialize the solver.

FemMatrix *	getMatrix ()

FemMatrix *	getMatrix ()

Public Member Functions inherited from plask::FemSolverWithMesh< SpaceT, MeshT >
	FemSolverWithMesh (const std::string &name="")

bool	parseFemConfiguration (XMLReader &reader, Manager &manager)

FemMatrix *	getMatrix ()

FemMatrix *	getMatrix ()

Public Member Functions inherited from plask::SolverWithMesh< SpaceT, MeshT >
	SolverWithMesh (const std::string &name="")

	~SolverWithMesh ()

void	loadConfiguration (XMLReader &source, Manager &manager) override
	Load configuration from given `source`.

void	parseStandardConfiguration (XMLReader &source, Manager &manager, const std::string &expected_msg="solver configuration element")

virtual void	onMeshChange (const typename MeshT::Event &PLASK_UNUSED(evt))
	This method is called just after the mesh has been changed.

void	onGeometryChange (const Geometry::Event &PLASK_UNUSED(evt)) override
	This method is called when the geometry is changed.

MeshT &	meshRef () const
	Get current module mesh.

shared_ptr< MeshT >	getMesh () const
	Get current solver mesh.

void	setMesh (const shared_ptr< MeshT > &mesh)
	Set new mesh for the solver.

void	setMesh (shared_ptr< MeshGeneratorD< MeshT::DIM > > generator)
	Set new mesh got from generator.

Public Member Functions inherited from plask::SolverOver< SpaceT >
	SolverOver (const std::string &name="")

	~SolverOver ()

void	parseStandardConfiguration (XMLReader &source, Manager &manager, const std::string &expected_msg="solver configuration element")

virtual void	onGeometryChange (const Geometry::Event &)
	This method is called when the geometry is changed.

shared_ptr< SpaceT >	getGeometry () const
	Get current solver geometry space.

void	setGeometry (const shared_ptr< SpaceT > &geometry)
	Set new geometry for the solver.

Public Member Functions inherited from plask::Solver
bool	initCalculation ()
	This should be called on beginning of each calculation method to ensure that solver will be initialized.

	Solver (const std::string &name="")
	Construct uninitialized solver.

virtual	~Solver ()
	Virtual destructor (for subclassing). Do nothing.

void	parseStandardConfiguration (XMLReader &source, Manager &manager, const std::string &expected_msg="solver configuration element")
	Load standard configuration (geometry, mesh) tags from `source`.

bool	isInitialized ()
	Check if solver is already initialized.

void	invalidate ()
	This method should be and is called if something important was changed: calculation space, mesh, etc.

std::string	getId () const
	Get solver id.

std::string	getName () const

virtual std::string	getClassDescription () const
	Get a description of this solver.

template<typename ArgT = double, typename ValT = double>
DataLog< ArgT, ValT >	dataLog (const std::string &chart_name, const std::string &axis_arg_name, const std::string &axis_val_name)

template<typename ArgT = double, typename ValT = double>
DataLog< ArgT, ValT >	dataLog (const std::string &axis_arg_name, const std::string &axis_val_name)

template<typename ... Args>
void	writelog (LogLevel level, std::string msg, Args &&... params) const
	Log a message for this solver.

Public Attributes
BoundaryConditions< RectangularMesh< 2 >::Boundary, double >	temperature_boundary
	Boundary condition of constant temperature (K)

ProviderFor< Temperature, Geometry2DType >::Delegate	outTemperature

ProviderFor< HeatFlux, Geometry2DType >::Delegate	outHeatFlux

ProviderFor< ThermalConductivity, Geometry2DType >::Delegate	outThermalConductivity

ReceiverFor< Heat, Geometry2DType >	inHeat

double	inittemp
	Initial temperature.

double	methodparam
	Initial parameter determining the calculation method (0.5 - Crank-Nicolson, 0 - explicit, 1 - implicit)

double	timestep
	Time step in nanoseconds.

double	elapstime
	Calculations elapsed time.

bool	lumping
	Whether use lumping for matrices?

size_t	rebuildfreq
	Frequency of mass matrix rebuilding.

size_t	logfreq
	Frequency of iteration progress reporting.

Public Attributes inherited from plask::FemSolverWithMesh< SpaceT, MeshT >
FemMatrixAlgorithm	algorithm = ALGORITHM_CHOLESKY
	Factorization algorithm to use.

IterativeMatrixParams	iter_params
	Parameters of iterative solver.

Protected Member Functions
void	setMatrix (FemMatrix &A, FemMatrix &B, DataVector< double > &F, const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &btemperature)
	Set stiffness matrix + load vector.

void	saveHeatFluxes ()
	Create 2D-vector with calculated heat fluxes.

void	onInitialize () override
	Initialize the solver.

void	onInvalidate () override
	Invalidate the data.

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 > > &dst_mesh, InterpolationMethod method)

void	setMatrix (FemMatrix &A, FemMatrix &B, DataVector< double > &F, const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &btemperature)

void	setMatrix (FemMatrix &A, FemMatrix &B, DataVector< double > &F, const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &btemperature)

Protected Member Functions inherited from plask::SolverOver< SpaceT >
template<typename Boundary , typename ConditionT >
void	readBoundaryConditions (Manager &manager, XMLReader &reader, BoundaryConditions< Boundary, ConditionT > &dest)
	Read boundary conditions using information about the geometry of this solver.

Protected Attributes
double	maxT
	Maximum temperature recorded.

DataVector< double >	temperatures
	Computed temperatures.

DataVector< double >	thickness
	Thicknesses of the layers.

DataVector< Vec< 2, double > >	fluxes
	Computed (only when needed) heat fluxes on our own mesh.

Protected Attributes inherited from plask::FemSolverWithMaskedMesh< Geometry2DType, RectangularMesh< 2 > >
plask::shared_ptr< RectangularMaskedMesh< MeshT::DIM > >	maskedMesh

EmptyElementsHandling	empty_elements
	Should we use full mesh?

Protected Attributes inherited from plask::SolverWithMesh< SpaceT, MeshT >
shared_ptr< MeshT >	mesh
	Mesh over which the calculations are performed.

boost::signals2::connection	mesh_signal_connection
	Connection of mesh to onMeshChange method, see http://www.boost.org/doc/libs/1_55_0/doc/html/signals2/tutorial.html#idp204830936.

Protected Attributes inherited from plask::SolverOver< SpaceT >
shared_ptr< SpaceT >	geometry
	Space in which the calculations are performed.

Protected Attributes inherited from plask::Solver
bool	initialized
	`true` only if solver is initialized

Additional Inherited Members
Public Types inherited from plask::SolverWithMesh< SpaceT, MeshT >
typedef MeshT	MeshType
	Type of the mesh for this solver.

Public Types inherited from plask::SolverOver< SpaceT >
typedef SpaceT	SpaceType
	of the space for this solver

Detailed Description

template<typename Geometry2DType>
struct plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >

Solver performing calculations in 2D Cartesian or Cylindrical space using finite element method.

Definition at line 27 of file femT2d.hpp.

Constructor & Destructor Documentation

◆ DynamicThermalFem2DSolver()

template<typename Geometry2DType >

plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::DynamicThermalFem2DSolver ( const std::string & name = "" )

Definition at line 21 of file femT2d.cpp.

◆ ~DynamicThermalFem2DSolver()

template<typename Geometry2DType >

plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::~DynamicThermalFem2DSolver ( )

Definition at line 41 of file femT2d.cpp.

Member Function Documentation

◆ compute()

template<typename Geometry2DType >

double plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::compute ( double time )

Run temperature calculations.

Returns: max correction of temperature against the last call

Definition at line 388 of file femT2d.cpp.

◆ getClassName() [1/3]

std::string plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DCartesian >::getClassName ( ) const

virtual

Get name of solver.

Returns: name of this solver

Implements plask::Solver.

Definition at line 558 of file femT2d.cpp.

◆ getClassName() [2/3]

std::string plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DCylindrical >::getClassName ( ) const

virtual

Get name of solver.

Returns: name of this solver

Implements plask::Solver.

Definition at line 559 of file femT2d.cpp.

◆ getClassName() [3/3]

template<typename Geometry2DType >

std::string plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::getClassName ( ) const

overridevirtual

Get name of solver.

Returns: name of this solver

Implements plask::Solver.

◆ getElapsTime()

template<typename Geometry2DType >

double plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::getElapsTime ( ) const

inline

Get calculations elapsed time.

Definition at line 81 of file femT2d.hpp.

◆ getHeatFluxes()

template<typename Geometry2DType >

const LazyData< Vec< 2 > > plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::getHeatFluxes	(	const shared_ptr< const MeshD< 2 > > &	dest_mesh,
		InterpolationMethod	method
	)

protected

Definition at line 505 of file femT2d.cpp.

◆ getTemperatures()

template<typename Geometry2DType >

const LazyData< double > plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::getTemperatures	(	const shared_ptr< const MeshD< 2 > > &	dest_mesh,
		InterpolationMethod	method
	)		const

protected

Definition at line 493 of file femT2d.cpp.

◆ getThermalConductivity()

template<typename Geometry2DType >

const LazyData< Tensor2< double > > plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::getThermalConductivity	(	const shared_ptr< const MeshD< 2 > > &	dst_mesh,
		InterpolationMethod	method
	)

protected

Definition at line 549 of file femT2d.cpp.

◆ loadConfiguration()

template<typename Geometry2DType >

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::loadConfiguration	(	XMLReader &	source,
		Manager &	manager
	)

overridevirtual

Load configuration from given source.

XML reader (source) point to opening of this solver tag and after return from this method should point to this solver closing tag.

Parameters

source	source of configuration
manager	manager from which information about geometry, meshes, materials, and so on can be get if needed

Reimplemented from plask::Solver.

Definition at line 46 of file femT2d.cpp.

◆ onInitialize()

template<typename Geometry2DType >

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::onInitialize ( )

overrideprotectedvirtual

Initialize the solver.

Reimplemented from plask::Solver.

Definition at line 75 of file femT2d.cpp.

◆ onInvalidate()

template<typename Geometry2DType >

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::onInvalidate ( )

overrideprotectedvirtual

Invalidate the data.

Reimplemented from plask::Solver.

Definition at line 122 of file femT2d.cpp.

◆ saveHeatFluxes()

template<typename Geometry2DType >

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::saveHeatFluxes ( )

protected

Create 2D-vector with calculated heat fluxes.

Definition at line 454 of file femT2d.cpp.

◆ setMatrix() [1/3]

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DCartesian >::setMatrix	(	FemMatrix &	A,
		FemMatrix &	B,
		DataVector< double > &	F,
		const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &	btemperature
	)

protected

Definition at line 129 of file femT2d.cpp.

◆ setMatrix() [2/3]

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DCylindrical >::setMatrix	(	FemMatrix &	A,
		FemMatrix &	B,
		DataVector< double > &	F,
		const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &	btemperature
	)

protected

Definition at line 259 of file femT2d.cpp.

◆ setMatrix() [3/3]

template<typename Geometry2DType >

void plask::thermal::dynamic::DynamicThermalFem2DSolver< Geometry2DType >::setMatrix	(	FemMatrix &	A,
		FemMatrix &	B,
		DataVector< double > &	F,
		const BoundaryConditionsWithMesh< RectangularMesh< 2 >::Boundary, double > &	btemperature
	)