https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Types | Protected Member Functions | Protected Attributes | Friends | List of all members
QuadSubChannelMesh Class Reference

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

#include <QuadSubChannelMesh.h>

Inheritance diagram for QuadSubChannelMesh:
[legend]

Public Types

enum  ParallelType { ParallelType::DEFAULT, ParallelType::REPLICATED, ParallelType::DISTRIBUTED }
 
typedef std::pair< const Node *, BoundaryIDPeriodicNodeInfo
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 QuadSubChannelMesh (const InputParameters &parameters)
 
 QuadSubChannelMesh (const QuadSubChannelMesh &other_mesh)
 
std::unique_ptr< MooseMeshsafeClone () const override
 
void buildMesh () override
 
void computeAssemblyHydraulicParameters ()
 Compute undeformed bundle-average inlet hydraulic quantities from generated mesh geometry. More...
 
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. More...
 
Real getSubchannelWettedPerimeter (unsigned int i_chan) const override
 Return undeformed wetted perimeter for a subchannel. More...
 
Node * getChannelNode (unsigned int i_chan, unsigned int iz) const override
 Get the subchannel mesh node for a given channel index and elevation index. More...
 
Node * getPinNode (unsigned int i_pin, unsigned int iz) const override
 Get the pin mesh node for a given pin index and elevation index. More...
 
unsigned int getNumOfChannels () const override
 Return the number of channels per layer. More...
 
unsigned int getNumOfGapsPerLayer () const override
 Return the number of gaps per layer. More...
 
unsigned int getNumOfPins () const override
 Return the number of pins. More...
 
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. More...
 
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. More...
 
const std::vector< unsigned int > & getChannelGaps (unsigned int i_chan) const override
 Return a vector of gap indices for a given channel index. More...
 
const std::vector< unsigned int > & getPinChannels (unsigned int i_pin) const override
 Return a vector of channel indices for a given Pin index. More...
 
const std::vector< unsigned int > & getChannelPins (unsigned int i_chan) const override
 Return a vector of pin indices for a given channel index. More...
 
const RealgetPitch () const override
 Return the undeformed pitch between 2 subchannels. More...
 
const RealgetCrossflowSign (unsigned int i_chan, unsigned int i_local) const override
 Return a sign for the crossflow given a subchannel index and local neighbor index. More...
 
const unsigned intgetNx () const
 Number of subchannels in the -x direction. More...
 
const unsigned intgetNy () const
 Number of subchannels in the -y direction. More...
 
const RealgetSideGap () const
 Returns the side gap, not to be confused with the gap between pins, this refers to the gap next to the duct. More...
 
unsigned int getSubchannelIndexFromPoint (const Point &p) const override
 Return a subchannel index for a given physical point p More...
 
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 More...
 
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. More...
 
Real getGapWidth (unsigned int axial_index, unsigned int gap_index) const override
 Return gap width for a given gap index. More...
 
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. 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 std::vector< Real > & getZBlockage () const
 Get axial location of blockage (in,out) [m]. More...
 
virtual const std::vector< unsigned int > & getIndexBlockage () const
 Get index of blocked subchannels. More...
 
virtual const std::vector< Real > & getReductionBlockage () const
 Get area reduction of blocked subchannels. More...
 
virtual const RealgetKij () const
 Return lateral loss coefficient. More...
 
virtual unsigned int getNumOfAxialCells () const
 Return the number of axial cells. More...
 
Node * getChannelNodeFromDuct (Node *duct_node) const
 Function that gets the channel node from the duct node. More...
 
Node * getDuctNodeFromChannel (Node *channel_node) const
 Function that gets the duct node from the channel node. More...
 
bool pinMeshExist () const
 Return if Pin Mesh exists or not. More...
 
bool ductMeshExist () const
 Return if Duct Mesh exists or not. More...
 
virtual const RealgetPinDiameter () const
 Return undeformed Pin diameter. More...
 
Real getAssemblyFlowArea () const
 Return undeformed bundle inlet flow area. More...
 
Real getAssemblyWettedPerimeter () const
 Return undeformed bundle inlet wetted perimeter. More...
 
Real getAssemblyHydraulicDiameter () const
 Return undeformed bundle-average hydraulic diameter. More...
 
virtual const RealgetHeatedLengthEntry () const
 Return unheated length at entry. More...
 
virtual const RealgetHeatedLength () const
 Return heated length. More...
 
virtual const RealgetHeatedLengthExit () const
 Return unheated length at exit. More...
 
const std::vector< Node * > & getDuctNodes () const
 Function that returns the vector with the duct nodes. More...
 
void setChannelToDuctMaps (const std::vector< Node *> &duct_nodes)
 Function that sets the channel-to-duct maps. More...
 
virtual MooseMeshclone () const
 
void determineUseDistributedMesh ()
 
std::unique_ptr< MeshBase > buildMeshBaseObject (unsigned int dim=libMesh::invalid_uint)
 
std::unique_ptr< TbuildTypedMesh (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< BoundaryIDgetBoundaryIDs (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< BoundaryIDgetBoundaryIDs (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 ()
 
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 ()
 
ConstElemPointerRangerefinedElementRange () const
 
ConstElemPointerRangecoarsenedElementRange () 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_typegetBoundaryActiveSemiLocalElemIds (BoundaryID bid) const
 
std::unordered_set< dof_id_typegetBoundaryActiveNeighborElemIds (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 RealVectorValuegetNormalByBoundaryID (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::PatchUpdateTypegetPatchUpdateStrategy () 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::MeshgetKokkosMesh ()
 
const Moose::Kokkos::MeshgetKokkosMesh () 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< SubdomainIDgetSubdomainIDs (const std::vector< SubdomainName > &subdomain_names) const
 
std::set< SubdomainIDgetSubdomainIDs (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< BoundaryIDgetSubdomainInterfaceBoundaryIds (const SubdomainID subdomain_id) const
 
std::set< SubdomainIDgetBoundaryConnectedBlocks (const BoundaryID bid) const
 
std::set< SubdomainIDgetBoundaryConnectedSecondaryBlocks (const BoundaryID bid) const
 
std::set< SubdomainIDgetInterfaceConnectedBlocks (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 MooseEnumpartitionerName () const
 
bool isPartitionerForced () const
 
void allowRecovery (bool allow)
 
void setCustomPartitioner (libMesh::Partitioner *partitioner)
 
bool isRegularOrthogonal ()
 
bool hasSecondOrderElements ()
 
virtual std::unique_ptr< libMesh::PointLocatorBasegetPointLocator () 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_typegetAllElemIDs (unsigned int elem_id_index) const
 
std::set< dof_id_typegetElemIDsOnBlocks (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
 
MooseAppCoordTransformcoordTransform ()
 
const MooseUnitslengthUnit () 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 ()
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject () const
 
MooseAppgetMooseApp () 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 InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool haveParameter (const std::string &name) const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
PerfGraphperfGraph ()
 
libMesh::ConstElemRangegetActiveLocalElementRange ()
 
libMesh::ConstElemRangegetActiveLocalElementRange ()
 
libMesh::NodeRangegetActiveNodeRange ()
 
libMesh::NodeRangegetActiveNodeRange ()
 
SemiLocalNodeRangegetActiveSemiLocalNodeRange () const
 
SemiLocalNodeRangegetActiveSemiLocalNodeRange () const
 
libMesh::ConstNodeRangegetLocalNodeRange ()
 
libMesh::ConstNodeRangegetLocalNodeRange ()
 
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 FaceInfofaceInfo (const Elem *elem, unsigned int side) const
 
const std::vector< const FaceInfo *> & faceInfo () const
 
const FaceInfofaceInfo (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 ElemInfoelemInfo (const dof_id_type id) const
 
const ElemInfoelemInfo (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. 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 unteaded 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 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 More...
 
 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

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
 
TdeclareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< TdeclareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const TgetRestartableData (const std::string &data_name) const
 
TdeclareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
TdeclareRecoverableData (const std::string &data_name, Args &&... args)
 
TdeclareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
TdeclareRestartableDataWithObjectNameWithContext (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 More...
 
unsigned int _ny
 number of subchannels in the y direction More...
 
unsigned int _n_channels
 number of subchannels in total More...
 
unsigned int _n_gaps
 Number of gaps per layer. More...
 
unsigned int _n_pins
 Number of pins. More...
 
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. More...
 
std::vector< std::vector< Node * > > _nodes
 vector of subchannel nodes More...
 
std::vector< std::vector< Node * > > _pin_nodes
 vector of fuel pin nodes More...
 
std::vector< std::vector< Node * > > _gapnodes
 vector of gap (interface between pairs of neighboring subchannels) nodes More...
 
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_chan_map
 map relating gap index to subchannel index More...
 
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_pin_map
 map relating gap index to fuel pin index More...
 
std::vector< std::vector< unsigned int > > _chan_to_gap_map
 map relating subchannel index to gap index More...
 
std::vector< std::vector< unsigned int > > _chan_to_pin_map
 map relating subchannel index to fuel pin index More...
 
std::vector< std::vector< unsigned int > > _pin_to_chan_map
 map relating fuel pin index to subchannel index More...
 
std::vector< std::vector< double > > _sign_id_crossflow_map
 Matrix used to give local sign to crossflow quantities. More...
 
std::vector< std::vector< Real > > _gij_map
 Vector to store gap size. More...
 
std::vector< EChannelType_subch_type
 Subchannel type. 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...
 
std::vector< Node * > _duct_nodes
 A list of all mesh nodes that form the (elements of) the duct mesh that surrounds the pins/subchannels. More...
 
std::map< Node *, Node * > _chan_to_duct_node_map
 Maps between channel nodes and duct nodes. More...
 
std::map< Node *, Node * > _duct_node_to_chan_map
 
std::vector< Real_spacer_z
 axial location of the spacers More...
 
std::vector< Real_spacer_k
 form loss coefficient of the spacers More...
 
std::vector< Real_z_blockage
 axial location of blockage (inlet, outlet) [m] More...
 
std::vector< unsigned int_index_blockage
 index of subchannels affected by blockage More...
 
std::vector< Real_reduction_blockage
 area reduction of subchannels affected by blockage More...
 
Real _kij
 Lateral form loss coefficient. More...
 
Real _pitch
 Distance between the neighbor fuel pins, pitch. More...
 
Real _pin_diameter
 fuel Pin diameter More...
 
Real _assembly_flow_area
 Undeformed bundle inlet flow area. More...
 
Real _assembly_wetted_perimeter
 Undeformed bundle inlet wetted perimeter. More...
 
Real _assembly_hydraulic_diameter
 Undeformed bundle-average hydraulic diameter. More...
 
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
 
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< 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::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
 

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

◆ QuadSubChannelMesh() [2/2]

QuadSubChannelMesh::QuadSubChannelMesh ( const QuadSubChannelMesh other_mesh)

Definition at line 29 of file QuadSubChannelMesh.C.

30  : SubChannelMesh(other_mesh),
31  _nx(other_mesh._nx),
32  _ny(other_mesh._ny),
33  _n_channels(other_mesh._n_channels),
34  _n_gaps(other_mesh._n_gaps),
35  _n_pins(other_mesh._n_pins),
36  _side_gap(other_mesh._side_gap),
37  _nodes(other_mesh._nodes),
38  _pin_nodes(other_mesh._pin_nodes),
39  _gapnodes(other_mesh._gapnodes),
41  _gap_to_pin_map(other_mesh._gap_to_pin_map),
46  _gij_map(other_mesh._gij_map),
47  _subch_type(other_mesh._subch_type)
48 {
49  if (_nx < 2 && _ny < 2)
50  mooseError(name(),
51  ": The number of subchannels cannot be less than 2 in both directions (x and y). "
52  "Smallest assembly allowed is either 2X1 or 1X2. ");
53 }
std::vector< std::vector< unsigned int > > _chan_to_gap_map
map relating subchannel index to gap index
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_pin_map
map relating gap index to fuel pin index
unsigned int _n_gaps
Number of gaps per layer.
std::vector< EChannelType > _subch_type
Subchannel type.
std::vector< std::vector< unsigned int > > _pin_to_chan_map
map relating fuel pin index to subchannel index
std::vector< std::vector< unsigned int > > _chan_to_pin_map
map relating subchannel index to fuel pin index
unsigned int _ny
number of subchannels in the y direction
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 o...
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_chan_map
map relating gap index to subchannel index
const std::string & name() const
std::vector< std::vector< double > > _sign_id_crossflow_map
Matrix used to give local sign to crossflow quantities.
unsigned int _n_channels
number of subchannels in total
std::vector< std::vector< Real > > _gij_map
Vector to store gap size.
std::vector< std::vector< Node * > > _nodes
vector of subchannel nodes
std::vector< std::vector< Node * > > _pin_nodes
vector of fuel pin nodes
void mooseError(Args &&... args) const
std::vector< std::vector< Node * > > _gapnodes
vector of gap (interface between pairs of neighboring subchannels) nodes
unsigned int _nx
number of subchannels in the x direction
SubChannelMesh(const InputParameters &parameters)

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.

152 {
153  // this is identical to getSubchannelIndexFromPoint, but when it is given a point "outside" the
154  // normal subchannel geometry (i.e. a point that lies in a gap around the lattice) we still report
155  // a valid subchannel index this is needed for transferring the solution onto a visualization mesh
156 
157  Real offset_x = (_nx - 1) * _pitch / 2.0;
158  Real offset_y = (_ny - 1) * _pitch / 2.0;
159  int i = (pt(0) + offset_x + 0.5 * _pitch) / _pitch;
160  int j = (pt(1) + offset_y + 0.5 * _pitch) / _pitch;
161 
162  i = std::max(0, i);
163  i = std::min(i, (int)(_nx - 1));
164 
165  j = std::max(0, j);
166  j = std::min(j, (int)(_ny - 1));
167 
168  return j * _nx + i;
169 }
unsigned int _ny
number of subchannels in the y direction
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
unsigned int _nx
number of subchannels in the x direction

◆ computeAssemblyHydraulicParameters()

void QuadSubChannelMesh::computeAssemblyHydraulicParameters ( )

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

Definition at line 67 of file QuadSubChannelMesh.C.

68 {
69  _assembly_flow_area = 0.0;
72 
73  const Real z = _z_grid.empty() ? 0.0 : _z_grid.front();
74 
75  for (unsigned int i_ch = 0; i_ch < _n_channels; i_ch++)
76  {
79  }
80 
81  if (_assembly_wetted_perimeter == 0.0)
82  mooseError(name(), ": Assembly wetted perimeter is zero; cannot compute hydraulic diameter.");
83 
85 }
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...
std::vector< Real > _z_grid
axial location of nodes
const std::string & name() const
Real _assembly_wetted_perimeter
Undeformed bundle inlet wetted perimeter.
Real _assembly_flow_area
Undeformed bundle inlet flow area.
unsigned int _n_channels
number of subchannels in total
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void mooseError(Args &&... args) const
Real getSubchannelWettedPerimeter(unsigned int i_chan) const override
Return undeformed wetted perimeter for a subchannel.
Real _assembly_hydraulic_diameter
Undeformed bundle-average hydraulic diameter.

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

◆ generatePinCenters()

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

Generate pin centers.

Parameters
nxnumber of channels in x-direction (must be more than 1)
nynumber of channels in y-direction (must be more than 1)
elevelevation in z-direction
pin_centersPositions in 3D space of pin centers

Definition at line 198 of file QuadSubChannelMesh.C.

Referenced by SCMDetailedQuadAssemblyMeshGenerator::generate().

200 {
201  mooseAssert(nx >= 2, "Number of channels in x-direction must be 2 or more.");
202  mooseAssert(ny >= 2, "Number of channels in y-direction must be 2 or more.");
203 
204  Real offset_x = (nx - 2) * pitch / 2.0;
205  Real offset_y = (ny - 2) * pitch / 2.0;
206  for (unsigned int iy = 0; iy < ny - 1; iy++)
207  for (unsigned int ix = 0; ix < nx - 1; ix++)
208  pin_centers.push_back(Point(pitch * ix - offset_x, pitch * iy - offset_y, elev));
209 }
static const std::string pitch
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

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

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

63 {
64  Real L = unheated_length_entry + heated_length + unheated_length_exit;
65  Real dz = L / n_cells;
66  for (unsigned int i = 0; i < n_cells + 1; i++)
67  z_grid.push_back(dz * i);
68 }
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
dof_id_type n_cells

◆ getAssemblyFlowArea()

Real SubChannelMesh::getAssemblyFlowArea ( ) const
inlineinherited

Return undeformed bundle inlet flow area.

Definition at line 152 of file SubChannelMesh.h.

Referenced by SubChannel1PhaseProblem::externalSolve().

152 { return _assembly_flow_area; }
Real _assembly_flow_area
Undeformed bundle inlet flow area.

◆ getAssemblyHydraulicDiameter()

Real SubChannelMesh::getAssemblyHydraulicDiameter ( ) const
inlineinherited

Return undeformed bundle-average hydraulic diameter.

Definition at line 162 of file SubChannelMesh.h.

Referenced by SubChannel1PhaseProblem::externalSolve().

Real _assembly_hydraulic_diameter
Undeformed bundle-average hydraulic diameter.

◆ 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; }
Real _assembly_wetted_perimeter
Undeformed bundle inlet 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.

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

63  {
64  return _chan_to_gap_map[i_chan];
65  }
std::vector< std::vector< unsigned int > > _chan_to_gap_map
map relating subchannel index to gap index

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

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

41  {
42  return _nodes[i_chan][iz];
43  }
std::vector< std::vector< Node * > > _nodes
vector of subchannel nodes

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

Referenced by SubChannel1PhaseProblem::externalSolve().

105 {
106  auto it = _duct_node_to_chan_map.find(duct_node);
107  return (it == _duct_node_to_chan_map.end()) ? nullptr : it->second;
108 }
std::map< Node *, Node * > _duct_node_to_chan_map

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

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

71  {
72  return _chan_to_pin_map[i_chan];
73  }
std::vector< std::vector< unsigned int > > _chan_to_pin_map
map relating subchannel index to fuel pin index

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

Referenced by QuadSubChannel1PhaseProblem::computeh().

76  {
77  return _sign_id_crossflow_map[i_chan][i_local];
78  }
std::vector< std::vector< double > > _sign_id_crossflow_map
Matrix used to give local sign to crossflow quantities.

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

Referenced by SubChannel1PhaseProblem::computeAddedHeatDuct().

98 {
99  auto it = _chan_to_duct_node_map.find(channel_node);
100  return (it == _chan_to_duct_node_map.end()) ? nullptr : it->second;
101 }
std::map< Node *, Node * > _chan_to_duct_node_map
Maps between channel nodes and duct nodes.

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

Referenced by SubChannel1PhaseProblem::externalSolve().

223 { return _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...

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

Referenced by QuadSubChannel1PhaseProblem::computeh().

55  {
56  return _gap_to_chan_map[i_gap];
57  }
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_chan_map
map relating gap index to subchannel index

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

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

59  {
60  return _gap_to_pin_map[i_gap];
61  }
std::vector< std::pair< unsigned int, unsigned int > > _gap_to_pin_map
map relating gap index to fuel pin index

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

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

100  {
101  return _gij_map[axial_index][gap_index];
102  }
std::vector< std::vector< Real > > _gij_map
Vector to store gap size.

◆ getHeatedLength()

virtual const Real& SubChannelMesh::getHeatedLength ( ) const
inlinevirtualinherited

◆ getHeatedLengthEntry()

virtual const Real& SubChannelMesh::getHeatedLengthEntry ( ) const
inlinevirtualinherited

◆ getHeatedLengthExit()

virtual const Real& SubChannelMesh::getHeatedLengthExit ( ) const
inlinevirtualinherited

Return unheated length at exit.

Definition at line 193 of file SubChannelMesh.h.

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

193 { return _unheated_length_exit; }
Real _unheated_length_exit
unheated length of the fuel Pin at the exit of the assembly

◆ getIndexBlockage()

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

Get index of blocked subchannels.

Definition at line 51 of file SubChannelMesh.h.

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

51 { return _index_blockage; }
std::vector< unsigned int > _index_blockage
index of subchannels affected by blockage

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

Referenced by SubChannel1PhaseProblem::computeDP().

41 { return _k_grid; }
std::vector< std::vector< Real > > _k_grid
axial form loss coefficient per computational cell

◆ getKij()

virtual const Real& SubChannelMesh::getKij ( ) const
inlinevirtualinherited

Return lateral loss coefficient.

Definition at line 61 of file SubChannelMesh.h.

61 { return _kij; }
Real _kij
Lateral form loss coefficient.

◆ getNumOfAxialCells()

virtual unsigned int SubChannelMesh::getNumOfAxialCells ( ) const
inlinevirtualinherited

◆ getNumOfChannels()

unsigned int QuadSubChannelMesh::getNumOfChannels ( ) const
inlineoverridevirtual

Return the number of channels per layer.

Implements SubChannelMesh.

Definition at line 50 of file QuadSubChannelMesh.h.

Referenced by QuadSubChannelNormalSliceValues::output().

50 { return processor_id() == 0 ? _n_channels : 0; }
unsigned int _n_channels
number of subchannels in total
processor_id_type processor_id() const

◆ 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; }
unsigned int _n_gaps
Number of gaps per layer.
processor_id_type processor_id() const

◆ getNumOfPins()

unsigned int QuadSubChannelMesh::getNumOfPins ( ) const
inlineoverridevirtual

Return the number of pins.

Implements SubChannelMesh.

Definition at line 52 of file QuadSubChannelMesh.h.

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

52 { return processor_id() == 0 ? _n_pins : 0; }
unsigned int _n_pins
Number of pins.
processor_id_type processor_id() const

◆ getNx()

const unsigned int& QuadSubChannelMesh::getNx ( ) const
inline

Number of subchannels in the -x direction.

Definition at line 81 of file QuadSubChannelMesh.h.

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

81 { return _nx; }
unsigned int _nx
number of subchannels in the x direction

◆ getNy()

const unsigned int& QuadSubChannelMesh::getNy ( ) const
inline

Number of subchannels in the -y direction.

Definition at line 83 of file QuadSubChannelMesh.h.

Referenced by QuadSubChannelNormalSliceValues::QuadSubChannelNormalSliceValues().

83 { return _ny; }
unsigned int _ny
number of subchannels in the y direction

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

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

67  {
68  return _pin_to_chan_map[i_pin];
69  }
std::vector< std::vector< unsigned int > > _pin_to_chan_map
map relating fuel pin index to subchannel index

◆ getPinDiameter()

virtual const Real& SubChannelMesh::getPinDiameter ( ) const
inlinevirtualinherited

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

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

173 {
174  if (_n_pins == 0)
175  mooseError(name(), ": Cannot compute a pin index because this mesh has no pins.");
176 
177  Real offset_x = (_nx - 2) * _pitch / 2.0;
178  Real offset_y = (_ny - 2) * _pitch / 2.0;
179  unsigned int i = (p(0) + offset_x) / _pitch;
180  unsigned int j = (p(1) + offset_y) / _pitch;
181  return j * (_nx - 1) + i;
182 }
unsigned int _ny
number of subchannels in the y direction
unsigned int _n_pins
Number of pins.
const std::string & name() const
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Real p
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
unsigned int _nx
number of subchannels in the x direction

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

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

46  {
47  return _pin_nodes[i_pin][iz];
48  }
std::vector< std::vector< Node * > > _pin_nodes
vector of fuel pin nodes

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

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

74 { return _pitch; }
Real _pitch
Distance between the neighbor fuel pins, pitch.

◆ getReductionBlockage()

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

Get area reduction of blocked subchannels.

Definition at line 56 of file SubChannelMesh.h.

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

56 { return _reduction_blockage; }
std::vector< Real > _reduction_blockage
area reduction of subchannels affected by blockage

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

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

89 { return _side_gap; }
Real _side_gap
The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct o...

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

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

89 {
90  Real standard_area = 0.0;
91  Real rod_area = 0.0;
92  Real additional_area = 0.0;
93 
94  const auto subch_type = getSubchannelType(i_chan);
95  if (subch_type == EChannelType::CORNER)
96  {
97  standard_area = 0.25 * _pitch * _pitch;
98  rod_area = 0.25 * 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
99  additional_area = _pitch * _side_gap + _side_gap * _side_gap;
100  }
101  else if (subch_type == EChannelType::EDGE)
102  {
103  standard_area = 0.5 * _pitch * _pitch;
104  rod_area = 0.5 * 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
105  additional_area = _pitch * _side_gap;
106  }
107  else
108  {
109  standard_area = _pitch * _pitch;
110  rod_area = 0.25 * libMesh::pi * _pin_diameter * _pin_diameter;
111  }
112 
113  Real flow_area = standard_area + additional_area - rod_area;
114 
115  unsigned int blockage_index = 0;
116  for (const auto & i_blockage : _index_blockage)
117  {
118  if (i_chan == i_blockage && z >= _z_blockage.front() && z <= _z_blockage.back())
119  flow_area *= _reduction_blockage[blockage_index];
120  blockage_index++;
121  }
122 
123  return flow_area;
124 }
EChannelType getSubchannelType(unsigned int index) const override
Return the type of the subchannel for given subchannel index.
Real _side_gap
The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct o...
std::vector< Real > _z_blockage
axial location of blockage (inlet, outlet) [m]
Real _pin_diameter
fuel Pin diameter
std::vector< Real > _reduction_blockage
area reduction of subchannels affected by blockage
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::vector< unsigned int > _index_blockage
index of subchannels affected by blockage
const Real pi

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

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

142 {
143  Real offset_x = (_nx - 1) * _pitch / 2.0;
144  Real offset_y = (_ny - 1) * _pitch / 2.0;
145  unsigned int i = (p(0) + offset_x + 0.5 * _pitch) / _pitch;
146  unsigned int j = (p(1) + offset_y + 0.5 * _pitch) / _pitch;
147  return j * _nx + i;
148 }
unsigned int _ny
number of subchannels in the y direction
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Real p
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
unsigned int _nx
number of subchannels in the x direction

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

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

97 { return _subch_type[index]; }
std::vector< EChannelType > _subch_type
Subchannel type.

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

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

128 {
129  const Real rod_circumference = libMesh::pi * _pin_diameter;
130  const auto subch_type = getSubchannelType(i_chan);
131 
132  if (subch_type == EChannelType::CORNER)
133  return 0.25 * rod_circumference + _pitch + 2.0 * _side_gap;
134  else if (subch_type == EChannelType::EDGE)
135  return 0.5 * rod_circumference + _pitch;
136  else
137  return rod_circumference;
138 }
EChannelType getSubchannelType(unsigned int index) const override
Return the type of the subchannel for given subchannel index.
Real _side_gap
The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct o...
Real _pin_diameter
fuel Pin diameter
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Real pi

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

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

46 { return _z_blockage; }
std::vector< Real > _z_blockage
axial location of blockage (inlet, outlet) [m]

◆ getZGrid()

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

◆ getZIndex()

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

Get axial index of point.

Definition at line 71 of file SubChannelMesh.C.

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

72 {
73  if (_z_grid.size() == 0)
74  mooseError("_z_grid is empty.");
75 
76  if (point(2) <= _z_grid[0])
77  return 0;
78  if (point(2) >= _z_grid[_z_grid.size() - 1])
79  return _z_grid.size() - 1;
80 
81  unsigned int lo = 0;
82  unsigned int hi = _z_grid.size();
83  while (lo < hi)
84  {
85  unsigned int mid = (lo + hi) / 2;
86  if (std::abs(_z_grid[mid] - point(2)) < 1e-5)
87  return mid;
88  else if (_z_grid[mid] < point(2))
89  lo = mid;
90  else
91  hi = mid;
92  }
93  return lo;
94 }
std::vector< Real > _z_grid
axial location of nodes
void mooseError(Args &&... args) const

◆ pinIndex()

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

Implements SubChannelMesh.

Definition at line 185 of file QuadSubChannelMesh.C.

186 {
187  if (_n_pins == 0)
188  mooseError(name(), ": Cannot compute a pin index because this mesh has no pins.");
189 
190  Real offset_x = (_nx - 2) * _pitch / 2.0;
191  Real offset_y = (_ny - 2) * _pitch / 2.0;
192  unsigned int i = (p(0) + offset_x) / _pitch;
193  unsigned int j = (p(1) + offset_y) / _pitch;
194  return j * (_nx - 1) + i;
195 }
unsigned int _ny
number of subchannels in the y direction
unsigned int _n_pins
Number of pins.
const std::string & name() const
Real _pitch
Distance between the neighbor fuel pins, pitch.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Real p
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
unsigned int _nx
number of subchannels in the x direction

◆ pinMeshExist()

bool SubChannelMesh::pinMeshExist ( ) const
inlineinherited

◆ safeClone()

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

Implements MooseMesh.

Definition at line 56 of file QuadSubChannelMesh.C.

57 {
58  return _app.getFactory().copyConstruct(*this);
59 }
std::unique_ptr< T > copyConstruct(const T &object)
Factory & getFactory()
MooseApp & _app

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

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

112 {
113  _duct_nodes.clear();
114  _chan_to_duct_node_map.clear();
115  _duct_node_to_chan_map.clear();
116 
117  if (_z_grid.empty())
118  mooseError("setChannelToDuctMaps: _z_grid is empty; cannot match duct nodes by z.");
119 
120  if (_subchannel_position.empty())
121  mooseError(
122  "setChannelToDuctMaps: _subchannel_position is empty; cannot map duct nodes to channels.");
123 
124  for (size_t i = 0; i < duct_nodes.size(); ++i)
125  {
126  Node * dn = duct_nodes[i];
127 
128  // 1) Find closest subchannel center in XY
129  unsigned int min_chan = 0;
130  Real min_dist = std::numeric_limits<Real>::max();
131 
132  const Point ductpos((*dn)(0), (*dn)(1), 0.0);
133 
134  for (unsigned int j = 0; j < _subchannel_position.size(); ++j)
135  {
136  const Point chanpos(_subchannel_position[j][0], _subchannel_position[j][1], 0.0);
137  const Real dist = (chanpos - ductpos).norm();
138 
139  if (dist < min_dist)
140  {
141  min_dist = dist;
142  min_chan = j;
143  }
144  }
145 
146  // 2) Find channel node at the same z using getZIndex + virtual accessor
147  const unsigned int iz = getZIndex(*dn);
148  Node * chan_node = getChannelNode(min_chan, iz);
149 
150  // 3) Store bidirectional mapping
151  _duct_node_to_chan_map[dn] = chan_node;
152  _chan_to_duct_node_map[chan_node] = dn;
153  }
154 
155  _duct_nodes = duct_nodes;
156  _duct_mesh_exist = true;
157 }
std::vector< Real > _z_grid
axial location of nodes
std::map< Node *, Node * > _duct_node_to_chan_map
std::vector< std::vector< Real > > _subchannel_position
x,y coordinates of the subchannel centroids
std::vector< Node * > _duct_nodes
A list of all mesh nodes that form the (elements of) the duct mesh that surrounds the pins/subchannel...
virtual Node * getChannelNode(unsigned int i_chan, unsigned int iz) const =0
Get the subchannel mesh node for a given channel index and elevation index.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
std::map< Node *, Node * > _chan_to_duct_node_map
Maps between channel nodes and duct nodes.
auto norm(const T &a)
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
virtual unsigned int getZIndex(const Point &point) const
Get axial index of point.

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

Referenced by QuadSubChannel1PhaseProblem::initializeSolution().

108  {
109  _gij_map[axial_index][gap_index] = gap_width;
110  }
std::vector< std::vector< Real > > _gij_map
Vector to store gap size.

◆ validParams()

InputParameters QuadSubChannelMesh::validParams ( )
static

Definition at line 19 of file QuadSubChannelMesh.C.

20 {
22  params.addClassDescription("Creates an subchannel mesh container for a square "
23  "lattice arrangement");
24  return params;
25 }
static InputParameters validParams()
void addClassDescription(const std::string &doc_string)

Friends And Related Function 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

◆ _assembly_hydraulic_diameter

Real SubChannelMesh::_assembly_hydraulic_diameter
protectedinherited

◆ _assembly_wetted_perimeter

Real SubChannelMesh::_assembly_wetted_perimeter
protectedinherited

◆ _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

◆ _duct_node_to_chan_map

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

◆ _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

◆ _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

◆ _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

◆ _ny

unsigned int QuadSubChannelMesh::_ny
protected

◆ _pin_diameter

Real SubChannelMesh::_pin_diameter
protectedinherited

◆ _pin_mesh_exist

bool SubChannelMesh::_pin_mesh_exist = false
inherited

◆ _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

◆ _reduction_blockage

std::vector<Real> SubChannelMesh::_reduction_blockage
protectedinherited

◆ _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

◆ _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

◆ _z_grid

std::vector<Real> SubChannelMesh::_z_grid
protectedinherited

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