1 / 49

Toshimasa Yoshiie Research Reactor Institute, Kyoto University

Nuclear Data and Materials Irradiation Effects - Analysis of irradiation damage structures and multiscale modeling -. Toshimasa Yoshiie Research Reactor Institute, Kyoto University. Comparison of irradiation effects between different facilities. Power reactors  Research reactors

mercia
Download Presentation

Toshimasa Yoshiie Research Reactor Institute, Kyoto University

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Nuclear Data and Materials Irradiation Effects- Analysis of irradiation damage structures and multiscale modeling - Toshimasa Yoshiie Research Reactor Institute, Kyoto University

  2. Comparison of irradiation effects between different facilities • Power reactors  Research reactors • Neutron irradiation Ion irradiation Electron irradiation DPA (Displacement per atom) The number of displacement of one atom during irradiation DPA dose not represent the effect of cascades

  3. Estimation of irradiation damage • Ion irradiation100MeVHe,1018ions/cm2 • Neutron irradiation1018n/cm2( ＞1MeV or ＞0.1MeV) • Displacement par atom (dpa) αｘ deposited energy in crystal lattice 2ｘ threshold energy of atomic displacement Kinchin-Peasemodel = DPA is the number of displacement for 1 atom during irradiation EP E T

  4. Frenkel pair formation Cascade formation Cascade

  5. Example of MD Simulation 40keV cascade of iron at 100K

  6. Clustering of Point Defects Stress field

  7. Interstitial Type Dislocation Loops in FCC Interstitial type dislocation loops in Al

  8. Voids Embrittlement Materials with no voids Materials with voids

  9. Stacking Fault Tetrahedra

  10. Stacking Fault Tetrahedra

  11. Cascade Damage High energy Particle Cascade Vacancy rich area Subcascades Nucleation of defect Interstitial rich area clusters

  12. Comparison of Subcascade Structuresby Thin Foil Irradiation 14MeV neutrons at 300K. Neutron PKA

  13. Au, KENS irradiation and 14 MeV irradiation at room temperature KENS neutrons 14MeVneutrons

  14. 14MeV neutron irradiated Cu

  15. 14MeV neutron irradiated Cu

  16. 14MeV neutron irradiated Cu

  17. 14MeV neutron irradiated Cu

  18. Temperature Effects of Cascade Damage Lower temperature Higher temperature Subcascades Subcascades fuse into a large cluster

  19. Thin foil Irradiated Au by Fusion Neutrons A group of SFTs A large SFT Subcascade Cascade fuses structures into one SFT

  20. 0.017dpa Large SFTs High number density of SFTs 0.044dpa Small SFTs Low number density of SFTs Fission-fusion Correlation of SFT in Au Fusion neutrons Fission neutrons 563K 573K

  21. NSFT: the concentration of SFTs observed α: the SFT formation efficiency ：the neutron fluence 　：the differential cross- section for PKA ETH :threshold energy for SFT formation PKA Energy Spectrum Analysis α= 0.05 ETH = 80keV

  22. MULTI-SCALE MODELING OF IRRADIATION EFFECTS IN SPALLATION NEUTRON SOURCE MATERIALS Toshimasa Yoshiie1, Takahiro Ito 2, Hiroshi Iwase 3, Yoshihisa Kaneko4, Masayoshi Kawai3, Ippei Kishida4, Satoshi Kunieda5, Koichi Sato1, Satoshi Shimakawa5, Futoshi Shimizu5, Satoshi Hashimoto4, Naoyuki Hashimoto6, Tokio Fukahori5, Yukinobu Watanabe7, Qiu Xu1, Shiori Ishino8 1Research Reactor Institute, Kyoto University 2 Department of Mechanical Engineering, Toyohashi University of Technology High Energy Accelerator Research Organization Osaka City University 5Japan Atomic Energy Agency 6Hokkaido University 7 Kyushu University 8Univerity of Tokyo

  23. Motivation • Spallation neutron source and Accelerator Driven System (ADS) are a coupling of a target and a proton accelerator. High energy protons of GeV order irradiated in the target produce a large number of neutrons. • The beam window and the target materials thus subjected to a very high irradiation load by source protons and spallation neutrons generated inside the target. • At present, there are no materials that enable the window to be operational for the desired period of time without deterioration of mechanical properties.

  24. Importance of GeV Order Proton Irradiation Effects in Materials • Spallation neutron source J-PARC (Japan) SNS (USA) • Accelerator driven system (in the planning stage ) 800MW ADS (Minor Actinide transmutation, Japan Atomic Energy Agency) 5MW Accelerator driven subcritical reactor (Kyoto University, neutron source)

  25. Purpose and outline • In this paper, mechanical property changes of nickel by 3 GeV protons were calculated by multi-scale modeling of irradiation effects. Nickel is considered to be a most simple model material of austenitic stainless steels used in beam window. The code consists of four parts. • Nuclear reaction, the interaction between high energy protons and nuclei in the target is calculated by PHITS code from 10-22 s. • Atomic collision by particles which do not cause nuclear reactions is calculated by molecular dynamics and k-Monte Carlo. As the energy of particles is high, subcascade analysis is employed. In each subcascade, the direct formation of clusters and the number of mobile defects are estimated. • Damage structure evolution is estimated by reaction kinetic analysis. • Mechanical property change is calculated by using 3D discrete dislocation dynamics (DDD). Stress-strain curves of high energy proton irradiated nickel are obtained.

  26. High energy particle Vacancy Size (m) PHITS 10-2010-1510-1010-51001051010

  27. Data flow between each code Nuclear reactions （PHITS code） Primary knock-on energy spectrum Atomic collisions （Molecular dynamics) Point defect distribution Damage structural evolutions (Reaction kinetic analysis) Concentration of defect clusters Mechanical property change (Three-dimensional discrete dislocation dynamics) Stress-strain curve

  28. 1. Nuclear Reaction Nucleation rate of neutrons, photons, charged particles and PKA energy spectrum by them

  29. Result of PHITS Simulation 3GeV protons Ni Number of PKA energy (left) and energy deposition by PKA (right) in 3 GeV proton irradiated Ni of 3 mm in thickness.

  30. Subcascade Analysis The number of subcascasdes are obtained from the result of PHITS. Large cascades are divided into subcascades. In the case of Ni, subcascade formation energy is calculated to be 10 keV. Cascade

  31. Number of subcascades by deposition of energy T TSC : Subcascade formation energy Total number of subcascades

  32. 2. Atomic Collision Simulation by MD Calculation Condition • Potential model by Daw and Baskes (1984) • NVE ensemble ( i.e., number of the atoms, cell volume and energy were kept constant) • 35x35x35 lattices (171500 atoms) in a simulation cell • Periodic boundary condition for the three directions • Initial condition : equilibrium for 50 ps at 300 K, 0MPa • PKA energy：10keV • MD runs with different initial directions of PKA ( none of which were parallel to the lattice vector.) z y x 35a0 a0 : lattice constant at 300 K

  33. Typical Distribution of Point Defects Marble : interstitial atoms，Violet ：Vacancy cites 0.006ps 0.025ps 1.132ps 18.40ps 96.38ps 7.395ps

  34. Results of MD Calculation Number of vacancies • Cascades terminated within 10~20ps • On average17 vacancies, 17 interstitials were produced. • Formation of defect clusters is calculated by k-Monte Carlo. • Clusters of three point defects are formed.

  35. 3.Damage Structure Evolutionby Reaction Kinetic Analysis In order to estimate damage structural evolution, the reaction kinetic analysis is used. Assumptions used in the calculation are as follows: (1) Mobile defects are interstitials, di-interstitials, tri-interstitials, vacancies and di-vacancies. (2) Thermal dissociation is considered for di-interstitials, tri-interstitials di-vacancies and tri-vacancies, and point defect clusters larger than 4 are set for stable clusters. (3) Time dependence of 10 variables, concentration of interstitials, di-interstitials, tri-interstitials, interstitial clusters (interstitial type dislocation loops), vacancies, di-vacancies, tri-vacancies, vacancy clusters (voids), total interstitials in interstitial clusters and total vacancies in vacancy clusters are calculated to 10 dpa. (4) Interstitial clusters (three interstitials) and vacancy clusters (three vacancies) are also formed directly in subcascades. (5) Materials temperature is 423 K during irradiation. The result is as follows, Formation of vacancy clusters of four vacancies, concentration: 0.59x10-3, Dislocation density: 1.1x10-9cm/cm3.

  36. damage productionI-I recombinationmutual annihilation absorption by loopsabsorption by voids Reaction Kinetic Analysis CI : Interstitial concentration (fractional unit). CV : Vacancy concentration Z : Cross section of reaction M: Mobility annihilation of interstitials at sink

  37. 4. Estimation of Mechanical Property Changes (Plastic Deformation of Metals) t t Plastic Deformation of Metals Dislocation Glide along Slip Planes s Work Hardening ↓ Motions of Dislocation Lines e

  38. Calculation of Dislocation Motions • Motions of Curved Dislocations • External Stress • Elastic Interaction between • Dislocations • Line Tension Division of a Dislocation Line ↓ Discrete Dislocation Dynamics (DDD) Simulation Very Complicated !

  39. b 3D-Discrete Dislocation Dynamics • Edge+Screw Dislocation Model • Devincre (1992～) External Stress Stresses from another segment • Mixed Dislocation Model • Zbib, et al (1998～) • Schwarz et al • Fivel et al ，Caiet al Line tension Peach-Koehler Force Dislocation Velocity Connection of Dislocation Segments with Mixed Characters

  40. Model Crystal and Mobile Dislocation Ni: shear modulus = 76.9 Gpa Poisson’s ratio = 0.31 Burgers vector = 0.24916 nm Mobile dislocation Stress axis ＝ [1 0 0] slip plane＝(1 1 1) slip vector＝[1 0 1] Obstacles（void） 0.593×10-3 void/ atom 0.14μm the densityalong the slip plane : 1.10x104 void / μm2 0.14μm 0.14μm

  41. Interaction between Dislocation and void Detrapping angle

  42. Determination of Detrapping Angle (Statistic Energy Calculation Method) [ l l l ] ( l l l) (Ī l 0) [ Ī Ī 2] void [ Ī l 0] (Ī Ī 2) Stress Size and orientation of model lattice used for Static energy calculation of dislocation movement. b is the Burgers vector. Stable atomic position is determined by an effective medium theory (EMT) potential for Ni fitted by Jacobsen et al. [K. W. Jacobsen, P. Stoltze, J. K. Norskov, Surf. Sci. 366 (1996) 394].

  43. Determination of Detrapping Angle between Edge Dislocations and 4 Vacancies 65° • Motion of dislocation near the slip plane under shear stress. X axis and Y axis are by atomic distance. In the left figure, a dislocation is separated into two partials. Four vacancies are at (16.0, 0.722, -0.408), (15.5, -0.144, -0.408), (16.0, 0.144, 0.408) and (15.5, 1.01, 0.408). Other figures indicate only right partial.

  44. Dislocation Motion by DDD Simulation (no void) 700MPa 400MPa 600MPa

  45. Dislocation motion in the crystal with randomly-distributed voids 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa 0MPa 200MPa 300MPa 400MPa 500MPa 600MPa 650MPa 700MPa 750MPa 800MPa 820MPa 840MPa 860MPa 880MPa 900MPa 920MPa 930MPa 935MPa Normal stress

  46. Stress-Strain Curves Calculated with DDD Simulation Plastic shear strain gp = A b / V Plastic strain ep = gp Sf A: area swept by a dislocation b: Burgers vector V: volume of model crystal Sf: Schmid factor (x10-6)

  47. Summary Importance of nuclear data for materials irradiation effects was shown. DPA is not good measure of irradiation damage. PKA energy spectrum is more useful to analyze damage structures. As an example of the analysis, the mechanical property change in Ni by 3 GeV proton irradiation was calculated and preliminary results were obtained. Nuclear reactions : PHITS PKA energy spectrum Atomic collision : MD, k-Monte Carlo, Subcascade analysis. The number of point defects and clusters in the subcascade Damage structure evolution : Reaction kinetic analysis Void density, Dislocation density Mechanical properties : Discrete dislocation dynamics Statistic energy calculation Stress strain curve

  48. Important data for materials irradiation effects High energy particles Primary knock on atom energy spectrum Formation rate of atoms by nuclear reactions High energy particle energy spectrum

More Related