Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Multi-approach assessment of the spatial distribution of the specific yield: application to the Crau plain aquifer, France

Spatially distributed values of the specific yield, a fundamental parameter for transient groundwater mass balance calculations, were obtained by means of three independent methods for the Crau plain, France. In contrast to its traditional use to assess recharge based on a given specific yield, the water-table fluctuation (WTF) method, applied using major recharging events, gave a first set of reference values. Then, large infiltration processes recorded by monitored boreholes and caused by major precipitation events were interpreted in terms of specific yield by means of a one-dimensional vertical numerical model solving Richards’ equations within the unsaturated zone. Finally, two gravity field campaigns, at low and high piezometric levels, were carried out to assess the groundwater mass variation and thus alternative specific yield values. The range obtained by the WTF method for this aquifer made of alluvial detrital material was 2.9– 26%, in line with the scarce data available so far. The average spatial value of specific yield by the WTF method (9.1%) is consistent with the aquifer scale value from the hydro-gravimetric approach. In this investigation, an estimate of the hitherto unknown spatial distribution of the specific yield over the Crau plain was obtained using the most reliable method (the WTF method). A groundwater mass balance calculation over the domain using this distribution yielded similar results to an independent quantification based on a stable isotope-mixing model. This agreement reinforces the relevance of such estimates, which can be used to build a more accurate transient hydrogeological model.

58 GEOSCIENCES↗
DOE OSTI2019

The transient response of atmospheric and oceanic heat transports to anthropogenic warming

Model projections of the near-future response to anthropogenic warming show compensation between meridional heat transports by the atmosphere (AHT) and ocean (OHT) that are largely symmetric about the equator1,2,3, the causes of which remain unclear. Here, using both the Coupled Model Intercomparison Project Phase 5 archive and Community Climate System Model version 4 simulations forced with Representative Concentration Pathway 8.5 to 2600, we show that this transient compensation—specifically during the initial stage of warming—is caused by combined changes in both atmospheric and oceanic circulations. In particular, it is caused by a southward OHT associated with a weakened Atlantic Meridional Overturning Circulation, a northward apparent OHT associated with an ocean heat storage maximum around the Southern Ocean, and a symmetric coupled response of the Hadley and Subtropical cells in the Indo-Pacific basin. It is further shown that the true advective OHT differs from the flux-inferred OHT in the initial warming due to the inhomogeneous responses of ocean heat storage. These results provide new insights to further our understanding of future heat transport responses, and thereby global climatic processes such as the redistribution of ocean heat.

He, Chengfei↗
DOE OSTI2023

Improved representation of laminar and turbulent sheet flow in subglacial drainage models

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants: Model-Predictive Control Library Development (Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2013

Gearbox and Drivetrain Models to Study Dynamic Effects of Modern Wind Turbines: Preprint

Wind turbine drivetrains consist of components that directly convert kinetic energy from the wind to electrical energy. Guaranteeing robust and reliable drivetrain designs is therefore important to minimize turbine downtime. Current drivetrain models often lack the ability to model both the impacts of electrical transients as well as wind turbulence and shear in one package. In this work, thecapability of the FAST wind turbine computer-aided engineering tool, developed by the National Renewable Energy Laboratory, is enhanced through integration of a dynamic model of the drivetrain. The dynamic drivetrain model is built using Simscape in the MATLAB/Simulink environment and incorporates detailed electrical generator models. This model can be used in the future to test advanced controlschemes to extend life of the gearbox.

17 WIND ENERGY↗
DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

The pros and cons of beyond standard model interpretations of ANITA events

The Antarctic Impulsive Transient Antenna (ANITA) experiment has observed two air shower events with energy ~500 PeV emerging from the Earth with exit angles ~30⁰ above the horizon. As was immediately noted by the ANITA Collaboration, these events (in principle) could originate in the atmospheric decay of an upgoing τ-lepton produced through a charged current interaction of a ντ inside the Earth. However, the relatively steep arrival angles of these perplexing events create tension with the standard model (SM) neutrino-nucleon interaction cross section. Deepening the conundrum, the IceCube neutrino telescope and the Pierre Auger Observatory with substantially larger exposures to cosmic ντ's in this energy range have not observed any events. This lack of observation implies that the messenger particle (MP) giving rise to ANITA events must produce an air shower event rate at least a factor of 40 larger than that produced by a flux of τ-neutrinos to avoid conflicts with the upper limits reported by the IceCube and the Pierre Auger collaborations. In addition, the sensitivity of ANITA to MP-induced events must be comparable to or larger than those of IceCube and Auger to avoid conflict with the non-observation of any signal at these facilities. Beyond SM interpretations of ANITA events can be classified according to whether the MPs: (i) live inside the Earth, (ii) originate in neutrino-nucleon collisions inside the Earth, (iii) come from cosmological distances. In this communication we investigate the positive and negative facets of these three classes of models.

None, None↗
DOE OSTI2025

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Divergence in land surface modeling: linking spread to structure

Divergence in land carbon cycle simulation is persistent and widespread. Regardless of model intercomparison project, results from individual models diverge significantly from each other and, in consequence, from reference datasets. Here we link model spread to structure using a 15-member ensemble of land surface models from the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) as a test case. Our analysis uses functional benchmarks and model structure as predicted by model skill in a machine learning framework to isolate discrete aspects of model structure associated with divergence. We also quantify how initial conditions prejudice present-day model outcomes after centennial-scale transient simulations. Overall, the functional benchmark and machine learning exercises emphasize the importance of ecosystem structure in correctly simulating carbon and water cycling, highlight uncertainties in the structure of carbon pools, and advise against hard parametric limits on ecosystem function. We also find that initial conditions explain 90% of variation in global satellite-era values—initial conditions largely predetermine transient endpoints, historical environmental change notwithstanding. As MsTMIP prescribes forcing data and spin-up protocol, the range in initial conditions and high levels of predetermination are also structural. Our results suggest that methodological tools linking divergence to discrete aspects of model structure would complement current community best practices in model development.

58 GEOSCIENCES↗
DOE OSTI2013

New Automotive Air Conditioning System Simulation Tool Developed in MATLAB/Simulink

Further improvements in vehicle fuel efficiency require accurate evaluation of the vehicle's transient total power requirement. When operated, the air conditioning (A/C) system is the largest auxiliary load on a vehicle; therefore, accurate evaluation of the load it places on the vehicle's engine and/or energy storage system is especially important. Vehicle simulation software, such as 'Autonomie,' has been used by OEMs to evaluate vehicles' energy performance. A transient A/C simulation tool incorporated into vehicle simulation models would also provide a tool for developing more efficient A/C systems through a thorough consideration of the transient A/C system performance. The dynamic system simulation software Matlab/Simulink was used to develop new and more efficient vehicle energy system controls. The various modeling methods used for the new simulation tool are described in detail. Comparison with measured data is provided to demonstrate the validity of the model.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Uncertainty Quantification of Extratropical Forest Biomass in CMIP5 Models over the Northern Hemisphere

Abstract Simplified representations of processes influencing forest biomass in Earth system models (ESMs) contribute to large uncertainty in projections. We evaluate forest biomass from eight ESMs outputs archived in the Coupled Model Intercomparison Project Phase 5 (CMIP5) using the biomass data synthesized from radar remote sensing and ground-based observations across northern extratropical latitudes. ESMs exhibit large biases in the forest distribution, forest fraction, and mass of carbon pools that contribute to uncertainty in forest total biomass (biases range from −20 Pg C to 135 Pg C). Forest total biomass is primarily positively correlated with precipitation variations, with surface temperature becoming equally important at higher latitudes, in both simulations and observations. Relatively small differences in forest biomass between the pre-industrial period and the contemporary period indicate uncertainties in forest biomass were introduced in the pre-industrial model equilibration (spin-up), suggesting parametric or structural model differences are a larger source of uncertainty than differences in transient responses. Our findings emphasize the importance of improved (1) models of carbon allocation to biomass compartments, (2) distribution of vegetation types in models, and (3) reproduction of pre-industrial vegetation conditions, in order to reduce the uncertainty in forest biomass simulated by ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗
DOE OSTI2025

Assessing High Burnup U-19Pu-10Zr Fuel Performance against Historical and Modeled Behavior

Advancing the deployment of sodium-cooled fast reactors (SFRs) requires thorough testing of metallic fuel pins under accident conditions to establish safe operational limits of high burnup fuel. To conduct transient testing, a comprehensive understanding of steady-state fuel behavior obtained through both experimental characterization and accurate predictive capabilities is needed. This study comparatively assesses the steady-state irradiation performance of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, irradiated under prototypic fast reactor conditions in preparation for planned safety testing at the Transient Reactor Test Facility. Since DP-40 was designated for use in the test and DP-36 serves as its sibling pin, non-destructive, engineering-scale post-irradiation examinations (PIE) were conducted on both pins while destructive examinations were performed exclusively on DP-36. The results were then assessed against historical performance data from similar fuel pins irradiated in the Experimental Breeder Reactor-II. Additionally, the steady-state irradiation of each pin was modeled using the BISON fuel performance code to assess the accuracy of current modeling capabilities in predicting the baseline irradiation behavior. Non-destructive examinations included neutron radiography to measure fuel column elongation, gamma scanning to verify pin integrity and fission product migration, and profilometry to assess dimensional changes. Benchmarking against existing PIE data revealed consistent patterns in axial fuel column growth and cladding diametral strain, though both pins exhibited longer low-density “fluff” structures, which can have implications for core reactivity and source term calculations. Destructive examinations on DP-36 included fission gas release analysis and sectioning for optical microscopy, which showed more complex constituent redistribution patterns than the traditionally accepted 3-ring model. The axial evolution of fractional areas and porosities of each of the redistributed zones were quantified and presented. Modeling comparisons showed agreement in fractional fission gas release but consistently overestimated axial and radial swelling and disagreed with measured axial porosity patterns. These conservative overpredictions suggested that the pins would appear closer to failure or operational limits at the start of transient tests, potentially leading to higher strain accumulation during the transient. While conservative estimates provide safety margins, they can negatively impact fuel economics. A review of the swelling models identified areas for improvement in the gaseous swelling, solid swelling, and fuel hot-pressing models when applied to ternary fuel. The results of this study highlight the critical importance of conducting pre-test characterization on both test and sibling pins to accurately capture steady-state fuel behavior, providing a precise baseline for post-test evaluations and essential inputs for transient modeling of the planned experiments. The analysis also revealed significant data gaps that require further investigation to enhance the understanding and prediction of fuel swelling and pore dynamics. Collecting comprehensive data across different irradiation conditions, burnup levels, and fuel compositions are essential for refining existing models and developing mechanistic models for both binary and ternary metallic fuels, ultimately improving the integration of modeling and experimental approaches in accident testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

RELAP5-3D validation studies based on the High Temperature Test facility

In the spring and summer of 2019, experiments were conducted at the High Temperature Test Facility (HTTF) that form the basis of an upcoming high-temperature gas-cooled reactor (HTGR) thermal hydraulics (T/H) benchmark. HTTF is an integral effects test facility for HTGR T/H modeling validation. This paper presents RELAP5-3D models of two of those experiments: PG-27, a pressurized conduction cooldown (PCC); and PG-29, a depressurized conduction cooldown (DCC). These models used the RELAP5-3D model of HTTF originally developed by Paul Bayless as a starting point. The sensitivity analysis and uncertainty quantification code, RAVEN was used to perform calibration studies for the steady-state portion of PG-27. Here we developed four PG-27 calibrations based on steady-state conditions. These calibrations all used an effective thermal conductivity equal to 36 % of the measured thermal conductivity, but they differed with respect to the frictional pressure drops and radial conduction models. These models all captured the trends in steady-state temperature distributions and transient temperature behavior well. All four calibrations show room for improvement in predicting the transient temperature rise. The smallest error in temperature rise during the transient was a 21 % underprediction, and the largest was a 48 % underprediction. The errors in transient temperature rise are largely a result of a mismatch in power density between the RELAP5-3D model and the experiment due to the location of active heater rods along the boundary between heat structures in the model. The best of these calibrations was applied to PG-29 to model the DCC. Once again, temperatures during the transient were underpredicted but trends in temperature were captured. The RELAP5-3D model captured trends in the data but could not reproduce measured temperatures exactly. This result is not attributed to deficiencies in the experimental data or to RELAP5–3D itself. Rather, this result likely arises due to the some of the assumptions and decisions made when the RELAP5-3D model was first developed, prior to the execution of HTTF experiments. An agreement in prediction of temperature trends but challenges reproducing HTTF temperatures within measurement uncertainty is consistent with previous analyses of HTTF in the literature. Future RELAP5-3D validation activities centered around HTTF may be able to provide greater insight into the code’s capabilities for HTGR modeling with a more finely nodalized model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Thermal Ratcheting Analysis of TEDS Packed-bed Thermocline Energy Storage Tank - Modeling Methodology and Data Validation

This report investigates numerical modeling methods for thermal ratcheting analysis of packed-bed thermal energy storage (TES) tank and discusses the validation results via comparison with experimental data. The experimental data obtained from various design characteristics of packed-bed thermocline tanks, including the Thermal Energy Distribution System (TEDS) TES tank at Idaho National Laboratory (INL), were used to validate thermal and mechanical models developed in this study to evaluate the thermal ratcheting potential. The thermal model was shown to predict the transient thermal propagation through the packed-bed thermocline tanks generally well. However, a larger discrepancy was observed during the comparison with the data from TEDS, presumably due to the uncertainty of boundary conditions given from the experiment. Based on the comparative study between the thermal model predictions and experimental data of various packed-bed thermocline tanks, potential improvements were suggested for the future TEDS experiments for more precise validation study. For mechanical (thermally induced stress) analysis, two different modeling approaches were tested to evaluate hoop stress applied to the packed-bed TES tank wall, which is a major cause of thermal ratcheting process: (i) infinite rigidity model and (ii) Drucker-Prager (DP) model. The ‘model (i)’ is a conservative method with infinite rigidity assumption of granular filler inside a TES tank, whereas the ‘model (ii)’ is a method that takes into account more realistic processes such as thermal expansion of filler and tank wall as well as inter-particle interactions during the cyclic operation of a packed-bed TES tank. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One pilot plant. Then, the effects of various model parameters were discussed to evaluate the thermal ratcheting potential of the TEDS TES tank. The preliminary thermal ratcheting analysis implies that the TEDS TES tank will hold its structural integrity during the normal operation cycles.

25 ENERGY STORAGE↗
DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗