OGS
|
Thermally induced deformation process in linear kinematics poro-mechanical/biphasic model.
The mixture momentum balance, the mixture mass balance and the mixture energy balance are solved under fully saturated conditions.
Definition at line 29 of file ThermoHydroMechanicsProcess.h.
#include <ThermoHydroMechanicsProcess.h>
Public Member Functions | |
ThermoHydroMechanicsProcess (std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const ¶meters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> &&process_variables, ThermoHydroMechanicsProcessData< DisplacementDim > &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme) | |
MathLib::MatrixSpecifications | getMatrixSpecifications (const int process_id) const override |
ODESystem interface | |
bool | isLinear () const override |
Public Member Functions inherited from ProcessLib::Process | |
Process (std::string name_, MeshLib::Mesh &mesh, std::unique_ptr< AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const ¶meters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> &&process_variables, SecondaryVariableCollection &&secondary_variables, const bool use_monolithic_scheme=true) | |
void | preTimestep (std::vector< GlobalVector * > const &x, const double t, const double delta_t, const int process_id) |
Preprocessing before starting assembly for new timestep. More... | |
void | postTimestep (std::vector< GlobalVector * > const &x, const double t, const double delta_t, int const process_id) |
Postprocessing after a complete timestep. More... | |
void | postNonLinearSolver (GlobalVector const &x, GlobalVector const &xdot, const double t, double const dt, int const process_id) |
void | preIteration (const unsigned iter, GlobalVector const &x) final |
void | computeSecondaryVariable (double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id) |
compute secondary variables for the coupled equations or for output. More... | |
NumLib::IterationResult | postIteration (GlobalVector const &x) final |
void | initialize () |
void | setInitialConditions (std::vector< GlobalVector * > &process_solutions, std::vector< GlobalVector * > const &process_solutions_prev, double const t, int const process_id) |
MathLib::MatrixSpecifications | getMatrixSpecifications (const int process_id) const override |
void | setCoupledSolutionsForStaggeredScheme (CoupledSolutionsForStaggeredScheme *const coupled_solutions) |
void | updateDeactivatedSubdomains (double const time, const int process_id) |
bool | isMonolithicSchemeUsed () const |
virtual void | setCoupledTermForTheStaggeredSchemeToLocalAssemblers (int const) |
virtual void | extrapolateIntegrationPointValuesToNodes (const double, std::vector< GlobalVector * > const &, std::vector< GlobalVector * > &) |
void | preAssemble (const double t, double const dt, GlobalVector const &x) final |
void | assemble (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) final |
void | assembleWithJacobian (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) final |
std::vector< NumLib::IndexValueVector< GlobalIndexType > > const * | getKnownSolutions (double const t, GlobalVector const &x, int const process_id) const final |
MeshLib::Mesh & | getMesh () const |
std::vector< std::reference_wrapper< ProcessVariable > > const & | getProcessVariables (const int process_id) const |
SecondaryVariableCollection const & | getSecondaryVariables () const |
std::vector< std::unique_ptr< IntegrationPointWriter > > const * | getIntegrationPointWriter (MeshLib::Mesh const &mesh) const |
virtual Eigen::Vector3d | getFlux (std::size_t, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &) const |
virtual void | solveReactionEquation (std::vector< GlobalVector * > &, std::vector< GlobalVector * > const &, double const, double const, NumLib::EquationSystem &, int const) |
Private Member Functions | |
void | constructDofTable () override |
void | initializeConcreteProcess (NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override |
Process specific initialization called by initialize(). More... | |
void | initializeBoundaryConditions () override |
void | assembleConcreteProcess (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override |
void | assembleWithJacobianConcreteProcess (const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) override |
void | preTimestepConcreteProcess (std::vector< GlobalVector * > const &x, double const t, double const dt, const int process_id) override |
void | postTimestepConcreteProcess (std::vector< GlobalVector * > const &x, const double t, const double dt, int const process_id) override |
void | postNonLinearSolverConcreteProcess (GlobalVector const &x, GlobalVector const &xdot, const double t, double const dt, int const process_id) override |
NumLib::LocalToGlobalIndexMap const & | getDOFTable (const int process_id) const override |
void | computeSecondaryVariableConcrete (double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, const int process_id) override |
std::tuple< NumLib::LocalToGlobalIndexMap *, bool > | getDOFTableForExtrapolatorData () const override |
bool | hasMechanicalProcess (int const process_id) const |
Private Attributes | |
std::vector< MeshLib::Node * > | _base_nodes |
std::unique_ptr< MeshLib::MeshSubset const > | _mesh_subset_base_nodes |
ThermoHydroMechanicsProcessData< DisplacementDim > | _process_data |
std::vector< std::unique_ptr< LocalAssemblerInterface > > | _local_assemblers |
std::unique_ptr< NumLib::LocalToGlobalIndexMap > | _local_to_global_index_map_single_component |
std::unique_ptr< NumLib::LocalToGlobalIndexMap > | _local_to_global_index_map_with_base_nodes |
GlobalSparsityPattern | _sparsity_pattern_with_linear_element |
MeshLib::PropertyVector< double > * | _nodal_forces = nullptr |
MeshLib::PropertyVector< double > * | _hydraulic_flow = nullptr |
MeshLib::PropertyVector< double > * | _heat_flux = nullptr |
ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::ThermoHydroMechanicsProcess | ( | std::string | name, |
MeshLib::Mesh & | mesh, | ||
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > && | jacobian_assembler, | ||
std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const & | parameters, | ||
unsigned const | integration_order, | ||
std::vector< std::vector< std::reference_wrapper< ProcessVariable >>> && | process_variables, | ||
ThermoHydroMechanicsProcessData< DisplacementDim > && | process_data, | ||
SecondaryVariableCollection && | secondary_variables, | ||
bool const | use_monolithic_scheme | ||
) |
Definition at line 27 of file ThermoHydroMechanicsProcess.cpp.
References ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::_heat_flux, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::_hydraulic_flow, ProcessLib::Process::_integration_point_writer, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::_local_assemblers, ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::_nodal_forces, MeshLib::Mesh::getDimension(), ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface::getEpsilon(), ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface::getSigma(), and MeshLib::Node.
|
overrideprivatevirtual |
Implements ProcessLib::Process.
Definition at line 308 of file ThermoHydroMechanicsProcess.cpp.
References ProcessLib::VectorMatrixAssembler::assemble(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), and ProcessLib::ProcessVariable::getActiveElementIDs().
|
overrideprivatevirtual |
Implements ProcessLib::Process.
Definition at line 327 of file ThermoHydroMechanicsProcess.cpp.
References ProcessLib::VectorMatrixAssembler::assembleWithJacobian(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and NumLib::transformVariableFromGlobalVector().
|
overrideprivatevirtual |
Reimplemented from ProcessLib::Process.
Definition at line 474 of file ThermoHydroMechanicsProcess.cpp.
References ProcessLib::LocalAssemblerInterface::computeSecondaryVariable(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), and ProcessLib::ProcessVariable::getActiveElementIDs().
|
overrideprivatevirtual |
This function is for general cases, in which all equations of the coupled processes have the same number of unknowns. For the general cases with the staggered scheme, all equations of the coupled processes share one DOF table hold by _local_to_global_index_map
. Other cases can be considered by overloading this member function in the derived class.
Reimplemented from ProcessLib::Process.
Definition at line 94 of file ThermoHydroMechanicsProcess.cpp.
References NumLib::BY_LOCATION, NumLib::computeSparsityPattern(), and MeshLib::getBaseNodes().
|
overrideprivatevirtual |
Reimplemented from ProcessLib::Process.
Definition at line 511 of file ThermoHydroMechanicsProcess.cpp.
|
overrideprivatevirtual |
Get the address of a LocalToGlobalIndexMap, and the status of its memory. If the LocalToGlobalIndexMap is created as new in this function, the function also returns a true boolean value to let Extrapolator manage the memory by the address returned by this function.
Reimplemented from ProcessLib::Process.
Definition at line 502 of file ThermoHydroMechanicsProcess.cpp.
|
override |
Get the size and the sparse pattern of the global matrix in order to create the global matrices and vectors for the system equations of this process.
process_id | Process ID. If the monolithic scheme is applied, process_id = 0. For the staggered scheme, process_id = 0 represents the hydraulic (H) process, while process_id = 1 represents the mechanical (M) process. |
Definition at line 75 of file ThermoHydroMechanicsProcess.cpp.
|
inlineprivate |
Check whether the process represented by process_id
is/has mechanical process. In the present implementation, the mechanical process has process_id == 1 in the staggered scheme.
Definition at line 137 of file ThermoHydroMechanicsProcess.h.
References ProcessLib::Process::_use_monolithic_scheme.
|
overrideprivatevirtual |
Member function to initialize the boundary conditions for all coupled processes. It is called by initialize().
Reimplemented from ProcessLib::Process.
Definition at line 280 of file ThermoHydroMechanicsProcess.cpp.
|
overrideprivatevirtual |
Process specific initialization called by initialize().
Implements ProcessLib::Process.
Definition at line 174 of file ThermoHydroMechanicsProcess.cpp.
References ProcessLib::ThermoHydroMechanics::createLocalAssemblers(), NumLib::SerialExecutor::executeMemberOnDereferenced(), MeshLib::Properties::existsPropertyVector(), MeshLib::Mesh::getDimension(), MeshLib::Mesh::getElements(), ProcessLib::getIntegrationPointMetaData(), ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface::getIntPtDarcyVelocity(), ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface::getIntPtEpsilon(), ProcessLib::ThermoHydroMechanics::LocalAssemblerInterface::getIntPtSigma(), MeshLib::Mesh::getProperties(), ProcessLib::LocalAssemblerInterface::initialize(), MeshLib::IntegrationPoint, MeshLib::Mesh::isAxiallySymmetric(), ProcessLib::makeExtrapolator(), MaterialPropertyLib::name, MeshLib::Node, and OGS_FATAL.
|
override |
Definition at line 68 of file ThermoHydroMechanicsProcess.cpp.
|
overrideprivatevirtual |
Reimplemented from ProcessLib::Process.
Definition at line 452 of file ThermoHydroMechanicsProcess.cpp.
References DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::LocalAssemblerInterface::postNonLinearSolver().
|
overrideprivatevirtual |
Reimplemented from ProcessLib::Process.
Definition at line 426 of file ThermoHydroMechanicsProcess.cpp.
References DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::LocalAssemblerInterface::postTimestep().
|
overrideprivatevirtual |
Reimplemented from ProcessLib::Process.
Definition at line 407 of file ThermoHydroMechanicsProcess.cpp.
References DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::LocalAssemblerInterface::preTimestep().
|
private |
Definition at line 106 of file ThermoHydroMechanicsProcess.h.
|
private |
Definition at line 144 of file ThermoHydroMechanicsProcess.h.
Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::ThermoHydroMechanicsProcess().
|
private |
Definition at line 143 of file ThermoHydroMechanicsProcess.h.
Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::ThermoHydroMechanicsProcess().
|
private |
Definition at line 110 of file ThermoHydroMechanicsProcess.h.
Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::ThermoHydroMechanicsProcess().
|
private |
Definition at line 113 of file ThermoHydroMechanicsProcess.h.
|
private |
Local to global index mapping for base nodes, which is used for linear interpolation for pressure in the staggered scheme.
Definition at line 118 of file ThermoHydroMechanicsProcess.h.
|
private |
Definition at line 107 of file ThermoHydroMechanicsProcess.h.
|
private |
Definition at line 142 of file ThermoHydroMechanicsProcess.h.
Referenced by ProcessLib::ThermoHydroMechanics::ThermoHydroMechanicsProcess< DisplacementDim >::ThermoHydroMechanicsProcess().
|
private |
Definition at line 108 of file ThermoHydroMechanicsProcess.h.
|
private |
Sparsity pattern for the flow equation, and it is initialized only if the staggered scheme is used.
Definition at line 122 of file ThermoHydroMechanicsProcess.h.