#include <transport_dg.hh>

Inheritance diagram for TransportDG< Model >::EqData:

Collaboration diagram for TransportDG< Model >::EqData:

[legend]

Public Member Functions
	EqData ()

void	set_DG_parameters_boundary (Side side, const int K_size, const std::vector< arma::mat33 > &K, const double flux, const arma::vec3 &normal_vector, const double alpha, double &gamma)
	Sets up parameters of the DG method on a given boundary edge. More...

double	elem_anisotropy (ElementAccessor< 3 > e) const
	Compute and return anisotropy of given element. More...

Public Attributes
MultiField< 3, FieldValue< 3 >::Scalar >	fracture_sigma
	Transition parameter for diffusive transfer on fractures (for each substance). More...

MultiField< 3, FieldValue< 3 >::Scalar >	dg_penalty
	Penalty enforcing inter-element continuity of solution (for each substance). More...

Field< 3, FieldValue< 3 >::Scalar >	region_id

Field< 3, FieldValue< 3 >::Scalar >	subdomain

EquationOutput	output_fields

std::vector< Vec >	ret_vec
	Auxiliary vectors for calculation of sources in balance due to retardation (e.g. sorption). More...

LinSys **	ls
	Linear algebra system for the transport equation. More...

LinSys **	ls_dt
	Linear algebra system for the time derivative (actually it is used only for handling the matrix structures). More...

std::shared_ptr< DOFHandlerMultiDim >	dh_
	Object for distribution of dofs. More...

GenericAssembly< MassAssemblyDim > *	mass_assembly_
	general assembly objects, hold assembly objects of appropriate dimension More...

GenericAssembly< StiffnessAssemblyDim > *	stiffness_assembly_

GenericAssembly< SourcesAssemblyDim > *	sources_assembly_

GenericAssembly< BdrConditionAssemblyDim > *	bdr_cond_assembly_

GenericAssembly< InitConditionAssemblyDim > *	init_cond_assembly_

FieldFEScalarVec	conc_fe

std::shared_ptr< DOFHandlerMultiDim >	dh_p0

Parameters of the numerical method
std::vector< std::vector< double > >	gamma
	Penalty parameters. More...

int	dg_variant
	DG variant ((non-)symmetric/incomplete. More...

unsigned int	dg_order
	Polynomial order of finite elements. More...

Auxiliary fields used during assembly
vector< vector< arma::vec3 > >	ad_coef
	Advection coefficients. More...

vector< vector< arma::mat33 > >	dif_coef
	Diffusion coefficients. More...

vector< vector< vector< arma::vec3 > > >	ad_coef_edg
	Advection coefficients on edges. More...

vector< vector< vector< arma::mat33 > > >	dif_coef_edg
	Diffusion coefficients on edges. More...

Detailed Description

template<class Model>
class TransportDG< Model >::EqData

Definition at line 144 of file transport_dg.hh.

Constructor & Destructor Documentation

template<class Model >

TransportDG< Model >::EqData::EqData ( )

Definition at line 107 of file transport_dg.cc.

Member Function Documentation

template<class Model >

double TransportDG< Model >::EqData::elem_anisotropy ( ElementAccessor< 3 > e ) const

Compute and return anisotropy of given element.

Definition at line 146 of file transport_dg.cc.

Here is the caller graph for this function:

template<class Model >

void TransportDG< Model >::EqData::set_DG_parameters_boundary	(	Side	side,
		const int	K_size,
		const std::vector< arma::mat33 > &	K,
		const double	flux,
		const arma::vec3 &	normal_vector,
		const double	alpha,
		double &	gamma
	)

Sets up parameters of the DG method on a given boundary edge.

Assumption is that the edge consists of only 1 side.

Parameters

side	The boundary side.
K_size	Size of vector of tensors K.
K	Dispersivity tensor.
ad_vector	Advection vector.
normal_vector	Normal vector (assumed constant along the edge).
alpha	Penalty parameter that influences the continuity of the solution (large value=more continuity).
gamma	Computed penalty parameters.