Creates the mesh of subchannels in a quadrilateral lattice. More...

#include <QuadSubChannelMesh.h>

Inheritance diagram for QuadSubChannelMesh:

Public Types
enum	ParallelType

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

typedef DataFileName	DataFileParameterType

Public Member Functions
	QuadSubChannelMesh (const InputParameters &parameters)

	QuadSubChannelMesh (const QuadSubChannelMesh &other_mesh)

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

void	buildMesh () override

void	computeAssemblyHydraulicParameters ()
	Compute undeformed bundle-average inlet hydraulic quantities from generated mesh geometry.

Real	getSubchannelFlowArea (unsigned int i_chan, Real z) const override
	Return undeformed flow area for a subchannel at an axial location, including any blockage reduction.

Real	getSubchannelWettedPerimeter (unsigned int i_chan) const override
	Return undeformed wetted perimeter for a subchannel.

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

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

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

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

unsigned int	getNumOfPins () const override
	Return the number of pins.

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

const std::pair< unsigned int, unsigned int > &	getGapPins (unsigned int i_gap) const override
	Return a pair of pin indices for a given gap index.

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

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

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

const Real &	getPitch () const override
	Return the undeformed pitch between 2 subchannels.

const Real &	getCrossflowSign (unsigned int i_chan, unsigned int i_local) const override
	Return a sign for the crossflow given a subchannel index and local neighbor index.

const unsigned int &	getNx () const
	Number of subchannels in the -x direction.

const unsigned int &	getNy () const
	Number of subchannels in the -y direction.

const Real &	getSideGap () const
	Returns the side gap, not to be confused with the gap between pins, this refers to the gap next to the duct.

unsigned int	getSubchannelIndexFromPoint (const Point &p) const override
	Return a subchannel index for a given physical point `p`

unsigned int	channelIndex (const Point &point) const override

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

unsigned int	pinIndex (const Point &p) const override

EChannelType	getSubchannelType (unsigned int index) const override
	Return the type of the subchannel for given subchannel index.

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

void	setGapWidth (unsigned int axial_index, unsigned int gap_index, Real gap_width)
	Set the gap width for a given axial cell and gap index.

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

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

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

virtual const std::vector< Real > &	getZBlockage () const
	Get axial location of blockage (in,out) [m].

virtual const std::vector< unsigned int > &	getIndexBlockage () const
	Get index of blocked subchannels.

virtual const std::vector< Real > &	getReductionBlockage () const
	Get area reduction of blocked subchannels.

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

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

Node *	getChannelNodeFromDuct (Node *duct_node) const
	Function that gets the channel node from the duct node.

Node *	getDuctNodeFromChannel (Node *channel_node) const
	Function that gets the duct node from the channel node.

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

bool	ductMeshExist () const
	Return if Duct Mesh exists or not.

virtual const Real &	getPinDiameter () const
	Return undeformed Pin diameter.

Real	getAssemblyFlowArea () const
	Return undeformed bundle inlet flow area.

Real	getAssemblyWettedPerimeter () const
	Return undeformed bundle inlet wetted perimeter.

Real	getAssemblyHydraulicDiameter () const
	Return undeformed bundle-average hydraulic diameter.

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

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

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

const std::vector< Node * > &	getDuctNodes () const
	Function that returns the vector with the duct nodes.

void	setChannelToDuctMaps (const std::vector< Node * > &duct_nodes)
	Function that sets the channel-to-duct maps.

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 std::vector< std::filesystem::path >	writeRecoveryFiles (const std::filesystem::path &file_base)

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

std::vector< std::vector< BoundaryID > >	getBoundaryIDs (const Elem *const elem) 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::unordered_map< dof_id_type, std::vector< dof_id_type > > &	nodeToElemMap ()

virtual bnd_node_iterator	bndNodesBegin ()

virtual bnd_node_iterator	bndNodesEnd ()

virtual bnd_elem_iterator	bndElemsBegin ()

virtual bnd_elem_iterator	bndElemsEnd ()

void	buildNodeListFromSideList ()

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

Moose::Kokkos::Mesh *	getKokkosMesh ()

const Moose::Kokkos::Mesh *	getKokkosMesh () 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 (const BoundaryID boundary_id) const

std::string	getBoundaryString (const BoundaryID boundary_id) const

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 (const unsigned int sys_num, const unsigned int var_num, const BoundaryID primary, const BoundaryID secondary)

const std::array< bool, 3 > &	queryPeriodicDimensions (const unsigned int sys_num, const unsigned int var_num) const

const std::array< bool, 3 > &	queryPeriodicDimensions (const MooseVariableBase &var) const

bool	isTranslatedPeriodic (const unsigned int sys_num, const unsigned int var_num, const unsigned int component) const

bool	isTranslatedPeriodic (const MooseVariableBase &var, const unsigned int component) const

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

RealVectorValue	minPeriodicVector (const unsigned int sys_num, const unsigned int var_num, Point p, Point q) const

RealVectorValue	minPeriodicVector (const MooseVariableBase &var, const Point &p, const Point &q) const

RealVectorValue	minPeriodicVector (const unsigned int var_num, const Point &p, const Point &q) const

Real	minPeriodicDistance (const unsigned int sys_num, const unsigned int var_num, const Point &p, const Point &q) const

Real	minPeriodicDistance (const MooseVariableBase &var, const Point &p, const Point &q) const

Real	minPeriodicDistance (const unsigned int var_num, const Point &p, const Point &q) const

void	detectPairedSidesets ()

bool	hasDetectedPairedSidesets () const

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

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

unsigned int	getMaxSidesPerElem () const

unsigned int	getMaxNodesPerElem () const

unsigned int	getMaxNodesPerSide () 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	hasLowerD () const

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

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

bool	getConstructNodeListFromSideList ()

bool	getDisplaceNodeListBySideList ()

bool	possiblyRebuildNodeToElemMap ()

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &name) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

std::string	messagePrefix (const bool hit_prefix=true) const

std::string	errorPrefix (const std::string &) const

void	mooseError (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecatedNoTrace (Args &&... args) const

void	mooseInfo (Args &&... args) const

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

PerfGraph &	perfGraph ()

const libMesh::ConstElemRange *	getActiveLocalElementRange ()

const 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	generatePinCenters (unsigned int nx, unsigned int ny, Real pitch, Real elev, std::vector< Point > &pin_centers)
	Generate pin centers.

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 unteaded lengths.

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 void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

static 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
bool	_pin_mesh_exist = false

bool	_duct_mesh_exist = false

std::vector< std::vector< Real > >	_subchannel_position
	x,y coordinates of the subchannel centroids

	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::array< bool, 3 >	periodic_dim_default

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

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
bool	allowRecovery () const

void	cacheInfo ()

void	freeBndNodes ()

void	freeBndElems ()

void	setPartitionerHelper (MeshBase *mesh=nullptr)

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) 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

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	_nx
	number of subchannels in the x direction

unsigned int	_ny
	number of subchannels in the y direction

unsigned int	_n_channels
	number of subchannels in total

unsigned int	_n_gaps
	Number of gaps per layer.

unsigned int	_n_pins
	Number of pins.

Real	_side_gap
	The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct or else the distance between the subchannel centroid to the duct wall.

std::vector< std::vector< Node * > >	_nodes
	vector of subchannel nodes

std::vector< std::vector< Node * > >	_pin_nodes
	vector of fuel pin nodes

std::vector< std::vector< Node * > >	_gapnodes
	vector of gap (interface between pairs of neighboring subchannels) nodes

std::vector< std::pair< unsigned int, unsigned int > >	_gap_to_chan_map
	map relating gap index to subchannel index

std::vector< std::pair< unsigned int, unsigned int > >	_gap_to_pin_map
	map relating gap index to fuel pin index

std::vector< std::vector< unsigned int > >	_chan_to_gap_map
	map relating subchannel index to gap index

std::vector< std::vector< unsigned int > >	_chan_to_pin_map
	map relating subchannel index to fuel pin index

std::vector< std::vector< unsigned int > >	_pin_to_chan_map
	map relating fuel pin index to subchannel index

std::vector< std::vector< double > >	_sign_id_crossflow_map
	Matrix used to give local sign to crossflow quantities.

std::vector< std::vector< Real > >	_gij_map
	Vector to store gap size.

std::vector< EChannelType >	_subch_type
	Subchannel type.

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

Real	_heated_length
	heated length of the fuel Pin

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

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

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

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

std::map< Node , Node >	_chan_to_duct_node_map
	Maps between channel nodes and duct nodes.

std::map< Node , Node >	_duct_node_to_chan_map

std::vector< Real >	_spacer_z
	axial location of the spacers

std::vector< Real >	_spacer_k
	form loss coefficient of the spacers

std::vector< Real >	_z_blockage
	axial location of blockage (inlet, outlet) [m]

std::vector< unsigned int >	_index_blockage
	index of subchannels affected by blockage

std::vector< Real >	_reduction_blockage
	area reduction of subchannels affected by blockage

Real	_kij
	Lateral form loss coefficient.

Real	_pitch
	Distance between the neighbor fuel pins, pitch.

Real	_pin_diameter
	fuel Pin diameter

Real	_assembly_flow_area
	Undeformed bundle inlet flow area.

Real	_assembly_wetted_perimeter
	Undeformed bundle inlet wetted perimeter.

Real	_assembly_hydraulic_diameter
	Undeformed bundle-average hydraulic diameter.

unsigned int	_n_cells
	number of axial cells

	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

std::unique_ptr< Moose::Kokkos::Mesh >	_kokkos_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< 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::unordered_map< dof_id_type, std::vector< dof_id_type > >	_node_to_elem_map

bool	_node_to_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::optional< std::vector< std::pair< BoundaryID, BoundaryID > > >	_paired_boundary

const bool	_is_split

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

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

Private Member Functions
std::unordered_map< dof_id_type, std::vector< dof_id_type > > &	internalNodeToElemMap ()

void	buildRefinementMap (const Elem &elem, libMesh::QBase &qrule, libMesh::QBase &qrule_face, int parent_side, int child, int child_side)

void	buildCoarseningMap (const Elem &elem, libMesh::QBase &qrule, libMesh::QBase &qrule_face, int input_side)

void	mapPoints (const std::vector< Point > &from, const std::vector< Point > &to, std::vector< QpMap > &qp_map)

void	findAdaptivityQpMaps (const Elem *template_elem, libMesh::QBase &qrule, libMesh::QBase &qrule_face, std::vector< std::vector< QpMap > > &refinement_map, std::vector< std::pair< unsigned int, QpMap > > &coarsen_map, int parent_side, int child, int child_side)

void	buildHRefinementAndCoarseningMaps (Assembly *assembly)

const std::vector< QpMap > &	getPRefinementMapHelper (const Elem &elem, const std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > > &) const

const std::vector< QpMap > &	getPCoarseningMapHelper (const Elem &elem, const std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > > &) const

void	updateCoordTransform ()

void	checkDuplicateSubdomainNames ()

void	computeMaxPerElemAndSide ()

void	buildElemIDInfo ()

void	buildLowerDMesh ()

RestartableDataValue &	registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const

void	registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)

RestartableData< T > &	declareRestartableDataHelper (const std::string &data_name, void *context, Args &&... args) const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

Private Attributes
std::unordered_map< dof_id_type, ElemInfo >	_elem_to_elem_info

std::vector< const ElemInfo * >	_elem_info

std::vector< FaceInfo >	_all_face_info

std::vector< const FaceInfo * >	_face_info

std::unordered_map< std::pair< const Elem , unsigned int >, FaceInfo >	_elem_side_to_face_info

bool	_finite_volume_info_dirty

bool	_linear_finite_volume_dofs_cached

std::map< std::pair< unsigned int, unsigned int >, std::array< bool, 3 > >	_periodic_dim

RealVectorValue	_half_range

std::vector< Node * >	_extreme_nodes

std::map< std::pair< int, libMesh::ElemType >, std::vector< std::vector< QpMap > > >	_elem_type_to_refinement_map

std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > >	_elem_type_to_p_refinement_map

std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > >	_elem_type_to_p_refinement_side_map

std::map< libMesh::ElemType, std::map< std::pair< int, int >, std::vector< std::vector< QpMap > > > >	_elem_type_to_child_side_refinement_map

std::map< std::pair< int, libMesh::ElemType >, std::vector< std::pair< unsigned int, QpMap > > >	_elem_type_to_coarsening_map

std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > >	_elem_type_to_p_coarsening_map

std::map< std::pair< libMesh::ElemType, unsigned int >, std::vector< QpMap > >	_elem_type_to_p_coarsening_side_map

std::unordered_map< SubdomainID, SubdomainData >	_sub_to_data

std::unordered_map< SubdomainID, std::set< BoundaryID > >	_neighbor_subdomain_boundary_ids

std::set< SubdomainID >	_lower_d_interior_blocks

std::set< SubdomainID >	_lower_d_boundary_blocks

std::unordered_map< std::pair< const Elem , unsigned short int >, const Elem >	_higher_d_elem_side_to_lower_d_elem

std::unordered_map< const Elem *, unsigned short int >	_lower_d_elem_to_higher_d_elem_side

bool	_allow_recovery

bool	_construct_node_list_from_side_list

bool	_displace_node_list_by_side_list

bool	_need_delete

bool	_allow_remote_element_removal

std::set< Elem * >	_ghost_elems_from_ghost_boundaries

bool	_need_ghost_ghosted_boundaries

std::vector< std::unordered_map< SubdomainID, std::set< dof_id_type > > >	_block_id_mapping

std::vector< dof_id_type >	_max_ids

std::vector< dof_id_type >	_min_ids

std::vector< std::vector< bool > >	_id_identical_flag

unsigned int	_max_sides_per_elem

unsigned int	_max_nodes_per_elem

unsigned int	_max_nodes_per_side

bool	_is_displaced

std::map< SubdomainID, Moose::CoordinateSystemType > &	_coord_sys

unsigned int	_rz_coord_axis

std::unordered_map< SubdomainID, std::pair< Point, RealVectorValue > >	_subdomain_id_to_rz_coord_axis

std::unique_ptr< MooseAppCoordTransform >	_coord_transform

bool	_coord_system_set

std::vector< SubdomainName >	_provided_coord_blocks

bool	_doing_p_refinement

unsigned int	_max_p_level

unsigned int	_max_h_level

const ParallelParamObject &	_parent

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

const RestartableDataMapName	_metaname

std::string	_restartable_name

Friends
class	SCMQuadAssemblyMeshGenerator

Detailed Description

Creates the mesh of subchannels in a quadrilateral lattice.

Definition at line 23 of file QuadSubChannelMesh.h.

Constructor & Destructor Documentation

◆ QuadSubChannelMesh() [1/2]

QuadSubChannelMesh::QuadSubChannelMesh ( const InputParameters & parameters )

Definition at line 27 of file QuadSubChannelMesh.C.

27: SubChannelMesh(params) {}

SubChannelMesh

Base class for subchannel meshes.

Definition SubChannelMesh.h:23

◆ QuadSubChannelMesh() [2/2]

QuadSubChannelMesh::QuadSubChannelMesh ( const QuadSubChannelMesh & other_mesh )

Definition at line 29 of file QuadSubChannelMesh.C.

  : SubChannelMesh(other_mesh),
    _nx(other_mesh._nx),
    _ny(other_mesh._ny),
    _n_channels(other_mesh._n_channels),
    _n_gaps(other_mesh._n_gaps),
    _n_pins(other_mesh._n_pins),
    _side_gap(other_mesh._side_gap),
    _nodes(other_mesh._nodes),
    _pin_nodes(other_mesh._pin_nodes),
    _gapnodes(other_mesh._gapnodes),
    _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),
    _chan_to_pin_map(other_mesh._chan_to_pin_map),
    _pin_to_chan_map(other_mesh._pin_to_chan_map),
    _sign_id_crossflow_map(other_mesh._sign_id_crossflow_map),
    _gij_map(other_mesh._gij_map),
    _subch_type(other_mesh._subch_type)
{
  if (_nx < 2 && _ny < 2)
    mooseError(name(),
               ": The number of subchannels cannot be less than 2 in both directions (x and y). "
               "Smallest assembly allowed is either 2X1 or 1X2. ");
}

Member Function Documentation

◆ buildMesh()

void QuadSubChannelMesh::buildMesh ( )

overridevirtual

Implements MooseMesh.

Definition at line 62 of file QuadSubChannelMesh.C.

63{

64}

◆ channelIndex()

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

overridevirtual

Implements SubChannelMesh.

Definition at line 151 of file QuadSubChannelMesh.C.

{
  // this is identical to getSubchannelIndexFromPoint, but when it is given a point "outside" the
  // normal subchannel geometry (i.e. a point that lies in a gap around the lattice) we still report
  // a valid subchannel index this is needed for transferring the solution onto a visualization mesh
 
  Real offset_x = (_nx - 1) * _pitch / 2.0;
  Real offset_y = (_ny - 1) * _pitch / 2.0;
  int i = (pt(0) + offset_x + 0.5 * _pitch) / _pitch;
  int j = (pt(1) + offset_y + 0.5 * _pitch) / _pitch;
 
  i = std::max(0, i);
  i = std::min(i, (int)(_nx - 1));
 
  j = std::max(0, j);
  j = std::min(j, (int)(_ny - 1));
 
  return j * _nx + i;
}

◆ computeAssemblyHydraulicParameters()

void QuadSubChannelMesh::computeAssemblyHydraulicParameters ( )

Compute undeformed bundle-average inlet hydraulic quantities from generated mesh geometry.

Definition at line 67 of file QuadSubChannelMesh.C.

{
  _assembly_flow_area = 0.0;
  _assembly_wetted_perimeter = 0.0;
  _assembly_hydraulic_diameter = 0.0;
 
  const Real z = _z_grid.empty() ? 0.0 : _z_grid.front();
 
  for (unsigned int i_ch = 0; i_ch < _n_channels; i_ch++)
  {
    _assembly_flow_area += getSubchannelFlowArea(i_ch, z);
    _assembly_wetted_perimeter += getSubchannelWettedPerimeter(i_ch);
  }
 
  if (_assembly_wetted_perimeter == 0.0)
    mooseError(name(), ": Assembly wetted perimeter is zero; cannot compute hydraulic diameter.");
 
  _assembly_hydraulic_diameter = 4.0 * _assembly_flow_area / _assembly_wetted_perimeter;
}

◆ ductMeshExist()

bool SubChannelMesh::ductMeshExist ( ) const

inlineinherited

Return if Duct Mesh exists or not.

Definition at line 101 of file SubChannelMesh.h.

101{ return _duct_mesh_exist; }

SubChannelMesh::_duct_mesh_exist

bool _duct_mesh_exist

Definition SubChannelMesh.h:231

◆ generatePinCenters()

void QuadSubChannelMesh::generatePinCenters	(	unsigned int	nx,
		unsigned int	ny,
		Real	pitch,
		Real	elev,
		std::vector< Point > &	pin_centers
	)

static

Generate pin centers.

Parameters

nx	number of channels in x-direction (must be more than 1)
ny	number of channels in y-direction (must be more than 1)
elev	elevation in z-direction
pin_centers	Positions in 3D space of pin centers

Definition at line 198 of file QuadSubChannelMesh.C.

{
  mooseAssert(nx >= 2, "Number of channels in x-direction must be 2 or more.");
  mooseAssert(ny >= 2, "Number of channels in y-direction must be 2 or more.");
 
  Real offset_x = (nx - 2) * pitch / 2.0;
  Real offset_y = (ny - 2) * pitch / 2.0;
  for (unsigned int iy = 0; iy < ny - 1; iy++)
    for (unsigned int ix = 0; ix < nx - 1; ix++)
      pin_centers.push_back(Point(pitch * ix - offset_x, pitch * iy - offset_y, elev));
}

Referenced by SCMDetailedQuadAssemblyMeshGenerator::generate().

◆ generateZGrid()

void SubChannelMesh::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 unteaded lengths.

Definition at line 58 of file SubChannelMesh.C.

{
  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);
}

Referenced by SCMQuadAssemblyMeshGenerator::SCMQuadAssemblyMeshGenerator(), SCMQuadDuctMeshGenerator::SCMQuadDuctMeshGenerator(), SCMTriAssemblyMeshGenerator::SCMTriAssemblyMeshGenerator(), and SCMTriDuctMeshGenerator::SCMTriDuctMeshGenerator().

◆ getAssemblyFlowArea()

Real SubChannelMesh::getAssemblyFlowArea ( ) const

inlineinherited

Return undeformed bundle inlet flow area.

Definition at line 152 of file SubChannelMesh.h.

152{ return _assembly_flow_area; }

Referenced by SubChannel1PhaseProblem::externalSolve().

◆ getAssemblyHydraulicDiameter()

Real SubChannelMesh::getAssemblyHydraulicDiameter ( ) const

inlineinherited

Return undeformed bundle-average hydraulic diameter.

Definition at line 162 of file SubChannelMesh.h.

162{ return _assembly_hydraulic_diameter; }

Referenced by SubChannel1PhaseProblem::externalSolve().

◆ getAssemblyWettedPerimeter()

Real SubChannelMesh::getAssemblyWettedPerimeter ( ) const

inlineinherited

Return undeformed bundle inlet wetted perimeter.

Definition at line 157 of file SubChannelMesh.h.

157{ return _assembly_wetted_perimeter; }

◆ getChannelGaps()

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

inlineoverridevirtual

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

Implements SubChannelMesh.

Definition at line 62 of file QuadSubChannelMesh.h.

  {
    return _chan_to_gap_map[i_chan];
  }

Referenced by QuadSubChannel1PhaseProblem::computeh(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getChannelNode()

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

inlineoverridevirtual

Get the subchannel mesh node for a given channel index and elevation index.

Implements SubChannelMesh.

Definition at line 40 of file QuadSubChannelMesh.h.

  {
    return _nodes[i_chan][iz];
  }

Referenced by QuadSubChannel1PhaseProblem::computeh(), QuadSubChannel1PhaseProblem::initializeSolution(), and QuadSubChannelNormalSliceValues::output().

◆ getChannelNodeFromDuct()

Node * SubChannelMesh::getChannelNodeFromDuct ( Node * duct_node ) const

inherited

Function that gets the channel node from the duct node.

Definition at line 104 of file SubChannelMesh.C.

{
  auto it = _duct_node_to_chan_map.find(duct_node);
  return (it == _duct_node_to_chan_map.end()) ? nullptr : it->second;
}

Referenced by SubChannel1PhaseProblem::externalSolve().

◆ getChannelPins()

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

inlineoverridevirtual

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

Implements SubChannelMesh.

Definition at line 70 of file QuadSubChannelMesh.h.

  {
    return _chan_to_pin_map[i_chan];
  }

Referenced by QuadSubChannel1PhaseProblem::computeAddedHeatPin(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getCrossflowSign()

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

inlineoverridevirtual

Return a sign for the crossflow given a subchannel index and local neighbor index.

Implements SubChannelMesh.

Definition at line 75 of file QuadSubChannelMesh.h.

  {
    return _sign_id_crossflow_map[i_chan][i_local];
  }

Referenced by QuadSubChannel1PhaseProblem::computeh().

◆ getDuctNodeFromChannel()

Node * SubChannelMesh::getDuctNodeFromChannel ( Node * channel_node ) const

inherited

Function that gets the duct node from the channel node.

Definition at line 97 of file SubChannelMesh.C.

{
  auto it = _chan_to_duct_node_map.find(channel_node);
  return (it == _chan_to_duct_node_map.end()) ? nullptr : it->second;
}

Referenced by SubChannel1PhaseProblem::computeAddedHeatDuct().

◆ getDuctNodes()

const std::vector< Node * > & SubChannelMesh::getDuctNodes ( ) const

inlineinherited

Function that returns the vector with the duct nodes.

Definition at line 223 of file SubChannelMesh.h.

223{ return _duct_nodes; }

SubChannelMesh::_duct_nodes

std::vector< Node * > _duct_nodes

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

Definition SubChannelMesh.h:248

Referenced by SubChannel1PhaseProblem::externalSolve().

◆ getGapChannels()

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

inlineoverridevirtual

Return a pair of subchannel indices for a given gap index.

Implements SubChannelMesh.

Definition at line 54 of file QuadSubChannelMesh.h.

  {
    return _gap_to_chan_map[i_gap];
  }

Referenced by QuadSubChannel1PhaseProblem::computeh().

◆ getGapPins()

const std::pair< unsigned int, unsigned int > & QuadSubChannelMesh::getGapPins ( unsigned int i_gap ) const

inlineoverridevirtual

Return a pair of pin indices for a given gap index.

Implements SubChannelMesh.

Definition at line 58 of file QuadSubChannelMesh.h.

  {
    return _gap_to_pin_map[i_gap];
  }

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

◆ getGapWidth()

Real QuadSubChannelMesh::getGapWidth	(	unsigned int	axial_index,
		unsigned int	gap_index
	)		const

inlineoverridevirtual

Return gap width for a given gap index.

Implements SubChannelMesh.

Definition at line 99 of file QuadSubChannelMesh.h.

  {
    return _gij_map[axial_index][gap_index];
  }

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

◆ getHeatedLength()

virtual const Real & SubChannelMesh::getHeatedLength ( ) const

inlinevirtualinherited

Return heated length.

Definition at line 188 of file SubChannelMesh.h.

188{ return _heated_length; }

SubChannelMesh::_heated_length

Real _heated_length

heated length of the fuel Pin

Definition SubChannelMesh.h:237

Referenced by QuadSubChannel1PhaseProblem::computeAddedHeatPin(), TriSubChannel1PhaseProblem::computeAddedHeatPin(), SCMQuadPowerAux::computeValue(), SCMTriPowerAux::computeValue(), SCMPinSurfaceTemperature::execute(), SCMPlanarMean::execute(), SCMQuadPowerAux::initialSetup(), SCMTriPowerAux::initialSetup(), SCMQuadPowerIC::initialSetup(), SCMTriPowerIC::initialSetup(), QuadSubChannelNormalSliceValues::output(), SCMQuadPowerIC::SCMQuadPowerIC(), SCMTriPowerIC::SCMTriPowerIC(), FCTFdisplacementIC::value(), SCMQuadPowerIC::value(), and SCMTriPowerIC::value().

◆ getHeatedLengthEntry()

virtual const Real & SubChannelMesh::getHeatedLengthEntry ( ) const

inlinevirtualinherited

Return unheated length at entry.

Definition at line 183 of file SubChannelMesh.h.

183{ return _unheated_length_entry; }

SubChannelMesh::_unheated_length_entry

Real _unheated_length_entry

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

Definition SubChannelMesh.h:235

Referenced by QuadSubChannel1PhaseProblem::computeAddedHeatPin(), TriSubChannel1PhaseProblem::computeAddedHeatPin(), SCMQuadPowerAux::computeValue(), SCMTriPowerAux::computeValue(), SCMPinSurfaceTemperature::execute(), SCMPlanarMean::execute(), SCMQuadPowerAux::initialSetup(), SCMTriPowerAux::initialSetup(), SCMQuadPowerIC::initialSetup(), SCMTriPowerIC::initialSetup(), QuadSubChannelNormalSliceValues::output(), FCTFdisplacementIC::value(), SCMQuadPowerIC::value(), and SCMTriPowerIC::value().

◆ getHeatedLengthExit()

virtual const Real & SubChannelMesh::getHeatedLengthExit ( ) const

inlinevirtualinherited

Return unheated length at exit.

Definition at line 193 of file SubChannelMesh.h.

193{ return _unheated_length_exit; }

SubChannelMesh::_unheated_length_exit

Real _unheated_length_exit

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

Definition SubChannelMesh.h:239

Referenced by SCMPinSurfaceTemperature::execute(), SCMPlanarMean::execute(), and QuadSubChannelNormalSliceValues::output().

◆ getIndexBlockage()

virtual const std::vector< unsigned int > & SubChannelMesh::getIndexBlockage ( ) const

inlinevirtualinherited

Get index of blocked subchannels.

Definition at line 51 of file SubChannelMesh.h.

51{ return _index_blockage; }

SubChannelMesh::_index_blockage

std::vector< unsigned int > _index_blockage

index of subchannels affected by blockage

Definition SubChannelMesh.h:262

Referenced by QuadSubChannel1PhaseProblem::initializeSolution(), TriSubChannel1PhaseProblem::initializeSolution(), and MarvelTriFlowAreaIC::value().

◆ getKGrid()

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

inlinevirtualinherited

Get axial cell location and value of loss coefficient.

Definition at line 41 of file SubChannelMesh.h.

41{ return _k_grid; }

SubChannelMesh::_k_grid

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

axial form loss coefficient per computational cell

Definition SubChannelMesh.h:244

Referenced by SubChannel1PhaseProblem::computeDP().

◆ getKij()

virtual const Real & SubChannelMesh::getKij ( ) const

inlinevirtualinherited

Return lateral loss coefficient.

Definition at line 61 of file SubChannelMesh.h.

61{ return _kij; }

SubChannelMesh::_kij

Real _kij

Lateral form loss coefficient.

Definition SubChannelMesh.h:267

◆ getNumOfAxialCells()

virtual unsigned int SubChannelMesh::getNumOfAxialCells ( ) const

inlinevirtualinherited

Return the number of axial cells.

Definition at line 66 of file SubChannelMesh.h.

66{ return processor_id() == 0 ? _n_cells : 0; }

SubChannelMesh::_n_cells

unsigned int _n_cells

number of axial cells

Definition SubChannelMesh.h:280

libMesh::ParallelObject::processor_id

processor_id_type processor_id() const

Referenced by TriSubChannel1PhaseProblem::computeh(), SCMDuctHeatRatePostprocessor::execute(), SCMPinPowerPostprocessor::execute(), SCMPinSurfaceTemperature::execute(), SCMPlanarMean::execute(), SCMTHPowerPostprocessor::execute(), SubChannelDelta::execute(), SCMQuadPowerAux::initialSetup(), SCMTriPowerAux::initialSetup(), SCMQuadPowerIC::initialSetup(), SCMTriPowerIC::initialSetup(), QuadSubChannelNormalSliceValues::output(), and SubChannel1PhaseProblem::SubChannel1PhaseProblem().

◆ getNumOfChannels()

unsigned int QuadSubChannelMesh::getNumOfChannels ( ) const

inlineoverridevirtual

Return the number of channels per layer.

Implements SubChannelMesh.

Definition at line 50 of file QuadSubChannelMesh.h.

50{ return processor_id() == 0 ? _n_channels : 0; }

Referenced by QuadSubChannelNormalSliceValues::output().

◆ getNumOfGapsPerLayer()

unsigned int QuadSubChannelMesh::getNumOfGapsPerLayer ( ) const

inlineoverridevirtual

Return the number of gaps per layer.

Implements SubChannelMesh.

Definition at line 51 of file QuadSubChannelMesh.h.

51{ return processor_id() == 0 ? _n_gaps : 0; }

◆ getNumOfPins()

unsigned int QuadSubChannelMesh::getNumOfPins ( ) const

inlineoverridevirtual

Return the number of pins.

Implements SubChannelMesh.

Definition at line 52 of file QuadSubChannelMesh.h.

52{ return processor_id() == 0 ? _n_pins : 0; }

Referenced by SCMQuadPowerAux::initialSetup(), SCMQuadPowerIC::initialSetup(), SCMQuadPowerAux::SCMQuadPowerAux(), and SCMQuadPowerIC::SCMQuadPowerIC().

◆ getNx()

const unsigned int & QuadSubChannelMesh::getNx ( ) const

inline

Number of subchannels in the -x direction.

Definition at line 81 of file QuadSubChannelMesh.h.

81{ return _nx; }

Referenced by QuadSubChannelNormalSliceValues::output(), and QuadSubChannelNormalSliceValues::QuadSubChannelNormalSliceValues().

◆ getNy()

const unsigned int & QuadSubChannelMesh::getNy ( ) const

inline

Number of subchannels in the -y direction.

Definition at line 83 of file QuadSubChannelMesh.h.

83{ return _ny; }

Referenced by QuadSubChannelNormalSliceValues::QuadSubChannelNormalSliceValues().

◆ getPinChannels()

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

inlineoverridevirtual

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

Implements SubChannelMesh.

Definition at line 66 of file QuadSubChannelMesh.h.

  {
    return _pin_to_chan_map[i_pin];
  }

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

◆ getPinDiameter()

virtual const Real & SubChannelMesh::getPinDiameter ( ) const

inlinevirtualinherited

Return undeformed Pin diameter.

Definition at line 147 of file SubChannelMesh.h.

147{ return _pin_diameter; }

SubChannelMesh::_pin_diameter

Real _pin_diameter

fuel Pin diameter

Definition SubChannelMesh.h:271

Referenced by SubChannelAddInitialConditionsAction::act(), TriSubChannel1PhaseProblem::computeh(), SCMMixingKimAndChung::computeLatticeMixingParameter(), SCMMixingChengTodreas::computeMixingParameter(), SCMHTCClosureBase::computeNusseltNumberPreInfo(), SCMFrictionUpdatedChengTodreas::computeQuadLatticeFrictionFactor(), SCMMixingChengTodreas::computeSweepFlowMixingParameter(), SCMFrictionUpdatedChengTodreas::computeTriLatticeFrictionFactor(), SubChannel1PhaseProblem::detectDeformation(), TriSubChannel1PhaseProblem::getSubChannelPeripheralDuctWidth(), TriSubChannel1PhaseProblem::initializeSolution(), SCMFrictionUpdatedChengTodreas::SCMFrictionUpdatedChengTodreas(), SCMMixingChengTodreas::SCMMixingChengTodreas(), MarvelTriFlowAreaIC::value(), and MarvelTriWettedPerimIC::value().

◆ getPinIndexFromPoint()

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

overridevirtual

Return a pin index for a given physical point p

Implements SubChannelMesh.

Definition at line 172 of file QuadSubChannelMesh.C.

{
  if (_n_pins == 0)
    mooseError(name(), ": Cannot compute a pin index because this mesh has no pins.");
 
  Real offset_x = (_nx - 2) * _pitch / 2.0;
  Real offset_y = (_ny - 2) * _pitch / 2.0;
  unsigned int i = (p(0) + offset_x) / _pitch;
  unsigned int j = (p(1) + offset_y) / _pitch;
  return j * (_nx - 1) + i;
}

Referenced by SCMQuadPowerAux::computeValue(), and SCMQuadPowerIC::value().

◆ getPinNode()

Node * QuadSubChannelMesh::getPinNode	(	unsigned int	i_pin,
		unsigned int	iz
	)		const

inlineoverridevirtual

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

Implements SubChannelMesh.

Definition at line 45 of file QuadSubChannelMesh.h.

  {
    return _pin_nodes[i_pin][iz];
  }

Referenced by QuadSubChannel1PhaseProblem::computeAddedHeatPin(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getPitch()

const Real & QuadSubChannelMesh::getPitch ( ) const

inlineoverridevirtual

Return the undeformed pitch between 2 subchannels.

Reimplemented from SubChannelMesh.

Definition at line 74 of file QuadSubChannelMesh.h.

74{ return _pitch; }

Referenced by QuadSubChannel1PhaseProblem::getSubChannelPeripheralDuctWidth(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getReductionBlockage()

virtual const std::vector< Real > & SubChannelMesh::getReductionBlockage ( ) const

inlinevirtualinherited

Get area reduction of blocked subchannels.

Definition at line 56 of file SubChannelMesh.h.

56{ return _reduction_blockage; }

SubChannelMesh::_reduction_blockage

std::vector< Real > _reduction_blockage

area reduction of subchannels affected by blockage

Definition SubChannelMesh.h:264

Referenced by QuadSubChannel1PhaseProblem::initializeSolution(), TriSubChannel1PhaseProblem::initializeSolution(), and MarvelTriFlowAreaIC::value().

◆ getSideGap()

const Real & QuadSubChannelMesh::getSideGap ( ) const

inline

Returns the side gap, not to be confused with the gap between pins, this refers to the gap next to the duct.

Edge Pitch W = (pitch/2 - pin_diameter/2 + gap) [m]

Definition at line 89 of file QuadSubChannelMesh.h.

89{ return _side_gap; }

Referenced by SCMFrictionUpdatedChengTodreas::computeQuadLatticeFrictionFactor(), QuadSubChannel1PhaseProblem::getSubChannelPeripheralDuctWidth(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getSubchannelFlowArea()

Real QuadSubChannelMesh::getSubchannelFlowArea	(	unsigned int	i_chan,
		Real	z
	)		const

overridevirtual

Return undeformed flow area for a subchannel at an axial location, including any blockage reduction.

Implements SubChannelMesh.

Definition at line 88 of file QuadSubChannelMesh.C.

{
  Real standard_area = 0.0;
  Real rod_area = 0.0;
  Real additional_area = 0.0;
 
  const auto subch_type = getSubchannelType(i_chan);
  if (subch_type == EChannelType::CORNER)
  {
    standard_area = 0.25 * _pitch * _pitch;
    rod_area = 0.25 * 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
    additional_area = _pitch * _side_gap + _side_gap * _side_gap;
  }
  else if (subch_type == EChannelType::EDGE)
  {
    standard_area = 0.5 * _pitch * _pitch;
    rod_area = 0.5 * 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
    additional_area = _pitch * _side_gap;
  }
  else
  {
    standard_area = _pitch * _pitch;
    rod_area = 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
  }
 
  Real flow_area = standard_area + additional_area - rod_area;
 
  unsigned int blockage_index = 0;
  for (const auto & i_blockage : _index_blockage)
  {
    if (i_chan == i_blockage && z >= _z_blockage.front() && z <= _z_blockage.back())
      flow_area *= _reduction_blockage[blockage_index];
    blockage_index++;
  }
 
  return flow_area;
}

Referenced by computeAssemblyHydraulicParameters(), and SCMQuadFlowAreaIC::value().

◆ getSubchannelIndexFromPoint()

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

overridevirtual

Return a subchannel index for a given physical point p

Implements SubChannelMesh.

Definition at line 141 of file QuadSubChannelMesh.C.

{
  Real offset_x = (_nx - 1) * _pitch / 2.0;
  Real offset_y = (_ny - 1) * _pitch / 2.0;
  unsigned int i = (p(0) + offset_x + 0.5 * _pitch) / _pitch;
  unsigned int j = (p(1) + offset_y + 0.5 * _pitch) / _pitch;
  return j * _nx + i;
}

Referenced by SCMQuadFlowAreaIC::value(), and SCMQuadWettedPerimIC::value().

◆ getSubchannelType()

EChannelType QuadSubChannelMesh::getSubchannelType ( unsigned int index ) const

inlineoverridevirtual

Return the type of the subchannel for given subchannel index.

Implements SubChannelMesh.

Definition at line 97 of file QuadSubChannelMesh.h.

97{ return _subch_type[index]; }

Referenced by getSubchannelFlowArea(), QuadSubChannel1PhaseProblem::getSubChannelPeripheralDuctWidth(), getSubchannelWettedPerimeter(), and QuadSubChannel1PhaseProblem::initializeSolution().

◆ getSubchannelWettedPerimeter()

Real QuadSubChannelMesh::getSubchannelWettedPerimeter ( unsigned int i_chan ) const

overridevirtual

Return undeformed wetted perimeter for a subchannel.

Implements SubChannelMesh.

Definition at line 127 of file QuadSubChannelMesh.C.

{
  const Real rod_circumference = libMesh::pi * _pin_diameter;
  const auto subch_type = getSubchannelType(i_chan);
 
  if (subch_type == EChannelType::CORNER)
    return 0.25 * rod_circumference + _pitch + 2.0 * _side_gap;
  else if (subch_type == EChannelType::EDGE)
    return 0.5 * rod_circumference + _pitch;
  else
    return rod_circumference;
}

Referenced by computeAssemblyHydraulicParameters(), and SCMQuadWettedPerimIC::value().

◆ getZBlockage()

virtual const std::vector< Real > & SubChannelMesh::getZBlockage ( ) const

inlinevirtualinherited

Get axial location of blockage (in,out) [m].

Definition at line 46 of file SubChannelMesh.h.

46{ return _z_blockage; }

Referenced by QuadSubChannel1PhaseProblem::initializeSolution(), TriSubChannel1PhaseProblem::initializeSolution(), and MarvelTriFlowAreaIC::value().

◆ getZGrid()

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

inlinevirtualinherited

Get axial location of layers.

Definition at line 31 of file SubChannelMesh.h.

31{ return _z_grid; }

Referenced by TriSubChannel1PhaseProblem::computeh(), SCMPinSurfaceTemperature::execute(), SCMPlanarMean::execute(), SCMQuadPowerAux::initialSetup(), SCMTriPowerAux::initialSetup(), SCMQuadPowerIC::initialSetup(), SCMTriPowerIC::initialSetup(), QuadSubChannelNormalSliceValues::output(), and SubChannel1PhaseProblem::SubChannel1PhaseProblem().

◆ getZIndex()

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

virtualinherited

Get axial index of point.

Definition at line 71 of file SubChannelMesh.C.

{
  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;
}

Referenced by SubChannel1PhaseProblem::externalSolve(), SCMSolutionTransfer::getFromNode(), and SubChannelMesh::setChannelToDuctMaps().

◆ pinIndex()

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

overridevirtual

Implements SubChannelMesh.

Definition at line 185 of file QuadSubChannelMesh.C.

{
  if (_n_pins == 0)
    mooseError(name(), ": Cannot compute a pin index because this mesh has no pins.");
 
  Real offset_x = (_nx - 2) * _pitch / 2.0;
  Real offset_y = (_ny - 2) * _pitch / 2.0;
  unsigned int i = (p(0) + offset_x) / _pitch;
  unsigned int j = (p(1) + offset_y) / _pitch;
  return j * (_nx - 1) + i;
}

◆ pinMeshExist()

bool SubChannelMesh::pinMeshExist ( ) const

inlineinherited

Return if Pin Mesh exists or not.

Definition at line 96 of file SubChannelMesh.h.

96{ return _pin_mesh_exist; }

SubChannelMesh::_pin_mesh_exist

bool _pin_mesh_exist

Definition SubChannelMesh.h:230

Referenced by SubChannelAddInitialConditionsAction::act(), SCMPinSurfaceTemperature::SCMPinSurfaceTemperature(), SCMQuadPowerAux::SCMQuadPowerAux(), SCMQuadPowerIC::SCMQuadPowerIC(), SCMTriPowerAux::SCMTriPowerAux(), and SCMTriPowerIC::SCMTriPowerIC().

◆ safeClone()

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

overridevirtual

Implements MooseMesh.

Definition at line 56 of file QuadSubChannelMesh.C.

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

◆ setChannelToDuctMaps()

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

inherited

Function that sets the channel-to-duct maps.

Definition at line 111 of file SubChannelMesh.C.

{
  _duct_nodes.clear();
  _chan_to_duct_node_map.clear();
  _duct_node_to_chan_map.clear();
 
  if (_z_grid.empty())
    mooseError("setChannelToDuctMaps: _z_grid is empty; cannot match duct nodes by z.");
 
  if (_subchannel_position.empty())
    mooseError(
        "setChannelToDuctMaps: _subchannel_position is empty; cannot map duct nodes to channels.");
 
  for (size_t i = 0; i < duct_nodes.size(); ++i)
  {
    Node * dn = duct_nodes[i];
 
    // 1) Find closest subchannel center in XY
    unsigned int min_chan = 0;
    Real min_dist = std::numeric_limits<Real>::max();
 
    const Point ductpos((*dn)(0), (*dn)(1), 0.0);
 
    for (unsigned int j = 0; j < _subchannel_position.size(); ++j)
    {
      const Point chanpos(_subchannel_position[j][0], _subchannel_position[j][1], 0.0);
      const Real dist = (chanpos - ductpos).norm();
 
      if (dist < min_dist)
      {
        min_dist = dist;
        min_chan = j;
      }
    }
 
    // 2) Find channel node at the same z using getZIndex + virtual accessor
    const unsigned int iz = getZIndex(*dn);
    Node * chan_node = getChannelNode(min_chan, iz);
 
    // 3) Store bidirectional mapping
    _duct_node_to_chan_map[dn] = chan_node;
    _chan_to_duct_node_map[chan_node] = dn;
  }
 
  _duct_nodes = duct_nodes;
  _duct_mesh_exist = true;
}

Referenced by SCMQuadDuctMeshGenerator::generate(), and SCMTriDuctMeshGenerator::generate().

◆ setGapWidth()

void QuadSubChannelMesh::setGapWidth	(	unsigned int	axial_index,
		unsigned int	gap_index,
		Real	gap_width
	)

inline

Set the gap width for a given axial cell and gap index.

Definition at line 107 of file QuadSubChannelMesh.h.

  {
    _gij_map[axial_index][gap_index] = gap_width;
  }

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

◆ validParams()

InputParameters QuadSubChannelMesh::validParams ( )

static

Definition at line 19 of file QuadSubChannelMesh.C.

{
  InputParameters params = SubChannelMesh::validParams();
  params.addClassDescription("Creates an subchannel mesh container for a square "
                             "lattice arrangement");
  return params;
}

Friends And Related Symbol Documentation

◆ SCMQuadAssemblyMeshGenerator

friend class SCMQuadAssemblyMeshGenerator

friend

Definition at line 170 of file QuadSubChannelMesh.h.

Member Data Documentation

◆ _assembly_flow_area

Real SubChannelMesh::_assembly_flow_area

protectedinherited

Undeformed bundle inlet flow area.

Definition at line 273 of file SubChannelMesh.h.

Referenced by computeAssemblyHydraulicParameters(), TriSubChannelMesh::computeAssemblyHydraulicParameters(), and SubChannelMesh::getAssemblyFlowArea().

◆ _assembly_hydraulic_diameter

Real SubChannelMesh::_assembly_hydraulic_diameter

protectedinherited

Undeformed bundle-average hydraulic diameter.

Definition at line 277 of file SubChannelMesh.h.

Referenced by computeAssemblyHydraulicParameters(), TriSubChannelMesh::computeAssemblyHydraulicParameters(), and SubChannelMesh::getAssemblyHydraulicDiameter().

◆ _assembly_wetted_perimeter

Real SubChannelMesh::_assembly_wetted_perimeter

protectedinherited

Undeformed bundle inlet wetted perimeter.

Definition at line 275 of file SubChannelMesh.h.

Referenced by computeAssemblyHydraulicParameters(), TriSubChannelMesh::computeAssemblyHydraulicParameters(), and SubChannelMesh::getAssemblyWettedPerimeter().

◆ _chan_to_duct_node_map

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

protectedinherited

Maps between channel nodes and duct nodes.

Definition at line 251 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getDuctNodeFromChannel(), and SubChannelMesh::setChannelToDuctMaps().

◆ _chan_to_gap_map

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

protected

map relating subchannel index to gap index

Definition at line 143 of file QuadSubChannelMesh.h.

Referenced by getChannelGaps(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _chan_to_pin_map

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

protected

map relating subchannel index to fuel pin index

Definition at line 145 of file QuadSubChannelMesh.h.

Referenced by getChannelPins(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _duct_mesh_exist

bool SubChannelMesh::_duct_mesh_exist = false

inherited

Definition at line 231 of file SubChannelMesh.h.

Referenced by SubChannelMesh::ductMeshExist(), SubChannelMesh::setChannelToDuctMaps(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _duct_node_to_chan_map

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

protectedinherited

Definition at line 252 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getChannelNodeFromDuct(), and SubChannelMesh::setChannelToDuctMaps().

◆ _duct_nodes

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

protectedinherited

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

Definition at line 248 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getDuctNodes(), and SubChannelMesh::setChannelToDuctMaps().

◆ _gap_to_chan_map

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

protected

map relating gap index to subchannel index

Definition at line 139 of file QuadSubChannelMesh.h.

Referenced by getGapChannels(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _gap_to_pin_map

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

protected

map relating gap index to fuel pin index

Definition at line 141 of file QuadSubChannelMesh.h.

Referenced by getGapPins(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _gapnodes

std::vector<std::vector<Node *> > QuadSubChannelMesh::_gapnodes

protected

vector of gap (interface between pairs of neighboring subchannels) nodes

Definition at line 136 of file QuadSubChannelMesh.h.

Referenced by SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _gij_map

std::vector<std::vector<Real> > QuadSubChannelMesh::_gij_map

protected

Vector to store gap size.

Definition at line 152 of file QuadSubChannelMesh.h.

Referenced by getGapWidth(), setGapWidth(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _heated_length

Real SubChannelMesh::_heated_length

protectedinherited

heated length of the fuel Pin

Definition at line 237 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getHeatedLength(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _index_blockage

std::vector<unsigned int> SubChannelMesh::_index_blockage

protectedinherited

index of subchannels affected by blockage

Definition at line 262 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getIndexBlockage(), getSubchannelFlowArea(), TriSubChannelMesh::getSubchannelFlowArea(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _k_grid

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

protectedinherited

axial form loss coefficient per computational cell

Definition at line 244 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getKGrid(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _kij

Real SubChannelMesh::_kij

protectedinherited

Lateral form loss coefficient.

Definition at line 267 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getKij(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _n_cells

unsigned int SubChannelMesh::_n_cells

protectedinherited

number of axial cells

Definition at line 280 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getNumOfAxialCells(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _n_channels

unsigned int QuadSubChannelMesh::_n_channels

protected

number of subchannels in total

Definition at line 118 of file QuadSubChannelMesh.h.

Referenced by computeAssemblyHydraulicParameters(), getNumOfChannels(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _n_gaps

unsigned int QuadSubChannelMesh::_n_gaps

protected

Number of gaps per layer.

Definition at line 120 of file QuadSubChannelMesh.h.

Referenced by getNumOfGapsPerLayer(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _n_pins

unsigned int QuadSubChannelMesh::_n_pins

protected

Number of pins.

Definition at line 122 of file QuadSubChannelMesh.h.

Referenced by getNumOfPins(), getPinIndexFromPoint(), pinIndex(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _nodes

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

protected

vector of subchannel nodes

Definition at line 132 of file QuadSubChannelMesh.h.

Referenced by getChannelNode(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _nx

unsigned int QuadSubChannelMesh::_nx

protected

number of subchannels in the x direction

Definition at line 114 of file QuadSubChannelMesh.h.

Referenced by channelIndex(), getNx(), getPinIndexFromPoint(), getSubchannelIndexFromPoint(), pinIndex(), QuadSubChannelMesh(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _ny

unsigned int QuadSubChannelMesh::_ny

protected

number of subchannels in the y direction

Definition at line 116 of file QuadSubChannelMesh.h.

Referenced by channelIndex(), getNy(), getPinIndexFromPoint(), getSubchannelIndexFromPoint(), pinIndex(), QuadSubChannelMesh(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _pin_diameter

Real SubChannelMesh::_pin_diameter

protectedinherited

fuel Pin diameter

Definition at line 271 of file SubChannelMesh.h.

Referenced by TriSubChannelMesh::channelIndex(), SubChannelMesh::getPinDiameter(), getSubchannelFlowArea(), TriSubChannelMesh::getSubchannelFlowArea(), getSubchannelWettedPerimeter(), TriSubChannelMesh::getSubchannelWettedPerimeter(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _pin_mesh_exist

bool SubChannelMesh::_pin_mesh_exist = false

inherited

Definition at line 230 of file SubChannelMesh.h.

Referenced by SubChannelMesh::pinMeshExist(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _pin_nodes

std::vector<std::vector<Node *> > QuadSubChannelMesh::_pin_nodes

protected

vector of fuel pin nodes

Definition at line 134 of file QuadSubChannelMesh.h.

Referenced by getPinNode(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _pin_to_chan_map

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

protected

map relating fuel pin index to subchannel index

Definition at line 147 of file QuadSubChannelMesh.h.

Referenced by getPinChannels(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _pitch

Real SubChannelMesh::_pitch

protectedinherited

Distance between the neighbor fuel pins, pitch.

Definition at line 269 of file SubChannelMesh.h.

Referenced by channelIndex(), TriSubChannelMesh::channelIndex(), getPinIndexFromPoint(), SubChannelMesh::getPitch(), getPitch(), getSubchannelFlowArea(), TriSubChannelMesh::getSubchannelFlowArea(), getSubchannelIndexFromPoint(), getSubchannelWettedPerimeter(), TriSubChannelMesh::getSubchannelWettedPerimeter(), pinIndex(), TriSubChannelMesh::pinIndex(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _reduction_blockage

std::vector<Real> SubChannelMesh::_reduction_blockage

protectedinherited

area reduction of subchannels affected by blockage

Definition at line 264 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getReductionBlockage(), getSubchannelFlowArea(), TriSubChannelMesh::getSubchannelFlowArea(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _side_gap

Real QuadSubChannelMesh::_side_gap

protected

The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct or else the distance between the subchannel centroid to the duct wall.

Edge Pitch W = (pitch/2 - pin_diameter/2 + gap) [m]

Definition at line 129 of file QuadSubChannelMesh.h.

Referenced by getSideGap(), getSubchannelFlowArea(), getSubchannelWettedPerimeter(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _sign_id_crossflow_map

std::vector<std::vector<double> > QuadSubChannelMesh::_sign_id_crossflow_map

protected

Matrix used to give local sign to crossflow quantities.

Definition at line 150 of file QuadSubChannelMesh.h.

Referenced by getCrossflowSign(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _spacer_k

std::vector<Real> SubChannelMesh::_spacer_k

protectedinherited

form loss coefficient of the spacers

Definition at line 257 of file SubChannelMesh.h.

Referenced by SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _spacer_z

std::vector<Real> SubChannelMesh::_spacer_z

protectedinherited

axial location of the spacers

Definition at line 255 of file SubChannelMesh.h.

Referenced by SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _subch_type

std::vector<EChannelType> QuadSubChannelMesh::_subch_type

protected

Subchannel type.

Definition at line 154 of file QuadSubChannelMesh.h.

Referenced by getSubchannelType(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _subchannel_position

std::vector<std::vector<Real> > SubChannelMesh::_subchannel_position

inherited

x,y coordinates of the subchannel centroids

Definition at line 284 of file SubChannelMesh.h.

Referenced by TriSubChannelMesh::channelIndex(), SubChannelMesh::setChannelToDuctMaps(), SubChannelMesh::SubChannelMesh(), SubChannelPointValue::SubChannelPointValue(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _unheated_length_entry

Real SubChannelMesh::_unheated_length_entry

protectedinherited

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

Definition at line 235 of file SubChannelMesh.h.

Referenced by SCMTriAssemblyMeshGenerator::generate(), SubChannelMesh::getHeatedLengthEntry(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _unheated_length_exit

Real SubChannelMesh::_unheated_length_exit

protectedinherited

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

Definition at line 239 of file SubChannelMesh.h.

Referenced by SubChannelMesh::getHeatedLengthExit(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _z_blockage

std::vector<Real> SubChannelMesh::_z_blockage

protectedinherited

axial location of blockage (inlet, outlet) [m]

Definition at line 260 of file SubChannelMesh.h.

Referenced by getSubchannelFlowArea(), TriSubChannelMesh::getSubchannelFlowArea(), SubChannelMesh::getZBlockage(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

◆ _z_grid

std::vector<Real> SubChannelMesh::_z_grid

protectedinherited

axial location of nodes

Definition at line 242 of file SubChannelMesh.h.

Referenced by computeAssemblyHydraulicParameters(), TriSubChannelMesh::computeAssemblyHydraulicParameters(), SubChannelMesh::getZGrid(), SubChannelMesh::getZIndex(), SubChannelMesh::setChannelToDuctMaps(), and SCMQuadAssemblyMeshGenerator::transferMetadata().

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

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

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Static Private Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ QuadSubChannelMesh() [1/2]

◆ QuadSubChannelMesh() [2/2]

Member Function Documentation

◆ buildMesh()

◆ channelIndex()

◆ computeAssemblyHydraulicParameters()

◆ ductMeshExist()

◆ generatePinCenters()

◆ generateZGrid()

◆ getAssemblyFlowArea()

◆ getAssemblyHydraulicDiameter()

◆ getAssemblyWettedPerimeter()

◆ getChannelGaps()

◆ getChannelNode()

◆ getChannelNodeFromDuct()

◆ getChannelPins()

◆ getCrossflowSign()

◆ getDuctNodeFromChannel()

◆ getDuctNodes()

◆ getGapChannels()

◆ getGapPins()

◆ getGapWidth()

◆ getHeatedLength()

◆ getHeatedLengthEntry()

◆ getHeatedLengthExit()

◆ getIndexBlockage()

◆ getKGrid()

◆ getKij()

◆ getNumOfAxialCells()

◆ getNumOfChannels()

◆ getNumOfGapsPerLayer()

◆ getNumOfPins()

◆ getNx()

◆ getNy()

◆ getPinChannels()

◆ getPinDiameter()

◆ getPinIndexFromPoint()

◆ getPinNode()

◆ getPitch()

◆ getReductionBlockage()

◆ getSideGap()

◆ getSubchannelFlowArea()

◆ getSubchannelIndexFromPoint()

◆ getSubchannelType()

◆ getSubchannelWettedPerimeter()

◆ getZBlockage()

◆ getZGrid()

◆ getZIndex()

◆ pinIndex()

◆ pinMeshExist()

◆ safeClone()

◆ setChannelToDuctMaps()

◆ setGapWidth()

◆ validParams()

Friends And Related Symbol Documentation

◆ SCMQuadAssemblyMeshGenerator

Member Data Documentation

◆ _assembly_flow_area

◆ _assembly_hydraulic_diameter

◆ _assembly_wetted_perimeter

◆ _chan_to_duct_node_map

◆ _chan_to_gap_map

◆ _chan_to_pin_map

◆ _duct_mesh_exist

◆ _duct_node_to_chan_map

◆ _duct_nodes