Flow123d  JS_before_hm-62-ge56f9d5
Classes | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
HM_Iterative Class Reference

Class for solution of fully coupled flow and mechanics using fixed-stress iterative splitting. More...

#include <hm_iterative.hh>

Inheritance diagram for HM_Iterative:
Inheritance graph
[legend]
Collaboration diagram for HM_Iterative:
Collaboration graph
[legend]

Classes

class  EqData
 

Public Member Functions

 HM_Iterative (Mesh &mesh, Input::Record in_record)
 
void initialize () override
 
void zero_time_step () override
 
void update_solution () override
 
const MH_DofHandlerget_mh_dofhandler () override
 
 ~HM_Iterative ()
 
- Public Member Functions inherited from DarcyFlowInterface
 DarcyFlowInterface (Mesh &mesh, const Input::Record in_rec)
 
virtual ~DarcyFlowInterface ()
 
- 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)
 
TimeGovernortime ()
 
virtual void set_time_governor (TimeGovernor &time)
 
double planned_time ()
 
double solved_time ()
 
Meshmesh ()
 
TimeMark::Type mark_type ()
 
FieldSetdata ()
 
virtual void get_solution_vector (FMT_UNUSED double *&vector, FMT_UNUSED unsigned int &size)
 
virtual void get_parallel_solution_vector (FMT_UNUSED Vec &vector)
 
virtual void output_data ()
 Write computed fields. More...
 

Static Public Member Functions

static const Input::Type::Recordget_input_type ()
 Define input record. More...
 
- Static Public Member Functions inherited from DarcyFlowInterface
static Input::Type::Abstractget_input_type ()
 

Private Member Functions

void update_potential ()
 
void update_flow_fields ()
 
void compute_iteration_error (double &difference, double &norm)
 

Private Attributes

std::shared_ptr< RichardsLMHflow_
 steady or unsteady water flow simulator based on MH scheme More...
 
std::shared_ptr< Elasticitymechanics_
 solute transport with chemistry through operator splitting More...
 
EqData data_
 
double beta_
 Tuning parameter for iterative splitting. More...
 
unsigned int min_it_
 Minimal number of iterations to perform. More...
 
unsigned int max_it_
 Maximal number of iterations. More...
 
double a_tol_
 Absolute tolerance for difference between two succeeding iterations. More...
 
double r_tol_
 Relative tolerance for difference between two succeeding iterations. More...
 

Static Private Attributes

static const int registrar
 

Additional Inherited Members

- Public Types inherited from DarcyFlowInterface
typedef DarcyFlowInterface FactoryBaseType
 Typedef for usage of Input::Factory in child classes. More...
 
- Protected Attributes inherited from EquationBase
bool equation_empty_
 flag is true if only default constructor was called More...
 
Meshmesh_
 
TimeGovernortime_
 
Input::Record input_record_
 
FieldSeteq_data_
 
std::shared_ptr< Balancebalance_
 object for calculation and writing the mass balance to file. More...
 

Detailed Description

Class for solution of fully coupled flow and mechanics using fixed-stress iterative splitting.

Flow and mechanics are solved separately and within each iteration the coupling terms are updated. Here we use the fixed-stress splitting [see Mikelic&Wheeler, Comput. Geosci. 17(3), 2013] which uses a tuning parameter "beta" to speed up the convergence.

TODO: The class is currently inherited from DarcyFlowInterface in order to provide MH_DofHandler for transport processes. This should be changed as soon as we replace MH_DofHandler by fields for velocity and pressure.

Definition at line 47 of file hm_iterative.hh.

Constructor & Destructor Documentation

HM_Iterative::HM_Iterative ( Mesh mesh,
Input::Record  in_record 
)

Definition at line 184 of file hm_iterative.cc.

HM_Iterative::~HM_Iterative ( )

Definition at line 465 of file hm_iterative.cc.

Member Function Documentation

void HM_Iterative::compute_iteration_error ( double &  difference,
double &  norm 
)
private

Definition at line 437 of file hm_iterative.cc.

Here is the caller graph for this function:

const it::Record & HM_Iterative::get_input_type ( )
static

Define input record.

Definition at line 96 of file hm_iterative.cc.

const MH_DofHandler & HM_Iterative::get_mh_dofhandler ( )
overridevirtual

Implements DarcyFlowInterface.

Definition at line 459 of file hm_iterative.cc.

void HM_Iterative::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 224 of file hm_iterative.cc.

void HM_Iterative::update_flow_fields ( )
private

Definition at line 371 of file hm_iterative.cc.

Here is the caller graph for this function:

void HM_Iterative::update_potential ( )
private

Definition at line 322 of file hm_iterative.cc.

Here is the caller graph for this function:

void HM_Iterative::update_solution ( void  )
overridevirtual

Calculation of the next time step and its output.

Reimplemented from EquationBase.

Definition at line 276 of file hm_iterative.cc.

void HM_Iterative::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 259 of file hm_iterative.cc.

Member Data Documentation

double HM_Iterative::a_tol_
private

Absolute tolerance for difference between two succeeding iterations.

Definition at line 114 of file hm_iterative.hh.

double HM_Iterative::beta_
private

Tuning parameter for iterative splitting.

Definition at line 105 of file hm_iterative.hh.

EqData HM_Iterative::data_
private

Definition at line 102 of file hm_iterative.hh.

std::shared_ptr<RichardsLMH> HM_Iterative::flow_
private

steady or unsteady water flow simulator based on MH scheme

Definition at line 97 of file hm_iterative.hh.

unsigned int HM_Iterative::max_it_
private

Maximal number of iterations.

Definition at line 111 of file hm_iterative.hh.

std::shared_ptr<Elasticity> HM_Iterative::mechanics_
private

solute transport with chemistry through operator splitting

Definition at line 100 of file hm_iterative.hh.

unsigned int HM_Iterative::min_it_
private

Minimal number of iterations to perform.

Definition at line 108 of file hm_iterative.hh.

double HM_Iterative::r_tol_
private

Relative tolerance for difference between two succeeding iterations.

Definition at line 117 of file hm_iterative.hh.

const int HM_Iterative::registrar
staticprivate
Initial value:
= Input::register_class< HM_Iterative, Mesh &, const Input::Record >("Coupling_Iterative")

Definition at line 94 of file hm_iterative.hh.


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