HydrogenSuperSolubility

Computes the supersolubility of hydrogen in Zr cladding.

note:Often Created by an Action

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

Description

This material computes the hydrogen solubility used by HydridePrecipitationRate.

Hydrogen Supersolubility

The supersolubility of hydrogen is based on the $var element = document.getElementById("moose-equation-7ba8b4b9-69b2-4742-bc09-3adcc6c3b825");katex.render("TSS_P", element, {displayMode:false,throwOnError:false});$ that follows an Arrhenius law:

$var element = document.getElementById("moose-equation-e8d03e54-6dd0-4d44-89df-a40e95202f2a");katex.render(" TSS_P = TSS_{P0} \\exp\\left(\\frac{-Q_P}{R T}\\right)", element, {displayMode:true,throwOnError:false});$

with $var element = document.getElementById("moose-equation-89cc3321-97c9-48f5-b24f-31f593e2371d");katex.render("Q_P = 35459", element, {displayMode:false,throwOnError:false});$ J/mol and $var element = document.getElementById("moose-equation-0ed43a95-c961-402d-b3d1-b8304a05709d");katex.render("TSS_{P0}=101999", element, {displayMode:false,throwOnError:false});$ wt.ppm as default values (Lacroix, 2019; Lacroix et al., 2021; Passelaigue et al., 2021).

It was further developed in (Passelaigue et al., 2022) to include the time dependency of the supersolubility. During temperature transients, the supersolubility (or effective $var element = document.getElementById("moose-equation-3920b72f-4896-4066-bfde-1c25add2e4e9");katex.render("TSS_P", element, {displayMode:false,throwOnError:false});$ ) follows the usual $var element = document.getElementById("moose-equation-e80e2d56-68f5-4941-9543-4e163dd8f17a");katex.render("TSS_P", element, {displayMode:false,throwOnError:false});$ :

$var element = document.getElementById("moose-equation-aad6a256-e22c-4e19-9768-1d91bb36e7bd");katex.render(" TSS_P^{eff} = TSS_P", element, {displayMode:true,throwOnError:false});$

During annealings, the supersolubility tends exponentially from $var element = document.getElementById("moose-equation-bad76f26-0bdc-46c8-8a7a-0cd4f49a2d37");katex.render("TSS_P", element, {displayMode:false,throwOnError:false});$ to the solubility (or effective $var element = document.getElementById("moose-equation-e4431479-8097-4c0f-ac11-df7b36abcc91");katex.render("TSS_D", element, {displayMode:false,throwOnError:false});$ ):

$var element = document.getElementById("moose-equation-a9b9aeeb-447f-4b01-8744-3e5d1cb9d228");katex.render(" TSS_P^{eff} = TSS_D^{eff} + (TSS_P - TSS_D^{eff}) \\exp\\left(-\\frac{t-t_0}{\\tau}\\right)", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-a03638a2-7d38-4f23-bf6b-f4d14069a9f0");katex.render("\\tau = 10^4", element, {displayMode:false,throwOnError:false});$ s is the characteristic time for the supersolubility decrease, $var element = document.getElementById("moose-equation-73e5eeda-eb03-417a-ac43-fbdda3c35494");katex.render("t_0", element, {displayMode:false,throwOnError:false});$ is the time when the annealing started (determined based on the time derivative of the temperature), and $var element = document.getElementById("moose-equation-2d8d4815-2ccf-4e41-ab93-80ca3518d06f");katex.render("TSS_D^{eff}", element, {displayMode:false,throwOnError:false});$ is computed in HydrogenSolubility.

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"
}

F. Passelaigue, P.-C. A. Simon, and A.T. Motta. Predicting the hydride rim by improving the solubility limits in the Hydride Nucleation-Growth-Dissolution (HNGD) model. Journal of Nuclear Materials, 558:153363, 2022.

@article{Passelaigue2022,
    author = "Passelaigue, F. and Simon, P.-C. A. and Motta, A.T.",
    journal = "Journal of Nuclear Materials",
    title = "{Predicting the hydride rim by improving the solubility limits in the Hydride Nucleation-Growth-Dissolution (HNGD) model}",
    volume = "558",
    year = "2022",
    pages = "153363"
}

(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
Hydrogen Supersolubility
Example Input Syntax
References