Mesh class for triangular, edge and corner inter_wrappers for hexagonal lattice fuel assemblies. More...

#include <TriInterWrapperMesh.h>

Inheritance diagram for TriInterWrapperMesh:

Public Types
enum	ParallelType { ParallelType::DEFAULT, ParallelType::REPLICATED, ParallelType::DISTRIBUTED }

typedef std::pair< const Node *, BoundaryID >	PeriodicNodeInfo

typedef DataFileName	DataFileParameterType

Public Member Functions
	TriInterWrapperMesh (const InputParameters &parameters)

	TriInterWrapperMesh (const TriInterWrapperMesh &other_mesh)

virtual std::unique_ptr< MooseMesh >	safeClone () const override

virtual void	buildMesh () override

virtual const unsigned int &	getNumOfAssemblies () const override
	Return the number of assemblies. More...

virtual Node *	getPinNode (unsigned int, unsigned) const override
	Get the pin mesh node for a given pin index and elevation index. More...

virtual bool	pinMeshExist () const override
	Return if Pin Mesh exists or not. More...

virtual bool	ductMeshExist () const override
	Return if duct Mesh exists or not. More...

virtual const Real &	getDuctToPinGap () const

virtual const unsigned int &	getPinIndex (const unsigned int channel_idx, const unsigned int neighbor_idx)
	Return Pin index given inter_wrapper index and local neighbor index. More...

virtual Node *	getChannelNode (unsigned int i_chan, unsigned iz) const override
	Get the inter-wrapper mesh node for a given channel index and elevation index. More...

virtual const unsigned int &	getNumOfChannels () const override
	Return the number of channels per layer. More...

virtual const unsigned int &	getNumOfGapsPerLayer () const override
	Return the number of gaps per layer. More...

virtual const std::pair< unsigned int, unsigned int > &	getGapChannels (unsigned int i_gap) const override
	Return a pair of inter-wrapper indices for a given gap index. More...

virtual const std::vector< unsigned int > &	getChannelGaps (unsigned int i_chan) const override
	Return a vector of gap indices for a given channel index. More...

virtual const std::vector< double > &	getGapMap () const override
	Return a map with gap sizes. More...

virtual const Real &	getCrossflowSign (unsigned int i_chan, unsigned int i_local) const override
	Return a signs for the cross flow given a inter-wrapper index and local neighbor index. More...

virtual unsigned int	getSubchannelIndexFromPoint (const Point &p) const override
	Return a inter-wrapper index for a given physical point `p` More...

virtual unsigned int	channelIndex (const Point &point) const override

virtual EChannelType	getSubchannelType (unsigned int index) const override
	Return the type of the inter-wrapper for given inter-wrapper index. More...

virtual Real	getGapWidth (unsigned int gap_index) const override
	Return gap width for a given gap index. More...

virtual const std::pair< unsigned int, unsigned int > &	getSweepFlowGaps (unsigned int i_chan) const

virtual const std::pair< unsigned int, unsigned int > &	getSweepFlowChans (unsigned int i_chan) const

virtual const std::vector< unsigned int > &	getPinChannels (unsigned int i_pin) const override
	Return a vector of channel indices for a given Pin index. More...

virtual const std::vector< unsigned int > &	getChannelPins (unsigned int i_chan) const override
	Return a vector of pin indices for a given channel index. More...

virtual const bool &	getIsTightSide () const

virtual unsigned int	getPinIndexFromPoint (const Point &p) const override
	Return a pin index for a given physical point `p` More...

virtual unsigned int	pinIndex (const Point &p) const override

void	setChannelToDuctMaps (const std::vector< Node *> &duct_nodes)
	Setup the internal maps when there is a outside duct present. More...

virtual const std::vector< Real > &	getZGrid () const
	Get axial location of layers. More...

virtual unsigned int	getZIndex (const Point &point) const
	Get axial index of point. More...

virtual const std::vector< std::vector< Real > > &	getKGrid () const
	Get axial cell location and value of loss coefficient. More...

virtual const Real &	getKij () const
	Return lateral loss coefficient. More...

virtual const unsigned int &	getNumOfAxialCells () const
	Return the number of axial cells. More...

virtual const Real &	getPitch () const
	Return the pitch between 2 inter-wrappers. More...

virtual const Real &	getSideX () const
	Return side lengths of the assembly. More...

virtual const Real &	getSideY () const

virtual const Real &	getHeatedLengthEntry () const
	Return unheated length at entry. More...

virtual const Real &	getHeatedLength () const
	Return heated length. More...

virtual const Real &	getHeatedLengthExit () const
	Return unheated length at exit. More...

virtual MooseMesh &	clone () const

void	determineUseDistributedMesh ()

std::unique_ptr< MeshBase >	buildMeshBaseObject (unsigned int dim=libMesh::invalid_uint)

std::unique_ptr< T >	buildTypedMesh (unsigned int dim=libMesh::invalid_uint)

void	setMeshBase (std::unique_ptr< MeshBase > mesh_base)

virtual void	init ()

virtual unsigned int	dimension () const

virtual unsigned int	spatialDimension () const

virtual unsigned int	effectiveSpatialDimension () const

unsigned int	getBlocksMaxDimension (const std::vector< SubdomainName > &blocks) const

std::vector< BoundaryID >	getBoundaryIDs (const Elem *const elem, const unsigned short int side) const

const std::set< BoundaryID > &	getBoundaryIDs () const

std::vector< BoundaryID >	getBoundaryIDs (const std::vector< BoundaryName > &boundary_name, bool generate_unknown=false) const

const Elem *	getLowerDElem (const Elem *, unsigned short int) const

unsigned int	getHigherDSide (const Elem *elem) const

void	buildNodeList ()

void	buildBndElemList ()

const std::map< dof_id_type, std::vector< dof_id_type > > &	nodeToElemMap ()

const std::map< dof_id_type, std::vector< dof_id_type > > &	nodeToActiveSemilocalElemMap ()

virtual bnd_node_iterator	bndNodesBegin ()

virtual bnd_node_iterator	bndNodesEnd ()

virtual bnd_elem_iterator	bndElemsBegin ()

virtual bnd_elem_iterator	bndElemsEnd ()

void	buildNodeListFromSideList ()

void	buildSideList (std::vector< dof_id_type > &el, std::vector< unsigned short int > &sl, std::vector< boundary_id_type > &il)

std::vector< std::tuple< dof_id_type, unsigned short int, boundary_id_type > >	buildSideList ()

std::vector< std::tuple< dof_id_type, unsigned short int, boundary_id_type > >	buildActiveSideList () const

unsigned int	sideWithBoundaryID (const Elem *const elem, const BoundaryID boundary_id) const

MeshBase::node_iterator	localNodesBegin ()

MeshBase::const_node_iterator	localNodesBegin () const

MeshBase::node_iterator	localNodesEnd ()

MeshBase::const_node_iterator	localNodesEnd () const

MeshBase::element_iterator	activeLocalElementsBegin ()

MeshBase::const_element_iterator	activeLocalElementsBegin () const

const MeshBase::element_iterator	activeLocalElementsEnd ()

const MeshBase::const_element_iterator	activeLocalElementsEnd () const

virtual dof_id_type	nNodes () const

virtual dof_id_type	nElem () const

virtual dof_id_type	nLocalNodes () const

virtual dof_id_type	nActiveElem () const

virtual dof_id_type	nActiveLocalElem () const

virtual SubdomainID	nSubdomains () const

virtual unsigned int	nPartitions () const

virtual bool	skipPartitioning () const

virtual bool	skipNoncriticalPartitioning () const

virtual dof_id_type	maxNodeId () const

virtual dof_id_type	maxElemId () const

virtual const Node &	node (const dof_id_type i) const

virtual Node &	node (const dof_id_type i)

virtual const Node &	nodeRef (const dof_id_type i) const

virtual Node &	nodeRef (const dof_id_type i)

virtual const Node *	nodePtr (const dof_id_type i) const

virtual Node *	nodePtr (const dof_id_type i)

virtual const Node *	queryNodePtr (const dof_id_type i) const

virtual Node *	queryNodePtr (const dof_id_type i)

virtual Elem *	elem (const dof_id_type i)

virtual const Elem *	elem (const dof_id_type i) const

virtual Elem *	elemPtr (const dof_id_type i)

virtual const Elem *	elemPtr (const dof_id_type i) const

virtual Elem *	queryElemPtr (const dof_id_type i)

virtual const Elem *	queryElemPtr (const dof_id_type i) const

bool	prepared () const

virtual void	prepared (bool state)

void	needsPrepareForUse ()

void	meshChanged ()

virtual void	onMeshChanged ()

void	cacheChangedLists ()

ConstElemPointerRange *	refinedElementRange () const

ConstElemPointerRange *	coarsenedElementRange () const

const std::vector< const Elem *> &	coarsenedElementChildren (const Elem *elem) const

void	updateActiveSemiLocalNodeRange (std::set< dof_id_type > &ghosted_elems)

bool	isSemiLocal (Node *const node) const

const std::unordered_map< boundary_id_type, std::unordered_set< dof_id_type > > &	getBoundariesToElems () const

const std::unordered_map< boundary_id_type, std::unordered_set< dof_id_type > > &	getBoundariesToActiveSemiLocalElemIds () const

std::unordered_set< dof_id_type >	getBoundaryActiveSemiLocalElemIds (BoundaryID bid) const

std::unordered_set< dof_id_type >	getBoundaryActiveNeighborElemIds (BoundaryID bid) const

bool	isBoundaryFullyExternalToSubdomains (BoundaryID bid, const std::set< SubdomainID > &blk_group) const

const std::set< SubdomainID > &	meshSubdomains () const

const std::set< BoundaryID > &	meshBoundaryIds () const

const std::set< BoundaryID > &	meshSidesetIds () const

const std::set< BoundaryID > &	meshNodesetIds () const

void	setBoundaryToNormalMap (std::unique_ptr< std::map< BoundaryID, RealVectorValue >> boundary_map)

void	setBoundaryToNormalMap (std::map< BoundaryID, RealVectorValue > *boundary_map)

void	setMeshBoundaryIDs (std::set< BoundaryID > boundary_IDs)

const RealVectorValue &	getNormalByBoundaryID (BoundaryID id) const

bool	prepare (const MeshBase *mesh_to_clone)

void	update ()

unsigned int	uniformRefineLevel () const

void	setUniformRefineLevel (unsigned int, bool deletion=true)

bool	skipDeletionRepartitionAfterRefine () const

bool	skipRefineWhenUseSplit () const

void	addGhostedBoundary (BoundaryID boundary_id)

void	setGhostedBoundaryInflation (const std::vector< Real > &inflation)

const std::set< unsigned int > &	getGhostedBoundaries () const

const std::vector< Real > &	getGhostedBoundaryInflation () const

void	ghostGhostedBoundaries ()

void	needGhostGhostedBoundaries (bool needghost)

unsigned int	getPatchSize () const

unsigned int	getGhostingPatchSize () const

unsigned int	getMaxLeafSize () const

void	setPatchUpdateStrategy (Moose::PatchUpdateType patch_update_strategy)

const Moose::PatchUpdateType &	getPatchUpdateStrategy () const

libMesh::BoundingBox	getInflatedProcessorBoundingBox (Real inflation_multiplier=0.01) const

	operator libMesh::MeshBase & ()

	operator const libMesh::MeshBase & () const

MeshBase &	getMesh ()

MeshBase &	getMesh (const std::string &name)

const MeshBase &	getMesh () const

const MeshBase &	getMesh (const std::string &name) const

const MeshBase *	getMeshPtr () const

void	printInfo (std::ostream &os=libMesh::out, const unsigned int verbosity=0) const

const std::set< SubdomainID > &	getNodeBlockIds (const Node &node) const

const std::vector< dof_id_type > &	getNodeList (boundary_id_type nodeset_id) const

const Node *	addUniqueNode (const Point &p, Real tol=1e-6)

Node *	addQuadratureNode (const Elem *elem, const unsigned short int side, const unsigned int qp, BoundaryID bid, const Point &point)

Node *	getQuadratureNode (const Elem *elem, const unsigned short int side, const unsigned int qp)

void	clearQuadratureNodes ()

BoundaryID	getBoundaryID (const BoundaryName &boundary_name) const

SubdomainID	getSubdomainID (const SubdomainName &subdomain_name) const

std::vector< SubdomainID >	getSubdomainIDs (const std::vector< SubdomainName > &subdomain_names) const

std::set< SubdomainID >	getSubdomainIDs (const std::set< SubdomainName > &subdomain_names) const

void	setSubdomainName (SubdomainID subdomain_id, const SubdomainName &name)

const std::string &	getSubdomainName (SubdomainID subdomain_id) const

std::vector< SubdomainName >	getSubdomainNames (const std::vector< SubdomainID > &subdomain_ids) const

void	setBoundaryName (BoundaryID boundary_id, BoundaryName name)

const std::string &	getBoundaryName (BoundaryID boundary_id)

void	buildPeriodicNodeMap (std::multimap< dof_id_type, dof_id_type > &periodic_node_map, unsigned int var_number, libMesh::PeriodicBoundaries *pbs) const

void	buildPeriodicNodeSets (std::map< BoundaryID, std::set< dof_id_type >> &periodic_node_sets, unsigned int var_number, libMesh::PeriodicBoundaries *pbs) const

Real	dimensionWidth (unsigned int component) const

bool	detectOrthogonalDimRanges (Real tol=1e-6)

void	addPeriodicVariable (unsigned int var_num, BoundaryID primary, BoundaryID secondary)

bool	isTranslatedPeriodic (unsigned int nonlinear_var_num, unsigned int component) const

RealVectorValue	minPeriodicVector (unsigned int nonlinear_var_num, Point p, Point q) const

Real	minPeriodicDistance (unsigned int nonlinear_var_num, Point p, Point q) const

const std::pair< BoundaryID, BoundaryID > *	getPairedBoundaryMapping (unsigned int component)

void	buildRefinementAndCoarseningMaps (Assembly *assembly)

const std::vector< std::vector< QpMap > > &	getRefinementMap (const Elem &elem, int parent_side, int child, int child_side)

const std::vector< std::pair< unsigned int, QpMap > > &	getCoarseningMap (const Elem &elem, int input_side)

void	changeBoundaryId (const boundary_id_type old_id, const boundary_id_type new_id, bool delete_prev)

const std::set< BoundaryID > &	getSubdomainBoundaryIds (const SubdomainID subdomain_id) const

std::set< BoundaryID >	getSubdomainInterfaceBoundaryIds (const SubdomainID subdomain_id) const

std::set< SubdomainID >	getBoundaryConnectedBlocks (const BoundaryID bid) const

std::set< SubdomainID >	getBoundaryConnectedSecondaryBlocks (const BoundaryID bid) const

std::set< SubdomainID >	getInterfaceConnectedBlocks (const BoundaryID bid) const

const std::set< SubdomainID > &	getBlockConnectedBlocks (const SubdomainID subdomain_id) const

bool	isBoundaryNode (dof_id_type node_id) const

bool	isBoundaryNode (dof_id_type node_id, BoundaryID bnd_id) const

bool	isBoundaryElem (dof_id_type elem_id) const

bool	isBoundaryElem (dof_id_type elem_id, BoundaryID bnd_id) const

void	errorIfDistributedMesh (std::string name) const

virtual bool	isDistributedMesh () const

bool	isParallelTypeForced () const

void	setParallelType (ParallelType parallel_type)

ParallelType	getParallelType () const

const MooseEnum &	partitionerName () const

bool	isPartitionerForced () const

void	allowRecovery (bool allow)

void	setCustomPartitioner (libMesh::Partitioner *partitioner)

bool	isRegularOrthogonal ()

bool	hasSecondOrderElements ()

virtual std::unique_ptr< libMesh::PointLocatorBase >	getPointLocator () const

virtual std::string	getFileName () const

void	needsRemoteElemDeletion (bool need_delete)

bool	needsRemoteElemDeletion () const

void	allowRemoteElementRemoval (bool allow_removal)

bool	allowRemoteElementRemoval () const

void	deleteRemoteElements ()

bool	hasMeshBase () const

bool	hasElementID (const std::string &id_name) const

unsigned int	getElementIDIndex (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::set< dof_id_type >	getAllElemIDs (unsigned int elem_id_index) const

std::set< dof_id_type >	getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const

std::unordered_map< dof_id_type, std::set< dof_id_type > >	getElemIDMapping (const std::string &from_id_name, const std::string &to_id_name) const

void	cacheFaceInfoVariableOwnership () const

void	cacheFVElementalDoFs () const

void	computeFiniteVolumeCoords () const

void	isDisplaced (bool is_displaced)

bool	isDisplaced () const

const std::map< boundary_id_type, std::vector< dof_id_type > > &	nodeSetNodes () const

Moose::CoordinateSystemType	getCoordSystem (SubdomainID sid) const

const std::map< SubdomainID, Moose::CoordinateSystemType > &	getCoordSystem () const

Moose::CoordinateSystemType	getUniqueCoordSystem () const

void	setCoordSystem (const std::vector< SubdomainName > &blocks, const MultiMooseEnum &coord_sys)

void	setAxisymmetricCoordAxis (const MooseEnum &rz_coord_axis)

void	setGeneralAxisymmetricCoordAxes (const std::vector< SubdomainName > &blocks, const std::vector< std::pair< Point, RealVectorValue >> &axes)

const std::pair< Point, RealVectorValue > &	getGeneralAxisymmetricCoordAxis (SubdomainID subdomain_id) const

bool	usingGeneralAxisymmetricCoordAxes () const

unsigned int	getAxisymmetricRadialCoord () const

void	checkCoordinateSystems ()

void	setCoordData (const MooseMesh &other_mesh)

void	markFiniteVolumeInfoDirty ()

bool	isFiniteVolumeInfoDirty () const

MooseAppCoordTransform &	coordTransform ()

const MooseUnits &	lengthUnit () const

const std::unordered_map< std::pair< const Elem , unsigned short int >, const Elem > &	getLowerDElemMap () const

bool	isSplit () const

void	buildFiniteVolumeInfo () const

void	setupFiniteVolumeMeshData () const

void	doingPRefinement (bool doing_p_refinement)

bool	doingPRefinement () const

unsigned int	maxPLevel () const

unsigned int	maxHLevel () const

const std::vector< QpMap > &	getPRefinementMap (const Elem &elem) const

const std::vector< QpMap > &	getPRefinementSideMap (const Elem &elem) const

const std::vector< QpMap > &	getPCoarseningMap (const Elem &elem) const

const std::vector< QpMap > &	getPCoarseningSideMap (const Elem &elem) const

void	buildPRefinementAndCoarseningMaps (Assembly *assembly)

bool	isLowerD (const SubdomainID subdomain_id) const

bool	hasLowerD () const

const std::set< SubdomainID > &	interiorLowerDBlocks () const

const std::set< SubdomainID > &	boundaryLowerDBlocks () const

bool	getConstructNodeListFromSideList ()

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () 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 InputParameters &	parameters () 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

T	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

PerfGraph &	perfGraph ()

libMesh::ConstElemRange *	getActiveLocalElementRange ()

libMesh::ConstElemRange *	getActiveLocalElementRange ()

libMesh::NodeRange *	getActiveNodeRange ()

libMesh::NodeRange *	getActiveNodeRange ()

SemiLocalNodeRange *	getActiveSemiLocalNodeRange () const

SemiLocalNodeRange *	getActiveSemiLocalNodeRange () const

libMesh::ConstNodeRange *	getLocalNodeRange ()

libMesh::ConstNodeRange *	getLocalNodeRange ()

libMesh::StoredRange< MooseMesh::const_bnd_node_iterator, const BndNode >	getBoundaryNodeRange ()

libMesh::StoredRange< MooseMesh::const_bnd_node_iterator, const BndNode >	getBoundaryNodeRange ()

libMesh::StoredRange< MooseMesh::const_bnd_elem_iterator, const BndElement >	getBoundaryElementRange ()

libMesh::StoredRange< MooseMesh::const_bnd_elem_iterator, const BndElement >	getBoundaryElementRange ()

virtual Real	getMinInDimension (unsigned int component) const

virtual Real	getMinInDimension (unsigned int component) const

virtual Real	getMaxInDimension (unsigned int component) const

virtual Real	getMaxInDimension (unsigned int component) const

bool	isCustomPartitionerRequested () const

bool	isCustomPartitionerRequested () const

void	setIsCustomPartitionerRequested (bool cpr)

void	setIsCustomPartitionerRequested (bool cpr)

unsigned int	nFace () const

unsigned int	nFace () const

const std::vector< const FaceInfo *> &	faceInfo () const

const FaceInfo *	faceInfo (const Elem *elem, unsigned int side) const

const std::vector< const FaceInfo *> &	faceInfo () const

const FaceInfo *	faceInfo (const Elem *elem, unsigned int side) const

face_info_iterator	ownedFaceInfoBegin ()

face_info_iterator	ownedFaceInfoBegin ()

face_info_iterator	ownedFaceInfoEnd ()

face_info_iterator	ownedFaceInfoEnd ()

elem_info_iterator	ownedElemInfoBegin ()

elem_info_iterator	ownedElemInfoBegin ()

elem_info_iterator	ownedElemInfoEnd ()

elem_info_iterator	ownedElemInfoEnd ()

const ElemInfo &	elemInfo (const dof_id_type id) const

const ElemInfo &	elemInfo (const dof_id_type id) const

const std::vector< const ElemInfo *> &	elemInfoVector () const

const std::vector< const ElemInfo *> &	elemInfoVector () const

const std::vector< FaceInfo > &	allFaceInfo () const

const std::vector< FaceInfo > &	allFaceInfo () 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	rodPositions (std::vector< Point > &positions, unsigned int nrings, Real assembly_pitch, Point center)
	Calculates and stores the Pin positions/centers for a hexagonal assembly containing the given number of rings in a triangular/alternating row grid spaced 'assembly_pitch' apart. More...

static void	generateZGrid (Real unheated_length_entry, Real heated_length, Real unheated_length_exit, unsigned int n_cells, std::vector< Real > &z_grid)
	Generate the spacing in z-direction using heated and unheated lengths. More...

static void	setSubdomainName (MeshBase &mesh, SubdomainID subdomain_id, const SubdomainName &name)

static void	changeBoundaryId (MeshBase &mesh, const boundary_id_type old_id, const boundary_id_type new_id, bool delete_prev)

static MooseEnum	partitioning ()

static MooseEnum	elemTypes ()

static void	setPartitioner (MeshBase &mesh_base, MooseEnum &partitioner, bool use_distributed_mesh, const InputParameters &params, MooseObject &context_obj)

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static const unsigned int	N_CORNERS = 6
	number of corners in the duct x-sec More...

Protected Types
typedef std::vector< BndNode *>::iterator	bnd_node_iterator_imp

typedef std::vector< BndNode *>::const_iterator	const_bnd_node_iterator_imp

typedef std::vector< BndElement *>::iterator	bnd_elem_iterator_imp

typedef std::vector< BndElement *>::const_iterator	const_bnd_elem_iterator_imp

Protected Member Functions
void	cacheInfo ()

void	freeBndNodes ()

void	freeBndElems ()

void	setPartitionerHelper (MeshBase *mesh=nullptr)

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

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

Protected Attributes
unsigned int	_n_rings
	number of rings of fuel pins More...

unsigned int	_n_channels
	number of subchannels More...

Real	_flat_to_flat
	the distance between flat surfaces of the duct facing each other More...

Real	_duct_to_pin_gap
	the gap thickness between the duct and peripheral fuel pins More...

std::vector< std::vector< Node * > >	_nodes
	nodes More...

std::vector< Node * >	_duct_nodes
	A list of all mesh nodes that form the (elements of) the hexagonal duct mesh that surrounds the pins/intertices. More...

std::map< Node , Node >	_chan_to_duct_node_map
	A map for providing the closest/corresponding duct node associated with each interstice node. More...

std::map< Node , Node >	_duct_node_to_chan_map
	A map for providing the closest/corresponding interstice node associated with each duct node. More...

std::vector< std::pair< unsigned int, unsigned int > >	_gap_to_chan_map
	stores the channel pairs for each gap More...

std::vector< std::pair< unsigned int, unsigned int > >	_gap_to_pin_map
	stores the fuel pins belonging to each gap More...

std::vector< std::vector< unsigned int > >	_chan_to_gap_map
	stores the gaps that forms each interstice More...

std::vector< std::vector< Real > >	_sign_id_crossflow_map
	Defines the global cross-flow direction -1 or 1 for each interstice and for all gaps that are belonging to the corresponding interstice. More...

std::vector< Real >	_gij_map
	gap size More...

std::vector< std::vector< Real > >	_subchannel_position
	x,y coordinates of the interstice More...

std::vector< Point >	_pin_position
	x,y coordinates of the fuel pins More...

std::vector< std::vector< Real > >	_pins_in_rings
	fuel pins that are belonging to each ring More...

std::vector< std::vector< unsigned int > >	_chan_to_pin_map
	stores the fuel pins belonging to each interstice More...

unsigned int	_n_assemblies
	number of assemblies More...

unsigned int	_n_gaps
	number of gaps More...

std::vector< EChannelType >	_subch_type
	interstice type More...

std::vector< EChannelType >	_gap_type
	gap type More...

std::vector< std::pair< unsigned int, unsigned int > >	_gap_pairs_sf
	sweeping flow model gap pairs per channel to specify directional edge flow More...

std::vector< std::pair< unsigned int, unsigned int > >	_chan_pairs_sf
	sweeping flow model channel pairs to specify directional edge flow More...

std::vector< std::vector< unsigned int > >	_pin_to_chan_map
	TODO: channel indices corresponding to a given pin index. More...

bool	_tight_side_bypass
	whether the side bypass shape follows the assemblies More...

Real	_unheated_length_entry
	unheated length of the fuel Pin at the entry of the assembly More...

Real	_heated_length
	heated length of the fuel Pin More...

Real	_unheated_length_exit
	unheated length of the fuel Pin at the exit of the assembly More...

std::vector< Real >	_z_grid
	axial location of nodes More...

std::vector< std::vector< Real > >	_k_grid
	axial form loss coefficient per computational cell More...

Real	_kij
	Lateral form loss coefficient. More...

Real	_assembly_pitch
	Distance between neighboring assemblies. More...

Real	_assembly_side_x
	Size of assembly sides. More...

Real	_assembly_side_y

unsigned int	_n_cells
	number of axial cells More...

	X

	Y

	Z

	MIN

	MAX

std::vector< std::unique_ptr< libMesh::GhostingFunctor > >	_ghosting_functors

std::vector< std::shared_ptr< RelationshipManager > >	_relationship_managers

bool	_built_from_other_mesh

ParallelType	_parallel_type

bool	_use_distributed_mesh

bool	_distribution_overridden

bool	_parallel_type_overridden

std::unique_ptr< libMesh::MeshBase >	_mesh

MooseEnum	_partitioner_name

bool	_partitioner_overridden

std::unique_ptr< libMesh::Partitioner >	_custom_partitioner

bool	_custom_partitioner_requested

unsigned int	_uniform_refine_level

bool	_skip_refine_when_use_split

bool	_skip_deletion_repartition_after_refine

bool	_is_changed

bool	_is_nemesis

bool	_moose_mesh_prepared

std::unique_ptr< ConstElemPointerRange >	_refined_elements

std::unique_ptr< ConstElemPointerRange >	_coarsened_elements

std::map< const Elem , std::vector< const Elem > >	_coarsened_element_children

std::set< Node *>	_semilocal_node_list

std::unique_ptr< libMesh::ConstElemRange >	_active_local_elem_range

std::unique_ptr< SemiLocalNodeRange >	_active_semilocal_node_range

std::unique_ptr< libMesh::NodeRange >	_active_node_range

std::unique_ptr< libMesh::ConstNodeRange >	_local_node_range

std::unique_ptr< libMesh::StoredRange< MooseMesh::const_bnd_node_iterator, const BndNode *> >	_bnd_node_range

std::unique_ptr< libMesh::StoredRange< MooseMesh::const_bnd_elem_iterator, const BndElement *> >	_bnd_elem_range

std::map< dof_id_type, std::vector< dof_id_type > >	_node_to_elem_map

bool	_node_to_elem_map_built

std::map< dof_id_type, std::vector< dof_id_type > >	_node_to_active_semilocal_elem_map

bool	_node_to_active_semilocal_elem_map_built

std::set< SubdomainID >	_mesh_subdomains

std::unique_ptr< std::map< BoundaryID, RealVectorValue > >	_boundary_to_normal_map

std::vector< BndNode *>	_bnd_nodes

std::map< boundary_id_type, std::set< dof_id_type > >	_bnd_node_ids

std::vector< BndElement *>	_bnd_elems

std::unordered_map< boundary_id_type, std::unordered_set< dof_id_type > >	_bnd_elem_ids

std::map< dof_id_type, Node *>	_quadrature_nodes

std::map< dof_id_type, std::map< unsigned int, std::map< dof_id_type, Node *> > >	_elem_to_side_to_qp_to_quadrature_nodes

std::vector< BndNode >	_extra_bnd_nodes

std::map< dof_id_type, std::set< SubdomainID > >	_block_node_list

std::map< boundary_id_type, std::vector< dof_id_type > >	_node_set_nodes

std::set< unsigned int >	_ghosted_boundaries

std::vector< Real >	_ghosted_boundaries_inflation

unsigned int	_patch_size

unsigned int	_ghosting_patch_size

unsigned int	_max_leaf_size

Moose::PatchUpdateType	_patch_update_strategy

std::vector< Node *>	_node_map

bool	_regular_orthogonal_mesh

std::vector< std::vector< Real > >	_bounds

std::vector< std::pair< BoundaryID, BoundaryID > >	_paired_boundary

const bool	_is_split

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

MooseApp &	_pg_moose_app

const std::string	_prefix

std::set< BoundaryID >	_mesh_boundary_ids

std::set< BoundaryID >	_mesh_sideset_ids

std::set< BoundaryID >	_mesh_nodeset_ids

const Parallel::Communicator &	_communicator

Friends
class	SCMTriInterWrapperMeshGenerator

Detailed Description

Mesh class for triangular, edge and corner inter_wrappers for hexagonal lattice fuel assemblies.

Definition at line 17 of file TriInterWrapperMesh.h.

Constructor & Destructor Documentation

◆ TriInterWrapperMesh() [1/2]

TriInterWrapperMesh::TriInterWrapperMesh ( const InputParameters & parameters )

Definition at line 25 of file TriInterWrapperMesh.C.

                                                                        : InterWrapperMesh(params)
 {
 }

◆ TriInterWrapperMesh() [2/2]

TriInterWrapperMesh::TriInterWrapperMesh ( const TriInterWrapperMesh & other_mesh )

Definition at line 29 of file TriInterWrapperMesh.C.

   : InterWrapperMesh(other_mesh),
     _n_rings(other_mesh._n_rings),
     _n_channels(other_mesh._n_channels),
     _flat_to_flat(other_mesh._flat_to_flat),
     _duct_to_pin_gap(other_mesh._duct_to_pin_gap),
     _nodes(other_mesh._nodes),
     _duct_nodes(other_mesh._duct_nodes),
     _chan_to_duct_node_map(other_mesh._chan_to_duct_node_map),
     _duct_node_to_chan_map(other_mesh._duct_node_to_chan_map),
     _gap_to_chan_map(other_mesh._gap_to_chan_map),
     _gap_to_pin_map(other_mesh._gap_to_pin_map),
     _chan_to_gap_map(other_mesh._chan_to_gap_map),
     _sign_id_crossflow_map(other_mesh._sign_id_crossflow_map),
     _gij_map(other_mesh._gij_map),
     _subchannel_position(other_mesh._subchannel_position),
     _pin_position(other_mesh._pin_position),
     _pins_in_rings(other_mesh._pins_in_rings),
     _chan_to_pin_map(other_mesh._chan_to_pin_map),
     _n_assemblies(other_mesh._n_assemblies),
     _n_gaps(other_mesh._n_gaps),
     _subch_type(other_mesh._subch_type),
     _gap_type(other_mesh._gap_type),
     _gap_pairs_sf(other_mesh._gap_pairs_sf),
     _chan_pairs_sf(other_mesh._chan_pairs_sf),
     _pin_to_chan_map(other_mesh._pin_to_chan_map),
     _tight_side_bypass(other_mesh._tight_side_bypass)
 {
 }

Member Function Documentation

◆ buildMesh()

void TriInterWrapperMesh::buildMesh ( )

overridevirtual

Implements MooseMesh.

Definition at line 94 of file TriInterWrapperMesh.C.

95 {

96 }

◆ channelIndex()

unsigned int TriInterWrapperMesh::channelIndex ( const Point & point ) const

overridevirtual

Implements InterWrapperMesh.

Definition at line 87 of file TriInterWrapperMesh.C.

 {
   mooseAssert(false, "This is not implemented");
   return 0;
 }

◆ ductMeshExist()

virtual bool TriInterWrapperMesh::ductMeshExist ( ) const

inlineoverridevirtual

Return if duct Mesh exists or not.

Implements InterWrapperMesh.

Definition at line 30 of file TriInterWrapperMesh.h.

30 { return false; }

◆ generateZGrid()

void InterWrapperMesh::generateZGrid	(	Real	unheated_length_entry,
		Real	heated_length,
		Real	unheated_length_exit,
		unsigned int	n_cells,
		std::vector< Real > &	z_grid
	)

staticinherited

Generate the spacing in z-direction using heated and unheated lengths.

Definition at line 38 of file InterWrapperMesh.C.

Referenced by SCMTriInterWrapperMeshGenerator::SCMTriInterWrapperMeshGenerator().

 {
   Real L = unheated_length_entry + heated_length + unheated_length_exit;
   Real dz = L / n_cells;
   for (unsigned int i = 0; i < n_cells + 1; i++)
     z_grid.push_back(dz * i);
 }

◆ getChannelGaps()

virtual const std::vector<unsigned int>& TriInterWrapperMesh::getChannelGaps ( unsigned int i_chan ) const

inlineoverridevirtual

Return a vector of gap indices for a given channel index.

Implements InterWrapperMesh.

Definition at line 54 of file TriInterWrapperMesh.h.

   {
     return _chan_to_gap_map[i_chan];
   }

◆ getChannelNode()

virtual Node* TriInterWrapperMesh::getChannelNode	(	unsigned int	i_chan,
		unsigned	iz
	)		const

inlineoverridevirtual

Get the inter-wrapper mesh node for a given channel index and elevation index.

Implements InterWrapperMesh.

Definition at line 43 of file TriInterWrapperMesh.h.

   {
     return _nodes[i_chan][iz];
   }

◆ getChannelPins()

virtual const std::vector<unsigned int>& TriInterWrapperMesh::getChannelPins ( unsigned int i_chan ) const

inlineoverridevirtual

Return a vector of pin indices for a given channel index.

Implements InterWrapperMesh.

Definition at line 89 of file TriInterWrapperMesh.h.

   {
     return _chan_to_pin_map[i_chan];
   }

◆ getCrossflowSign()

virtual const Real& TriInterWrapperMesh::getCrossflowSign	(	unsigned int	i_chan,
		unsigned int	i_local
	)		const

inlineoverridevirtual

Return a signs for the cross flow given a inter-wrapper index and local neighbor index.

Implements InterWrapperMesh.

Definition at line 59 of file TriInterWrapperMesh.h.

   {
     return _sign_id_crossflow_map[i_chan][i_local];
   }

◆ getDuctToPinGap()

virtual const Real& TriInterWrapperMesh::getDuctToPinGap ( ) const

inlinevirtual

Definition at line 32 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), TriInterWrapperWettedPerimIC::value(), and TriInterWrapperFlowAreaIC::value().

32 { return _duct_to_pin_gap; }

TriInterWrapperMesh::_duct_to_pin_gap

Real _duct_to_pin_gap

the gap thickness between the duct and peripheral fuel pins

Definition: TriInterWrapperMesh.h:112

◆ getGapChannels()

virtual const std::pair<unsigned int, unsigned int>& TriInterWrapperMesh::getGapChannels ( unsigned int i_gap ) const

inlineoverridevirtual

Return a pair of inter-wrapper indices for a given gap index.

Implements InterWrapperMesh.

Definition at line 50 of file TriInterWrapperMesh.h.

   {
     return _gap_to_chan_map[i_gap];
   }

◆ getGapMap()

virtual const std::vector<double>& TriInterWrapperMesh::getGapMap ( ) const

inlineoverridevirtual

Return a map with gap sizes.

Implements InterWrapperMesh.

Definition at line 58 of file TriInterWrapperMesh.h.

58 { return _gij_map; }

TriInterWrapperMesh::_gij_map

std::vector< Real > _gij_map

gap size

Definition: TriInterWrapperMesh.h:138

◆ getGapWidth()

virtual Real TriInterWrapperMesh::getGapWidth ( unsigned int gap_index ) const

inlineoverridevirtual

Return gap width for a given gap index.

Implements InterWrapperMesh.

Definition at line 72 of file TriInterWrapperMesh.h.

72 { return _gij_map[gap_index]; }

TriInterWrapperMesh::_gij_map

std::vector< Real > _gij_map

gap size

Definition: TriInterWrapperMesh.h:138

◆ getHeatedLength()

virtual const Real& InterWrapperMesh::getHeatedLength ( ) const

inlinevirtualinherited

Return heated length.

Definition at line 135 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), InterWrapper1PhaseProblem::computeh(), QuadInterWrapperPowerIC::initialSetup(), TriInterWrapperPowerIC::initialSetup(), QuadInterWrapperPowerIC::QuadInterWrapperPowerIC(), TriInterWrapperPowerIC::TriInterWrapperPowerIC(), QuadInterWrapperPowerIC::value(), and TriInterWrapperPowerIC::value().

135 { return _heated_length; }

InterWrapperMesh::_heated_length

Real _heated_length

heated length of the fuel Pin

Definition: InterWrapperMesh.h:170

◆ getHeatedLengthEntry()

virtual const Real& InterWrapperMesh::getHeatedLengthEntry ( ) const

inlinevirtualinherited

Return unheated length at entry.

Definition at line 130 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), InterWrapper1PhaseProblem::computeh(), QuadInterWrapperPowerIC::initialSetup(), TriInterWrapperPowerIC::initialSetup(), QuadInterWrapperPowerIC::value(), and TriInterWrapperPowerIC::value().

130 { return _unheated_length_entry; }

InterWrapperMesh::_unheated_length_entry

Real _unheated_length_entry

unheated length of the fuel Pin at the entry of the assembly

Definition: InterWrapperMesh.h:168

◆ getHeatedLengthExit()

virtual const Real& InterWrapperMesh::getHeatedLengthExit ( ) const

inlinevirtualinherited

Return unheated length at exit.

Definition at line 140 of file InterWrapperMesh.h.

140 { return _unheated_length_exit; }

InterWrapperMesh::_unheated_length_exit

Real _unheated_length_exit

unheated length of the fuel Pin at the exit of the assembly

Definition: InterWrapperMesh.h:172

◆ getIsTightSide()

virtual const bool& TriInterWrapperMesh::getIsTightSide ( ) const

inlinevirtual

Definition at line 94 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapperWettedPerimIC::value(), and TriInterWrapperFlowAreaIC::value().

94 { return _tight_side_bypass; }

TriInterWrapperMesh::_tight_side_bypass

bool _tight_side_bypass

whether the side bypass shape follows the assemblies

Definition: TriInterWrapperMesh.h:162

◆ getKGrid()

virtual const std::vector<std::vector<Real> >& InterWrapperMesh::getKGrid ( ) const

inlinevirtualinherited

Get axial cell location and value of loss coefficient.

Definition at line 38 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeDP(), and InterWrapper1PhaseProblem::computeDP().

38 { return _k_grid; }

InterWrapperMesh::_k_grid

std::vector< std::vector< Real > > _k_grid

axial form loss coefficient per computational cell

Definition: InterWrapperMesh.h:176

◆ getKij()

virtual const Real& InterWrapperMesh::getKij ( ) const

inlinevirtualinherited

Return lateral loss coefficient.

Definition at line 43 of file InterWrapperMesh.h.

43 { return _kij; }

InterWrapperMesh::_kij

Real _kij

Lateral form loss coefficient.

Definition: InterWrapperMesh.h:178

◆ getNumOfAssemblies()

virtual const unsigned int& TriInterWrapperMesh::getNumOfAssemblies ( ) const

inlineoverridevirtual

Return the number of assemblies.

Implements InterWrapperMesh.

Definition at line 25 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapperPowerIC::initialSetup(), and TriInterWrapperPowerIC::TriInterWrapperPowerIC().

25 { return _n_assemblies; }

TriInterWrapperMesh::_n_assemblies

unsigned int _n_assemblies

number of assemblies

Definition: TriInterWrapperMesh.h:148

◆ getNumOfAxialCells()

virtual const unsigned int& InterWrapperMesh::getNumOfAxialCells ( ) const

inlinevirtualinherited

Return the number of axial cells.

Definition at line 48 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), QuadInterWrapperPowerIC::initialSetup(), and TriInterWrapperPowerIC::initialSetup().

48 { return _n_cells; }

InterWrapperMesh::_n_cells

unsigned int _n_cells

number of axial cells

Definition: InterWrapperMesh.h:185

◆ getNumOfChannels()

virtual const unsigned int& TriInterWrapperMesh::getNumOfChannels ( ) const

inlineoverridevirtual

Return the number of channels per layer.

Implements InterWrapperMesh.

Definition at line 47 of file TriInterWrapperMesh.h.

47 { return _n_channels; }

TriInterWrapperMesh::_n_channels

unsigned int _n_channels

number of subchannels

Definition: TriInterWrapperMesh.h:108

◆ getNumOfGapsPerLayer()

virtual const unsigned int& TriInterWrapperMesh::getNumOfGapsPerLayer ( ) const

inlineoverridevirtual

Return the number of gaps per layer.

Implements InterWrapperMesh.

Definition at line 48 of file TriInterWrapperMesh.h.

48 { return _n_gaps; }

TriInterWrapperMesh::_n_gaps

unsigned int _n_gaps

number of gaps

Definition: TriInterWrapperMesh.h:150

◆ getPinChannels()

virtual const std::vector<unsigned int>& TriInterWrapperMesh::getPinChannels ( unsigned int i_pin ) const

inlineoverridevirtual

Return a vector of channel indices for a given Pin index.

Implements InterWrapperMesh.

Definition at line 84 of file TriInterWrapperMesh.h.

   {
     return _pin_to_chan_map[i_pin];
   }

◆ getPinIndex()

virtual const unsigned int& TriInterWrapperMesh::getPinIndex	(	const unsigned int	channel_idx,
		const unsigned int	neighbor_idx
	)

inlinevirtual

Return Pin index given inter_wrapper index and local neighbor index.

Definition at line 37 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapperPowerIC::value().

   {
     return _chan_to_pin_map[channel_idx][neighbor_idx];
   }

◆ getPinIndexFromPoint()

unsigned int TriInterWrapperMesh::getPinIndexFromPoint ( const Point & p ) const

overridevirtual

Return a pin index for a given physical point p

Implements InterWrapperMesh.

Definition at line 99 of file TriInterWrapperMesh.C.

 {
   mooseAssert(false, "This is not implemented");
   return 0;
 }

◆ getPinNode()

virtual Node* TriInterWrapperMesh::getPinNode	(	unsigned	i_pin,
		unsigned	iz
	)		const

inlineoverridevirtual

Get the pin mesh node for a given pin index and elevation index.

Implements InterWrapperMesh.

Definition at line 27 of file TriInterWrapperMesh.h.

27 { return nullptr; }

◆ getPitch()

virtual const Real& InterWrapperMesh::getPitch ( ) const

inlinevirtualinherited

Return the pitch between 2 inter-wrappers.

Reimplemented in QuadInterWrapperMesh.

Definition at line 114 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), InterWrapper1PhaseProblem::computeWij(), and TriInterWrapperFlowAreaIC::value().

114 { return _assembly_pitch; }

InterWrapperMesh::_assembly_pitch

Real _assembly_pitch

Distance between neighboring assemblies.

Definition: InterWrapperMesh.h:180

◆ getSideX()

virtual const Real& InterWrapperMesh::getSideX ( ) const

inlinevirtualinherited

Return side lengths of the assembly.

Definition at line 119 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh(), InterWrapper1PhaseProblem::externalSolve(), QuadInterWrapperWettedPerimIC::value(), TriInterWrapperWettedPerimIC::value(), QuadInterWrapperFlowAreaIC::value(), and TriInterWrapperFlowAreaIC::value().

119 { return _assembly_side_x; }

InterWrapperMesh::_assembly_side_x

Real _assembly_side_x

Size of assembly sides.

Definition: InterWrapperMesh.h:182

◆ getSideY()

virtual const Real& InterWrapperMesh::getSideY ( ) const

inlinevirtualinherited

Definition at line 120 of file InterWrapperMesh.h.

Referenced by InterWrapper1PhaseProblem::externalSolve(), QuadInterWrapperWettedPerimIC::value(), and QuadInterWrapperFlowAreaIC::value().

120 { return _assembly_side_y; }

InterWrapperMesh::_assembly_side_y

Real _assembly_side_y

Definition: InterWrapperMesh.h:183

◆ getSubchannelIndexFromPoint()

unsigned int TriInterWrapperMesh::getSubchannelIndexFromPoint ( const Point & p ) const

overridevirtual

Return a inter-wrapper index for a given physical point p

Implements InterWrapperMesh.

Definition at line 66 of file TriInterWrapperMesh.C.

Referenced by TriInterWrapperWettedPerimIC::value(), TriInterWrapperFlowAreaIC::value(), and TriInterWrapperPowerIC::value().

 {
   Real distance0 = 1.0e+8;
   Real distance1;
   unsigned int j = 0;
 
   for (unsigned int i = 0; i < _n_channels; i++)
   {
     distance1 = std::sqrt(std::pow((p(0) - _subchannel_position[i][0]), 2.0) +
                           std::pow((p(1) - _subchannel_position[i][1]), 2.0));
 
     if (distance1 < distance0)
     {
       j = i;
       distance0 = distance1;
     } // if
   }   // for
   return j;
 }

◆ getSubchannelType()

virtual EChannelType TriInterWrapperMesh::getSubchannelType ( unsigned int index ) const

inlineoverridevirtual

Return the type of the inter-wrapper for given inter-wrapper index.

Implements InterWrapperMesh.

Definition at line 67 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapperWettedPerimIC::value(), TriInterWrapperFlowAreaIC::value(), and TriInterWrapperPowerIC::value().

   {
     return _subch_type[index];
   }

◆ getSweepFlowChans()

virtual const std::pair<unsigned int, unsigned int>& TriInterWrapperMesh::getSweepFlowChans ( unsigned int i_chan ) const

inlinevirtual

Definition at line 79 of file TriInterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeh().

   {
     return _chan_pairs_sf[i_chan];
   }

◆ getSweepFlowGaps()

virtual const std::pair<unsigned int, unsigned int>& TriInterWrapperMesh::getSweepFlowGaps ( unsigned int i_chan ) const

inlinevirtual

Definition at line 74 of file TriInterWrapperMesh.h.

   {
     return _gap_pairs_sf[i_chan];
   }

◆ getZGrid()

virtual const std::vector<Real>& InterWrapperMesh::getZGrid ( ) const

inlinevirtualinherited

Get axial location of layers.

Definition at line 28 of file InterWrapperMesh.h.

Referenced by TriInterWrapper1PhaseProblem::computeDP(), TriInterWrapper1PhaseProblem::computeh(), QuadInterWrapperPowerIC::initialSetup(), and TriInterWrapperPowerIC::initialSetup().

28 { return _z_grid; }

InterWrapperMesh::_z_grid

std::vector< Real > _z_grid

axial location of nodes

Definition: InterWrapperMesh.h:174

◆ getZIndex()

unsigned int InterWrapperMesh::getZIndex ( const Point & point ) const

virtualinherited

Get axial index of point.

Definition at line 51 of file InterWrapperMesh.C.

Referenced by InterWrapperSolutionTransfer::getFromNode().

 {
   if (_z_grid.size() == 0)
     mooseError("_z_grid is empty.");
 
   if (point(2) <= _z_grid[0])
     return 0;
   if (point(2) >= _z_grid[_z_grid.size() - 1])
     return _z_grid.size() - 1;
 
   unsigned int lo = 0;
   unsigned int hi = _z_grid.size();
   while (lo < hi)
   {
     unsigned int mid = (lo + hi) / 2;
     if (std::abs(_z_grid[mid] - point(2)) < 1e-5)
       return mid;
     else if (_z_grid[mid] < point(2))
       lo = mid;
     else
       hi = mid;
   }
   return lo;
 }

◆ pinIndex()

unsigned int TriInterWrapperMesh::pinIndex ( const Point & p ) const

overridevirtual

Implements InterWrapperMesh.

Definition at line 106 of file TriInterWrapperMesh.C.

 {
   mooseAssert(false, "This is not implemented");
   return 0;
 }

◆ pinMeshExist()

virtual bool TriInterWrapperMesh::pinMeshExist ( ) const

inlineoverridevirtual

Return if Pin Mesh exists or not.

Implements InterWrapperMesh.

Definition at line 29 of file TriInterWrapperMesh.h.

29 { return false; }

◆ rodPositions()

void TriInterWrapperMesh::rodPositions	(	std::vector< Point > &	positions,
		unsigned int	nrings,
		Real	assembly_pitch,
		Point	center
	)

static

Calculates and stores the Pin positions/centers for a hexagonal assembly containing the given number of rings in a triangular/alternating row grid spaced 'assembly_pitch' apart.

The points are generated such that the duct is centered at the given center point.

Defining parameters

Definition at line 113 of file TriInterWrapperMesh.C.

Referenced by SCMTriInterWrapperMeshGenerator::SCMTriInterWrapperMeshGenerator().

 {
   // distance: it is the distance to the next Pin
   //
   Real theta = 0.0;
   Real dtheta = 0.0;
   Real distance = 0.0;
   Real theta1 = 0.0;
   Real theta_corrected = 0.0;
   Real pi = libMesh::pi;
   unsigned int k = 0;
   positions.emplace_back(0.0, 0.0);
   for (unsigned int i = 1; i < nrings; i++)
   {
     dtheta = 2.0 * pi / (i * 6);
     theta = 0.0;
     for (unsigned int j = 0; j < i * 6; j++)
     {
       k = k + 1;
       theta1 = fmod(theta + 1.0e-10, pi / 3.0);
       distance = std::sqrt(
           (pow(i * assembly_pitch, 2) + pow(theta1 / dtheta * assembly_pitch, 2) -
            2.0 * i * assembly_pitch * (theta1 / dtheta * assembly_pitch) * std::cos(pi / 3.0)));
       theta_corrected = std::acos(1.0 / (i * assembly_pitch) / distance / 2.0 *
                                   (std::pow(i * assembly_pitch, 2) + std::pow(distance, 2) -
                                    std::pow(theta1 / dtheta * assembly_pitch, 2)));
       if (theta1 < 1.0e-6)
       {
         theta_corrected = theta;
       }
       else
       {
         if (theta > pi / 3.0 && theta <= 2.0 / 3.0 * pi)
           theta_corrected = theta_corrected + pi / 3.0;
         else if (theta > 2.0 / 3.0 * pi && theta <= pi)
           theta_corrected = theta_corrected + 2.0 / 3.0 * pi;
         else if (theta > pi && theta <= 4.0 / 3.0 * pi)
           theta_corrected = theta_corrected + pi;
         else if (theta > 4.0 / 3.0 * pi && theta <= 5.0 / 3.0 * pi)
           theta_corrected = theta_corrected + 4.0 / 3.0 * pi;
         else if (theta > 5.0 / 3.0 * pi && theta <= 2.0 * pi)
           theta_corrected = theta_corrected + 5.0 / 3.0 * pi;
       }
       positions.emplace_back(center(0) + distance * std::cos(theta_corrected),
                              center(1) + distance * std::sin(theta_corrected));
       theta = theta + dtheta;
     } // j
   }   // i
 }

◆ safeClone()

std::unique_ptr< MooseMesh > TriInterWrapperMesh::safeClone ( ) const

overridevirtual

Implements MooseMesh.

Definition at line 60 of file TriInterWrapperMesh.C.

 {
   return _app.getFactory().copyConstruct(*this);
 }

◆ setChannelToDuctMaps()

void TriInterWrapperMesh::setChannelToDuctMaps ( const std::vector< Node *> & duct_nodes )

Setup the internal maps when there is a outside duct present.

Definition at line 168 of file TriInterWrapperMesh.C.

 {
   const Real tol = 1e-10;
   for (size_t i = 0; i < duct_nodes.size(); i++)
   {
     int min_chan = 0;
     Real min_dist = std::numeric_limits<double>::max();
     Point ductpos((*duct_nodes[i])(0), (*duct_nodes[i])(1), 0);
     for (size_t j = 0; j < _subchannel_position.size(); j++)
     {
       Point chanpos(_subchannel_position[j][0], _subchannel_position[j][1], 0);
       auto dist = (chanpos - ductpos).norm();
       if (dist < min_dist)
       {
         min_dist = dist;
         min_chan = j;
       }
     }
 
     Node * chan_node = nullptr;
     for (auto cn : _nodes[min_chan])
     {
       if (std::abs((*cn)(2) - (*duct_nodes[i])(2)) < tol)
       {
         chan_node = cn;
         break;
       }
     }
 
     if (chan_node == nullptr)
       mooseError("failed to find matching channel node for duct node");
 
     _duct_node_to_chan_map[duct_nodes[i]] = chan_node;
     _chan_to_duct_node_map[chan_node] = duct_nodes[i];
   }
 
   _duct_nodes = duct_nodes;
 }

◆ validParams()

InputParameters TriInterWrapperMesh::validParams ( )

static

Definition at line 17 of file TriInterWrapperMesh.C.

 {
   InputParameters params = InterWrapperMesh::validParams();
   params.addClassDescription("Creates an inter-wrappper mesh container for a triangular "
                              "lattice subchannel arrangement");
   return params;
 }

Friends And Related Function Documentation

◆ SCMTriInterWrapperMeshGenerator

friend class SCMTriInterWrapperMeshGenerator

friend

Definition at line 178 of file TriInterWrapperMesh.h.

Member Data Documentation

◆ _assembly_pitch

Real InterWrapperMesh::_assembly_pitch

protectedinherited

Distance between neighboring assemblies.

Definition at line 180 of file InterWrapperMesh.h.

Referenced by QuadInterWrapperMesh::channelIndex(), QuadInterWrapperMesh::getPinIndexFromPoint(), QuadInterWrapperMesh::getPitch(), InterWrapperMesh::getPitch(), QuadInterWrapperMesh::getSubchannelIndexFromPoint(), and QuadInterWrapperMesh::pinIndex().

◆ _assembly_side_x

Real InterWrapperMesh::_assembly_side_x

protectedinherited

Size of assembly sides.

Definition at line 182 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getSideX().

◆ _assembly_side_y

Real InterWrapperMesh::_assembly_side_y

protectedinherited

Definition at line 183 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getSideY().

◆ _chan_pairs_sf

std::vector<std::pair<unsigned int, unsigned int> > TriInterWrapperMesh::_chan_pairs_sf

protected

sweeping flow model channel pairs to specify directional edge flow

Definition at line 158 of file TriInterWrapperMesh.h.

Referenced by getSweepFlowChans().

◆ _chan_to_duct_node_map

std::map<Node *, Node *> TriInterWrapperMesh::_chan_to_duct_node_map

protected

A map for providing the closest/corresponding duct node associated with each interstice node.

i.e. a map of interstice mesh nodes to duct mesh nodes.

Definition at line 121 of file TriInterWrapperMesh.h.

Referenced by setChannelToDuctMaps().

◆ _chan_to_gap_map

std::vector<std::vector<unsigned int> > TriInterWrapperMesh::_chan_to_gap_map

protected

stores the gaps that forms each interstice

Definition at line 131 of file TriInterWrapperMesh.h.

Referenced by getChannelGaps().

◆ _chan_to_pin_map

std::vector<std::vector<unsigned int> > TriInterWrapperMesh::_chan_to_pin_map

protected

stores the fuel pins belonging to each interstice

Definition at line 146 of file TriInterWrapperMesh.h.

Referenced by getChannelPins(), and getPinIndex().

◆ _duct_node_to_chan_map

std::map<Node *, Node *> TriInterWrapperMesh::_duct_node_to_chan_map

protected

A map for providing the closest/corresponding interstice node associated with each duct node.

i.e. a map of duct mesh nodes to interstice mesh nodes.

Definition at line 124 of file TriInterWrapperMesh.h.

Referenced by setChannelToDuctMaps().

◆ _duct_nodes

std::vector<Node *> TriInterWrapperMesh::_duct_nodes

protected

A list of all mesh nodes that form the (elements of) the hexagonal duct mesh that surrounds the pins/intertices.

Definition at line 118 of file TriInterWrapperMesh.h.

Referenced by setChannelToDuctMaps().

◆ _duct_to_pin_gap

Real TriInterWrapperMesh::_duct_to_pin_gap

protected

the gap thickness between the duct and peripheral fuel pins

Definition at line 112 of file TriInterWrapperMesh.h.

Referenced by getDuctToPinGap().

◆ _flat_to_flat

Real TriInterWrapperMesh::_flat_to_flat

protected

the distance between flat surfaces of the duct facing each other

Definition at line 110 of file TriInterWrapperMesh.h.

◆ _gap_pairs_sf

std::vector<std::pair<unsigned int, unsigned int> > TriInterWrapperMesh::_gap_pairs_sf

protected

sweeping flow model gap pairs per channel to specify directional edge flow

Definition at line 156 of file TriInterWrapperMesh.h.

Referenced by getSweepFlowGaps().

◆ _gap_to_chan_map

std::vector<std::pair<unsigned int, unsigned int> > TriInterWrapperMesh::_gap_to_chan_map

protected

stores the channel pairs for each gap

Definition at line 127 of file TriInterWrapperMesh.h.

Referenced by getGapChannels().

◆ _gap_to_pin_map

std::vector<std::pair<unsigned int, unsigned int> > TriInterWrapperMesh::_gap_to_pin_map

protected

stores the fuel pins belonging to each gap

Definition at line 129 of file TriInterWrapperMesh.h.

◆ _gap_type

std::vector<EChannelType> TriInterWrapperMesh::_gap_type

protected

gap type

Definition at line 154 of file TriInterWrapperMesh.h.

◆ _gij_map

std::vector<Real> TriInterWrapperMesh::_gij_map

protected

gap size

Definition at line 138 of file TriInterWrapperMesh.h.

Referenced by getGapMap(), and getGapWidth().

◆ _heated_length

Real InterWrapperMesh::_heated_length

protectedinherited

heated length of the fuel Pin

Definition at line 170 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getHeatedLength().

◆ _k_grid

std::vector<std::vector<Real> > InterWrapperMesh::_k_grid

protectedinherited

axial form loss coefficient per computational cell

Definition at line 176 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getKGrid().

◆ _kij

Real InterWrapperMesh::_kij

protectedinherited

Lateral form loss coefficient.

Definition at line 178 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getKij().

◆ _n_assemblies

unsigned int TriInterWrapperMesh::_n_assemblies

protected

number of assemblies

Definition at line 148 of file TriInterWrapperMesh.h.

Referenced by getNumOfAssemblies().

◆ _n_cells

unsigned int InterWrapperMesh::_n_cells

protectedinherited

number of axial cells

Definition at line 185 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getNumOfAxialCells().

◆ _n_channels

unsigned int TriInterWrapperMesh::_n_channels

protected

number of subchannels

Definition at line 108 of file TriInterWrapperMesh.h.

Referenced by getNumOfChannels(), and getSubchannelIndexFromPoint().

◆ _n_gaps

unsigned int TriInterWrapperMesh::_n_gaps

protected

number of gaps

Definition at line 150 of file TriInterWrapperMesh.h.

Referenced by getNumOfGapsPerLayer().

◆ _n_rings

unsigned int TriInterWrapperMesh::_n_rings

protected

number of rings of fuel pins

Definition at line 106 of file TriInterWrapperMesh.h.

◆ _nodes

std::vector<std::vector<Node *> > TriInterWrapperMesh::_nodes

protected

nodes

Definition at line 114 of file TriInterWrapperMesh.h.

Referenced by getChannelNode(), and setChannelToDuctMaps().

◆ _pin_position

std::vector<Point> TriInterWrapperMesh::_pin_position

protected

x,y coordinates of the fuel pins

Definition at line 142 of file TriInterWrapperMesh.h.

◆ _pin_to_chan_map

std::vector<std::vector<unsigned int> > TriInterWrapperMesh::_pin_to_chan_map

protected

TODO: channel indices corresponding to a given pin index.

Definition at line 160 of file TriInterWrapperMesh.h.

Referenced by getPinChannels().

◆ _pins_in_rings

std::vector<std::vector<Real> > TriInterWrapperMesh::_pins_in_rings

protected

fuel pins that are belonging to each ring

Definition at line 144 of file TriInterWrapperMesh.h.

◆ _sign_id_crossflow_map

std::vector<std::vector<Real> > TriInterWrapperMesh::_sign_id_crossflow_map

protected

Defines the global cross-flow direction -1 or 1 for each interstice and for all gaps that are belonging to the corresponding interstice.

Given a interstice and a gap, if the neighbor interstice index belonging to the same gap is lower, set it to -1, otherwise set it to 1.

Definition at line 136 of file TriInterWrapperMesh.h.

Referenced by getCrossflowSign().

◆ _subch_type

std::vector<EChannelType> TriInterWrapperMesh::_subch_type

protected

interstice type

Definition at line 152 of file TriInterWrapperMesh.h.

Referenced by getSubchannelType().

◆ _subchannel_position

std::vector<std::vector<Real> > TriInterWrapperMesh::_subchannel_position

protected

x,y coordinates of the interstice

Definition at line 140 of file TriInterWrapperMesh.h.

Referenced by getSubchannelIndexFromPoint(), and setChannelToDuctMaps().

◆ _tight_side_bypass

bool TriInterWrapperMesh::_tight_side_bypass

protected

whether the side bypass shape follows the assemblies

Definition at line 162 of file TriInterWrapperMesh.h.

Referenced by getIsTightSide().

◆ _unheated_length_entry

Real InterWrapperMesh::_unheated_length_entry

protectedinherited

unheated length of the fuel Pin at the entry of the assembly

Definition at line 168 of file InterWrapperMesh.h.

Referenced by SCMQuadInterWrapperMeshGenerator::generate(), SCMTriInterWrapperMeshGenerator::generate(), and InterWrapperMesh::getHeatedLengthEntry().

◆ _unheated_length_exit

Real InterWrapperMesh::_unheated_length_exit

protectedinherited

unheated length of the fuel Pin at the exit of the assembly

Definition at line 172 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getHeatedLengthExit().

◆ _z_grid

std::vector<Real> InterWrapperMesh::_z_grid

protectedinherited

axial location of nodes

Definition at line 174 of file InterWrapperMesh.h.

Referenced by InterWrapperMesh::getZGrid(), and InterWrapperMesh::getZIndex().

◆ N_CORNERS

const unsigned int TriInterWrapperMesh::N_CORNERS = 6

static

number of corners in the duct x-sec

Definition at line 181 of file TriInterWrapperMesh.h.

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

subchannel/include/mesh/TriInterWrapperMesh.h
subchannel/src/mesh/TriInterWrapperMesh.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ TriInterWrapperMesh() [1/2]

◆ TriInterWrapperMesh() [2/2]

Member Function Documentation

◆ buildMesh()

◆ channelIndex()

◆ ductMeshExist()

◆ generateZGrid()

◆ getChannelGaps()

◆ getChannelNode()

◆ getChannelPins()

◆ getCrossflowSign()

◆ getDuctToPinGap()

◆ getGapChannels()

◆ getGapMap()

◆ getGapWidth()

◆ getHeatedLength()

◆ getHeatedLengthEntry()

◆ getHeatedLengthExit()

◆ getIsTightSide()

◆ getKGrid()

◆ getKij()

◆ getNumOfAssemblies()

◆ getNumOfAxialCells()

◆ getNumOfChannels()

◆ getNumOfGapsPerLayer()

◆ getPinChannels()

◆ getPinIndex()

◆ getPinIndexFromPoint()

◆ getPinNode()

◆ getPitch()

◆ getSideX()

◆ getSideY()

◆ getSubchannelIndexFromPoint()

◆ getSubchannelType()

◆ getSweepFlowChans()

◆ getSweepFlowGaps()

◆ getZGrid()

◆ getZIndex()

◆ pinIndex()

◆ pinMeshExist()

◆ rodPositions()

◆ safeClone()

◆ setChannelToDuctMaps()

◆ validParams()

Friends And Related Function Documentation

◆ SCMTriInterWrapperMeshGenerator

Member Data Documentation

◆ _assembly_pitch

◆ _assembly_side_x

◆ _assembly_side_y

◆ _chan_pairs_sf

◆ _chan_to_duct_node_map

◆ _chan_to_gap_map

◆ _chan_to_pin_map

◆ _duct_node_to_chan_map

◆ _duct_nodes

◆ _duct_to_pin_gap

◆ _flat_to_flat

◆ _gap_pairs_sf

◆ _gap_to_chan_map

◆ _gap_to_pin_map

◆ _gap_type

◆ _gij_map

◆ _heated_length

◆ _k_grid

◆ _kij

◆ _n_assemblies

◆ _n_cells

◆ _n_channels