CrackMeshCut3DUserObject Class Reference

CrackMeshCut3DUserObject: (1) reads in a mesh describing the crack surface, (2) uses the mesh to do initial cutting of 3D elements, and (3) grows the mesh based on prescribed growth functions. More...

#include <CrackMeshCut3DUserObject.h>

Inheritance diagram for CrackMeshCut3DUserObject:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CrackMeshCut3DUserObject (const InputParameters &parameters)
virtual void	initialSetup () override
virtual void	initialize () override
virtual const std::vector< Point >	getCrackFrontPoints (unsigned int num_crack_front_points) const override
	get a set of points along a crack front from a XFEM GeometricCutUserObject More...
virtual const std::vector< RealVectorValue >	getCrackPlaneNormals (unsigned int num_crack_front_points) const override
	get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject More...
virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutEdge > &cut_edges, std::vector< Xfem::CutNode > &cut_nodes) const override
virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutFace > &cut_faces) const override
virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_edges, std::vector< Xfem::CutEdge > &cut_edges) const override
virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_faces, std::vector< Xfem::CutFace > &cut_faces) const override
void	findActiveBoundaryNodes ()
	Find all active boundary nodes in the cutter mesh Find boundary nodes that will grow; nodes outside of the structural mesh are inactive. More...
std::vector< int >	getFrontPointsIndex () const
	Get crack front points in the active segment -1 means inactive; positive is the point's index in the Crack Front Definition starting from 0. More...
virtual unsigned int	getNumberOfCrackFrontPoints () const override
	Return the total number of crack front points. More...
MeshBase &	getCutterMesh () const
	Get a reference to the cutter mesh. More...
virtual void	execute () override
virtual void	threadJoin (const UserObject &y) override
virtual void	finalize () override
virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutEdge > &cut_edges, std::vector< Xfem::CutNode > &cut_nodes) const =0
	Check to see whether a specified 2D element should be cut based on geometric conditions. More...
virtual bool	cutElementByGeometry (const Elem *elem, std::vector< Xfem::CutFace > &cut_faces) const =0
	Check to see whether a specified 3D element should be cut based on geometric conditions. More...
virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_edges, std::vector< Xfem::CutEdge > &cut_edges) const =0
	Check to see whether a fragment of a 2D element should be cut based on geometric conditions. More...
virtual bool	cutFragmentByGeometry (std::vector< std::vector< Point >> &frag_faces, std::vector< Xfem::CutFace > &cut_faces) const =0
	Check to see whether a fragment of a 3D element should be cut based on geometric conditions. More...
unsigned int	getInterfaceID () const
	Get the interface ID for this cutting object. More...
void	setInterfaceID (unsigned int interface_id)
	Set the interface ID for this cutting object. More...
bool	shouldHealMesh () const
	Should the elements cut by this cutting object be healed in the current time step? More...
virtual CutSubdomainID	getCutSubdomainID (const Node *) const
	Get CutSubdomainID telling which side the node belongs to relative to the cut. More...
CutSubdomainID	getCutSubdomainID (const Elem *elem) const
	Get the CutSubdomainID for the given element. More...
bool	usesMesh () const
	Getter for if a cutter mesh is used in a derived class. More...
virtual Real	spatialValue (const Point &) const
virtual const std::vector< Point >	spatialPoints () const
void	setPrimaryThreadCopy (UserObject *primary)
UserObject *	primaryThreadCopy ()
SubProblem &	getSubProblem () const
bool	shouldDuplicateInitialExecution () const
void	gatherSum (T &value)
void	gatherMax (T &value)
void	gatherMin (T &value)
void	gatherProxyValueMax (T1 &proxy, T2 &value)
void	gatherProxyValueMin (T1 &proxy, T2 &value)
std::set< UserObjectName >	getDependObjects () const
const std::set< std::string > &	getRequestedItems () override
const std::set< std::string > &	getSuppliedItems () override
unsigned int	systemNumber () const
virtual bool	needThreadedCopy () const
virtual bool	enabled () const
std::shared_ptr< MooseObject >	getSharedPtr ()
std::shared_ptr< const MooseObject >	getSharedPtr () const
bool	isKokkosObject () const
MooseApp &	getMooseApp () const
const std::string &	type () const
const std::string &	name () const
std::string	typeAndName () const
MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const
MooseObjectName	uniqueName () const
const InputParameters &	parameters () const
const hit::Node *	getHitNode () const
bool	hasBase () const
const std::string &	getBase () const
const T &	getParam (const std::string &name) const
std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const
const T *	queryParam (const std::string &name) const
const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const
T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
bool	haveParameter (const std::string &name) const
bool	isParamValid (const std::string &name) const
bool	isParamSetByUser (const std::string &name) const
void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
void	paramError (const std::string &param, Args... args) const
void	paramWarning (const std::string &param, Args... args) const
void	paramWarning (const std::string &param, Args... args) const
void	paramInfo (const std::string &param, Args... args) const
std::string	messagePrefix (const bool hit_prefix=true) const
std::string	errorPrefix (const std::string &) const
void	mooseError (Args &&... args) const
void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
void	mooseErrorNonPrefixed (Args &&... args) const
void	mooseWarning (Args &&... args) const
void	mooseWarning (Args &&... args) const
void	mooseWarningNonPrefixed (Args &&... args) const
void	mooseWarningNonPrefixed (Args &&... args) const
void	mooseDeprecated (Args &&... args) const
void	mooseDeprecated (Args &&... args) const
void	mooseDeprecatedNoTrace (Args &&... args) const
void	mooseInfo (Args &&... args) const
void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
std::string	getDataFileName (const std::string &param) const
std::string	getDataFileNameByName (const std::string &relative_path) const
std::string	getDataFilePath (const std::string &relative_path) const
virtual void	timestepSetup ()
virtual void	jacobianSetup ()
virtual void	residualSetup ()
virtual void	subdomainSetup ()
virtual void	customSetup (const ExecFlagType &)
const ExecFlagEnum &	getExecuteOnEnum () const
UserObjectName	getUserObjectName (const std::string &param_name) const
const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const
const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
const UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const
const UserObjectBase &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
const std::vector< SubdomainName > &	blocks () const
unsigned int	numBlocks () const
virtual const std::set< SubdomainID > &	blockIDs () const
unsigned int	blocksMaxDimension () const
bool	hasBlocks (const SubdomainName &name) const
bool	hasBlocks (const std::vector< SubdomainName > &names) const
bool	hasBlocks (const std::set< SubdomainName > &names) const
bool	hasBlocks (SubdomainID id) const
bool	hasBlocks (const std::vector< SubdomainID > &ids) const
bool	hasBlocks (const std::set< SubdomainID > &ids) const
bool	isBlockSubset (const std::set< SubdomainID > &ids) const
bool	isBlockSubset (const std::vector< SubdomainID > &ids) const
bool	hasBlockMaterialProperty (const std::string &prop_name)
const std::set< SubdomainID > &	meshBlockIDs () const
virtual bool	blockRestricted () const
virtual void	checkVariable (const MooseVariableFieldBase &variable) const
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)
const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)
const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)
const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)
const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name_in)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)
Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)
std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > >	getKokkosBlockMaterialProperty (const MaterialPropertyName &name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()
const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)
const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)
std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)
std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)
std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)
std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)
void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)
void	statefulPropertiesAllowed (bool)
virtual bool	getMaterialPropertyCalled () const
virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const
virtual void	resolveOptionalProperties ()
const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
virtual const VariableValue &	coupledValueByName (const std::string &var_name)
virtual const ArrayVariableValue &	coupledArrayValueByName (const std::string &var_name)
const std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > &	getCoupledVars () const
const std::vector< MooseVariableFieldBase *> &	getCoupledMooseVars () const
const std::vector< MooseVariable *> &	getCoupledStandardMooseVars () const
const std::vector< VectorMooseVariable *> &	getCoupledVectorMooseVars () const
const std::vector< ArrayMooseVariable *> &	getCoupledArrayMooseVars () const
void	addFEVariableCoupleableVectorTag (TagID tag)
void	addFEVariableCoupleableMatrixTag (TagID tag)
std::set< TagID > &	getFEVariableCoupleableVectorTags ()
const std::set< TagID > &	getFEVariableCoupleableVectorTags () const
std::set< TagID > &	getFEVariableCoupleableMatrixTags ()
const std::set< TagID > &	getFEVariableCoupleableMatrixTags () const
auto &	getWritableCoupledVariables () const
bool	hasWritableCoupledVariables () const
const ADVariableValue *	getADDefaultValue (const std::string &var_name) const
const ADVectorVariableValue *	getADDefaultVectorValue (const std::string &var_name) const
const ADVariableGradient &	getADDefaultGradient () const
const ADVectorVariableGradient &	getADDefaultVectorGradient () const
const ADVariableSecond &	getADDefaultSecond () const
const ADVectorVariableCurl &	getADDefaultCurl () const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::VariableGradient	kokkosCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::VariableValue	kokkosCoupledValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledValues (const std::string &var_name) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradient (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradients (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValues (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValues (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledValuesOld (const std::string &var_name) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOld (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOld (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValuesOld (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledValuesOlder (const std::string &var_name) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableGradient	kokkosCoupledGradientsOlder (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalValuesOlder (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledDofValuesOlder (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledDot (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledDots (const std::string &var_name) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosCoupledNodalDots (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValues (const std::string &var_name) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradient (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradients (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValues (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValues (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValuesOld (const std::string &var_name) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientsOld (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValuesOld (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValuesOld (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledValuesOlder (const std::string &var_name) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableGradient	kokkosADCoupledGradientsOlder (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalValuesOlder (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDofValuesOlder (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDot (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledDots (const std::string &var_name) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::ADVariableValue	kokkosADCoupledNodalDots (const std::string &var_name) const
Moose::Kokkos::Scalar< const Real >	kokkosCoupledDotDu (const std::string &var_name, unsigned int comp=0) const
Moose::Kokkos::VariableValue	kokkosZeroValue () const
Moose::Kokkos::VariableGradient	kokkosZeroGradient () const
Moose::Kokkos::VariableValue	kokkosZeroNodalValue () const
const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()
const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const
const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const
const std::set< MooseVariableFieldBase *> &	getMooseVariableDependencies () const
std::set< MooseVariableFieldBase *>	checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})
std::set< MooseVariableFieldBase *>	checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
void	addMooseVariableDependency (MooseVariableFieldBase *var)
void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
bool	isImplicit ()
Moose::StateArg	determineState () const
void	setRandomResetFrequency (ExecFlagType exec_flag)
unsigned long	getRandomLong () const
Real	getRandomReal () const
unsigned int	getSeed (std::size_t id)
unsigned int	getMasterSeed () const
bool	isNodal () const
ExecFlagType	getResetOnTime () const
void	setRandomDataPointer (RandomData *random_data)
virtual unsigned int	getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const
virtual unsigned int	getElementIDIndexByName (const std::string &id_name) const
virtual const dof_id_type &	getElementID (const std::string &id_parameter_name, unsigned int comp=0) const
dof_id_type	getElementID (const Elem *elem, unsigned int elem_id_index) const
virtual const dof_id_type &	getElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const
virtual const dof_id_type &	getElementIDByName (const std::string &id_name) const
virtual const dof_id_type &	getElementIDNeighborByName (const std::string &id_name) const
bool	hasElementID (const std::string &id_name) const
dof_id_type	maxElementID (unsigned int elem_id_index) const
dof_id_type	minElementID (unsigned int elem_id_index) const
bool	areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const
std::unordered_map< dof_id_type, std::set< dof_id_type > >	getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const
std::set< dof_id_type >	getAllElemIDs (unsigned int elem_id_index) const
std::set< dof_id_type >	getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObject (const std::string &param_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
bool	hasUserObjectByName (const UserObjectName &object_name) const
const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)
const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)
const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)
MaterialBase &	getMaterial (const std::string &name)
MaterialBase &	getMaterial (const std::string &name)
MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)
MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)
bool	hasMaterialProperty (const std::string &name)
bool	hasMaterialProperty (const std::string &name)
bool	hasMaterialPropertyByName (const std::string &name)
bool	hasMaterialPropertyByName (const std::string &name)
bool	hasADMaterialProperty (const std::string &name)
bool	hasADMaterialProperty (const std::string &name)
bool	hasADMaterialPropertyByName (const std::string &name)
bool	hasADMaterialPropertyByName (const std::string &name)
bool	hasKokkosMaterialProperty (const std::string &name)
bool	hasKokkosMaterialProperty (const std::string &name)
bool	hasKokkosMaterialPropertyByName (const std::string &name)
bool	hasKokkosMaterialPropertyByName (const std::string &name)
bool	hasGenericMaterialProperty (const std::string &name)
bool	hasGenericMaterialProperty (const std::string &name)
bool	hasGenericMaterialPropertyByName (const std::string &name)
bool	hasGenericMaterialPropertyByName (const std::string &name)
const Function &	getFunction (const std::string &name) const
const Function &	getFunctionByName (const FunctionName &name) const
bool	hasFunction (const std::string &param_name) const
bool	hasFunctionByName (const FunctionName &name) const
Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const
const T &	getKokkosFunction (const std::string &name) const
Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const
const T &	getKokkosFunctionByName (const FunctionName &name) const
bool	hasKokkosFunction (const std::string &param_name) const
bool	hasKokkosFunctionByName (const FunctionName &name) const
bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
bool	hasPostprocessorByName (const PostprocessorName &name) const
std::size_t	coupledPostprocessors (const std::string &param_name) const
const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
bool	hasVectorPostprocessor (const std::string &param_name) const
bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const
T &	getSampler (const std::string &name)
Sampler &	getSampler (const std::string &name)
T &	getSamplerByName (const SamplerName &name)
Sampler &	getSamplerByName (const SamplerName &name)
virtual void	meshChanged ()
virtual void	meshDisplaced ()
PerfGraph &	perfGraph ()
const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const
virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const
const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const
bool	isVectorPostprocessorDistributed (const std::string &param_name) const
bool	isVectorPostprocessorDistributed (const std::string &param_name) const
bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
const Distribution &	getDistribution (const std::string &name) const
const T &	getDistribution (const std::string &name) const
const Distribution &	getDistribution (const std::string &name) const
const T &	getDistribution (const std::string &name) const
const Distribution &	getDistributionByName (const DistributionName &name) const
const T &	getDistributionByName (const std::string &name) const
const Distribution &	getDistributionByName (const DistributionName &name) const
const T &	getDistributionByName (const std::string &name) const
const Parallel::Communicator &	comm () const
processor_id_type	n_processors () const
processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()
static void	callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
static void	sort (typename std::vector< T > &vector)
static void	sortDFS (typename std::vector< T > &vector)
static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header, NameFunc &&name_func)
static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

Public Attributes
	usingCombinedWarningSolutionWarnings
const ConsoleStream	_console

Static Public Attributes
static const std::string	type_param
static const std::string	name_param
static const std::string	unique_name_param
static const std::string	app_param
static const std::string	moose_base_param
static const std::string	kokkos_object_param
static constexpr PropertyValue::id_type	default_property_id
static constexpr PropertyValue::id_type	zero_property_id
static constexpr auto	SYSTEM
static constexpr auto	NAME

Protected Types
enum	GrowthDirectionEnum { GrowthDirectionEnum::MAX_HOOP_STRESS, GrowthDirectionEnum::FUNCTION }
	Enum to for crack growth direction. More...
enum	GrowthRateEnum { GrowthRateEnum::REPORTER, GrowthRateEnum::FUNCTION }
	Enum to for crack growth rate. More...

Protected Member Functions
virtual bool	intersectWithEdge (const Point &p1, const Point &p2, const std::vector< Point > &_vertices, Point &point) const
	Check if a line intersects with an element. More...
bool	findIntersection (const Point &p1, const Point &p2, const std::vector< Point > &vertices, Point &point) const
	Find directional intersection along the positive extension of the vector from p1 to p2. More...
bool	isInsideEdge (const Point &p1, const Point &p2, const Point &p) const
	Check if point p is inside the edge p1-p2. More...
Real	getRelativePosition (const Point &p1, const Point &p2, const Point &p) const
	Get the relative position of p from p1. More...
bool	isInsideCutPlane (const std::vector< Point > &_vertices, const Point &p) const
	Check if point p is inside a plane. More...
void	findBoundaryNodes ()
	Find boundary nodes of the cutter mesh This is a simple algorithm simply based on the added angle = 360 degrees Works fine for planar cutting surface for curved cutting surface, need to re-work this subroutine to make it more general. More...
void	findBoundaryEdges ()
	Find boundary edges of the cutter mesh. More...
void	sortBoundaryNodes ()
	Sort boundary nodes to be in the right order along the boundary. More...
Real	findDistance (dof_id_type node1, dof_id_type node2)
	Find distance between two nodes. More...
void	refineBoundary ()
	If boundary nodes are too sparse, add nodes in between. More...
void	findActiveBoundaryDirection ()
	Find growth direction at each active node. More...
void	growFront ()
	Grow the cutter mesh. More...
void	sortFrontNodes ()
	Sort the front nodes. More...
void	findFrontIntersection ()
	Find front-structure intersections. More...
void	refineFront ()
	Refine the mesh at the front. More...
void	triangulation ()
	Create tri3 elements between the new front and the old front. More...
void	joinBoundary ()
	Join active boundaries and inactive boundaries to be the new boundary. More...
void	initializeCrackFrontNodes ()
	Determine initial crack front nodes from cutter mesh. More...
virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override
virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
virtual void	addUserObjectDependencyHelper (const UserObjectBase &uo) const override
void	addReporterDependencyHelper (const ReporterName &reporter_name) override
void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const
const ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const
const ReporterName &	getReporterName (const std::string &param_name) const
T &	declareRestartableData (const std::string &data_name, Args &&... args)
ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
const T &	getRestartableData (const std::string &data_name) const
T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
T &	declareRecoverableData (const std::string &data_name, Args &&... args)
T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
std::string	restartableName (const std::string &data_name) const
const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)
const T &	getMeshProperty (const std::string &data_name)
bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const
bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const
bool	hasMeshProperty (const std::string &data_name) const
bool	hasMeshProperty (const std::string &data_name) const
std::string	meshPropertyName (const std::string &data_name) const
PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const
PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
std::string	timedSectionName (const std::string &section_name) const
virtual bool	hasBlockMaterialPropertyHelper (const std::string &prop_name)
void	initializeBlockRestrictable (const MooseObject *moose_object)
void	initializeKokkosBlockRestrictable ()
Moose::CoordinateSystemType	getBlockCoordSystem ()
KOKKOS_FUNCTION dof_id_type	numKokkosBlockElements () const
KOKKOS_FUNCTION dof_id_type	numKokkosBlockNodes () const
KOKKOS_FUNCTION dof_id_type	numKokkosBlockSides () const
KOKKOS_FUNCTION ContiguousElementID	kokkosBlockElementID (Moose::Kokkos::ThreadID tid) const
KOKKOS_FUNCTION ContiguousElementID	kokkosBlockNodeID (Moose::Kokkos::ThreadID tid) const
KOKKOS_FUNCTION auto	kokkosBlockElementSideID (Moose::Kokkos::ThreadID tid) const
virtual void	checkMaterialProperty (const std::string &name, const unsigned int state)
virtual void	getKokkosMaterialPropertyHook (const std::string &, const unsigned int)
void	markMatPropRequested (const std::string &)
MaterialPropertyName	getMaterialPropertyName (const std::string &name) const
void	checkExecutionStage ()
virtual void	coupledCallback (const std::string &, bool) const
virtual bool	isCoupled (const std::string &var_name, unsigned int i=0) const
virtual bool	isCoupledConstant (const std::string &var_name) const
unsigned int	coupledComponents (const std::string &var_name) const
VariableName	coupledName (const std::string &var_name, unsigned int comp=0) const
std::vector< VariableName >	coupledNames (const std::string &var_name) const
virtual unsigned int	coupled (const std::string &var_name, unsigned int comp=0) const
std::vector< unsigned int >	coupledIndices (const std::string &var_name) const
virtual const VariableValue &	coupledValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValues (const std::string &var_name) const
std::vector< const VectorVariableValue *>	coupledVectorValues (const std::string &var_name) const
const GenericVariableValue< is_ad > &	coupledGenericValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericValue (const std::string &var_name, unsigned int comp) const
const GenericVectorVariableValue< is_ad > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const
const GenericVectorVariableValue< false > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
const GenericVectorVariableValue< true > &	coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
std::vector< const GenericVariableValue< is_ad > *>	coupledGenericValues (const std::string &var_name) const
std::vector< const GenericVariableValue< false > *>	coupledGenericValues (const std::string &var_name) const
std::vector< const GenericVariableValue< true > *>	coupledGenericValues (const std::string &var_name) const
const GenericVariableValue< is_ad > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< is_ad > &	coupledGenericDot (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< is_ad > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
virtual const VariableValue &	coupledValueLower (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableValue *>	adCoupledValues (const std::string &var_name) const
const ADVariableValue &	adCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableValue &	adCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVectorVariableValue *>	adCoupledVectorValues (const std::string &var_name) const
virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const
virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const
std::vector< const ArrayVariableValue *>	coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const
virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableGradient &	coupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, TagID tag) const
std::vector< const VariableGradient *>	coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const
virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const
std::vector< const ArrayVariableGradient *>	coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const
virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
const ArrayVariableValue &	coupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const
virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
virtual const VariableValue &	coupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, TagID tag) const
std::vector< const VariableValue *>	coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const
virtual const VectorVariableValue &	coupledVectorValue (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValue (const std::string &var_name, unsigned int comp=0) const
std::vector< const ArrayVariableValue *>	coupledArrayValues (const std::string &var_name) const
MooseWritableVariable &	writableVariable (const std::string &var_name, unsigned int comp=0)
virtual VariableValue &	writableCoupledValue (const std::string &var_name, unsigned int comp=0)
void	checkWritableVar (MooseWritableVariable *var)
virtual const VariableValue &	coupledValueOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValuesOld (const std::string &var_name) const
std::vector< const VectorVariableValue *>	coupledVectorValuesOld (const std::string &var_name) const
virtual const VariableValue &	coupledValueOlder (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledValuesOlder (const std::string &var_name) const
virtual const VariableValue &	coupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradient (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableGradient *>	coupledGradients (const std::string &var_name) const
const ADVariableGradient &	adCoupledGradient (const std::string &var_name, unsigned int comp=0) const
const ADVariableGradient &	adCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableGradient *>	adCoupledGradients (const std::string &var_name) const
const GenericVariableGradient< is_ad > &	coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const
const GenericVariableGradient< false > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const
const GenericVariableGradient< true > &	coupledGenericGradient (const std::string &var_name, unsigned int comp) const
std::vector< const GenericVariableGradient< is_ad > *>	coupledGenericGradients (const std::string &var_name) const
std::vector< const GenericVariableGradient< false > *>	coupledGenericGradients (const std::string &var_name) const
std::vector< const GenericVariableGradient< true > *>	coupledGenericGradients (const std::string &var_name) const
const ADVectorVariableGradient &	adCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
const ADVariableSecond &	adCoupledSecond (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableSecond &	adCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableGradient *>	coupledGradientsOld (const std::string &var_name) const
virtual const VariableGradient &	coupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableGradient &	coupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableGradient &	coupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableGradient &	coupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurl (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurlOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableCurl &	coupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableCurl &	adCoupledCurl (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDiv (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDivOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableDivergence &	coupledDivOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecond (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
virtual const VariableSecond &	coupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledDots (const std::string &var_name) const
virtual const VariableValue &	coupledDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledDot (const std::string &var_name, unsigned int comp=0) const
std::vector< const ADVariableValue *>	adCoupledDots (const std::string &var_name) const
const ADVariableValue &	adCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
const ADVectorVariableValue &	adCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VectorVariableValue &	coupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDot (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const ArrayVariableValue &	coupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDu (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValue (const std::string &var_name, unsigned int comp=0) const
const Moose::ADType< T >::type &	adCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
const T &	coupledNodalDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
virtual const VariableValue &	coupledDofValues (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValues (const std::string &var_name) const
virtual const VariableValue &	coupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValuesOld (const std::string &var_name) const
virtual const VariableValue &	coupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
std::vector< const VariableValue *>	coupledAllDofValuesOlder (const std::string &var_name) const
virtual const ArrayVariableValue &	coupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
virtual const ADVariableValue &	adCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adZeroValue () const
const ADVariableGradient &	adZeroGradient () const
const ADVariableSecond &	adZeroSecond () const
const GenericVariableValue< is_ad > &	genericZeroValue ()
const GenericVariableValue< false > &	genericZeroValue ()
const GenericVariableValue< true > &	genericZeroValue ()
const GenericVariableGradient< is_ad > &	genericZeroGradient ()
const GenericVariableGradient< false > &	genericZeroGradient ()
const GenericVariableGradient< true > &	genericZeroGradient ()
const GenericVariableSecond< is_ad > &	genericZeroSecond ()
const GenericVariableSecond< false > &	genericZeroSecond ()
const GenericVariableSecond< true > &	genericZeroSecond ()
bool	checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const
const MooseVariableFieldBase *	getFEVar (const std::string &var_name, unsigned int comp) const
const MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp) const
MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp)
std::vector< const MooseVariableFieldBase *>	getFieldVars (const std::string &var_name) const
const T *	getVarHelper (const std::string &var_name, unsigned int comp) const
T *	getVarHelper (const std::string &var_name, unsigned int comp)
std::vector< const T *>	getVarsHelper (const std::string &var_name) const
MooseVariable *	getVar (const std::string &var_name, unsigned int comp)
const MooseVariable *	getVar (const std::string &var_name, unsigned int comp) const
VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp)
const VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp) const
ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp)
const ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp) const
void	validateExecutionerType (const std::string &name, const std::string &fn_name) const
std::vector< T >	coupledVectorHelper (const std::string &var_name, const Func &func) const
bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const
unsigned int	coupledScalarComponents (const std::string &var_name) const
unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const
libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const
const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const
const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)
const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)
const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
bool	hasReporterValue (const std::string &param_name) const
bool	hasReporterValue (const std::string &param_name) const
bool	hasReporterValue (const std::string &param_name) const
bool	hasReporterValue (const std::string &param_name) const
bool	hasReporterValueByName (const ReporterName &reporter_name) const
bool	hasReporterValueByName (const ReporterName &reporter_name) const
bool	hasReporterValueByName (const ReporterName &reporter_name) const
bool	hasReporterValueByName (const ReporterName &reporter_name) const
const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)
const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)
const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)
const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)
const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)
const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)
void	serialize (std::string &serialized_buffer)
	Methods to pack/unpack the _marked_elems_2d and _marked_elems_3d data into a structure suitable for parallel communication. More...
void	deserialize (std::vector< std::string > &serialized_buffers)

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const unsigned int	_cut_elem_nnode = 3
	The cutter mesh has triangluar elements only. More...
const unsigned int	_cut_elem_dim = 2
MooseMesh &	_mesh
	The structural mesh. More...
CrackFrontDefinition *	_crack_front_definition
	The crack front definition. More...
std::vector< dof_id_type >	_crack_front_points
	updated crack front definition they are in the same order as defined in the input but the number of nodes may increase its difference from _front is that: _front does not necessarily follow the order of crack front definition therefore, _crack_front_points is generated from _front with the order of crack front definition limitation: this approach does not currently support the growth of one crack front into two More...
const GrowthDirectionEnum	_growth_dir_method
	The direction method for growing mesh at the front. More...
const GrowthRateEnum	_growth_increment_method
	The growth increment method for growing mesh at the front. More...
const unsigned int	_elem_dim = 3
	The structural mesh must be 3D only. More...
const Real	_const_intersection = 0.01
	Used to define intersection points. More...
Real	_size_control
	Used for cutter mesh refinement and front advancement. More...
unsigned int	_n_step_growth
	Number of steps to grow the mesh. More...
bool	_stop
	Variables to help control the work flow. More...
bool	_grow
std::vector< dof_id_type >	_boundary
	Boundary nodes of the cutter mesh. More...
std::vector< std::vector< dof_id_type > >	_active_boundary
	Active boundary nodes where growth is allowed. More...
std::vector< unsigned int >	_inactive_boundary_pos
	Inactive boundary. More...
std::vector< dof_id_type >	_tracked_crack_front_points
	Front nodes that are grown from the crack front definition defined in the input therefore, they are (1) in the same order as defined in the input and (2) the number of nodes does not change. More...
bool	_cfd
	is it using the crack_front_definition More...
std::set< Xfem::CutEdge >	_boundary_edges
	Edges at the boundary. More...
std::map< dof_id_type, std::vector< dof_id_type > >	_boundary_map
	A map of boundary nodes and their neighbors. More...
std::vector< std::vector< Point > >	_active_direction
	Growth direction for active boundaries. More...
std::vector< std::vector< dof_id_type > >	_front
	New boundary after growth. More...
bool	_is_mesh_modified
	Indicator that shows if the cutting mesh is modified or not in this calculation step. More...
unsigned int	_num_crack_front_points
	Total number of crack front points in the mesh cutter. More...
const Function *	_func_x
	Parsed functions of front growth. More...
const Function *	_func_y
const Function *	_func_z
const Function *	_func_v
const std::vector< Real > *const	_ki_vpp
	Pointer to fracture integral ki if available. More...
const std::vector< Real > *const	_kii_vpp
	Pointer to fracture integral kii if available. More...
const std::vector< Real > *const	_growth_inc_reporter
	Pointer to reporter with growth increment if available. More...
std::unique_ptr< MeshBase >	_cutter_mesh
	The xfem cutter mesh. More...
std::shared_ptr< XFEM >	_xfem
	Pointer to the XFEM controller object. More...
unsigned int	_interface_id
	Associated interface id. More...
bool	_heal_always
	Heal the mesh. More...
int	_last_step_initialized
	Time step information needed to advance a 3D crack only at the real beginning of a time step. More...
const bool	_uses_mesh
	bool to set if CrackFrontPointsProvider derived objects use a cutter mesh More...
const Elem *const &	_current_elem
const Real &	_current_elem_volume
const MooseArray< Point > &	_q_point
const QBase *const &	_qrule
const MooseArray< Real > &	_JxW
const MooseArray< Real > &	_coord
const Moose::CoordinateSystemType &	_coord_sys
const THREAD_ID	_tid
SubProblem &	_subproblem
FEProblemBase &	_fe_problem
SystemBase &	_sys
Assembly &	_assembly
const bool	_duplicate_initial_execution
std::set< std::string >	_depend_uo
const bool &	_enabled
MooseApp &	_app
Factory &	_factory
ActionFactory &	_action_factory
const std::string &	_type
const std::string &	_name
const InputParameters &	_pars
const ExecFlagEnum &	_execute_enum
const ExecFlagType &	_current_execute_flag
MooseApp &	_restartable_app
const std::string	_restartable_system_name
const THREAD_ID	_restartable_tid
const bool	_restartable_read_only
FEProblemBase &	_mci_feproblem
FEProblemBase &	_mdi_feproblem
MooseApp &	_pg_moose_app
const std::string	_prefix
const MaterialData *	_blk_material_data
const InputParameters &	_mi_params
const std::string	_mi_name
const MooseObjectName	_mi_moose_object_name
FEProblemBase &	_mi_feproblem
SubProblem &	_mi_subproblem
const THREAD_ID	_mi_tid
const bool	_is_kokkos_object
const Moose::MaterialDataType	_material_data_type
MaterialData &	_material_data
bool	_stateful_allowed
bool	_get_material_property_called
std::vector< std::unique_ptr< PropertyValue > >	_default_properties
std::unordered_set< unsigned int >	_material_property_dependencies
const MaterialPropertyName	_get_suffix
const bool	_use_interpolated_state
const InputParameters &	_c_parameters
const std::string &	_c_name
const std::string &	_c_type
FEProblemBase &	_c_fe_problem
const SystemBase *const	_c_sys
std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> >	_coupled_vars
std::vector< MooseVariableFieldBase *>	_coupled_moose_vars
std::vector< MooseVariable *>	_coupled_standard_moose_vars
std::vector< VectorMooseVariable *>	_coupled_vector_moose_vars
std::vector< ArrayMooseVariable *>	_coupled_array_moose_vars
std::vector< MooseVariableField< Real > *>	_coupled_fv_moose_vars
const std::unordered_map< std::string, std::string > &	_new_to_deprecated_coupled_vars
bool	_c_nodal
bool	_c_is_implicit
const bool	_c_allow_element_to_nodal_coupling
THREAD_ID	_c_tid
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > >	_ad_default_value
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > >	_default_vector_value
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > >	_default_array_value
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > >	_ad_default_vector_value
VariableValue	_default_value_zero
VariableGradient	_default_gradient
MooseArray< ADRealVectorValue >	_ad_default_gradient
MooseArray< ADRealTensorValue >	_ad_default_vector_gradient
VariableSecond	_default_second
MooseArray< ADRealTensorValue >	_ad_default_second
MooseArray< ADRealVectorValue >	_ad_default_curl
const VariableValue &	_zero
const VariablePhiValue &	_phi_zero
const MooseArray< ADReal > &	_ad_zero
const VariableGradient &	_grad_zero
const MooseArray< ADRealVectorValue > &	_ad_grad_zero
const VariablePhiGradient &	_grad_phi_zero
const VariableSecond &	_second_zero
const MooseArray< ADRealTensorValue > &	_ad_second_zero
const VariablePhiSecond &	_second_phi_zero
const VectorVariableValue &	_vector_zero
const VectorVariableCurl &	_vector_curl_zero
VectorVariableValue	_default_vector_value_zero
VectorVariableGradient	_default_vector_gradient
VectorVariableCurl	_default_vector_curl
VectorVariableDivergence	_default_div
ArrayVariableValue	_default_array_value_zero
ArrayVariableGradient	_default_array_gradient
bool	_coupleable_neighbor
FEProblemBase &	_sc_fe_problem
const THREAD_ID	_sc_tid
const Real &	_real_zero
const VariableValue &	_scalar_zero
const Point &	_point_zero
const InputParameters &	_ti_params
FEProblemBase &	_ti_feproblem
bool	_is_implicit
Real &	_t
const Real &	_t_old
int &	_t_step
Real &	_dt
Real &	_dt_old
bool	_is_transient
const Parallel::Communicator &	_communicator
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo2D > >	_marked_elems_2d
	Containers with information about all 2D and 3D elements marked for cutting by this object. More...
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > >	_marked_elems_3d

Static Protected Attributes
static const std::string	_interpolated_old
static const std::string	_interpolated_older

Detailed Description

CrackMeshCut3DUserObject: (1) reads in a mesh describing the crack surface, (2) uses the mesh to do initial cutting of 3D elements, and (3) grows the mesh based on prescribed growth functions.

Definition at line 25 of file CrackMeshCut3DUserObject.h.

Member Enumeration Documentation

GrowthDirectionEnum

enum CrackMeshCut3DUserObject::GrowthDirectionEnum

strongprotected

Enum to for crack growth direction.

Enumerator
MAX_HOOP_STRESS
FUNCTION

Definition at line 88 of file CrackMeshCut3DUserObject.h.

  {
    MAX_HOOP_STRESS,
    FUNCTION
  };

GrowthRateEnum

enum CrackMeshCut3DUserObject::GrowthRateEnum

strongprotected

Enum to for crack growth rate.

Enumerator
REPORTER
FUNCTION

Definition at line 97 of file CrackMeshCut3DUserObject.h.

  {
    REPORTER,
    FUNCTION
  };

Constructor & Destructor Documentation

CrackMeshCut3DUserObject()

CrackMeshCut3DUserObject::CrackMeshCut3DUserObject

(

const InputParameters &

parameters

)

Definition at line 60 of file CrackMeshCut3DUserObject.C.

  : MeshCutUserObjectBase(parameters),
    _mesh(_subproblem.mesh()),
    _growth_dir_method(getParam<MooseEnum>("growth_dir_method").getEnum<GrowthDirectionEnum>()),
    _growth_increment_method(
        getParam<MooseEnum>("growth_increment_method").getEnum<GrowthRateEnum>()),
    _n_step_growth(getParam<unsigned int>("n_step_growth")),
    _is_mesh_modified(false),
    _func_x(parameters.isParamValid("growth_direction_x") ? &getFunction("growth_direction_x")
                                                          : nullptr),
    _func_y(parameters.isParamValid("growth_direction_y") ? &getFunction("growth_direction_y")
                                                          : nullptr),
    _func_z(parameters.isParamValid("growth_direction_z") ? &getFunction("growth_direction_z")
                                                          : nullptr),
    _func_v(parameters.isParamValid("growth_rate") ? &getFunction("growth_rate") : nullptr),
    _ki_vpp((_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS)
                ? &getVectorPostprocessorValue(
                      "ki_vectorpostprocessor",
                      getParam<VectorPostprocessorName>("ki_vectorpostprocessor"))
                : nullptr),
    _kii_vpp((_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS)
                 ? &getVectorPostprocessorValue(
                       "kii_vectorpostprocessor",
                       getParam<VectorPostprocessorName>("kii_vectorpostprocessor"))
                 : nullptr),
    _growth_inc_reporter((_growth_increment_method == GrowthRateEnum::REPORTER)
                             ? &getReporterValueByName<std::vector<Real>>(
                                   getParam<ReporterName>("growth_reporter"), REPORTER_MODE_ROOT)
                             : nullptr)
{
  _grow = (_n_step_growth == 0 ? 0 : 1);

  if (_grow)
  {
    if (!isParamValid("size_control"))
      paramError("size_control", "Crack growth needs size control.");

    _size_control = getParam<Real>("size_control");

    if (_growth_dir_method == GrowthDirectionEnum::FUNCTION &&
        (_func_x == nullptr || _func_y == nullptr || _func_z == nullptr))
      mooseError("function is not specified for the function method that defines growth direction");
  }

  if (_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS ||
      _growth_increment_method == GrowthRateEnum::REPORTER)
    _cfd = true;
  else
    _cfd = false;

  if (_grow)
  {
    findBoundaryNodes();
    findBoundaryEdges();
    sortBoundaryNodes();
    findActiveBoundaryNodes();
    initializeCrackFrontNodes();
  }
  // test element type; only tri3 elements are allowed
  for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
  {
    if (cut_elem->n_nodes() != _cut_elem_nnode)
      mooseError("The input cut mesh should include tri elements only!");
    if (cut_elem->dim() != _cut_elem_dim)
      mooseError("The input cut mesh must be 2D elements only!");
  }
}

Member Function Documentation

cutElementByGeometry() [1/4]

bool CrackMeshCut3DUserObject::cutElementByGeometry ( const Elem *  elem, std::vector< Xfem::CutEdge > &  cut_edges, std::vector< Xfem::CutNode > &  cut_nodes  ) const

overridevirtual

Definition at line 175 of file CrackMeshCut3DUserObject.C.

{
  mooseError("invalid method for 3D mesh cutting");
  return false;
}

cutElementByGeometry() [2/4]

bool CrackMeshCut3DUserObject::cutElementByGeometry ( const Elem *  elem, std::vector< Xfem::CutFace > &  cut_faces  ) const

overridevirtual

Definition at line 184 of file CrackMeshCut3DUserObject.C.

{
  bool elem_cut = false;

  if (elem->dim() != _elem_dim)
    mooseError("The structural mesh to be cut by a surface mesh must be 3D!");

  for (unsigned int i = 0; i < elem->n_sides(); ++i)
  {
    // This returns the lowest-order type of side.
    std::unique_ptr<const Elem> curr_side = elem->side_ptr(i);
    if (curr_side->dim() != 2)
      mooseError("In cutElementByGeometry dimension of side must be 2, but it is ",
                 curr_side->dim());
    unsigned int n_edges = curr_side->n_sides();

    std::vector<unsigned int> cut_edges;
    std::vector<Real> cut_pos;

    for (unsigned int j = 0; j < n_edges; j++)
    {
      // This returns the lowest-order type of side.
      std::unique_ptr<const Elem> curr_edge = curr_side->side_ptr(j);
      if (curr_edge->type() != EDGE2)
        mooseError("In cutElementByGeometry face edge must be EDGE2, but type is: ",
                   libMesh::Utility::enum_to_string(curr_edge->type()),
                   " base element type is: ",
                   libMesh::Utility::enum_to_string(elem->type()));
      const Node * node1 = curr_edge->node_ptr(0);
      const Node * node2 = curr_edge->node_ptr(1);

      for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
      {
        std::vector<Point> vertices;

        for (auto & node : cut_elem->node_ref_range())
        {
          Point & this_point = node;
          vertices.push_back(this_point);
        }

        Point intersection;
        if (intersectWithEdge(*node1, *node2, vertices, intersection))
        {
          cut_edges.push_back(j);
          cut_pos.emplace_back(getRelativePosition(*node1, *node2, intersection));
        }
      }
    }

    // if two edges of an element are cut, it is considered as an element being cut
    if (cut_edges.size() == 2)
    {
      elem_cut = true;
      Xfem::CutFace mycut;
      mycut._face_id = i;
      mycut._face_edge.push_back(cut_edges[0]);
      mycut._face_edge.push_back(cut_edges[1]);
      mycut._position.push_back(cut_pos[0]);
      mycut._position.push_back(cut_pos[1]);
      cut_faces.push_back(mycut);
    }
  }
  return elem_cut;
}

cutElementByGeometry() [3/4]

virtual bool GeometricCutUserObject::cutElementByGeometry ( const Elem *  elem, std::vector< Xfem::CutEdge > &  cut_edges, std::vector< Xfem::CutNode > &  cut_nodes  ) const

pure virtualinherited

Check to see whether a specified 2D element should be cut based on geometric conditions.

Parameters

elem	Pointer to the libMesh element to be considered for cutting
cut_edges	Data structure filled with information about edges to be cut
cut_nodes	Data structure filled with information about nodes to be cut

Returns

bool true if element is to be cut

Implemented in GeometricCut3DUserObject.

Referenced by GeometricCutUserObject::execute().

cutElementByGeometry() [4/4]

virtual bool GeometricCutUserObject::cutElementByGeometry ( const Elem *  elem, std::vector< Xfem::CutFace > &  cut_faces  ) const

pure virtualinherited

Check to see whether a specified 3D element should be cut based on geometric conditions.

Parameters

elem	Pointer to the libMesh element to be considered for cutting
cut_faces	Data structure filled with information about edges to be cut

Returns

bool true if element is to be cut

Implemented in GeometricCut3DUserObject.

cutFragmentByGeometry() [1/4]

bool CrackMeshCut3DUserObject::cutFragmentByGeometry ( std::vector< std::vector< Point >> &  frag_edges, std::vector< Xfem::CutEdge > &  cut_edges  ) const

overridevirtual

Definition at line 255 of file CrackMeshCut3DUserObject.C.

{
  mooseError("invalid method for 3D mesh cutting");
  return false;
}

cutFragmentByGeometry() [2/4]

bool CrackMeshCut3DUserObject::cutFragmentByGeometry ( std::vector< std::vector< Point >> &  frag_faces, std::vector< Xfem::CutFace > &  cut_faces  ) const

overridevirtual

Definition at line 263 of file CrackMeshCut3DUserObject.C.

{
  // TODO: Need this for branching in 3D
  mooseError("cutFragmentByGeometry not yet implemented for 3D mesh cutting");
  return false;
}

cutFragmentByGeometry() [3/4]

virtual bool GeometricCutUserObject::cutFragmentByGeometry ( std::vector< std::vector< Point >> &  frag_edges, std::vector< Xfem::CutEdge > &  cut_edges  ) const

pure virtualinherited

Check to see whether a fragment of a 2D element should be cut based on geometric conditions.

Parameters

frag_edges	Data structure defining the current fragment to be considered
cut_edges	Data structure filled with information about fragment edges to be cut

Returns

bool true if fragment is to be cut

Implemented in GeometricCut3DUserObject.

Referenced by GeometricCutUserObject::execute().

cutFragmentByGeometry() [4/4]

virtual bool GeometricCutUserObject::cutFragmentByGeometry ( std::vector< std::vector< Point >> &  frag_faces, std::vector< Xfem::CutFace > &  cut_faces  ) const

pure virtualinherited

Check to see whether a fragment of a 3D element should be cut based on geometric conditions.

Parameters

frag_faces	Data structure defining the current fragment to be considered
cut_faces	Data structure filled with information about fragment faces to be cut

Returns

bool true if fragment is to be cut

Implemented in GeometricCut3DUserObject.

deserialize()

void GeometricCutUserObject::deserialize ( std::vector< std::string > &  serialized_buffers )

protectedinherited

Definition at line 219 of file GeometricCutUserObject.C.

Referenced by GeometricCutUserObject::finalize().

{
  mooseAssert(serialized_buffers.size() == _app.n_processors(),
              "Unexpected size of serialized_buffers: " << serialized_buffers.size());

  // The input string stream used for deserialization
  std::istringstream iss;

  // Loop over all datastructures for all processors to perfrom the gather operation
  for (unsigned int rank = 0; rank < serialized_buffers.size(); ++rank)
  {
    // skip the current processor (its data is already in the structures)
    if (rank == processor_id())
      continue;

    // populate the stream with a new buffer and reset stream state
    iss.clear();
    iss.str(serialized_buffers[rank]);

    // Load the communicated data into temporary structures
    std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo2D>> other_marked_elems_2d;
    std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo3D>> other_marked_elems_3d;
    dataLoad(iss, other_marked_elems_2d, this);
    dataLoad(iss, other_marked_elems_3d, this);

    // merge the data in with the current processor's data
    _marked_elems_2d.insert(other_marked_elems_2d.begin(), other_marked_elems_2d.end());
    _marked_elems_3d.insert(other_marked_elems_3d.begin(), other_marked_elems_3d.end());
  }
}

execute()

void GeometricCutUserObject::execute ( )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 64 of file GeometricCutUserObject.C.

{
  if (_current_elem->dim() == 2)
  {
    std::vector<Xfem::CutEdge> elem_cut_edges;
    std::vector<Xfem::CutNode> elem_cut_nodes;
    std::vector<Xfem::CutEdge> frag_cut_edges;
    std::vector<std::vector<Point>> frag_edges;

    EFAElement2D * EFAElem = _xfem->getEFAElem2D(_current_elem);

    // Don't cut again if elem has been already cut twice
    if (!EFAElem->isFinalCut())
    {
      // get fragment edges
      _xfem->getFragmentEdges(_current_elem, EFAElem, frag_edges);

      // mark cut edges for the element and its fragment
      bool cut = cutElementByGeometry(_current_elem, elem_cut_edges, elem_cut_nodes);
      if (EFAElem->numFragments() > 0)
        cut |= cutFragmentByGeometry(frag_edges, frag_cut_edges);

      if (cut)
      {
        Xfem::GeomMarkedElemInfo2D gmei2d;
        gmei2d._elem_cut_edges = elem_cut_edges;
        gmei2d._elem_cut_nodes = elem_cut_nodes;
        gmei2d._frag_cut_edges = frag_cut_edges;
        gmei2d._frag_edges = frag_edges;
        _marked_elems_2d[_current_elem->id()].push_back(gmei2d);
      }
    }
  }
  else if (_current_elem->dim() == 3)
  {
    std::vector<Xfem::CutFace> elem_cut_faces;
    std::vector<Xfem::CutFace> frag_cut_faces;
    std::vector<std::vector<Point>> frag_faces;

    EFAElement3D * EFAElem = _xfem->getEFAElem3D(_current_elem);

    // Don't cut again if elem has been already cut twice
    if (!EFAElem->isFinalCut())
    {
      // get fragment edges
      _xfem->getFragmentFaces(_current_elem, EFAElem, frag_faces);

      // mark cut faces for the element and its fragment
      bool cut = cutElementByGeometry(_current_elem, elem_cut_faces);
      // TODO: This would be done for branching, which is not yet supported in 3D
      // if (EFAElem->numFragments() > 0)
      //  cut |= cutFragmentByGeometry(frag_faces, frag_cut_faces, _t);

      if (cut)
      {
        Xfem::GeomMarkedElemInfo3D gmei3d;
        gmei3d._elem_cut_faces = elem_cut_faces;
        gmei3d._frag_cut_faces = frag_cut_faces;
        gmei3d._frag_faces = frag_faces;
        _marked_elems_3d[_current_elem->id()].push_back(gmei3d);
      }
    }
  }
}

finalize()

void GeometricCutUserObject::finalize ( )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 251 of file GeometricCutUserObject.C.

{
  // for single processor runs we do not need to do anything here
  if (_app.n_processors() > 1)
  {
    // create send buffer
    std::string send_buffer;

    // create byte buffers for the streams received from all processors
    std::vector<std::string> recv_buffers;

    // pack the complex datastructures into the string stream
    serialize(send_buffer);

    // broadcast serialized data to and receive from all processors
    _communicator.allgather(send_buffer, recv_buffers);

    // unpack the received data and merge it into the local data structures
    deserialize(recv_buffers);
  }

  for (const auto & it : _marked_elems_2d)
    for (const auto & gmei : it.second)
      _xfem->addGeomMarkedElem2D(it.first, gmei, _interface_id);

  for (const auto & it : _marked_elems_3d)
    for (const auto & gmei : it.second)
      _xfem->addGeomMarkedElem3D(it.first, gmei, _interface_id);

  _marked_elems_2d.clear();
  _marked_elems_3d.clear();
}

findActiveBoundaryDirection()

void CrackMeshCut3DUserObject::findActiveBoundaryDirection ( )

protected

Find growth direction at each active node.

Definition at line 708 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
  mooseAssert(!(_cfd && _active_boundary.size() != 1),
              "crack-front-definition using the cutter mesh only supports one active crack front "
              "segment for now");

  _active_direction.clear();

  for (unsigned int i = 0; i < _active_boundary.size(); ++i)
  {
    std::vector<Point> temp;
    Point dir;

    if (_inactive_boundary_pos.size() != 0)
    {
      for (unsigned int j = 0; j < 3; ++j)
        dir(j) = 0;
      temp.push_back(dir);
    }

    unsigned int i1 = 1;
    unsigned int i2 = _active_boundary[i].size() - 1;
    if (_inactive_boundary_pos.size() == 0)
    {
      i1 = 0;
      i2 = _active_boundary[i].size();
    }

    if (_growth_dir_method == GrowthDirectionEnum::FUNCTION)
      // loop over active front points
      for (unsigned int j = i1; j < i2; ++j)
      {
        Node * this_node = _cutter_mesh->node_ptr(_active_boundary[i][j]);
        mooseAssert(this_node, "Node is NULL");
        Point & this_point = *this_node;
        dir(0) = _func_x->value(0, this_point);
        dir(1) = _func_y->value(0, this_point);
        dir(2) = _func_z->value(0, this_point);

        temp.push_back(dir);
      }
    // determine growth direction based on KI and KII at the crack front
    else if (_growth_dir_method == GrowthDirectionEnum::MAX_HOOP_STRESS)
    {
      mooseAssert(_ki_vpp->size() == _kii_vpp->size(), "KI and KII VPPs must be the same size");
      mooseAssert(_ki_vpp->size() == _active_boundary[0].size(),
                  "the number of crack front nodes in the self-similar method should equal to the "
                  "size of VPP defined at the crack front");
      mooseAssert(_crack_front_points.size() == _active_boundary[0].size(),
                  "the number of crack front nodes should be the same in _crack_front_points and "
                  "_active_boundary[0]");

      // the node order in _active_boundary[0] and _crack_front_points may be the same or opposite,
      // their correspondence is needed
      std::vector<int> index = getFrontPointsIndex();

      for (unsigned int j = i1; j < i2; ++j)
      {
        // growth direction in crack front coord (cfc) system based on the max hoop stress criterion
        // Jiang, Wen, Benjamin W.Spencer, and John E.Dolbow.
        // "Ceramic nuclear fuel fracture modeling with the extended finite "
        // "element method." Engineering Fracture Mechanics 223(2020):106713.
        // https://doi.org/10.1016/j.engfracmech.2019.106713
        // Equation 6
        int ind = index[j];
        Real ki = _ki_vpp->at(ind);
        Real kii = _kii_vpp->at(ind);
        Real sqrt_k = std::sqrt(ki * ki + 8 * kii * kii);

        Real theta_m = 0;
        Real theta_p = 0;
        if (std::abs(kii) > libMesh::TOLERANCE)
        {
          theta_m = 2 * std::atan((ki - sqrt_k) / (4 * kii));
          theta_p = 2 * std::atan((ki + sqrt_k) / (4 * kii));
        }

        // Equation 5 check relative sigma_tt
        Real sigma_tt_m = ki * (3 * std::cos(theta_m / 2) + std::cos(3 * theta_m / 2)) +
                          kii * (-3 * std::sin(theta_m / 2) - 3 * std::sin(3 * theta_m / 2));
        Real sigma_tt_p = ki * (3 * std::cos(theta_p / 2) + std::cos(3 * theta_p / 2)) +
                          kii * (-3 * std::sin(theta_p / 2) - 3 * std::sin(3 * theta_p / 2));
        Real theta;
        if (sigma_tt_m > sigma_tt_p)
          theta = theta_m;
        else
          theta = theta_p;

        // growth direction in crack front coord (cfc) system based on the max hoop stress criterion
        RealVectorValue dir_cfc;

        // growth direction in global coord system based on the max hoop stress criterion
        RealVectorValue dir;

        dir_cfc(0) = std::cos(theta);
        dir_cfc(1) = std::sin(theta);
        dir_cfc(2) = 0;
        dir = _crack_front_definition->rotateFromCrackFrontCoordsToGlobal(dir_cfc, ind);

        temp.push_back(dir);
      }
    }
    else
      mooseError("This growth_dir_method is not pre-defined!");

    if (_inactive_boundary_pos.size() != 0)
    {
      for (unsigned int j = 0; j < 3; ++j)
        dir(j) = 0;
      temp.push_back(dir);
    }

    _active_direction.push_back(temp);
  }

  // normalize the directional vector
  Real maxl = 0;

  for (unsigned int i = 0; i < _active_direction.size(); ++i)
    for (unsigned int j = 0; j < _active_direction[i].size(); ++j)
    {
      Point pt = _active_direction[i][j];
      Real length = std::sqrt(pt * pt);
      if (length > maxl)
        maxl = length;
    }

  for (unsigned int i = 0; i < _active_direction.size(); ++i)
    for (unsigned int j = 0; j < _active_direction[i].size(); ++j)
      _active_direction[i][j] /= maxl;
}

findActiveBoundaryNodes()

void CrackMeshCut3DUserObject::findActiveBoundaryNodes

(

)

Find all active boundary nodes in the cutter mesh Find boundary nodes that will grow; nodes outside of the structural mesh are inactive.

Definition at line 648 of file CrackMeshCut3DUserObject.C.

Referenced by CrackMeshCut3DUserObject(), CrackGrowthReporterBase::getCutterMeshIndices(), and initialize().

{
  _active_boundary.clear();
  _inactive_boundary_pos.clear();

  std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
  pl->enable_out_of_mesh_mode();

  unsigned int n_boundary = _boundary.size();

  // if the node is outside of the structural model, store its position in _boundary to
  // _inactive_boundary_pos
  for (unsigned int j = 0; j < n_boundary; ++j)
  {
    Node * this_node = _cutter_mesh->node_ptr(_boundary[j]);
    mooseAssert(this_node, "Node is NULL");
    Point & this_point = *this_node;

    const Elem * elem = (*pl)(this_point);
    if (elem == nullptr)
      _inactive_boundary_pos.push_back(j);
  }

  unsigned int n_inactive_boundary = _inactive_boundary_pos.size();

  // all nodes are inactive, stop
  if (n_inactive_boundary == n_boundary)
    _stop = 1;

  // find and store active boundary segments in "_active_boundary"
  if (n_inactive_boundary == 0)
    _active_boundary.push_back(_boundary);
  else
  {
    for (unsigned int i = 0; i < n_inactive_boundary - 1; ++i)
    {
      if (_inactive_boundary_pos[i + 1] - _inactive_boundary_pos[i] != 1)
      {
        std::vector<dof_id_type> temp;
        for (unsigned int j = _inactive_boundary_pos[i]; j <= _inactive_boundary_pos[i + 1]; ++j)
        {
          temp.push_back(_boundary[j]);
        }
        _active_boundary.push_back(temp);
      }
    }
    if (_inactive_boundary_pos[n_inactive_boundary - 1] - _inactive_boundary_pos[0] <
        n_boundary - 1)
    {
      std::vector<dof_id_type> temp;
      for (unsigned int j = _inactive_boundary_pos[n_inactive_boundary - 1]; j < n_boundary; ++j)
        temp.push_back(_boundary[j]);
      for (unsigned int j = 0; j <= _inactive_boundary_pos[0]; ++j)
        temp.push_back(_boundary[j]);
      _active_boundary.push_back(temp);
    }
  }
}

findBoundaryEdges()

void CrackMeshCut3DUserObject::findBoundaryEdges ( )

protected

Find boundary edges of the cutter mesh.

Definition at line 393 of file CrackMeshCut3DUserObject.C.

Referenced by CrackMeshCut3DUserObject().

{
  _boundary_edges.clear();

  std::vector<dof_id_type> corner_elem_id;
  unsigned int counter = 0;

  std::vector<dof_id_type> node_id(_cut_elem_nnode);
  std::vector<bool> is_node_on_boundary(_cut_elem_nnode);

  for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
  {
    for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
    {
      node_id[i] = cut_elem->node_ptr(i)->id();
      is_node_on_boundary[i] = (_boundary_map.find(node_id[i]) != _boundary_map.end());
    }

    if (is_node_on_boundary[0] && is_node_on_boundary[1] && is_node_on_boundary[2])
    {
      // this is an element at the corner; all nodes are on the boundary but not all edges are on
      // the boundary
      corner_elem_id.push_back(counter);
    }
    else
    {
      // for other elements, find and store boundary edges
      for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
      {
        // if both nodes on an edge are on the boundary, it is a boundary edge.
        if (is_node_on_boundary[i] && is_node_on_boundary[(i + 1 <= 2) ? i + 1 : 0])
        {
          dof_id_type node1 = node_id[i];
          dof_id_type node2 = node_id[(i + 1 <= 2) ? i + 1 : 0];
          if (node1 > node2)
            std::swap(node1, node2);

          Xfem::CutEdge ce;

          if (node1 > node2)
            std::swap(node1, node2);
          ce._id1 = node1;
          ce._id2 = node2;

          _boundary_edges.insert(ce);
        }
      }
    }
    ++counter;
  }

  // loop over edges in corner elements
  // if an edge is shared by two elements, it is not an boundary edge (is_edge_inside = 1)
  for (unsigned int i = 0; i < corner_elem_id.size(); ++i)
  {
    auto elem_it = _cutter_mesh->elements_begin();

    for (dof_id_type j = 0; j < corner_elem_id[i]; ++j)
      ++elem_it;
    Elem * cut_elem = *elem_it;

    for (unsigned int j = 0; j < _cut_elem_nnode; ++j)
    {
      bool is_edge_inside = 0;

      dof_id_type node1 = cut_elem->node_ptr(j)->id();
      dof_id_type node2 = cut_elem->node_ptr((j + 1 <= 2) ? j + 1 : 0)->id();
      if (node1 > node2)
        std::swap(node1, node2);

      unsigned int counter = 0;
      for (const auto & cut_elem2 : _cutter_mesh->element_ptr_range())
      {
        if (counter != corner_elem_id[i])
        {
          for (unsigned int k = 0; k < _cut_elem_nnode; ++k)
          {
            dof_id_type node3 = cut_elem2->node_ptr(k)->id();
            dof_id_type node4 = cut_elem2->node_ptr((k + 1 <= 2) ? k + 1 : 0)->id();
            if (node3 > node4)
              std::swap(node3, node4);

            if (node1 == node3 && node2 == node4)
            {
              is_edge_inside = 1;
              goto endloop;
            }
          }
        }
        ++counter;
      }
    endloop:
      if (is_edge_inside == 0)
      {
        // store boundary edges
        Xfem::CutEdge ce;

        if (node1 > node2)
          std::swap(node1, node2);
        ce._id1 = node1;
        ce._id2 = node2;

        _boundary_edges.insert(ce);
      }
      else
      {
        // this is not a boundary edge; remove it from existing edge list
        for (auto it = _boundary_edges.begin(); it != _boundary_edges.end();)
        {
          if ((*it)._id1 == node1 && (*it)._id2 == node2)
            it = _boundary_edges.erase(it);
          else
            ++it;
        }
      }
    }
  }
}

findBoundaryNodes()

void CrackMeshCut3DUserObject::findBoundaryNodes ( )

protected

Find boundary nodes of the cutter mesh This is a simple algorithm simply based on the added angle = 360 degrees Works fine for planar cutting surface for curved cutting surface, need to re-work this subroutine to make it more general.

Definition at line 381 of file CrackMeshCut3DUserObject.C.

Referenced by CrackMeshCut3DUserObject().

{
  auto boundary_node_ids = MeshTools::find_boundary_nodes(*_cutter_mesh);
  for (auto it = boundary_node_ids.cbegin(); it != boundary_node_ids.cend(); it++)
  {
    dof_id_type id = *it;
    std::vector<dof_id_type> neighbors;
    _boundary_map[id] = neighbors;
  }
}

findDistance()

Real CrackMeshCut3DUserObject::findDistance ( dof_id_type  node1, dof_id_type  node2  )

protected

Find distance between two nodes.

Definition at line 583 of file CrackMeshCut3DUserObject.C.

Referenced by refineBoundary(), refineFront(), and triangulation().

{
  Node * n1 = _cutter_mesh->node_ptr(node1);
  mooseAssert(n1 != nullptr, "Node is NULL");
  Node * n2 = _cutter_mesh->node_ptr(node2);
  mooseAssert(n2 != nullptr, "Node is NULL");
  Real distance = (*n1 - *n2).norm();
  return distance;
}

findFrontIntersection()

void CrackMeshCut3DUserObject::findFrontIntersection ( )

protected

Find front-structure intersections.

Definition at line 922 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
  ConstBndElemRange & range = *_mesh.getBoundaryElementRange();

  for (unsigned int i = 0; i < _front.size(); ++i)
  {
    if (_front[i].size() >= 2)
    {
      std::vector<Point> pint1;
      std::vector<Point> pint2;
      std::vector<Real> length1;
      std::vector<Real> length2;

      Real node_id = _front[i][0];
      Node * this_node = _cutter_mesh->node_ptr(node_id);
      mooseAssert(this_node, "Node is NULL");
      Point & p2 = *this_node;

      if (_front[i].size() >= 4)
        node_id = _front[i][2];
      else
        node_id = _front[i][1];

      this_node = _cutter_mesh->node_ptr(node_id);
      mooseAssert(this_node, "Node is NULL");
      Point & p1 = *this_node;

      node_id = _front[i].back();
      this_node = _cutter_mesh->node_ptr(node_id);
      mooseAssert(this_node, "Node is NULL");
      Point & p4 = *this_node;

      if (_front[i].size() >= 4)
        node_id = _front[i][_front[i].size() - 3];
      else
        node_id = _front[i][_front[i].size() - 2];

      this_node = _cutter_mesh->node_ptr(node_id);
      mooseAssert(this_node, "Node is NULL");
      Point & p3 = *this_node;

      bool do_inter1 = 1;
      bool do_inter2 = 1;

      std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
      pl->enable_out_of_mesh_mode();
      const Elem * elem = (*pl)(p1);
      if (elem == nullptr)
        do_inter1 = 0;
      elem = (*pl)(p4);
      if (elem == nullptr)
        do_inter2 = 0;

      for (const auto & belem : range)
      {
        Point pt;
        std::vector<Point> vertices;

        elem = belem->_elem;
        std::unique_ptr<const Elem> curr_side = elem->side_ptr(belem->_side);
        for (unsigned int j = 0; j < curr_side->n_nodes(); ++j)
        {
          const Node * node = curr_side->node_ptr(j);
          const Point & this_point = *node;
          vertices.push_back(this_point);
        }

        if (findIntersection(p1, p2, vertices, pt))
        {
          pint1.push_back(pt);
          length1.push_back((pt - p1) * (pt - p1));
        }
        if (findIntersection(p3, p4, vertices, pt))
        {
          pint2.push_back(pt);
          length2.push_back((pt - p3) * (pt - p3));
        }
      }

      if (length1.size() != 0 && do_inter1)
      {
        auto it1 = std::min_element(length1.begin(), length1.end());
        Point inter1 = pint1[std::distance(length1.begin(), it1)];
        inter1 += (inter1 - p1) * _const_intersection;

        Node * this_node = Node::build(inter1, _cutter_mesh->n_nodes()).release();
        _cutter_mesh->add_node(this_node);

        mooseAssert(_cutter_mesh->n_nodes() - 1 > 0, "The cut mesh must be at least one element.");
        unsigned int n = _cutter_mesh->n_nodes() - 1;

        auto it = _front[i].begin();
        _front[i].insert(it, n);

        if (_cfd)
          _tracked_crack_front_points[_tracked_crack_front_points.size() - 1] = n;
      }

      if (length2.size() != 0 && do_inter2)
      {
        auto it2 = std::min_element(length2.begin(), length2.end());
        Point inter2 = pint2[std::distance(length2.begin(), it2)];
        inter2 += (inter2 - p2) * _const_intersection;

        Node * this_node = Node::build(inter2, _cutter_mesh->n_nodes()).release();
        _cutter_mesh->add_node(this_node);

        dof_id_type n = _cutter_mesh->n_nodes() - 1;

        auto it = _front[i].begin();
        unsigned int m = _front[i].size();
        _front[i].insert(it + m, n);

        if (_cfd)
          _tracked_crack_front_points[0] = n;
      }
    }
  }
}

findIntersection()

bool CrackMeshCut3DUserObject::findIntersection ( const Point &  p1, const Point &  p2, const std::vector< Point > &  vertices, Point &  point  ) const

protected

Find directional intersection along the positive extension of the vector from p1 to p2.

Definition at line 303 of file CrackMeshCut3DUserObject.C.

Referenced by findFrontIntersection().

{
  bool has_intersection = false;

  Plane elem_plane(vertices[0], vertices[1], vertices[2]);
  Point point = vertices[0];
  Point normal = elem_plane.unit_normal(point);

  std::array<Real, 3> plane_point = {{point(0), point(1), point(2)}};
  std::array<Real, 3> planenormal = {{normal(0), normal(1), normal(2)}};
  std::array<Real, 3> p_begin = {{p1(0), p1(1), p1(2)}};
  std::array<Real, 3> p_end = {{p2(0), p2(1), p2(2)}};
  std::array<Real, 3> cut_point = {{0.0, 0.0, 0.0}};

  if (Xfem::plane_normal_line_exp_int_3d(
          &plane_point[0], &planenormal[0], &p_begin[0], &p_end[0], &cut_point[0]) == 1)
  {
    Point p(cut_point[0], cut_point[1], cut_point[2]);
    Real dotp = ((p - p1) * (p2 - p1)) / ((p2 - p1) * (p2 - p1));
    if (isInsideCutPlane(vertices, p) && dotp > 1)
    {
      pint = p;
      has_intersection = true;
    }
  }
  return has_intersection;
}

getCrackFrontPoints()

const std::vector< Point > CrackMeshCut3DUserObject::getCrackFrontPoints ( unsigned int  int ) const

overridevirtual

get a set of points along a crack front from a XFEM GeometricCutUserObject

Returns

A vector which contains all crack front points

Implements CrackFrontPointsProvider.

Definition at line 1254 of file CrackMeshCut3DUserObject.C.

{
  std::vector<Point> crack_front_points(number_crack_front_points);
  // number_crack_front_points is updated via
  // _crack_front_definition->updateNumberOfCrackFrontPoints(_crack_front_points.size())
  if (number_crack_front_points != _crack_front_points.size())
    mooseError("Number of nodes in CrackFrontDefinition does not match the number of nodes in the "
               "cutter_mesh.\nCrackFrontDefinition nodes = " +
               Moose::stringify(number_crack_front_points) +
               "\ncutter_mesh nodes = " + Moose::stringify(_crack_front_points.size()));

  for (unsigned int i = 0; i < number_crack_front_points; ++i)
  {
    dof_id_type id = _crack_front_points[i];
    Node * this_node = _cutter_mesh->node_ptr(id);
    mooseAssert(this_node, "Node is NULL");
    Point & this_point = *this_node;
    crack_front_points[i] = this_point;
  }
  return crack_front_points;
}

getCrackPlaneNormals()

const std::vector< RealVectorValue > CrackMeshCut3DUserObject::getCrackPlaneNormals ( unsigned int  int ) const

overridevirtual

get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject

Returns

A vector which contains all crack front normals

Implements CrackFrontPointsProvider.

Definition at line 1277 of file CrackMeshCut3DUserObject.C.

{
  std::vector<RealVectorValue> crack_plane_normals(number_crack_front_points);

  // build the node-to-elems map
  std::unordered_map<dof_id_type, std::vector<dof_id_type>> node_to_elems_map;
  node_to_elems_map.clear();
  for (const auto & elem : _cutter_mesh->element_ptr_range())
    for (auto & node : elem->node_ref_range())
      node_to_elems_map[node.id()].push_back(elem->id());

  // build the elem-to-normal map
  std::unordered_map<dof_id_type, RealVectorValue> elem_to_normal_map;
  elem_to_normal_map.clear();
  for (const auto & elem : _cutter_mesh->element_ptr_range())
  {
    Point & p1 = *elem->node_ptr(0);
    Point & p2 = *elem->node_ptr(1);
    Point & p3 = *elem->node_ptr(2);
    Plane elem_plane(p3, p2, p1); // to match the current normal of 0,0,-1;
    RealVectorValue normal = elem_plane.unit_normal(p1);
    elem_to_normal_map[elem->id()] = normal;
  }

  // for any front node, the normal is averaged based on the normals of all elements sharing this
  // node this code may fail when the front node has no element connected to it, e.g. refinement at
  // step 1 has to be disabled
  for (unsigned int i = 0; i < number_crack_front_points; ++i)
  {
    dof_id_type id = _crack_front_points[i];
    std::vector<dof_id_type> elems = node_to_elems_map[id];
    unsigned int n_elem = elems.size();

    RealVectorValue normal_avr = 0;
    for (unsigned int j = 0; j < n_elem; ++j)
      normal_avr += elem_to_normal_map[elems[j]];
    normal_avr = normal_avr / n_elem;
    crack_plane_normals[i] = normal_avr;
  }
  return crack_plane_normals;
}

getCutSubdomainID() [1/2]

virtual CutSubdomainID GeometricCutUserObject::getCutSubdomainID ( const Node *  ) const

inlinevirtualinherited

Get CutSubdomainID telling which side the node belongs to relative to the cut.

The returned ID contains no physical meaning, but should be consistent throughout the simulation.

Parameters

node	Pointer to the node

Returns

An unsigned int indicating the side

Definition at line 189 of file GeometricCutUserObject.h.

Referenced by CutElementSubdomainModifier::computeSubdomainID(), CutSubdomainIDAux::computeValue(), and XFEM::getCutSubdomainID().

  {
    mooseError("Objects that inherit from GeometricCutUserObject should override the "
               "getCutSubdomainID method");
    return 0;
  }

getCutSubdomainID() [2/2]

CutSubdomainID GeometricCutUserObject::getCutSubdomainID ( const Elem *  elem ) const

inherited

Get the CutSubdomainID for the given element.

Parameters

node	Pointer to the element

Returns

The CutSubdomainID

Definition at line 285 of file GeometricCutUserObject.C.

{
  return _xfem->getCutSubdomainID(this, elem);
}

getCutterMesh()

MeshBase & MeshCutUserObjectBase::getCutterMesh ( ) const

inherited

Get a reference to the cutter mesh.

Returns

reference to the cutter mesh

Definition at line 58 of file MeshCutUserObjectBase.C.

Referenced by XFEMCutMeshOutput::output().

{
  mooseAssert(_cutter_mesh, "MeshCutUserObjectBase::getCutterMesh _cutter_mesh is nullptr");
  return *_cutter_mesh;
}

getFrontPointsIndex()

std::vector< int > CrackMeshCut3DUserObject::getFrontPointsIndex

(

)

const

Get crack front points in the active segment -1 means inactive; positive is the point's index in the Crack Front Definition starting from 0.

Definition at line 1320 of file CrackMeshCut3DUserObject.C.

Referenced by findActiveBoundaryDirection(), CrackGrowthReporterBase::getCutterMeshIndices(), and growFront().

{
  // Crack front definition using the cutter mesh currently only supports one active crack front
  // segment
  unsigned int ibnd = 0;
  unsigned int size_this_segment = _active_boundary[ibnd].size();
  unsigned int n_inactive_nodes = _inactive_boundary_pos.size();

  std::vector<int> index(size_this_segment, -1);

  unsigned int i1 = n_inactive_nodes == 0 ? 0 : 1;
  unsigned int i2 = n_inactive_nodes == 0 ? size_this_segment : size_this_segment - 1;

  // loop over active front points
  for (unsigned int j = i1; j < i2; ++j)
  {
    dof_id_type id = _active_boundary[ibnd][j];
    auto it = std::find(_crack_front_points.begin(), _crack_front_points.end(), id);
    index[j] = std::distance(_crack_front_points.begin(), it);
  }

  return index;
}

getInterfaceID()

unsigned int GeometricCutUserObject::getInterfaceID ( ) const

inlineinherited

Get the interface ID for this cutting object.

Returns

the interface ID

Definition at line 167 of file GeometricCutUserObject.h.

{ return _interface_id; };

getNumberOfCrackFrontPoints()

unsigned int CrackMeshCut3DUserObject::getNumberOfCrackFrontPoints ( ) const

overridevirtual

Return the total number of crack front points.

Returns the number of crack front points for use by CrackFrontDefinition. Overrides the virtual function defined in CrackFrontPointsProvider.

Reimplemented from GeometricCutUserObject.

Definition at line 1345 of file CrackMeshCut3DUserObject.C.

{
  return _num_crack_front_points;
}

getRelativePosition()

Real CrackMeshCut3DUserObject::getRelativePosition ( const Point &  p1, const Point &  p2, const Point &  p  ) const

protected

Get the relative position of p from p1.

Definition at line 343 of file CrackMeshCut3DUserObject.C.

{
  Real full_len = (p2 - p1).norm();
  Real len_p1_p = (p - p1).norm();
  return len_p1_p / full_len;
}

growFront()

void CrackMeshCut3DUserObject::growFront ( )

protected

Grow the cutter mesh.

Definition at line 841 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
  _front.clear();

  for (unsigned int i = 0; i < _active_boundary.size(); ++i)
  {
    std::vector<dof_id_type> temp;

    unsigned int i1 = 1;
    unsigned int i2 = _active_boundary[i].size() - 1;
    if (_inactive_boundary_pos.size() == 0)
    {
      i1 = 0;
      i2 = _active_boundary[i].size();
    }

    std::vector<int> index = getFrontPointsIndex();
    for (unsigned int j = i1; j < i2; ++j)
    {
      Node * this_node = _cutter_mesh->node_ptr(_active_boundary[i][j]);
      mooseAssert(this_node, "Node is NULL");
      Point & this_point = *this_node;
      Point dir = _active_direction[i][j];

      Point x;
      Real growth_increment = 0;
      switch (_growth_increment_method)
      {
        case GrowthRateEnum::FUNCTION:
        {
          growth_increment = _func_v->value(0, Point(0, 0, 0));
          break;
        }
        case GrowthRateEnum::REPORTER:
        {
          int ind = index[j];
          if (index[j] == -1)
            growth_increment = 0;
          else
            growth_increment = _growth_inc_reporter->at(ind);
          break;
        }
        default:
        {
          mooseError("This growth_increment_method is not pre-defined!");
          break;
        }
      }
      for (unsigned int k = 0; k < 3; ++k)
        x(k) = this_point(k) + dir(k) * growth_increment;

      this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
      _cutter_mesh->add_node(this_node);

      dof_id_type id = _cutter_mesh->n_nodes() - 1;
      temp.push_back(id);

      if (_cfd)
      {
        auto it = std::find(_tracked_crack_front_points.begin(),
                            _tracked_crack_front_points.end(),
                            _active_boundary[0][j]);
        if (it != _tracked_crack_front_points.end())
        {
          unsigned int pos = std::distance(_tracked_crack_front_points.begin(), it);
          _tracked_crack_front_points[pos] = id;
        }
      }
    }

    _front.push_back(temp);
  }
}

initialize()

void CrackMeshCut3DUserObject::initialize ( )

overridevirtual

Reimplemented from GeometricCutUserObject.

Definition at line 137 of file CrackMeshCut3DUserObject.C.

{
  _is_mesh_modified = false;

  if (_grow)
  {
    if (_t_step == 1)
      _last_step_initialized = 1;

    _stop = 0;

    if (_t_step > 1 && _t_step != _last_step_initialized)
    {
      _last_step_initialized = _t_step;

      for (unsigned int i = 0; i < _n_step_growth; ++i)
      {
        if (_stop != 1)
        {
          findActiveBoundaryNodes();
          findActiveBoundaryDirection();
          _is_mesh_modified = true;
          growFront();
          sortFrontNodes();
          if (_inactive_boundary_pos.size() != 0)
            findFrontIntersection();
          refineFront();
          triangulation();
          joinBoundary();
        }
      }
    }
  }
  if (_cfd)
    _crack_front_definition->isCutterModified(_is_mesh_modified);
}

initializeCrackFrontNodes()

void CrackMeshCut3DUserObject::initializeCrackFrontNodes ( )

protected

Determine initial crack front nodes from cutter mesh.

Definition at line 573 of file CrackMeshCut3DUserObject.C.

Referenced by CrackMeshCut3DUserObject().

{
  _tracked_crack_front_points.clear();
  _tracked_crack_front_points.assign(_active_boundary[0].rbegin(), _active_boundary[0].rend());

  _num_crack_front_points = _tracked_crack_front_points.size();
  _crack_front_points = _tracked_crack_front_points;
}

initialSetup()

void CrackMeshCut3DUserObject::initialSetup ( )

overridevirtual

Reimplemented from ElementUserObject.

Definition at line 129 of file CrackMeshCut3DUserObject.C.

{
  if (_cfd)
    _crack_front_definition =
        &_fe_problem.getUserObject<CrackFrontDefinition>("crackFrontDefinition");
}

intersectWithEdge()

bool CrackMeshCut3DUserObject::intersectWithEdge ( const Point &  p1, const Point &  p2, const std::vector< Point > &  _vertices, Point &  point  ) const

protectedvirtual

Check if a line intersects with an element.

Definition at line 272 of file CrackMeshCut3DUserObject.C.

{
  bool has_intersection = false;

  Plane elem_plane(vertices[0], vertices[1], vertices[2]);
  Point point = vertices[0];
  Point normal = elem_plane.unit_normal(point);

  std::array<Real, 3> plane_point = {{point(0), point(1), point(2)}};
  std::array<Real, 3> planenormal = {{normal(0), normal(1), normal(2)}};
  std::array<Real, 3> edge_point1 = {{p1(0), p1(1), p1(2)}};
  std::array<Real, 3> edge_point2 = {{p2(0), p2(1), p2(2)}};
  std::array<Real, 3> cut_point = {{0.0, 0.0, 0.0}};

  if (Xfem::plane_normal_line_exp_int_3d(
          &plane_point[0], &planenormal[0], &edge_point1[0], &edge_point2[0], &cut_point[0]) == 1)
  {
    Point temp_p(cut_point[0], cut_point[1], cut_point[2]);
    if (isInsideCutPlane(vertices, temp_p) && isInsideEdge(p1, p2, temp_p))
    {
      pint = temp_p;
      has_intersection = true;
    }
  }
  return has_intersection;
}

isInsideCutPlane()

bool CrackMeshCut3DUserObject::isInsideCutPlane ( const std::vector< Point > &  _vertices, const Point &  p  ) const

protected

Check if point p is inside a plane.

Definition at line 353 of file CrackMeshCut3DUserObject.C.

Referenced by findIntersection(), and intersectWithEdge().

{
  unsigned int n_node = vertices.size();

  Plane elem_plane(vertices[0], vertices[1], vertices[2]);
  Point normal = elem_plane.unit_normal(vertices[0]);

  bool inside = false;
  unsigned int counter = 0;

  for (unsigned int i = 0; i < n_node; ++i)
  {
    unsigned int iplus1 = (i < n_node - 1 ? i + 1 : 0);
    Point middle2p = p - 0.5 * (vertices[i] + vertices[iplus1]);
    const Point side_tang = vertices[iplus1] - vertices[i];
    Point side_norm = side_tang.cross(normal);
    Xfem::normalizePoint(middle2p);
    Xfem::normalizePoint(side_norm);
    if (middle2p * side_norm <= 0.0)
      counter += 1;
  }
  if (counter == n_node)
    inside = true;
  return inside;
}

isInsideEdge()

bool CrackMeshCut3DUserObject::isInsideEdge ( const Point &  p1, const Point &  p2, const Point &  p  ) const

protected

Check if point p is inside the edge p1-p2.

Definition at line 335 of file CrackMeshCut3DUserObject.C.

Referenced by intersectWithEdge().

{
  Real dotp1 = (p1 - p) * (p2 - p1);
  Real dotp2 = (p2 - p) * (p2 - p1);
  return (dotp1 * dotp2 <= 0.0);
}

joinBoundary()

void CrackMeshCut3DUserObject::joinBoundary ( )

protected

Join active boundaries and inactive boundaries to be the new boundary.

Definition at line 1214 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
  if (_inactive_boundary_pos.size() == 0)
  {
    _boundary = _front[0];
    _boundary.pop_back();
    return;
  }

  std::vector<dof_id_type> full_front;

  unsigned int size1 = _active_boundary.size();

  for (unsigned int i = 0; i < size1; ++i)
  {
    unsigned int size2 = _active_boundary[i].size();

    dof_id_type bd1 = _active_boundary[i][size2 - 1];
    dof_id_type bd2 = _active_boundary[i + 1 < size1 ? i + 1 : 0][0];

    full_front.insert(full_front.end(), _front[i].begin(), _front[i].end());

    auto it1 = std::find(_boundary.begin(), _boundary.end(), bd1);
    unsigned int pos1 = std::distance(_boundary.begin(), it1);
    auto it2 = std::find(_boundary.begin(), _boundary.end(), bd2);
    unsigned int pos2 = std::distance(_boundary.begin(), it2);

    if (pos1 <= pos2)
      full_front.insert(full_front.end(), _boundary.begin() + pos1, _boundary.begin() + pos2 + 1);
    else
    {
      full_front.insert(full_front.end(), _boundary.begin() + pos1, _boundary.end());
      full_front.insert(full_front.end(), _boundary.begin(), _boundary.begin() + pos2 + 1);
    }
  }

  _boundary = full_front;
}

refineBoundary()

void CrackMeshCut3DUserObject::refineBoundary ( )

protected

If boundary nodes are too sparse, add nodes in between.

Definition at line 594 of file CrackMeshCut3DUserObject.C.

{
  std::vector<dof_id_type> new_boundary_order(_boundary.begin(), _boundary.end());

  mooseAssert(_boundary.size() >= 2, "Boundary must be at least two nodes");

  for (unsigned int i = _boundary.size() - 1; i >= 1; --i)
  {
    dof_id_type node1 = _boundary[i - 1];
    dof_id_type node2 = _boundary[i];

    Real distance = findDistance(node1, node2);

    if (distance > _size_control)
    {
      unsigned int n = static_cast<unsigned int>(distance / _size_control);
      std::array<Real, 3> x1;
      std::array<Real, 3> x2;

      Node * n1 = _cutter_mesh->node_ptr(node1);
      mooseAssert(n1 != nullptr, "Node is NULL");
      Point & p1 = *n1;
      Node * n2 = _cutter_mesh->node_ptr(node2);
      mooseAssert(n2 != nullptr, "Node is NULL");
      Point & p2 = *n2;

      for (unsigned int j = 0; j < 3; ++j)
      {
        x1[j] = p1(j);
        x2[j] = p2(j);
      }

      for (unsigned int j = 0; j < n; ++j)
      {
        Point x;
        for (unsigned int k = 0; k < 3; ++k)
          x(k) = x2[k] - (x2[k] - x1[k]) * (j + 1) / (n + 1);

        Node * this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
        _cutter_mesh->add_node(this_node);

        dof_id_type id = _cutter_mesh->n_nodes() - 1;
        auto it = new_boundary_order.begin();
        new_boundary_order.insert(it + i, id);
      }
    }
  }

  _boundary = new_boundary_order;
  mooseAssert(_boundary.size() > 0, "Boundary should not have zero size");
  _boundary.pop_back();
}

refineFront()

void CrackMeshCut3DUserObject::refineFront ( )

protected

Refine the mesh at the front.

Definition at line 1043 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
  std::vector<std::vector<dof_id_type>> new_front(_front.begin(), _front.end());

  for (unsigned int ifront = 0; ifront < _front.size(); ++ifront)
  {
    unsigned int i1 = _front[ifront].size() - 1;
    if (_inactive_boundary_pos.size() == 0)
      i1 = _front[ifront].size();

    for (unsigned int i = i1; i >= 1; --i)
    {
      unsigned int i2 = i;
      if (_inactive_boundary_pos.size() == 0)
        i2 = (i <= _front[ifront].size() - 1 ? i : 0);

      dof_id_type node1 = _front[ifront][i - 1];
      dof_id_type node2 = _front[ifront][i2];
      Real distance = findDistance(node1, node2);

      if (distance > _size_control)
      {
        unsigned int n = static_cast<int>(distance / _size_control);
        std::array<Real, 3> x1;
        std::array<Real, 3> x2;

        Node * this_node = _cutter_mesh->node_ptr(node1);
        mooseAssert(this_node, "Node is NULL");
        Point & p1 = *this_node;
        this_node = _cutter_mesh->node_ptr(node2);
        mooseAssert(this_node, "Node is NULL");
        Point & p2 = *this_node;

        for (unsigned int j = 0; j < 3; ++j)
        {
          x1[j] = p1(j);
          x2[j] = p2(j);
        }

        for (unsigned int j = 0; j < n; ++j)
        {
          Point x;
          for (unsigned int k = 0; k < 3; ++k)
            x(k) = x2[k] - (x2[k] - x1[k]) * (j + 1) / (n + 1);

          Node * this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
          _cutter_mesh->add_node(this_node);

          dof_id_type id = _cutter_mesh->n_nodes() - 1;

          auto it = new_front[ifront].begin();
          new_front[ifront].insert(it + i, id);
        }
      }
    }
  }

  _front = new_front;

  if (_cfd)
  {
    if (_front[0][0] == _tracked_crack_front_points[0] &&
        _front[0].back() == _tracked_crack_front_points.back())
      _crack_front_points = _front[0];
    else if (_front[0][0] == _tracked_crack_front_points.back() &&
             _front[0].back() == _tracked_crack_front_points[0])
    {
      _crack_front_points = _front[0];
      std::reverse(_crack_front_points.begin(), _crack_front_points.end());
    }
    else
      mooseError("the crack front and the tracked crack front definition must match in terms of "
                 "their end nodes\n _front[0][0]= " +
                 Moose::stringify(_front[0][0]) + "\n _tracked_crack_front_points[0]= " +
                 Moose::stringify(_tracked_crack_front_points[0]) +
                 "\n _tracked_crack_front_points.back()=" +
                 Moose::stringify(_tracked_crack_front_points.back()));

    _num_crack_front_points = _crack_front_points.size();
    _crack_front_definition->updateNumberOfCrackFrontPoints(_num_crack_front_points);
  }
}

serialize()

void GeometricCutUserObject::serialize ( std::string &  serialized_buffer )

protectedinherited

Methods to pack/unpack the _marked_elems_2d and _marked_elems_3d data into a structure suitable for parallel communication.

Definition at line 206 of file GeometricCutUserObject.C.

Referenced by GeometricCutUserObject::finalize().

{
  // stream for serializing the _marked_elems_2d and _marked_elems_3d data structures to a byte
  // stream
  std::ostringstream oss;
  dataStore(oss, _marked_elems_2d, this);
  dataStore(oss, _marked_elems_3d, this);

  // Populate the passed in string pointer with the string stream's buffer contents
  serialized_buffer.assign(oss.str());
}

setInterfaceID()

void GeometricCutUserObject::setInterfaceID ( unsigned int  interface_id )

inlineinherited

Set the interface ID for this cutting object.

Parameters

the	interface ID

Definition at line 173 of file GeometricCutUserObject.h.

Referenced by XFEM::addGeometricCut().

{ _interface_id = interface_id; };

shouldHealMesh()

bool GeometricCutUserObject::shouldHealMesh ( ) const

inlineinherited

Should the elements cut by this cutting object be healed in the current time step?

Returns

true if the cut element should be healed

Definition at line 180 of file GeometricCutUserObject.h.

Referenced by XFEM::cutMeshWithEFA().

{ return _heal_always; };

sortBoundaryNodes()

void CrackMeshCut3DUserObject::sortBoundaryNodes ( )

protected

Sort boundary nodes to be in the right order along the boundary.

Definition at line 513 of file CrackMeshCut3DUserObject.C.

Referenced by CrackMeshCut3DUserObject().

{
  _boundary.clear();

  for (auto it = _boundary_edges.begin(); it != _boundary_edges.end(); ++it)
  {
    dof_id_type node1 = (*it)._id1;
    dof_id_type node2 = (*it)._id2;

    mooseAssert(_boundary_map.find(node1) != _boundary_map.end(),
                "_boundary_map does not have this key");
    mooseAssert(_boundary_map.find(node2) != _boundary_map.end(),
                "_boundary_map does not have this key");

    _boundary_map.find(node1)->second.push_back(node2);
    _boundary_map.find(node2)->second.push_back(node1);
  }

  auto it = _boundary_map.begin();
  while (it != _boundary_map.end())
  {
    if (it->second.size() != 2)
      mooseError(
          "Boundary nodes in the cutter mesh must have exactly two neighbors; this one has: ",
          it->second.size());
    ++it;
  }

  auto it2 = _boundary_edges.begin();
  dof_id_type node1 = (*it2)._id1;
  dof_id_type node2 = (*it2)._id2;
  _boundary.push_back(node1);
  _boundary.push_back(node2);

  for (unsigned int i = 0; i < _boundary_edges.size() - 1; ++i)
  {
    mooseAssert(_boundary_map.find(node2) != _boundary_map.end(),
                "_boundary_map does not have this key");

    dof_id_type node3 = _boundary_map.find(node2)->second[0];
    dof_id_type node4 = _boundary_map.find(node2)->second[1];

    if (node3 == node1)
    {
      _boundary.push_back(node4);
      node1 = node2;
      node2 = node4;
    }
    else if (node4 == node1)
    {
      _boundary.push_back(node3);
      node1 = node2;
      node2 = node3;
    }
    else
      mooseError("Discontinuity in cutter boundary");
  }
}

sortFrontNodes()

void CrackMeshCut3DUserObject::sortFrontNodes ( )

protected

Sort the front nodes.

Definition at line 916 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{
}

threadJoin()

void GeometricCutUserObject::threadJoin ( const UserObject &  y )

overridevirtualinherited

Implements ElementUserObject.

Definition at line 130 of file GeometricCutUserObject.C.

{
  const auto & gcuo = static_cast<const GeometricCutUserObject &>(y);

  for (const auto & it : gcuo._marked_elems_2d)
  {
    mooseAssert(_marked_elems_2d.find(it.first) == _marked_elems_2d.end(),
                "Element already inserted in map from a different thread");
    _marked_elems_2d[it.first] = it.second;
  }
  for (const auto & it : gcuo._marked_elems_3d)
  {
    mooseAssert(_marked_elems_3d.find(it.first) == _marked_elems_3d.end(),
                "Element already inserted in map from a different thread");
    _marked_elems_3d[it.first] = it.second;
  }
}

triangulation()

void CrackMeshCut3DUserObject::triangulation ( )

protected

Create tri3 elements between the new front and the old front.

Definition at line 1127 of file CrackMeshCut3DUserObject.C.

Referenced by initialize().

{

  mooseAssert(_active_boundary.size() == _front.size(),
              "_active_boundary and _front must be the same size!");

  if (_inactive_boundary_pos.size() == 0)
  {
    _active_boundary[0].push_back(_active_boundary[0][0]);
    _front[0].push_back(_front[0][0]);
  }

  // loop over active segments
  for (unsigned int k = 0; k < _front.size(); ++k)
  {
    unsigned int n1 = _active_boundary[k].size();
    unsigned int n2 = _front[k].size();

    unsigned int i1 = 0;
    unsigned int i2 = 0;

    // stop when all nodes are associated with an element
    while (!(i1 == n1 - 1 && i2 == n2 - 1))
    {
      std::vector<dof_id_type> elem;

      dof_id_type p1 = _active_boundary[k][i1]; // node in the old front
      dof_id_type p2 = _front[k][i2];           // node in the new front

      if (i1 != n1 - 1 && i2 != n2 - 1)
      {
        dof_id_type p3 = _active_boundary[k][i1 + 1]; // next node in the old front
        dof_id_type p4 = _front[k][i2 + 1];           // next node in the new front

        elem.push_back(p1);
        elem.push_back(p2);

        Real d1 = findDistance(p1, p4);
        Real d2 = findDistance(p3, p2);

        if (d1 < d2)
        {
          elem.push_back(p4);
          i2++;
        }

        else
        {
          elem.push_back(p3);
          i1++;
        }
      }

      else if (i1 == n1 - 1)
      {
        dof_id_type p4 = _front[k][i2 + 1]; // next node in the new front

        elem.push_back(p1);
        elem.push_back(p2);
        elem.push_back(p4);
        i2++;
      }

      else if (i2 == n2 - 1)
      {
        dof_id_type p3 = _active_boundary[k][i1 + 1]; // next node in the old front

        elem.push_back(p1);
        elem.push_back(p2);
        elem.push_back(p3);
        i1++;
      }

      Elem * new_elem = Elem::build(TRI3).release();

      for (unsigned int i = 0; i < _cut_elem_nnode; ++i)
      {
        mooseAssert(_cutter_mesh->node_ptr(elem[i]) != nullptr, "Node is NULL");
        new_elem->set_node(i, _cutter_mesh->node_ptr(elem[i]));
      }

      _cutter_mesh->add_elem(new_elem);
    }
  }
}

usesMesh()

bool CrackFrontPointsProvider::usesMesh ( ) const

inlineinherited

Getter for if a cutter mesh is used in a derived class.

Returns

bool indicating if a cutter mesh is used in the derived class

Definition at line 45 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontDefinition::initialSetup().

{ return _uses_mesh; }

validParams()

InputParameters CrackMeshCut3DUserObject::validParams ( )

static

Definition at line 27 of file CrackMeshCut3DUserObject.C.

{
  InputParameters params = MeshCutUserObjectBase::validParams();
  MooseEnum growthDirection("MAX_HOOP_STRESS FUNCTION", "FUNCTION");
  params.addParam<MooseEnum>(
      "growth_dir_method", growthDirection, "choose from FUNCTION, MAX_HOOP_STRESS");
  MooseEnum growthRate("REPORTER FUNCTION", "FUNCTION");
  params.addParam<MooseEnum>(
      "growth_increment_method", growthRate, "choose from FUNCTION, REPORTER");
  params.addParam<FunctionName>("growth_direction_x",
                                "Function defining x-component of crack growth direction");
  params.addParam<FunctionName>("growth_direction_y",
                                "Function defining y-component of crack growth direction");
  params.addParam<FunctionName>("growth_direction_z",
                                "Function defining z-component of crack growth direction");

  params.addParam<VectorPostprocessorName>(
      "ki_vectorpostprocessor", "II_KI_1", "Name of the VectorPostprocessor that computes K_I");
  params.addParam<VectorPostprocessorName>("kii_vectorpostprocessor",
                                           "II_KII_1",
                                           "The name of the vectorpostprocessor that contains KII");
  params.addParam<ReporterName>("growth_reporter",
                                "The name of the Reporter that computes the growth increment");
  params.addParam<FunctionName>("growth_rate", "Function defining crack growth rate");
  params.addParam<Real>(
      "size_control", 0, "Criterion for refining elements while growing the crack");
  params.addParam<unsigned int>("n_step_growth", 0, "Number of steps for crack growth");
  params.addClassDescription("Creates a UserObject for a mesh cutter in 3D problems");
  return params;
}

Member Data Documentation

_active_boundary

std::vector<std::vector<dof_id_type> > CrackMeshCut3DUserObject::_active_boundary

protected

Active boundary nodes where growth is allowed.

Definition at line 125 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), findActiveBoundaryNodes(), getFrontPointsIndex(), growFront(), initializeCrackFrontNodes(), joinBoundary(), and triangulation().

_active_direction

std::vector<std::vector<Point> > CrackMeshCut3DUserObject::_active_direction

protected

Growth direction for active boundaries.

Definition at line 144 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), and growFront().

_boundary

std::vector<dof_id_type> CrackMeshCut3DUserObject::_boundary

protected

Boundary nodes of the cutter mesh.

Definition at line 122 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), joinBoundary(), refineBoundary(), and sortBoundaryNodes().

_boundary_edges

std::set<Xfem::CutEdge> CrackMeshCut3DUserObject::_boundary_edges

protected

Edges at the boundary.

Definition at line 138 of file CrackMeshCut3DUserObject.h.

Referenced by findBoundaryEdges(), and sortBoundaryNodes().

_boundary_map

std::map<dof_id_type, std::vector<dof_id_type> > CrackMeshCut3DUserObject::_boundary_map

protected

A map of boundary nodes and their neighbors.

Definition at line 141 of file CrackMeshCut3DUserObject.h.

Referenced by findBoundaryEdges(), findBoundaryNodes(), and sortBoundaryNodes().

_cfd

bool CrackMeshCut3DUserObject::_cfd

protected

is it using the crack_front_definition

Definition at line 135 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), findActiveBoundaryDirection(), findFrontIntersection(), growFront(), initialize(), initialSetup(), and refineFront().

_const_intersection

const Real CrackMeshCut3DUserObject::_const_intersection = 0.01

protected

Used to define intersection points.

Definition at line 109 of file CrackMeshCut3DUserObject.h.

Referenced by findFrontIntersection().

_crack_front_definition

CrackFrontDefinition* CrackMeshCut3DUserObject::_crack_front_definition

protected

The crack front definition.

Definition at line 78 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), initialize(), initialSetup(), and refineFront().

_crack_front_points

std::vector<dof_id_type> CrackMeshCut3DUserObject::_crack_front_points

protected

updated crack front definition they are in the same order as defined in the input but the number of nodes may increase its difference from _front is that: _front does not necessarily follow the order of crack front definition therefore, _crack_front_points is generated from _front with the order of crack front definition limitation: this approach does not currently support the growth of one crack front into two

Definition at line 85 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), getCrackFrontPoints(), getCrackPlaneNormals(), getFrontPointsIndex(), initializeCrackFrontNodes(), and refineFront().

_cut_elem_dim

const unsigned int CrackMeshCut3DUserObject::_cut_elem_dim = 2

protected

Definition at line 72 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject().

_cut_elem_nnode

const unsigned int CrackMeshCut3DUserObject::_cut_elem_nnode = 3

protected

The cutter mesh has triangluar elements only.

Definition at line 71 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), findBoundaryEdges(), and triangulation().

_cutter_mesh

std::unique_ptr<MeshBase> MeshCutUserObjectBase::_cutter_mesh

protectedinherited

The xfem cutter mesh.

Definition at line 32 of file MeshCutUserObjectBase.h.

Referenced by MeshCut2DUserObjectBase::addNucleatedCracksToMesh(), CrackMeshCut3DUserObject(), MeshCut2DUserObjectBase::cutElementByGeometry(), MeshCut2DUserObjectBase::cutFragmentByGeometry(), findActiveBoundaryDirection(), MeshCut2DFunctionUserObject::findActiveBoundaryGrowth(), findActiveBoundaryNodes(), findBoundaryEdges(), findBoundaryNodes(), findDistance(), findFrontIntersection(), MeshCut2DUserObjectBase::findOriginalCrackFrontNodes(), getCrackFrontPoints(), MeshCut2DUserObjectBase::getCrackFrontPoints(), getCrackPlaneNormals(), MeshCut2DUserObjectBase::getCrackPlaneNormals(), MeshCutUserObjectBase::getCutterMesh(), MeshCut2DUserObjectBase::growFront(), growFront(), MeshCut2DUserObjectBase::MeshCut2DUserObjectBase(), MeshCutUserObjectBase::MeshCutUserObjectBase(), refineBoundary(), refineFront(), MeshCut2DUserObjectBase::removeNucleatedCracksTooCloseToExistingCracks(), and triangulation().

_elem_dim

const unsigned int CrackMeshCut3DUserObject::_elem_dim = 3

protected

The structural mesh must be 3D only.

Definition at line 106 of file CrackMeshCut3DUserObject.h.

_front

std::vector<std::vector<dof_id_type> > CrackMeshCut3DUserObject::_front

protected

New boundary after growth.

Definition at line 147 of file CrackMeshCut3DUserObject.h.

Referenced by findFrontIntersection(), growFront(), joinBoundary(), refineFront(), and triangulation().

_func_v

const Function* CrackMeshCut3DUserObject::_func_v

protected

Definition at line 260 of file CrackMeshCut3DUserObject.h.

Referenced by growFront().

_func_x

const Function* CrackMeshCut3DUserObject::_func_x

protected

Parsed functions of front growth.

Definition at line 257 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

_func_y

const Function* CrackMeshCut3DUserObject::_func_y

protected

Definition at line 258 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

_func_z

const Function* CrackMeshCut3DUserObject::_func_z

protected

Definition at line 259 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

_grow

bool CrackMeshCut3DUserObject::_grow

protected

Definition at line 119 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and initialize().

_growth_dir_method

const GrowthDirectionEnum CrackMeshCut3DUserObject::_growth_dir_method

protected

The direction method for growing mesh at the front.

Definition at line 94 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and findActiveBoundaryDirection().

_growth_inc_reporter

const std::vector<Real>* const CrackMeshCut3DUserObject::_growth_inc_reporter

protected

Pointer to reporter with growth increment if available.

Definition at line 267 of file CrackMeshCut3DUserObject.h.

Referenced by growFront().

_growth_increment_method

const GrowthRateEnum CrackMeshCut3DUserObject::_growth_increment_method

protected

The growth increment method for growing mesh at the front.

Definition at line 103 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and growFront().

_heal_always

bool GeometricCutUserObject::_heal_always

protectedinherited

Heal the mesh.

Definition at line 211 of file GeometricCutUserObject.h.

_inactive_boundary_pos

std::vector<unsigned int> CrackMeshCut3DUserObject::_inactive_boundary_pos

protected

Inactive boundary.

Definition at line 128 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection(), findActiveBoundaryNodes(), getFrontPointsIndex(), growFront(), initialize(), joinBoundary(), refineFront(), and triangulation().

_interface_id

unsigned int GeometricCutUserObject::_interface_id

protectedinherited

Associated interface id.

Definition at line 208 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::finalize(), and GeometricCutUserObject::GeometricCutUserObject().

_is_mesh_modified

bool CrackMeshCut3DUserObject::_is_mesh_modified

protected

Indicator that shows if the cutting mesh is modified or not in this calculation step.

Definition at line 150 of file CrackMeshCut3DUserObject.h.

Referenced by initialize().

_ki_vpp

const std::vector<Real>* const CrackMeshCut3DUserObject::_ki_vpp

protected

Pointer to fracture integral ki if available.

Definition at line 263 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection().

_kii_vpp

const std::vector<Real>* const CrackMeshCut3DUserObject::_kii_vpp

protected

Pointer to fracture integral kii if available.

Definition at line 265 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryDirection().

_last_step_initialized

int GeometricCutUserObject::_last_step_initialized

protectedinherited

Time step information needed to advance a 3D crack only at the real beginning of a time step.

Definition at line 214 of file GeometricCutUserObject.h.

Referenced by initialize().

_marked_elems_2d

std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo2D> > GeometricCutUserObject::_marked_elems_2d

protectedinherited

Containers with information about all 2D and 3D elements marked for cutting by this object.

Definition at line 217 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::deserialize(), GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::initialize(), GeometricCutUserObject::serialize(), and GeometricCutUserObject::threadJoin().

_marked_elems_3d

std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo3D> > GeometricCutUserObject::_marked_elems_3d

protectedinherited

Definition at line 218 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::deserialize(), GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::initialize(), GeometricCutUserObject::serialize(), and GeometricCutUserObject::threadJoin().

_mesh

MooseMesh& CrackMeshCut3DUserObject::_mesh

protected

The structural mesh.

Definition at line 75 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), and findFrontIntersection().

_n_step_growth

unsigned int CrackMeshCut3DUserObject::_n_step_growth

protected

Number of steps to grow the mesh.

Definition at line 115 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), and initialize().

_num_crack_front_points

unsigned int CrackMeshCut3DUserObject::_num_crack_front_points

protected

Total number of crack front points in the mesh cutter.

Definition at line 153 of file CrackMeshCut3DUserObject.h.

Referenced by getNumberOfCrackFrontPoints(), initializeCrackFrontNodes(), and refineFront().

_size_control

Real CrackMeshCut3DUserObject::_size_control

protected

Used for cutter mesh refinement and front advancement.

Definition at line 112 of file CrackMeshCut3DUserObject.h.

Referenced by CrackMeshCut3DUserObject(), refineBoundary(), and refineFront().

_stop

bool CrackMeshCut3DUserObject::_stop

protected

Variables to help control the work flow.

Definition at line 118 of file CrackMeshCut3DUserObject.h.

Referenced by findActiveBoundaryNodes(), and initialize().

_tracked_crack_front_points

std::vector<dof_id_type> CrackMeshCut3DUserObject::_tracked_crack_front_points

protected

Front nodes that are grown from the crack front definition defined in the input therefore, they are (1) in the same order as defined in the input and (2) the number of nodes does not change.

Definition at line 132 of file CrackMeshCut3DUserObject.h.

Referenced by findFrontIntersection(), growFront(), initializeCrackFrontNodes(), and refineFront().

_uses_mesh

const bool CrackFrontPointsProvider::_uses_mesh

protectedinherited

bool to set if CrackFrontPointsProvider derived objects use a cutter mesh

Definition at line 49 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontPointsProvider::usesMesh().

_xfem

std::shared_ptr<XFEM> GeometricCutUserObject::_xfem

protectedinherited

Pointer to the XFEM controller object.

Definition at line 205 of file GeometricCutUserObject.h.

Referenced by GeometricCutUserObject::execute(), GeometricCutUserObject::finalize(), GeometricCutUserObject::GeometricCutUserObject(), and GeometricCutUserObject::getCutSubdomainID().

---

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

• xfem/include/userobjects/CrackMeshCut3DUserObject.h
• xfem/src/userobjects/CrackMeshCut3DUserObject.C