OGS
LiquidFlowProcess.h
Go to the documentation of this file.
1
13#pragma once
14
15#include <memory>
16
17#include "LiquidFlowData.h"
20#include "ProcessLib/Process.h"
22
23namespace MeshLib
24{
25class Element;
26class Mesh;
27template <typename PROP_VAL_TYPE>
28class PropertyVector;
29}
30
31namespace ProcessLib
32{
33namespace LiquidFlow
34{
55class LiquidFlowProcess final : public Process
56{
57public:
59 std::string name, MeshLib::Mesh& mesh,
60 std::unique_ptr<AbstractJacobianAssembler>&& jacobian_assembler,
61 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
62 parameters,
63 unsigned const integration_order,
64 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
65 process_variables,
66 LiquidFlowData&& process_data,
67 SecondaryVariableCollection&& secondary_variables,
68 std::unique_ptr<ProcessLib::SurfaceFluxData>&& surfaceflux,
69 bool const is_linear);
70
71 void computeSecondaryVariableConcrete(double const t, double const dt,
72 std::vector<GlobalVector*> const& x,
73 GlobalVector const& x_prev,
74 int const process_id) override;
75
76 bool isLinear() const override { return _is_linear; }
77
78 Eigen::Vector3d getFlux(std::size_t const element_id,
79 MathLib::Point3d const& p,
80 double const t,
81 std::vector<GlobalVector*> const& x) const override;
82
83 void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
84 std::vector<GlobalVector*> const& x_prev,
85 const double t,
86 const double dt,
87 int const process_id) override;
88
89private:
91 NumLib::LocalToGlobalIndexMap const& dof_table,
92 MeshLib::Mesh const& mesh, unsigned const integration_order) override;
93
94 void assembleConcreteProcess(const double t, double const dt,
95 std::vector<GlobalVector*> const& x,
96 std::vector<GlobalVector*> const& x_prev,
97 int const process_id, GlobalMatrix& M,
98 GlobalMatrix& K, GlobalVector& b) override;
99
101 const double t, double const dt, std::vector<GlobalVector*> const& x,
102 std::vector<GlobalVector*> const& x_prev, int const process_id,
103 GlobalVector& b, GlobalMatrix& Jac) override;
104
106
107 std::vector<std::unique_ptr<LiquidFlowLocalAssemblerInterface>>
109
110 std::unique_ptr<ProcessLib::SurfaceFluxData> _surfaceflux;
112 bool _is_linear = false;
113};
114
115} // namespace LiquidFlow
116} // namespace ProcessLib
Global vector based on Eigen vector.
Definition EigenVector.h:25
A class to simulate the liquid flow process in porous media described by.
Eigen::Vector3d getFlux(std::size_t const element_id, MathLib::Point3d const &p, double const t, std::vector< GlobalVector * > const &x) const override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalVector &b, GlobalMatrix &Jac) override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, const double t, const double dt, int const process_id) override
std::vector< std::unique_ptr< LiquidFlowLocalAssemblerInterface > > _local_assemblers
void computeSecondaryVariableConcrete(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id) override
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
std::unique_ptr< ProcessLib::SurfaceFluxData > _surfaceflux
LiquidFlowProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< 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, LiquidFlowData &&process_data, SecondaryVariableCollection &&secondary_variables, std::unique_ptr< ProcessLib::SurfaceFluxData > &&surfaceflux, bool const is_linear)
MeshLib::PropertyVector< double > * _hydraulic_flow
void assembleConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override
std::string const name
Definition Process.h:362
Handles configuration of several secondary variables from the project file.