Finite element modeling of hydrogen atom diffusion and distribution at corrosion defect on aged pipelines transporting hydrogen

Repurposing existing natural gas pipelines for hydrogen transport has attracted wide interests. However, the corrosion defect present on these aged pipelines can affect hydrogen (H) atom accumulation, potentially causing hydrogen embrittlement. In this work, a finite element-based model was developed by coupling solid mechanics and H atom diffusion to investigate the distribution of H atoms at a corrosion defect on a steel pipe segment under applied longitudinal tensile strains. The applied strain causes local stress (both Mises stress and hydrostatic stress) and strain concentrations at the corrosion defect, affecting the H atom diffusion and distribution. In the absence of the tensile strain, the H atoms, once entering the interior of pipe, diffuse uniformly into the pipe body along the radial direction driven by a concentration gradient. When a strain is applied on the pipe, the H atom diffusion is driven by hydrostatic stress. The maximum H atom concentration exceeds the initial concentration of H atoms entering the steel pipe, indicating the H atom accumulation at the corrosion defect. The applied tensile strain also affects the location where the H atoms accumulate. For both internal and external corrosion defects, more H atoms will be concentrated at the defect center when the defect length reduces and the depth increases.     

Evolution of microstructure and texture in pipeline steels at different TMCP procedures with regard to hydrogen induced cracking

SEM and EBSD techniques are used to evaluate hydrogen induced cracking susceptibility in API X70 pipeline steels produced by thermo-mechanical controlled process (TMCP) in laboratory scale. Based on the observations, there is no dominant texture in the specimens and the grains are randomly distributed. Different TMCP parameters and rolling processes generates different grain size, and grains are often elongated along the rolling direction. The results also show that cooling rate is another factor affecting the grain size. A high cooling rate does not allow the grains to grow. The reason for the transgranular type of cracking might be the strong grain boundaries in ambient temperatures which prevents the intergranular cracking. Based on experiments, the hydrogen environment does not have permanent effects on the mechanical properties of the investigated specimens. The electrochemical hydrogen charging experiment shows that the grain refinement improves the resistance to hydrogen embrittlement.     

Finite element analysis of hydrogen diffusion/plasticity coupled behaviors of low-alloy ferritic steel at large strain

Hydrogen embrittlement is commonly considered as a typical failure mechanism for low-alloy ferritic steel under high pressure hydrogen environment. Currently, the hydrogen enhanced localized plasticity theory has been largely recognized for studying the hydrogen embrittlement mechanism by introducing the localized plastic flow and the hydrogen induced strain concept. However, the hydrogen induced strain and the plastic strain are often solved respectively in this theory, which may weaken the effect of hydrogen on the plastic deformation. The purpose of this paper is to propose a modified theoretical model from the microstructural level by emphasizing the coupling mechanism between the hydrogen diffusion and the plastic deformation at large strain, where the hydrogen induced strain is superimposed on the equivalent plastic strain instead of on the strain components. Fully implicit backward Euler algorithm by finite element analysis (FEA) under the corotational configuration is used to implement the proposed model, where the hydrogen induced strain is involved in the stress return process within each iteration, indicating a more direct interaction between them than existing works. FEA by using finite element software ABAQUS-UMAT subroutine is performed for the smooth tensile specimen and the notch specimen respectively under slow tensile strain rate loading and different hydrogen pressure. Developed direct coupling model is expected to further gain insight into the hydrogen embrittlement effect on the plastic deformation, especially at the trapping sites.     

Effects of vanadium precipitates on hydrogen trapping efficiency and hydrogen induced cracking resistance in X80 pipeline steel

In this study, the number and size distribution of vanadium precipitates and their effects on hydrogen trapping efficiency and hydrogen-induced cracking (HIC) susceptibility were investigated in X80 pipeline steel. The results showed that as the vanadium content increased, the number of nanoscale vanadium precipitates clearly increased. Furthermore, the amount of hydrogen atoms trapped by vanadium precipitates gradually increased and the hydrogen diffusion coefficient decreased from 4.74 × 10^?6 cm² s^?1 in the vanadium-free V0 steel to 8.48 × 10^?7 cm² s^?1 in the V4 steel with 0.16% V, according to hydrogen permeation results. It also reduced the possibility of hydrogen atoms diffusing into the sites of harmful defects such as large-size oxides and elongated MnS inclusions, where cracks were caused more easily. In addition, the V3 steel with 0.12% V, containing the largest number of vanadium carbide particles of less than 60 nm, had the lowest HIC susceptibility.     

Fracture toughness properties of HIC susceptible carbon steels in sour service conditions

Hydrogen Induced Cracking (HIC) in carbon steels is a well-studied mechanism, where diffusing hydrogen atoms accumulates at the steel imperfections/laminations to create gaseous hydrogen with very high pressure, leading to initiation and growth of internal cavities, so-called HIC. Measurements of relevant fracture toughness properties of non-HIC resistant steels in hydrogen environment is critical to predict and assess the initiation and growth of HIC. The present work attempts to quantify the effect of hydrogen on the fracture toughness properties (K_Q and CTOD) of an API X42 pipeline steel under simulated H₂S in-service conditions. The fracture toughness properties are measured in TL and SL directions: perpendicular and parallel to the pipeline wall thickness, respectively, following ASTM E1820, standard. Since the X42 is a non-HIC resistant steel, the measurement of the fracture toughness properties in the SL direction is more relevant in terms of HIC initiation and growth than fracture toughness properties in the TL direction. Indeed, parallel to the thickness of the pipeline wall, X42 steel shows microstructural features prone to HIC formation and growth. Steady state H₂S in-service conditions were simulated by charging the specimen for 48 h using a special electrolytic solution and then tested (ex-situ) to evaluate the fracture toughness properties. The steady state H₂S environment was obtained by measuring the Hydrogen Concentration (CH) in the bulk of the specimen, using Thermal desorption Spectroscopy at three levels of CH. It was observed that the K_Q was not affected in the SL direction, while it was reduced in the TL direction for 1.5 ppmw of CH. The CTOD showed mixed results in the TL direction while it was significantly reduced in the SL direction reaching a saturation at 1 ppmw of CH. Besides, microstructural analyses showed that the presence of inclusions coalescence in form of dimples promote the early failure, which is more pronounced in the hydrogen environment especially at higher levels of CH.     

Effect of Ce content on the hydrogen induced cracking of X80 pipeline steel

In this study, the effect of Ce content on hydrogen induced cracking (HIC) of X80 pipeline steel has been investigated. The results show that as the Ce content increased from 0 wt% to 0.0042 wt%, 0.016 wt% and 0.024 wt%, the HIC susceptibility of tested steels decreased first and then increased. The steel containing 0.016 wt% Ce possessed the lowest HIC susceptibility because Ce modified inclusions, promoted the formation of acicular ferrite, and decreased the number of hydrogen traps and intergranular cracks.     

Hydrogen trapping and hydrogen induced cracking of welded X100 pipeline steel in H2S environments

Hydrogen induced cracking (HIC) susceptibility of the welded X100 pipeline steel was evaluated in NACE “A” solution at room temperature according to the NACE TM0284-2011 standard. Both the kinetic parameters of the permeability $\left(J_{\infty }L\right)$, the apparent diffusivity (D_app) and the concentration of reversible and irreversible hydrogen in the base metal and welded joint of X100 pipeline steel were quantitatively investigated by hydrogen permeation test. The results showed that the welded joint with an inhomogeneous microstructure had a higher trap density and more susceptible to HIC due to two orders of magnitude larger in the concentration of irreversible hydrogen than that of base metal, though all presenting poor HIC resistance for both base metal and the welded joint. The HIC cracks initiated from the inclusions enriching in Al, Ca, Si, Mn. The cracks are primarily transgranular, accompanying with limited intergranular ones.     

Comprehensive study on hydrogen induced cracking of electrical resistance welded API X52 pipeline steel

Different heat treatment cycles were designed in order to investigate the effect of microstructural changes on hydrogen induced cracking resistance (HIC) and mechanical properties of the electric resistance welded steel. The heat treating of the as-welded specimen improved the ductility and impact toughness. After heat treatment, the uniform hardness profile was obtained for the welded specimens. The removal of local hard zones reduced the risk of HIC. The chemical composition and clustering of inclusions have a deleterious effect on cracking resistance in the H₂S environment. Aluminosilicate compounds and MnS inclusions were favorite sites for HIC. The most promising post weld heat treatment for improving mechanical properties and the resistance to HIC was the application of two-cycle quenching followed by tempering.     

Investigation of the hydrogen induced cracking behaviour of API 5L X65 pipeline steel

The effect of microstructural features on the hydrogen induced cracking (HIC) susceptibility of two API 5L X65 pipeline steels were investigated by cathodic charging, hydrogen permeation and hydrogen microprint experiments. Microstructural evaluation after hydrogen charging revealed cracks at the mid-thickness (segregation zone) of both plates. However, more severe cracks were observed in the plate with higher dislocation density and residual stresses. The plate with lower plastic strain and more {111}-oriented grains had less severe cracks. Inclusions found along the crack path, comprising of Si-enriched oxides and carbides contributed to the initiation and propagation of cracks. The variation of the trapping behaviour and hydrogen diffusion through the plates were examined. The results confirmed that a higher ratio of reversible to irreversible traps contributes to increasing HIC severity in steels. Additionally, hydrogen transport through the steels was most prominent along the grain boundaries, indicating the importance of grain boundary character to HIC.     

An effective finite element model for the prediction of hydrogen induced cracking in steel pipelines

This paper presents a comprehensive finite element model for the numerical simulation of Hydrogen Induced Cracking (HIC) in steel pipelines exposed to sulphurous compounds, such as hydrogen sulphide (H₂S). The model is able to mimic the pressure build-up mechanism related to the recombination of atomic hydrogen into hydrogen gas within the crack cavity. In addition, the strong couplings between non-Fickian hydrogen diffusion, pressure build-up and crack extension are accounted for. In order to enhance the predictive capabilities of the proposed model, problem boundary conditions are based on actual in-field operating parameters, such as pH and partial pressure of H₂S. The computational results reported herein show that, during the extension phase, the propagating crack behaves like a trap attracting more hydrogen, and that the hydrostatic stress field at the crack tip speed-up HIC related crack initiation and growth. In addition, HIC is reduced when the pH increases and the partial pressure of H₂S decreases. Furthermore, the relation between the crack growth rate and (i) the initial crack radius and position, (ii) the pipe wall thickness and (iii) the fracture toughness, is also evaluated. Numerical results agree well with experimental data retrieved from the literature.     

Finite element model for charge and discharge cycle of activated carbon hydrogen storage

One of the main challenges to introduce hydrogen on the energy market is to improve on-board hydrogen storage and develop more efficient distribution technologies to increase the amount of stored gas while lowering the storage pressure. The physisorption of hydrogen on activated carbons (AC) is being investigated as a possible route for hydrogen storage. The objective of this work is to study the performance of adsorption-based hydrogen storage units from a "systems" point of view. A realistic two-dimensional axisymmetric geometric model which couples mass, momentum and energy balances is established based on the thermodynamic conservation laws using finite element method as implemented in COMSOL Multiphysics™. We consider the charging and discharging of the storage unit at a rated pressure of 9 MPa, and at an initial temperature of 302 K. The results are compared with experimental data obtained at the Hydrogen Research Institute of the University of Quebec at Trois-Rivieres. The storage tank is cooled by ice water. Research results show that both the simulated variations of pressure and temperatures during charge and discharge processes are in good agreement with the experimental data. The temperatures in the central region of tank are higher than those at the entrance and near the wall at the end of charge time while they are lower than those at the entrance and near the wall at the end of discharge time. The velocities are largest at the entrance, and decrease gradually along the axis of the tank. Owing to thermal effects, the larger flow rates result in less amount of adsorption in the condition of the same charging pressure. Hence measures of increasing heat transfer should be adopted, such as increasing the thermal conductivity of the storage bed. From the point of view of storage capacity, it is therefore possible to realize rapid hydrogenation, which is conducive to the use of such systems for on-board hydrogen storage based on activated carbon adsorption.     

Using a pulsed electric current to reduce the susceptibility of high-strength steel to the stress corrosion cracking caused by hydrogen

Under the tensile loading, the damage of metals in the corrosive medium is the most destructive and harmful. In this study, the stress corrosion cracking behavior of H-charged high-strength steel in 3.5 wt% NaCl solution after electropulsing treatment was investigated. The experimental results from elongation, yield strength, fracture morphology, and polarization curves all demonstrate the positive effect of the pulsed processing, as it reduced the susceptibility of steel to stress corrosion cracking by removing hydrogen by electropulsing. The reduction in hydrogen content of the pulsed high–strength steels was attributed to electromigration and increased system free energy, which drove the hydrogen atoms in the steel to de–trap and reduced the susceptibility to stress corrosion cracking.     

Mechanism of hydrogen-restrained crack propagation and practical application research of thermohydrogen treatment in a TiAl-based alloy

The mechanism of hydrogen-restrained crack propagation and practical application of thermohydrogen treatment in a TiAl-based alloy was investigated in this study. Hydrogenated and unhydrogenated alloys were subjected to high-temperature compression test, with a temperature range 1050–1200 °C and strain rate range 0.001–1 s⁻¹. The results showed that crack propagation was restrained due to hydrogen addition. The main mechanism of hydrogen-restrained crack propagation of such alloy was revealed that hydrogen-promoted lamella bending and hydrogen-decreased Young's modulus induced inter-lamellar cracks transforming into trans-lamellar cracks, decreasing cracks in the hydrogenated alloy. Additionally, hydrogen-induced mechanical twinning in γ-phase lamellae partly restrained inter-lamellar crack propagation. In the two-step forging process, the optimum forging parameters were determined. It was found that hydrogen could effectively restrain crack propagation during the two-step forging process. Hydrogen refined grains of the forged billets, which improved toughness of such billets. The hydrogen content of the forged hydrogenated billets could be decreased to a desired value, and the phase composition and content were basically identical to those of the initial unhydrogenated alloy.     

Dependence of hydrogen embrittlement on hydrogen in the surface layer in type 304 stainless steel

Hydrogen embrittlement (HE) together with the hydrogen transport behavior in hydrogen-charged type 304 stainless steel was investigated by combined tension and outgassing experiments. The hydrogen release rate and HE of hydrogen-charged 304 specimens increase with the hydrogen pressure for hydrogen-charging (or hydrogen content) and almost no HE is observed below the hydrogen content of 8.5 mass ppm. Baking at 433 K for 48 h can eliminate HE of the hydrogen-charged 304 specimen, while removing the surface layer will restore HE, which indicates that hydrogen in the surface layer plays the primary role in HE. Scanning electron microscopy (SEM) and scanning tunnel microscopy (STM) observations show that particles attributed to the strain-induced α′ martensite formation break away from the matrix and the small holes form during deformation on the specimen surface. With increasing strain, the connection among small holes along {111} slip planes of austenite will cause crack initiation on the surface, and then the hydrogen induced crack propagates from the surface to interior.     

Finite element simulation of welding residual stresses in martensitic steel pipes

Abstract

Finite element (FE) simulations of the welding of two high grade steel pipes are described. The first is a P91 steel pipe welded with a similar P91 weld consumable, and the second is a P92 steel pipe welded with dissimilar nickel–chromium based weld consumables. Both welds are multipass circumferential butt welds, having 73 weld beads in the P91 pipe and 36 beads in the P92 pipe. Since the pipes and welds are symmetric around their axes, the FE simulations are axisymmetric, allowing high FE mesh refinement and residual stress prediction accuracy. The FE simulations of the welding of the P91 and P92 pipes comprise thermal and sequentially coupled structural analyses. The thermal analyses model the heat evolution produced by the welding arc, determining the temperature history throughout the FE models. Structural analyses use the computed temperature history as input data to predict the residual stress fields throughout the models. Post-weld heat treatment (PWHT) of both pipes has also been numerically simulated by assuming that the FE models obey the Norton creep law during the hold time period at 760°C. The residual stresses presented here have all been validated by corresponding experimental measurements. Before PWHT, it has been found that, at certain locations in the weld region and heat affected zone (HAZ) in the pipes, tensile hoop and axial residual stresses approach the tensile strength of the material, presenting a high risk of failure. It has also been found that PWHT substantially reduces the magnitude of residual stresses by varying degrees depending on the material.     

Coupled computing for reactive hydrogen leakage phenomena with crack propagation in a pressure tank

This research focuses on reactive hydrogen leakage due to the fracturing of high-pressure hydrogen tanks to develop a coupled computing method that can simultaneously perform reactive fluid-structure interaction analyses. The integrated computational approach to reactive hydrogen leakage with combustion due to wall crack propagation was implemented using a hybrid of the coupled particle and Eulerian methods. This computational method provides valuable safety information for predicting crack propagation and hydrogen leakage with combustion reaction in pressure tanks as an essential part of assessing hydrogen as an energy vector. As a result, the effect of crack formation and wall boundary conditions on hydrogen turbulent diffusion and concentration distribution and the thermodynamic behavior of the chemical reaction were predicted computationally. It was found that a combustion reaction does not occur near the streamwise axis at the duct center because there is an excess of hydrogen fuel, which increases scalar dissipation, and that the combustion reaction separated into upper and lower regions as it proceeded.     

An integrated experimental and modeling approach to determine hydrogen diffusion and trapping in a high-strength steel

The hydrogen diffusion and trapping in AISI 4330 V high-strength steel is investigated by repetitive electrochemical hydrogen permeation tests, thermal desorption analysis, and hot and melt extraction. The analysis is coupled with a numerical model based on McNabb and Foster's kinetics with varying degrees of trap occupancy. The trapping parameters are obtained by fitting the numerical model to the experimental data, which permits to describe the diffusion and trapping processes for all tested conditions. In addition, the predictions calculated by the model are critically discussed and compared with those derived from usual approaches based on analytical solutions from Fick's laws and from Choo and Lee's method. An important difference is observed, indicating that the use of general analytical methods may not be adequate in the case of the studied steel. The use of more rigorous analysis provides a better understanding of trapping phenomenon and improved predictions of charging times and contents.     

Influence of grain boundary misorientation on hydrogen embrittlement in bi-crystal nickel

Computational techniques and tools have been developed to understand hydrogen embrittlement and hydrogen induced intergranular cracking based on grain boundary (GB) engineering with the help of computational materials engineering. This study can help to optimize GB misorientation configurations by identifying the cases that would improve the material properties increasing resistance to hydrogen embrittlement. In order to understand and optimize, it is important to understand the influence of misorientation angle on the atomic clustered hydrogen distribution under the impact of dilatational stress distributions. In this study, a number of bi-crystal models with tilt grain boundary (TGB) misorientation angles (θ) ranging between 0°≤ θ ≤ 90° were developed, with rotation performed about the [001] axis, using numerical microstructural finite element analysis. Subsequently, local stress and strain concentrations generated along the TGB (due to the difference in individual neighbouring crystals elastic anisotropy response as functions of misorientation angles) were evaluated when bi-crystals were subjected to overall uniform applied traction. Finally, the hydrogen distribution and segregations as a function of misorientation angles were studied. In real nickel, as opposed to the numerical model, geometrically necessary dislocations are generated due to GB misorientation. The generated dislocation motion along TGBs in response to dilatational mismatch varies depending on the misorientation angles. These generated dislocation motions affect the stress, strain and hydrogen distribution. Hydrogen segregates along these dislocations acting as traps and since the dislocation distribution varies depending on misorientation angles the hydrogen traps are also influenced by misorientation angles. From the results of numerical modelling it has been observed that the local stress, strain and hydrogen distributions are inhomogeneous, affected by the misorientation angles, orientations of neighbouring crystal and boundary conditions. In real material, as opposed to the numerical model, the clustered atomic hydrogens are segregated in traps near to the TGB due to the influence of dislocations developed under the effects of applied mechanical stress. The numerical model predicts maximum hydrogen concentrations are accumulated on the TGB with misorientation angles ranging between 15°< θ < 45°. This investigation reinforces the importance of GB engineering for designing and optimizing these materials to decrease hydrogen segregation arising from TGB misorientation angles.     

Finite element analysis of internal penny-shaped crack problem in an unbounded thermoelastic medium

In this article, the generalized thermoelastic theory with one relaxation time is used to investigate the internal penny-shaped crack problem in an unbounded thermoelastic material. The boundary of the crack is subjected to a prescribed stress distribution and temperature. The problem has been solved numerically using a finite element method. The results of the temperature, displacement components, and stress components have been verified numerically and are represented graphically. Comparisons are made with the results predicted by the presence and absence of the relaxation times.     

The synergistic effects of hydrogen embrittlement and transient gas flow conditions on integrity assessment of a precracked steel pipeline

We are reporting in this study the hydrogen permeation in the lattice structure of a steel pipeline designed for natural gas transportation by investigating the influence of blending gaseous hydrogen into natural gas flow and resulted internal pressure values on the structural integrity of cracked pipes. The presence of cracks may provoke pipeline failure and hydrogen leakage. The auto-ignition of hydrogen leaks, although been small, leads to a flame difficult to be seen. The latter makes such a phenomenon extremely dangerous as explosions became very likely to happen. In this paper, a reliable method is presented that can be used to predict the acceptable defect in order to reduce risks caused by pipe failure due to hydrogen embrittlement. The presented model takes into account the synergistic effects of transient gas flow conditions in pipelines and hydrogen embrittlement of steel material due to pressurized hydrogen gas permeation. It is found that blending hydrogen gas into natural gas pipelines increases the internal load on the pipeline walls due to overpressure values that may be reached in a transient gas flow regime. Also, the interaction between transient hydrogen gas flow and embrittlement of API 5L X52 steel pipeline was investigated using Failure Assessment Diagram (FAD) and the results have shown that transient flow enhances pipeline failure due to hydrogen permeation. It was shown that hydrogen embrittlement of steel pipelines in contact with the hydrogen environment, together with the transient gas flow and significantly increased transient pressure values, also increases the probability of failure of a cracked pipeline. Such a situation threatens the integrity of high stress pipelines, especially under the real working conditions of hydrogen gas transportation.