HydrideVolumeFraction

Computes the hydride volume fraction as a degree 3 polynomial of temperature.

note:Often Created by an Action

This object can be set up automatically by using the CladdingHydrides action.

Description

This material computes the hydride formation energy used by HydrideNucleationKinetics, HydrideGrowthKinetics and HydrogenSolubility.

Volume fraction computation

The volume fraction of hydrides is given by

$var element = document.getElementById("moose-equation-b5597144-ddb6-410a-a698-5a152a66213f");katex.render(" v_\\delta = \\dfrac{x_{Prec}}{x_{\\delta} - TSS_D^{at}}.", element, {displayMode:true,throwOnError:false});$

$var element = document.getElementById("moose-equation-ebd51627-815b-4e01-afae-3b7ef951056b");katex.render("x_{Prec}", element, {displayMode:false,throwOnError:false});$ is the atomic fraction of hydrides (Passelaigue et al., 2021)

$var element = document.getElementById("moose-equation-1b3396eb-c82c-44bb-b250-600dba5af483");katex.render("x_{Prec} = \\dfrac{\\dfrac{C_{Prec}}{M_H}}{\\dfrac{C_{tot}}{M_H}+\\dfrac{10^6-C_{tot}}{M_{Zr}}},", element, {displayMode:true,throwOnError:false});$

$var element = document.getElementById("moose-equation-d4a71c76-51f4-45dc-9e29-fa1ff3a384f1");katex.render("TSS_D^{at}", element, {displayMode:false,throwOnError:false});$ is the hydrogen solubility expressed as an atomic fraction

$var element = document.getElementById("moose-equation-18709cbf-bc20-4682-b08a-be7ea28737ef");katex.render(" TSS_D^{at} = \\dfrac{\\dfrac{TSS_D^{eff}}{M_H}}{\\dfrac{TSS_D^{eff}}{M_H} + \\dfrac{10^6 - TSS_D^{eff}}{M_{Zr}}}", element, {displayMode:true,throwOnError:false});$

and $var element = document.getElementById("moose-equation-05ee5896-0c72-4c5f-912a-f4559ec730ae");katex.render("x_{\\delta}", element, {displayMode:false,throwOnError:false});$ is the atomic fraction of hydrogen at the $var element = document.getElementById("moose-equation-a58938b4-3ef6-4a39-822d-536e3860fdcb");katex.render("(\\delta / \\alpha+\\delta)", element, {displayMode:false,throwOnError:false});$ boundary, fitted with a degree 3 polynomial of the temperature as

$var element = document.getElementById("moose-equation-659fdf6a-518d-4684-81a7-dac275ad5f32");katex.render("x_\\delta(T) = x_0 + x_1 T + x_2 T^2 + x_3 T^3,", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-f21ebecd-8a89-4cd3-af67-b3432c02c974");katex.render("x_0 = 0.623", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-aa9f2b39-358f-438c-94b7-a8d29f1f96f4");katex.render("x_1 = -5.73\\times 10^{-5} K^{-1}", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-4646df56-5ce8-432d-b3a5-997d1201231b");katex.render("x_2 = 8.48\\times 10^{-8} K^{-2}", element, {displayMode:false,throwOnError:false});$ , and $var element = document.getElementById("moose-equation-d5b6ae95-c242-4ff8-9065-67b643a10516");katex.render("x_3 = -9.93\\times 10^{-11} K^{-3}", element, {displayMode:false,throwOnError:false});$ (Lacroix, 2019; Lacroix et al., 2021).

This class also exports the value of $var element = document.getElementById("moose-equation-ea4e2a9e-2804-403c-ae36-eac824a6717c");katex.render("x_{\\delta}", element, {displayMode:false,throwOnError:false});$ used by HydrideGrowthKinetics.

Example Input Syntax

[CladdingHydrides<<<{"href": "../../syntax/CladdingHydrides/index.html"}>>>]
  [hydrides]
    block<<<{"description": "The list of block ids for the cladding."}>>> = 0
    temperature<<<{"description": "Temperature (K)."}>>> = temp
    hydrogen_in_solution_ppm<<<{"description": "Concentration of dissolved hydrogen (wt.ppm)."}>>> = Css
    hydrogen_as_hydride_ppm<<<{"description": "Concentration of hydrogen as hydride (wt.ppm)."}>>> = Cprec
    solubility_frequency_factor<<<{"description": "Frequency factor for Arrhenius TSSd (s-1)."}>>> = 67116
    solubility_activation_energy<<<{"description": "Activation energy for Arrhenius TSSd (J/mol)."}>>> = 32294
    # solubility_g = 0
    # solubility_delta = 1
  []
[]

References

E. Lacroix. Modeling Zirconium Hydride Precipitation and Dissolution in Zirconium Alloys. PhD thesis, The Pennsylvania State University, 2019.

@phdthesis{lacroixPhd,
    Author = "Lacroix, E.",
    Title = "Modeling Zirconium Hydride Precipitation and Dissolution in Zirconium Alloys",
    Year = "2019",
    School = "The Pennsylvania State University"
}

E. Lacroix, P.-C. A. Simon, A. T. Motta, and J.D. Almer. Zirconium hydride precipitation and dissolution kinetics in the hysteresis region in zirconium alloys. Zirconium in the Nuclear Industry: 19th International Symposium, ASTM STP 1597, pages 67–91, 2021.

@article{Lacroix2019,
    author = "Lacroix, E. and Simon, P.-C. A. and Motta, A. T. and Almer, J.D.",
    title = "{Zirconium hydride precipitation and dissolution kinetics in the hysteresis region in zirconium alloys}",
    journal = "Zirconium in the Nuclear Industry: 19th International Symposium, ASTM STP 1597",
    year = "2021",
    pages = "67-91"
}

F. Passelaigue, E. Lacroix, G. Pastore, and A.T. Motta. Implementation and Validation of the Hydride Nucleation-Growth-Dissolution (HNGD) model in BISON. Journal of Nuclear Materials, 544:152683, 2021.

@article{Passelaigue2021,
    author = "Passelaigue, F. and Lacroix, E. and Pastore, G. and Motta, A.T.",
    journal = "Journal of Nuclear Materials",
    title = "{Implementation and Validation of the Hydride Nucleation-Growth-Dissolution (HNGD) model in BISON}",
    volume = "544",
    year = "2021",
    pages = "152683"
}

(test/tests/hydrogen/benchmark_transients.i)

# Benchmark : Precipitation and dissolution kinetics (E. Lacroix, PhD Dissertation (2019), 5.2.3)
# In this simulation, the sample is subjected to a uniform thermal treatment that shows the kinetics
# of Nucleation, Growth and Dissolution of hydrides in a zircaloy cladding.

[GlobalParams]
  order = FIRST
  family = LAGRANGE
[]

[Mesh]
  [generated_mesh]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = 25.4e-5
    ymin = 0
    ymax = 25e-5
    nx = 1
    ny = 1
    bias_x = 1
  []
[]

[Variables]
  [temp]
  []

  [Css]
    scaling=1e6
  []

  [Cprec]
    order = FIRST
    family = LAGRANGE
    scaling=1e6
  []
[]

[AuxVariables]
  # Used to prevent negative concentrations
  [bounds_dummy]
    order = FIRST
    family = LAGRANGE
  []
[]

[Bounds]
  # To prevent negative concentrations
  [Cprec_lower_bound]
    type = ConstantBounds
    variable = bounds_dummy
    bounded_variable = Cprec
    bound_type = lower
    bound_value = 0.
  []
[]

[ICs]
  [temp_init]
    variable=temp
    value=300
    type=ConstantIC
 []
 [Css_init]
   variable=Css
   value=0.06
   type=ConstantIC
 []
 [Cprec_init]
   variable=Cprec
   value=253.94
   type=ConstantIC
 []
[]

[Kernels]
  [heatflow]
    type = CoefDiffusion
    variable = temp
    coef = 1e-4
  []
  [dTdt]
    type = TimeDerivative
    variable = temp
  []
[]

[CladdingHydrides]
  [hydrides]
    block=0
    temperature = temp
    hydrogen_in_solution_ppm = Css
    hydrogen_as_hydride_ppm = Cprec
    solubility_frequency_factor = 67116
    solubility_activation_energy = 32294
    # solubility_g = 0
    # solubility_delta = 1
  []
[]

[Functions]
  [fct_temp]
    type = PiecewiseLinear
    x = '0    1969  2594    3877  4411  5499  6096  6609  7705  8419  8988  1.02E+04 1.03E+04 1.04E+04 1.05E+04 1.06E+04 1.07E+04'
    y = '300  697.4 697.4   442.7 442.7 656.6 656.6 557.7 558.1 697.4 697.8 451.2    406.4    374.2    350.89   330.3    323.13'
  []
[]

[BCs]
  [T]
    boundary= '0'
    type=FunctionDirichletBC
    variable=temp
    function = fct_temp
  []
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]

  type = Transient

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -snes_type'
  petsc_options_value = 'lu vinewtonrsls' # This petsc option help prevent negative concentrations

  line_search = 'none'

  l_max_its = 100
  l_tol = 8e-3

  nl_max_its = 15
  nl_rel_tol = 1e-6
  nl_abs_tol = 1e-10

  start_time = 0
  end_time = 1.07E+04
  n_startup_steps = 1

  dtmax = 20
  dtmin = 1

  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 20
    timestep_limiting_function = fct_temp
    force_step_every_function_point = true
    optimal_iterations = 60
    iteration_window = 0
    linear_iteration_ratio = 100
  []
[]

[Postprocessors]
  # [temp]
  #   type = ElementAverageValue
  #   variable = 'temp'
  #   block = 0
  # []
  [Css]
    type = ElementAverageValue
    variable = 'Css'
    block = 0
  []
[]

[Outputs]
  # exodus = true
  [console]
    type = Console
    max_rows = 25
  []
  execute_on = 'timestep_end'
  csv = true
[]

Description
Volume fraction computation
Example Input Syntax
References