A stabilized element for the incompressible Navier-Stokes equations, utilizing lagrangian_Eulerian approach for droplet dynamics. More...

#include <surface_tension.h>

Inheritance diagram for Kratos::SurfaceTension< TDim, TNumNodes >:

Collaboration diagram for Kratos::SurfaceTension< TDim, TNumNodes >:

Public Member Functions
Life Cycle
	SurfaceTension (IndexType NewId=0)
	Default constuctor. More...

	SurfaceTension (IndexType NewId, const NodesArrayType &ThisNodes)
	Constructor using an array of nodes. More...

	SurfaceTension (IndexType NewId, GeometryType::Pointer pGeometry)
	Constructor using a geometry object. More...

	SurfaceTension (IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties)
	Constuctor using geometry and properties. More...

	~SurfaceTension () override
	Destructor. More...

Operations
Element::Pointer	Create (IndexType NewId, NodesArrayType const &ThisNodes, PropertiesType::Pointer pProperties) const override
	Create a new element of this type. More...

Element::Pointer	Create (IndexType NewId, GeometryType::Pointer pGeom, PropertiesType::Pointer pProperties) const override
	It creates a new element pointer. More...

void	CalculateLocalSystem (MatrixType &rLeftHandSideMatrix, VectorType &rRightHandSideVector, ProcessInfo &rCurrentProcessInfo) override
	Provides local contributions from body forces and OSS projection terms. More...

void	CalculateLeftHandSide (MatrixType &rLeftHandSideMatrix, ProcessInfo &rCurrentProcessInfo) override
	Returns a zero matrix of appropiate size (provided for compatibility with scheme) More...

void	CalculateRightHandSide (VectorType &rRightHandSideVector, ProcessInfo &rCurrentProcessInfo) override
	Provides local contributions from body forces and projections to the RHS. More...

void	CalculateMassMatrix (MatrixType &rMassMatrix, ProcessInfo &rCurrentProcessInfo) override
	Computes local contributions to the mass matrix. More...

void	CalculateLocalVelocityContribution (MatrixType &rDampingMatrix, VectorType &rRightHandSideVector, ProcessInfo &rCurrentProcessInfo) override
	Computes the local contribution associated to 'new' velocity and pressure values. More...

void	FinalizeNonLinearIteration (ProcessInfo &rCurrentProcessInfo) override

void	Calculate (const Variable< double > &rVariable, double &rOutput, const ProcessInfo &rCurrentProcessInfo) override
	Implementation of Calculate to compute an error estimate. More...

void	Calculate (const Variable< array_1d< double, 3 > > &rVariable, array_1d< double, 3 > &rOutput, const ProcessInfo &rCurrentProcessInfo) override
	Implementation of Calculate to compute the local OSS projections. More...

void	EquationIdVector (EquationIdVectorType &rResult, ProcessInfo &rCurrentProcessInfo) override
	Provides the global indices for each one of this element's local rows. More...

void	GetDofList (DofsVectorType &rElementalDofList, ProcessInfo &rCurrentProcessInfo) override
	Returns a list of the element's Dofs. More...

void	GetFirstDerivativesVector (Vector &Values, int Step=0) override
	Returns VELOCITY_X, VELOCITY_Y, (VELOCITY_Z,) PRESSURE for each node. More...

void	GetSecondDerivativesVector (Vector &Values, int Step=0) override
	Returns ACCELERATION_X, ACCELERATION_Y, (ACCELERATION_Z,) 0 for each node. More...

void	GetValueOnIntegrationPoints (const Variable< array_1d< double, 3 > > &rVariable, std::vector< array_1d< double, 3 > > &rOutput, const ProcessInfo &rCurrentProcessInfo) override
	Obtain an array_1d<double,3> elemental variable, evaluated on gauss points. More...

void	GetValueOnIntegrationPoints (const Variable< double > &rVariable, std::vector< double > &rValues, const ProcessInfo &rCurrentProcessInfo) override
	Obtain a double elemental variable, evaluated on gauss points. More...

void	GetValueOnIntegrationPoints (const Variable< array_1d< double, 6 > > &rVariable, std::vector< array_1d< double, 6 > > &rValues, const ProcessInfo &rCurrentProcessInfo) override
	Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning. More...

void	GetValueOnIntegrationPoints (const Variable< Vector > &rVariable, std::vector< Vector > &rValues, const ProcessInfo &rCurrentProcessInfo) override
	Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning. More...

void	GetValueOnIntegrationPoints (const Variable< Matrix > &rVariable, std::vector< Matrix > &rValues, const ProcessInfo &rCurrentProcessInfo) override
	Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning. More...

Elemental Data
int	Check (const ProcessInfo &rCurrentProcessInfo) override
	Checks the input and that all required Kratos variables have been registered. More...

Input and output
std::string	Info () const override
	Turn back information as a string. More...

void	PrintInfo (std::ostream &rOStream) const override
	Print information about this object. More...

Public Member Functions inherited from Kratos::Element
	Element (IndexType NewId=0)

	Element (IndexType NewId, const NodesArrayType &ThisNodes)

	Element (IndexType NewId, GeometryType::Pointer pGeometry)

	Element (IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties)

	Element (Element const &rOther)
	Copy constructor. More...

	~Element () override
	Destructor. More...

Element &	operator= (Element const &rOther)
	Assignment operator. More...

virtual Pointer	Clone (IndexType NewId, NodesArrayType const &ThisNodes) const
	It creates a new element pointer and clones the previous element data. More...

virtual void	EquationIdVector (EquationIdVectorType &rResult, const ProcessInfo &rCurrentProcessInfo) const

virtual void	GetDofList (DofsVectorType &rElementalDofList, const ProcessInfo &rCurrentProcessInfo) const

virtual IntegrationMethod	GetIntegrationMethod () const

virtual void	GetValuesVector (Vector &values, int Step=0) const

virtual void	GetFirstDerivativesVector (Vector &values, int Step=0) const

virtual void	GetSecondDerivativesVector (Vector &values, int Step=0) const

virtual void	Initialize (const ProcessInfo &rCurrentProcessInfo)

virtual void	ResetConstitutiveLaw ()

virtual void	InitializeSolutionStep (const ProcessInfo &rCurrentProcessInfo)

virtual void	InitializeNonLinearIteration (const ProcessInfo &rCurrentProcessInfo)

virtual void	FinalizeNonLinearIteration (const ProcessInfo &rCurrentProcessInfo)

virtual void	FinalizeSolutionStep (const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateLocalSystem (MatrixType &rLeftHandSideMatrix, VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateLeftHandSide (MatrixType &rLeftHandSideMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateRightHandSide (VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateFirstDerivativesContributions (MatrixType &rLeftHandSideMatrix, VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateFirstDerivativesLHS (MatrixType &rLeftHandSideMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateFirstDerivativesRHS (VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateSecondDerivativesContributions (MatrixType &rLeftHandSideMatrix, VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateSecondDerivativesLHS (MatrixType &rLeftHandSideMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateSecondDerivativesRHS (VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateMassMatrix (MatrixType &rMassMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateDampingMatrix (MatrixType &rDampingMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateLumpedMassVector (VectorType &rLumpedMassVector, const ProcessInfo &rCurrentProcessInfo) const

virtual void	AddExplicitContribution (const ProcessInfo &rCurrentProcessInfo)

virtual void	AddExplicitContribution (const VectorType &rRHSVector, const Variable< VectorType > &rRHSVariable, const Variable< double > &rDestinationVariable, const ProcessInfo &rCurrentProcessInfo)
	This function is designed to make the element to assemble an rRHS vector identified by a variable rRHSVariable by assembling it to the nodes on the variable rDestinationVariable. (This is the double version) More...

virtual void	AddExplicitContribution (const VectorType &rRHSVector, const Variable< VectorType > &rRHSVariable, const Variable< array_1d< double, 3 > > &rDestinationVariable, const ProcessInfo &rCurrentProcessInfo)
	This function is designed to make the element to assemble an rRHS vector identified by a variable rRHSVariable by assembling it to the nodes on the variable rDestinationVariable. (This is the vector version) More...

virtual void	AddExplicitContribution (const MatrixType &rLHSMatrix, const Variable< MatrixType > &rLHSVariable, const Variable< Matrix > &rDestinationVariable, const ProcessInfo &rCurrentProcessInfo)
	This function is designed to make the element to assemble an rRHS vector identified by a variable rRHSVariable by assembling it to the nodes on the variable rDestinationVariable. (This is the matrix version) More...

virtual void	Calculate (const Variable< Vector > &rVariable, Vector &Output, const ProcessInfo &rCurrentProcessInfo)

virtual void	Calculate (const Variable< Matrix > &rVariable, Matrix &Output, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< bool > &rVariable, std::vector< bool > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< int > &rVariable, std::vector< int > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< double > &rVariable, std::vector< double > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< array_1d< double, 3 >> &rVariable, std::vector< array_1d< double, 3 >> &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< array_1d< double, 4 >> &rVariable, std::vector< array_1d< double, 4 >> &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< array_1d< double, 6 >> &rVariable, std::vector< array_1d< double, 6 >> &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< array_1d< double, 9 >> &rVariable, std::vector< array_1d< double, 9 >> &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< Vector > &rVariable, std::vector< Vector > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< Matrix > &rVariable, std::vector< Matrix > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateOnIntegrationPoints (const Variable< ConstitutiveLaw::Pointer > &rVariable, std::vector< ConstitutiveLaw::Pointer > &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< bool > &rVariable, const std::vector< bool > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< int > &rVariable, const std::vector< int > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< double > &rVariable, const std::vector< double > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< array_1d< double, 3 >> &rVariable, const std::vector< array_1d< double, 3 >> &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< array_1d< double, 4 >> &rVariable, const std::vector< array_1d< double, 4 >> &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< array_1d< double, 6 >> &rVariable, const std::vector< array_1d< double, 6 >> &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< array_1d< double, 9 >> &rVariable, const std::vector< array_1d< double, 9 >> &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< Vector > &rVariable, const std::vector< Vector > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< Matrix > &rVariable, const std::vector< Matrix > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual void	SetValuesOnIntegrationPoints (const Variable< ConstitutiveLaw::Pointer > &rVariable, const std::vector< ConstitutiveLaw::Pointer > &rValues, const ProcessInfo &rCurrentProcessInfo)

virtual int	Check (const ProcessInfo &rCurrentProcessInfo) const

virtual void	MassMatrix (MatrixType &rMassMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	AddMassMatrix (MatrixType &rLeftHandSideMatrix, double coeff, const ProcessInfo &rCurrentProcessInfo)

virtual void	DampMatrix (MatrixType &rDampMatrix, const ProcessInfo &rCurrentProcessInfo)

virtual void	AddInertiaForces (VectorType &rRightHandSideVector, double coeff, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateLocalVelocityContribution (MatrixType &rDampingMatrix, VectorType &rRightHandSideVector, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateSensitivityMatrix (const Variable< double > &rDesignVariable, Matrix &rOutput, const ProcessInfo &rCurrentProcessInfo)

virtual void	CalculateSensitivityMatrix (const Variable< array_1d< double, 3 > > &rDesignVariable, Matrix &rOutput, const ProcessInfo &rCurrentProcessInfo)

PropertiesType::Pointer	pGetProperties ()
	returns the pointer to the property of the element. Does not throw an error, to allow copying of elements which don't have any property assigned. More...

const PropertiesType::Pointer	pGetProperties () const

PropertiesType &	GetProperties ()

PropertiesType const &	GetProperties () const

void	SetProperties (PropertiesType::Pointer pProperties)

bool	HasProperties () const
	Check that the Element has a correctly initialized pointer to a Properties instance. More...

virtual const Parameters	GetSpecifications () const
	This method provides the specifications/requirements of the element. More...

void	PrintData (std::ostream &rOStream) const override
	Print object's data. More...

	KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION (Element)

Public Member Functions inherited from Kratos::GeometricalObject
	GeometricalObject (IndexType NewId=0)
	Default constructor. More...

	GeometricalObject (IndexType NewId, GeometryType::Pointer pGeometry)
	Default constructor. More...

	~GeometricalObject () override
	Destructor. More...

	GeometricalObject (GeometricalObject const &rOther)
	Copy constructor. More...

GeometricalObject &	operator= (GeometricalObject const &rOther)
	Assignment operator. More...

virtual void	SetGeometry (GeometryType::Pointer pGeometry)
	Sets the pointer to the geometry. More...

GeometryType::Pointer	pGetGeometry ()
	Returns the pointer to the geometry. More...

const GeometryType::Pointer	pGetGeometry () const
	Returns the pointer to the geometry (const version) More...

GeometryType &	GetGeometry ()
	Returns the reference of the geometry. More...

GeometryType const &	GetGeometry () const
	Returns the reference of the geometry (const version) More...

Flags &	GetFlags ()
	Returns the flags of the object. More...

Flags const &	GetFlags () const
	Returns the flags of the object (const version) More...

void	SetFlags (Flags const &rThisFlags)
	Sets the flags of the object. More...

DataValueContainer &	Data ()

DataValueContainer &	GetData ()

DataValueContainer const &	GetData () const

void	SetData (DataValueContainer const &rThisData)

template<class TDataType >
bool	Has (const Variable< TDataType > &rThisVariable) const

template<class TVariableType >
void	SetValue (const TVariableType &rThisVariable, typename TVariableType::Type const &rValue)

template<class TVariableType >
TVariableType::Type &	GetValue (const TVariableType &rThisVariable)

template<class TVariableType >
TVariableType::Type const &	GetValue (const TVariableType &rThisVariable) const

void	PrintData (std::ostream &rOStream) const override
	Print object's data. More...

unsigned int	use_count () const noexcept

	KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION (GeometricalObject)
	Pointer definition of GeometricalObject. More...

bool	IsActive () const
	Checks if the GeometricalObject is active. More...

Public Member Functions inherited from Kratos::IndexedObject
	IndexedObject (IndexType NewId=0)
	Default constructor. More...

virtual	~IndexedObject ()
	Destructor. More...

	IndexedObject (IndexedObject const &rOther)
	Copy constructor. More...

IndexedObject &	operator= (IndexedObject const &rOther)
	Assignment operator. More...

template<class TObjectType >
IndexType	operator() (TObjectType const &rThisObject) const

IndexType	Id () const

IndexType	GetId () const

virtual void	SetId (IndexType NewId)

IndexType &	DepricatedIdAccess ()
	TODO: remove this function when removing data_file_io object. More...

	KRATOS_CLASS_POINTER_DEFINITION (IndexedObject)
	Pointer definition of IndexedObject. More...

Public Member Functions inherited from Kratos::Flags
Flags &	operator= (Flags const &rOther)
	Assignment operator. More...

	operator bool () const

Flags	operator~ () const

bool	operator! () const

void	AssignFlags (Flags const &rOther)

void	Set (const Flags ThisFlag)

void	Set (const Flags ThisFlag, bool Value)

void	Reset (const Flags ThisFlag)

void	Flip (const Flags ThisFlag)

void	SetPosition (IndexType Position, bool Value=true)

bool	GetPosition (IndexType Position) const

void	FlipPosition (IndexType Position)

void	ClearPosition (IndexType Position)

void	Clear ()

Flags	AsFalse () const

bool	Is (Flags const &rOther) const

bool	IsDefined (Flags const &rOther) const

bool	IsNot (Flags const &rOther) const

bool	IsNotDefined (Flags const &rOther) const

	KRATOS_CLASS_POINTER_DEFINITION (Flags)
	Pointer definition of Flags. More...

const Flags &	operator\|= (const Flags &Other)

const Flags &	operator&= (const Flags &Other)

	Flags ()
	Default constructor. More...

	Flags (Flags const &rOther)
	Copy constructor. More...

virtual	~Flags ()
	Destructor. More...

Protected Member Functions
Protected Operations
virtual void	CalculateTau (double &TauOne, double &TauTwo, const array_1d< double, 3 > &rAdvVel, const double ElemSize, const double Density, const double Viscosity, const ProcessInfo &rCurrentProcessInfo)
	Calculate Stabilization parameters. More...

virtual void	CalculateStaticTau (double &TauOne, const array_1d< double, 3 > &rAdvVel, const double ElemSize, const double Density, const double Viscosity)
	Calculate momentum stabilization parameter (without time term). More...

virtual void	AddMomentumRHS (VectorType &F, const double Density, const array_1d< double, TNumNodes > &rShapeFunc, const double Weight)
	Add the momentum equation contribution to the RHS (body forces) More...

virtual void	AddProjectionToRHS (VectorType &RHS, const array_1d< double, 3 > &rAdvVel, const double Density, const double TauOne, const double TauTwo, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight, const double DeltaTime=1.0)
	Add OSS projection terms to the RHS. More...

void	CalculateLumpedMassMatrix (MatrixType &rLHSMatrix, const double Mass)
	Add lumped mass matrix. More...

void	AddConsistentMassMatrixContribution (MatrixType &rLHSMatrix, const array_1d< double, TNumNodes > &rShapeFunc, const double Density, const double Weight)

void	AddMassStabTerms (MatrixType &rLHSMatrix, const double Density, const array_1d< double, 3 > &rAdvVel, const double TauOne, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Add mass-like stabilization terms to LHS. More...

void	AddIntegrationPointVelocityContribution (MatrixType &rDampingMatrix, VectorType &rDampRHS, const double Density, const double Viscosity, const array_1d< double, 3 > &rAdvVel, const double TauOne, const double TauTwo, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Add a the contribution from a single integration point to the velocity contribution. More...

void	ApplySurfaceTensionContribution (MatrixType &rDampingMatrix, VectorType &rRightHandSideVector, const array_1d< double, 3 > &node_indx, const int &k, const ProcessInfo &rCurrentProcessInfo)
	Add the surface tension term to the velocity contribution. More...

void	ApplySurfaceTensionContribution3D (MatrixType &rDampingMatrix, VectorType &rRightHandSideVector, const array_1d< double, 4 > &node_indx, const int &k, const ProcessInfo &rCurrentProcessInfo)

void	AddViscousStress2D ()
	Add the viscous stress to air domain in two dimensions. More...

void	AddProjectionResidualContribution (const array_1d< double, 3 > &rAdvVel, const double Density, array_1d< double, 3 > &rElementalMomRes, double &rElementalMassRes, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Assemble the contribution from an integration point to the element's residual. More...

void	ASGSMomResidual (const array_1d< double, 3 > &rAdvVel, const double Density, array_1d< double, 3 > &rElementalMomRes, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Assemble the contribution from an integration point to the element's residual. More...

void	OSSMomResidual (const array_1d< double, 3 > &rAdvVel, const double Density, array_1d< double, 3 > &rElementalMomRes, const array_1d< double, TNumNodes > &rShapeFunc, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Assemble the contribution from an integration point to the element's residual. More...

virtual double	EffectiveViscosity (double Density, const array_1d< double, TNumNodes > &rN, const BoundedMatrix< double, TNumNodes, TDim > &rDN_DX, double ElemSize, const ProcessInfo &rProcessInfo)
	EffectiveViscosity Calculate the viscosity at given integration point, using Smagorinsky if enabled. More...

double	EquivalentStrainRate (const BoundedMatrix< double, TNumNodes, TDim > &rDN_DX) const
	EquivalentStrainRate Calculate the second invariant of the strain rate tensor GammaDot = (2SijSij)^0.5. More...

void	GetAdvectiveVel (array_1d< double, 3 > &rAdvVel, const array_1d< double, TNumNodes > &rShapeFunc)
	Write the advective velocity evaluated at this point to an array. More...

virtual void	GetAdvectiveVel (array_1d< double, 3 > &rAdvVel, const array_1d< double, TNumNodes > &rShapeFunc, const std::size_t Step)
	Write the advective velocity evaluated at this point to an array. More...

void	GetConvectionOperator (array_1d< double, TNumNodes > &rResult, const array_1d< double, 3 > &rVelocity, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv)
	Write the convective operator evaluated at this point (for each nodal funciton) to an array. More...

virtual void	EvaluateInPoint (double &rResult, const Variable< double > &rVariable, const array_1d< double, TNumNodes > &rShapeFunc)
	Write the value of a variable at a point inside the element to a double. More...

virtual void	EvaluateInPoint (array_1d< double, 3 > &rResult, const Variable< array_1d< double, 3 > > &rVariable, const array_1d< double, TNumNodes > &rShapeFunc)
	Write the value of a variable at a point inside the element to a double. More...

double	ElementSize (const double)
	Return an estimate for the element size h, used to calculate the stabilization parameters. More...

virtual void	AddViscousTerm (MatrixType &rDampingMatrix, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Adds the contribution of the viscous term to the momentum equation. More...

void	AddBTransCB (MatrixType &rDampingMatrix, BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv, const double Weight)
	Adds the contribution of the viscous term to the momentum equation (alternate). More...

void	ModulatedGradientDiffusion (MatrixType &rDampingMatrix, const BoundedMatrix< double, TNumNodes, TDim > &rDN_DX, const double Weight)

void	CalculateB (BoundedMatrix< double,(TDim TNumNodes)/2, TDim TNumNodes > &rB, const BoundedMatrix< double, TNumNodes, TDim > &rShapeDeriv)
	Calculate the strain rate matrix. More...

virtual void	CalculateC (BoundedMatrix< double,(TDim TNumNodes)/2,(TDim TNumNodes)/2 > &rC, const double Viscosity)
	Calculate a matrix that provides the stress given the strain rate. More...

double	ConsistentMassCoef (const double Area)

double	SubscaleErrorEstimate (const ProcessInfo &rProcessInfo)

array_1d< double, 2 >	Vector2D (const double x0, const double y0, const double x1, const double y1)

array_1d< double, 3 >	Vector3D (const double x0, const double y0, const double z0, const double x1, const double y1, const double z1)

array_1d< double, 3 >	CrossProduct3D (const array_1d< double, 3 > &a, const array_1d< double, 3 > &b)

double	Norm2D (const array_1d< double, 2 > &a)

double	Norm3D (const array_1d< double, 3 > &a)

void	NormalizeVec2D (array_1d< double, 2 > &input)

void	NormalizeVec3D (array_1d< double, 3 > &input)

double	DotProduct2D (const array_1d< double, 2 > &a, const array_1d< double, 2 > &b)

double	DotProduct3D (const array_1d< double, 3 > &a, const array_1d< double, 3 > &b)

Type Definitions
typedef IndexedObject	BaseType
	base type: an IndexedObject that automatically has a unique number More...

typedef Node	NodeType
	definition of node type (default is: Node) More...

typedef Properties	PropertiesType

typedef Geometry< NodeType >	GeometryType
	definition of the geometry type with given NodeType More...

typedef Geometry< NodeType >::PointsArrayType	NodesArrayType
	definition of nodes container type, redefined from GeometryType More...

typedef Vector	VectorType

typedef Matrix	MatrixType

typedef std::size_t	IndexType

typedef std::size_t	SizeType

typedef std::vector< std::size_t >	EquationIdVectorType

typedef std::vector< Dof< double >::Pointer >	DofsVectorType

typedef PointerVectorSet< Dof< double >, IndexedObject >	DofsArrayType

typedef array_1d< double, TNumNodes >	ShapeFunctionsType

typedef BoundedMatrix< double, TNumNodes, TDim >	ShapeFunctionDerivativesType

	KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION (SurfaceTension)
	Pointer definition of SurfaceTension. More...

Serialization
class	Serializer

Specialized implementation of SurfaceTension for functions that depend on TDim
void	EquationIdVector (EquationIdVectorType &rResult, ProcessInfo &rCurrentProcessInfo)

void	EquationIdVector (EquationIdVectorType &rResult, ProcessInfo &rCurrentProcessInfo)

void	GetDofList (DofsVectorType &rElementalDofList, ProcessInfo &rCurrentProcessInfo)

void	GetDofList (DofsVectorType &rElementalDofList, ProcessInfo &rCurrentProcessInfo)

void	GetFirstDerivativesVector (Vector &Values, int Step)

void	GetFirstDerivativesVector (Vector &Values, int Step)

void	GetSecondDerivativesVector (Vector &Values, int Step)

void	GetSecondDerivativesVector (Vector &Values, int Step)

void	GetValueOnIntegrationPoints (const Variable< array_1d< double, 3 > > &rVariable, std::vector< array_1d< double, 3 > > &rOutput, const ProcessInfo &rCurrentProcessInfo)

void	GetValueOnIntegrationPoints (const Variable< array_1d< double, 3 > > &rVariable, std::vector< array_1d< double, 3 > > &rOutput, const ProcessInfo &rCurrentProcessInfo)

double	ElementSize (const double Area)

double	ElementSize (const double Volume)

double	EquivalentStrainRate (const boost::numeric::ublas::bounded_matrix< double, 3, 2 > &rDN_DX) const

double	EquivalentStrainRate (const boost::numeric::ublas::bounded_matrix< double, 4, 3 > &rDN_DX) const

void	CalculateB (boost::numeric::ublas::bounded_matrix< double, 3, 6 > &rB, const boost::numeric::ublas::bounded_matrix< double, 3, 2 > &rShapeDeriv)

void	CalculateB (boost::numeric::ublas::bounded_matrix< double, 6, 12 > &rB, const boost::numeric::ublas::bounded_matrix< double, 4, 3 > &rShapeDeriv)

void	CalculateC (boost::numeric::ublas::bounded_matrix< double, 3, 3 > &rC, const double Viscosity)

void	CalculateC (boost::numeric::ublas::bounded_matrix< double, 6, 6 > &rC, const double Viscosity)

void	AddViscousTerm (MatrixType &rDampingMatrix, const boost::numeric::ublas::bounded_matrix< double, 3, 2 > &rShapeDeriv, const double Weight)

void	AddViscousTerm (MatrixType &rDampingMatrix, const boost::numeric::ublas::bounded_matrix< double, 4, 3 > &rShapeDeriv, const double Weight)

double	ConsistentMassCoef (const double Area)

double	ConsistentMassCoef (const double Volume)

Additional Inherited Members
Public Types inherited from Kratos::Element
typedef Element	ElementType
	definition of element type More...

typedef GeometricalObject	BaseType
	base type: an GeometricalObject that automatically has a unique number More...

typedef Node	NodeType
	definition of node type (default is: Node) More...

typedef Properties	PropertiesType

typedef Geometry< NodeType >	GeometryType
	definition of the geometry type with given NodeType More...

typedef Geometry< NodeType >::PointsArrayType	NodesArrayType
	definition of nodes container type, redefined from GeometryType More...

typedef Vector	VectorType

typedef Matrix	MatrixType

typedef std::size_t	IndexType

typedef std::size_t	SizeType

typedef Dof< double >	DofType

typedef std::vector< std::size_t >	EquationIdVectorType

typedef std::vector< DofType::Pointer >	DofsVectorType

typedef PointerVectorSet< DofType >	DofsArrayType

typedef GeometryData::IntegrationMethod	IntegrationMethod
	Type definition for integration methods. More...

typedef GeometryData	GeometryDataType

Public Types inherited from Kratos::GeometricalObject
typedef Node	NodeType
	Definition of the node type. More...

typedef Geometry< NodeType >	GeometryType
	The geometry type definition. More...

typedef std::size_t	IndexType
	Defines the index type. More...

typedef std::size_t	result_type
	Defines the result type. More...

Public Types inherited from Kratos::IndexedObject
typedef std::size_t	IndexType
	The definition of the index type. More...

typedef std::size_t	result_type
	The definition of the result_type. More...

Public Types inherited from Kratos::Flags
enum	FlagsList { Flag0 = BlockType(1) , Flag1 = BlockType(1) << 1 , Flag2 = BlockType(1) << 2 , Flag3 = BlockType(1) << 3 , Flag4 = BlockType(1) << 4 , Flag5 = BlockType(1) << 5 , Flag6 = BlockType(1) << 6 , Flag7 = BlockType(1) << 7 , Flag8 = BlockType(1) << 8 , Flag9 = BlockType(1) << 9 , Flag10 = BlockType(1) << 10 , Flag11 = BlockType(1) << 11 , Flag12 = BlockType(1) << 12 , Flag13 = BlockType(1) << 13 , Flag14 = BlockType(1) << 14 , Flag15 = BlockType(1) << 15 , Flag16 = BlockType(1) << 16 , Flag17 = BlockType(1) << 17 , Flag18 = BlockType(1) << 18 , Flag19 = BlockType(1) << 19 , Flag20 = BlockType(1) << 20 , Flag21 = BlockType(1) << 21 , Flag22 = BlockType(1) << 22 , Flag23 = BlockType(1) << 23 , Flag24 = BlockType(1) << 24 , Flag25 = BlockType(1) << 25 , Flag26 = BlockType(1) << 26 , Flag27 = BlockType(1) << 27 , Flag28 = BlockType(1) << 28 , Flag29 = BlockType(1) << 29 , Flag30 = BlockType(1) << 30 }

typedef int64_t	BlockType

typedef int64_t	FlagType

typedef std::size_t	IndexType

Static Public Member Functions inherited from Kratos::GeometricalObject
static bool	HasSameType (const GeometricalObject &rLHS, const GeometricalObject &rRHS)
	Checks if two GeometricalObject have the same type. More...

static bool	HasSameType (const GeometricalObject rLHS, const GeometricalObject rRHS)
	Checks if two GeometricalObject have the same type (pointer version) More...

static bool	HasSameGeometryType (const GeometricalObject &rLHS, const GeometricalObject &rRHS)
	Checks if two GeometricalObject have the same geometry type. More...

static bool	HasSameGeometryType (const GeometricalObject rLHS, const GeometricalObject rRHS)
	Checks if two GeometricalObject have the same geometry type (pointer version) More...

static bool	IsSame (const GeometricalObject &rLHS, const GeometricalObject &rRHS)
	Checks if two GeometricalObject are the same. More...

static bool	IsSame (const GeometricalObject rLHS, const GeometricalObject rRHS)
	Checks if two GeometricalObject are the same (pointer version) More...

Static Public Member Functions inherited from Kratos::Flags
static const Flags	AllDefined ()

static const Flags	AllTrue ()

static Flags	Create (IndexType ThisPosition, bool Value=true)

Detailed Description

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>
class Kratos::SurfaceTension< TDim, TNumNodes >

A stabilized element for the incompressible Navier-Stokes equations, utilizing lagrangian_Eulerian approach for droplet dynamics.

This class is based on the VMS element in the fluid dynamics, It is extended to solve for the droplet dynamics, and surface tensions. for more information about the the governer equations and thier discretization, please find it in the below references: Alex Jarauta, Pavel Ryzhakov, Marc Secanell, Prashant R Waghmare, and Jordi Pons-Prats. Numerical study of droplet dynamics in a polymer electrolyte fuel cell gas channel using an embedded eulerian-lagrangian approach. Journal of Power Sources, 323:201–212, 2016.

Member Typedef Documentation

◆ BaseType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef IndexedObject Kratos::SurfaceTension< TDim, TNumNodes >::BaseType

base type: an IndexedObject that automatically has a unique number

◆ DofsArrayType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef PointerVectorSet<Dof<double>, IndexedObject> Kratos::SurfaceTension< TDim, TNumNodes >::DofsArrayType

◆ DofsVectorType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef std::vector< Dof<double>::Pointer > Kratos::SurfaceTension< TDim, TNumNodes >::DofsVectorType

◆ EquationIdVectorType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef std::vector<std::size_t> Kratos::SurfaceTension< TDim, TNumNodes >::EquationIdVectorType

◆ GeometryType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Geometry<NodeType> Kratos::SurfaceTension< TDim, TNumNodes >::GeometryType

definition of the geometry type with given NodeType

◆ IndexType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef std::size_t Kratos::SurfaceTension< TDim, TNumNodes >::IndexType

◆ MatrixType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Matrix Kratos::SurfaceTension< TDim, TNumNodes >::MatrixType

◆ NodesArrayType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Geometry<NodeType>::PointsArrayType Kratos::SurfaceTension< TDim, TNumNodes >::NodesArrayType

definition of nodes container type, redefined from GeometryType

◆ NodeType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Node Kratos::SurfaceTension< TDim, TNumNodes >::NodeType

definition of node type (default is: Node)

◆ PropertiesType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Properties Kratos::SurfaceTension< TDim, TNumNodes >::PropertiesType

Properties are used to store any parameters related to the constitutive law

◆ ShapeFunctionDerivativesType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef BoundedMatrix<double, TNumNodes, TDim> Kratos::SurfaceTension< TDim, TNumNodes >::ShapeFunctionDerivativesType

◆ ShapeFunctionsType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef array_1d<double, TNumNodes> Kratos::SurfaceTension< TDim, TNumNodes >::ShapeFunctionsType

◆ SizeType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef std::size_t Kratos::SurfaceTension< TDim, TNumNodes >::SizeType

◆ VectorType

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

typedef Vector Kratos::SurfaceTension< TDim, TNumNodes >::VectorType

Constructor & Destructor Documentation

◆ SurfaceTension() [1/4]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::SurfaceTension ( IndexType NewId = 0 )

inline

Default constuctor.

Parameters

NewId Index number of the new element (optional)

◆ SurfaceTension() [2/4]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::SurfaceTension	(	IndexType	NewId,
		const NodesArrayType &	ThisNodes
	)

inline

Constructor using an array of nodes.

Parameters

NewId	Index of the new element
ThisNodes	An array containing the nodes of the new element

◆ SurfaceTension() [3/4]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::SurfaceTension	(	IndexType	NewId,
		GeometryType::Pointer	pGeometry
	)

inline

Constructor using a geometry object.

Parameters

NewId	Index of the new element
pGeometry	Pointer to a geometry object

◆ SurfaceTension() [4/4]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::SurfaceTension	(	IndexType	NewId,
		GeometryType::Pointer	pGeometry,
		PropertiesType::Pointer	pProperties
	)

inline

Constuctor using geometry and properties.

Parameters

NewId	Index of the new element
pGeometry	Pointer to a geometry object
pProperties	Pointer to the element's properties

◆ ~SurfaceTension()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::~SurfaceTension ( )

inlineoverride

Destructor.

Member Function Documentation

◆ AddBTransCB()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddBTransCB	(	MatrixType &	rDampingMatrix,
		BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Adds the contribution of the viscous term to the momentum equation (alternate).

This function is an alternate implementation of ST::AddViscousTerm. This version works with ublas matrices, using the relationship between stress and rate of strain given by ST::CalculateC. It is currently unused (as ST::AddViscousTerm is a more efficient implementation of the Cauchy equation) but it is left here so derived classes can use it to implement other constitutive equations.

Parameters

rDampingMatrix	Elemental Damping matrix
rShapeDeriv	Elemental shape function derivatives
Weight	Effective viscosity, in dynamic units, weighted by the integration point area

◆ AddConsistentMassMatrixContribution()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddConsistentMassMatrixContribution	(	MatrixType &	rLHSMatrix,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const double	Density,
		const double	Weight
	)

inlineprotected

◆ AddIntegrationPointVelocityContribution()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddIntegrationPointVelocityContribution	(	MatrixType &	rDampingMatrix,
		VectorType &	rDampRHS,
		const double	Density,
		const double	Viscosity,
		const array_1d< double, 3 > &	rAdvVel,
		const double	TauOne,
		const double	TauTwo,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Add a the contribution from a single integration point to the velocity contribution.

◆ AddMassStabTerms()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddMassStabTerms	(	MatrixType &	rLHSMatrix,
		const double	Density,
		const array_1d< double, 3 > &	rAdvVel,
		const double	TauOne,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Add mass-like stabilization terms to LHS.

This function is only used in ASGS. For OSS, we avoid computing these terms, as they shoud cancel out with the dynamic part of the projection (which is not computed either)

Parameters

rLHSMatrix	Left hand side of the velocity-pressure system
Density	Density on integration point
rAdvVel	Advective velocity on integration point
TauOne	Stabilization parameter for momentum equation
rShapeFunc	Shape funcitions evaluated on integration point
rShapeDeriv	Shape function derivatives evaluated on integration point
Weight	Area (or volume) times integration point weight

◆ AddMomentumRHS()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::AddMomentumRHS	(	VectorType &	F,
		const double	Density,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const double	Weight
	)

inlineprotectedvirtual

Add the momentum equation contribution to the RHS (body forces)

◆ AddProjectionResidualContribution()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddProjectionResidualContribution	(	const array_1d< double, 3 > &	rAdvVel,
		const double	Density,
		array_1d< double, 3 > &	rElementalMomRes,
		double &	rElementalMassRes,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Assemble the contribution from an integration point to the element's residual.

Note that the dynamic term is not included in the momentum equation. If OSS_SWITCH = 1, we don't take into account the 'dynamic' stabilization terms, as it they belong to the finite element space.

◆ AddProjectionToRHS()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::AddProjectionToRHS	(	VectorType &	RHS,
		const array_1d< double, 3 > &	rAdvVel,
		const double	Density,
		const double	TauOne,
		const double	TauTwo,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight,
		const double	DeltaTime = `1.0`
	)

inlineprotectedvirtual

Add OSS projection terms to the RHS.

◆ AddViscousStress2D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::AddViscousStress2D ( )

inlineprotected

Add the viscous stress to air domain in two dimensions.

◆ AddViscousTerm() [1/3]

void Kratos::SurfaceTension< 2, 3 >::AddViscousTerm	(	MatrixType &	rDampingMatrix,
		const boost::numeric::ublas::bounded_matrix< double, 3, 2 > &	rShapeDeriv,
		const double	Weight
	)

protected

See also: SurfaceTension::AddViscousTerm

◆ AddViscousTerm() [2/3]

void Kratos::SurfaceTension< 3, 4 >::AddViscousTerm	(	MatrixType &	rDampingMatrix,
		const boost::numeric::ublas::bounded_matrix< double, 4, 3 > &	rShapeDeriv,
		const double	Weight
	)

protected

See also: SurfaceTension::AddViscousTerm

◆ AddViscousTerm() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::AddViscousTerm	(	MatrixType &	rDampingMatrix,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

protectedvirtual

Adds the contribution of the viscous term to the momentum equation.

The viscous term is written in stress-divergence (Cauchy) form.

Parameters

rDampingMatrix	Elemental Damping matrix
rShapeDeriv	Elemental shape function derivatives
Weight	Effective viscosity, in dynamic units, weighted by the integration point area

◆ ApplySurfaceTensionContribution()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::ApplySurfaceTensionContribution	(	MatrixType &	rDampingMatrix,
		VectorType &	rRightHandSideVector,
		const array_1d< double, 3 > &	node_indx,
		const int &	k,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineprotected

Add the surface tension term to the velocity contribution.

◆ ApplySurfaceTensionContribution3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::ApplySurfaceTensionContribution3D	(	MatrixType &	rDampingMatrix,
		VectorType &	rRightHandSideVector,
		const array_1d< double, 4 > &	node_indx,
		const int &	k,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineprotected

◆ ASGSMomResidual()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::ASGSMomResidual	(	const array_1d< double, 3 > &	rAdvVel,
		const double	Density,
		array_1d< double, 3 > &	rElementalMomRes,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Assemble the contribution from an integration point to the element's residual.

ASGS version. Note that rElementalMomRes should be initialized before calling this.

Parameters

rAdvVel	Convection velocity (not including subscale)
Density	Fluid density evaluated at integration point
rElementalMomRes	Result
rShapeFunc	Shape functions evaluated at integration point
rShapeDeriv	Shape function derivatives evaluated at integration point
Weight	Integration point weight (as a fraction of area or volume)

◆ Calculate() [1/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::Calculate	(	const Variable< array_1d< double, 3 > > &	rVariable,
		array_1d< double, 3 > &	rOutput,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverridevirtual

Implementation of Calculate to compute the local OSS projections.

If rVariable == ADVPROJ, This function computes the OSS projection terms using pressure and velocity values from the previous iteration. The projections are then added to the nodal variables ADVPROJ (Momentum residual) and DIVPROJ (Mass continuity residual). It is assumed that the scheme will divide the result by the assembled NODAL_AREA, which is equivalent to a nodal interpolation using a lumped mass matrix.

Parameters

rVariable	Use ADVPROJ
Output	Will be overwritten with the elemental momentum error
rCurrentProcessInfo	Process info instance (unused)

Return output

Reimplemented from Kratos::Element.

◆ Calculate() [2/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::Calculate	(	const Variable< double > &	rVariable,
		double &	rOutput,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverridevirtual

Implementation of Calculate to compute an error estimate.

If rVariable == ERROR_RATIO, this function will provide an a posteriori estimate of the norm of the subscale velocity, calculated as TauOne*||MomentumResidual||. Note that the residual of the momentum equation is evaluated at the element center and that the result has units of velocity (L/T). The error estimate both saved as the elemental ERROR_RATIO variable and returned as rOutput. If rVARIABLE == NODAL_AREA, the element's contribution to nodal area is added to its nodes.

Parameters

rVariable	Use ERROR_RATIO or NODAL_AREA
rOutput	Returns the error estimate for ERROR_RATIO, unused for NODAL_AREA
rCurrentProcessInfo	Process info instance (will be checked for OSS_SWITCH)

See also: MarkForRefinement for a use of the error ratio

Reimplemented from Kratos::Element.

◆ CalculateB() [1/3]

void Kratos::SurfaceTension< 2, 3 >::CalculateB	(	boost::numeric::ublas::bounded_matrix< double, 3, 6 > &	rB,
		const boost::numeric::ublas::bounded_matrix< double, 3, 2 > &	rShapeDeriv
	)

protected

See SurfaceTension::CalculateB

◆ CalculateB() [2/3]

void Kratos::SurfaceTension< 3, 4 >::CalculateB	(	boost::numeric::ublas::bounded_matrix< double, 6, 12 > &	rB,
		const boost::numeric::ublas::bounded_matrix< double, 4, 3 > &	rShapeDeriv
	)

protected

See SurfaceTension::CalculateB

◆ CalculateB() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateB	(	BoundedMatrix< double,(TDim TNumNodes)/2, TDim TNumNodes > &	rB,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv
	)

protected

Calculate the strain rate matrix.

Unused, left to support derived classes.

See also: SurfaceTension::AddBTransCB

Parameters

rB	Strain rate matrix
rShapeDeriv	Nodal shape funcion derivatives

◆ CalculateC() [1/3]

void Kratos::SurfaceTension< 2, 3 >::CalculateC	(	boost::numeric::ublas::bounded_matrix< double, 3, 3 > &	rC,
		const double	Viscosity
	)

protected

See SurfaceTension::CalculateC

◆ CalculateC() [2/3]

void Kratos::SurfaceTension< 3, 4 >::CalculateC	(	boost::numeric::ublas::bounded_matrix< double, 6, 6 > &	rC,
		const double	Viscosity
	)

protected

See SurfaceTension::CalculateC

◆ CalculateC() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateC	(	BoundedMatrix< double,(TDim TNumNodes)/2,(TDim TNumNodes)/2 > &	rC,
		const double	Viscosity
	)

protectedvirtual

Calculate a matrix that provides the stress given the strain rate.

Unused, left to support derived classes.

See also: SurfaceTension::AddBTransCB. Note that only non-zero terms are written, so the output matrix should be initialized before calling this.

Parameters

rC	Matrix representation of the stress tensor (output)
Viscosity	Effective viscosity, in dynamic units, weighted by the integration point area

◆ CalculateLeftHandSide()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateLeftHandSide	(	MatrixType &	rLeftHandSideMatrix,
		ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Returns a zero matrix of appropiate size (provided for compatibility with scheme)

Parameters

rLeftHandSideMatrix	Local matrix, will be filled with zeros
rCurrentProcessInfo	Process info instance

◆ CalculateLocalSystem()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateLocalSystem	(	MatrixType &	rLeftHandSideMatrix,
		VectorType &	rRightHandSideVector,
		ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Provides local contributions from body forces and OSS projection terms.

This is called during the assembly process and provides the terms of the system that are either constant or computed explicitly (from the 'old' iteration variables). In this case this means the body force terms and the OSS projections, that are treated explicitly.

Parameters

rLeftHandSideMatrix	the elemental left hand side matrix. Not used here, required for compatibility purposes only.
rRightHandSideVector	the elemental right hand side
rCurrentProcessInfo	the current process info

◆ CalculateLocalVelocityContribution()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateLocalVelocityContribution	(	MatrixType &	rDampingMatrix,
		VectorType &	rRightHandSideVector,
		ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Computes the local contribution associated to 'new' velocity and pressure values.

Provides local contributions to the system associated to the velocity and pressure terms (convection, diffusion, pressure gradient/velocity divergence and stabilization).

Parameters

rDampingMatrix	Will be filled with the velocity-proportional "damping" matrix
rRightHandSideVector	the elemental right hand side vector
rCurrentProcessInfo	the current process info instance

◆ CalculateLumpedMassMatrix()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateLumpedMassMatrix	(	MatrixType &	rLHSMatrix,
		const double	Mass
	)

inlineprotected

Add lumped mass matrix.

Adds the lumped mass matrix to an elemental LHS matrix. Note that the time factor (typically 1/(k*Dt) ) is added by the scheme outside the element.

Parameters

rLHSMatrix	The local matrix where the result will be added
Mass	The weight assigned to each node (typically Density * Area / NumNodes or Density*Volume / NumNodes)

◆ CalculateMassMatrix()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateMassMatrix	(	MatrixType &	rMassMatrix,
		ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Computes local contributions to the mass matrix.

Provides the local contributions to the mass matrix, which is defined here as the matrix associated to velocity derivatives. Note that the mass matrix implemented here is lumped.

Parameters

rMassMatrix	Will be filled with the elemental mass matrix
rCurrentProcessInfo	the current process info instance

◆ CalculateRightHandSide()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateRightHandSide	(	VectorType &	rRightHandSideVector,
		ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Provides local contributions from body forces and projections to the RHS.

This is called during the assembly process and provides the RHS terms of the system that are either constant or computed explicitly (from the 'old' iteration variables). In this case this means the body force terms and the OSS projections, that are treated explicitly.

Parameters

rRightHandSideVector	Will be filled with the elemental right hand side
rCurrentProcessInfo	ProcessInfo instance from the ModelPart. It is expected to contain values for OSS_SWITCH, DYNAMIC_TAU and DELTA_TIME

◆ CalculateStaticTau()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateStaticTau	(	double &	TauOne,
		const array_1d< double, 3 > &	rAdvVel,
		const double	ElemSize,
		const double	Density,
		const double	Viscosity
	)

inlineprotectedvirtual

Calculate momentum stabilization parameter (without time term).

Calculates the momentum tau parameter based on a given advective velocity. The dynamic term is not taken into account. This function is intended for error estimation only. In other cases use CalculateTau

Parameters

TauOne	First stabilization parameter (momentum equation)
rAdvVel	advection velocity
ElemSize	Characteristic element length
Density	Density on integrartion point
Viscosity	Dynamic viscosity (mu) on integrartion point

◆ CalculateTau()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::CalculateTau	(	double &	TauOne,
		double &	TauTwo,
		const array_1d< double, 3 > &	rAdvVel,
		const double	ElemSize,
		const double	Density,
		const double	Viscosity,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineprotectedvirtual

Calculate Stabilization parameters.

Calculates both tau parameters based on a given advective velocity. Takes time step and dynamic coefficient from given ProcessInfo instance. ProcessInfo variables DELTA_TIME and DYNAMIC_TAU will be used.

Parameters

TauOne	First stabilization parameter (momentum equation)
TauTwo	Second stabilization parameter (mass equation)
rAdvVel	advection velocity
ElemSize	Characteristic element length
Density	Density on integrartion point
Viscosity	Dynamic viscosity (mu) on integrartion point
rCurrentProcessInfo	Process info instance

◆ Check()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

int Kratos::SurfaceTension< TDim, TNumNodes >::Check ( const ProcessInfo & rCurrentProcessInfo )

inlineoverride

Checks the input and that all required Kratos variables have been registered.

This function provides the place to perform checks on the completeness of the input. It is designed to be called only once (or anyway, not often) typically at the beginning of the calculations, so to verify that nothing is missing from the input or that no common error is found.

Parameters

rCurrentProcessInfo The ProcessInfo of the ModelPart that contains this element.

Returns: 0 if no errors were found.

◆ ConsistentMassCoef() [1/3]

double Kratos::SurfaceTension< 2, 3 >::ConsistentMassCoef ( const double Area )

protected

◆ ConsistentMassCoef() [2/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::ConsistentMassCoef ( const double Area )

protected

◆ ConsistentMassCoef() [3/3]

double Kratos::SurfaceTension< 3, 4 >::ConsistentMassCoef ( const double Volume )

protected

◆ Create() [1/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Element::Pointer Kratos::SurfaceTension< TDim, TNumNodes >::Create	(	IndexType	NewId,
		GeometryType::Pointer	pGeom,
		PropertiesType::Pointer	pProperties
	)		const

inlineoverridevirtual

It creates a new element pointer.

Parameters

NewId	the ID of the new element
pGeom	the geometry to be employed
pProperties	the properties assigned to the new element

Returns: a Pointer to the new element

Reimplemented from Kratos::Element.

◆ Create() [2/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Element::Pointer Kratos::SurfaceTension< TDim, TNumNodes >::Create	(	IndexType	NewId,
		NodesArrayType const &	ThisNodes,
		PropertiesType::Pointer	pProperties
	)		const

inlineoverridevirtual

Create a new element of this type.

Returns a pointer to a new SurfaceTension element, created using given input

Parameters

NewId	the ID of the new element
ThisNodes	the nodes of the new element
pProperties	the properties assigned to the new element

Returns: a Pointer to the new element

Reimplemented from Kratos::Element.

◆ CrossProduct3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

array_1d<double,3> Kratos::SurfaceTension< TDim, TNumNodes >::CrossProduct3D	(	const array_1d< double, 3 > &	a,
		const array_1d< double, 3 > &	b
	)

inlineprotected

◆ DotProduct2D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::DotProduct2D	(	const array_1d< double, 2 > &	a,
		const array_1d< double, 2 > &	b
	)

inlineprotected

◆ DotProduct3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::DotProduct3D	(	const array_1d< double, 3 > &	a,
		const array_1d< double, 3 > &	b
	)

inlineprotected

◆ EffectiveViscosity()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual double Kratos::SurfaceTension< TDim, TNumNodes >::EffectiveViscosity	(	double	Density,
		const array_1d< double, TNumNodes > &	rN,
		const BoundedMatrix< double, TNumNodes, TDim > &	rDN_DX,
		double	ElemSize,
		const ProcessInfo &	rProcessInfo
	)

inlineprotectedvirtual

EffectiveViscosity Calculate the viscosity at given integration point, using Smagorinsky if enabled.

The Smagorinsky model is used only if the C_SMAGORINSKY is defined on the elemental data container.

Note: : This function is redefined when using Non-Newtonian constitutive models. It is important to keep its signature, otherwise non-Newtonian models will stop working.

Parameters

Density	The fluid's density at the integration point.
rN	Nodal shape functions evaluated at the integration points (area coordinates for the point).
rDN_DX	Shape function derivatives at the integration point.
ElemSize	Representative length of the element (used only for Smagorinsky).
rProcessInfo	ProcessInfo instance passed from the ModelPart.

Returns: Effective viscosity, in dynamic units (Pa*s or equivalent).

◆ ElementSize() [1/3]

double Kratos::SurfaceTension< 2, 3 >::ElementSize ( const double Area )

protected

The size of the 2D element is estimated as the diameter of a circle of the same area. Area = Pi * (h/2)^2

See also: SurfaceTension::ElementSize

◆ ElementSize() [2/3]

double Kratos::SurfaceTension< 3, 4 >::ElementSize ( const double Volume )

protected

The size of the 3D element is estimated as the diameter of the sphere circumscribed to a regular tetrahedron with the same volume.

See also: SurfaceTension::ElementSize

◆ ElementSize() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::ElementSize ( const double )

protected

Return an estimate for the element size h, used to calculate the stabilization parameters.

Estimate the element size from its area or volume, required to calculate stabilization parameters. Note that its implementation is different for 2D or 3D elements.

See also: ST2D, ST3D for actual implementation

Parameters

Volume (in 3D) or Area (in 2D) of the element

Returns: Element size h

◆ EquationIdVector() [1/3]

void Kratos::SurfaceTension< 2 >::EquationIdVector	(	EquationIdVectorType &	rResult,
		ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::EquationIdVector

◆ EquationIdVector() [2/3]

void Kratos::SurfaceTension< 3 >::EquationIdVector	(	EquationIdVectorType &	rResult,
		ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::EquationIdVector

◆ EquationIdVector() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::EquationIdVector	(	EquationIdVectorType &	rResult,
		ProcessInfo &	rCurrentProcessInfo
	)

override

Provides the global indices for each one of this element's local rows.

this determines the elemental equation ID vector for all elemental DOFs

Parameters

rResult	A vector containing the global Id of each row
rCurrentProcessInfo	the current process info object (unused)

◆ EquivalentStrainRate() [1/3]

double Kratos::SurfaceTension< 2, 3 >::EquivalentStrainRate ( const boost::numeric::ublas::bounded_matrix< double, 3, 2 > & rDN_DX ) const

protected

◆ EquivalentStrainRate() [2/3]

double Kratos::SurfaceTension< 3, 4 >::EquivalentStrainRate ( const boost::numeric::ublas::bounded_matrix< double, 4, 3 > & rDN_DX ) const

protected

◆ EquivalentStrainRate() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::EquivalentStrainRate ( const BoundedMatrix< double, TNumNodes, TDim > & rDN_DX ) const

protected

EquivalentStrainRate Calculate the second invariant of the strain rate tensor GammaDot = (2SijSij)^0.5.

Note: Our implementation of non-Newtonian consitutive models such as Bingham relies on this funcition being defined on all fluid elements.

Parameters

rDN_DX Shape function derivatives at the integration point.

Returns: GammaDot = (2SijSij)^0.5.

◆ EvaluateInPoint() [1/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::EvaluateInPoint	(	array_1d< double, 3 > &	rResult,
		const Variable< array_1d< double, 3 > > &	rVariable,
		const array_1d< double, TNumNodes > &	rShapeFunc
	)

inlineprotectedvirtual

Write the value of a variable at a point inside the element to a double.

Evaluate a scalar variable in the point where the form functions take the values given by rShapeFunc and write the result to rResult. This is an auxiliary function used to compute values in integration points.

Parameters

rResult	The double where the value will be added to
rVariable	The nodal variable to be read
rShapeFunc	The values of the form functions in the point

◆ EvaluateInPoint() [2/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::EvaluateInPoint	(	double &	rResult,
		const Variable< double > &	rVariable,
		const array_1d< double, TNumNodes > &	rShapeFunc
	)

inlineprotectedvirtual

Write the value of a variable at a point inside the element to a double.

Evaluate a scalar variable in the point where the form functions take the values given by rShapeFunc and write the result to rResult. This is an auxiliary function used to compute values in integration points.

Parameters

rResult	The double where the value will be added to
rVariable	The nodal variable to be read
rShapeFunc	The values of the form functions in the point
Step	The time Step (Defaults to 0 = Current)

◆ FinalizeNonLinearIteration()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::FinalizeNonLinearIteration ( ProcessInfo & rCurrentProcessInfo )

inlineoverride

◆ GetAdvectiveVel() [1/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetAdvectiveVel	(	array_1d< double, 3 > &	rAdvVel,
		const array_1d< double, TNumNodes > &	rShapeFunc
	)

inlineprotected

Write the advective velocity evaluated at this point to an array.

Writes the value of the advective velocity evaluated at a point inside the element to an array_1d

Parameters

rAdvVel	Output array
rShapeFunc	Shape functions evaluated at the point of interest

◆ GetAdvectiveVel() [2/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

virtual void Kratos::SurfaceTension< TDim, TNumNodes >::GetAdvectiveVel	(	array_1d< double, 3 > &	rAdvVel,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const std::size_t	Step
	)

inlineprotectedvirtual

Write the advective velocity evaluated at this point to an array.

Writes the value of the advective velocity evaluated at a point inside the element to an array_1d

Parameters

rAdvVel	Output array
rShapeFunc	Shape functions evaluated at the point of interest
Step	The time Step

◆ GetConvectionOperator()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetConvectionOperator	(	array_1d< double, TNumNodes > &	rResult,
		const array_1d< double, 3 > &	rVelocity,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv
	)

inlineprotected

Write the convective operator evaluated at this point (for each nodal funciton) to an array.

Evaluate the convective operator for each node's shape function at an arbitrary point

Parameters

rResult	Output vector
rVelocity	Velocity evaluated at the integration point
rShapeDeriv	Derivatives of shape functions evaluated at the integration point

See also: GetAdvectiveVel provides rVelocity

◆ GetDofList() [1/3]

void Kratos::SurfaceTension< 2 >::GetDofList	(	DofsVectorType &	rElementalDofList,
		ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::GetDofList

◆ GetDofList() [2/3]

void Kratos::SurfaceTension< 3 >::GetDofList	(	DofsVectorType &	rElementalDofList,
		ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::GetDofList

◆ GetDofList() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetDofList	(	DofsVectorType &	rElementalDofList,
		ProcessInfo &	rCurrentProcessInfo
	)

override

Returns a list of the element's Dofs.

Parameters

ElementalDofList	the list of DOFs
rCurrentProcessInfo	the current process info instance

◆ GetFirstDerivativesVector() [1/3]

void Kratos::SurfaceTension< 2 >::GetFirstDerivativesVector	(	Vector &	Values,
		int	Step
	)

See also: SurfaceTension::GetFirstDerivativesVector

◆ GetFirstDerivativesVector() [2/3]

void Kratos::SurfaceTension< 3 >::GetFirstDerivativesVector	(	Vector &	Values,
		int	Step
	)

See also: SurfaceTension::GetFirstDerivativesVector

◆ GetFirstDerivativesVector() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetFirstDerivativesVector	(	Vector &	Values,
		int	Step = `0`
	)

override

Returns VELOCITY_X, VELOCITY_Y, (VELOCITY_Z,) PRESSURE for each node.

Parameters

Values	Vector of nodal unknowns
Step	Get result from 'Step' steps back, 0 is current step. (Must be smaller than buffer size)

◆ GetSecondDerivativesVector() [1/3]

void Kratos::SurfaceTension< 2 >::GetSecondDerivativesVector	(	Vector &	Values,
		int	Step
	)

See also: SurfaceTension::GetSecondDerivativesVector

◆ GetSecondDerivativesVector() [2/3]

void Kratos::SurfaceTension< 3 >::GetSecondDerivativesVector	(	Vector &	Values,
		int	Step
	)

See also: SurfaceTension::GetSecondDerivativesVector

◆ GetSecondDerivativesVector() [3/3]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetSecondDerivativesVector	(	Vector &	Values,
		int	Step = `0`
	)

override

Returns ACCELERATION_X, ACCELERATION_Y, (ACCELERATION_Z,) 0 for each node.

Parameters

Values	Vector of nodal second derivatives
Step	Get result from 'Step' steps back, 0 is current step. (Must be smaller than buffer size)

◆ GetValueOnIntegrationPoints() [1/7]

void Kratos::SurfaceTension< 2 >::GetValueOnIntegrationPoints	(	const Variable< array_1d< double, 3 > > &	rVariable,
		std::vector< array_1d< double, 3 > > &	rOutput,
		const ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::GetValueOnIntegrationPoints

◆ GetValueOnIntegrationPoints() [2/7]

void Kratos::SurfaceTension< 3 >::GetValueOnIntegrationPoints	(	const Variable< array_1d< double, 3 > > &	rVariable,
		std::vector< array_1d< double, 3 > > &	rOutput,
		const ProcessInfo &	rCurrentProcessInfo
	)

See also: SurfaceTension::GetValueOnIntegrationPoints

◆ GetValueOnIntegrationPoints() [3/7]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetValueOnIntegrationPoints	(	const Variable< array_1d< double, 3 > > &	rVariable,
		std::vector< array_1d< double, 3 > > &	rOutput,
		const ProcessInfo &	rCurrentProcessInfo
	)

override

Obtain an array_1d<double,3> elemental variable, evaluated on gauss points.

If the variable is VORTICITY, computes the vorticity (rotational of the velocity) based on the current velocity values. Otherwise, it assumes that the input variable is an elemental value and retrieves it. Implemented for a single gauss point only.

Parameters

rVariable	Kratos vector variable to get
Output	Will be filled with the values of the variable on integrartion points
rCurrentProcessInfo	Process info instance

◆ GetValueOnIntegrationPoints() [4/7]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetValueOnIntegrationPoints	(	const Variable< array_1d< double, 6 > > &	rVariable,
		std::vector< array_1d< double, 6 > > &	rValues,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning.

◆ GetValueOnIntegrationPoints() [5/7]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetValueOnIntegrationPoints	(	const Variable< double > &	rVariable,
		std::vector< double > &	rValues,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Obtain a double elemental variable, evaluated on gauss points.

If the variable is TAUONE or TAUTWO, calculates the corresponding stabilization parameter for the element, based on rCurrentProcessInfo's DELTA_TIME and DYNAMIC_TAU. If the variable is MU, calculates the effective viscosity at the element center due to Smagorinsky (in 'dynamic' units). Otherwise, it assumes that the input variable is an elemental value and retrieves it. Implemented for a single gauss point only.

Parameters

rVariable	Kratos vector variable to compute
Output	Will be filled with the values of the variable on integrartion points
rCurrentProcessInfo	Process info instance

◆ GetValueOnIntegrationPoints() [6/7]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetValueOnIntegrationPoints	(	const Variable< Matrix > &	rVariable,
		std::vector< Matrix > &	rValues,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning.

◆ GetValueOnIntegrationPoints() [7/7]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::GetValueOnIntegrationPoints	(	const Variable< Vector > &	rVariable,
		std::vector< Vector > &	rValues,
		const ProcessInfo &	rCurrentProcessInfo
	)

inlineoverride

Empty implementation of unused CalculateOnIntegrationPoints overloads to avoid compilation warning.

◆ Info()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

std::string Kratos::SurfaceTension< TDim, TNumNodes >::Info ( ) const

inlineoverridevirtual

Turn back information as a string.

Reimplemented from Kratos::Element.

◆ KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

Kratos::SurfaceTension< TDim, TNumNodes >::KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION ( SurfaceTension< TDim, TNumNodes > )

Pointer definition of SurfaceTension.

◆ ModulatedGradientDiffusion()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::ModulatedGradientDiffusion	(	MatrixType &	rDampingMatrix,
		const BoundedMatrix< double, TNumNodes, TDim > &	rDN_DX,
		const double	Weight
	)

inlineprotected

◆ Norm2D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::Norm2D ( const array_1d< double, 2 > & a )

inlineprotected

◆ Norm3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::Norm3D ( const array_1d< double, 3 > & a )

inlineprotected

◆ NormalizeVec2D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::NormalizeVec2D ( array_1d< double, 2 > & input )

inlineprotected

◆ NormalizeVec3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::NormalizeVec3D ( array_1d< double, 3 > & input )

inlineprotected

◆ OSSMomResidual()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::OSSMomResidual	(	const array_1d< double, 3 > &	rAdvVel,
		const double	Density,
		array_1d< double, 3 > &	rElementalMomRes,
		const array_1d< double, TNumNodes > &	rShapeFunc,
		const BoundedMatrix< double, TNumNodes, TDim > &	rShapeDeriv,
		const double	Weight
	)

inlineprotected

Assemble the contribution from an integration point to the element's residual.

OSS version. Note that rElementalMomRes should be initialized before calling this.

Parameters

rAdvVel	Convection velocity (not including subscale)
Density	Fluid density evaluated at integration point
rElementalMomRes	Result
rShapeFunc	Shape functions evaluated at integration point
rShapeDeriv	Shape function derivatives evaluated at integration point
Weight	Integration point weight (as a fraction of area or volume)

◆ PrintInfo()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

void Kratos::SurfaceTension< TDim, TNumNodes >::PrintInfo ( std::ostream & rOStream ) const

inlineoverridevirtual

Print information about this object.

Reimplemented from Kratos::Element.

◆ SubscaleErrorEstimate()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

double Kratos::SurfaceTension< TDim, TNumNodes >::SubscaleErrorEstimate ( const ProcessInfo & rProcessInfo )

inlineprotected

◆ Vector2D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

array_1d<double,2> Kratos::SurfaceTension< TDim, TNumNodes >::Vector2D	(	const double	x0,
		const double	y0,
		const double	x1,
		const double	y1
	)

inlineprotected

◆ Vector3D()

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

array_1d<double,3> Kratos::SurfaceTension< TDim, TNumNodes >::Vector3D	(	const double	x0,
		const double	y0,
		const double	z0,
		const double	x1,
		const double	y1,
		const double	z1
	)

inlineprotected

Friends And Related Function Documentation

◆ Serializer

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>

friend class Serializer

friend

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

/home/runner/work/Documentation/Documentation/master/applications/ULFapplication/custom_elements/surface_tension.h

Public Member Functions

Protected Member Functions

Type Definitions

Serialization

Specialized implementation of SurfaceTension for functions that depend on TDim

Additional Inherited Members

Detailed Description

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1> class Kratos::SurfaceTension< TDim, TNumNodes >

Member Typedef Documentation

◆ BaseType

◆ DofsArrayType

◆ DofsVectorType

◆ EquationIdVectorType

◆ GeometryType

◆ IndexType

◆ MatrixType

◆ NodesArrayType

◆ NodeType

◆ PropertiesType

◆ ShapeFunctionDerivativesType

◆ ShapeFunctionsType

◆ SizeType

◆ VectorType

Constructor & Destructor Documentation

◆ SurfaceTension() [1/4]

◆ SurfaceTension() [2/4]

◆ SurfaceTension() [3/4]

◆ SurfaceTension() [4/4]

◆ ~SurfaceTension()

Member Function Documentation

◆ AddBTransCB()

◆ AddConsistentMassMatrixContribution()

◆ AddIntegrationPointVelocityContribution()

◆ AddMassStabTerms()

◆ AddMomentumRHS()

◆ AddProjectionResidualContribution()

◆ AddProjectionToRHS()

◆ AddViscousStress2D()

◆ AddViscousTerm() [1/3]

◆ AddViscousTerm() [2/3]

◆ AddViscousTerm() [3/3]

◆ ApplySurfaceTensionContribution()

◆ ApplySurfaceTensionContribution3D()

◆ ASGSMomResidual()

◆ Calculate() [1/2]

◆ Calculate() [2/2]

◆ CalculateB() [1/3]

◆ CalculateB() [2/3]

◆ CalculateB() [3/3]

◆ CalculateC() [1/3]

◆ CalculateC() [2/3]

◆ CalculateC() [3/3]

◆ CalculateLeftHandSide()

◆ CalculateLocalSystem()

◆ CalculateLocalVelocityContribution()

◆ CalculateLumpedMassMatrix()

◆ CalculateMassMatrix()

◆ CalculateRightHandSide()

◆ CalculateStaticTau()

◆ CalculateTau()

◆ Check()

◆ ConsistentMassCoef() [1/3]

◆ ConsistentMassCoef() [2/3]

◆ ConsistentMassCoef() [3/3]

◆ Create() [1/2]

◆ Create() [2/2]

◆ CrossProduct3D()

◆ DotProduct2D()

◆ DotProduct3D()

◆ EffectiveViscosity()

◆ ElementSize() [1/3]

◆ ElementSize() [2/3]

◆ ElementSize() [3/3]

◆ EquationIdVector() [1/3]

◆ EquationIdVector() [2/3]

◆ EquationIdVector() [3/3]

◆ EquivalentStrainRate() [1/3]

◆ EquivalentStrainRate() [2/3]

◆ EquivalentStrainRate() [3/3]

◆ EvaluateInPoint() [1/2]

template<unsigned int TDim, unsigned int TNumNodes = TDim + 1>
class Kratos::SurfaceTension< TDim, TNumNodes >