armi.reactor.composites module

This module contains the basic composite pattern underlying the reactor package.

This follows the principles of the Composite Design Pattern to allow the construction of a part/whole hierarchy representing a physical nuclear reactor. The composite objects act somewhat like lists: they can be indexed, iterated over, appended, extended, inserted, etc. Each member of the hierarchy knows its children and its parent, so full access to the hierarchy is available from everywhere. This design was chosen because of the close analogy of the model to the physical nature of nuclear reactors.

Warning

Because each member of the hierarchy is linked to the entire tree, it is often unsafe to save references to individual members; it can cause large and unexpected memory inefficiencies.

See Also: Developer Docs.

class armi.reactor.composites.FlagSerializer

Bases: armi.reactor.parameters.parameterDefinitions.Serializer

Serializer implementation for Flags.

This operates by converting each set of Flags (too large to fit in a uint64) into a sequence of enough uint8 elements to represent all flags. These constitute a dimension of a 2-D numpy array containing all Flags for all objects provided to the pack() function.

version: Optional[str] = '1'

static pack(data)

Flags are represented as a 2-D numpy array of uint8 (single-byte, unsigned integers), where each row contains the bytes representing a single Flags instance. We also store the list of field names so that we can verify that the reader and the writer can agree on the meaning of each bit.

Under the hood, this calls the private implementation providing the armi.reactor.flags.Flags class as the target output class.

classmethod unpack(data, version, attrs)

Reverse the pack operation.

This will allow for some degree of conversion from old flags to a new set of flags, as long as all of the source flags still exist in the current set of flags.

Under the hood, this calls the private implementation providing the armi.reactor.flags.Flags class as the target output class.

class armi.reactor.composites.CompositeModelType(name, bases, attrs)

Bases: armi.reactor.parameters.resolveCollections.ResolveParametersMeta

Metaclass for tracking subclasses of ArmiObject subclasses.

It is often useful to have an easily-accessible collection of all classes that participate in the ARMI composite reactor model. This metaclass maintains a collection of all defined subclasses, called TYPES.

TYPES: Dict[str, Type] = {'ArmiObject': <class 'armi.reactor.composites.ArmiObject'>, 'Assembly': <class 'armi.reactor.assemblies.Assembly'>, 'AssemblyList': <class 'armi.reactor.assemblyLists.AssemblyList'>, 'Block': <class 'armi.reactor.blocks.Block'>, 'CartesianAssembly': <class 'armi.reactor.assemblies.CartesianAssembly'>, 'CartesianBlock': <class 'armi.reactor.blocks.CartesianBlock'>, 'Circle': <class 'armi.reactor.components.basicShapes.Circle'>, 'Component': <class 'armi.reactor.components.component.Component'>, 'Composite': <class 'armi.reactor.composites.Composite'>, 'Core': <class 'armi.reactor.reactors.Core'>, 'Cube': <class 'armi.reactor.components.volumetricShapes.Cube'>, 'DerivedShape': <class 'armi.reactor.components.DerivedShape'>, 'DifferentialRadialSegment': <class 'armi.reactor.components.volumetricShapes.DifferentialRadialSegment'>, 'Helix': <class 'armi.reactor.components.complexShapes.Helix'>, 'HexAssembly': <class 'armi.reactor.assemblies.HexAssembly'>, 'HexBlock': <class 'armi.reactor.blocks.HexBlock'>, 'HexHoledCircle': <class 'armi.reactor.components.complexShapes.HexHoledCircle'>, 'Hexagon': <class 'armi.reactor.components.basicShapes.Hexagon'>, 'HoledHexagon': <class 'armi.reactor.components.complexShapes.HoledHexagon'>, 'HoledRectangle': <class 'armi.reactor.components.complexShapes.HoledRectangle'>, 'HoledSquare': <class 'armi.reactor.components.complexShapes.HoledSquare'>, 'MockBlock': <class 'armi.nuclearDataIO.tests.test_xsCollections.MockBlock'>, 'NullComponent': <class 'armi.reactor.components.NullComponent'>, 'PositiveOrNegativeVolumeComponent': <class 'armi.reactor.components.PositiveOrNegativeVolumeComponent'>, 'RZAssembly': <class 'armi.reactor.assemblies.RZAssembly'>, 'RadialSegment': <class 'armi.reactor.components.volumetricShapes.RadialSegment'>, 'Reactor': <class 'armi.reactor.reactors.Reactor'>, 'Rectangle': <class 'armi.reactor.components.basicShapes.Rectangle'>, 'ShapedComponent': <class 'armi.reactor.components.component.ShapedComponent'>, 'SolidRectangle': <class 'armi.reactor.components.basicShapes.SolidRectangle'>, 'SpentFuelPool': <class 'armi.reactor.assemblyLists.SpentFuelPool'>, 'Sphere': <class 'armi.reactor.components.volumetricShapes.Sphere'>, 'Square': <class 'armi.reactor.components.basicShapes.Square'>, 'ThRZAssembly': <class 'armi.reactor.assemblies.ThRZAssembly'>, 'ThRZBlock': <class 'armi.reactor.blocks.ThRZBlock'>, 'Torus': <class 'armi.reactor.components.volumetricShapes.Torus'>, 'Triangle': <class 'armi.reactor.components.basicShapes.Triangle'>, 'UnshapedComponent': <class 'armi.reactor.components.UnshapedComponent'>, 'UnshapedVolumetricComponent': <class 'armi.reactor.components.UnshapedVolumetricComponent'>, 'ZeroMassComponent': <class 'armi.reactor.components.ZeroMassComponent'>}

class armi.reactor.composites.ArmiObject(name)

Bases: object

The abstract base class for all composites and leaves.

This:

• declares the interface for objects in the composition

• implements default behavior for the interface common to all classes

• Declares an interface for accessing and managing child objects

• Defines an interface for accessing parents.

Called “component” in gang of four, this is an ArmiObject here because the word component was already taken in ARMI.

The armi.reactor.parameters.ResolveParametersMeta metaclass is used to automatically create ParameterCollection subclasses for storing parameters associated with any particular subclass of ArmiObject. Defining a pDefs class attribute in the definition of a subclass of ArmiObject will lead to the creation of a new subclass of py:class:armi.reactor.parameters.ParameterCollection, which will contain the definitions from that class’s pDefs as well as the definitions for all of its parents. A new paramCollectionType class attribute will be added to the ArmiObject subclass to reflect which type of parameter collection should be used.

Warning

This class has far too many public methods. We are in the midst of a composite tree cleanup that will likely break these out onto a number of separate functional classes grouping things like composition, location, shape/dimensions, and various physics queries. Methods are being collected here from the various specialized subclasses (Block, Assembly) in preparation for this next step. As a result, the public API on this method should be considered unstable.

Variables

• ~ArmiObject.name (str) – Object name

• ~ArmiObject.parent (ArmiObject) – The object’s parent in a hierarchical tree

• ~ArmiObject.cached (dict) – Some cached values for performance

• ~ArmiObject.p (ParameterCollection) – The state variables

• ~ArmiObject.spatialGrid (grids.Grid) – The spatial grid that this object contains

• ~ArmiObject.spatialLocator (grids.LocationBase) – The location of this object in its parent grid, or global space

See also

armi.reactor.parameters

paramCollectionType

alias of armi.reactor.parameters.parameterCollections.ArmiObjectParameterCollection

pDefs = <armi.reactor.parameters.parameterDefinitions.ParameterDefinitionCollection object>

duplicate()

Make a clean copy of this object.

Warning

Be careful with inter-object dependencies. If one object contains a reference to another object which contains links to the entire hierarchical tree, memory can fill up rather rapidly. Weak references are designed to help with this problem.

clearCache()

Clear the cache so all new values are recomputed.

copyParamsFrom(other)

Overwrite this object’s params with other object’s.

Parameters

other (ArmiObject) – The object to copy params from

updateParamsFrom(new)

Update this object’s params with a new object’s.

Parameters

new (ArmiObject) – The object to copy params from

getChildren(deep=False, generationNum=1, includeMaterials=False)

Return the children of this object.

getChildrenWithFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exactMatch=True)

Get all children that have given flags.

getComponents(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Return all armi.reactor.component.Component within this Composite.

Parameters

• typeSpec (TypeSpec) – Component flags. Will restrict Components to specific ones matching the flags specified.

• exact (bool, optional) – Only match exact component labels (names). If True, ‘coolant’ will not match ‘interCoolant’. This has no impact if compLabel is None.

Returns

items matching compLabel and exact criteria

Return type

list of Component

iterComponents(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Yield components one by one in a generator.

doChildrenHaveFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], deep=False)

Generator that yields True if the next child has given flags.

Parameters

typeSpec (TypeSpec) – Requested type of the child

containsAtLeastOneChildWithFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]])

Return True if any of the children are of a given type.

Parameters

typeSpec (TypeSpec) – Requested type of the children

See also

self.doChildrenHaveFlags, self.containsOnlyChildrenWithFlags

containsOnlyChildrenWithFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]])

Return True if all of the children are of a given type.

Parameters

typeSpec (TypeSpec) – Requested type of the children

See also

self.doChildrenHaveFlags, self.containsAtLeastOneChildWithFlags

copyParamsToChildren(paramNames)

Copy param values in paramNames to all children.

Parameters

paramNames (list) – List of param names to copy to children

classmethod getParameterCollection()

Return a new instance of the specific ParameterCollection type associated with this object.

This has the same effect as obj.paramCollectionType(). Getting a new instance through a class method like this is useful in situations where the paramCollectionType is not a top-level object and therefore cannot be trivially pickled. Since we know that by the time we want to make any instances of/unpickle a given ArmiObject, such a class attribute will have been created and associated. So, we use this top-level class method to dig dynamically down to the underlying parameter collection type.

See also

armi.reactor.parameters.parameterCollections.ParameterCollection.__reduce__()

getParamNames()

Get a list of parameters keys that are available on this object.

Will not have any corner, edge, or timenode dependence.

nameContains(s)

True if s is in this object’s name (eg. nameContains(‘fuel’)==True for ‘testfuel’.

Notes

Case insensitive (all gets converted to lower)

getName()

setName(name)

hasFlags(typeID: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exact=False)

Determine if this object is of a certain type.

Parameters

• typeID (TypeSpec) – Flags to test the object against, to see if it contains them. If a list is provided, each element is treated as a “candidate” set of flags. Return True if any of candidates match. When exact is True, the object must match one of the candidates exactly. If exact is False, the object must have at least the flags contained in a candidate for that candidate to be a match; extra flags on the object are permitted. None matches all objects if exact is False, or no objects if exact is True.

• exact (bool, optional) – Require the type of the object to fully match the provided typeID(s)

Returns

hasFlags – True if this object is in the typeID list.

Return type

bool

Notes

Type comparisons use bitwise comparisons using valid flags.

If you have an ‘inner control’ assembly, then this will evaluate True for the INNER | CONTROL flag combination. If you just want all FUEL, simply use FUEL with no additional qualifiers. For more complex comparisons, use bitwise operations.

Always returns true if typeID is none and exact is False, allowing for default parameters to be passed in when the method does not care about the object type. If the typeID is none and exact is True, this will always return False.

Examples

If you have an object with the INNER, DRIVER, and FUEL flags, then

>>> obj.getType()
[some integer]

>>> obj.hasFlags(Flags.FUEL)
True

>>> obj.hasFlags(Flags.INNER | Flags.DRIVER | Flags.FUEL)
True

>>> obj.hasFlags(Flags.OUTER | Flags.DRIVER | Flags.FUEL)
False

>>> obj.hasFlags(Flags.INNER | Flags.FUEL)
True

>>> obj.hasFlags(Flags.INNER | Flags.FUEL, exact=True)
False

>>> obj.hasFlags([Flags.INNER | Flags.DRIVER | Flags.FUEL,
... Flags.OUTER | Flags.DRIVER | Flags.FUEL], exact=True)
False

getType()

Return the object type.

setType(typ, flags: Optional[armi.reactor.flags.Flags] = None)

Set the object type.

Parameters

• typ (str) – The desired “type” for the object. Type describes the general class of the object, and typically corresponds to the name of the blueprint that created it.

• flags (Flags, optional) – The set of Flags to apply to the object. If these are omitted, then Flags will be derived from the typ.

Warning

We are in the process of developing more robust definitions for things like “name” and “type”. “type” will generally refer to the name of the blueprint that created a particular object. When present, a “name” will refer to a specific instance of an object of a particular “type”. Think unique names for each assembly in a core, even if they are all created from the same blueprint and therefore have the same “type”. When this work is complete, it will be strongly discouraged, or even disallowed to change the type of an object after it has been created, and setType() may be removed entirely.

getVolume()

getArea(cold=False)

getVolumeFractions()

Return volume fractions of each child.

Sets volume or area of missing piece (like coolant) if it exists. Caching would be nice here.

Returns

fracs – list of (component, volFrac) tuples

Return type

list

See also

test_block.Block_TestCase.test_consistentAreaWithOverlappingComponents

Notes

void areas can be negative in gaps between fuel/clad/liner(s), but these negative areas are intended to account for overlapping positive areas to insure the total area of components inside the clad is accurate. See test_block.Block_TestCase.test_consistentAreaWithOverlappingComponents

getVolumeFraction()

Return the volume fraction that this object takes up in its parent.

getMaxArea()

” The maximum area of this object if it were totally full.

See also

armi.reactor.blocks.HexBlock.getMaxArea

getMaxVolume()

The maximum volume of this object if it were totally full.

Returns

vol – volume in cm^3.

Return type

float

getMass(nuclideNames=None)

Determine the mass in grams of nuclide(s) and/or elements in this object.

Parameters

nuclideNames (str, optional) – The nuclide/element specifier to get the mass of in the object. If omitted, total mass is returned.

Returns

mass – The mass in grams.

Return type

float

getMassFrac(nucName)

Get the mass fraction of a nuclide.

Notes

If you need multiple mass fractions, use getMassFracs.

getMicroSuffix()

getMassFracs()

Get mass fractions of all nuclides in object.

Ni [1/cm3] * Ai [g/mole] ~ mass

setMassFrac(nucName, val)

Adjust the composition of this object so the mass fraction of nucName is val.

See also

setMassFracs

efficiently set multiple mass fractions at the same time.

setMassFracs(massFracs)

Apply one or more adjusted mass fractions.

This will adjust the total mass of the object, as the mass of everything designated will change, while anything else will not.

\[ \begin{align}\begin{aligned}m_i = \frac{M_i}{\sum_j(M_j)}\\(M_{j \ne i} + M_i) m_i = M_i\\\frac{m_i M_{j \ne i}}{1-m_i} = M_i\\\frac{m_i M_{j \ne i}}{V(1-m_i)} = M_i/V = m_i \rho\\N_i = \frac{m_i \rho N_A}{A_i}\\N_i = \frac{m_i M_{j \ne i} N_A}{V (1-m_i) {A_i}}\\\frac{M_{j \ne i}}{V} = m_{j \ne i} \rho\\m_{j \ne i} = 1 - m_i\end{aligned}\end{align} \]

Notes

You can’t just change one mass fraction though, you have scale all others to fill the remaining frac.

Parameters

massFracs (dict) – nucName : new mass fraction pairs.

adjustMassFrac(nuclideToAdjust=None, elementToAdjust=None, nuclideToHoldConstant=None, elementToHoldConstant=None, val=0.0)

Set the initial Zr mass fraction while maintaining Uranium enrichment, but general purpose.

Parameters

• nuclideToAdjust (str, optional) – The nuclide name to adjust

• elementToAdjust (str, optional) – The element to adjust. All isotopes in this element will adjust

• nuclideToHoldconstant (str, optional) – A nuclide to hold constant

• elementToHoldConstant (str) – Same

• val (float) – The value to set the adjust mass fraction to be.

Notes

If you use this for two elements one after the other, you will probably get something wrong. For instance, if you have U-10Zr and add Pu at 10% mass fraction, the Zr fraction will drop below 10% of the total. The U-Zr fractions will remain constant though. So this is mostly useful if you have U-10Zr and want to change it to U-5Zr.

Theory:

Mass fraction of each nuclide to be adjusted = Ai where A1+A2+A…+AI = A Mass fraction of nuclides to be held constant = Ci where sum = C Mass fraction of other nuclides is Oi, sum = O new value for A is v

A+C+O = 1.0 A’=v. If A>0, then A’=A*f1=v where f1 = v/A If A=0, then Ai’ = v/len(A), distributing the value evenly among isotopes

Now, to adjust the other nuclides, we know A’+C+O’ = 1.0 , or v+C+O’ = 1.0 So, O’= 1.0-v-C We can scale each Oi evenly by multiplying by the factor f2 Oi’ = Oi * (1-C-v)/O = Oi * f2 where f2= (1-C-v)

See also

setMassFrac, getMassFrac

adjustMassEnrichment(massFraction)

Adjust the enrichment of this object.

If it’s Uranium, enrichment means U-235 fraction. If it’s Boron, enrichment means B-10 fraction, etc.

Parameters

newEnrich (float) – The new enrichment as a fraction.

getNumberDensity(nucName)

Return the number density of a nuclide in atoms/barn-cm.

Notes

This can get called very frequently and has to do volume computations so should use some kind of caching that is invalidated by any temperature, composition, etc. changes. Even with caching the volume calls are still somewhat expensive so prefer the methods in see also.

See also

ArmiObject.getNuclideNumberDensities

More efficient for >1 specific nuc density is needed.

ArmiObject.getNumberDensities

More efficient for when all nucs in object is needed.

getNuclideNumberDensities(nucNames)

Return a list of number densities in atoms/barn-cm for the nuc names requested.

getNumberDensities(expandFissionProducts=False)

Retrieve the number densities in atoms/barn-cm of all nuclides (or those requested) in the object.

Parameters

expandFissionProducts (bool (optional)) – expand the fission product number densities

Returns

numberDensities – nucName keys, number density values (atoms/bn-cm)

Return type

dict

getNeutronEnergyDepositionConstants()

Get the neutron energy deposition group constants for a composite.

Returns

energyDepConstants – Neutron energy generation group constants (in Joules/cm)

Return type

numpy.array

Raises

RuntimeError: – Reports if a cross section library is not assigned to a reactor.

getGammaEnergyDepositionConstants()

Get the gamma energy deposition group constants for a composite.

Returns

energyDepConstants – Energy generation group constants (in Joules/cm)

Return type

numpy.array

Raises

RuntimeError: – Reports if a cross section library is not assigned to a reactor.

getTotalEnergyGenerationConstants()

Get the total energy generation group constants for a composite.

Gives the total energy generation rates when multiplied by the multigroup flux.

Returns

totalEnergyGenConstant – Total (fission + capture) energy generation group constants (Joules/cm)

Return type

numpy.array

getFissionEnergyGenerationConstants()

Get the fission energy generation group constants for a composite.

Gives the fission energy generation rates when multiplied by the multigroup flux.

Returns

fissionEnergyGenConstant – Energy generation group constants (Joules/cm)

Return type

numpy.array

Raises

RuntimeError: – Reports if a cross section library is not assigned to a reactor.

getCaptureEnergyGenerationConstants()

Get the capture energy generation group constants for a composite.

Gives the capture energy generation rates when multiplied by the multigroup flux.

Returns

fissionEnergyGenConstant – Energy generation group constants (Joules/cm)

Return type

numpy.array

Raises

RuntimeError: – Reports if a cross section library is not assigned to a reactor.

getChildrenWithNuclides(nucNames)

Return children that contain any nuclides in nucNames.

getAncestor(fn)

Return the first ancestor that satisfies the supplied predicate.

Parameters

fn (Function-like object) – The predicate used to test the validity of an ancestor. Should return true if the ancestor satisfies the caller’s requirements

getAncestorAndDistance(fn, _distance=0) → Optional[Tuple[armi.reactor.composites.ArmiObject, int]]

Return the first ancestor that satisfies the supplied predicate, along with how many levels above self the ancestor lies.

Parameters

fn (Function-like object) – The predicate used to test the validity of an ancestor. Should return true if the ancestor satisfies the caller’s requirements

getAncestorWithFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exactMatch=False)

Return the first ancestor that matches the passed flags.

Parameters

• typeSpec (TypeSpec) – A collection of flags to match on candidate parents

• exactMatch (bool) – Whether the flags match should be exact

Returns

the first ancestor up the chain of parents that matches the passed flags

Return type

armi.composites.ArmiObject

See also

ArmiObject.hasFlags

getTotalNDens()

Return the total number density of all atoms in this object.

Returns

nTot – Total ndens of all nuclides in atoms/bn-cm. Not homogenized.

Return type

float

setNumberDensity(nucName, val)

Set the number density of this nuclide to this value.

This distributes atom density evenly across all children that contain nucName. If the nuclide doesn’t exist in any of the children, then that’s actually an error. This would only happen if some unnatural nuclide like Pu239 built up in fresh UZr. That should be anticipated and dealt with elsewhere.

setNumberDensities(numberDensities)

Set one or more multiple number densities. Reset any non-listed nuclides to 0.0.

Parameters

numberDensities (dict) – nucName: ndens pairs.

Notes

We’d like to not have to call setNumberDensity for each nuclide because we don’t want to call getVolumeFractions for each nuclide (it’s inefficient).

updateNumberDensities(numberDensities)

Set one or more multiple number densities. Leaves unlisted number densities alone.

This changes a nuclide number density only on children that already have that nuclide, thereby allowing, for example, actinides to stay in the fuel component when setting block-level values.

The complication is that various number densities are distributed among various components. This sets the number density for each nuclide evenly across all components that contain it.

Parameters

numberDensities (dict) – nucName: ndens pairs.

changeNDensByFactor(factor)

Change the number density of all nuclides within the object by a multiplicative factor.

clearNumberDensities()

Reset all the number densities to nearly zero.

Set to almost zero, so components remember which nuclides are where.

density()

Returns the mass density of the object in g/cc.

getNumberOfAtoms(nucName)

Return the number of atoms of nucName in this object.

getLumpedFissionProductCollection()

Get collection of LFP objects. Will work for global or block-level LFP models.

Returns

lfps – lfpName keys , lfp object values

Return type

LumpedFissionProduct

See also

armi.physics.neutronics.fissionProductModel.lumpedFissionProduct.LumpedFissionProduct

LFP object

setLumpedFissionProducts(lfpCollection)

setChildrenLumpedFissionProducts(lfpCollection)

getFissileMassEnrich()

Returns the fissile mass enrichment.

getBoronMassEnrich()

Return B-10 mass fraction.

getUraniumMassEnrich()

Returns U-235 mass fraction assuming U-235 and U-238 only.

getUraniumNumEnrich()

Returns U-235 number fraction.

getPuN()

Returns total number density of Pu isotopes.

getPuMoles()

Returns total number of moles of Pu isotopes.

calcTotalParam(param, objs=None, volumeIntegrated=False, addSymmetricPositions=False, typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, generationNum=1, calcBasedOnFullObj=False)

Sums up a parameter throughout the object’s children or list of objects.

Parameters

• param (str) – Name of the block parameter to sum

• objs (iterable, optional) – A list of objects to sum over. If none, all children in object will be used

• volumeIntegrated (bool, optional) – Integrate over volume

• addSymmetricPositions (bool, optional) – If True, will multiply by the symmetry factor of the core (3 for 1/3 models, 1 for full core models)

• typeSpec (TypeSpec) – object types to restrict to

• generationNum (int, optional) – Which generation to consider. 1 means direct children, 2 means children of children. Default: Just return direct children.

• calcBasedOnFullObj (bool, optional) – Some assemblies or blocks, such as the center assembly in a third core model, are not modeled as full assemblies or blocks. In the third core model objects at these postions are modeled as having 1/3 the volume and thus 1/3 the power. Setting this argument to True will apply the full value of the parameter as if it was a full block or assembly.

calcAvgParam(param, typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, weightingParam=None, volumeAveraged=True, absolute=True, generationNum=1)

Calculate the child-wide average of a parameter.

Parameters

• param (str) – The ARMI block parameter that you want the average from

• typeSpec (TypeSpec) – The child types that should be included in the calculation. Restrict average to a certain child type with this parameter.

• weightingParam (None or str, optional) – An optional block param that the average will be weighted against

• volumeAveraged (bool, optional) – volume (or height, or area) average this param

• absolute (bool, optional) – Returns the average of the absolute value of param

• generationNum (int, optional) – Which generation to average over (1 for children, 2 for grandchildren)

The weighted sum is:

\[\left<\text{x}\right> = \frac{\sum_{i} x_i w_i}{\sum_i w_i}\]

where \(i\) is each child, \(x_i\) is the param value of the i-th child, and \(w_i\) is the weighting param value of the i-th child.

Warning

If a param is unset/zero on any of the children, this will be included in the average and may significantly perturb results.

Returns

The average parameter value.

Return type

float

getMaxParam(param, typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, absolute=True, generationNum=1, returnObj=False)

Find the maximum value for the parameter in this container.

Parameters

• param (str) – block parameter that will be sought.

• typeSpec (TypeSpec) – restricts the search to cover a variety of block types.

• absolute (bool) – looks for the largest magnitude value, regardless of sign, default: true

• returnObj (bool, optional) – If true, returns the child object as well as the value.

Returns

• maxVal (float) – The maximum value of the parameter asked for

• obj (child object) – The object that has the max (only returned if returnObj==True)

getMinParam(param, typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, absolute=True, generationNum=1, returnObj=False)

Find the minimum value for the parameter in this container.

See also

getMaxParam

details

getChildParamValues(param)

Get the child parameter values in a numpy array.

isFuel()

True if this is a fuel block.

containsHeavyMetal()

True if this has HM.

getNuclides()

Determine which nuclides are present in this armi object.

Returns

List of nuclide names that exist in this

Return type

list

getFissileMass()

Returns fissile mass in grams.

getHMMass()

Returns heavy metal mass in grams.

getHMMoles()

Get the number of moles of heavy metal in this object in full symmetry.

Notes

If an object is on a symmetry line, the number of moles will be scaled up by the symmetry factor. This is done because this is typically used for tracking burnup, and BOL moles are computed in full objects too so there are no complications as things move on and off of symmetry lines.

Warning

getHMMoles is different than every other get mass call since it multiplies by symmetry factor but getVolume() on the block level divides by symmetry factor causing them to cancel out.

This was needed so that HM moles mass did not change based on if the block/assembly was on a symmetry line or not.

getHMDens()

Compute the total heavy metal density of this object.

Returns

hmDens – The total heavy metal number (atom) density in atoms/bn-cm.

Return type

float

getPuMass()

Get the mass of Pu in this object in grams.

getPuFrac()

Compute the Pu/HM mass fraction in this object.

Returns

puFrac – The pu mass fraction in heavy metal in this assembly

Return type

float

getZrFrac()

Return the total zr/(hm+zr) fraction in this assembly.

getMaxUraniumMassEnrich()

getFPMass()

Returns mass of fission products in this block in grams.

getFuelMass()

returns mass of fuel in grams.

constituentReport()

A print out of some pertinent constituent information.

getAtomicWeight()

Calculate the atomic weight of this object in g/mole of atoms.

Warning

This is not the molecular weight, which is grams per mole of molecules (grams/gram-molecule). That requires knowledge of the chemical formula. Don’t be surprised when you run this on UO2 and find it to be 90; there are a lot of Oxygen atoms in UO2.

\[A = \frac{\sum_i N_i A_i }{\sum_i N_i}\]

getMasses()

Return a dictionary of masses indexed by their nuclide names.

Notes

Implemented to get number densities and then convert to mass because getMass is too slow on a large tree.

getIntegratedMgFlux(adjoint=False, gamma=False)

getMgFlux(adjoint=False, average=False, volume=None, gamma=False)

Return the multigroup neutron flux in [n/cm^2/s].

The first entry is the first energy group (fastest neutrons). Each additional group is the next energy group, as set in the ISOTXS library.

On blocks, it is stored integrated over volume on <block>.p.mgFlux

Parameters

• adjoint (bool, optional) – Return adjoint flux instead of real

• average (bool, optional) – If true, will return average flux between latest and previous. Doesn’t work for pin detailed yet

• volume (float, optional) – If average=True, the volume-integrated flux is divided by volume before being returned. The user may specify a volume here, or the function will obtain the block volume directly.

• gamma (bool, optional) – Whether to return the neutron flux or the gamma flux.

Returns

flux – multigroup neutron flux in [n/cm^2/s]

Return type

numpy.array

removeMass(nucName, mass)

addMass(nucName, mass)

Add mass to a particular nuclide.

Parameters

• nucName (str) – nuclide name e.g. ‘U235’

• mass (float) – mass in grams of nuclide to be added to this armi Object

addMasses(masses)

Adds a vector of masses.

Parameters

masses (dict) – a dictionary of masses (g) indexed by nucNames (string)

setMass(nucName, mass)

Set the mass in an object by adjusting the ndens of the nuclides.

Parameters

• nucName (str) – Nuclide name to set mass of

• mass (float) – Mass in grams to set.

setMasses(masses)

Set a vector of masses.

Parameters

masses (dict) – a dictionary of masses (g) indexed by nucNames (string)

getSymmetryFactor()

Return a scaling factor due to symmetry on the area of the object or its children.

See also

armi.reactor.blocks.HexBlock.getSymmetryFactor

concrete implementation

getBoundingIndices()

Find the 3-D index bounds (min, max) of all children in the spatial grid of this object.

Returns

bounds – ((minI, maxI), (minJ, maxJ), (minK, maxK))

Return type

tuple

getComponentNames()

Get all unique component names of this Composite.

Returns

A set of all unique component names found in this Composite.

Return type

set or str

getComponentsOfShape(shapeClass)

Return list of components in this block of a particular shape.

Parameters

shapeClass (Component) – The class of component, e.g. Circle, Helix, Hexagon, etc.

Returns

param – List of components in this block that are of the given shape.

Return type

list

getComponentsOfMaterial(material=None, materialName=None)

Return list of components in this block that are made of a particular material.

Only one of the selectors may be used

Parameters

• material (Material object, optional) – The material to match

• materialName (str, optional) – The material name to match.

Returns

componentsWithThisMat

Return type

list

hasComponents(typeSpec: Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]], None, List[Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]]]], exact=False)

Return true if components matching all TypeSpec exist in this object.

Parameters

typeSpec (Flags or iterable of Flags) – Component flags to check for

getComponentByName(name)

Gets a particular component from this object, based on its name.

Parameters

name (str) – The blueprint name of the component to return

getComponent(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exact=False, quiet=False)

Get a particular component from this object.

Parameters

• typeSpec (flags.Flags or list of Flags) – The type specification of the component to return

• exact (boolean, optional) – Demand that the component flags be exactly equal to the typespec. Default: False

• quiet (boolean, optional) – Warn if the component is not found. Default: False

• object (Careful with multiple similar names in one) –

Returns

Component

Return type

The component that matches the critera or None

getNumComponents(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exact=False)

Get the number of components that have these flags, taking into account multiplicity. Useful for getting nPins even when there are pin detailed cases.

Parameters

typeSpec (Flags) – Expected flags of the component to get. e.g. Flags.FUEL

Returns

total – the number of components of this type in this object, including multiplicity.

Return type

int

setComponentDimensionsReport()

Makes a summary of the dimensions of the components in this object.

printDensities(expandFissionProducts=False)

Get lines that have the number densities of a object.

expandAllElementalsToIsotopics()

expandElementalToIsotopics(elementalNuclide)

Expands the density of a specific elemental nuclides to its natural isotopics.

Parameters

elementalNuclide (armi.nucDirectory.nuclideBases.NaturalNuclide) – natural nuclide to replace.

getAverageTempInC(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Return the average temperature of the ArmiObject in C by averaging all components.

resolveLinkedDims(components)

Resolve link strings to links on all child components.

getDominantMaterial(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Return the first sample of the most dominant material (by volume) in this object.

Parameters

• typeSpec (Flags or iterable of Flags, optional) – The types of components to consider (e.g. [Flags.FUEL, Flags.CONTROL])

• exact (bool, optional) – Whether or not the TypeSpec is exact

Returns

mat – the first instance of the most dominant material (by volume) in this object.

Return type

Material

See also

getComponentsOfMaterial

Gets components that are made of a particular material

gatherMaterialsByVolume

Classifies all materials by volume

class armi.reactor.composites.Composite(name)

Bases: armi.reactor.composites.ArmiObject

An ArmiObject that has children.

This is a fundamental ARMI state object that generally represents some piece of the nuclear reactor that is made up of other smaller pieces. This object can cache information about its children to help performance.

Details about spatial representation

Spatial representation of a Composite is handled through a combination of the spatialLocator and spatialGrid parameters. The spatialLocator is a numpy triple representing either:

1. Indices in the parent’s spatialGrid (for lattices, etc.), used when the dtype is int.

2. Coordinates in the parent’s universe in cm, used when the dtype is float.

The top parent of any composite must have a coordinate-based spatialLocator. For example, a Reactor an a Pump should both have coordinates based on how far apart they are.

The traversal of indices and grids is recursive. The Reactor/Core/Assembly/Block model is handled by putting a 2-D grid (either Theta-R, Hex, or Cartesian) on the Core and individual 1-D Z-meshes on the assemblies. Then, Assemblies have 2-D spatialLocators (i,j,0) and Blocks have 1-D spatiaLocators (0,0,k). These get added to form the global indices. This way, if an assembly is moved, all the blocks immediately and naturally move with it. Individual children may have coordinate-based spatialLocators mixed with siblings in a grid. This allows mixing grid-representation with explicit representation, often useful in advanced assemblies and thermal reactors.

The traversal of indices and grids is recursive. The Reactor/Core/Assembly/Block model is handled by putting a 2-D grid (either Theta-R, Hex, or Cartesian) on the Core and individual 1-D Z-meshes on the assemblies. Then, Assemblies have 2-D spatialLocators (i,j,0) and Blocks have 1-D spatiaLocators (0,0,k). These get added to form the global indices. This way, if an assembly is moved, all the blocks immediately and naturally move with it. Individual children may have coordinate-based spatialLocators mixed with siblings in a grid. This allows mixing grid-representation with explicit representation, often useful in advanced assemblies and thermal reactors.

sort()

Sort the children of this object.

index(obj)

Obtain the list index of a particular child.

append(obj)

Append a child to this object.

extend(seq)

Add a list of children to this object.

add(obj)

Add one new child.

remove(obj)

Remove a particular child.

moveTo(locator)

Move to specific location in parent. Often in a grid.

insert(index, obj)

Insert an object into the list of children at a particular index.

removeAll()

Remove all children.

setChildren(items)

Clear this container and fills it with new children.

getChildren(deep=False, generationNum=1, includeMaterials=False, predicate=None)

Return the children objects of this composite.

Parameters

• deep (boolean, optional) – Return all children of all levels.

• generationNum (int, optional) – Which generation to return. 1 means direct children, 2 means children of children. Setting this parameter will only return children of this generation, not their parents. Default: Just return direct children.

• includeMaterials (bool, optional) – Include the material properties

• predicate (callable, optional) – An optional unary predicate to use for filtering results. This can be used to request children of specific types, or with desired attributes. Not all ArmiObjects have the same methods and members, so care should be taken to make sure that the predicate executes gracefully in all cases (e.g., use getattr(obj, "attribute", None) to access instance attributes). Failure to meet the predicate only affects the object in question; children will still be considered.

Examples

>>> obj.getChildren()
[child1, child2, child3]

>>> obj.getChildren(generationNum=2)
[grandchild1, grandchild2, grandchild3]

>>> obj.getChildren(deep=True)
[child1, child2, child3, grandchild1, grandchild2, grandchild3]

# Assuming that grandchild1 and grandchild3 are Component objects >>> obj.getChildren(deep=True, predicate=lambda o: isinstance(o, Component)) [grandchild1, grandchild3]

getChildrenWithFlags(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]], exactMatch=False)

Get all children of a specific type.

getChildrenOfType(typeName)

Get children that have a specific input type name.

getComponents(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

iterComponents(typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Return an iterator of armi.reactor.component.Component objects within this Composite.

Parameters

• typeSpec (TypeSpec) – Component flags. Will restrict Components to specific ones matching the flags specified.

• exact (bool, optional) – Only match exact component labels (names). If True, ‘coolant’ will not match ‘interCoolant’. This has no impact if typeSpec is None.

Returns

items matching typeSpec and exact criteria

Return type

iterator of Component

syncMpiState()

Synchronize all parameters of this object and all children to all worker nodes over the network using MPI.

In parallelized runs, if each process has its own copy of the entire reactor hierarchy, this method synchronizes the state of all parameters on all objects.

Returns

number of parameters synchronized over all components

Return type

int

retainState(paramsToApply=None)

Restores a state before and after some operation.

Parameters

paramsToApply (iterable) – Parameters that should be applied to the state after existing the state retainer. All others will be reverted to their values upon entering.

Notes

This should be used in a with statement.

backUp()

Create and store a backup of the state.

This needed to be overridden due to linked components which actually have a parameter value of another ARMI component.

restoreBackup(paramsToApply)

Restore the parameters from previously created backup.

Parameters

paramsToApply (list of ParmeterDefinitions) – restores the state of all parameters not in paramsToApply

getLumpedFissionProductsIfNecessary(nuclides=None)

Return Lumped Fission Product objects that belong to this object or any of its children.

getLumpedFissionProductCollection()

Get collection of LFP objects. Will work for global or block-level LFP models.

Returns

lfps – lfpName keys, lfp object values

Return type

object

See also

armi.physics.neutronics.fissionProductModel.lumpedFissionProduct.LumpedFissionProduct

LFP object

requiresLumpedFissionProducts(nuclides=None)

True if any of the nuclides in this object are Lumped nuclides.

getIntegratedMgFlux(adjoint=False, gamma=False)

Returns the multigroup neutron tracklength in [n-cm/s].

The first entry is the first energy group (fastest neutrons). Each additional group is the next energy group, as set in the ISOTXS library.

Parameters

• adjoint (bool, optional) – Return adjoint flux instead of real

• gamma (bool, optional) – Whether to return the neutron flux or the gamma flux.

Returns

integratedFlux – multigroup neutron tracklength in [n-cm/s]

Return type

numpy.array

getReactionRates(nucName, nDensity=None)

Get the reaction rates of a certain nuclide on this object.

Parameters

• nucName (str) – nuclide name – e.g. ‘U235’

• nDensity (float) – number Density

Returns

rxnRates – reaction rates (1/s) for nG, nF, n2n, nA and nP

Return type

dict

Notes

This is volume integrated NOT (1/cm3-s)

If you set nDensity to 1 this makes 1-group cross section generation easier

printContents(includeNuclides=True)

Display information about all the comprising children in this object.

isOnWhichSymmetryLine()

getBoundingCircleOuterDiameter(Tc=None, cold=False)

Get sum circle bound.

Used to roughly approximate relative size vs. other objects

paramCollectionType

alias of armi.reactor.parameters.parameterCollections.ArmiObjectParameterCollection

class armi.reactor.composites.StateRetainer(composite, paramsToApply=None)

Bases: object

Retains state during some operations.

This can be used to temporarily cache state, perform an operation, extract some info, and then revert back to the original state.

• A state retainer is faster than restoring state from a database as it reduces the number of IO reads; however, it does use more memory.

• This can be used on any object within the composite pattern via with [rabc].retainState([list], [of], [parameters], [to], [retain]):. Use on an object up in the hierarchy applies to all objects below as well.

• This is intended to work across MPI, so that if you were to broadcast the reactor the state would be correct; however the exact implication on parameters may be unclear.

Create an instance of a StateRetainer.

Parameters

• composite (Composite) – composite object to retain state (recursively)

• paramsToApply (iterable of parameters.Parameter) – Iterable of parameters.Parameter to retain updated values after __exit__. All other parameters are reverted to the original state, i.e. retained at the original value.

armi.reactor.composites.gatherMaterialsByVolume(objects: List[armi.reactor.composites.ArmiObject], typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Compute the total volume of each material in a set of objects and give samples.

Parameters

• objects (list of ArmiObject) – Objects to look within. This argument allows clients to search though some subset of the three (e.g. when you’re looking for all CLADDING components within FUEL blocks)

• typeSpec (TypeSpec) – Flags for the components to look at

• exact (bool) – Whether or not the TypeSpec is exact

Notes

This helper method is outside the main ArmiObject tree for the special clients that need to filter both by container type (e.g. Block type) with one set of flags, and Components with another set of flags.

Warning

This is a composition related helper method that will likely be filed into classes/modules that deal specifically with the composition of things in the data model. Thus clients that use it from here should expect to need updates soon.

armi.reactor.composites.getDominantMaterial(objects: List[armi.reactor.composites.ArmiObject], typeSpec: Optional[Union[armi.utils.flags.Flag, armi.utils.flags.auto, Sequence[Union[armi.utils.flags.Flag, armi.utils.flags.auto]]]] = None, exact=False)

Return the first sample of the most dominant material (by volume) in a set of objects.

Warning

This is a composition related helper method that will likely be filed into classes/modules that deal specifically with the composition of things in the data model. Thus clients that use it from here should expect to need updates soon.

armi.reactor.composites.getReactionRateDict(nucName, lib, xsSuffix, mgFlux, nDens)

Parameters

• nucName (str) – nuclide name – e.g. ‘U235’, ‘PU239’, etc. Not to be confused with the nuclide _label_, see the nucDirectory module for a description of the difference.

• lib (isotxs) – cross section library

• xsSuffix (str) – cross section suffix, consisting of the type followed by the burnup group, e.g. ‘AB’ for the second burnup group of type A

• mgFlux (numpy.nArray) – integrated mgFlux (n-cm/s)

• nDens (float) – number density (atom/bn-cm)

Returns

dictionary of reaction rates (rxn/s) for nG, nF, n2n, nA and nP

Return type

rxnRates - dict

Notes

Assume there is no n3n cross section in ISOTXS