Detailed Description

A class to simulate the liquid flow process in porous media described by.

\[ \frac{\partial n \rho_l}{\partial T} \frac{\partial T}{\partial t}/\rho_l + (\frac{\partial n \rho_l}{\partial p}/\rho_l + \beta_s) \frac{\partial p}{\partial t} -\nabla (\frac{K}{\mu}(\nabla p + \rho_l g \nabla z) ) = Q \]

where

\begin{eqnarray*} &p:& \mbox{pore pressure,}\\ &T: & \mbox{Temperature,}\\ &\rho_l:& \mbox{liquid density,}\\ &\beta_s:& \mbox{specific storage,}\\ &K:& \mbox{permeability,}\\ &\mu:& \mbox{viscosity,}\\ \end{eqnarray*}

Definition at line 55 of file LiquidFlowProcess.h.

#include <LiquidFlowProcess.h>

Inheritance diagram for ProcessLib::LiquidFlow::LiquidFlowProcess:

[legend]

Collaboration diagram for ProcessLib::LiquidFlow::LiquidFlowProcess:

[legend]

Public Member Functions
	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)

void	computeSecondaryVariableConcrete (double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id) override

bool	isLinear () const override

Eigen::Vector3d	getFlux (std::size_t const element_id, MathLib::Point3d const &p, double const t, std::vector< GlobalVector * > const &x) const override

void	postTimestepConcreteProcess (std::vector< GlobalVector * > const &x, const double t, const double dt, int const process_id) 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 &parameters, 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

virtual NumLib::LocalToGlobalIndexMap const &	getDOFTable (const int) const

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 void	solveReactionEquation (std::vector< GlobalVector * > &, std::vector< GlobalVector * > const &, double const, double const, NumLib::EquationSystem &, int const)

Private Member Functions
void	initializeConcreteProcess (NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
	Process specific initialization called by initialize(). More...

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

Private Attributes
LiquidFlowData	_process_data

std::vector< std::unique_ptr< LiquidFlowLocalAssemblerInterface > >	_local_assemblers

std::unique_ptr< ProcessLib::SurfaceFluxData >	_surfaceflux

MeshLib::PropertyVector< double > *	_hydraulic_flow = nullptr

Additional Inherited Members
Public Types inherited from ProcessLib::Process
using	NonlinearSolver = NumLib::NonlinearSolverBase

using	TimeDiscretization = NumLib::TimeDiscretization

Public Attributes inherited from ProcessLib::Process
std::string const	name

Protected Member Functions inherited from ProcessLib::Process
NumLib::Extrapolator &	getExtrapolator () const

NumLib::LocalToGlobalIndexMap const &	getSingleComponentDOFTable () const

void	initializeProcessBoundaryConditionsAndSourceTerms (const NumLib::LocalToGlobalIndexMap &dof_table, const int process_id)

virtual void	constructDofTable ()

void	constructMonolithicProcessDofTable ()

void	constructDofTableOfSpecifiedProcessStaggeredScheme (const int specified_prosess_id)

virtual std::tuple< NumLib::LocalToGlobalIndexMap *, bool >	getDOFTableForExtrapolatorData () const

Protected Attributes inherited from ProcessLib::Process
MeshLib::Mesh &	_mesh

std::unique_ptr< MeshLib::MeshSubset const >	_mesh_subset_all_nodes

std::unique_ptr< NumLib::LocalToGlobalIndexMap >	_local_to_global_index_map

SecondaryVariableCollection	_secondary_variables

VectorMatrixAssembler	_global_assembler

const bool	_use_monolithic_scheme

CoupledSolutionsForStaggeredScheme *	_coupled_solutions

unsigned const	_integration_order

std::vector< std::unique_ptr< IntegrationPointWriter > >	_integration_point_writer

GlobalSparsityPattern	_sparsity_pattern

std::vector< std::vector< std::reference_wrapper< ProcessVariable > > >	_process_variables

std::vector< BoundaryConditionCollection >	_boundary_conditions

Constructor & Destructor Documentation

◆ LiquidFlowProcess()

ProcessLib::LiquidFlow::LiquidFlowProcess::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
	)

Definition at line 28 of file LiquidFlowProcess.cpp.

     : Process(std::move(name), mesh, std::move(jacobian_assembler), parameters,
               integration_order, std::move(process_variables),
               std::move(secondary_variables)),
       _process_data(std::move(process_data)),
       _surfaceflux(std::move(surfaceflux))
 {
     DBUG("Create Liquid flow process.");
  
     _hydraulic_flow = MeshLib::getOrCreateMeshProperty<double>(
         mesh, "HydraulicFlow", MeshLib::MeshItemType::Node, 1);
 }

References _hydraulic_flow, DBUG(), and MeshLib::Node.

Member Function Documentation

◆ assembleConcreteProcess()

void ProcessLib::LiquidFlow::LiquidFlowProcess::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
	)

overrideprivatevirtual

Implements ProcessLib::Process.

Definition at line 72 of file LiquidFlowProcess.cpp.

 {
     DBUG("Assemble LiquidFlowProcess.");
  
     std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
         dof_table = {std::ref(*_local_to_global_index_map)};
  
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     // Call global assembler for each local assembly item.
     GlobalExecutor::executeSelectedMemberDereferenced(
         _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
         pv.getActiveElementIDs(), dof_table, t, dt, x, xdot, process_id, M, K,
         b);
  
     MathLib::finalizeMatrixAssembly(M);
     MathLib::finalizeMatrixAssembly(K);
     MathLib::finalizeVectorAssembly(b);
  
     auto const residuum = computeResiduum(*x[0], *xdot[0], M, K, b);
     transformVariableFromGlobalVector(residuum, 0 /*variable id*/,
                                       *_local_to_global_index_map,
                                       *_hydraulic_flow, std::negate<double>());
 }

References ProcessLib::Process::_global_assembler, _hydraulic_flow, _local_assemblers, ProcessLib::Process::_local_to_global_index_map, ProcessLib::VectorMatrixAssembler::assemble(), ProcessLib::computeResiduum(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), MathLib::finalizeMatrixAssembly(), MathLib::finalizeVectorAssembly(), ProcessLib::ProcessVariable::getActiveElementIDs(), ProcessLib::Process::getProcessVariables(), and NumLib::transformVariableFromGlobalVector().

◆ assembleWithJacobianConcreteProcess()

void ProcessLib::LiquidFlow::LiquidFlowProcess::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
	)

overrideprivatevirtual

Implements ProcessLib::Process.

Definition at line 100 of file LiquidFlowProcess.cpp.

 {
     DBUG("AssembleWithJacobian LiquidFlowProcess.");
  
     std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
         dof_table = {std::ref(*_local_to_global_index_map)};
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     // Call global assembler for each local assembly item.
     GlobalExecutor::executeSelectedMemberDereferenced(
         _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
         _local_assemblers, pv.getActiveElementIDs(), dof_table, t, dt, x, xdot,
         dxdot_dx, dx_dx, process_id, M, K, b, Jac);
 }

References ProcessLib::Process::_global_assembler, _local_assemblers, ProcessLib::Process::_local_to_global_index_map, ProcessLib::VectorMatrixAssembler::assembleWithJacobian(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::Process::getProcessVariables().

◆ computeSecondaryVariableConcrete()

void ProcessLib::LiquidFlow::LiquidFlowProcess::computeSecondaryVariableConcrete	(	double const	t,
		double const	dt,
		std::vector< GlobalVector * > const &	x,
		GlobalVector const &	x_dot,
		int const	process_id
	)

overridevirtual

Reimplemented from ProcessLib::Process.

Definition at line 119 of file LiquidFlowProcess.cpp.

 {
     DBUG("Compute the velocity for LiquidFlowProcess.");
     std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
     dof_tables.reserve(x.size());
     std::generate_n(std::back_inserter(dof_tables), x.size(),
                     [&]() { return _local_to_global_index_map.get(); });
  
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LiquidFlowLocalAssemblerInterface::computeSecondaryVariable,
         _local_assemblers, pv.getActiveElementIDs(), dof_tables, t, dt, x,
         x_dot, process_id);
 }

References _local_assemblers, ProcessLib::LocalAssemblerInterface::computeSecondaryVariable(), DBUG(), NumLib::SerialExecutor::executeSelectedMemberOnDereferenced(), ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::Process::getProcessVariables().

◆ getFlux()

Eigen::Vector3d ProcessLib::LiquidFlow::LiquidFlowProcess::getFlux	(	std::size_t const	element_id,
		MathLib::Point3d const &	p,
		double const	t,
		std::vector< GlobalVector * > const &	x
	)		const

overridevirtual

Reimplemented from ProcessLib::Process.

Definition at line 136 of file LiquidFlowProcess.cpp.

 {
     // fetch local_x from primary variable
     std::vector<GlobalIndexType> indices_cache;
     auto const r_c_indices = NumLib::getRowColumnIndices(
         element_id, *_local_to_global_index_map, indices_cache);
     constexpr int process_id = 0;  // monolithic scheme.
     std::vector<double> local_x(x[process_id]->get(r_c_indices.rows));
  
     return _local_assemblers[element_id]->getFlux(p, t, local_x);
 }

References _local_assemblers, ProcessLib::Process::_local_to_global_index_map, and NumLib::getRowColumnIndices().

◆ initializeConcreteProcess()

void ProcessLib::LiquidFlow::LiquidFlowProcess::initializeConcreteProcess	(	NumLib::LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh,
		unsigned const	integration_order
	)

overrideprivatevirtual

Process specific initialization called by initialize().

Implements ProcessLib::Process.

Definition at line 51 of file LiquidFlowProcess.cpp.

 {
     const int process_id = 0;
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
  
     int const mesh_space_dimension = _process_data.mesh_space_dimension;
     ProcessLib::createLocalAssemblers<LiquidFlowLocalAssembler>(
         mesh_space_dimension, mesh.getElements(), dof_table,
         pv.getShapeFunctionOrder(), _local_assemblers,
         mesh.isAxiallySymmetric(), integration_order, _process_data);
  
     _secondary_variables.addSecondaryVariable(
         "darcy_velocity",
         makeExtrapolator(
             mesh_space_dimension, getExtrapolator(), _local_assemblers,
             &LiquidFlowLocalAssemblerInterface::getIntPtDarcyVelocity));
 }

References _local_assemblers, _process_data, ProcessLib::Process::_secondary_variables, ProcessLib::SecondaryVariableCollection::addSecondaryVariable(), MeshLib::Mesh::getElements(), ProcessLib::Process::getExtrapolator(), ProcessLib::LiquidFlow::LiquidFlowLocalAssemblerInterface::getIntPtDarcyVelocity(), ProcessLib::Process::getProcessVariables(), ProcessLib::ProcessVariable::getShapeFunctionOrder(), MeshLib::Mesh::isAxiallySymmetric(), ProcessLib::makeExtrapolator(), and ProcessLib::LiquidFlow::LiquidFlowData::mesh_space_dimension.

◆ isLinear()

bool ProcessLib::LiquidFlow::LiquidFlowProcess::isLinear ( ) const

inlineoverride

Definition at line 76 of file LiquidFlowProcess.h.

76 { return true; }

◆ postTimestepConcreteProcess()

void ProcessLib::LiquidFlow::LiquidFlowProcess::postTimestepConcreteProcess	(	std::vector< GlobalVector * > const &	x,
		const double	t,
		const double	dt,
		int const	process_id
	)

overridevirtual

Reimplemented from ProcessLib::Process.

Definition at line 153 of file LiquidFlowProcess.cpp.

 {
     if (!_surfaceflux)  // computing the surfaceflux is optional
     {
         return;
     }
  
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
     _surfaceflux->integrate(x, t, *this, process_id, _integration_order, _mesh,
                             pv.getActiveElementIDs());
 }

References ProcessLib::Process::_integration_order, ProcessLib::Process::_mesh, _surfaceflux, ProcessLib::ProcessVariable::getActiveElementIDs(), and ProcessLib::Process::getProcessVariables().

Member Data Documentation

◆ _hydraulic_flow

MeshLib::PropertyVector<double>* ProcessLib::LiquidFlow::LiquidFlowProcess::_hydraulic_flow = nullptr

private

Definition at line 111 of file LiquidFlowProcess.h.

Referenced by LiquidFlowProcess(), and assembleConcreteProcess().

◆ _local_assemblers

std::vector<std::unique_ptr<LiquidFlowLocalAssemblerInterface> > ProcessLib::LiquidFlow::LiquidFlowProcess::_local_assemblers

private

Definition at line 108 of file LiquidFlowProcess.h.

Referenced by assembleConcreteProcess(), assembleWithJacobianConcreteProcess(), computeSecondaryVariableConcrete(), getFlux(), and initializeConcreteProcess().

◆ _process_data

LiquidFlowData ProcessLib::LiquidFlow::LiquidFlowProcess::_process_data

private

Definition at line 105 of file LiquidFlowProcess.h.

Referenced by initializeConcreteProcess().

◆ _surfaceflux

std::unique_ptr<ProcessLib::SurfaceFluxData> ProcessLib::LiquidFlow::LiquidFlowProcess::_surfaceflux

private

Definition at line 110 of file LiquidFlowProcess.h.

Referenced by postTimestepConcreteProcess().

The documentation for this class was generated from the following files:

ProcessLib/LiquidFlow/LiquidFlowProcess.h
ProcessLib/LiquidFlow/LiquidFlowProcess.cpp

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ LiquidFlowProcess()

Member Function Documentation

◆ assembleConcreteProcess()

◆ assembleWithJacobianConcreteProcess()

◆ computeSecondaryVariableConcrete()

◆ getFlux()

◆ initializeConcreteProcess()

◆ isLinear()

◆ postTimestepConcreteProcess()

Member Data Documentation

◆ _hydraulic_flow

◆ _local_assemblers

◆ _process_data

◆ _surfaceflux