42 using namespace
Input::Type;
47 return Record(
"DualPorosity",
48 "Dual porosity model in transport problems.\n" 49 "Provides computing the concentration of substances in mobile and immobile zone.\n" 53 "Containes region specific data necessary to construct dual porosity model.")
55 "Tolerance according to which the explicit Euler scheme is used or not." 56 "Set 0.0 to use analytic formula only (can be slower).")
61 EqData().output_fields.make_output_type(
"DualPorosity",
""),
62 IT::Default(
"{ \"fields\": [ \"conc_immobile\" ] }"),
63 "Setting of the fields output.")
68 Input::register_class< DualPorosity, Mesh &, Input::Record >(
"DualPorosity") +
74 .
name(
"diffusion_rate_immobile")
75 .
description(
"Diffusion coefficient of non-equilibrium linear exchange between mobile and immobile zone.")
80 .
name(
"porosity_immobile")
87 .
name(
"init_conc_immobile")
88 .
description(
"Initial concentration of substances in the immobile zone.")
99 .
name(
"conc_immobile")
187 .set_time_governor(*
time_);
196 .set_time_governor(*
time_);
251 std::stringstream ss;
307 conc_mob, conc_immob,
308 previous_conc_mob, previous_conc_immob,
325 temp_exponent = (por_mob + por_immob) / (por_mob * por_immob) *
time_->
dt();
329 exponent = diff_vec[sbi] * temp_exponent;
335 conc_average = ((por_mob * previous_conc_mob) + (por_immob * previous_conc_immob))
336 / (por_mob + por_immob);
338 conc_max = std::max(previous_conc_mob-conc_average, previous_conc_immob-conc_average);
343 if(
time_->
dt() <= por_mob*por_immob/(max(diff_vec)*(por_mob+por_immob)) &&
346 double temp = diff_vec[sbi]*(previous_conc_immob - previous_conc_mob) *
time_->
dt();
348 conc_mob = temp / por_mob + previous_conc_mob;
351 conc_immob = -temp / por_immob + previous_conc_immob;
355 double temp = exp(-exponent);
357 conc_mob = (previous_conc_mob - conc_average) * temp + conc_average;
360 conc_immob = (previous_conc_immob - conc_average) * temp + conc_average;
379 for (sbi = 0; sbi < n_subst; sbi++) {
427 bool cfl_changed =
false;
438 time_constraint = std::min(time_constraint, cfl_immobile);
void output_type(OutputTime::DiscreteSpace rt)
unsigned int size() const
get global size
void set_input_list(Input::Array input_list)
Sorption model in immobile zone in case dual porosity is considered.
void set_initial_condition()
Sets initial condition from input.
VecScatter vconc_out_scatter
Output vector scatter.
std::shared_ptr< ReactionTerm > reaction_mobile
Reaction running in mobile zone.
void output_data(void) override
Main output routine.
double ** concentration_matrix_
const std::vector< std::string > & names()
MultiField< 3, FieldValue< 3 >::Scalar > conc_immobile
Calculated concentrations in the immobile zone.
void output(TimeStep step)
SubstanceList substances_
Names belonging to substances.
EquationOutput output_fields
Fields indended for output, i.e. all input fields plus those representing solution.
RegionSet get_region_set(const string &set_name) const
void zero_time_step() override
Field< 3, FieldValue< 3 >::Scalar > porosity_immobile
Immobile porosity field.
~DualPorosity(void)
Destructor.
std::vector< VectorSeqDouble > conc_immobile_out
concentration array output for immobile phase (gathered - sequential)
const RegionDB & region_db() const
IdxInt * el_4_loc_
Indices of elements belonging to local dofs.
ElementAccessor< 3 > element_accessor(unsigned int idx, bool boundary=false)
void allocate_output_mpi(void)
Allocates petsc vectors, prepares them for output and creates output vector scatter.
const TimeStep & step(int index=-1) const
static Input::Type::Abstract & it_abstract_mobile_term()
FieldCommon & units(const UnitSI &units)
Set basic units of the field.
unsigned int size() const
std::shared_ptr< OutputTime > output_stream_
Pointer to a transport output stream.
void update_solution(void) override
ReactionTerm(Mesh &init_mesh, Input::Record in_rec)
Constructor.
Class ReactionTerm is an abstract class representing reaction term in transport.
Class Dual_por_exchange implements the model of dual porosity.
double scheme_tolerance_
Dual porosity computational scheme tolerance.
std::shared_ptr< ReactionTerm > reaction_immobile
Reaction running in immobile zone.
void make_reactions()
Resolves construction of following reactions.
unsigned int n_elements() const
void initialize() override
Prepares the object to usage.
Distribution * distribution_
Pointer to reference to distribution of elements between processors.
static constexpr Mask equation_result
Match result fields. These are never given by input or copy of input.
static constexpr Mask input_copy
FieldCommon & input_default(const string &input_default)
#define START_TIMER(tag)
Starts a timer with specified tag.
#define OLD_ASSERT_LESS(a, b)
bool evaluate_time_constraint(double &time_constraint) override
Computes a constraint for time step.
IdxInt * row_4_el_
Indices of rows belonging to elements.
void initialize_fields()
Initializes field sets.
static Input::Type::Abstract & it_abstract_immobile_term()
virtual Value::return_type const & value(const Point &p, const ElementAccessor< spacedim > &elm) const
EqData()
Collect all fields.
Class representing dual porosity model in transport.
static Input::Type::Abstract & it_abstract_term()
Support classes for parallel programing.
Sorption model in mobile zone in case dual porosity is considered.
void output_vector_gather(void) override
Gathers all the parallel vectors to enable them to be output.
FieldCommon & description(const string &description)
void initialize(std::shared_ptr< OutputTime > stream, Mesh *mesh, Input::Record in_rec, const TimeGovernor &tg)
Input::Record input_record_
void set_components(const std::vector< string > &names)
double ** compute_reaction(double **concentrations, int loc_el) override
bool set_time(const TimeStep &time, LimitSide limit_side)
#define WarningOut()
Macro defining 'warning' record of log.
FieldCommon & name(const string &name)
#define END_TIMER(tag)
Ends a timer with specified tag.
MultiField< 3, FieldValue< 3 >::Scalar > diffusion_rate_immobile
Mass transfer coefficients between mobile and immobile pores.
arma::vec::fixed< spacedim > centre() const
static const int registrar
Registrar of class to factory.
void set_mesh(const Mesh &mesh)
MultiField< 3, FieldValue< 3 >::Scalar > init_conc_immobile
Initial concentrations in the immobile zone.
bool is_field_output_time(const FieldCommon &field, TimeStep step) const
FieldCommon & set_limits(double min, double max=std::numeric_limits< double >::max())
Field< 3, FieldValue< 3 >::Scalar > porosity
Porosity field.
FieldCommon & flags(FieldFlag::Flags::Mask mask)
Class for representation SI units of Fields.
static const Input::Type::Record & get_input_type()
static UnitSI & dimensionless()
Returns dimensionless unit.
static bool print_message_table(ostream &stream, std::string equation_name)
Vec * vconc_immobile
PETSC concentration vector for immobile phase (parallel).
This file contains classes representing sorption model. Sorption model can be computed both in case t...
#define FLOW123D_FORCE_LINK_IN_CHILD(x)
ElementVector element
Vector of elements of the mesh.
unsigned int lsize(int proc) const
get local size