Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations
Abstract
Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$$$N_\text {V}$$$$ ), He atoms ( $$$$N_\text {He}$$$$ ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$$$N_\text {V}$$$$ vacancies and $$$$N_\text {He}$$$$ helium atoms, the $$$$N_\text {He}/N_\text {V}$$$$ ratio that causes the event, the resulting increase in $$$$N_\text {V}$$$$ , and the associated shift of the bubble depth are formulated as a function of $$$$N_\text {V}$$$$ and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$$$N_\text {He}/N_\text {V}$$$$ ratio in order to burst. The burst depth and $$$$N_\text {He}/N_\text {V}$$$$ are also modeled as a function of $$$$N_\text {V}$$$$ and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.
- Authors:
- Publication Date:
- Research Org.:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Fusion Energy Sciences (FES)
- OSTI Identifier:
- 1984907
- Alternate Identifier(s):
- OSTI ID: 1986526
- Report Number(s):
- PNNL-SA-183858
Journal ID: ISSN 2045-2322; 9601; PII: 35803
- Grant/Contract Number:
- AC05-76RL0-1830; 4000159718; 4000159698; SC0018421; AC05-76RL01830; SC0006661
- Resource Type:
- Published Article
- Journal Name:
- Scientific Reports
- Additional Journal Information:
- Journal Name: Scientific Reports Journal Volume: 13 Journal Issue: 1; Journal ID: ISSN 2045-2322
- Publisher:
- Nature Publishing Group
- Country of Publication:
- United Kingdom
- Language:
- English
- Subject:
- 36 MATERIALS SCIENCE; fusion energy materials; materials science; physics
Citation Formats
Setyawan, Wahyu, Dasgupta, Dwaipayan, Blondel, Sophie, Nandipati, Giridhar, Hammond, Karl D., Maroudas, Dimitrios, and Wirth, Brian D. Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations. United Kingdom: N. p., 2023.
Web. doi:10.1038/s41598-023-35803-3.
Setyawan, Wahyu, Dasgupta, Dwaipayan, Blondel, Sophie, Nandipati, Giridhar, Hammond, Karl D., Maroudas, Dimitrios, & Wirth, Brian D. Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations. United Kingdom. https://doi.org/10.1038/s41598-023-35803-3
Setyawan, Wahyu, Dasgupta, Dwaipayan, Blondel, Sophie, Nandipati, Giridhar, Hammond, Karl D., Maroudas, Dimitrios, and Wirth, Brian D. Tue .
"Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations". United Kingdom. https://doi.org/10.1038/s41598-023-35803-3.
@article{osti_1984907,
title = {Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations},
author = {Setyawan, Wahyu and Dasgupta, Dwaipayan and Blondel, Sophie and Nandipati, Giridhar and Hammond, Karl D. and Maroudas, Dimitrios and Wirth, Brian D.},
abstractNote = {Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.},
doi = {10.1038/s41598-023-35803-3},
journal = {Scientific Reports},
number = 1,
volume = 13,
place = {United Kingdom},
year = {Tue Jun 13 00:00:00 EDT 2023},
month = {Tue Jun 13 00:00:00 EDT 2023}
}
https://doi.org/10.1038/s41598-023-35803-3
Works referenced in this record:
Tungsten surface evolution by helium bubble nucleation, growth and rupture
journal, June 2013
- Sefta, Faiza; Hammond, Karl D.; Juslin, Niklas
- Nuclear Fusion, Vol. 53, Issue 7
Helium segregation and transport behavior near ⟨100⟩ and ⟨110⟩ symmetric tilt grain boundaries in tungsten
journal, June 2018
- Yang, Zhuocen; Hu, Lin; Maroudas, Dimitrios
- Journal of Applied Physics, Vol. 123, Issue 22
Equation of state and melting properties of from measurements to 20 kbar
journal, June 1980
- Mills, R. L.; Liebenberg, D. H.; Bronson, J. C.
- Physical Review B, Vol. 21, Issue 11
Operating temperature windows for fusion reactor structural materials
journal, November 2000
- Zinkle, S. J.; Ghoniem, N. M.
- Fusion Engineering and Design, Vol. 51-52, p. 55-71
Formation of helium induced nanostructure ‘fuzz’ on various tungsten grades
journal, September 2010
- Baldwin, M. J.; Doerner, R. P.
- Journal of Nuclear Materials, Vol. 404, Issue 3
Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations
journal, May 2019
- Hammond, Karl D.; Naeger, Ian V.; Widanagamaachchi, Wathsala
- Nuclear Fusion, Vol. 59, Issue 6
The European power plant conceptual study
journal, November 2005
- Maisonnier, D.; Cook, I.; Pierre, Sardain
- Fusion Engineering and Design, Vol. 75-79
Helium segregation on surfaces of plasma-exposed tungsten
journal, January 2016
- Maroudas, Dimitrios; Blondel, Sophie; Hu, Lin
- Journal of Physics: Condensed Matter, Vol. 28, Issue 6
Helium in-plane migration behavior on 〈1 0 0〉 symmetric tilt grain boundaries in tungsten
journal, July 2018
- Yang, Zhuocen; Hammond, Karl D.
- Journal of Physics: Condensed Matter, Vol. 30, Issue 32
Dynamics of small mobile helium clusters near tungsten surfaces
journal, August 2014
- Hu, Lin; Hammond, Karl D.; Wirth, Brian D.
- Surface Science, Vol. 626
Formation of Nanostructured Tungsten with Arborescent Shape due to Helium Plasma Irradiation
journal, January 2006
- Takamura, Shuichi; Ohno, Noriyasu; Nishijima, Dai
- Plasma and Fusion Research, Vol. 1
Thermal stability of helium–vacancy clusters in iron
journal, April 2003
- Morishita, K.; Sugano, R.; Wirth, B. D.
- Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, Vol. 202
Effect of rhenium addition on fracture toughness of tungsten at elevated temperatures
journal, February 1995
- Mutoh, Y.; Ichikawa, K.; Nagata, K.
- Journal of Materials Science, Vol. 30, Issue 3
Interactions of mobile helium clusters with surfaces and grain boundaries of plasma-exposed tungsten
journal, May 2014
- Hu, Lin; Hammond, Karl D.; Wirth, Brian D.
- Journal of Applied Physics, Vol. 115, Issue 17
Tungsten ‘fuzz’ growth re-examined: the dependence on ion fluence in non-erosive and erosive helium plasma
journal, August 2015
- Petty, T. J.; Baldwin, M. J.; Hasan, M. I.
- Nuclear Fusion, Vol. 55, Issue 9
An improved N -body semi-empirical model for body-centred cubic transition metals
journal, July 1987
- Ackland, G. J.; Thetford, R.
- Philosophical Magazine A, Vol. 56, Issue 1
Crystal orientation effects on helium ion depth distributions and adatom formation processes in plasma-facing tungsten
journal, October 2014
- Hammond, Karl D.; Wirth, Brian D.
- Journal of Applied Physics, Vol. 116, Issue 14
Continuum-scale modeling of helium bubble bursting under plasma-exposed tungsten surfaces
journal, November 2018
- Blondel, Sophie; Bernholdt, David E.; Hammond, Karl D.
- Nuclear Fusion, Vol. 58, Issue 12
A review of recent experiments on W and high Z materials as plasma-facing components in magnetic fusion devices
journal, February 1997
- Noda, Nobuaki; Philipps, Volker; Neu, Rudolf
- Journal of Nuclear Materials, Vol. 241-243
A new interatomic potential function for helium
journal, January 1968
- Beck, D. E.
- Molecular Physics, Vol. 14, Issue 4
Challenges and opportunities of modeling plasma–surface interactions in tungsten using high-performance computing
journal, August 2015
- Wirth, Brian D.; Hammond, K. D.; Krasheninnikov, S. I.
- Journal of Nuclear Materials, Vol. 463
Strain-rate dependence of the brittle-to-ductile transition temperature in tungsten
journal, June 2007
- Giannattasio, A.; Roberts, S. G.
- Philosophical Magazine, Vol. 87, Issue 17
Formation process of tungsten nanostructure by the exposure to helium plasma under fusion relevant plasma conditions
journal, August 2009
- Kajita, Shin; Sakaguchi, Wataru; Ohno, Noriyasu
- Nuclear Fusion, Vol. 49, Issue 9
Large-scale atomistic simulations of low-energy helium implantation into tungsten single crystals
journal, February 2018
- Hammond, Karl D.; Blondel, Sophie; Hu, Lin
- Acta Materialia, Vol. 144
LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales
journal, February 2022
- Thompson, Aidan P.; Aktulga, H. Metin; Berger, Richard
- Computer Physics Communications, Vol. 271
Molecular dynamics simulations on the effect of sub-surface helium bubbles on the sputtering yield of tungsten
journal, July 2013
- Sefta, F.; Juslin, N.; Hammond, K. D.
- Journal of Nuclear Materials, Vol. 438
Fracture toughness of polycrystalline tungsten alloys
journal, November 2010
- Gludovatz, B.; Wurster, S.; Hoffmann, A.
- International Journal of Refractory Metals and Hard Materials, Vol. 28, Issue 6
A comparison of interatomic potentials for modeling tungsten–hydrogen–helium plasma–surface interactions
journal, August 2015
- Cusentino, Mary Alice; Hammond, Karl D.; Sefta, Faiza
- Journal of Nuclear Materials, Vol. 463
Fracture-mechanical properties of neutron irradiated ITER specification tungsten
journal, April 2021
- Gaganidze, E.; Chauhan, A.; Schneider, H. -C.
- Journal of Nuclear Materials, Vol. 547
Theoretical Model of Helium Bubble Growth and Density in Plasma-Facing Metals
journal, February 2020
- Hammond, Karl D.; Maroudas, Dimitrios; Wirth, Brian D.
- Scientific Reports, Vol. 10, Issue 1
Molecular-dynamics analysis of mobile helium cluster reactions near surfaces of plasma-exposed tungsten
journal, October 2015
- Hu, Lin; Hammond, Karl D.; Wirth, Brian D.
- Journal of Applied Physics, Vol. 118, Issue 16
On the origin of ‘fuzz’ formation in plasma-facing materials
journal, July 2019
- Dasgupta, Dwaipayan; Kolasinski, Robert D.; Friddle, Raymond W.
- Nuclear Fusion, Vol. 59, Issue 8
Interatomic potentials for simulation of He bubble formation in W
journal, January 2013
- Juslin, N.; Wirth, B. D.
- Journal of Nuclear Materials, Vol. 432, Issue 1-3
Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten
journal, October 2020
- Weerasinghe, Asanka; Hu, Lin; Hammond, Karl D.
- Journal of Applied Physics, Vol. 128, Issue 16
Limitations of W and W–1%La2O3 for use as structural materials
journal, June 2005
- Rieth, Michael; Dafferner, Bernhard
- Journal of Nuclear Materials, Vol. 342, Issue 1-3
Morphologies of tungsten nanotendrils grown under helium exposure
journal, February 2017
- Wang, Kun; Doerner, R. P.; Baldwin, M. J.
- Scientific Reports, Vol. 7, Issue 1
Anisotropic fracture behaviour and brittle-to-ductile transition of polycrystalline tungsten
journal, August 2010
- Rupp, D.; Weygand, S. M.
- Philosophical Magazine, Vol. 90, Issue 30
Competing Kinetics and He Bubble Morphology in W
journal, March 2015
- Sandoval, Luis; Perez, Danny; Uberuaga, Blas P.
- Physical Review Letters, Vol. 114, Issue 10
Ductile to brittle transition in ITER specification tungsten assessed by combined fracture toughness and bending tests analysis
journal, March 2019
- Yin, Chao; Terentyev, Dmitry; Pardoen, Thomas
- Materials Science and Engineering: A, Vol. 750
Helium induced nanoscopic morphology on tungsten under fusion relevant plasma conditions
journal, January 2008
- Baldwin, M. J.; Doerner, R. P.
- Nuclear Fusion, Vol. 48, Issue 3
An improved xenon equation of state for nanobubbles in UO2
journal, December 2022
- Yang, L.; Wirth, B. D.
- Journal of Nuclear Materials, Vol. 572