Source code for armi.materials.thU

# Copyright 2019 TerraPower, LLC
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#     http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

"""
Thorium Uranium metal

Data is from [#IAEA-TECDOCT-1450]_.

.. [#IAEA-TECDOCT-1450] Thorium fuel cycle -- Potential benefits and challenges, IAEA-TECDOC-1450 (2005).
    https://www-pub.iaea.org/mtcd/publications/pdf/te_1450_web.pdf
"""

from armi.utils.units import getTk
from armi.materials.material import Material
from armi import runLog


[docs]class ThU(Material):
    name = "ThU"
    enrichedNuclide = "U233"

[docs]    def getEnrichment(self):
        return self.getMassFrac("U233") / (
            self.getMassFrac("U233") + self.getMassFrac("TH232")
        )

[docs]    def applyInputParams(self, U233_wt_frac=None):
        runLog.warning("Material {} has not yet been tested for accuracy".format("ThU"))

        if U233_wt_frac:
            self.adjustMassEnrichment(U233_wt_frac)

[docs]    def setDefaultMassFracs(self):
        self.setMassFrac("TH232", 1.0)
        self.setMassFrac("U233", 0.0)

[docs]    def density(self, Tk=None, Tc=None):
        Tk = getTk(Tc, Tk)
        """g/cc from IAEA TE 1450"""
        return 11.68

[docs]    def linearExpansion(self, Tk=None, Tc=None):
        r"""m/m/K from IAEA TE 1450"""
        Tk = getTk(Tc, Tk)
        self.checkTempRange(30, 600, Tk, "linear expansionn")
        return 11.9e-6

[docs]    def thermalConductivity(self, Tk=None, Tc=None):
        r"""W/m-K from IAEA TE 1450"""
        Tk = getTk(Tc, Tk)
        return 43.1

[docs]    def meltingPoint(self):
        r"""melting point in K from IAEA TE 1450"""
        return 2025.0