OGS
ComponentTransportProcess.h
Go to the documentation of this file.
1
11#pragma once
12
16#include "ProcessLib/Process.h"
17
18namespace ChemistryLib
19{
20class ChemicalSolverInterface;
21}
22
23namespace ProcessLib
24{
25struct SurfaceFluxData;
26
27namespace ComponentTransport
28{
95{
96public:
98 std::string name,
99 MeshLib::Mesh& mesh,
100 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&&
101 jacobian_assembler,
102 std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&
103 parameters,
104 unsigned const integration_order,
105 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
106 process_variables,
107 ComponentTransportProcessData&& process_data,
108 SecondaryVariableCollection&& secondary_variables,
109 bool const use_monolithic_scheme,
110 std::unique_ptr<ProcessLib::SurfaceFluxData>&& surfaceflux,
111 std::unique_ptr<ChemistryLib::ChemicalSolverInterface>&&
112 chemical_solver_interface);
113
116
117 bool isLinear() const override { return false; }
119
120 Eigen::Vector3d getFlux(std::size_t const element_id,
121 MathLib::Point3d const& p, double const t,
122 std::vector<GlobalVector*> const& x) const override;
123
124 void solveReactionEquation(std::vector<GlobalVector*>& x,
125 std::vector<GlobalVector*> const& x_prev,
126 double const t, double const dt,
127 NumLib::EquationSystem& ode_sys,
128 int const process_id) override;
129
130 void computeSecondaryVariableConcrete(double const /*t*/,
131 double const /*dt*/,
132 std::vector<GlobalVector*> const& x,
133 GlobalVector const& /*x_dot*/,
134 int const /*process_id*/) override;
135
136 void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
137 const double t,
138 const double dt,
139 int const process_id) override;
140
141private:
143 NumLib::LocalToGlobalIndexMap const& dof_table,
144 MeshLib::Mesh const& mesh,
145 unsigned const integration_order) override;
146
147 void setInitialConditionsConcreteProcess(std::vector<GlobalVector*>& x,
148 double const t,
149 int const process_id) override;
150
151 void assembleConcreteProcess(const double t, double const dt,
152 std::vector<GlobalVector*> const& x,
153 std::vector<GlobalVector*> const& xdot,
154 int const process_id, GlobalMatrix& M,
155 GlobalMatrix& K, GlobalVector& b) override;
156
158 const double t, double const dt, std::vector<GlobalVector*> const& x,
159 std::vector<GlobalVector*> const& xdot, int const process_id,
161 GlobalMatrix& Jac) override;
162
164
165 std::vector<std::unique_ptr<ComponentTransportLocalAssemblerInterface>>
167
168 std::unique_ptr<ProcessLib::SurfaceFluxData> _surfaceflux;
169
170 std::unique_ptr<ChemistryLib::ChemicalSolverInterface>
172
174};
175
176} // namespace ComponentTransport
177} // namespace ProcessLib
Global vector based on Eigen vector.
Definition: EigenVector.h:28
void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &x, double const t, 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< ChemistryLib::ChemicalSolverInterface > _chemical_solver_interface
void assembleWithJacobianConcreteProcess(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, GlobalMatrix &Jac) 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 computeSecondaryVariableConcrete(double const, double const, std::vector< GlobalVector * > const &x, GlobalVector const &, int const) override
void solveReactionEquation(std::vector< GlobalVector * > &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, NumLib::EquationSystem &ode_sys, int const process_id) override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, const double t, const double dt, int const process_id) override
ComponentTransportProcess(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, ComponentTransportProcessData &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme, std::unique_ptr< ProcessLib::SurfaceFluxData > &&surfaceflux, std::unique_ptr< ChemistryLib::ChemicalSolverInterface > &&chemical_solver_interface)
std::vector< std::unique_ptr< ComponentTransportLocalAssemblerInterface > > _local_assemblers
std::unique_ptr< ProcessLib::SurfaceFluxData > _surfaceflux
Eigen::Vector3d getFlux(std::size_t const element_id, MathLib::Point3d const &p, double const t, std::vector< GlobalVector * > const &x) const override
std::string const name
Definition: Process.h:325
Handles configuration of several secondary variables from the project file.
static const double t