Results for “INTERNAL STRESS” — Page 14

DOE OSTI2021

“Seeing Is Believing”—In-Depth Analysis by Co-Imaging of Periodically-Poled X-Cut Lithium Niobate Thin Films

Nonlinear and quantum optical devices based on periodically-poled thin film lithium niobate (PP-TFLN) have gained considerable interest lately, due to their significantly improved performance as compared to their bulk counterparts. Nevertheless, performance parameters such as conversion efficiency, minimum pump power, and spectral bandwidth strongly depend on the quality of the domain structure in these PP-TFLN samples, e.g., their homogeneity and duty cycle, as well as on the overlap and penetration depth of domains with the waveguide mode. Hence, in order to propose improved fabrication protocols, a profound quality control of domain structures is needed that allows quantifying and thoroughly analyzing these parameters. In this paper, we propose to combine a set of nanometer-to-micrometer-scale imaging techniques, i.e., piezoresponse force microscopy (PFM), second-harmonic generation (SHG), and Raman spectroscopy (RS), to access the relevant and crucial sample properties through cross-correlating these methods. Based on our findings, we designate SHG to be the best-suited standard imaging technique for this purpose, in particular when investigating the domain poling process in x-cut TFLNs. While PFM is excellently recommended for near-surface high-resolution imaging, RS provides thorough insights into stress and/or defect distributions, as associated with these domain structures. In this context, our work here indicates unexpectedly large signs for internal fields occurring in x-cut PP-TFLNs that are substantially larger as compared to previous observations in bulk LN.

Reitzig, Sven (ORCID:0000000286970379)↗

DOE OSTI2021

On the origin of internal obstacles to dislocation glide in single-phase NiFe random alloys

Over the past decade, a new class of unconventional alloys, where several constitutive elements are presented in close to equal concentrations, is considered as promising structural materials in many applications. Properties and behavior of many concentrated solid solution alloys (CSSAs) under different conditions, including mechanical loading, differs from that of conventional alloys and has been the subject of intensive studies by different techniques. A key feature of these materials is that there is no clear distinction between solute and solvent species of atoms, which is the conventional approach to considering dislocation motion in alloys. To understand this regime, we have recently reported on the glide of a screw dislocation glide in a simple equiatomic NiFe alloy and reported the two-mode character. At low applied stress, the dislocation motion is similar to a glide through a field of obstacles. In this work, we report the results of a detailed study on the nature of these internal obstacles. We demonstrate that these are not conventional localized “obstacles” but is the result of significantly different energies of instantaneous dislocation configurations over the alloy bulk. The strong pinning positions are associated with composition changes in the partial dislocation cores, while observing no significant composition changes along the stacking faults. Local energies can vary significantly, as demonstrated using simulations of short-periodicity dislocation lines. Therefore, the critical stress to move a dislocation under such conditions is not a constant Peierls stress but represents a wide spectrum of barriers between different non-straight valleys that the dislocation line accommodates in the bulk at particular stress/strain conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A method for assessing economic, environmental, and reliability tradeoffs of interregional transmission connecting ERCOT (the Texas grid) to the eastern and western grids

Reliable development of the power grid is an evolving concern for humanity due to extreme weather that frequently threatens power sector infrastructure. The state of Texas is a uniquely structured testbed for grid planners to study when looking for solutions to development, innovation, and overcoming such challenges. Because of its size and islanded structure, Texas is small enough to model, but big enough to matter. Texas is a global leader in energy production, energy consumption, and maintains an unusually diverse fuel mix. In addition, the state has experienced winter freezes, heat waves, wind storms, droughts and floods that have threatened power sector infrastructure or caused recent blackouts and calls for demand side conservation. One of the most devastating of these events was the North American winter storm, dubbed “Winter Storm Uri” by the Weather Channel, that froze the region in February 2021 and led to an extended power outage event that put the majority of Texan residents in darkness for days. While preparing to avoid such outage events in the future, various tools have been proposed to improve grid reliability, including energy efficiency, demand response, and distributed energy resources. An additional option would be to develop interregional transmission that connects the Texas grid to other national grids. To assess the merits of this idea, we developed a novel, universally-applicable and internationally-relevant framework to study how the Texas grid would evolve alongside access to various interregional ties. This method allows us to stress the synthetic grid structure and analyze how it would respond to the shock of a simulated winter storm event. Our method leverages open-source modeling tools, such as PowerGenome, pyGRETA, and GenX to synthesize unique zonal grid data, construct a consolidated network of model regions, and simulate different developmental pathways of capacity expansion and operational dispatch. We demonstrate our method with an analysis connecting the Electric Reliability Council of Texas (ERCOT), the grid that serves most of Texas, the Western Electricity Coordinating Council (WECC), the grid that serves the western half of the contiguous U.S., and the Eastern Interconnect, the grid that serves the eastern half of the contiguous U.S. Our results indicate that the cost-optimal capacity of interregional transmission connecting the ERCOT grid to other grids lies between 9–13 GW assuming baseline conditions. Building this amount of connecting capacity in one or multiple directions lowers the costs and emissions of development and operation by up to $16 billion and 257 million metric tonnes (MMT) respectively. Additionally, our results show that the interregional connections between ERCOT and other national grids reduce the amount of total load shed required through mild winter storm events. However, our results also show that there is a threshold of very extreme winter storm conditions, spanning multiple service areas, above which the connections exacerbate resource adequacy problems. Therefore, the results indicate that the connections need to be carefully planned alongside the rest of the grid infrastructure to avoid over-reliance on specific resources or technology options.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Mechanical Properties of Nominally Ductile and Brittle Zones within the Caney Shale Formation

We present a case study from the Caney Shale, which is an emerging hydrocarbon-bearing formation in Oklahoma. Open hole logs lead to a nominal distinction of five zones in an exploration well, two of which are designated “ductile” and three of which are designated as less ductile “reservoir” zones. Subsequently, a detailed mechanical characterization was carried out on core plugs taken from each of these zones. Results show nominally “ductile” zones are weaker in terms of UCS, tensile strength, fracture toughness, and confined compressive strength. However, by a variety of brittleness indices, it is found the nominal “ductile” regions do not behave as expected. They are weaker, but not consistently more ductile. Furthermore, some brittleness indices lead to contradictory conclusions. For example, zones which would be considered the most brittle according to the proportion of nonlinear stress-strain behavior prior to failure and/or by indices based on elastic properties would be considered the least brittle according to their internal friction angle.

Benge, Margaret↗

DOE OSTI2023

Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity

This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.

42 ENGINEERING↗

DOE OSTI2026

A comparison of wrought and powder metallurgical FeCrAl claddings under simulated LWR accident transients

Iron-chromium-aluminum (FeCrAl) alloys are potential accident tolerant fuel (ATF) cladding candidates for light-water reactors but are difficult to fabricate as thin-walled tubes via conventional cast-and-wrought routes. Powder metallurgy (PM) offers a manufacturing alternative with improved compositional control, but its transient accident performance has not been directly benchmarked against wrought variants. This study evaluates the burst behavior of commercially developed PM-processed FeCrAl alloys, PM-C26M (Fe-12Cr-6Al-2Mo) and the high precipitate density FA-SMT (Fe-22Cr-5Al-3Mo), under simulated light-water reactor accident transient conditions. Burst testing was conducted using the Severe Accident Test Station with heating rates of 5 °C/s and 50 °C/s and internal pressures ranging from 25 MPa to 100 MPa. PM-C26M reproduced wrought C26M burst behavior within 7–37 °C across the stress range, indicating that PM processing does not compromise transient strength. FA-SMT exhibited markedly higher burst temperatures and reduced heating-rate sensitivity, consistent with its engineered precipitate strengthening. FA-SMT rupture exhibited axial "unzipping" rather than the lateral tearing characteristic of PM- and wrought C26M. Post-test EBSD and fractography indicate that this behavior is strongly correlated with strain-gated intergranular void nucleation associated with the dense precipitate architecture of FA-SMT, a response absent in the comparatively clean PM-C26M matrix and consistent with rupture morphologies reported for oxide-dispersion strengthened (ODS) FeCrAl of similar base-matrix chemistry to PM-C26M. These findings highlight the potential of powder metallurgy as a viable fabrication route for ATF claddings from an accident performance standpoint.

Bell, Sam [ORNL] (ORCID:0000000251905657)↗

DOE OSTI2021

Magneto mitochondrial dysfunction mediated cancer cell death using intracellular magnetic nano-transducers

Mitochondria are crucial regulators of the intrinsic pathway of cancer cell death. The high sensitivity of cancer cells to mitochondrial dysfunction offers opportunities for emerging targets in cancer therapy. Herein, magnetic nano-transducers, which convert external magnetic fields into physical stress, are designed to induce mitochondrial dysfunction to remotely kill cancer cells. Spindle-shaped iron oxide nanoparticles were synthesized to maximize cellular internalization and magnetic transduction. The magneto-mechanical transduction of nano-transducers in mitochondria enhances cancer cell apoptosis by promoting a mitochondrial quality control mechanism, referred to as mitophagy. Furthermore, in the liver cancer animal model, nano-transducers are infused into the local liver tumor via the hepatic artery. After treatment with a magnetic field, in vivo mitophagy-mediated cancer cell death was also confirmed by mitophagy markers, mitochondrial DNA damage assay, and TUNEL staining of tissues. This study is expected to contribute to the development of nanoparticle-mediated mitochondria-targeting cancer therapy and biological tools, such as magneto-genetics.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Irreversible thermodynamics of curved lipid membranes. II. Permeability and osmosis

We present a theory that combines the framework of irreversible thermodynamics with modified integral theorems to model arbitrarily curved and deforming membranes immersed in bulk fluid solutions. We study the coupling between the mechanics and permeability of a viscous and elastically bendable membrane, and a multicomponent bulk fluid solution. An equation for the internal entropy production for irreversibilities at the membrane is derived, determining the generalized thermodynamic forces and fluxes, from which we identify the deviatoric stress as a novel driving force for permeability. A complete set of equations of motion, constitutive laws, and boundary conditions to model the lipid membrane and bulk fluid system are provided.

Alkadri, Ahmad M↗

DOE OSTI2023

Low temperature absolute internal quantum efficiency of InGaN-based light-emitting diodes

A method to reveal low temperature absolute internal quantum efficiency of an InGaN-based light-emitting diode (LED) is reported. Structural (via scanning electron and atomic force microscopies, x-ray diffraction, and Raman spectroscopy) and optical (via temperature- and power-dependent photoluminescence) properties of InGaN-based blue LEDs are coupled with a channel-based recombination model, which reveals the low temperature absolute internal quantum efficiency of the LED on traditional sapphire and emerging Si (111) substrates as 27.5% and 71.1%, respectively. High low-temperature absolute internal quantum efficiency of the former is attributed to its lower defect density and less strain. Furthermore, our results show assuming perfect (i.e., 100%) internal quantum efficiency in an InGaN-based LED at low temperatures could be misleading, and using negative thermal quenching properties of defect luminesce bands with a channel-based recombination model can be a method to quantify LED's less than ideal low temperature internal quantum efficiencies.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Thermal stress in degraded forests in the Brazilian Amazon Arc of Deforestation

Understanding thermal stress in tropical forests has taken on new urgency in light of accelerating climate change and expansion of deforestation and forest degradation. Degraded tropical forests in particular may be approaching critical temperature thresholds even more rapidly than intact forests, with implications for tree survival and ecosystem recovery. We investigate thermal stress in degraded tropical forests within the Brazilian Amazon Arc of Deforestation. Using land surface temperature data from the ECOsystem Spaceborne Thermal Radiometer Experiment on the international Space Station (ECOSTRESS), we compared canopy temperatures of intact, selectively logged, and burned forests in Feliz Natal, Mato Grosso, Brazil. Upper canopy temperatures in previously burned forests were 4.1% higher (mean = 36.5 °C) and 50.9% more variable compared to intact and logged forests, which showed remarkably similar temperature distributions (means of 34.9 °C and 35.1 °C, respectively). Modeled leaf temperature distributions based on canopy temperature measurements from one of the warmest days in a 2 year dry season record indicated that 87% of leaves in the warmest burned forest patches exceeded the temperature threshold where respiration surpasses photosynthesis, compared to approximately 72%–74% in intact and logged forests, respectively, in the same time period. After controlling for environmental factors, burned forests were predicted to be 2.6 °C warmer on average [1.39 °C–3.96 °C, 95% credible interval] than intact forests across a 5–40 m canopy height range. Burned forests showed modest thermal recovery over time, with temperatures decreasing by approximately 1.2 °C over a 30 year recovery period. In contrast, logged forests showed minimal thermal differences from intact forests (−0.353 on average, [−1.125–0.274, 95% credible interval]) and negligible change (0.15 °C) over the same timeframe. While the absolute probabilities of exceeding damaging thermal thresholds remain low across all forest types under current climate conditions, the probability of leaves exceeding temperatures that cause permanent leaf damage was ten times higher in burned forests, with implications for the future of burned forest regeneration in water-limited regions of the Amazon basin. In particular, the combination of higher mean temperatures, greater temperature variability, and more frequent exposure to damaging thermal thresholds implies that burned tropical forests will experience substantially higher mortality rates and slower biomass recovery compared to intact and selectively logged forests, especially in water-limited regions where trees cannot rely on evaporative cooling to moderate canopy temperatures.

Brazilian Amazon Arc of Deforestation↗

DOE OSTI2024

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Controllable Solar Flux Heating for Freeze Recovery in Molten Salt Parabolic Trough Collectors

In molten-salt parabolic trough plants, the melting process is particularly important for freeze recovery of salt that is solidified in a collector loop, should such an event occur. Currently impedance heating is expected for freeze recovery of the collector loops, but this method can be expensive. Here, a lower-cost alternative is proposed to use controllable concentrated solar flux directly from the parabolic mirrors to thaw salt that is frozen in the collector. A computational fluid dynamics model was developed to explore the solidification and melting processes of molten salt in a parabolic trough receiver and to assess the viability of this concept. Results indicate that concentrated solar heating has the potential to melt frozen salt in 5.6 h, compared to that in 8.8 h for a 300 W m –1 impedance heating system. At the same time, controllable solar flux heating introduces nonuniform solar fluxes on the receiver surface, which can induce significant thermal stress on the receiver tube. A preliminary stress analysis indicates that the temperature difference across the receiver tube should be maintained below about 70 °C for heating up to 300 °C at internal pressures ≤10 bar. At these conditions, freeze recovery using solar flux heating will not significantly affect receiver lifetime. These results suggest that controllable solar flux heating could effectively supplement or replace impedance heating in the freeze recovery system. Incorporating this methodology in future parabolic trough concentrating solar power plants is an opportunity for capital and operational cost-savings.

14 SOLAR ENERGY↗

DOE OSTI2023

Hydrogen Extremely Low Probability of Rupture Version 1

SAND2024-08511O The Hydrogen Extremely Low Probability of Rupture (HELPR) software toolkit provides probabilistic fracture assessments of natural gas infrastructure for transporting hydrogen. Using probabilistic sampling and sensitivity analysis—along with engineering fatigue and fracture mechanics— the software develops models that can predict structural integrity measures for natural gas pipelines under stress. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Schroeder, Benjamin↗

DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗

DOE OSTI2020

Coupling between TOUGH3 and FLAC3D

Sandia National Laboratories has hired Itasca Consulting Group, Inc., the authors of the FLAC3D geomechanics software, to couple FLAC3D with TOUGH3, the porous media flow solver. The work is being done to enable a coupled mechanical-thermal-hydraulic analysis of a potential criticality event in a dual purpose cannister (DPC). The U.S. Department of Energy Office of Spent Fuel and Waste Science & Technology is investigating the performance of DPCs for direct geological disposal of spent nuclear fuel. Post closure criticality control is an important aspect of this investigation. Over geological timescales, it is envisioned that the canister and canister overpack will develop fractures due to stress corrosion processes. A breach in the canister could allow groundwater to fill the canister. Fresh water is a neutron moderator; thus, if the canister internals and fuel assemblies have been sufficiently degraded, a criticality event could occur. Such an event would release enough energy to boil the water between the fuel rods and pressurize the cannister. This internal pressurization may cause the initial fractures in the canister and overpack to grow. It is important to understand the change in hydraulic transmissivity between the canister and surroundings for two reasons: first, because it may control the potential for and frequency of subsequent criticality events; second, because it will control the release of radionuclides from the canister. The motivation for this work is to better understand the potential for periodic criticality events, cannister damage, and release of radionuclides during a criticality event in a DPC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Multi physics simulation on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) by using Cardinal

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗