OGS
List of all members | Classes | Public Member Functions | Private Types | Private Member Functions | Private Attributes
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim > Class Template Reference

Detailed Description

template<typename ShapeFunction, int GlobalDim>
class ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >

Definition at line 67 of file LiquidFlowLocalAssembler.h.

#include <LiquidFlowLocalAssembler.h>

Inheritance diagram for ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >:
[legend]
Collaboration diagram for ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >:
[legend]

Classes

struct  AnisotropicCalculator
 
struct  IsotropicCalculator
 

Public Member Functions

 LiquidFlowLocalAssembler (MeshLib::Element const &element, std::size_t const, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, LiquidFlowData const &process_data)
 
void assemble (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data) override
 
Eigen::Vector3d getFlux (MathLib::Point3d const &p_local_coords, double const t, std::vector< double > const &local_x) const override
 
Eigen::Map< const Eigen::RowVectorXd > getShapeMatrix (const unsigned integration_point) const override
 Provides the shape matrix at the given integration point.
 
std::vector< double > const & getIntPtDarcyVelocity (const double t, std::vector< GlobalVector * > const &x, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_table, std::vector< double > &velocity_cache) const override
 
- Public Member Functions inherited from ProcessLib::LocalAssemblerInterface
virtual ~LocalAssemblerInterface ()=default
 
virtual void setInitialConditions (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, int const process_id)
 
virtual void initialize (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
 
virtual void preAssemble (double const, double const, std::vector< double > const &)
 
virtual void assembleForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
 
virtual void assembleWithJacobian (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > const &local_x_prev, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
 
virtual void assembleWithJacobianForStaggeredScheme (double const t, double const dt, Eigen::VectorXd const &local_x, Eigen::VectorXd const &local_x_prev, int const process_id, std::vector< double > &local_b_data, std::vector< double > &local_Jac_data)
 
virtual void computeSecondaryVariable (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
 
virtual void preTimestep (std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, double const t, double const delta_t)
 
virtual void postTimestep (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
void postNonLinearSolver (std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, double const t, double const dt, int const process_id)
 
virtual Eigen::Vector3d getFlux (MathLib::Point3d const &, double const, std::vector< std::vector< double > > const &) const
 Fits to staggered scheme.
 
virtual std::optional< VectorSegment > getVectorDeformationSegment () const
 
- Public Member Functions inherited from NumLib::ExtrapolatableElement
virtual ~ExtrapolatableElement ()=default
 

Private Types

using ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>
 
using ShapeMatrices = typename ShapeMatricesType::ShapeMatrices
 
using LocalAssemblerTraits
 
using NodalMatrixType = typename LocalAssemblerTraits::LocalMatrix
 
using NodalVectorType = typename LocalAssemblerTraits::LocalVector
 
using NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType
 
using GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType
 
using GlobalDimMatrixType = typename ShapeMatricesType::GlobalDimMatrixType
 
using GlobalDimNodalMatrixType
 

Private Member Functions

template<typename LaplacianGravityVelocityCalculator >
void assembleMatrixAndVector (double const t, double const dt, std::vector< double > const &local_x, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
 
template<typename LaplacianGravityVelocityCalculator , typename VelocityCacheType >
void computeProjectedDarcyVelocity (const double t, const double dt, std::vector< double > const &local_x, ParameterLib::SpatialPosition const &pos, VelocityCacheType &darcy_velocity_at_ips) const
 
template<typename VelocityCacheType >
void computeDarcyVelocity (bool const is_scalar_permeability, const double t, const double dt, std::vector< double > const &local_x, ParameterLib::SpatialPosition const &pos, VelocityCacheType &darcy_velocity_at_ips) const
 

Private Attributes

MeshLib::Element const & _element
 
NumLib::GenericIntegrationMethod const & _integration_method
 
std::vector< IntegrationPointData< GlobalDimNodalMatrixType > > _ip_data
 
const LiquidFlowData & _process_data
 

Member Typedef Documentation

◆ GlobalDimMatrixType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::GlobalDimMatrixType = typename ShapeMatricesType::GlobalDimMatrixType
private

Definition at line 79 of file LiquidFlowLocalAssembler.h.

◆ GlobalDimNodalMatrixType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::GlobalDimNodalMatrixType
private
Initial value:
typename ShapeMatricesType::GlobalDimNodalMatrixType
EigenFixedShapeMatrixPolicy::GlobalDimNodalMatrixType
MatrixType< GlobalDim, ShapeFunction::NPOINTS > GlobalDimNodalMatrixType
Definition ShapeMatrixPolicy.h:82

Definition at line 80 of file LiquidFlowLocalAssembler.h.

◆ GlobalDimVectorType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::GlobalDimVectorType = typename ShapeMatricesType::GlobalDimVectorType
private

Definition at line 78 of file LiquidFlowLocalAssembler.h.

◆ LocalAssemblerTraits

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LocalAssemblerTraits
private
Initial value:
ProcessLib::LocalAssemblerTraits<
ShapeMatricesType, ShapeFunction::NPOINTS, NUM_NODAL_DOF, GlobalDim>
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::ShapeMatricesType
ShapeMatrixPolicyType< ShapeFunction, GlobalDim > ShapeMatricesType
Definition LiquidFlowLocalAssembler.h:69

Definition at line 72 of file LiquidFlowLocalAssembler.h.

◆ NodalMatrixType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::NodalMatrixType = typename LocalAssemblerTraits::LocalMatrix
private

Definition at line 75 of file LiquidFlowLocalAssembler.h.

◆ NodalRowVectorType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType
private

Definition at line 77 of file LiquidFlowLocalAssembler.h.

◆ NodalVectorType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::NodalVectorType = typename LocalAssemblerTraits::LocalVector
private

Definition at line 76 of file LiquidFlowLocalAssembler.h.

◆ ShapeMatrices

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::ShapeMatrices = typename ShapeMatricesType::ShapeMatrices
private

Definition at line 70 of file LiquidFlowLocalAssembler.h.

◆ ShapeMatricesType

template<typename ShapeFunction , int GlobalDim>
using ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>
private

Definition at line 69 of file LiquidFlowLocalAssembler.h.

Constructor & Destructor Documentation

◆ LiquidFlowLocalAssembler()

template<typename ShapeFunction , int GlobalDim>
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LiquidFlowLocalAssembler ( MeshLib::Element const & element,
std::size_t const ,
NumLib::GenericIntegrationMethod const & integration_method,
bool const is_axially_symmetric,
LiquidFlowData const & process_data )
inline

Definition at line 84 of file LiquidFlowLocalAssembler.h.

90 : _element(element),
91 _integration_method(integration_method),
92 _process_data(process_data)
93 {
94 unsigned const n_integration_points =
95 _integration_method.getNumberOfPoints();
96 _ip_data.reserve(n_integration_points);
97
98 auto const& shape_matrices =
99 NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
100 GlobalDim>(element, is_axially_symmetric,
101 _integration_method);
102
103 ParameterLib::SpatialPosition pos;
104 pos.setElementID(_element.getID());
105
106 double const aperture_size =
107 (_element.getDimension() == 3u)
108 ? 1.0
109 : _process_data.aperture_size(0.0, pos)[0];
110
111 for (unsigned ip = 0; ip < n_integration_points; ip++)
112 {
113 _ip_data.emplace_back(
114 shape_matrices[ip].dNdx,
115 _integration_method.getWeightedPoint(ip).getWeight() *
116 shape_matrices[ip].integralMeasure *
117 shape_matrices[ip].detJ * aperture_size);
118 }
119 }
MathLib::WeightedPoint::getWeight
constexpr double getWeight() const
Definition WeightedPoint.h:82
MeshLib::Element::getDimension
virtual constexpr unsigned getDimension() const =0
Get dimension of the mesh element.
MeshLib::Element::getID
std::size_t getID() const
Returns the ID of the element.
Definition Element.h:91
NumLib::GenericIntegrationMethod::getWeightedPoint
MathLib::WeightedPoint const & getWeightedPoint(unsigned const igp) const
Definition GenericIntegrationMethod.h:36
ParameterLib::SpatialPosition::setElementID
void setElementID(std::size_t element_id)
Definition SpatialPosition.h:76
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::_ip_data
std::vector< IntegrationPointData< GlobalDimNodalMatrixType > > _ip_data
Definition LiquidFlowLocalAssembler.h:156
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::_integration_method
NumLib::GenericIntegrationMethod const & _integration_method
Definition LiquidFlowLocalAssembler.h:155
NumLib::initShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
Definition InitShapeMatrices.h:54

References ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_element, ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_integration_method, ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_ip_data, ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_process_data, ProcessLib::LiquidFlow::LiquidFlowData::aperture_size, MeshLib::Element::getDimension(), MeshLib::Element::getID(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), NumLib::initShapeMatrices(), and ParameterLib::SpatialPosition::setElementID().

Member Function Documentation

◆ assemble()

template<typename ShapeFunction , int GlobalDim>
void ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::assemble ( double const t,
double const dt,
std::vector< double > const & local_x,
std::vector< double > const & ,
std::vector< double > & local_M_data,
std::vector< double > & local_K_data,
std::vector< double > & local_b_data )
overridevirtual

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 27 of file LiquidFlowLocalAssembler-impl.h.

32{
33 // Dummy integration point
34 unsigned const ip = 0;
35 // auto const& ip_data = _ip_data[ip];
36 auto const& Ns = _process_data.shape_matrix_cache
37 .NsHigherOrder<typename ShapeFunction::MeshElement>();
38 auto const& N = Ns[ip];
39
40 ParameterLib::SpatialPosition const pos{
41 std::nullopt, _element.getID(),
42 MathLib::Point3d(
43 NumLib::interpolateCoordinates<ShapeFunction, ShapeMatricesType>(
44 _element, N))};
45
46 auto const& medium = *_process_data.media_map.getMedium(_element.getID());
47 MaterialPropertyLib::VariableArray vars;
48 vars.temperature =
49 medium[MaterialPropertyLib::PropertyType::reference_temperature]
50 .template value<double>(vars, pos, t, dt);
51 vars.liquid_phase_pressure = std::numeric_limits<double>::quiet_NaN();
52 GlobalDimMatrixType const permeability =
53 MaterialPropertyLib::formEigenTensor<GlobalDim>(
54 medium[MaterialPropertyLib::PropertyType::permeability].value(
55 vars, pos, t, dt));
56 // Note: For Inclined 1D in 2D/3D or 2D element in 3D, the first item in
57 // the assert must be changed to permeability.rows() ==
58 // _element->getDimension()
59 assert(permeability.rows() == GlobalDim || permeability.rows() == 1);
60
61 if (permeability.size() == 1)
62 { // isotropic or 1D problem.
63 assembleMatrixAndVector<IsotropicCalculator>(
64 t, dt, local_x, local_M_data, local_K_data, local_b_data);
65 }
66 else
67 {
68 assembleMatrixAndVector<AnisotropicCalculator>(
69 t, dt, local_x, local_M_data, local_K_data, local_b_data);
70 }
71}
NumLib::ShapeMatrixCache::NsHigherOrder
auto const & NsHigherOrder() const
Definition ShapeMatrixCache.h:99
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::assembleMatrixAndVector
void assembleMatrixAndVector(double const t, double const dt, std::vector< double > const &local_x, std::vector< double > &local_M_data, std::vector< double > &local_K_data, std::vector< double > &local_b_data)
Definition LiquidFlowLocalAssembler-impl.h:122
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::GlobalDimMatrixType
typename ShapeMatricesType::GlobalDimMatrixType GlobalDimMatrixType
Definition LiquidFlowLocalAssembler.h:79
MaterialPropertyLib::formEigenTensor
Eigen::Matrix< double, GlobalDim, GlobalDim > formEigenTensor(MaterialPropertyLib::PropertyDataType const &values)
Definition FormEigenTensor.cpp:126
NumLib::interpolateCoordinates
std::array< double, 3 > interpolateCoordinates(MeshLib::Element const &e, typename ShapeMatricesType::ShapeMatrices::ShapeType const &N)
Definition InitShapeMatrices.h:102
ProcessLib::LiquidFlow::LiquidFlowData::shape_matrix_cache
NumLib::ShapeMatrixCache shape_matrix_cache
caches for each mesh element type the shape matrix
Definition LiquidFlowData.h:57
ProcessLib::LiquidFlow::LiquidFlowData::media_map
MaterialPropertyLib::MaterialSpatialDistributionMap media_map
Definition LiquidFlowData.h:42

References MaterialPropertyLib::formEigenTensor(), NumLib::interpolateCoordinates(), MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::permeability, MaterialPropertyLib::reference_temperature, and MaterialPropertyLib::VariableArray::temperature.

◆ assembleMatrixAndVector()

template<typename ShapeFunction , int GlobalDim>
template<typename LaplacianGravityVelocityCalculator >
void ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::assembleMatrixAndVector ( double const t,
double const dt,
std::vector< double > const & local_x,
std::vector< double > & local_M_data,
std::vector< double > & local_K_data,
std::vector< double > & local_b_data )
private

Definition at line 121 of file LiquidFlowLocalAssembler-impl.h.

127{
128 auto const local_matrix_size = local_x.size();
129 assert(local_matrix_size == ShapeFunction::NPOINTS);
130
131 auto local_M = MathLib::createZeroedMatrix<NodalMatrixType>(
132 local_M_data, local_matrix_size, local_matrix_size);
133 auto local_K = MathLib::createZeroedMatrix<NodalMatrixType>(
134 local_K_data, local_matrix_size, local_matrix_size);
135 auto local_b = MathLib::createZeroedVector<NodalVectorType>(
136 local_b_data, local_matrix_size);
137 const auto local_p_vec =
138 MathLib::toVector<NodalVectorType>(local_x, local_matrix_size);
139
140 unsigned const n_integration_points =
141 _integration_method.getNumberOfPoints();
142
143 ParameterLib::SpatialPosition pos;
144 pos.setElementID(_element.getID());
145
146 auto const& medium = *_process_data.media_map.getMedium(_element.getID());
147 auto const& fluid_phase = fluidPhase(medium);
148
149 MaterialPropertyLib::VariableArray vars;
150 auto& phase_pressure = medium.hasPhase("Gas") ? vars.gas_phase_pressure
151 : vars.liquid_phase_pressure;
152
153 GlobalDimVectorType const projected_body_force_vector =
154 _process_data.element_rotation_matrices[_element.getID()] *
155 _process_data.element_rotation_matrices[_element.getID()].transpose() *
156 _process_data.specific_body_force;
157
158 auto const& Ns = _process_data.shape_matrix_cache
159 .NsHigherOrder<typename ShapeFunction::MeshElement>();
160
161 for (unsigned ip = 0; ip < n_integration_points; ip++)
162 {
163 auto const& ip_data = _ip_data[ip];
164 auto const& N = Ns[ip];
165
166 phase_pressure = N.dot(local_p_vec);
167 ParameterLib::SpatialPosition const pos{
168 std::nullopt, _element.getID(),
169 MathLib::Point3d(NumLib::interpolateCoordinates<ShapeFunction,
170 ShapeMatricesType>(
171 _element, N))};
172
173 vars.temperature =
174 medium[MaterialPropertyLib::PropertyType::reference_temperature]
175 .template value<double>(vars, pos, t, dt);
176
177 auto const [fluid_density, viscosity] =
178 MaterialPropertyLib::getFluidDensityAndViscosity(t, dt, pos,
179 fluid_phase, vars);
180
181 auto const porosity =
182 medium[MaterialPropertyLib::PropertyType::porosity]
183 .template value<double>(vars, pos, t, dt);
184 auto const specific_storage =
185 medium[MaterialPropertyLib::PropertyType::storage]
186 .template value<double>(vars, pos, t, dt);
187
188 auto const get_drho_dp = [&]()
189 {
190 return fluid_phase[MaterialPropertyLib::PropertyType::density]
191 .template dValue<double>(vars, _process_data.phase_variable,
192 pos, t, dt);
193 };
194 auto const storage =
195 _process_data.equation_balance_type == EquationBalanceType::volume
196 ? specific_storage
197 : specific_storage + porosity * get_drho_dp() / fluid_density;
198
199 double const scaling_factor =
200 _process_data.equation_balance_type == EquationBalanceType::volume
201 ? 1.0
202 : fluid_density;
203 // Assemble mass matrix, M
204 local_M.noalias() += scaling_factor * storage * N.transpose() * N *
205 ip_data.integration_weight;
206
207 GlobalDimMatrixType const permeability =
208 MaterialPropertyLib::formEigenTensor<GlobalDim>(
209 medium[MaterialPropertyLib::PropertyType::permeability].value(
210 vars, pos, t, dt));
211
212 // Assemble Laplacian, K, and RHS by the gravitational term
213 LaplacianGravityVelocityCalculator::calculateLaplacianAndGravityTerm(
214 local_K, local_b, ip_data, scaling_factor * permeability, viscosity,
215 fluid_density, projected_body_force_vector,
216 _process_data.has_gravity);
217 }
218}
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::GlobalDimVectorType
typename ShapeMatricesType::GlobalDimVectorType GlobalDimVectorType
Definition LiquidFlowLocalAssembler.h:78
MaterialPropertyLib::getFluidDensityAndViscosity
std::tuple< double, double > getFluidDensityAndViscosity(double const t, double const dt, ParameterLib::SpatialPosition const &pos, MaterialPropertyLib::Phase const &fluid_phase, MaterialPropertyLib::VariableArray &vars)
It computes fluid density and viscosity for single phase flow model.
Definition GetFluidDensityAndViscosity.cpp:42
MaterialPropertyLib::fluidPhase
Phase const & fluidPhase(Medium const &medium)
Returns a gas or aqueous liquid phase of the given medium.
Definition Medium.cpp:100
MathLib::createZeroedVector
Eigen::Map< Vector > createZeroedVector(std::vector< double > &data, Eigen::VectorXd::Index size)
Definition EigenMapTools.h:153
MathLib::toVector
Eigen::Map< const Vector > toVector(std::vector< double > const &data, Eigen::VectorXd::Index size)
Creates an Eigen mapped vector from the given data vector.
Definition EigenMapTools.h:167
MathLib::createZeroedMatrix
Eigen::Map< Matrix > createZeroedMatrix(std::vector< double > &data, Eigen::MatrixXd::Index rows, Eigen::MatrixXd::Index cols)
Definition EigenMapTools.h:32
ProcessLib::LiquidFlow::LiquidFlowData::phase_variable
MaterialPropertyLib::Variable const phase_variable
Definition LiquidFlowData.h:59
ProcessLib::LiquidFlow::LiquidFlowData::element_rotation_matrices
std::vector< Eigen::MatrixXd > const element_rotation_matrices
A vector of the rotation matrices for all elements.
Definition LiquidFlowData.h:45
ProcessLib::LiquidFlow::LiquidFlowData::specific_body_force
Eigen::VectorXd const specific_body_force
Definition LiquidFlowData.h:49
ProcessLib::LiquidFlow::LiquidFlowData::equation_balance_type
EquationBalanceType const equation_balance_type
Definition LiquidFlowData.h:40

References MathLib::createZeroedMatrix(), MathLib::createZeroedVector(), MaterialPropertyLib::density, MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::VariableArray::gas_phase_pressure, MaterialPropertyLib::getFluidDensityAndViscosity(), NumLib::interpolateCoordinates(), MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::permeability, MaterialPropertyLib::porosity, MaterialPropertyLib::reference_temperature, ParameterLib::SpatialPosition::setElementID(), MaterialPropertyLib::storage, MaterialPropertyLib::VariableArray::temperature, MathLib::toVector(), and ProcessLib::LiquidFlow::volume.

◆ computeDarcyVelocity()

template<typename ShapeFunction , int GlobalDim>
template<typename VelocityCacheType >
void ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::computeDarcyVelocity ( bool const is_scalar_permeability,
const double t,
const double dt,
std::vector< double > const & local_x,
ParameterLib::SpatialPosition const & pos,
VelocityCacheType & darcy_velocity_at_ips ) const
private

Definition at line 222 of file LiquidFlowLocalAssembler-impl.h.

227{
228 if (is_scalar_permeability)
229 { // isotropic or 1D problem.
230 computeProjectedDarcyVelocity<IsotropicCalculator>(
231 t, dt, local_x, pos, darcy_velocity_at_ips);
232 }
233 else
234 {
235 computeProjectedDarcyVelocity<AnisotropicCalculator>(
236 t, dt, local_x, pos, darcy_velocity_at_ips);
237 }
238}
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::computeProjectedDarcyVelocity
void computeProjectedDarcyVelocity(const double t, const double dt, std::vector< double > const &local_x, ParameterLib::SpatialPosition const &pos, VelocityCacheType &darcy_velocity_at_ips) const
Definition LiquidFlowLocalAssembler-impl.h:302

◆ computeProjectedDarcyVelocity()

template<typename ShapeFunction , int GlobalDim>
template<typename LaplacianGravityVelocityCalculator , typename VelocityCacheType >
void ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::computeProjectedDarcyVelocity ( const double t,
const double dt,
std::vector< double > const & local_x,
ParameterLib::SpatialPosition const & pos,
VelocityCacheType & darcy_velocity_at_ips ) const
private

Definition at line 301 of file LiquidFlowLocalAssembler-impl.h.

306{
307 auto const local_matrix_size = local_x.size();
308 assert(local_matrix_size == ShapeFunction::NPOINTS);
309
310 const auto local_p_vec =
311 MathLib::toVector<NodalVectorType>(local_x, local_matrix_size);
312
313 unsigned const n_integration_points =
314 _integration_method.getNumberOfPoints();
315
316 auto const& medium = *_process_data.media_map.getMedium(_element.getID());
317 auto const& fluid_phase = fluidPhase(medium);
318
319 MaterialPropertyLib::VariableArray vars;
320 auto& phase_pressure = medium.hasPhase("Gas") ? vars.gas_phase_pressure
321 : vars.liquid_phase_pressure;
322
323 GlobalDimVectorType const projected_body_force_vector =
324 _process_data.element_rotation_matrices[_element.getID()] *
325 _process_data.element_rotation_matrices[_element.getID()].transpose() *
326 _process_data.specific_body_force;
327
328 auto const& Ns = _process_data.shape_matrix_cache
329 .NsHigherOrder<typename ShapeFunction::MeshElement>();
330
331 for (unsigned ip = 0; ip < n_integration_points; ip++)
332 {
333 auto const& ip_data = _ip_data[ip];
334 auto const& N = Ns[ip];
335
336 phase_pressure = N.dot(local_p_vec);
337 ParameterLib::SpatialPosition const pos{
338 std::nullopt, _element.getID(),
339 MathLib::Point3d(NumLib::interpolateCoordinates<ShapeFunction,
340 ShapeMatricesType>(
341 _element, N))};
342
343 vars.temperature =
344 medium[MaterialPropertyLib::PropertyType::reference_temperature]
345 .template value<double>(vars, pos, t, dt);
346 auto const [fluid_density, viscosity] =
347 MaterialPropertyLib::getFluidDensityAndViscosity(t, dt, pos,
348 fluid_phase, vars);
349
350 GlobalDimMatrixType const permeability =
351 MaterialPropertyLib::formEigenTensor<GlobalDim>(
352 medium[MaterialPropertyLib::PropertyType::permeability].value(
353 vars, pos, t, dt));
354
355 darcy_velocity_at_ips.col(ip) =
356 LaplacianGravityVelocityCalculator::calculateVelocity(
357 local_p_vec, ip_data, permeability, viscosity, fluid_density,
358 projected_body_force_vector, _process_data.has_gravity);
359 }
360}

References MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::VariableArray::gas_phase_pressure, MaterialPropertyLib::getFluidDensityAndViscosity(), NumLib::interpolateCoordinates(), MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::permeability, MaterialPropertyLib::reference_temperature, MaterialPropertyLib::VariableArray::temperature, and MathLib::toVector().

◆ getFlux()

template<typename ShapeFunction , int GlobalDim>
Eigen::Vector3d ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::getFlux ( MathLib::Point3d const & p_local_coords,
double const t,
std::vector< double > const & local_x ) const
overridevirtual

Computes the flux in the point p_local_coords that is given in local coordinates using the values from local_x.

Reimplemented from ProcessLib::LocalAssemblerInterface.

Definition at line 74 of file LiquidFlowLocalAssembler-impl.h.

77{
78 // TODO (tf) Temporary value not used by current material models. Need
79 // extension of getFlux interface
80 double const dt = std::numeric_limits<double>::quiet_NaN();
81
82 // Note: Axial symmetry is set to false here, because we only need dNdx
83 // here, which is not affected by axial symmetry.
84 auto const shape_matrices =
85 NumLib::computeShapeMatrices<ShapeFunction, ShapeMatricesType,
86 GlobalDim>(_element,
87 false /*is_axially_symmetric*/,
88 std::array{p_local_coords})[0];
89
90 // create pos object to access the correct media property
91 ParameterLib::SpatialPosition pos;
92 pos.setElementID(_element.getID());
93
94 auto const& medium = *_process_data.media_map.getMedium(_element.getID());
95 auto const& fluid_phase = fluidPhase(medium);
96
97 MaterialPropertyLib::VariableArray vars;
98
99 double pressure = 0.0;
100 NumLib::shapeFunctionInterpolate(local_x, shape_matrices.N, pressure);
101 vars.liquid_phase_pressure = pressure;
102
103 GlobalDimMatrixType const intrinsic_permeability =
104 MaterialPropertyLib::formEigenTensor<GlobalDim>(
105 medium[MaterialPropertyLib::PropertyType::permeability].value(
106 vars, pos, t, dt));
107 auto const viscosity =
108 fluid_phase[MaterialPropertyLib::PropertyType::viscosity]
109 .template value<double>(vars, pos, t, dt);
110
111 Eigen::Vector3d flux(0.0, 0.0, 0.0);
112 flux.head<GlobalDim>() =
113 -intrinsic_permeability / viscosity * shape_matrices.dNdx *
114 Eigen::Map<const NodalVectorType>(local_x.data(), local_x.size());
115
116 return flux;
117}
NumLib::detail::shapeFunctionInterpolate
void shapeFunctionInterpolate(const NodalValues &, const ShapeMatrix &)
Definition Interpolation.h:27
NumLib::computeShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > computeShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, PointContainer const &points)
Definition InitShapeMatrices.h:26

References NumLib::computeShapeMatrices(), MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::permeability, ParameterLib::SpatialPosition::setElementID(), NumLib::detail::shapeFunctionInterpolate(), and MaterialPropertyLib::viscosity.

◆ getIntPtDarcyVelocity()

template<typename ShapeFunction , int GlobalDim>
std::vector< double > const & ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::getIntPtDarcyVelocity ( const double t,
std::vector< GlobalVector * > const & x,
std::vector< NumLib::LocalToGlobalIndexMap const * > const & dof_table,
std::vector< double > & velocity_cache ) const
overridevirtual

Implements ProcessLib::LiquidFlow::LiquidFlowLocalAssemblerInterface.

Definition at line 242 of file LiquidFlowLocalAssembler-impl.h.

247{
248 // TODO (tf) Temporary value not used by current material models. Need
249 // extension of secondary variable interface.
250 double const dt = std::numeric_limits<double>::quiet_NaN();
251
252 constexpr int process_id = 0;
253 auto const indices =
254 NumLib::getIndices(_element.getID(), *dof_tables[process_id]);
255 auto const local_x = x[process_id]->get(indices);
256 auto const n_integration_points = _integration_method.getNumberOfPoints();
257 velocity_cache.clear();
258
259 // Dummy integration point
260 unsigned const ip = 0;
261 // auto const& ip_data = _ip_data[ip];
262 auto const& Ns = _process_data.shape_matrix_cache
263 .NsHigherOrder<typename ShapeFunction::MeshElement>();
264 auto const& N = Ns[ip];
265
266 ParameterLib::SpatialPosition const pos{
267 std::nullopt, _element.getID(),
268 MathLib::Point3d(
269 NumLib::interpolateCoordinates<ShapeFunction, ShapeMatricesType>(
270 _element, N))};
271
272 auto const& medium = *_process_data.media_map.getMedium(_element.getID());
273 MaterialPropertyLib::VariableArray vars;
274 vars.temperature =
275 medium[MaterialPropertyLib::PropertyType::reference_temperature]
276 .template value<double>(vars, pos, t, dt);
277 vars.liquid_phase_pressure = std::numeric_limits<double>::quiet_NaN();
278
279 GlobalDimMatrixType const permeability =
280 MaterialPropertyLib::formEigenTensor<GlobalDim>(
281 medium[MaterialPropertyLib::PropertyType::permeability].value(
282 vars, pos, t, dt));
283
284 assert(permeability.rows() == GlobalDim || permeability.rows() == 1);
285
286 bool const is_scalar_permeability = (permeability.size() == 1);
287
288 auto velocity_cache_vectors = MathLib::createZeroedMatrix<
289 Eigen::Matrix<double, GlobalDim, Eigen::Dynamic, Eigen::RowMajor>>(
290 velocity_cache, GlobalDim, n_integration_points);
291
292 computeDarcyVelocity(is_scalar_permeability, t, dt, local_x, pos,
293 velocity_cache_vectors);
294
295 return velocity_cache;
296}
ProcessLib::LiquidFlow::LiquidFlowLocalAssembler::computeDarcyVelocity
void computeDarcyVelocity(bool const is_scalar_permeability, const double t, const double dt, std::vector< double > const &local_x, ParameterLib::SpatialPosition const &pos, VelocityCacheType &darcy_velocity_at_ips) const
Definition LiquidFlowLocalAssembler-impl.h:222
NumLib::getIndices
std::vector< GlobalIndexType > getIndices(std::size_t const mesh_item_id, NumLib::LocalToGlobalIndexMap const &dof_table)
Definition DOFTableUtil.cpp:112

References MathLib::createZeroedMatrix(), MaterialPropertyLib::formEigenTensor(), NumLib::getIndices(), NumLib::interpolateCoordinates(), MaterialPropertyLib::VariableArray::liquid_phase_pressure, MaterialPropertyLib::permeability, MaterialPropertyLib::reference_temperature, and MaterialPropertyLib::VariableArray::temperature.

◆ getShapeMatrix()

template<typename ShapeFunction , int GlobalDim>
Eigen::Map< const Eigen::RowVectorXd > ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::getShapeMatrix ( const unsigned integration_point) const
inlineoverridevirtual

Provides the shape matrix at the given integration point.

Implements NumLib::ExtrapolatableElement.

Definition at line 134 of file LiquidFlowLocalAssembler.h.

136 {
137 auto const& N =
138 _process_data.shape_matrix_cache
139 .NsHigherOrder<typename ShapeFunction::MeshElement>();
140
141 // assumes N is stored contiguously in memory
142 return Eigen::Map<const Eigen::RowVectorXd>(
143 N[integration_point].data(), N[integration_point].size());
144 }
MathLib::VectorizedTensor::size
constexpr int size(int const displacement_dim)
Vectorized tensor size for given displacement dimension.
Definition VectorizedTensor.h:25

References ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_process_data, NumLib::ShapeMatrixCache::NsHigherOrder(), and ProcessLib::LiquidFlow::LiquidFlowData::shape_matrix_cache.

Member Data Documentation

◆ _element

template<typename ShapeFunction , int GlobalDim>
MeshLib::Element const& ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_element
private

Definition at line 153 of file LiquidFlowLocalAssembler.h.

Referenced by ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LiquidFlowLocalAssembler().

◆ _integration_method

template<typename ShapeFunction , int GlobalDim>
NumLib::GenericIntegrationMethod const& ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_integration_method
private

Definition at line 155 of file LiquidFlowLocalAssembler.h.

Referenced by ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LiquidFlowLocalAssembler().

◆ _ip_data

template<typename ShapeFunction , int GlobalDim>
std::vector<IntegrationPointData<GlobalDimNodalMatrixType> > ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_ip_data
private

Definition at line 156 of file LiquidFlowLocalAssembler.h.

Referenced by ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LiquidFlowLocalAssembler().

◆ _process_data

template<typename ShapeFunction , int GlobalDim>
const LiquidFlowData& ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::_process_data
private

Definition at line 225 of file LiquidFlowLocalAssembler.h.

Referenced by ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::LiquidFlowLocalAssembler(), and ProcessLib::LiquidFlow::LiquidFlowLocalAssembler< ShapeFunction, GlobalDim >::getShapeMatrix().

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