https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | List of all members
MeshCut2DUserObjectBase Class Referenceabstract

MeshCut2DUserObjectBase: (1) reads in a mesh describing the crack surface, (2) Fills xfem cut element ojbects. More...

#include <MeshCut2DUserObjectBase.h>

Inheritance diagram for MeshCut2DUserObjectBase:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 MeshCut2DUserObjectBase (const InputParameters &parameters)
 
virtual void initialSetup () override final
 
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
 
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< RealVectorValuegetCrackPlaneNormals (unsigned int num_crack_front_points) const override
 get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject CrackFrontDefinition wants the normal so this implementation of getCrackPlaneNormals gives the CrackFrontDefinition a normal for a line element with a tangent direction in the [001] direction. More...
 
MeshBase & getCutterMesh () const
 
virtual void initialize () override
 
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...
 
SubProblemgetSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () 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)
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObjectprimaryThreadCopy ()
 
std::set< UserObjectName > getDependObjects () const
 
virtual bool needThreadedCopy () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParametersparameters () const
 
MooseObjectName uniqueName () 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
 
getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (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
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) 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 ExecFlagEnumgetExecuteOnEnum () 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 UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (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)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (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< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (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)
 
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 VariableValuecoupledValueByName (const std::string &var_name)
 
virtual const ArrayVariableValuecoupledArrayValueByName (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 ADVariableValuegetADDefaultValue (const std::string &var_name) const
 
const ADVectorVariableValuegetADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradientgetADDefaultGradient () const
 
const ADVectorVariableGradientgetADDefaultVectorGradient () const
 
const ADVariableSecondgetADDefaultSecond () const
 
const ADVectorVariableCurlgetADDefaultCurl () 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_typegetElementID (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_typegetElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const
 
virtual const dof_id_typegetElementIDByName (const std::string &id_name) const
 
virtual const dof_id_typegetElementIDNeighborByName (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_typegetAllElemIDs (unsigned int elem_id_index) const
 
std::set< dof_id_typegetElemIDsOnBlocks (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)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (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 hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (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 VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (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)
 
SamplergetSampler (const std::string &name)
 
T & getSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (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 DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (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 sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 

Public Attributes

const ConsoleStream _console
 

Static Public Attributes

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

virtual void findActiveBoundaryGrowth ()=0
 
Find growth direction at each active node More...
 
void findOriginalCrackFrontNodes ()
 Find the original crack front nodes in the cutter mesh and use to populate _original_and_current_front_node_ids. More...
 
void growFront ()
 grow the cutter mesh More...
 
void addNucleatedCracksToMesh ()
 Calls into MeshCutNucleation UO to add cracks. More...
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObject &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
const ReporterNamegetReporterName (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)
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
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 intcoupledIndices (const std::string &var_name) const
 
virtual const VariableValuecoupledValue (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 VariableValuecoupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledValues (const std::string &var_name) const
 
const ADVariableValueadCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue *> adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagValue (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 ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (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 VariableGradientcoupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradientcoupledVectorTagGradient (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 ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (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 VariableValuecoupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValuecoupledVectorTagArrayDofValue (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 VariableValuecoupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledMatrixTagValue (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 VectorVariableValuecoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue *> coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariablewritableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValuewritableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValuecoupledValueOld (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 VariableValuecoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValuecoupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradients (const std::string &var_name) const
 
const ADVariableGradientadCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradientadCoupledGradientDot (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 ADVectorVariableGradientadCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecondadCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecondadCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradientcoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurladCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledDots (const std::string &var_name) const
 
virtual const VariableValuecoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledDots (const std::string &var_name) const
 
const ADVariableValueadCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledArrayDotDu (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 >::typeadCoupledNodalValue (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 VariableValuecoupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValues (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOlder (const std::string &var_name) const
 
virtual const ArrayVariableValuecoupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValueadCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadZeroValue () const
 
const ADVariableGradientadZeroGradient () const
 
const ADVariableSecondadZeroSecond () 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 MooseVariableFieldBasegetFEVar (const std::string &var_name, unsigned int comp) const
 
const MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp) const
 
MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp)
 
const T * getVarHelper (const std::string &var_name, unsigned int comp) const
 
T * getVarHelper (const std::string &var_name, unsigned int comp)
 
MooseVariablegetVar (const std::string &var_name, unsigned int comp)
 
const MooseVariablegetVar (const std::string &var_name, unsigned int comp) const
 
VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp)
 
const VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp) const
 
ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp)
 
const ArrayMooseVariablegetArrayVar (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 VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (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 VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (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

MooseMesh_mesh
 The FE solution mesh. More...
 
std::unique_ptr< MeshBase > _cutter_mesh
 The xfem cutter mesh. More...
 
const MeshCut2DNucleationBase_nucleate_uo
 2D UO for nucleating cracks More...
 
bool _is_mesh_modified
 Indicator that shows if the cutting mesh is modified or not in this calculation step. More...
 
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
 This vector of pairs orders crack tips to make the order used in this class the same as those for the fracture integrals VectorPostprocessors created by CrackFrontDefinition. More...
 
std::vector< std::pair< dof_id_type, Point > > _active_front_node_growth_vectors
 contains the active node ids and their growth vectors More...
 
CrackFrontDefinition_crack_front_definition
 user object for communicating between solid_mechanics interaction integrals and 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
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
const THREAD_ID _tid
 
Assembly_assembly
 
const Moose::CoordinateSystemType_coord_sys
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp_app
 
const std::string _type
 
const std::string _name
 
const InputParameters_pars
 
Factory_factory
 
ActionFactory_action_factory
 
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 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< MooseVariableFV< Real > *> _coupled_standard_fv_moose_vars
 
std::vector< MooseLinearVariableFV< Real > *> _coupled_standard_linear_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
 

Private Member Functions

void removeNucleatedCracksTooCloseToEachOther (std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &nucleated_elems_map, Real nucleationRadius)
 Remove nucleated cracks that are too close too each other. More...
 
void removeNucleatedCracksTooCloseToExistingCracks (std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &nucleated_elems_map, Real nucleationRadius)
 Remove nucleated cracks that are too close to a pre-existing crack in the mesh. More...
 

Detailed Description

MeshCut2DUserObjectBase: (1) reads in a mesh describing the crack surface, (2) Fills xfem cut element ojbects.

Derived classes modify the class and grow the mesh

Definition at line 21 of file MeshCut2DUserObjectBase.h.

Constructor & Destructor Documentation

◆ MeshCut2DUserObjectBase()

MeshCut2DUserObjectBase::MeshCut2DUserObjectBase ( const InputParameters parameters)

Definition at line 36 of file MeshCut2DUserObjectBase.C.

39  _nucleate_uo(isParamValid("nucleate_uo")
40  ? &getUserObject<MeshCut2DNucleationBase>("nucleate_uo")
41  : nullptr),
42  _is_mesh_modified(false)
43 {
44  _depend_uo.insert(getParam<UserObjectName>("crack_front_definition"));
45 
46  // only the Exodus type is currently supported
47  MeshFileName cutterMeshFileName = getParam<MeshFileName>("mesh_file");
48  _cutter_mesh = std::make_unique<ReplicatedMesh>(_communicator);
49  _cutter_mesh->read(cutterMeshFileName);
50  // test element type; only line elements are allowed
51  for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
52  {
53  if (cut_elem->n_nodes() != 2)
54  mooseError("The input cut mesh should include EDGE2 elements only!");
55  if (cut_elem->dim() != 1)
56  mooseError("The input cut mesh should have 1D elements (in a 2D space) only!");
57  }
58 
59  // find node fronts of the original cutmesh. This is used to order EVERYTHING.
61 }
virtual MooseMesh & mesh()=0
void findOriginalCrackFrontNodes()
Find the original crack front nodes in the cutter mesh and use to populate _original_and_current_fron...
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
GeometricCutUserObject(const InputParameters &parameters, const bool uses_mesh=false)
const Parallel::Communicator & _communicator
SubProblem & _subproblem
bool isParamValid(const std::string &name) const
MooseMesh & _mesh
The FE solution mesh.
const MeshCut2DNucleationBase * _nucleate_uo
2D UO for nucleating cracks
std::set< std::string > _depend_uo
bool _is_mesh_modified
Indicator that shows if the cutting mesh is modified or not in this calculation step.
void mooseError(Args &&... args) const
const InputParameters & parameters() const

Member Function Documentation

◆ addNucleatedCracksToMesh()

void MeshCut2DUserObjectBase::addNucleatedCracksToMesh ( )
protected

Calls into MeshCutNucleation UO to add cracks.

Definition at line 363 of file MeshCut2DUserObjectBase.C.

Referenced by MeshCut2DFractureUserObject::initialize(), and MeshCut2DFunctionUserObject::initialize().

364 {
365  if (_nucleate_uo)
366  {
367  std::map<unsigned int, std::pair<RealVectorValue, RealVectorValue>> nucleated_elems_map =
369  const Real nucleationRadius = _nucleate_uo->getNucleationRadius();
370 
371  removeNucleatedCracksTooCloseToEachOther(nucleated_elems_map, nucleationRadius);
372  removeNucleatedCracksTooCloseToExistingCracks(nucleated_elems_map, nucleationRadius);
373 
374  std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
375  pl->enable_out_of_mesh_mode();
376  for (const auto & elem_nodes : nucleated_elems_map)
377  {
378  std::pair<RealVectorValue, RealVectorValue> nodes = elem_nodes.second;
379  // add nodes for the elements that define the nucleated cracks
380  Node * node_0 = Node::build(nodes.first, _cutter_mesh->n_nodes()).release();
381  _cutter_mesh->add_node(node_0);
382  dof_id_type node_id_0 = _cutter_mesh->n_nodes() - 1;
383  Node * node_1 = Node::build(nodes.second, _cutter_mesh->n_nodes()).release();
384  _cutter_mesh->add_node(node_1);
385  dof_id_type node_id_1 = _cutter_mesh->n_nodes() - 1;
386  // add elements that define nucleated cracks
387  std::vector<dof_id_type> elem;
388  elem.push_back(node_id_0);
389  elem.push_back(node_id_1);
390  Elem * new_elem = Elem::build(EDGE2).release();
391  for (unsigned int i = 0; i < new_elem->n_nodes(); ++i)
392  {
393  mooseAssert(_cutter_mesh->node_ptr(elem[i]) != nullptr, "Node is NULL");
394  new_elem->set_node(i, _cutter_mesh->node_ptr(elem[i]));
395  }
396  _cutter_mesh->add_elem(new_elem);
397  // now add the nucleated nodes to the crack id data struct
398  // edge nucleated cracks will add one node to _original_and_current_front_node_ids
399  // bulk nucleated cracks will add two nodes to _original_and_current_front_node_ids
400  Point & point_0 = *node_0;
401  const Elem * crack_front_elem_0 = (*pl)(point_0);
402  if (crack_front_elem_0 != NULL)
403  _original_and_current_front_node_ids.push_back(std::make_pair(node_id_0, node_id_0));
404 
405  Point & point_1 = *node_1;
406  const Elem * crack_front_elem_1 = (*pl)(point_1);
407  if (crack_front_elem_1 != NULL)
408  _original_and_current_front_node_ids.push_back(std::make_pair(node_id_1, node_id_1));
409 
410  _is_mesh_modified = true;
411  }
412  _cutter_mesh->prepare_for_use();
413  }
414 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
This vector of pairs orders crack tips to make the order used in this class the same as those for the...
Real getNucleationRadius() const
Provide getter to MeshCut2DUserObjectBase for member data set in input.
MooseMesh & _mesh
The FE solution mesh.
std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue > > getNucleatedElemsMap() const
Provide getter to MeshCut2DUserObjectBase for a map of nucleated cracks.
const MeshCut2DNucleationBase * _nucleate_uo
2D UO for nucleating cracks
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
bool _is_mesh_modified
Indicator that shows if the cutting mesh is modified or not in this calculation step.
virtual std::unique_ptr< libMesh::PointLocatorBase > getPointLocator() const
void removeNucleatedCracksTooCloseToExistingCracks(std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &nucleated_elems_map, Real nucleationRadius)
Remove nucleated cracks that are too close to a pre-existing crack in the mesh.
void removeNucleatedCracksTooCloseToEachOther(std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &nucleated_elems_map, Real nucleationRadius)
Remove nucleated cracks that are too close too each other.
uint8_t dof_id_type

◆ cutElementByGeometry() [1/4]

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

Definition at line 71 of file MeshCut2DUserObjectBase.C.

74 {
75  // With the crack defined by a line, this method cuts a 2D elements by a line
76  // Fixme lynn Copy and paste from InterfaceMeshCut2DUserObject::cutElementByGeometry
77  mooseAssert(elem->dim() == 2, "Dimension of element to be cut must be 2");
78 
79  bool elem_cut = false;
80 
81  for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
82  {
83  unsigned int n_sides = elem->n_sides();
84 
85  for (unsigned int i = 0; i < n_sides; ++i)
86  {
87  std::unique_ptr<const Elem> curr_side = elem->side_ptr(i);
88 
89  mooseAssert(curr_side->type() == EDGE2, "Element side type must be EDGE2.");
90 
91  const Node * node1 = curr_side->node_ptr(0);
92  const Node * node2 = curr_side->node_ptr(1);
93  Real seg_int_frac = 0.0;
94 
95  const std::pair<Point, Point> elem_endpoints(cut_elem->node_ref(0), cut_elem->node_ref(1));
96 
97  if (Xfem::intersectSegmentWithCutLine(*node1, *node2, elem_endpoints, 1.0, seg_int_frac))
98  {
99  if (seg_int_frac > Xfem::tol && seg_int_frac < 1.0 - Xfem::tol)
100  {
101  elem_cut = true;
102  Xfem::CutEdge mycut;
103  mycut._id1 = node1->id();
104  mycut._id2 = node2->id();
105  mycut._distance = seg_int_frac;
106  mycut._host_side_id = i;
107  cut_edges.push_back(mycut);
108  }
109  else if (seg_int_frac < Xfem::tol)
110  {
111  elem_cut = true;
112  Xfem::CutNode mycut;
113  mycut._id = node1->id();
114  mycut._host_id = i;
115  cut_nodes.push_back(mycut);
116  }
117  }
118  }
119  }
120  return elem_cut;
121 }
Data structure defining a cut through a node.
unsigned int _host_id
Local ID of this node in the host element.
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
Data structure defining a cut on an element edge.
unsigned int _id1
ID of the first node on the edge.
static const double tol
Definition: XFEMFuncs.h:21
bool intersectSegmentWithCutLine(const Point &segment_point1, const Point &segment_point2, const std::pair< Point, Point > &cutting_line_points, const Real &cutting_line_fraction, Real &segment_intersection_fraction)
Determine whether a line segment is intersected by a cutting line, and compute the fraction along tha...
Definition: XFEMFuncs.C:774
Real _distance
Fractional distance along the edge (from node 1 to 2) where the cut is located.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
EDGE2
unsigned int _id2
ID of the second node on the edge.
unsigned int _host_side_id
Local ID of this side in the host element.
unsigned int _id
ID of the cut node.

◆ cutElementByGeometry() [2/4]

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

Definition at line 124 of file MeshCut2DUserObjectBase.C.

126 {
127  mooseError("Invalid method for 2D mesh cutting.");
128  return false;
129 }
void mooseError(Args &&... args) const

◆ 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
elemPointer to the libMesh element to be considered for cutting
cut_edgesData structure filled with information about edges to be cut
cut_nodesData 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
elemPointer to the libMesh element to be considered for cutting
cut_facesData 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 MeshCut2DUserObjectBase::cutFragmentByGeometry ( std::vector< std::vector< Point >> &  frag_edges,
std::vector< Xfem::CutEdge > &  cut_edges 
) const
overridevirtual

Definition at line 132 of file MeshCut2DUserObjectBase.C.

134 {
135  bool cut_frag = false;
136 
137  for (const auto & cut_elem : _cutter_mesh->element_ptr_range())
138  {
139  const std::pair<Point, Point> elem_endpoints(cut_elem->node_ref(0), cut_elem->node_ref(1));
140  unsigned int n_sides = frag_edges.size();
141  for (unsigned int i = 0; i < n_sides; ++i)
142  {
143  Real seg_int_frac = 0.0;
145  frag_edges[i][0], frag_edges[i][1], elem_endpoints, 1, seg_int_frac))
146  {
147  cut_frag = true;
148  Xfem::CutEdge mycut;
149  mycut._id1 = i;
150  mycut._id2 = (i < (n_sides - 1) ? (i + 1) : 0);
151  mycut._distance = seg_int_frac;
152  mycut._host_side_id = i;
153  cut_edges.push_back(mycut);
154  }
155  }
156  }
157  return cut_frag;
158 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
Data structure defining a cut on an element edge.
unsigned int _id1
ID of the first node on the edge.
bool intersectSegmentWithCutLine(const Point &segment_point1, const Point &segment_point2, const std::pair< Point, Point > &cutting_line_points, const Real &cutting_line_fraction, Real &segment_intersection_fraction)
Determine whether a line segment is intersected by a cutting line, and compute the fraction along tha...
Definition: XFEMFuncs.C:774
Real _distance
Fractional distance along the edge (from node 1 to 2) where the cut is located.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
unsigned int _id2
ID of the second node on the edge.
unsigned int _host_side_id
Local ID of this side in the host element.

◆ cutFragmentByGeometry() [2/4]

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

Definition at line 161 of file MeshCut2DUserObjectBase.C.

163 {
164  mooseError("Invalid method for 2D mesh fragment cutting.");
165  return false;
166 }
void mooseError(Args &&... args) const

◆ 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_edgesData structure defining the current fragment to be considered
cut_edgesData 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_facesData structure defining the current fragment to be considered
cut_facesData 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().

220 {
221  mooseAssert(serialized_buffers.size() == _app.n_processors(),
222  "Unexpected size of serialized_buffers: " << serialized_buffers.size());
223 
224  // The input string stream used for deserialization
225  std::istringstream iss;
226 
227  // Loop over all datastructures for all processors to perfrom the gather operation
228  for (unsigned int rank = 0; rank < serialized_buffers.size(); ++rank)
229  {
230  // skip the current processor (its data is already in the structures)
231  if (rank == processor_id())
232  continue;
233 
234  // populate the stream with a new buffer and reset stream state
235  iss.clear();
236  iss.str(serialized_buffers[rank]);
237 
238  // Load the communicated data into temporary structures
239  std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo2D>> other_marked_elems_2d;
240  std::map<unsigned int, std::vector<Xfem::GeomMarkedElemInfo3D>> other_marked_elems_3d;
241  dataLoad(iss, other_marked_elems_2d, this);
242  dataLoad(iss, other_marked_elems_3d, this);
243 
244  // merge the data in with the current processor's data
245  _marked_elems_2d.insert(other_marked_elems_2d.begin(), other_marked_elems_2d.end());
246  _marked_elems_3d.insert(other_marked_elems_3d.begin(), other_marked_elems_3d.end());
247  }
248 }
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...
void dataLoad(std::istream &stream, Xfem::CutFace &cf, void *context)
processor_id_type n_processors() const
processor_id_type processor_id() const
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d

◆ execute()

void GeometricCutUserObject::execute ( )
overridevirtualinherited

Implements ElementUserObject.

Definition at line 64 of file GeometricCutUserObject.C.

65 {
66  if (_current_elem->dim() == 2)
67  {
68  std::vector<Xfem::CutEdge> elem_cut_edges;
69  std::vector<Xfem::CutNode> elem_cut_nodes;
70  std::vector<Xfem::CutEdge> frag_cut_edges;
71  std::vector<std::vector<Point>> frag_edges;
72 
73  EFAElement2D * EFAElem = _xfem->getEFAElem2D(_current_elem);
74 
75  // Don't cut again if elem has been already cut twice
76  if (!EFAElem->isFinalCut())
77  {
78  // get fragment edges
79  _xfem->getFragmentEdges(_current_elem, EFAElem, frag_edges);
80 
81  // mark cut edges for the element and its fragment
82  bool cut = cutElementByGeometry(_current_elem, elem_cut_edges, elem_cut_nodes);
83  if (EFAElem->numFragments() > 0)
84  cut |= cutFragmentByGeometry(frag_edges, frag_cut_edges);
85 
86  if (cut)
87  {
89  gmei2d._elem_cut_edges = elem_cut_edges;
90  gmei2d._elem_cut_nodes = elem_cut_nodes;
91  gmei2d._frag_cut_edges = frag_cut_edges;
92  gmei2d._frag_edges = frag_edges;
93  _marked_elems_2d[_current_elem->id()].push_back(gmei2d);
94  }
95  }
96  }
97  else if (_current_elem->dim() == 3)
98  {
99  std::vector<Xfem::CutFace> elem_cut_faces;
100  std::vector<Xfem::CutFace> frag_cut_faces;
101  std::vector<std::vector<Point>> frag_faces;
102 
103  EFAElement3D * EFAElem = _xfem->getEFAElem3D(_current_elem);
104 
105  // Don't cut again if elem has been already cut twice
106  if (!EFAElem->isFinalCut())
107  {
108  // get fragment edges
109  _xfem->getFragmentFaces(_current_elem, EFAElem, frag_faces);
110 
111  // mark cut faces for the element and its fragment
112  bool cut = cutElementByGeometry(_current_elem, elem_cut_faces);
113  // TODO: This would be done for branching, which is not yet supported in 3D
114  // if (EFAElem->numFragments() > 0)
115  // cut |= cutFragmentByGeometry(frag_faces, frag_cut_faces, _t);
116 
117  if (cut)
118  {
120  gmei3d._elem_cut_faces = elem_cut_faces;
121  gmei3d._frag_cut_faces = frag_cut_faces;
122  gmei3d._frag_faces = frag_faces;
123  _marked_elems_3d[_current_elem->id()].push_back(gmei3d);
124  }
125  }
126  }
127 }
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...
Data structure describing geometrically described cut through 3D element.
virtual unsigned int numFragments() const
Definition: EFAElement2D.C:207
std::vector< CutNode > _elem_cut_nodes
Container for data about all cut nodes in this element.
std::vector< CutEdge > _elem_cut_edges
Container for data about all cut edges in this element.
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.
std::vector< std::vector< Point > > _frag_edges
Container for data about all cut edges in cut fragments in this element.
virtual bool isFinalCut() const
Definition: EFAElement2D.C:796
std::vector< CutEdge > _frag_cut_edges
Container for data about all cut fragments in this element.
Data structure describing geometrically described cut through 2D element.
std::shared_ptr< XFEM > _xfem
Pointer to the XFEM controller object.
const Elem *const & _current_elem
std::vector< CutFace > _frag_cut_faces
Container for data about all faces this element&#39;s fragment.
std::vector< std::vector< Point > > _frag_faces
Container for data about all cut faces in cut fragments in this element.
std::vector< CutFace > _elem_cut_faces
Container for data about all cut faces in this element.
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.
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d
virtual bool isFinalCut() const
Definition: EFAElement3D.C:825

◆ finalize()

void GeometricCutUserObject::finalize ( )
overridevirtualinherited

Implements ElementUserObject.

Definition at line 251 of file GeometricCutUserObject.C.

252 {
253  // for single processor runs we do not need to do anything here
254  if (_app.n_processors() > 1)
255  {
256  // create send buffer
257  std::string send_buffer;
258 
259  // create byte buffers for the streams received from all processors
260  std::vector<std::string> recv_buffers;
261 
262  // pack the complex datastructures into the string stream
263  serialize(send_buffer);
264 
265  // broadcast serialized data to and receive from all processors
266  _communicator.allgather(send_buffer, recv_buffers);
267 
268  // unpack the received data and merge it into the local data structures
269  deserialize(recv_buffers);
270  }
271 
272  for (const auto & it : _marked_elems_2d)
273  for (const auto & gmei : it.second)
274  _xfem->addGeomMarkedElem2D(it.first, gmei, _interface_id);
275 
276  for (const auto & it : _marked_elems_3d)
277  for (const auto & gmei : it.second)
278  _xfem->addGeomMarkedElem3D(it.first, gmei, _interface_id);
279 
280  _marked_elems_2d.clear();
281  _marked_elems_3d.clear();
282 }
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...
void allgather(const T &send_data, std::vector< T, A > &recv_data) const
unsigned int _interface_id
Associated interface id.
const Parallel::Communicator & _communicator
processor_id_type n_processors() const
void deserialize(std::vector< std::string > &serialized_buffers)
void serialize(std::string &serialized_buffer)
Methods to pack/unpack the _marked_elems_2d and _marked_elems_3d data into a structure suitable for p...
std::shared_ptr< XFEM > _xfem
Pointer to the XFEM controller object.
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d

◆ findActiveBoundaryGrowth()

virtual void MeshCut2DUserObjectBase::findActiveBoundaryGrowth ( )
protectedpure virtual


Find growth direction at each active node

Implemented in MeshCut2DFunctionUserObject, and MeshCut2DFractureUserObject.

◆ findOriginalCrackFrontNodes()

void MeshCut2DUserObjectBase::findOriginalCrackFrontNodes ( )
protected

Find the original crack front nodes in the cutter mesh and use to populate _original_and_current_front_node_ids.

Definition at line 284 of file MeshCut2DUserObjectBase.C.

Referenced by MeshCut2DUserObjectBase().

285 {
286  std::unique_ptr<PointLocatorBase> pl = _mesh.getPointLocator();
287  pl->enable_out_of_mesh_mode();
288  std::unordered_set boundary_nodes = MeshTools::find_boundary_nodes(*_cutter_mesh);
289  for (const auto & node : boundary_nodes)
290  {
291  auto node_id = node;
292  Node * this_node = _cutter_mesh->node_ptr(node_id);
293  mooseAssert(this_node, "Node is NULL");
294  Point & this_point = *this_node;
295 
296  const Elem * elem = (*pl)(this_point);
297  if (elem != NULL)
298  _original_and_current_front_node_ids.push_back(std::make_pair(node, node));
299  }
300  std::sort(_original_and_current_front_node_ids.begin(),
302 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
This vector of pairs orders crack tips to make the order used in this class the same as those for the...
MooseMesh & _mesh
The FE solution mesh.
virtual std::unique_ptr< libMesh::PointLocatorBase > getPointLocator() const

◆ getCrackFrontPoints()

const std::vector< Point > MeshCut2DUserObjectBase::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 176 of file MeshCut2DUserObjectBase.C.

177 {
178  std::vector<Point> crack_front_points(number_crack_front_points);
179  // number_crack_front_points is updated via
180  // _crack_front_definition->updateNumberOfCrackFrontPoints(_crack_front_points.size())
181  if (number_crack_front_points != _original_and_current_front_node_ids.size())
182  mooseError("MeshCut2DFractureUserObject::getCrackFrontPoints: number_crack_front_points=" +
183  Moose::stringify(number_crack_front_points) +
184  " does not match the number of nodes given in "
185  "_original_and_current_front_node_ids=" +
187 
188  for (unsigned int i = 0; i < number_crack_front_points; ++i)
189  {
191  Node * this_node = _cutter_mesh->node_ptr(id);
192  mooseAssert(this_node, "Node is NULL");
193  Point & this_point = *this_node;
194  crack_front_points[i] = this_point;
195  }
196  return crack_front_points;
197 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
This vector of pairs orders crack tips to make the order used in this class the same as those for the...
std::string stringify(const T &t)
void mooseError(Args &&... args) const
uint8_t dof_id_type

◆ getCrackPlaneNormals()

const std::vector< RealVectorValue > MeshCut2DUserObjectBase::getCrackPlaneNormals ( unsigned int  num_crack_front_points) const
overridevirtual

get a set of normal vectors along a crack front from a XFEM GeometricCutUserObject CrackFrontDefinition wants the normal so this implementation of getCrackPlaneNormals gives the CrackFrontDefinition a normal for a line element with a tangent direction in the [001] direction.

Returns
A vector which contains all crack front normals

Implements CrackFrontPointsProvider.

Definition at line 200 of file MeshCut2DUserObjectBase.C.

201 {
202  if (number_crack_front_points != _original_and_current_front_node_ids.size())
203  mooseError("MeshCut2DFractureUserObject::getCrackPlaneNormals: number_crack_front_points=" +
204  Moose::stringify(number_crack_front_points) +
205  " does not match the number of nodes given in "
206  "_original_and_current_front_node_ids=" +
208  ". This will happen if a crack front exits the boundary because the number of "
209  "points in the CrackFrontDefinition is never updated.");
210 
211  std::vector<std::pair<dof_id_type, RealVectorValue>> crack_plane_normals;
212  for (const auto & elem : _cutter_mesh->element_ptr_range())
213  {
214  dof_id_type id0 = elem->node_id(0);
215  dof_id_type id1 = elem->node_id(1);
216  dof_id_type id;
217 
218  auto it0 = std::find_if(_original_and_current_front_node_ids.begin(),
220  [&id0](const std::pair<dof_id_type, dof_id_type> & element)
221  { return element.second == id0; });
222  auto it1 = std::find_if(_original_and_current_front_node_ids.begin(),
224  [&id1](const std::pair<dof_id_type, dof_id_type> & element)
225  { return element.second == id1; });
226 
227  bool found_it0 = (it0 != _original_and_current_front_node_ids.end());
228  bool found_it1 = (it1 != _original_and_current_front_node_ids.end());
229 
230  // Newly nucleated crack elements can have one normal if they are on the edge OR
231  // two normals if they are in the bulk.
232  if (found_it0)
233  {
234  Point end_pt, connecting_pt;
235 
236  end_pt = elem->node_ref(0);
237  connecting_pt = elem->node_ref(1);
238  id = it0->first; // sort by original crack front node ids
239 
240  Point fracture_dir = end_pt - connecting_pt;
241  // The crack normal is orthogonal to the crack extension direction (fracture_dir),
242  // and is defined in this implementation as the cross product of the direction of crack
243  // extension with the tangent direction, which is always (0, 0, 1) in 2D.
244  RealVectorValue normal_dir{fracture_dir(1), -fracture_dir(0), 0};
245  normal_dir /= normal_dir.norm();
246  crack_plane_normals.push_back(std::make_pair(id, normal_dir));
247  }
248 
249  if (found_it1)
250  {
251  Point end_pt, connecting_pt;
252 
253  end_pt = elem->node_ref(1);
254  connecting_pt = elem->node_ref(0);
255  id = it1->first; // sort by original crack front node ids
256 
257  Point fracture_dir = end_pt - connecting_pt;
258  // The crack normal is orthogonal to the crack extension direction (fracture_dir),
259  // and is defined in this implementation as the cross product of the direction of crack
260  // extension with the tangent direction, which is always (0, 0, 1) in 2D.
261  RealVectorValue normal_dir{fracture_dir(1), -fracture_dir(0), 0};
262  normal_dir /= normal_dir.norm();
263  crack_plane_normals.push_back(std::make_pair(id, normal_dir));
264  }
265  }
266  mooseAssert(
267  _original_and_current_front_node_ids.size() == crack_plane_normals.size(),
268  "Boundary nodes are attached to more than one element. This should not happen for a 1D "
269  "cutter mesh."
270  "\n Number of _original_and_current_front_node_ids=" +
272  "\n Number of crack_plane_normals=" + Moose::stringify(crack_plane_normals.size()));
273 
274  // the crack_plane_normals are now sorted by the ORIGINAL crack front ids
275  std::sort(crack_plane_normals.begin(), crack_plane_normals.end());
276  std::vector<RealVectorValue> sorted_crack_plane_normals;
277  for (auto & crack : crack_plane_normals)
278  sorted_crack_plane_normals.push_back(crack.second);
279 
280  return sorted_crack_plane_normals;
281 }
auto norm() const -> decltype(std::norm(Real()))
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
This vector of pairs orders crack tips to make the order used in this class the same as those for the...
std::string stringify(const T &t)
void mooseError(Args &&... args) const
uint8_t dof_id_type

◆ 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
nodePointer to the node
Returns
An unsigned int indicating the side

Definition at line 183 of file GeometricCutUserObject.h.

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

184  {
185  mooseError("Objects that inherit from GeometricCutUserObject should override the "
186  "getCutSubdomainID method");
187  return 0;
188  }
void mooseError(Args &&... args) const

◆ getCutSubdomainID() [2/2]

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

Get the CutSubdomainID for the given element.

Parameters
nodePointer to the element
Returns
The CutSubdomainID

Definition at line 285 of file GeometricCutUserObject.C.

286 {
287  return _xfem->getCutSubdomainID(this, elem);
288 }
std::shared_ptr< XFEM > _xfem
Pointer to the XFEM controller object.

◆ getCutterMesh()

MeshBase & MeshCut2DUserObjectBase::getCutterMesh ( ) const

Definition at line 169 of file MeshCut2DUserObjectBase.C.

Referenced by XFEMCutMeshOutput::output().

170 {
171  mooseAssert(_cutter_mesh, "MeshCut2DUserObjectBase::getCutterMesh _cutter_mesh is nullptr");
172  return *_cutter_mesh;
173 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.

◆ getInterfaceID()

unsigned int GeometricCutUserObject::getInterfaceID ( ) const
inlineinherited

Get the interface ID for this cutting object.

Returns
the interface ID

Definition at line 161 of file GeometricCutUserObject.h.

161 { return _interface_id; };
unsigned int _interface_id
Associated interface id.

◆ growFront()

void MeshCut2DUserObjectBase::growFront ( )
protected

grow the cutter mesh

Definition at line 305 of file MeshCut2DUserObjectBase.C.

Referenced by MeshCut2DFractureUserObject::initialize(), and MeshCut2DFunctionUserObject::initialize().

306 {
307  dof_id_type current_front_node_id;
308  for (std::size_t i = 0; i < _original_and_current_front_node_ids.size(); ++i)
309  {
310  current_front_node_id = _original_and_current_front_node_ids[i].second;
311  // check if node front node id is active
312  auto direction_iter =
313  std::find_if(_active_front_node_growth_vectors.begin(),
315  [&current_front_node_id](const std::pair<dof_id_type, Point> & element)
316  { return element.first == current_front_node_id; });
317  // only add an element for active node front ids
318  if (direction_iter != _active_front_node_growth_vectors.end())
319  {
320  Node * this_node = _cutter_mesh->node_ptr(current_front_node_id);
321  mooseAssert(this_node, "Node is NULL");
322  Point & this_point = *this_node;
323 
324  Point new_node_offset = direction_iter->second;
325  Point x = this_point + new_node_offset;
326 
327  // TODO: Should check if cut line segment created between "this_point" and "x" crosses
328  // another line element in the cutter mesh or solid mesh boundary.
329  // Crossing another line element would be a special case that still needs to be handled,
330  // however, it doesnot cause an error, it will just ignore the other line segment and recut
331  // the solid mesh element.
332  // Crossing a solid mesh boundary would be for aesthetics reasons so
333  // that element was trimmed close to the boundary but would have not effect on the simulation.
334  // Crossing a solid mesh boundary should be handled by something like
335  // MeshCut2DRankTwoTensorNucleation::lineLineIntersect2D
336 
337  // add node to front
338  this_node = Node::build(x, _cutter_mesh->n_nodes()).release();
339  _cutter_mesh->add_node(this_node);
340  dof_id_type new_front_node_id = _cutter_mesh->n_nodes() - 1;
341 
342  // add element to front
343  std::vector<dof_id_type> elem;
344  elem.push_back(current_front_node_id);
345  elem.push_back(new_front_node_id);
346  Elem * new_elem = Elem::build(EDGE2).release();
347  for (unsigned int i = 0; i < new_elem->n_nodes(); ++i)
348  {
349  mooseAssert(_cutter_mesh->node_ptr(elem[i]) != nullptr, "Node is NULL");
350  new_elem->set_node(i, _cutter_mesh->node_ptr(elem[i]));
351  }
352  _cutter_mesh->add_elem(new_elem);
353  // now push to the end of _original_and_current_front_node_ids for tracking and fracture
354  // integrals
355  _original_and_current_front_node_ids[i].second = new_front_node_id;
356  _is_mesh_modified = true;
357  }
358  }
359  _cutter_mesh->prepare_for_use();
360 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
std::vector< std::pair< dof_id_type, dof_id_type > > _original_and_current_front_node_ids
This vector of pairs orders crack tips to make the order used in this class the same as those for the...
const std::vector< double > x
bool _is_mesh_modified
Indicator that shows if the cutting mesh is modified or not in this calculation step.
std::vector< std::pair< dof_id_type, Point > > _active_front_node_growth_vectors
contains the active node ids and their growth vectors
uint8_t dof_id_type

◆ initialize()

void GeometricCutUserObject::initialize ( )
overridevirtualinherited

Implements ElementUserObject.

Reimplemented in InterfaceMeshCutUserObjectBase, CrackMeshCut3DUserObject, MeshCut2DFunctionUserObject, and MeshCut2DFractureUserObject.

Definition at line 57 of file GeometricCutUserObject.C.

58 {
59  _marked_elems_2d.clear();
60  _marked_elems_3d.clear();
61 }
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...
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d

◆ initialSetup()

void MeshCut2DUserObjectBase::initialSetup ( )
finaloverridevirtual

Reimplemented from ElementUserObject.

Definition at line 64 of file MeshCut2DUserObjectBase.C.

65 {
66  const auto uo_name = getParam<UserObjectName>("crack_front_definition");
68 }
T & getUserObject(const std::string &name, unsigned int tid=0) const
Class used in fracture integrals to define geometric characteristics of the crack front...
FEProblemBase & _fe_problem
CrackFrontDefinition * _crack_front_definition
user object for communicating between solid_mechanics interaction integrals and xfem cutter mesh ...

◆ removeNucleatedCracksTooCloseToEachOther()

void MeshCut2DUserObjectBase::removeNucleatedCracksTooCloseToEachOther ( std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &  nucleated_elems_map,
Real  nucleationRadius 
)
private

Remove nucleated cracks that are too close too each other.

Lowest map key wins

Parameters
nucleated_elems_mapmap from nucleation userObject with key for mesh element id and two nodes of nucleated crack
nucleationRadiusexclusion distance between cracks

Definition at line 417 of file MeshCut2DUserObjectBase.C.

Referenced by addNucleatedCracksToMesh().

420 {
421  // remove nucleated elements that are too close too each other. Lowest key wins
422  for (auto it1 = nucleated_elems_map.begin(); it1 != nucleated_elems_map.end(); ++it1)
423  {
424  std::pair<RealVectorValue, RealVectorValue> nodes = it1->second;
425  Point p2 = nodes.first;
426  Point p1 = nodes.second;
427  Point p = p1 + (p2 - p1) / 2;
428  for (auto it2 = nucleated_elems_map.begin(); it2 != nucleated_elems_map.end();)
429  {
430  if (it1 == it2)
431  {
432  ++it2;
433  continue;
434  }
435 
436  nodes = it2->second;
437  p2 = nodes.first;
438  p1 = nodes.second;
439  Point q = p1 + (p2 - p1) / 2;
440  Point pq = q - p;
441  if (pq.norm() <= nucleationRadius)
442  it2 = nucleated_elems_map.erase(it2);
443  else
444  ++it2;
445  }
446  }
447 }

◆ removeNucleatedCracksTooCloseToExistingCracks()

void MeshCut2DUserObjectBase::removeNucleatedCracksTooCloseToExistingCracks ( std::map< unsigned int, std::pair< RealVectorValue, RealVectorValue >> &  nucleated_elems_map,
Real  nucleationRadius 
)
private

Remove nucleated cracks that are too close to a pre-existing crack in the mesh.

Parameters
nucleated_elems_mapmap from nucleation userObject with key for mesh element id and two nodes of nucleated crack
nucleationRadiusexclusion distance between cracks

Definition at line 450 of file MeshCut2DUserObjectBase.C.

Referenced by addNucleatedCracksToMesh().

453 {
454  for (auto it = nucleated_elems_map.begin(); it != nucleated_elems_map.end();)
455  {
456  std::pair<RealVectorValue, RealVectorValue> nodes = it->second;
457  Point p2 = nodes.first;
458  Point p1 = nodes.second;
459  Point p = p1 + (p2 - p1) / 2;
460  bool removeNucleatedElem = false;
461  for (const auto & cutter_elem :
462  as_range(_cutter_mesh->active_elements_begin(), _cutter_mesh->active_elements_end()))
463  {
464  const Node * const * cutter_elem_nodes = cutter_elem->get_nodes();
465  Point m = *cutter_elem_nodes[1] - *cutter_elem_nodes[0];
466  Real t = m * (p - *cutter_elem_nodes[0]) / m.norm_sq();
467  Real d = std::numeric_limits<Real>::max();
468  if (t <= 0)
469  {
470  Point j = p - *cutter_elem_nodes[0];
471  d = j.norm();
472  }
473  else if (t >= 0)
474  {
475  Point j = p - *cutter_elem_nodes[1];
476  d = j.norm();
477  }
478  else
479  {
480  Point j = p - (*cutter_elem_nodes[0] + t * m);
481  d = j.norm();
482  }
483  if (d <= nucleationRadius)
484  {
485  removeNucleatedElem = true;
486  break;
487  }
488  }
489  if (removeNucleatedElem)
490  it = nucleated_elems_map.erase(it);
491  else
492  ++it;
493  }
494 }
std::unique_ptr< MeshBase > _cutter_mesh
The xfem cutter mesh.
SimpleRange< IndexType > as_range(const std::pair< IndexType, IndexType > &p)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")

◆ 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().

207 {
208  // stream for serializing the _marked_elems_2d and _marked_elems_3d data structures to a byte
209  // stream
210  std::ostringstream oss;
211  dataStore(oss, _marked_elems_2d, this);
212  dataStore(oss, _marked_elems_3d, this);
213 
214  // Populate the passed in string pointer with the string stream's buffer contents
215  serialized_buffer.assign(oss.str());
216 }
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...
void dataStore(std::ostream &stream, Xfem::CutFace &cf, void *context)
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d

◆ setInterfaceID()

void GeometricCutUserObject::setInterfaceID ( unsigned int  interface_id)
inlineinherited

Set the interface ID for this cutting object.

Parameters
theinterface ID

Definition at line 167 of file GeometricCutUserObject.h.

Referenced by XFEM::addGeometricCut().

167 { _interface_id = interface_id; };
unsigned int _interface_id
Associated 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 174 of file GeometricCutUserObject.h.

Referenced by XFEM::cutMeshWithEFA().

174 { return _heal_always; };
bool _heal_always
Heal the mesh.

◆ threadJoin()

void GeometricCutUserObject::threadJoin ( const UserObject y)
overridevirtualinherited

Implements ElementUserObject.

Definition at line 130 of file GeometricCutUserObject.C.

131 {
132  const auto & gcuo = static_cast<const GeometricCutUserObject &>(y);
133 
134  for (const auto & it : gcuo._marked_elems_2d)
135  {
136  mooseAssert(_marked_elems_2d.find(it.first) == _marked_elems_2d.end(),
137  "Element already inserted in map from a different thread");
138  _marked_elems_2d[it.first] = it.second;
139  }
140  for (const auto & it : gcuo._marked_elems_3d)
141  {
142  mooseAssert(_marked_elems_3d.find(it.first) == _marked_elems_3d.end(),
143  "Element already inserted in map from a different thread");
144  _marked_elems_3d[it.first] = it.second;
145  }
146 }
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...
const std::vector< double > y
std::map< unsigned int, std::vector< Xfem::GeomMarkedElemInfo3D > > _marked_elems_3d

◆ 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 40 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontDefinition::initialSetup().

40 { return _uses_mesh; }
const bool _uses_mesh
bool to set if CrackFrontPointsProvider derived objects use a cutter mesh

◆ validParams()

InputParameters MeshCut2DUserObjectBase::validParams ( )
static

Definition at line 22 of file MeshCut2DUserObjectBase.C.

Referenced by MeshCut2DFractureUserObject::validParams(), and MeshCut2DFunctionUserObject::validParams().

23 {
25  params.addRequiredParam<MeshFileName>(
26  "mesh_file",
27  "Mesh file for the XFEM geometric cut; currently only the Exodus type is supported");
28  params.addParam<UserObjectName>("nucleate_uo", "The MeshCutNucleation UO for nucleating cracks.");
29  params.addParam<UserObjectName>("crack_front_definition",
30  "crackFrontDefinition",
31  "The CrackFrontDefinition user object name");
32  params.addClassDescription("Creates a UserObject base class for a mesh cutter in 2D problems");
33  return params;
34 }
static InputParameters validParams()
Factory constructor, takes parameters so that all derived classes can be built using the same constru...
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addRequiredParam(const std::string &name, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _active_front_node_growth_vectors

std::vector<std::pair<dof_id_type, Point> > MeshCut2DUserObjectBase::_active_front_node_growth_vectors
protected

contains the active node ids and their growth vectors

Definition at line 78 of file MeshCut2DUserObjectBase.h.

Referenced by MeshCut2DFractureUserObject::findActiveBoundaryGrowth(), MeshCut2DFunctionUserObject::findActiveBoundaryGrowth(), and growFront().

◆ _crack_front_definition

CrackFrontDefinition* MeshCut2DUserObjectBase::_crack_front_definition
protected

user object for communicating between solid_mechanics interaction integrals and xfem cutter mesh

Definition at line 81 of file MeshCut2DUserObjectBase.h.

Referenced by MeshCut2DFractureUserObject::findActiveBoundaryGrowth(), MeshCut2DFractureUserObject::initialize(), MeshCut2DFunctionUserObject::initialize(), and initialSetup().

◆ _cutter_mesh

std::unique_ptr<MeshBase> MeshCut2DUserObjectBase::_cutter_mesh
protected

◆ _heal_always

bool GeometricCutUserObject::_heal_always
protectedinherited

Heal the mesh.

Definition at line 205 of file GeometricCutUserObject.h.

◆ _interface_id

unsigned int GeometricCutUserObject::_interface_id
protectedinherited

Associated interface id.

Definition at line 202 of file GeometricCutUserObject.h.

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

◆ _is_mesh_modified

bool MeshCut2DUserObjectBase::_is_mesh_modified
protected

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

Definition at line 65 of file MeshCut2DUserObjectBase.h.

Referenced by addNucleatedCracksToMesh(), growFront(), MeshCut2DFractureUserObject::initialize(), and MeshCut2DFunctionUserObject::initialize().

◆ _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 208 of file GeometricCutUserObject.h.

Referenced by CrackMeshCut3DUserObject::initialize().

◆ _marked_elems_2d

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

◆ _marked_elems_3d

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

◆ _mesh

MooseMesh& MeshCut2DUserObjectBase::_mesh
protected

The FE solution mesh.

Definition at line 56 of file MeshCut2DUserObjectBase.h.

Referenced by addNucleatedCracksToMesh(), and findOriginalCrackFrontNodes().

◆ _nucleate_uo

const MeshCut2DNucleationBase* MeshCut2DUserObjectBase::_nucleate_uo
protected

2D UO for nucleating cracks

Definition at line 62 of file MeshCut2DUserObjectBase.h.

Referenced by addNucleatedCracksToMesh().

◆ _original_and_current_front_node_ids

std::vector<std::pair<dof_id_type, dof_id_type> > MeshCut2DUserObjectBase::_original_and_current_front_node_ids
protected

This vector of pairs orders crack tips to make the order used in this class the same as those for the fracture integrals VectorPostprocessors created by CrackFrontDefinition.

The original crack front node ids found in the cutter mesh are put in pair.first and the assosciated current crack front node id that grew from the original crack front node id is in pair.second. This vector is sorted on pair.first which makes the ordering of this vector the same as that used in the CrackFrontDefinition

Definition at line 75 of file MeshCut2DUserObjectBase.h.

Referenced by addNucleatedCracksToMesh(), MeshCut2DFractureUserObject::findActiveBoundaryGrowth(), MeshCut2DFunctionUserObject::findActiveBoundaryGrowth(), findOriginalCrackFrontNodes(), getCrackFrontPoints(), getCrackPlaneNormals(), growFront(), and MeshCut2DFractureUserObject::initialize().

◆ _uses_mesh

const bool CrackFrontPointsProvider::_uses_mesh
protectedinherited

bool to set if CrackFrontPointsProvider derived objects use a cutter mesh

Definition at line 44 of file CrackFrontPointsProvider.h.

Referenced by CrackFrontPointsProvider::usesMesh().

◆ _xfem

std::shared_ptr<XFEM> GeometricCutUserObject::_xfem
protectedinherited

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