Dual flow hydraulic 2D & 3D for Quad, Tri, Hex and Tet. More...

#include <gmpHydraulicDualFlow.h>

Inheritance diagram for GmpHydraulicDualFlow:

[legend]

Collaboration diagram for GmpHydraulicDualFlow:

[legend]

Public Member Functions
	GmpHydraulicDualFlow (const char pluginType, GmSimulationData simulation, QString id, QString description, const GmpFemPhysicsCommonMaterialFactory *matFactory, const GmLogCategory &logger)

virtual	~GmpHydraulicDualFlow ()
	Destructor.

virtual const QVariantMap *	physicsMetaDataMap ()
	Returns a reference for the Hydraulics dual flow attribute map, built when the function is called for the first time.

virtual GmpFemPhysics::FemResultType	fillElementData (const GmElement *e, GmpFemMatrixSet &elemMatrices, GmpFemVectorSet &elemVectors)
	Fills the element matrices and vectors for a generic dual flow element.

virtual bool	fixedNodalForcesBc (QVector< int > &nodes, QVector< int > &dof, QVector< double > &values) const
	See comments on base class. Fills vectors with prescribed nodal forces - flow.

virtual bool	fixedNodalDofsBc (QVector< int > &nodes, QVector< int > &dof, QVector< double > &values, bool *constantValues) const
	See comments on base class. Fills vectors with prescribed node pore pressures.

virtual void	fillFractureStorageCompMatrix (const GmElement e, const GmVector coord, int ip, GmMatrix &Sfc, const GmVector &N, double c) const
	Calculates the stiffness matrix of "storage / compressibility of fluid" for element 'e' at the given integration point (coord, ip) storing the result at the Sfc matrix. N is the shape function vector and c a constant equal to w * h * detJ.

virtual void	fillTransferMatrix (const GmElement e, const GmVector coord, int ip, GmMatrix &Wfm, const GmVector &N, double c) const
	Fills the Transfer matrix with the computed scalar "shape factor" Returns Wfm = ws * c * N * N.t, where: ws is the shape factor.

virtual void	fillFluidFlowMatrix (const GmElement e, const GmVector coord, int ip, GmMatrix &Kfc, const GmMatrix &Bp, double c) const
	Calculates the flow stiffness matrix of fluid 2 due to permeability contribution or permeability matrix for element 'e' at the given integration point (coord, ip) storing the result at the Kfc matrix. Bp is the strain pore pressure matrix and c a constant equal to w * h * detJ.

virtual void	fillElementFluidPorePressure (const GmElement *e, GmVector &Pfe)
	Given an element, fills the vector pfe with nodal pore pressureof the seconf fluid. The vector should have size equal to the number of calculated nodes.

virtual void	fillDampingMatrix (const GmElement *e, GmMatrix &C, const GmMatrix &Sc, const GmMatrix &Sfc, const QVector< int > &iPe, const QVector< int > &iPfe) const
	Fills the Damping matrix with the computed "storage/compressibility" Matrizes Returns C = Sc + Sfc, where: Sc is the "storage/compressibility" Matrix of the medium Sfc is the "storage/compressibility" Matrix of the fluid.

virtual void	fillStiffnessMatrix (const GmElement *e, GmMatrix &K, const GmMatrix &Kc, const GmMatrix &Kfc, const GmMatrix &Wfm, const QVector< int > &iPe, const QVector< int > &iPfe) const
	Fills the Stiffness matrix with the computed "permeability" Matrizes Returns K = Kc + Kfc, where: Kc is the permeability Matrix of the medium Kfc is the permeability Matrix of the fluid.

virtual void	fillInternalForceVector (const GmElement *e, GmVector &Fi, const GmMatrix &Kc, const GmMatrix &Kfc, const GmMatrix &Wfm, const GmVector &pe, const GmVector &pfe, const QVector< int > &iPe, const QVector< int > &iPfe) const
	Fills the Internal force vector for both fluids Returns fi = Kcpe + Kfcpfe, where: Kc is the permeability Matrix of the medium Kfc is the permeability Matrix of the fluid.

virtual void	FillDofIndexPosition (const GmElement *e, QVector< int > &iPe, QVector< int > &iPfe) const
	Fills the position index of the two hydraulic degrees of freedom. the implementation supports order of the defined state variables i.e {P, Pf} or {Pf, P}.

virtual bool	calcDerivedResults (bool nonLinearSolver)
	Calc velocity and flow on nodes and/or Gauss points.

virtual GmPropertySet *	fracturePropertySet (int id) const

Public Member Functions inherited from GmpHydraulic
	GmpHydraulic (const char pluginType, GmSimulationData simulation, QString id, QString description, const GmpFemPhysicsCommonMaterialFactory *matFactory, const GmLogCategory &logger)
	Constructor. Will be called by the plugin loading code.

virtual	~GmpHydraulic ()
	Destructor.

virtual const char *	pluginName () const

virtual const char *	pluginType () const

virtual bool	supportsStateDumping ()

virtual FemResultType	fillElementDataForBc (const GmElement e, const GmBoundaryCondition bc, int bcIndex, int bcListIndex, int border, GmpFemMatrixSet &elemMatrices, GmpFemVectorSet &elemVectors)

virtual FemResultType	fillElementDataBc (const GmElement e, const GmBoundaryCondition bc, int bcIndex, int bcListIndex, int border, GmMatrix &elemMat, GmVector &elemVec)

virtual void	fillStorageCompMatrix (const GmElement e, const GmVector coord, int ip, GmMatrix &Sc, const GmVector &N, double c)
	Calculates the stiffness matrix of "storage / compressibility" for element 'e' at the given integration point (coord, ip) storing the result at the Sc matrix. N is the shape function vector and c a constant equal to w * h * detJ.

virtual void	fillFlowMatrix (const GmElement e, const GmVector coord, int ip, GmMatrix &Kc, const GmMatrix &Bp, double c) const
	Calculates the flow stiffness matrix due to permeability contribution or permeability matrix for element 'e' at the given integration point (coord, ip) storing the result at the Kc matrix. Bp is the strain pore pressure matrix and c a constant equal to w * h * detJ.

virtual void	fillFluidVelocity (const GmElement e, const GmVector coord, int ip, GmVector &Vip, const GmVector &Pe, const GmMatrix &Bp) const
	Evaluates velocity at integration points (it should be improved to consider gravity effect)

virtual void	fillElementPorePressure (const GmElement *e, GmVector &pe)
	Given an element, fills the vector pe with nodal pore pressure. The vector should have size equal to the number of calculated nodes.

virtual void	gravitySupplyFluid (const GmElement e, const GmVector coord, int ip, GmVector &qext_gsf, GmMatrix &Bp, double c)

virtual void	fillSourceTermVector (const GmElement e, const GmVector coord, int ip, GmVector &St, const GmVector &N, double c)
	Calculates the source term due the porosity rate of porous eskeleton for element 'e' at the given integration point (coord, ip) storing the result at external vector. N is the shape function vector and c a constant equal to w * h * detJ.

virtual double	fillBpMatrix (const GmShape *shape, const GmVector &ncoord, const GmMatrix &X, const GmVector &N, const GmMatrix &J, GmMatrix &Bp)=0
	Given a point 'ncoord', the element coordinate matrix 'X', the point shape functions 'N' and the transformation Jacobian 'J', fills the matrix with cartesian partial derivatives Bp using the provided shape function and returns the scaled jacobian determinant.

virtual void	gravityVector (const GmElement e, const GmVector coord, int ip, GmVector &grav) const
	Returns the gravity accelaration vector.

virtual bool	isAxisymmetric ()
	Returns TRUE only for axisymmetric models.

virtual double	axisymmetricFactor (const GmElement *e, const GmMatrix &X, const GmVector &N) const
	Returns the factor that should be applied when using axis symetric elements. 1.0 for other etypes.

HMaterialPointAccessors &	hydraulicMaterialPoint ()
	Returns the hydraulic material point accessor.

Protected Types
enum	DualFlowStateVarIds { PF_ID = GmpHydraulic::NUM_STATEVAR_IDS , NUM_STATEVAR_IDS }
	IDs fo physics state vars. More...

enum	DualFlowNodeAttributeIds { VFF_NA_ID = GmpHydraulic::NUM_NA_IDS , VEF_NA_ID , NUM_NA_IDS }
	IDs for physics node attributes. More...

enum	DualFlowGaussAttributeIds { VEF_GA_ID = GmpHydraulic::NUM_GA_IDS , NUM_GA_IDS }
	IDs for physics Gauss attributes. More...

enum	DualFlowBoundaryConditionIds { FIXED_PF_BC_ID = GmpHydraulic::NUM_BC_IDS , FIXED_QF_BC_ID , NUM_BC_IDS }
	IDs for accepted boundary condition types. More...

enum	DualFlowBoundaryConditionValueIds { BC_PF_ID = GmpHydraulic::NUM_BCV_IDS , BC_QWF_ID , NUM_BCV_IDS }
	IDs for property values from accepted boundary condition types. More...

enum	DualFlowPhysicsAttributeIds { FRACPROPSET_ID = GmpHydraulic::NUM_ATTRI_IDS , NUM_ATTRI_IDS }
	IDs for physics attributes. More...

Protected Member Functions
virtual const GmElementDof *	dofMapping (GmCellType type) const

virtual bool	checkAndLoadPrivateData (LuaTable &table)

virtual void	calcElementDerivedResultAtPoints (const GmElement *e, int resultId, const GmMatrix &evalPoints, bool evalAtIp, GmMatrix &result)
	The virtual function used by calcDerivedResults() to effectivelly calculate the velocity (resultId = 0) or the flow (resultId = 1) at the given evaluation points.

void	calcElementVelocityAtPoints (const GmElement *e, const GmMatrix &evalPoints, GmMatrix &v)
	The function responsible for effectivelly calculating the velocity at the given evaluation points, filling v.

void	calcElementFlowAtPoints (const GmElement *e, const GmMatrix &evalPoints, bool evalAtIp, GmMatrix &q)
	The function responsible for effectivelly calculating the flow at the given evaluation points, filling q.

Protected Member Functions inherited from GmpHydraulic
virtual bool	checkAndLoadAttributeAccessors (LuaTable &nodeTable, LuaTable &cellTable, LuaTable &gaussTable)
	Checks the loaded data.

void	calcElementVelocityAtPoints (const GmElement *e, const GmMatrix &evalPoints, GmMatrix &v)
	The function responsible for effectivelly calculating the velocity at the given evaluation points, filling v.

void	calcElementFlowAtPoints (const GmElement *e, const GmMatrix &evalPoints, bool evalAtIp, GmMatrix &q)
	The function responsible for effectivelly calculating the flow at the given evaluation points, filling q.

Protected Attributes
GmpFemPhysicsData	_fractureData
	fracture data

GmValueAccessor *	_anAcc
	fracture aperture accessor

Protected Attributes inherited from GmpHydraulic
HMaterialPointAccessors	_hmaterialPointAccessor
	Access to hydraulic material point.

bool	_setInitialCond
	Sets the initial conditions required by the physics.

Additional Inherited Members
Public Types inherited from GmpHydraulic
enum	HMaterialPointStates { saturated , unsaturated }

enum	PhysicsAttributeIds { ISOPARAMETRIC_ID , POROSITY_UPDATE_ID , PERMEABILITY_UPDATE_ID , USE_CONSISTENT_MATRIX_ID , VELOCITY_MODE_ID , FLOW_MODE_ID , DIRECT_ID , CGRAV_ID , VF_MODE_ID , NUM_ATTRI_IDS }
	IDs for physics attributes. More...

enum	StateVarIds { P_ID , NUM_STATEVAR_IDS }
	IDs fo physics state vars. More...

enum	NodeAttributeIds { V_NA_ID , Q_NA_ID , VF_NA_ID , NUM_NA_IDS }
	IDs for physics node attributes. More...

enum	cellAttributeIds { NUM_CA_IDS }
	IDs for physics Cell attributes. More...

enum	GaussAttributeIds { V_GA_ID , Q_GA_ID , FT_GA_ID , FTOLD_GA_ID , NUM_GA_IDS }
	IDs for physics Gauss attributes. More...

enum	ElementPropertyIds { THICKNESS_ID , GRAV_ID , SOURTERM_ID , NUM_PROPER_IDS }
	IDs for physics element properties. More...

enum	BoundaryConditionIds { FIXED_P_BC_ID , FIXED_Q_BC_ID , FIXED_SPF_BC_ID , NUM_BC_IDS }
	IDs for accepted boundary condition types. More...

enum	BoundaryConditionValueIds { BC_P_ID , BC_QW_ID , BC_IF_ID , NUM_BCV_IDS }
	IDs for property values from accepted boundary condition types. More...

Detailed Description

Dual flow hydraulic 2D & 3D for Quad, Tri, Hex and Tet.

Member Enumeration Documentation

◆ DualFlowBoundaryConditionIds

enum GmpHydraulicDualFlow::DualFlowBoundaryConditionIds

protected

IDs for accepted boundary condition types.

Enumerator
FIXED_PF_BC_ID	Id for fixed node fracture pore pressure boundary condition.
FIXED_QF_BC_ID	Id for fixed node fracture flow boundary condition.
NUM_BC_IDS	The number of boundary conditions ids above.

◆ DualFlowBoundaryConditionValueIds

enum GmpHydraulicDualFlow::DualFlowBoundaryConditionValueIds

protected

IDs for property values from accepted boundary condition types.

Enumerator
BC_PF_ID	Nodal pore pressure for fixed node fracture pore pressure boundary condition.
BC_QWF_ID	Nodal pore flow for fixed node fracture pore flow boundary condition.
NUM_BCV_IDS	The number of boundary conditions value ids above.

◆ DualFlowGaussAttributeIds

enum GmpHydraulicDualFlow::DualFlowGaussAttributeIds

protected

IDs for physics Gauss attributes.

Enumerator
VEF_GA_ID	Base Id for Gauss attribute(s) used to store the calculated fracture velocity.
NUM_GA_IDS	The number of Gauss attribute ids above.

◆ DualFlowNodeAttributeIds

enum GmpHydraulicDualFlow::DualFlowNodeAttributeIds

protected

IDs for physics node attributes.

Enumerator
VFF_NA_ID	Base Id for node attribute(s) used to store the calculated reaction volume flux.
VEF_NA_ID	Base Id for node attribute(s) used to store the calculated fracture velocity.
NUM_NA_IDS	The number of node attribute ids above.

◆ DualFlowPhysicsAttributeIds

enum GmpHydraulicDualFlow::DualFlowPhysicsAttributeIds

protected

IDs for physics attributes.

Enumerator
FRACPROPSET_ID	Id for retrieving the fracture property attribute.
NUM_ATTRI_IDS	The number of physics attribute ids above.

◆ DualFlowStateVarIds

enum GmpHydraulicDualFlow::DualFlowStateVarIds

protected

IDs fo physics state vars.

Enumerator
PF_ID	Id for retrieving the accessor to the fracture pore pressure state var.
NUM_STATEVAR_IDS	The number of property ids above.

Member Function Documentation

◆ calcDerivedResults()

bool GmpHydraulicDualFlow::calcDerivedResults ( bool nonLinearSolver )

virtual

Calc velocity and flow on nodes and/or Gauss points.

Reimplemented from GmpHydraulic.

◆ calcElementDerivedResultAtPoints()

void GmpHydraulicDualFlow::calcElementDerivedResultAtPoints	(	const GmElement *	e,
		int	resultId,
		const GmMatrix &	evalPoints,
		bool	evalAtIp,
		GmMatrix &	result )

protectedvirtual

The virtual function used by calcDerivedResults() to effectivelly calculate the velocity (resultId = 0) or the flow (resultId = 1) at the given evaluation points.

See the base class for a parameter description.

Reimplemented from GmpHydraulic.

◆ calcElementFlowAtPoints()

void GmpHydraulicDualFlow::calcElementFlowAtPoints	(	const GmElement *	e,
		const GmMatrix &	evalPoints,
		bool	evalAtIp,
		GmMatrix &	q )

protected

The function responsible for effectivelly calculating the flow at the given evaluation points, filling q.

This function result is a matrix with size d x m where 'm' is the number of evaluation points and 'd' the result dimension (2 or 3 for representing qxx, qyy and qzz)

Parameters

e	The element where we are calculating the result matrix
evalPoints	Matrix filled with natural coordinates for the 'm' points where results will be evaluated. Size = natural_coord_dimension x m.
evalAtIp	Flag that when set to true means that column 'i' in evalPoints is the coordinate of the 'i'th integration point.
q	Calculated results for each of the evaluation points. For each point of the 'm' points, a matrix column should be filled with the given point result. Size = d x m, where 'd' is the result dimension.

◆ calcElementVelocityAtPoints()

void GmpHydraulicDualFlow::calcElementVelocityAtPoints	(	const GmElement *	e,
		const GmMatrix &	evalPoints,
		GmMatrix &	v )

protected

The function responsible for effectivelly calculating the velocity at the given evaluation points, filling v.

This function result is a matrix with size d x m where 'm' is the number of evaluation points and 'd' the result dimension (2 or 3 for representing vxx, vyy and vzz)

IMPORTANT: This function currently requires that evalPoints are equal to the elements integration points, and in the same order.

Parameters

e	The element where we are calculating the result matrix
evalPoints	Matrix filled with natural coordinates for the 'm' points where results will be evaluated. Size = natural_coord_dimension x m.
v	Calculated results for each of the evaluation points. For each point of the 'm' points, a matrix column should be filled with the given point result. Size = d x m, where 'd' is the result dimension.

◆ fillElementData()

GmpFemPhysics::FemResultType GmpHydraulicDualFlow::fillElementData	(	const GmElement *	e,
		GmpFemMatrixSet &	elemMatrices,
		GmpFemVectorSet &	elemVectors )

virtual

Fills the element matrices and vectors for a generic dual flow element.

Reimplemented from GmpHydraulic.

Reimplemented in GmpHydraulicEnhancedDualFlow.

◆ fillFluidFlowMatrix()

void GmpHydraulicDualFlow::fillFluidFlowMatrix	(	const GmElement *	e,
		const GmVector *	coord,
		int	ip,
		GmMatrix &	Kfc,
		const GmMatrix &	Bp,
		double	c ) const

virtual

Calculates the flow stiffness matrix of fluid 2 due to permeability contribution or permeability matrix for element 'e' at the given integration point (coord, ip) storing the result at the Kfc matrix. Bp is the strain pore pressure matrix and c a constant equal to w * h * detJ.

This function can operate either in an iso parametric setting (Kc = n x n) or in a super parametric setting (Kfc = nf x nf) where n is the number of element nodes and nf the number of element corner nodes.

Returns Kfc = Bp.t() * Kaux * Bp * c, where c should be equal to w * h * detJ

◆ fillFractureStorageCompMatrix()

void GmpHydraulicDualFlow::fillFractureStorageCompMatrix	(	const GmElement *	e,
		const GmVector *	coord,
		int	ip,
		GmMatrix &	Sfc,
		const GmVector &	N,
		double	c ) const

virtual

Calculates the stiffness matrix of "storage / compressibility of fluid" for element 'e' at the given integration point (coord, ip) storing the result at the Sfc matrix. N is the shape function vector and c a constant equal to w * h * detJ.

This function can operate either in an iso parametric setting (Sfc = n x n) or in a super parametric setting (Sfc = nf x nf) where n is the number of element nodes and nf the number of element corner nodes.

Returns Sfc = N * (Bp + porosity/Kw) * N.t() * c, where c should be equal to w * h * detJ

◆ fixedNodalDofsBc()

bool GmpHydraulicDualFlow::fixedNodalDofsBc	(	QVector< int > &	nodes,
		QVector< int > &	dof,
		QVector< double > &	values,
		bool *	constantValues ) const

virtual

See comments on base class. Fills vectors with prescribed node pore pressures.

Reimplemented from GmpHydraulic.

◆ fixedNodalForcesBc()

bool GmpHydraulicDualFlow::fixedNodalForcesBc	(	QVector< int > &	nodes,
		QVector< int > &	dof,
		QVector< double > &	values ) const

virtual

See comments on base class. Fills vectors with prescribed nodal forces - flow.

Reimplemented from GmpHydraulic.

◆ physicsMetaDataMap()

const QVariantMap * GmpHydraulicDualFlow::physicsMetaDataMap ( )

virtual

Returns a reference for the Hydraulics dual flow attribute map, built when the function is called for the first time.

Reimplemented from GmpHydraulic.

Reimplemented in GmpHydraulicEnhancedDualFlow.

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

inc/gmpHydraulicDualFlow.h
src/gmpHydraulicDualFlow.cpp

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ DualFlowBoundaryConditionIds

◆ DualFlowBoundaryConditionValueIds

◆ DualFlowGaussAttributeIds

◆ DualFlowNodeAttributeIds

◆ DualFlowPhysicsAttributeIds

◆ DualFlowStateVarIds

Member Function Documentation

◆ calcDerivedResults()

◆ calcElementDerivedResultAtPoints()

◆ calcElementFlowAtPoints()

◆ calcElementVelocityAtPoints()

◆ fillElementData()

◆ fillFluidFlowMatrix()

◆ fillFractureStorageCompMatrix()

◆ fixedNodalDofsBc()

◆ fixedNodalForcesBc()

◆ physicsMetaDataMap()