OGS
ComponentTransportProcess.h
Go to the documentation of this file.
1 
11 #pragma once
12 
13 #include "ComponentTransportFEM.h"
16 #include "ProcessLib/Process.h"
17 
18 namespace ChemistryLib
19 {
20 class ChemicalSolverInterface;
21 }
22 
23 namespace ProcessLib
24 {
25 struct SurfaceFluxData;
26 
27 namespace ComponentTransport
28 {
94 class ComponentTransportProcess final : public Process
95 {
96 public:
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 
125  int const process_id) override;
126 
127  void solveReactionEquation(std::vector<GlobalVector*>& x,
128  std::vector<GlobalVector*> const& x_prev,
129  double const t, double const dt,
130  NumLib::EquationSystem& ode_sys,
131  int const process_id) override;
132 
134  const double t,
135  std::vector<GlobalVector*> const& integration_point_values_vectors,
136  std::vector<GlobalVector*>& nodal_values_vectors) override;
137 
138  void computeSecondaryVariableConcrete(double const /*t*/,
139  double const /*dt*/,
140  std::vector<GlobalVector*> const& x,
141  GlobalVector const& /*x_dot*/,
142  int const /*process_id*/) override;
143 
144  void postTimestepConcreteProcess(std::vector<GlobalVector*> const& x,
145  const double t,
146  const double dt,
147  int const process_id) override;
148 
149 private:
151  NumLib::LocalToGlobalIndexMap const& dof_table,
152  MeshLib::Mesh const& mesh,
153  unsigned const integration_order) override;
154 
155  void setInitialConditionsConcreteProcess(std::vector<GlobalVector*>& x,
156  double const t,
157  int const process_id) override;
158 
159  void assembleConcreteProcess(const double t, double const dt,
160  std::vector<GlobalVector*> const& x,
161  std::vector<GlobalVector*> const& xdot,
162  int const process_id, GlobalMatrix& M,
163  GlobalMatrix& K, GlobalVector& b) override;
164 
166  const double t, double const dt, std::vector<GlobalVector*> const& x,
167  std::vector<GlobalVector*> const& xdot, const double dxdot_dx,
168  const double dx_dx, int const process_id, GlobalMatrix& M,
169  GlobalMatrix& K, GlobalVector& b, GlobalMatrix& Jac) override;
170 
172 
173  std::vector<std::unique_ptr<ComponentTransportLocalAssemblerInterface>>
175 
176  std::unique_ptr<ProcessLib::SurfaceFluxData> _surfaceflux;
177 
178  std::unique_ptr<ChemistryLib::ChemicalSolverInterface>
180 };
181 
182 } // namespace ComponentTransport
183 } // namespace ProcessLib
Global vector based on Eigen vector.
Definition: EigenVector.h:26
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 extrapolateIntegrationPointValuesToNodes(const double t, std::vector< GlobalVector * > const &integration_point_values_vectors, std::vector< GlobalVector * > &nodal_values_vectors) 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 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 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
void setCoupledTermForTheStaggeredSchemeToLocalAssemblers(int const process_id) override
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
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::string const name
Definition: Process.h:323
Handles configuration of several secondary variables from the project file.