Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Yamaguchi, Masatake; Ebihara, Kenichi; Itakura, Mitsuhiro; Tsuru, Tomohito
Scripta Materialia, 255, p.116366_1 - 116366_5, 2025/01
Times Cited Count:0 Percentile:0.00(Nanoscience & Nanotechnology)Reduction of grain boundary cohesive energy by hydrogen has been considered as one of the candidate causes of grain boundary fracture in steels and aluminum alloys. Recently, the effects of grain boundary segregation elements have been investigated by first-principles calculations, but there are few studies that quantitatively estimate the grain boundary cohesive energy. In this study, we describe a quantitative evaluation method using first-principles calculation results and show some examples of experimental calculations.
Okada, Kazuho*; Shibata, Akinobu*; Kimura, Yuji*; Yamaguchi, Masatake; Ebihara, Kenichi; Tsuji, Nobuhiro*
Acta Materialia, 280, p.120288_1 - 120288_14, 2024/11
Times Cited Count:0 Percentile:0.00(Materials Science, Multidisciplinary)Higa, Ryota*; Fujihara, Hiro*; Toda, Hiroyuki*; Kobayashi, Masakazu*; Ebihara, Kenichi; Takeuchi, Akihisa*
Materials Transactions, 65(8), p.899 - 906, 2024/08
Times Cited Count:0 Percentile:0.00(Materials Science, Multidisciplinary)It is indispensable to suppress hydrogen embrittlement (HE) to develop the strength of the Al-Zn-Mg alloy. Because intergranular fracture (IGF) is mainly observed when HE occurs in the alloy, we need to understand the IGF initiation to suppress HE. In the present study, we investigated the stress, strain, and H concentration, which influence the IGF initiation, in actual fractured regions by simulation of a crystal plasticity finite element method and H diffusion analysis in a 3D image-based model, which was created based on 3D polycrystalline microstructure data obtained from X-ray imaging technique. Combining the simulation and in-situ observation of the tensile test sample by X-ray CT, we examined the stress, strain, and H concentration, and discussed the IG crack initiation condition. As a result, it is revealed that stress normal to grain boundary induced by crystal plasticity dominates IG crack initiation while the accumulation of H due to stress has little impact on it.
Fujihara, Hiro*; Toda, Hiroyuki*; Ebihara, Kenichi; Kobayashi, Masakazu*; Mayama, Tsuyoshi*; Hirayama, Kyosuke*; Shimizu, Kazuyuki*; Takeuchi, Akihisa*; Uesugi, Masayuki*
International Journal of Plasticity, 174, p.103897_1 - 103897_22, 2024/03
Times Cited Count:7 Percentile:95.75(Engineering, Mechanical)Hydrogen(H) embrittlement in high-strength aluminum(Al) alloys is a crucial problem. H accumulation at the interface of precipitates in Al alloy is considered to cause embrittlement. However, there is no quantitative knowledge regarding the interaction between H distribution and stress field near cracks. In this study, using a multi-modal three-dimensional image-based simulation combining the crystal plasticity finite element method and H diffusion analysis, we tried to capture the stress distribution near the crack, its influence on the H distribution, and the probability of crack initiation in the experimental condition. As a result, it was found that grain boundary cracks transition to quasi-cleavage cracks in the region where the cohesive energy of the semi-coherent interface of MgZn precipitates decreases due to H accumulation near the tip. We believe the present simulation method successfully bridges nanoscale delamination and macroscale brittle fracture.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Journal of Applied Physics, 135(7), p.075102_1 - 075102_7, 2024/02
Times Cited Count:0 Percentile:0.00(Physics, Applied)Fracture of body centred cubic (bcc) metals and alloys below the ductile-to-brittle transition temperature is brittle. This is theoretically explained by the notion that the critical stress intensity factor of a given crack front for brittle fracture is smaller than that for plasticdeformation; hence, brittle fracture is chosen over plastic deformation. Although this view is true from a macroscopic point of view, such brittle fracture is always accompanied by small-scale plastic deformation in the vicinity of the crack tip, i.e. crack tip plasticity. This short paper investigates the origin of this plasticity using atomistic modeling with a recently developed machine-learning interatomic potential of -Fe. The computational results identified the precursor of crack tip plasticity, i.e. the group of activated atoms dynamically nucleated by fast crack propagation.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Zairyo, 73(2), p.129 - 135, 2024/02
Body-centered-cubic transition metals, such as Fe and W, cleave along the {100} plane. To find out the mechanism of this response, atomistic simulations of curved crack-fronts of bcc Fe were conducted at 0 K using an interatomic potential created by an artificial neural network (ANN) technique. We discovered that dislocations can be emitted from the curved crack fronts along the {110} crack plane, and this phenomenon explains why the cleavage is observed only along the {100} plane. In addition, the cleavage simulations along {100} at the elevated temperature were found to be accompanied by plasticity; namely, they represented more realistic fracture.
Tang, J.*; Wang, Y.*; Fujihara, Hiro*; Shimizu, Kazuyuki*; Hirayama, Kyosuke*; Ebihara, Kenichi; Takeuchi, Akihisa*; Uesugi, Masayuki*; Toda, Hiroyuki*
Scripta Materialia, 239, p.115804_1 - 115804_5, 2024/01
Times Cited Count:9 Percentile:72.83(Nanoscience & Nanotechnology)Stress corrosion cracking (SCC) behaviors induced by the combination of external and internal hydrogen (H) in an Al-Zn-Mg-Cu alloy were systematically investigated via in situ 3D characterization techniques. SCC of the Al-Zn-Mg-Cu alloy could initiate and propagate in the potential crack region where the H concentration exceeded a critical value, in which the nanoscopic H-induced decohesion of -MgZn
precipitates resulted in macroscopic cracking. External H that penetrated the alloy from the environment played a crucial role during the SCC of the Al-Zn-Mg-Cu alloy by generating gradient-distributed H-affected zones near the crack tips, which made Al alloys in water environment more sensitive to SCC. Additionally, the pre-existing internal H was driven toward the crack tips during plastic deformation. It was involved in the SCC and made contributions to both the cracks initiation and propagation.
Higa, Ryota*; Fujihara, Hiro*; Toda, Hiroyuki*; Kobayashi, Masakazu*; Ebihara, Kenichi; Takeuchi, Akihisa*
Keikinzoku, 73(11), p.530 - 536, 2023/11
In Al-Zn-Mg alloys, suppression of hydrogen embrittlement is necessary to improve their strength. In this study, the distribution of stress, strain, and hydrogen concentration in the actual fracture region was investigated using the crystal plasticity finite element method and hydrogen diffusion analysis based on a model derived from three-dimensional polycrystalline microstructural data obtained by X-ray CT. In addition, the distributions of stress, strain, and hydrogen concentration were compared with the actual crack initiation behavior by combining in-situ observation of tensile tests using X-ray CT and simulation. The results show that stress loading perpendicular to the grain boundary due to crystal plasticity dominates grain boundary crack initiation. It was also found that internal hydrogen accumulation due to crystal plasticity has little effect on crack initiation.
Yamaguchi, Masatake; Ebihara, Kenichi; Tsuru, Tomohito; Itakura, Mitsuhiro
Materials Transactions, 64(11), p.2553 - 2559, 2023/11
Times Cited Count:5 Percentile:58.10(Materials Science, Multidisciplinary)We attempted to calculate the hydrogen trapping energies on the incoherent interfaces of MgZn precipitates and Mg
Si crystallites in aluminum alloys from first-principles calculations. Since the unit cell containing the incoherent interface does not satisfy the periodic boundary condition, resulting in a discontinuity of crystal blocks, the hydrogen trapping energy was calculated in a region far from the discontinuity (vacuum) region. We found considerable trapping energies for hydrogen atoms at the incoherent interfaces consisting of assumed atomistic arrangement. We also conducted preliminary calculations of the reduction in the cohesive energy by hydrogen trapping on the incoherent interfaces of Mg
Si in the aluminum matrix.
Ebihara, Kenichi; Sekine, Daiki*; Sakiyama, Yuji*; Takahashi, Jun*; Takai, Kenichi*; Omura, Tomohiko*
International Journal of Hydrogen Energy, 48(79), p.30949 - 30962, 2023/09
Times Cited Count:0 Percentile:0.00(Chemistry, Physical)To understand hydrogen embrittlement (HE), which is one of the stress corrosion cracking of steel materials, it is necessary to know the H distribution in steel, which can be effectively interpreted by numerical simulation of thermal desorption spectra. In weld metals and TRIP steels, residual austenite significantly influences the spectra, but a clear H distribution is not well known. In this study, an originally coded two-dimensional model was used to numerically simulate the previously reported spectra of high-carbon ferritic-austenitic duplex stainless steels, and it was found that H is mainly trapped at the carbide surface when the amount of H in the steel is low and at the duplex interface when the amount of H is high. It was also found that the thickness dependence of the H desorption peak for the interface trap site is caused by a different reason than the conventional one.
Tsuru, Tomohito; Nishimura, Katsuhiko*; Matsuda, Kenji*; Nunomura, Norio*; Namiki, Takahiro*; Lee, S.*; Higemoto, Wataru; Matsuzaki, Teiichiro*; Yamaguchi, Masatake; Ebihara, Kenichi; et al.
Metallurgical and Materials Transactions A, 54(6), p.2374 - 2383, 2023/06
Times Cited Count:2 Percentile:33.88(Materials Science, Multidisciplinary)Although hydrogen embrittlement susceptibility of high-strength Al alloys is recognized as a critical issue in the practical use of Al alloys, identifying the hydrogen trapping or distribution has been challenging. In the present study, an effective approach based on experiment and simulation is proposed to explore the potential trap sites in Al alloys. Zero-field muon spin relaxation experiments were carried out for Al-0.5%Mg, Al-0.2%Cu, Al-0.15%Ti, Al-0.011%Ti, Al-0.28%V, and Al-0.015%V (at.%) in the temperature range from 5 to 300 K. The temperature variations of the dipole field widths have revealed three peaks for Al-0.5%Mg, four peaks for Al-0.2%Cu, three peaks for Al-0.011%Ti and Al-0.015%V. Atomic configurations of the muon trapping sites corresponding to the observed peaks are well assigned using the first-principles calculations for the trap energies of hydrogen around a solute and solute-vacancy pair. The extracted linear relationship between the muon
peak temperature and the trap energy enables us to explore the potential alloying elements and their complex that have strong binding energies with hydrogen in Al alloys.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Scientific Reports (Internet), 12, p.19701_1 - 19701_10, 2022/11
Times Cited Count:9 Percentile:56.77(Multidisciplinary Sciences)Body-centered-cubic (bcc) transition metals, such as -Fe and W, cleave along the {100} plane, even though the surface energy is the lowest along the {110} plane. To unravel the mechanism of this odd response, large-scale atomistic simulations of curved cleavage cracks of
-Fe were conducted in association with stress intensity factor analyses of straight crack fronts using an interatomic potential created by an artificial neural network technique. The study provides novel findings: Dislocations are emitted from the crack fronts along the {110} cleavage plane, and this phenomenon explains why the {100} plane can be the cleavage plane. However, the simple straight crack-front analyses did not yield the same conclusion. It is suggested that atomistic modeling, at sufficiently large scales to capture the inherent complexities of materials using highly accurate potentials, is necessary to correctly predict the mechanical strength. The method adopted in this study is generally applicable to the cleavage problem of bcc transition metals and alloys.
Ebihara, Kenichi
Zairyo, 71(5), p.481 - 487, 2022/05
no abstracts in English
Ebihara, Kenichi; Suzudo, Tomoaki
Metals, 12(4), p.662_1 - 662_10, 2022/04
Times Cited Count:4 Percentile:36.22(Materials Science, Multidisciplinary)Phosphorus atoms in steels accumulate at grain boundaries via thermal and/or irradiation effects and induce grain boundary embrittlement. Quantitative prediction of phosphorus segregation at grain boundaries under various temperature and irradiation conditions is therefore essential for preventing embrittlement. To develop a model of grain boundary phosphorus segregation in -iron, we studied the migration of a phosphorus atom in two types of symmetrical tilt grain boundaries (
3[1-10](111) and
5[100](0-13) grain boundaries) using molecular dynamics simulations with an embedded atom method potential. The results revealed that, in the
3 grain boundary, phosphorus atoms migrate three-dimensionally mainly in the form of interstitial atoms, whereas in the
5 grain boundary, these atoms migrate one-dimensionally mainly via vacancy-atom exchanges. Moreover, de-trapping of phosphorus atoms and vacancies was investigated.
Yamaguchi, Masatake; Itakura, Mitsuhiro; Tsuru, Tomohito; Ebihara, Kenichi
Materials Transactions, 62(5), p.582 - 589, 2021/05
Times Cited Count:16 Percentile:71.72(Materials Science, Multidisciplinary)no abstracts in English
Saito, Shimpei*; De Rosis, A.*; Fei, L.*; Luo, K. H.*; Ebihara, Kenichi; Kaneko, Akiko*; Abe, Yutaka*
Physics of Fluids, 33(2), p.023307_1 - 023307_21, 2021/02
Times Cited Count:43 Percentile:98.23(Mechanics)A Boiling phenomenon in a liquid flow field is known as forced-convection boiling. We numerically investigated the boiling system on a cylinder in a flow at a saturated condition. To deal with such a phenomenon, we developed a numerical scheme based on the pseudopotential lattice Boltzmann method. The collision was performed in the space of central moments (CMs) to enhance stability for high Reynolds numbers. Furthermore, additional terms for thermodynamic consistency were derived in a CMs framework. The effectiveness of the model was tested against some boiling processes, including nucleation, growth, and departure of a vapor bubble for high Reynolds numbers. Our model can reproduce all the boiling regimes without the artificial initial vapor phase. We found that the Nukiyama curve appears even though the focused system is the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid-liquid contact even in the film-boiling regime.
Hata, Kuniki; Takamizawa, Hisashi; Hojo, Tomohiro*; Ebihara, Kenichi; Nishiyama, Yutaka; Nagai, Yasuyoshi*
Journal of Nuclear Materials, 543, p.152564_1 - 152564_10, 2021/01
Times Cited Count:16 Percentile:90.59(Materials Science, Multidisciplinary)Reactor pressure vessel (RPV) steels for pressurized water reactors (PWRs) with bulk P contents ranging from 0.007 to 0.012wt.% were subjected to neutron irradiation at fluences ranging from 0.3 to 1.210
n/cm
(E
1 MeV) in PWRs or a materials testing reactor (MTR). Grain-boundary P segregation was analyzed using Auger electron spectroscopy (AES) on intergranular facets and found to increase with increasing neutron fluence. A rate theory model was also used to simulate the increase in grain-boundary P segregation for RPV steels with a bulk P content up to 0.020wt.%. The increase in grain-boundary P segregation in RPV steel with a bulk P content of 0.015wt.% (the maximum P concentration found in RPV steels used in Japanese nuclear power plants intended for restart) was estimated to be less than 0.1 in monolayer coverage at 1.0
10
n/cm
(E
1 MeV). A comparison of the PWR data with the MTR data showed that neutron flux had no effect upon grain-boundary P segregation. The effects of grain-boundary P segregation upon changes in irradiation hardening and ductile-brittle transition temperature (DBTT) shifts were also discussed. A linear relationship between irradiation hardening and the DBTT shift with a slope of 0.63 obtained for RPV steels with a bulk P content up to 0.026wt.%, which is higher than that of most U.S. A533B steels. It is concluded that the intergranular embrittlement is unlikely to occur for RPV steels irradiated in PWRs.
Ebihara, Kenichi; Sugiyama, Yuri*; Matsumoto, Ryosuke*; Takai, Kenichi*; Suzudo, Tomoaki
Metallurgical and Materials Transactions A, 52(1), p.257 - 269, 2021/01
Times Cited Count:9 Percentile:43.14(Materials Science, Multidisciplinary)We simulated the thermal desorption spectra of a small-size iron specimen to which was applied during charging with hydrogen atoms using a model incorporating the behavior of vacancies and vacancy clusters. The model considered up to vacancy clusters , which is composed of nine vacancies and employed the parameters based on atomistic calculations, including the H trapping energy of vacancies and vacancy clusters that we estimated using the molecular static calculation. As a result, we revealed that the model could, on the whole, reproduced the experimental spectra except two characteristic differences, and also the dependence of the spectra on the aging temperature. By examining the cause of the differences, the possibilities that the diffusion of clusters of
and
is slower than the model and that vacancy clusters are generated by applying strain and H charging concurrently were indicated.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito
AIP Advances (Internet), 10(11), p.115209_1 - 115209_8, 2020/11
Times Cited Count:14 Percentile:61.66(Nanoscience & Nanotechnology)The mechanism of their brittle fracture of BCC metals is not fully understood. In this study, we conduct a series of three-dimensional molecular dynamics simulations of cleavage fracture of -iron. In particular, we focus on mode-I loading starting from curved crack fronts. In the simulations, brittle fractures are observed at cleavages on the {100} plane, while the initial cracks become blunted on other planes as a result of dislocation emissions. Our modeling results agreed with a common experimental observation, that is, {100} is the preferential cleavage plane in bcc transition metals.
Yamaguchi, Masatake; Tsuru, Tomohito; Ebihara, Kenichi; Itakura, Mitsuhiro; Matsuda, Kenji*; Shimizu, Kazuyuki*; Toda, Hiroyuki*
Materials Transactions, 61(10), p.1907 - 1911, 2020/10
Times Cited Count:18 Percentile:66.74(Materials Science, Multidisciplinary)no abstracts in English