From 27060c7bcc0744a99a638fafc77ef183a0d1c8b2 Mon Sep 17 00:00:00 2001
From: RemDelaporteMathurin <rdelaportemathurin@gmail.com>
Date: Tue, 17 Sep 2024 11:32:28 -0400
Subject: [PATCH] added Okana properties

---
 h_transport_materials/material.py             |  5 +-
 .../property_database/palladium_copper.py     | 53 +++++++++++++------
 h_transport_materials/references.bib          | 16 ++++++
 3 files changed, 58 insertions(+), 16 deletions(-)

diff --git a/h_transport_materials/material.py b/h_transport_materials/material.py
index 7aff11f..e1cbd9f 100644
--- a/h_transport_materials/material.py
+++ b/h_transport_materials/material.py
@@ -113,9 +113,11 @@ class MoltenSalt(Material):
 class FeCrAl(Steel):
     family = "fecral"
 
+
 class PalladiumAlloy(Alloy):
     family = "palladium alloy"
 
+
 TUNGSTEN = Tungsten()
 BERYLLIUM = Beryllium()
 CARBON = Carbon()
@@ -165,4 +167,5 @@ class PalladiumAlloy(Alloy):
 
 TZM = Alloy("tzm")
 
-PD52CU = PalladiumAlloy("pd52cu")
\ No newline at end of file
+PD60CU40 = PalladiumAlloy("60pd40cu")
+PD52CU = PalladiumAlloy("pd52cu")
diff --git a/h_transport_materials/property_database/palladium_copper.py b/h_transport_materials/property_database/palladium_copper.py
index ff5ee9e..ac7da2c 100644
--- a/h_transport_materials/property_database/palladium_copper.py
+++ b/h_transport_materials/property_database/palladium_copper.py
@@ -2,19 +2,22 @@
 from h_transport_materials import (
     Permeability,
     Diffusivity,
-) 
+)
 import numpy as np
 
 u = htm.ureg
 
-li_data_T = np.array(
-    [
-        350,
-        375,
-        400,
-        425,
-    ]
-) * u.degC
+li_data_T = (
+    np.array(
+        [
+            350,
+            375,
+            400,
+            425,
+        ]
+    )
+    * u.degC
+)
 
 li_data_y = (
     np.array(
@@ -25,9 +28,9 @@
             1.60e-8,
         ]
     )
-    * u.mol 
-    * u.m**-1 
-    * u.Pa**-0.5 
+    * u.mol
+    * u.m**-1
+    * u.Pa**-0.5
     * u.s**-1
 )
 
@@ -36,7 +39,8 @@
     data_y=li_data_y,
     source="li_low_2023",
     isotope="H",
-    note="SI Table 1 (supporting information)"
+    note="SI Table 1 (supporting information)",
+    material=htm.PD52CU,
 )
 
 piper_diffusivity_h = Diffusivity(
@@ -46,14 +50,33 @@
     isotope="H",
     source="piper_diffusion_2004",
     note="Equation 6 - this Arrhenius fit holds for when the alloy is in the beta phase for 52.5 percent copper. These range values were found outside of the paper.",
+    material=htm.PD52CU,
+)
+
+onaka_diffusivity = Diffusivity(
+    D_0=2.9e-8 * u.m**2 * u.s**-1,
+    E_D=0.048 * u.eV * u.particle**-1,
+    isotope="D",
+    range=(u.Quantity(350, u.K), u.Quantity(600, u.K)),
+    source="onaka_characteristic_2016",
+    material=htm.PD60CU40,
+)
+
+onaka_recombination = htm.RecombinationCoeff(
+    pre_exp=1.43e-26 * u.m**4 * u.s**-1 * u.particle**-1,
+    act_energy=0.40 * u.eV * u.particle**-1,
+    isotope="D",
+    range=(u.Quantity(350, u.K), u.Quantity(600, u.K)),
+    source="onaka_characteristic_2016",
+    material=htm.PD60CU40,
 )
 
 properties = [
     li_permeability_h,
     piper_diffusivity_h,
+    onaka_recombination,
+    onaka_diffusivity,
 ]
 
-for prop in properties:
-    prop.material = htm.PD52CU
 
 htm.database += properties
diff --git a/h_transport_materials/references.bib b/h_transport_materials/references.bib
index 9203cc3..69fa529 100644
--- a/h_transport_materials/references.bib
+++ b/h_transport_materials/references.bib
@@ -2708,6 +2708,22 @@ @article{shirasu_thermodynamic_1993
 	pages = {218--222},
 }
 
+@article{onaka_characteristic_2016,
+	title = {Characteristic of a {pdCu} membrane as atomic hydrogen probe for {QUEST}},
+	volume = {9},
+	issn = {2352-1791},
+	url = {https://www.sciencedirect.com/science/article/pii/S235217911530048X},
+	doi = {10.1016/j.nme.2016.09.001},
+	abstract = {A permeation probe is a useful device for detecting the atomic hydrogen flux to plasma-facing walls. Recently, we developed a new type of probe using 60Pd-40Cu alloy (PdCu) as the permeation material. The deuterium behaviors in PdCu samples were investigated using nuclear reaction analysis (NRA) and permeation observations, and the diffusion coefficient and recombination coefficients were determined from these observations. It was found that the sensitivity of a 0.02-mm-thick probe was as high as 0.5 below 473K and was independent of the incident flux. The response time at 473K was 0.41s and 1.3s under an incident flux of 1020m−2s−1 and 1019m−2s−1, respectively. Thus, we concluded that the new PdCu probe can effectively detect the incident atomic hydrogen with high sensitivity and a suitable response time.},
+	urldate = {2024-09-17},
+	journal = {Nuclear Materials and Energy},
+	author = {Onaka, Masayuki and Takagi, Ikuji and Kobayashi, Taishi and Sasaki, Takayuki and Kuzmin, Arseniy and Zushi, Hideki},
+	month = dec,
+	year = {2016},
+	keywords = {Diffusion, Hydrogen, Nuclear reaction analysis, Palladium alloy, Permeation, Plasma-wall interaction},
+	pages = {104--108},
+}
+
 @article{fuerst_hastelloyn_2024,
 	title = {Hydrogen and deuterium permeation in Hastelloy N},
 	journal = {Journal of Nuclear Materials},