News and Views, Volume 54 | Benefits of High-Performance Computing for the Seismic Analysis of Critical Facilities

By Julio Garcia, PhD, PE

Abstract
Continued advances in the fields of structural and geotechnical engineering and increasing expectations for fidelity and refinement of seismic analysis models have resulted in a computational bottleneck emerging within the structural design activities for critical infrastructure. The Soil-Structure Interaction (SSI) analysis is key to characterize a realistic seismic response. Unfortunately, this is computationally demanding, given the need to characterize both the dynamic response of the structure and the site. The traditional SSI framework, System for Analysis of Soil Structure Interaction (SASSI), is a trusted analysis approach with decades of precedence. Our proprietary SSI analysis software, SC-SASSI, incorporates modern enhanced solution algorithms into SASSI and takes advantage of state-of-the-art computing, thereby drastically reducing analysis run-time.

The expanded SSI capabilities coded in SC-SASSI effectively eliminate the model size limitations of previous SASSI platforms and allow more accurate and streamlined treatment of behavior that have traditionally been neglected or often treated conservatively.

This article provides a variety of recent project examples highlighting how the enhanced solutions algorithms of SC-SASSI coupled with HPC capabilities have been leveraged, resulting in gains such as improved accuracy, reduced uncertainty, and accelerated schedules.

Introduction
Seismic analysis of critical industrial facilities requires the inclusion of the soil into the analysis model to evaluate the influence of the subgrade’s flexibility on the structure’s dynamic behavior. As a result, a soil-structure system is used for this purpose, and the consideration of the soil in the structural behavior is referred to as SSI. Industries other than nuclear also benefit from SSI consideration, both for seismic as well as for other dynamic loading conditions, like machinery vibration. For decades, the SASSI [1] program and its improved derivatives have been an accepted and favored analysis tool for the treatment of seismic SSI effects. With continued advances in the fields of structural and geotechnical engineering, and increasing detail and refinement of analysis models, SASSI has become a computational bottleneck in the design and evaluation of nuclear plant structures, and may not provide the fidelity and detail expected in modern engineering. Compared to the current state of practice for structural finite element (FE) analysis, gross model simplifications are too often perceived to be necessary to ensure reasonable analysis schedules, necessitating numerous assumptions, sensitivity studies, and additional qualitative evaluations to address specific aspects of the real seismic response.

SC-SASSI
Enhanced solution algorithms have become available since the development and widespread adoption of the original SASSI framework in the early 1980s. The commercial code, SC-SASSI [2], leverages these algorithms to efficiently take advantage of modern high-performance computing (HPC) with dramatic results. Figure 1 shows a flowchart highlighting some of the SC-SASSI improvements to the traditional SASSI. The gold boxes represent areas of key improvement and include the all-new RESPONSE module, which provides advanced postprocessing and is also a “living” module that continues to grow and adapt to the user’s needs.

Figure 1. Key SC-SASSI Improvements (Shown in Gold) to Traditional SASSI Modules.

Expanded SASSI capabilities via HPC functionality allow two primary advantages that lead to project gains: (1) large increases in permissible model sizes (often measured as a larger number of “interaction nodes”) and (2) drastically reduced analysis run-time. The large models now possibly allow direct treatment of SSI phenomena that previously required gross simplification and/or compromise. Key advantages of the HPC enhancements for large models are the handling of features like deep embedment, high-frequency hazard, structure-soil-structure interaction, and high-fidelity models for coupled stress analysis and seismic analysis. The example outlined below demonstrates some of the gains achieved using large models.

Figure 2. Example Nuclear Power Facility FE Model with Direct Consideration of SSSI.

Structure-Soil-Structure Interaction (SSSI)
An ideal example to demonstrate SSSI is that of a Nuclear power generating station. Typically, these sites have multiple structures located near one another, each with components critical to the facility’s safe operation.  The SSI behavior of one structure can affect the SSI behavior of an adjacent structure, referred to as structure-soil-structure interaction (SSSI).  Direct consideration of SSSI for multiple large structures has been impractical with traditional SASSI, necessitating the use of 2D “slice” SSSI models or “cascade analysis” techniques. These simplifications seek to reduce the size of the model(s) used to evaluate SSSI effects and/or avoid excessive computer run-times. However, they often insert significant inefficiencies (by way of required benchmarking analysis or decoupled, multi-step analyses) and thus remove realism.   

In contrast, SASSI with HPC provides efficiency to include multiple structures in the same SSI model, thus directly capturing SSSI effects.  In these direct SSSI models, entire multi-building nuclear island complexes have been incorporated into a single analysis run regardless of separation distance, soil properties, foundation depths, etc.  The benefits of direct treatment of SSSI effects include: (a) captures full 3D effects of SSSI behavior on adjacent structures and spatially variable amplitude of SSSI effects; (b) captures two-way feedback between adjacent structures; (c) minimizes the number of separate sensitivity study models to generate and maintain; and (d) allows more streamlined calculation of the relative displacements between adjacent structures.

An example multi-building FE model for direct consideration of SSSI effects in SC-SASSI is illustrated in Figure 2. This pressurized water reactor (PWR) complex [3] combines multiple buildings including the containment building, auxiliary building, radwaste building, turbine building, and intake structure into a single model.  In addition to capturing traditional SSSI effects, the combined model allowed consideration for partially shared load paths between adjacent structures, shared foundations despite separate super-structures, and relative movement between adjacent structures (to directly assess the potential for building impacts and failure of commodities spanning across seismic joints).

Benefits for Short Computational Times
HPC enhancements in SC-SASSI also allow much faster SSI analysis run-times than were possible with traditional SASSI.  The rapid analysis now allows more rigorous treatment of considerations that previously were handled via a sensitivity study or by the enveloping response of bounding cases. Key advantages of the HPC enhancements for shortening the computational times are the handling of features like probabilistic SSI analysis, multiple configurations, and/or alternate boundary conditions (like spent fuel installation facilities that may be analyzed under different cask loading configurations), as well as multiple hazard levels (i.e., design basis and beyond design basis earthquakes). The example outlined below demonstrates gains achieved using rapid analysis.

Probabilistic SSI Analysis
Probabilistic SSI analysis is advantageous over more traditional deterministic methods since it can directly address the inherent uncertainties related to the seismic input motion, soil and structural material behavior, and modeling assumptions. These uncertainties are directly addressed in the probabilistic SSI through randomization of the seismic input motions as well as soil and structural stiffness and damping. However, a drawback in probabilistic SSI analysis is the number of SSI models that are needed to obtain stable and reliable results, and the perceived computational effort associated with analysis of these various models. In practice, at least 30 SSI models are needed using the Latin Hypercube Sampling (LHS) method for the generation of SSI models. For this reason, many structures have traditionally been either analyzed using a deterministic approach or overly simplified for probabilistic analysis to achieve reasonable computational run-time. The use of HPC can mitigate this drawback by significantly reducing the analysis run-time, thus making it feasible to utilize detailed models for probabilistic analysis.

By covering realistic variability of key inputs like soil, structure, and ground motion, probabilistic SSI analyses provide a range of expected structural response results that allow easy identification of either realistic or conservative responses, according to the objective of the evaluation; this represents a significant advantage over deterministic analysis.

The case described below provides an example where probabilistic SSI analyses were used to obtain realistic responses for the purpose of fragility analysis of structures, systems, and components (SSCs) in a nuclear power plant [3]. Each of the 30 probabilistic SSI cases, which included variation of the soil, structural stiffness, damping, and earthquake ground motion, was run for each model. All 30 SSI cases were run in less than one day for each model. These analyses generated 30 In-Structure Response Spectra (ISRS) for each of the selected locations, and the median and 80th percentile spectra were directly extracted. The response spectra variability (βRS), which often must be assumed in fragility calculations, could be directly and accurately calculated for each set of ISRS. Results post-processing is efficiently automated in SC-SASSI. Figure 3 illustrates each resulting ISRS and their post-processing for a typical location.

Figure 3. Example Calculation of Median and 80th% ISRS from Suite of Probabilistic Results.

Ground Motion Incoherence
SC-SASSI offers additional useful features, like ground motion incoherence effects. Ground motion incoherence refers to the spatial variation of the ground motion. In other words, there are differences in the ground motion experienced by two points separated by some distance; the larger the distance, the larger the variability. Considering the effects of ground motion incoherency in SSI analysis enhances the realism of predicted structural response and generally reduces high-frequency response. The relative significance of this refinement depends on several factors, including the foundation input motion’s high-frequency energy content and the foundation footprint’s size. The inclusion of ground motion incoherence in SASSI requires additional computational effort versus coherent ground motion. A common approach in the industry to minimize the additional computational effort is to limit the number of coherency modes in the analysis to some minimum number of modes that can still capture realistic SSI behavior. However, the industry debates on how many modes are sufficient under different conditions. HPC capabilities in SC-SASSI mitigate the computational penalty for maintaining many coherency modes, thus allowing consideration of many modes, even for large and complex structures.

Figure 4. Synchronous Condenser on Pile Foundations and Retrofit Evaluation.

Industrial Facilities on Pile Foundations
Pile foundations are frequently used in industrial facilities. While their static analysis and design are widely understood and applied in practice, the contribution to the dynamic response of piles is often neglected, leading to a mischaracterization of the response, like in the case portrayed below.

At the subject hydroelectric generation facility, the synchronous condenser, responsible for adjusting the electric power before it gets input into the transmission grid, experienced undesirable vibration levels during its initial stages of operation. The source of the high vibration levels is attributed to the resonant dynamic interaction between the condenser, its foundation, and the supporting soil media, which went undetected during the initial design process. The foundation of the condenser consists of a massive reinforced concrete block supported on steel casing piles with reinforced concrete infill.

A detailed three-dimensional finite element SSI model was developed to accurately represent the interaction between forces, structure, soil, and foundation, as displayed in Figure 4a. The model was analyzed using the computer software SC-SASSI, including the recently implemented and validated pile elements. The SSI analysis was able to diagnose the problem, and it also served to evaluate feasible remediation solutions (see Figure 4b). The selection of the retrofit solution was performed with the objectives of (1) reducing vibration amplitudes, (2) shifting the system frequency away from the operation frequency, both (1) and (2) being achieved by increasing the system stiffness and therefore minimizing potential resonance conditions, and (3) minimizing the disruption of the surrounding critical equipment already installed.

Conclusions
Since SASSI’s development and widespread adoption, enhanced solution algorithms have become available. SC-SASSI has incorporated these into the trusted SASSI framework to efficiently take advantage of modern high-performance computing (HPC). Leveraging SC-SASSI to fully utilize HPC can lead to dramatic project gains with reduced schedules, improved realism, and often more targeted, less conservative solutions.

Expanded SASSI capabilities via HPC functionality allow two primary advantages that lead to project gains: (1) large increases in permissible model sizes and (2) drastically reduced analysis run-time.  These two primary advantages allow direct and rigorous treatment of SSI phenomena that previously required gross simplification, compromise, sensitivity studies, and/or overly conservative assumptions.  More specifically, SASSI with HPC allows project improvements in areas such as (a) deep embedment with soft soil and high-frequency hazard, (b) structure-soil-structure interaction, (c) combined high-fidelity models for coupled stress analysis and seismic analysis, (d) probabilistic SSI analysis, (e) multiple configurations and/or alternate boundary conditions, (f) multiple hazard levels, (g) pile foundations, and (h) ground motion incoherency effects.

Project examples illustrated in this article demonstrate that the HPC capabilities in SC-SASSI allow treatment of complex SSI phenomena more rigorously and realistically than traditional SASSI approaches, while at the same time offering faster and more efficient analyses. When implemented to their maximum potential, HPC enhancements in SC-SASSI can lead to project gains by reducing cost, schedule, and risk.

References

  1. Lysmer, J., Tabatabaie, M., Tajirian, F.,  Vahdani, S., and Ostadan, F.  “SASSI – A system for analysis of soil-structure interaction”. Technical Report UCB/GT/81-02, Geotechnical Engineering, University of California, Berkeley, 1981.
  2. SC Solutions, Inc. “SC-SASSI Manual, Version 2.3.2”. Sunnyvale, California, 2018.
  3. Kosbab, B., Garcia, J., Li, W.,  Tran, H., Talebinejad, I., Bassam, A., and Tilow, K. “Soil-Structure Interaction: Project Gains through use of High-Performance Computing”. Transactions, SMiRT-24, Busan, Korea, 2017.

Get News & Views, Volume 54

News and Views, Volume 54 | Advanced NDE for Hydroelectric Penstock Inspection

By:  Jason Van Velsor and Jeff Milligan

At SI, we regularly combine advanced NDE inspections with fitness for service evaluations to provide value-added solutions for our clients.

Hydroelectric power plants harness two of the most powerful forces on earth, water and gravity.  The integrity of the penstocks that flow water to and away from the turbines in these plants is paramount to safe operation and the safety of the surrounding population. With many hydro plants approaching 100 years of service, critical issues can arise with these penstocks, which may have little to no fabrication documentation, may have significant fabrication imperfections, and may have significant accumulated damage from the environment and many years of service.

Our talented, highly experienced NDE staff can complete in-depth hydroelectric penstock inspections for our clients, providing peace of mind that their assets are safe to continue operating long into the future.

Example of Inspections Completed By SI

Figure 1. (Top) C-Scan image of corrosion spots on the ID of a penstock with the color scheme relative to material thickness; (Bottom) B-Scan image showing the thinnest thickness reading, 0.140 inch.

Phased Array Ultrasonic (PAUT) Corrosion Mapping of Penstocks

SI works with hydroelectric utility companies to provide detailed thickness measurements on specific areas of penstocks. Phased array ultrasonic corrosion mapping provides hundreds of thousands of thickness readings that can produce a detailed image of the inside surface of a penstock, which may have experienced corrosion or erosion that limits penstock life. This highly detailed scan produces the data required for accurate Fitness for Service (FFS) calculations. PAUT corrosion mapping can be performed from the outside diameter of penstocks that are above ground to detect internal wall thinning, or it can be performed from the inside of a penstock, as long as safe access and confined space requirements are met.  An example of the phased array ultrasonic data can be seen in Figure 1 Thickness measurement values, often taken every 0.04 inch (1 mm) along the scan area can be saved and exported using industry standard formats (e.g., CSV, excel, etc.) to support further statistical analysis of the ultrasonic data.

Corrosion Surface Profilometry

Figure 2. Laser scan image of corrosion on the OD of a penstock with pit depths identified.

For a pitted or corroded surface that is accessible to an inspector (ID or OD), the use of a laser profilometry device can be a valuable tool to map and depth size corrosion.  Conventional pit gauge measurement of corrosion on penstock surfaces can be a time consuming and inaccurate process.  The results rely heavily on the experience of the technician as well as the surrounding surface condition.  A more efficient and accurate method is to create an exact image of the surface using a handheld laser scanner.  The laser scanning process is very fast and the results can be displayed as a 3D surface or unrolled to a 2D view.  Reconstruction of the surface is real-time, with color coding used to provide a visual relevance for material loss.  Off-line analysis can be used to make discrete readings of wall loss, or the full map data file can be exported using standard file formats to allow other subsequent analysis to be conducted.

Short Range Guided Wave Technique to Inspect Riveted Joints for Crevice Corrosion

Figure 3. Drawing of potential location for crevice corrosion on a riveted lap-joint.

Figure 4. Mock-up drawing and actual examination data from SR-GWT test.

Some penstocks with riveted joints may suffer from crevice corrosion due to the construction geometry of the riveted lap joints.  The lap-joint area is generally covered in concrete and buried so the external surface is not easily accessible. The concrete can disbond over time and water collects and runs along the crevice of the plate and butt strap, causing corrosion. The area of crevice corrosion may not be assessed from the internal surface with traditional ultrasonic methods due to the obstruction from the internal butt strap, therefore SI has developed a short range guided wave testing (SR-GWT) technique to inspect for this crevice corrosion. Figure 3 shows an image to help visualize the issue.  For this technique electromagnetic acoustic transducers (EMAT) are utilized to propagate sound waves along the volume of the penstock plate to detect a change in the cross-sectional area.  Figure 4 shows an example of inspection data that was collected on a calibration plate to prove this technique.

Lap-Welded Longitudinal Seam Inspection

Our staff can complete lap-welded seam identification and inspections. Lap-welded or forge-weld longitudinal seam pipes and penstocks were manufactured in the 1920s (Figure 5). Oftentimes, little to no fabrication documentation exists that will tell a hydroelectric utility if their penstock cans were made with lap-welded seams. SI developed a phased array ultrasonic examination technique to identify areas of lap-welded seams and look for lack of fusion and service damage.

SI developed two different PAUT techniques to identify and examine lap-welded penstocks. A refracted shear wave technique is preferred for ID-connected indications, while a longitudinal wave technique produces sound that is more perpendicular to the weld bond line of the lap-welded joint.  Figure 6 shows PAUT scan data along with an explanation of typical features from three different lap-welded seam examinations.

ID/OD Girth Weld Inspection

Figure 5. Illustrated cross-section of a lap weld.

Phased array ultrasonic examination of girth welds is also a common inspection for hydro penstocks. SI’s vast experience inspecting high energy piping systems and pressure vessels translate perfectly to penstock girth weld applications. SI uses ultrasonic software simulation programs to create scan plans that calculate the necessary beam angles and focusing to ensure 100% weld coverage during a PAUT examination. Encoding PAUT data with an automated, semi-automated, or manual encoding device allows for off-line analysis and a permanent record of inspection data. SI also has the ability to create custom scanners, probes, and wedges, when required for difficult inspection applications.  The advanced NDE equipment and experience of SI is unmatched. 

Figure 6. PAUT scans at 0° incidence (image is analogous to a cross-section of the weld) from hammer-welded pipe showing seam indications (labelled “1”) increased wall thickness at the seam (labelled “2”), and relatively clean areas adjacent to the seam (labelled “3”). The top image and center images show indications from the bond line (likely indicating lack-of-fusion), while the bottom image shows no indications from the bond line. Note that the base material on the left and right ends of the image shows dense areas of small indications, likely attributable to numerous impurities and inclusions in low-quality skelp from which this pipe was manufactured.

SI Can Help

The situations discussed previously demonstrate solutions that SI has developed for hydroelectric penstocks. With our extensive expertise, SI can work with our clients to develop and execute inspection solutions that are customized to their specific needs and provide various engineering analyses based on the examination results.

 

Get News & Views, Volume 54

News & Views, Volume 54

The next issue of Structural Integrity’s bi-annual technical magazine, News and Views, is now available. In this issue, experts from our Nuclear, Energy Services, Nondestructive Examination (NDE), Pipeline Integrity Compliance Solutions, and Metallurgical Laboratory groups share project case studies and details on technical advancements across a range of industries and disciplines.

Get News & Views, Volume 54

News and Views, Volume 54 | Materials Lab Featured Damage Mechanism

STRAIN-INDUCED PRECIPITATION HARDENING (SIPH) IN AUSTENITIC STAINLESS BOILER TUBES

By:  Wendy Weiss

Structural Integrity’s Metallurgical Laboratory offers comprehensive metallurgical laboratory services to support client material issues.

Strain-Induced Precipitation Hardening, also known as SIPH, is a commonly misinterpreted boiler tube failure mechanism that occurs when austenitic stainless steel tubing is cold or warm worked during fabrication and then is installed with either improper or no solution annealing heat treatment. While the basic mechanism and the root cause are understood, the complex interaction between heat chemistry, quantity of cold or warm work, and subsequent thermal history makes it very difficult to predict under precisely what circumstances damage due to SIPH will result in failure of a boiler tube.

Longitudinally oriented crack at the extrados of a bend in a stainless steel superheater tube.

Mechanism
SIPH occurs when a heat of austenitic stainless steel containing certain precipitate-forming elements (e.g., niobium, titanium, vanadium, etc.) either intentionally or as residuals is cold or warm worked during subsequent material processing. The cold or warm working creates excess defects in the sub-structure of the material, which serve as preferred sites for precipitation of temper-resistant carbides or carbo-nitrides. Precipitation occurs when the material is heated to a sufficiently high temperature that is well below the solution annealing temperature. This can occur rapidly during a poorly executed heat treatment if the material does not reach the proper solution annealing temperature, or it can occur more slowly at typical operating temperatures for superheater or reheater tubing in utility-type boilers.

Once formed, precipitates anchor to the defects, resulting in a substantial increase in the elevated temperature or creep strength of the interior of the grains. At the same time, there is a narrow zone of material immediately adjoining the grain boundaries that remains largely precipitate-free due to the diffusional characteristics of the grain boundary itself. Ultimately, the interior portion of the individual grains becomes very strong at elevated temperatures while the material immediately adjoining the grain boundaries becomes comparatively creep weak. In addition, any surface-active elements that may be present in the material, such as arsenic, tin, antimony, etc., will tend to concentrate at the grain boundaries, further reducing their strength.

Regardless of when the precipitation occurs, once the interior of the grains has been strengthened, the grain boundary regions are weakened. Any strain imposed on the material in response to an applied or residual stress is forced to concentrate in the grain boundary region, which substantially magnifies its effect. For example, suppose the bulk strain experienced by a cold-worked stainless superheater tube segment is very small – a fraction of a percent – and the material has undergone the SIPH reaction. The strengthened grain interiors will undergo no strain. Conversely, within the much smaller volume of the comparatively weak grain boundaries, the accumulated strain will be orders of magnitude higher than the bulk level.

Typical Locations

  • Bends
  • Offsets
  • Swages
  • Welded attachments

Features

  • Can initiate at midwall
  • Early-stage damage consists largely of grain boundary cavities and microfissures
  • Intergranular cracking

Root Causes
The single root cause of SIPH is the failure to properly solution anneal susceptible heats of austenitic stainless steel tubing that has been either cold or warm worked during fabrication.

An overall, cross-sectional view of the intergranular crack (image A) with higher magnification views of some of the secondary grain boundary microfissures and voids (image C). Slip bands (parallel lines in image B) indicate local deformation.

For any situation involving material property characterization, Structural Integrity has an experienced group of materials specialists and a full-service metallurgical testing laboratory that can help.

Visit Our Metallurgical Laboratory Services Page

SEE THESE OTHER RELATED MATERIALS FROM THE SIA METALLURGICAL LABORATORY

Featured Articles


Case Studies


Metallurgical Laboratory Product & Services Information


News & Views, Volume 53 | Serviceability Assessment of an L-Grade Stainless Steel Pipe Fitting

By: Terry Totemeier

A client recently ordered a Type 316 stainless steel pipe coupling fitting for use in a high-pressure, high-temperature steam line operating at 1005°F.  The fitting that was received was so-called dual grade Type 316/316L stainless steel.  Given the limitations on using “L” grades of stainless steel at high temperatures, the client requested that SI perform a serviceability assessment for the fitting to determine if it could be safely used until the next scheduled outage when a replacement non-L grade fitting would be available.

BACKGROUND
The fitting ordered was a ½” nominal diameter (NPS ½), 6000# (Class 6000) full coupling socket-welding fitting in accordance with the ASME B16.11 specification, material ASME SA-182 forging, Type 316 stainless steel (designated as F316 in SA-182).  The fitting supplied was dual grade F316/316L material with a carbon content of 0.023% per the material test certificate.  The designation of this material as “dual grade” means that it meets the requirements of both F316 and F316L material grades.  This is possible because the chemical composition requirements of these two grades overlap, with the primary difference between them being carbon content.  For F316 the carbon content is specified to be 0.08% maximum (no minimum), while for F316L the carbon content is specified to be 0.030% maximum.  Therefore, material with carbon content less than 0.030% will meet the requirements for both grades.  It is worth noting that the carbon content of “H” grade of 316 stainless steel (F316H per SA-182) is specified to be 0.04-0.10%.  The H grade is intended for use at high temperatures.

The received fitting was installed in a main steam valve pressure equalizing line with a steam temperature/pressure of 2750 psia/1015°F at design conditions and 2520 psia/1005°F at operating conditions.  The fitting was welded to Grade P11 pipe on one side and Grade P22 pipe on the other side.  The applicable code was stated to be ASME BPVC Section I.

With a reported carbon content of less than 0.04%, the fitting is technically not permitted for use in ASME Section I construction above a temperature of 1000°F.  Per the ASME Boiler and Pressure Vessel Code (BPVC) Section II, Part D, Table 1A, the allowable stresses for SA-182, F316 material are valid at or above 1000°F only when the carbon content is greater than 0.04% (Note G12).  Per the same table, SA-182, F316L material is only permitted for use in Section I construction up to 850°F.  The reason for this temperature limitation is that the long-term creep-rupture strength of Type 316 stainless steel with lower carbon content is reduced compared to material with higher carbon content because fewer carbides form during service to strengthen the grain boundaries.  There are no other adverse impacts of the lower carbon content, e.g., on fatigue strength or oxidation resistance.

The short-term serviceability of the fitting with low carbon content was assessed by comparing bounding pressure stresses in the fitting with the reported creep-rupture strength for Type 316L material.  Per the ASME B16.11 specification, Class 6000 socket-welding fittings are compatible with NPS Schedule 160 pipe, meaning that pressure stresses in the fitting will be less than those in Sch 160 pipe with minimum wall thickness according to ASME B36.10 (pipe dimension specification), in other words, the fitting will be at least as strong as the pipe.  

ASSESSMENT
The dimensions of NPS ½, Schedule 160 pipe per the ASME B36.10 pipe specification are 0.84” outer diameter (OD), 0.165” minimum wall thickness (MWT).  For an operating steam pressure of 2,520 psi, the reference hoop stress per the equation in ASME BPVC Section I, Appendix A-317 is 5.05 ksi.  Per the general design guidance in ASME B16.11 (Section 2.1.1) the pressure stresses in the fitting must be less than this.  

Figure 1. Schematic diagram for a socket-welding coupling fitting. Per ASME B16.11, an NPS ½, Class 6000 fitting has relevant dimensions B = 0.875” maximum, C = 0.204” minimum, and D = 0.434” minimum.

Since the fitting in question is cylindrical, comparative hoop stresses can also be calculated from dimensions given in ASME B16.11, although these may not be exact due to the varied wall thickness in the fitting.  According to Table I-1 of ASME B16.11, the central body of the fitting is 1.283” OD and 0.395” MWT (Figure 1).  The reference hoop stress calculated using the A-317 equation at 2,520 psi stream pressure and these dimensions is 2.63 ksi, considerably less than 5.05 ksi.  In the female socket ends of the fitting, the OD is also 1.283”, but the minimum wall thickness is 0.204”, leading to a calculated reference hoop pressure stress of 6.58 ksi.  Note that the actual stresses in the socket ends will be much less than this because the pipe will be inserted and welded into the socket, taking up the pressure loading, but the calculated stress can be taken as a bounding value.

Creep-rupture strengths for Type 316L stainless steel have been reported in ASTM Data Series DS 5S2 publication, “An Evaluation of the Yield, Tensile, Creep, and Rupture Strengths of Wrought 304, 316, 321, and 347 Stainless Steels at Elevated Temperatures” (ASTM, 1969).  According to Table 7 in this report, the average 10,000 hour creep-rupture strengths for Type 316L at 1000°F and 1050°F are 34.5 and 25 ksi, respectively.  Minimum creep-rupture strengths are typically taken as 80% of the average strength, so the inferred minimum strengths at 1000°F and 1050°F are 27.6 and 20 ksi, respectively.  

The reported 10,000 hour creep-rupture strengths in the temperature range of interest are more than twice the calculated bounding pressure stresses in the fitting, so it was judged that there is very little risk of failure of the fitting by creep-rupture in the next 10,000 hours of service.

This result is unsurprising since the 1005°F is barely into the creep range for Type 316 regardless of carbon content.  The carbon content effects become more pronounced at higher temperatures (approximately 1100°F and above).

CONCLUSION
Based on the above assessment, it was SI’s opinion that the Type 316L fitting with carbon content less than 0.03% was suitable for a limited period of service (less than 10,000 hours) until it can be replaced.  Given that the fitting is reportedly welded to low-alloy steel pipe on either side, SI also recommended that a Grade 22 (2.25Cr-1Mo) low-alloy steel fitting be considered as a replacement, which would eliminate dissimilar metal welds (DMWs) between the fitting and pipes.  DMWs are prone to premature failure due to thermal fatigue, weld fusion line cracking, and decarburization of the ferritic material. This voluntary recommendation made by SI, was not part of the original scope of work, but may have been just as critical a finding as it shed light upon a failure risk previously unknown by the client. 

Get News & Views, Volume 53

News & Views, Volume 53 | Encoded Phased Array Ultrasonic Examination Services for Cast Austenitic Stainless Steel (CASS) Piping Welds

IN PRESSURIZED WATER REACTOR (PWR) COOLANT SYSTEMS

By:  John Hayden and Jason Van Velsor

The CASS piping welds present in many PWR plants provide numerous and complicated challenges to their effective ultrasonic examinations. To this point, a viable ultrasonic examination solution for the inspection of these piping components, as required by ASME Code Section IX,  had previously not been available. By leveraging our technical expertise in materials, technology development, and advanced NDE deployment, Structural Integrity Associates, Inc (SI) has developed a new system that will provide a meaningful solution for the examination of CASS piping components. The result of this program will be the first commercial offering for the volumetric examination of CASS components in the nuclear industry.

BACKGROUND INFORMATION
ASME Section XI Class 1 RCS piping system welds fabricated using CASS materials pose serious and well-understood challenges to their effective ultrasonic examination. For decades, utilities and regulators have struggled with the administrative and financial burdens of Relief Requests, which were, and still are, based on the inability to perform meaningful volumetric examinations of welds in CASS components. 

Many years of futility and frustration may have fostered the belief that technology allowing effective and meaningful examination of CASS materials would never be achievable. This is no longer the case.

The failure mechanism for CASS material occurs through the loss of fracture toughness due to thermal aging embrittlement. The susceptibility of CASS material to thermal aging embrittlement is strongly affected by several factors, primary of which are system operating time and temperature, the casting method used during component manufacture, and molybdenum and ferrite content. In addition to the existing ASME Section XI requirements for the examination of welds in CASS materials, the susceptibility to thermal aging embrittlement drives the requirement for additional examinations (including ultrasonic examinations) as directed by several NRC-published NUREGs required for plant license renewal. The existence of a viable, effective examination capability for CASS materials plays a very important part in both currently required Inservice Inspections (ISI) and plant license renewal.

Figure 1. An example of the widely-varying microstructure of a centrifugally cast piping segment. False-color imaging is used to aid visualizing grain variations. (Image from NUREG/CR-6933 PNNL-16292)

CASS MATERIAL PROPERTIES AND EFFECT ON ULTRASONIC EXAMINATION
Metallurgical studies have revealed that the microstructure of CASS piping can vary drastically in the radial (through-wall) direction, as well as around the circumference and along the length of any given piping segment. Large and small equiaxed, columnar and mixed (combinations of equiaxed and columnar grains), and banding (layers of substantially different grain structures) are commonly observed in CASS piping materials. None of these conditions favor the performance of effective ultrasonic examinations.

Figure 2. PWR RCS Major Components

The very large and widely varying types (equiaxed, columnar, and randomly mixed), sizes and orientations of the anisotropic grains in CASS material are very problematic. Anisotropic is defined as an object or substance having a physical property that has a different value when measured in different directions. Such physical properties strongly affect the propagation of ultrasound in CASS material by causing severe attenuation (loss of energy through beam scattering and absorption), beam redirection, and unpredictable changes in ultrasonic wave velocity. These factors are responsible for the inability of ultrasonic examination to completely and reliably interrogate the Code-required volume (inner 1/3 Tnom) of welds in CASS piping material. Interestingly, CASS materials less than 1.6” Tnom (Pressurizer Surge Piping) can be effectively examined, while CASS materials over 2.00” (Main RCS Coolant Loop Piping) are less effectively examined.  Consequently, an ASME Section XI, Appendix VIII qualification program for CASS piping components has not been established and remains in the course of preparation. Nonetheless, ASME Section XI requirements to conduct inservice examinations of RCS piping welds fabricated from CASS components remain fully in force.

ASME CODE ACTIONS AFFECTING CASS PIPING EXAMINATIONS
ASME Section XI Code Case N-824, “Ultrasonic Examination of Cast Austenitic Piping Welds From the Outside Surface,” was approved by ASME in October 2012 and by the NRC in October 2019. This Code Case provides the first approved direction for the ultrasonic examination of welds joining CASS piping components. The ASME B&PV Code, Section XI, 2015 Edition, incorporates Code Case N 824 into Mandatory Appendix III in the form of Mandatory Supplement 2. To date, these two ASME Section XI Code documents remain the sole sources approved by ASME and NRC that provide specific direction for the examination of CASS RCS piping system welds and, therefore, form the foundation of SI’s approach for the development of our CASS ultrasonic examination solution.

SI’S CASS PROGRAM DESCRIPTION
SI is developing the industry’s most well-conceived and capable ultrasonic system for the examination of welds in CASS piping components. To accomplish this objective, SI has drawn upon our internal knowledge and experience, supplemented by a careful study of numerous authoritative bodies of knowledge relating to the examination of CASS components. The development of the SI examination system has been guided by both SI’s industry-leading 17 years of experience conducting phased array examinations in nuclear power plants and the knowledge acquired through the careful study of the topical information contained within industry-recognized publications. These published results of extensive industry research provided both guidance for the selection of phased array system components and CASS-specific material insights that strengthen the technical content of our Appendix III-based procedure. 

Figure 3. RCS Coolant Pump and Crossover Piping

CASS PROGRAM ELEMENTS
SI believes that the procedure, equipment and personnel featured in this program will be equivalent or superior to those that will form the industry-consensus approach for CASS ultrasonic examinations needed to successfully achieve Appendix VIII, (future) Supplement 9, “Qualification Requirements for Cast Austenitic Piping Welds.”

Ultrasonic Procedure – SI has crafted an ultrasonic examination procedure framework that is fully compliant with ASME Section XI, Mandatory Appendix III, Supplement 2, along with referenced Section XI Appendices as modified by the applicable regulatory documents.

Ultrasonic Equipment – SI has acquired and assembled the ultrasonic system components required by Code Case N-824 and Appendix III, Supplement 2, which includes the following:

  • Ultrasonic instrumentation capable of functioning over the entire expected range of examination frequencies. The standard examination frequency range extends from low-frequency, 500 KHz operation for RCS main loop piping welds through 1.0 MHz for pressurizer surge piping. 

SI has designed and acquired additional phased array transducers that meet the physical requirements of frequency, wave mode, and aperture size and are capable of generating the prescribed examination angles with the required focal properties. SI has designed and fabricated an assortment of wedge assemblies that will be mated with our phased array probes to provide effective sound field coupling to the CASS components being examined. SI’s wedge designs consider the CASS pipe outside diameter and thickness dimensions and employ natural wedge-to-material refraction to assure optimal energy transmission and sound field focusing.

SI also possesses several data encoding options that are necessary to acquire ultrasonic data over the expected range of component access and surface conditions. The encoding options will include:

  • Fully-automated scanning system, capable of driving the relatively large and heavy 500KHz phased array probes
  • The SI-developed Latitude manually-driven encoding system, which has been deployed during PDI-qualified dissimilar metal DM weld examinations in nuclear power plants

    Figure 4. Steam Generator Details

Examination Personnel – SI’s ultrasonic examination personnel are thoroughly trained and experienced in all elements of encoded ultrasonic data acquisition and analysis in nuclear plants. SI’s examiners have a minimum of 10 years of experience and hold multiple PDI qualifications in manual and encoded techniques. SI recognizes the challenges that exist with the examination of CASS piping welds and has developed a comprehensive program of specialized, mandatory training for personnel involved with CASS examinations. This training includes descriptions of coarse grain structures, their effect on the ultrasonic beam, and the expected ultrasonic response characteristics of metallurgical and flaw reflectors, as well as the evaluation of CASS component surface conditions.

ULTRASONIC TECHNQUE VALIDATION
Although not required by the ASME Code, SI has arranged for access to CASS piping system specimens from reputable sources to validate the efficiency of our data acquisition process and the performance of our ultrasonic examination techniques. The specimens represent various pipe sizes and wall thicknesses and contain flaws of known location and size to permit the validation and optimization of SI’s data acquisition and analysis processes. SI will thoroughly analyze, document, and publish the results of our system performance during the examination of the subject CASS specimens.

Figure 5. Pressurizer and Surge Line Details

CASS PIPING SYSTEM APPLICATIONS
Typical CASS Piping Weld Locations in PWR Reactor Coolant Systems
The following graphic illustrates the location and extent of CASS materials in the RCS of many PWR plants.

RCS Main Loop Piping Welds: This portion of the RCS contains large diameter butt welds that join centrifugally cast stainless steel (CCSS) piping segments to statically cast stainless steel (SCSS) elbows and reactor coolant pump (RCP) casings. RCS main loop piping includes the following subassemblies:

  • Hot leg piping from the Reactor Vessel Outlet to the SG Inlet
  • Cross-over piping from the SG Outlet to the RCP Inlet
  • Cold leg piping from the RCP Outlet to the RPV Inlet

Steam Generator Inlet / Outlet Nozzle DM Welds: These terminal end DM butt welds are present in PWR plants, both with and without safe ends between the SCSS elbows and the ferritic steel nozzle forgings. 

Pressurizer Surge Piping Welds: This portion of the RCS contains a series of butt welds fabricated using CCSS piping segments to SCSS elbows between the Pressurizer Surge nozzle end and the Hot Leg Surge nozzle. 

SUMMARY
The CASS piping welds present in many PWR plants provide numerous and complicated challenges to their effective ultrasonic examinations. SI’s new CASS ultrasonic examination system will provide a new and meaningful solution.

PROJECT TIMELINE
SI is working to complete the development, integration and capability demonstrations of the CASS ultrasonic examination system described in this document for limited (emergent) fall 2023 and scheduled deployments beginning in spring 2024.

Get News & Views, Volume 53

News & Views, Volume 53 | Phased Array Ultrasonic Testing (PAUT) Monitoring with Ultrasonic Thick-Film Arrays

Traditional nondestructive examination (NDE) activities are planned based on hours of service, number of load cycles, time elapsed since previous inspections, or after the emergence of clear and obvious damage in a component. While engineering judgment and risk analysis can, and should, be used to prioritize inspections, these prioritizations are not based on the actual physical condition of the component or material it is constructed from but on precursory conditions that may or may not lead to eventual damage. Alternatively, continuous monitoring approaches can facilitate advanced planning and the optimization of Operations and Maintenance (O&M) spending by enabling the prioritization of inspections based on a component’s actual current condition. Furthermore, continuous monitoring enables earlier detection, which allows the extension of the component’s remaining useful life through modified operation. 

SI’s recent advances with thick-film are breakthrough technologies for long-term monitoring and imaging of crack growth in critical components.

Figure 1. Installed thick-film UT sensors for thickness monitoring of elbows.

Given the trend of fewer on-site resources and tighter O&M budgets, the energy industry has a strong motivation to progress toward condition-based inspection and maintenance. To facilitate this evolution in asset management strategy, new monitoring sensor technologies are needed, ones that provide meaningful monitoring data directly correlated to the condition of the material or asset. To support this need Structural Integrity has developed a novel thick-film ultrasonic sensor solution. Initially developed for basic applications, such as thickness monitoring, SI’s recent advances with this technology make long-term monitoring and imaging of crack growth in critical components possible.

BACKGROUND
Ultrasonic thick-films are comprised of a piezoelectric ceramic coating that is deposited on the surface of the component that will be monitored. A conductive layer is then placed over the ceramic layer, and the ceramic layer deforms when an electric potential is applied across the film. When a sinusoidal excitation pulse in the ultrasonic frequency range is applied across the film, the vibration of the film is transferred into the test component as an ultrasonic stress wave.   

Structural Integrity initially developed our thick-film ultrasonic sensors for real-time thickness monitoring and has demonstrated the performance and longevity of this technology through laboratory testing and installation in industrial power plant environments, as seen in the photograph in Figure 1, where the sensors have been installed on multiple high-temperature piping components that are susceptible to wall thinning from erosion. In this application, the sensors are fabricated directly on the pipe’s external surface, covered with a protective coating, and then covered with the original piping insulation. Following installation, data can either be collected and transferred automatically using an installed data acquisition instrument, or a connection panel can be installed that permits users to acquire data periodically using a traditional off-the-shelf ultrasonic instrument. Example ultrasonic datasets are shown in Figure 2.

Figure 2. Ultrasonic datasets from an installed thick-film UT sensor at two different points in time.

TECHNOLOGY ADVANCEMENTS
Recently, SI has demonstrated the ability to create thick-film sensors with complex element arrays that can be individually controlled to steer and focus the sound field, as with traditional phased-array ultrasonic testing (PAUT). Moreover, data from individual array elements can be acquired and post-processed using full-matrix capture (FMC) techniques. FMC is a data acquisition technique where all elements in the array are used to both transmit and receive ultrasonic waves. The result is a large data matrix that can be used for further processing with various post-processing techniques. Compared to more traditional active focusing, FMC is well-suited for a fixed transducer array, as scanning speed is not a concern. Another advantage is that the electronics needed for data acquisition can be simplified – requiring only a single pulsing channel.

A thick-film Linear-Phased Array (LPA) installed on a standard calibration block is shown in Figure 3. The two images shown on the right were generated using the Total Focusing Method (TFM) post-processing algorithm, with the image on the far right having an adjusted color scale to highlight the imaging of the notches toward the bottom of the calibration block. TFM is an amplitude-based image reconstruction algorithm where the A-scans from the FMC dataset are used to synthetically focus on every point in a defined region of interest.

Figure 3. FMC TFM results from a thick-film linear phased array installed on a calibration block.

Using other information from the FMC dataset, such as the phase of the waveforms, has proven to be beneficial in certain cases. At each focal point in the region of interest, a large phase coherence among all the waveforms can be indicative of a focused reflector. This can then be applied to the TFM image at each focal point as a weighting factor (also known as the Phase Coherence Factor (PCF)) to improve the signal-to-noise ratio. 

Figures 4 and 5 illustrate the results of applying the phase coherence imaging technique to the FMC datasets collected with thick film transducer arrays. The sample is a section of high-energy piping approximately 1.7 inches thick with cracking at various positions along a girth weld. The sample has a counterbore with ID-initiated cracks up to approximately 0.5 inches in length coming from the taper of the counterbore. The thick film transducer arrays were located at different positions along the weld.

SUMMARY
The energy industry is moving away from traditional scheduled-based planning for inspection and maintenance activities and toward “smart plant” concepts that rely more heavily on data correlated to actual component conditions. To accomplish this, there is a need for new and novel monitoring technologies that are both unobtrusive and able to withstand the harsh conditions of industrial facilities. Collecting robust and meaningful monitoring data will be critical in ensuring that safety and asset reliability are maintained and even improved. Structural Integrity’s thick-film UT technology has been developed to achieve this goal and continues to evolve for higher-temperature components and more advanced applications. We are ready to support a variety of in-field applications, contact one of SI’s experts if you have questions or a potential application that could benefit from installed thick-film UT sensors.

Figure 4. Phase coherence imaging result from a thick film transducer array on a cracked weld sample.

Figure 5. Phase coherence imaging result from a thick film transducer array on a cracked weld sample.

News & Views, Volume 53 | An ECA Process for the Impact of Hydrogen Blending on Girth Weld Defects

By:  Scott Riccardella, Owen Malinowski and Chris Tipple

Several pipeline operators have established pilot demonstration programs to blend hydrogen with natural gas (hydrogen blending) in their gas transmission pipelines.  Structural Integrity Associates (SI) has been providing clients technical consulting support to complete engineering critical assessment (ECA) projects to help evaluate the potential impact to pipeline integrity and help ensure the safety of the public, customers, employees, and the natural gas pipeline infrastructure. 

In a recent study, girth weld defects were identified as a key threat to pipeline integrity, particularly when the pipeline is exposed to large axial strain due to soil movement (which can be experienced from landslides, underwater erosion, storm surge, ground settlement and lateral spreading).  The impact to girth weld defects combined with large strain can pose a significant threat that is further exacerbated with hydrogen blending.  SI developed and implemented a program to complete a detailed ECA using probabilistic risk modeling to assess the probability of rupture (POR) to an offshore pipeline that had experienced significant strain due to erosion of the channel area, pipeline movement, and sand waves in the sea channel.  

To complete the ECA, a probabilistic analysis was performed consisting of the following activities:

REVIEW OF IN-LINE INSPECTION RESULTS

  • Recent strain data collected from an Inertial Mapping Unit (IMU) In-Line Inspection (ILI) tool were reviewed and analyzed to create a map of applicable strain at each girth weld in the study. 

MATERIAL PROPERTY, DEFECT AND OPERATING DATA ANALYSIS

  • Pipe populations were developed with specific characteristics that make them more compatible with hydrogen blending, or less compatible due to the respective susceptibility to hydrogen-related threats under different operating conditions.
  • SI developed Statistical distributions for key material properties (strength, toughness, wall thickness, etc.) and girth weld defect characteristics (length, depth, etc) using client specific and industry databases.
  • SI reviewed and incorporated relevant material tests performed to evaluate the effects of targeted hydrogen blend levels on the materials of interest (carbon steel base metal, longitudinal seam welds and girth welds).

DETERMINISTIC ANALYSIS USING FINITE ELEMENT MODELING (FEM)

  • A finite element analysis was utilized to determine the stress intensity factor of a circumferentially oriented crack subjected to high bending loads resulting in large axial strain.  The elastic-plastic analysis was used to determine the stress intensity factor as a function of strain, for a circumferentially oriented, externally breaking crack subject to a bending stress.

DEVELOPMENT OF A FRACTURE MECHANICS MODEL 9for probabilistic modeling)

  • From the FEA results a simplified elastic model was developed relating the stress intensity factor to the peak tensile axial strain resulting from bending.
  • SI incorporated the stress intensity factor from this model into an API 579 FAD based evaluation of girth weld, crack-like defects.

REVISIONS TO SI SYNTHESIS™ SOFTWARE

  • SI has developed specialized risk analysis software tools to evaluate pipeline POR which were applied to evaluate the impact or hydrogen blending to the POR. 
  • The software was specifically enhanced for this analysis to incorporate the following items:
    • Evaluation of flaws associated with circumferential cracking (such as those that may be encountered in vintage girth welds).
    • Incorporation of secondary loads and stresses (such as those encountered through land/soil movement).

PROBABILISTIC ANALYSIS

  • SI applied the probabilistic framework to evaluate the increased susceptibility to failure imposed from hydrogen blending with special consideration for ground movement and girth weld defects.  
  • This framework used Probabilistic Fracture Mechanics (PFM) and addressed the following phenomena associated with hydrogen blending:
    • Accelerated crack growth rates and 
    • Hydrogen embrittlement of the pipeline steel.
  • The POR was then evaluated for each active threat on the pipeline, comparing the risks associated with pure natural gas service to natural gas with hydrogen blending, considering various assessment options (hydrotest or ILI) prior to hydrogen injection.

CONCLUSION

Key challenges have been identified with blending hydrogen in gas transmission pipelines.  The susceptibility to failure of girth weld defects exposed to significant strain can be further exacerbated by the presence of hydrogen.  SI has developed a probabilistic framework and supporting tools to complete an ECA and provide a better understanding of the threats and subsequent impact to risk posed by cracks and crack-like defects in a hydrogen blending environment.

Get News & Views, Volume 53

News & Views, Volume 53 | Materials Lab Featured Damage Mechanism

CIRCUMFERENTIAL THERMAL FATIGUE IN CONVENTIONAL WATERWALL TUBES

By:  Wendy Weiss

Circumferential Thermal Fatigue damage in Conventional Waterwall Tubes most commonly appears as circumferentially oriented cracking in the waterwalls of coal-fired supercritical units. Initially, the formation of ripple magnetite was a significant factor in the formation of this damage. Later, the introduction of oxygenated treatment controlled the formation of ripple magnetite, thus greatly reducing this damage mechanism.  In the early 2000s, however, this type of thermal fatigue began occurring more frequently as low NOx burners and separated overfire air systems were introduced. 

Figure 1. Tube with a series of circumferential cracks

Mechanism

Three basic factors contribute to this type of thermal fatigue damage. 

  1. The first factor is the starting tube temperature (i.e., the temperature under normal operating conditions). The higher the starting temperature, the greater the accumulation of damage in the affected tubing. For example, tubes subjected to higher heat flux or tubes with thick weld overlays will have higher average metal temperatures and accumulate damage more quickly. 
  2. The second factor is the extent of gradually increasing tube temperature caused by reasons such as internal deposit buildup, flame impingement, or unstable flow. 
  3. The third factor is the contribution of thermal transients due to slag shedding or using sootblowers or water cannons. 

Essentially, the thermal fatigue cracking results from the combination of increasing tube metal temperature and thermal transients and is aggravated by high starting tube temperatures. 

Figure 2. The external surface of the tube after the external deposits were removed

TYPICAL LOCATIONS

Figure 3. Cross-sectional views of the cracking in the etched (Top) and unetched (BOTTOM) conditions

  • Tubes with slag buildup and shedding
  • Areas affected by wall blow quenching
  • High heat flux locations
  • Areas affected by flame impingement
  • Cracking can be localized or widespread
  • Tends to be contained within a relatively narrow range of elevations

FEATURES

  • Circumferentially oriented, multiple, parallel cracks along the hot side of the tubes.
  • Notch shaped, oxide filled cracks in cross-section.
  • Adjacent tubes can exhibit variability in crack density.

ROOT CAUSES

  • High Initial Waterwall Tube Temperatures
    • Thick weld overlays
    • Higher heat flux
    • Flame impingement
  • Increasing Waterwall Tube Temperatures
    • Internal deposits including ripple magnetite, thick oxide layers, or feedwater corrosion products
    • Reduced internal flow rate
    • Formation of external oxides and deposits
  • Severe Thermal Transients
    • Natural or forced slag removal, including slag shedding and sootblowing
    • Use of water cannons or improper sootblowing
    • Flame instabilities
    • Unit operation, including forced fan cooling, rapid startups, frequent load cycling

Get News & Views, Volume 53

News & Views, Volume 53 | PEGASUSTM Nuclear Fuel Code

ADVANCED FUEL MODELING DEVELOPMENT STATUS

By:  Bill Lyon

INTRODUCTION
The PEGASUS nuclear fuel behavior code features a robust 3D, finite element modeling (FEM) computational foundation capable of performing both thermo-mechanical and structural non-linear analyses within a highly versatile and customizable computational platform. The first applications of PEGASUS were for light water reactor (LWR) fuels and materials. Development work on PEGASUS has been extended to advanced fuel designs such as those proposed for Advanced Technology Fuel (ATF) LWR applications and Gen IV reactor designs, including gas and liquid metal-cooled reactors (GCRs and LMRs). 

The versatility and adaptability of PEGASUS is key in enabling extensions to non-conventional operating environments, materials, fuel forms, and geometries.

LWR APPLICATIONS
SiC Cladding

Figure 1. SiGA cladding is a multi-layered composite design composed of SiC fiber in a SiC matrix.

A project is underway to further the development and irradiation testing of a composite silicon carbide matrix as an ATF cladding material. This research is supported through a DOE Funding Opportunity award (DE-FOA-0002308) for the irradiation of a composite silicon-carbide (SiC) ceramic matrix material in an existing U.S. commercial LWR. This work is led by General Atomics – Electromagnetic Systems (GA-EMS) with Structural Integrity Associates (SIA) as a primary subcontractor. For this work, PEGASUS is being adapted to model monolithic and composite SiC manufactured by GA-EMS, SiGA [1], through the incorporation of proprietary material constitutive models. PEGASUS will then be used to provide independent test performance analyses aiding in the design of the irradiation vehicle and predicted material performance. The goal of the testing is to gather irradiation data under prototypic LWR operating conditions and to inform and confirm material performance models for the SiGA-based cladding. A follow-on activity is planned to evaluate the predicted performance compared to data gathered during the post-irradiation examination phase of the project.

Figure 2. Lightbridge Fuel Design PEGASUS Models

Cruciform Metallic Fuel
An additional fuel concept that has been explored using PEGASUS is a cruciform, extruded metallic fuel design proposed by Lightbridge Corporation [2]. This fuel is characterized by a unique multi-lobed fuel cross-section and features a U-50Zr fuel composition. Recent work has been published on fabrication testing of this proposed fuel design by Pacific Northwest National Laboratory (PNNL) [3]. PEGASUS has been used previously to prototype 2D and 3D geometric models and meshes of Lightbridge fuel and to perform fundamental temperature and stress distributions for this fuel under prototypic LWR conditions. PEGASUS has specific modeling tools designed to facilitate “extruded” 3D fuel designs that automate the meshing of these geometries. More work in this area is planned as a proposal has recently been awarded under the DOE NEUP program (DE-FOA-0002732) funding a collaborative project led by Texas A&M University along with Lightbridge, NuScale, and Structural Integrity Associates, Inc. (SI) for modeling this type of fuel for application in a LWR SMR.

 

URANIUM METAL ALLOY FUELS FOR SODIUM-COOLED FAST REACTORS
The initial implementation of metallic alloy fuel and stainless-steel cladding material constitutive models for prototypic fast reactor fuel designs in PEGASUS has been completed. Material properties and behavioral models for U-Pu-Zr fuel and HT-9 (high Chromium, martensitic stainless steel) cladding have been added. Ongoing work includes the implementation of a gaseous swelling and fission gas release behavior model for U-Pu-Zr fuel, a Zr-redistribution model, and a fuel-cladding chemical interaction (FCCI) model that includes the effect on cladding wall thinning.

To test the implementation of these models, benchmark tests were prepared that provided comparative data for assessment of the models’ performance. Test cases were chosen from two experimental series irradiated in EBR-II: X430, a 37-pin hexagonal sub-assembly, and X441, a 61-pin bundle. These experiments were designed to test numerous fuel rod design variables and fuel response as a function of fuel alloy composition, smear density, plenum-to-fuel volume ratio, power, and coolant conditions [4]. The general experimental fuel rod design corresponds to the typical driver fuel configuration shown in Figure 3.

Figure 3. Typical EBR-II Mark-III or Mark-IIIA Fuel Element [5]

Figure 4. Left: 2D Computational Model of Rod DP2, Right: Temperature Contour Plot of the Fuel Stack Region for Rod DP21 at Peak Power (plenum region removed for detail)

An illustration of the model and selected results from the initial analysis of rod DP21, assembly X441 are shown in the figures above. Figure 4 provides a diagram of the computational model showing the primary components of the model and a plot of the temperature distribution throughout the fueled region of the rod at peak power. Figure 5 provides the radial temperature profile across the fuel rod from the center to the cladding outer surface at peak power near the end of the irradiation period. Temperatures vary from just ~900 K at the pellet center to ~650 K at the cladding surface. The temperature differential is fairly low at ~250 K, as would be expected from a high-conductivity metal fuel rod with a Na-bonded fuel cladding gap. These results are consistent with published experimental observations.

TRISO FUEL MODELING DEVELOPMENT
Several advanced fuel material models have been implemented specifically for TRISO fuel in PEGASUS, including thermal and mechanical models for UCO or UO2 kernels, PyC, SiC materials, and a fission gas release model for computing the release of gaseous fission products such as Xe and Kr. In addition to the standard 3D and 2D axisymmetric modeling FEM capabilities in the code, PEGASUS contains several unique tools designed specifically to support TRISO fuel modeling and analysis. These include a “spherical mesh object” tool that can automate the process of generating 2D/3D TRISO spheres, meshing them, and embedding them into a fuel matrix to allow modeling of individual TRISO kernels or fully encapsulated TRISO fuel forms. An example of models generated using the spherical mesh object tool is shown in Figure 6. The spherical mesh object capability is, to our knowledge, unique to PEGASUS and not found in any other fuel performance or general-purpose FEM code. PEGASUS also has a “reshape” function that can automate the process of meshing and modeling deformed TRISO particles to increase user efficiency. Figure 7 illustrates particle meshes that were created using the reshape meshing tool.

These modeling capabilities allow PEGASUS to be used to investigate very detailed mechanical and structural effects in TRISO fuel forms. For example, enabling the detailed analysis of the mechanical interaction between TRISO fuel layers explicitly examining the effects of cracking, debonding, and asphericity within whole or damaged particles.

Figure 5. Radial Temperature Distribution Across the Fuel Rod Model at ~ 486 Days of Irradiation

Planned future development work includes the integration of damage-mechanics modeling and fission product diffusion in the TRISO particle, fuel compact, ad matrix. One failure mode of particular interest that has been identified is cracking of the IPyC layer which propagates through the SiC outer layer. This can create a pathway for enhanced fission product release from the TRISO particle to the surrounding fuel matrix. This failure mechansim appears to occur when the buffer layer remains bonded to the IPyC layer providing the conditions for a synergistic mechanical and chemical failure mechanism that combines cracking, stress concentration, and chemical corrosion (localized Pd-induced corrosion in the SiC [6]. This failure mode is of interest because it can have a strong impact on fuel source term determination for operational TRISO fuel.

Figure 6. Temperature distribution in a cross-section of a 3D slab of a TRISO compact matrix model with a “sparse”, random kernel distribution under prototypic gas-cooled reactor conditions. (Generated using the “spherical mesh object” tool.)

SUMMARY
PEGASUS is an advanced analysis tool developed for industry applications that can provide a complimentary and independent capability for nuclear fuel performance. Recent development work on PEGASUS has focused on expanding the applicability of the code to the advanced fuel (ATF) and advanced reactor arena. Future development is planned for PEGASUS that will continue along multiple avenues with an emphasis on advanced fuels and specific thermo-mechanical issues within the industry, such as deterministic failure model development. One example of this is the aforementioned Pd-induced failure mechanism identified for TRISO fuel. SI is actively seeking partners within the advanced fuel community to collaborate with on this work and would welcome inquiries and proposals for expanded application of PEGASUS.

 

Figure 7. Deformed 3D TRISO particle meshes generated using the “reshape” function tool in PEGASUS.

References

  1.   C. P. Deck et al., Overview of General Atomics SiGA™ SiC-SiC Composite Development for Accident Tolerant Fuel, Transactions of the American Nuclear Society, Vol. 120, Minneapolis, Minnesota, June 9–13, 2019.
  2.   J. Malone, A. Totemeier, N. Shapiro, and S.  Vaidyanathan, Lightbridge Corporation’s Advanced Metallic Fuel, Nuclear Technology, Vol 180, Dec. 2012.
  3.  Z. Huber and E. Conte, Casting and Characterization of U-50Zr, PNNL-33873, Pacific Northwest National Laboratory, Richland, Washington 99354, Jan. 2023.
  4.  C. E. Lahm, J. F. Koenig, R. G. Pahl, D. L. Porter, and D. C. Crawford, “Experience with Advanced Driver Fuels in EBR-II,” J. Nucl. Mat. 204 (1993) 119-123.
  5. G. L. Hofman, M. C. Billone, J. F. Koenig, J. M. Kramer, J. D. B. Lambert, L. Leibowitz, Y. Orechwa, D. R. Pedersen, D. L. Porter, H. Tsai, and A. E. Wright. Metallic fuels handbook, Technical Report ANL-NSE-3, Argonne National Laboratory, 2019.
  6. J.D Hunn, C.A. Baldwin, T.J. Gerzak, F.C. Montgomery, R.N. Morris, C.M. Silva, P.A Demkowicz, J.M. Harp, and S.A. Ploger, Detection and Analysis of Particles with Failed SiC in AGR-1 Fuel Compacts, Nuclear Engineering and Design, April 2016.

Get News & Views, Volume 53