33 #ifndef DARCY_FLOW_MH_HH
34 #define DARCY_FLOW_MH_HH
130 <<
"Diverged nonlinear solver. Reason: " << EI_Reason::val
133 <<
"Missing the key 'time', obligatory for the transient problems.");
218 std::shared_ptr<DOFHandlerMultiDim>
dh_;
219 std::shared_ptr<SubDOFHandlerMultiDim>
dh_cr_;
272 {
eq_fields_->extra_storativity = extra_stor; }
403 {
return "Flow_Darcy_MH";}
410 void dofs_range(
unsigned int n_dofs,
unsigned int &min,
unsigned int &max,
unsigned int component);
MortarMethod
Type of experimental Mortar-like method for non-compatible 1d-2d interaction.
std::shared_ptr< SubDOFHandlerMultiDim > dh_cr_
DOF handler represents DOFs of edges.
MortarMethod mortar_method_
int is_linear
Hack fo BDDC solver.
MultidimAssembly multidim_assembler
unsigned int n_schur_compls
std::shared_ptr< DOFHandlerMultiDim > dh_
full DOF handler represents DOFs of sides, elements and edges
std::vector< char > bc_switch_dirichlet
Idicator of dirichlet or neumann type of switch boundary conditions.
std::shared_ptr< DOFHandlerMultiDim > dh_cr_disc_
DOF handler represents DOFs of sides.
std::shared_ptr< Balance > balance
bool force_no_neumann_bc
auxiliary flag for switchting Dirichlet like BC
Field< 3, FieldValue< 3 >::Scalar > init_piezo_head
Same as previous but used in boundary fields.
BCField< 3, FieldValue< 3 >::Scalar > bc_piezo_head
BCField< 3, FieldValue< 3 >::Enum > bc_type
BCField< 3, FieldValue< 3 >::Scalar > bc_switch_piezo_head
BCField< 3, FieldValue< 3 >::VectorFixed > bc_gravity
Holds gravity vector acceptable in FieldModel.
Field< 3, FieldValue< 3 >::Scalar > extra_source
Externally added storativity.
Field< 3, FieldValue< 3 >::Scalar > field_ele_piezo_head
Field< 3, FieldValue< 3 >::TensorFixed > anisotropy
Field< 3, FieldValue< 3 >::Scalar > conductivity
Field< 3, FieldValue< 3 >::Scalar > field_ele_pressure
Externally added water source.
EqFields()
Creation of all fields.
Field< 3, FieldValue< 3 >::VectorFixed > gravity_field
Field< 3, FieldValue< 3 >::Scalar > field_edge_pressure
Field< 3, FieldValue< 3 >::Scalar > sigma
Field< 3, FieldValue< 3 >::Scalar > init_pressure
Field< 3, FieldValue< 3 >::Scalar > cross_section
Field< 3, FieldValue< 3 >::Scalar > water_source_density
FieldCoords & X()
Return coords field.
BCField< 3, FieldValue< 3 >::Scalar > bc_pressure
BCField< 3, FieldValue< 3 >::Scalar > bc_switch_pressure
BCField< 3, FieldValue< 3 >::Scalar > bc_robin_sigma
BCField< 3, FieldValue< 3 >::Scalar > bc_flux
static const Input::Type::Selection & get_bc_type_selection()
Return a Selection corresponding to enum BC_Type.
Field< 3, FieldValue< 3 >::VectorFixed > field_ele_velocity
Field< 3, FieldValue< 3 >::VectorFixed > flux
Field< 3, FieldValue< 3 >::Scalar > extra_storativity
Field< 3, FieldValue< 3 >::Scalar > storativity
Mixed-hybrid model of linear Darcy flow, possibly unsteady.
DarcyFlowMHOutput * output_object
void set_mesh_data_for_bddc(LinSys_BDDC *bddc_ls)
virtual double solved_time() override
void prepare_new_time_step()
DECLARE_INPUT_EXCEPTION(ExcMissingTimeGovernor,<< "Missing the key 'time', obligatory for the transient problems.")
void assembly_mh_matrix(MultidimAssembly &assembler)
static const Input::Type::Record & type_field_descriptor()
virtual void read_initial_condition()
void print_matlab_matrix(string matlab_file)
Print darcy flow matrix in matlab format into a file.
virtual void initialize_specific()
std::shared_ptr< Balance > balance_
std::vector< int > get_component_indices_vec(unsigned int component) const
Get vector of all DOF indices of given component (0..side, 1..element, 2..edge)
void allocate_mh_matrix()
void initialize() override
static std::string equation_name()
void solve_nonlinear()
Solve method common to zero_time_step and update solution.
virtual ~DarcyMH() override
void set_extra_source(const Field< 3, FieldValue< 3 >::Scalar > &extra_src)
std::shared_ptr< EqData > eq_data_
bool use_steady_assembly_
void zero_time_step() override
std::shared_ptr< EqFields > eq_fields_
static const int registrar
Registrar of class to factory.
virtual double solution_precision() const
unsigned int nonlinear_iteration_
void set_extra_storativity(const Field< 3, FieldValue< 3 >::Scalar > &extra_stor)
virtual void output_data() override
Write computed fields.
virtual void assembly_linear_system()
static const Input::Type::Selection & get_mh_mortar_selection()
Selection for enum MortarMethod.
void reconstruct_solution_from_schur(MultidimAssembly &assembler)
DECLARE_EXCEPTION(ExcSolverDiverge,<< "Diverged nonlinear solver. Reason: "<< EI_Reason::val)
virtual bool zero_time_term(bool time_global=false)
DarcyMH(Mesh &mesh, const Input::Record in_rec, TimeGovernor *tm=nullptr)
CREATE AND FILL GLOBAL MH MATRIX OF THE WATER MODEL.
virtual void postprocess()
postprocess velocity field (add sources)
TYPEDEF_ERR_INFO(EI_Reason, string)
static const Input::Type::Record & get_input_type()
void update_solution() override
void create_linear_system(Input::AbstractRecord rec)
void solve_time_step(bool output=true)
Solve the problem without moving to next time and without output.
virtual void setup_time_term()
virtual void assembly_source_term()
Source term is implemented differently in LMH version.
Container for various descendants of FieldCommonBase.
FieldCoords X_
Field holds coordinates for computing of FieldFormulas.
Class template representing a field with values dependent on: point, element, and region.
Basic time management functionality for unsteady (and steady) solvers (class Equation).
std::vector< std::shared_ptr< AssemblyFlowBase > > MultidimAssembly
void dofs_range(unsigned int n_dofs, unsigned int &min, unsigned int &max, unsigned int component)
Helper method fills range (min and max) of given component.
void mat_count_off_proc_values(Mat m, Vec v)
Basic time management class.