Coupling of a transport model with a reaction model by operator splitting. More...

#include <transport_operator_splitting.hh>

Inheritance diagram for TransportOperatorSplitting:

Collaboration diagram for TransportOperatorSplitting:

Public Types
typedef AdvectionProcessBase	FactoryBaseType

Public Member Functions
	TransportOperatorSplitting (Mesh &init_mesh, const Input::Record in_rec)
	Constructor. More...

virtual	~TransportOperatorSplitting ()
	Destructor. More...

virtual void	set_velocity_field (const MH_DofHandler &dh) override

void	initialize () override

void	zero_time_step () override

void	update_solution () override

void	compute_until_save_time ()

void	compute_internal_step ()

void	output_data () override
	Write computed fields. More...

Public Member Functions inherited from AdvectionProcessBase
	AdvectionProcessBase (Mesh &mesh, const Input::Record in_rec)

Public Member Functions inherited from EquationBase
	EquationBase ()

	EquationBase (Mesh &mesh, const Input::Record in_rec)

virtual	~EquationBase ()

virtual void	choose_next_time ()

virtual void	set_time_upper_constraint (double dt, std::string message)

virtual void	set_time_lower_constraint (double dt, std::string message)

TimeGovernor &	time ()

virtual void	set_time_governor (TimeGovernor &time)

double	planned_time ()

double	solved_time ()

Mesh &	mesh ()

TimeMark::Type	mark_type ()

FieldSet &	data ()

virtual void	get_solution_vector (double *&vector, unsigned int &size)

virtual void	get_parallel_solution_vector (Vec &vector)

Static Public Member Functions
static const Input::Type::Record &	get_input_type ()
	Declare input record type for the equation TransportOperatorSplittiong. More...

Static Public Member Functions inherited from AdvectionProcessBase
static Input::Type::Abstract &	get_input_type ()
	Common specification of the input record for secondary equations. More...

Private Attributes
std::shared_ptr< ConcentrationTransportBase >	convection

std::shared_ptr< ReactionTerm >	reaction

double	cfl_convection
	Time restriction due to transport. More...

double	cfl_reaction
	Time restriction due to reactions. More...

Static Private Attributes
static const int	registrar
	Registrar of class to factory. More...

Additional Inherited Members
Protected Attributes inherited from EquationBase
bool	equation_empty_
	flag is true if only default constructor was called More...

Mesh *	mesh_

TimeGovernor *	time_

Input::Record	input_record_

FieldSet *	eq_data_

std::shared_ptr< Balance >	balance_
	object for calculation and writing the mass balance to file. More...

Detailed Description

Coupling of a transport model with a reaction model by operator splitting.

Outline: Transport model is any descendant of TransportBase (even TransportOperatorSplitting itself). This should perform the transport possibly with diffusion and usually without coupling between substances and phases.

Reaction is any descendant of the ReactionBase class. This represents reactions in general way of any coupling that happens between substances and phases on one element or more generally on one DoF.

Definition at line 202 of file transport_operator_splitting.hh.

Member Typedef Documentation

typedef AdvectionProcessBase TransportOperatorSplitting::FactoryBaseType

Definition at line 204 of file transport_operator_splitting.hh.

Constructor & Destructor Documentation

TransportOperatorSplitting::TransportOperatorSplitting	(	Mesh &	init_mesh,
		const Input::Record	in_rec
	)

Constructor.

Definition at line 146 of file transport_operator_splitting.cc.

TransportOperatorSplitting::~TransportOperatorSplitting ( )

virtual

Destructor.

Definition at line 208 of file transport_operator_splitting.cc.

Member Function Documentation

void TransportOperatorSplitting::compute_internal_step ( )

void TransportOperatorSplitting::compute_until_save_time ( )

const Record & TransportOperatorSplitting::get_input_type ( )

static

Declare input record type for the equation TransportOperatorSplittiong.

TODO: The question is if this should be a general coupling class (e.g. allow coupling TranportDG with reactions even if it is not good idea for numerical reasons.) To make this a coupling class we should modify all main input files for transport problems.

Definition at line 70 of file transport_operator_splitting.cc.

void TransportOperatorSplitting::initialize ( )

overridevirtual

This method should initialize fields of the equation. All members (e.g. number of components) that are necessary for the field initialization must be set between construction and call of initialize. After this method the upper level coupling may set sharing of some fields between equations.

Reimplemented from EquationBase.

Definition at line 216 of file transport_operator_splitting.cc.

void TransportOperatorSplitting::output_data ( void )

overridevirtual

Write computed fields.

Reimplemented from EquationBase.

Definition at line 231 of file transport_operator_splitting.cc.

void TransportOperatorSplitting::set_velocity_field ( const MH_DofHandler & dh )

overridevirtual

This method takes sequential PETSc vector of side velocities and update transport matrix. The ordering is same as ordering of sides in the mesh. We just keep the pointer, but do not destroy the object.

TODO: We should pass whole velocity field object (description of base functions and dof numbering) and vector.

Implements AdvectionProcessBase.

Definition at line 349 of file transport_operator_splitting.cc.

void TransportOperatorSplitting::update_solution ( void )

overridevirtual

Calculation of the next time step and its output.

Reimplemented from EquationBase.

Definition at line 260 of file transport_operator_splitting.cc.

void TransportOperatorSplitting::zero_time_step ( )

overridevirtual

Initialization of the solution in the zero time.

There may be fields that can not be initialized in the initialize method as they are provided by the coupling. Fields coming from coupling has to be set after the initialize method and before zero_time_step.

Reimplemented from EquationBase.

Definition at line 245 of file transport_operator_splitting.cc.

Member Data Documentation

double TransportOperatorSplitting::cfl_convection

private

Time restriction due to transport.

Definition at line 242 of file transport_operator_splitting.hh.

double TransportOperatorSplitting::cfl_reaction

private

Time restriction due to reactions.

Definition at line 243 of file transport_operator_splitting.hh.

std::shared_ptr<ConcentrationTransportBase> TransportOperatorSplitting::convection

private

Definition at line 236 of file transport_operator_splitting.hh.

std::shared_ptr<ReactionTerm> TransportOperatorSplitting::reaction

private

Definition at line 237 of file transport_operator_splitting.hh.

const int TransportOperatorSplitting::registrar

staticprivate

Initial value:

=
        Input::register_class< TransportOperatorSplitting, Mesh &, const Input::Record>("Coupling_OperatorSplitting") +
        TransportOperatorSplitting::get_input_type().size()

Registrar of class to factory.

Definition at line 234 of file transport_operator_splitting.hh.

The documentation for this class was generated from the following files:

/opt/flow123d/flow123d/src/transport/transport_operator_splitting.hh
/opt/flow123d/flow123d/src/transport/transport_operator_splitting.cc

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation