Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transient load analysis”

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 595 records · Page 33

Analysis of a nuclear hybrid energy system using absorption chillers and stratified chilled-water storage with an mPower reactor - 187

The intermittency of renewable energy technology, such as wind and solar power, puts strain on electric grids, often forcing carbon-based and nuclear sources of energy to operate in a load follow mode. Operating nuclear reactors in a load-follow fashion is undesirable due to the associated thermal and mechanical stresses placed on the fuel and reactor components. Various methods of Thermal Energy Storage (TES) can be coupled to nuclear (or renewable) power sources to help absorb grid instabilities caused by daily load demand changes and renewable intermittency. This is known as a Nuclear Hybrid Energy System (NHES). Previous research has found that a stratified chilled-water storage tank can help displace peak cooling loads to off-peak hours. In this study, a dynamic single effect, lithium bromide absorption chiller model was developed as a callable FORTRAN subroutine. Cooling towers, a stratified chilled-water storage tank, and cooling loads characteristic of a large campus were included to form a comprehensive TES model. This TES model was coupled to a high-fidelity mPower sized Small Modular Reactor (SMR) FORTRAN code that provides the transient behavior of primary and secondary-side components. Control algorithms based on excess capacity and conditions in the stratified chilled-water storage tank were used to determine when to divert low-pressure steam from the turbines to the absorption chillers in order to keep reactor power near steady-state. Simulation results depict how the TES system and SMR interact, and demonstrate that stratified chilled-water storage is a viable TES technology for an SMR in a NHES. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Climate Shifts Within Major Agricultural Seasons for +1.5 and +2.0 °C Worlds: HAPPI Projections and AgMIP Modeling Scenarios

This study compares climate changes in major agricultural regions and current agricultural seasons associated with global warming of +1.5 or +2.0 °C above pre-industrial conditions. It describes the generation of climate scenarios for agricultural modeling applications conducted as part of the Agricultural Model Intercomparison and Improvement Project (AgMIP) Coordinated Global and Regional Assessments. Climate scenarios from the Half a degree Additional warming, Projections, Prognosis and Impacts project (HAPPI) are largely consistent with transient scenarios extracted from RCP4.5 simulations of the Coupled Model Intercomparison Project phase 5 (CMIP5). Focusing on food and agricultural systems and top-producing breadbaskets in particular, we distinguish maize, rice, wheat, and soy season changes from global annual mean climate changes. Many agricultural regions warm at a rate that is faster than the global mean surface temperature (including oceans) but slower than the mean land surface temperature, leading to regional warming that exceeds 0.5 °C between the +1.5 and +2.0 °C Worlds. Agricultural growing seasons warm at a pace slightly behind the annual temperature trends in most regions, while precipitation increases slightly ahead of the annual rate. Rice cultivation regions show reduced warming as they are concentrated where monsoon rainfall is projected to intensify, although projections are influenced by Asian aerosol loading in climate mitigation scenarios. Compared to CMIP5, HAPPI slightly underestimates the CO2 concentration that corresponds to the +1.5 °C World but overestimates the CO2 concentration for the +2.0 °C World, which means that HAPPI scenarios may also lead to an overestimate in the beneficial effects of CO2 on crops in the +2.0 °C World. HAPPI enables detailed analysis of the shifting distribution of extreme growing season temperatures and precipitation, highlighting widespread increases in extreme heat seasons and heightened skewness toward hot seasons in the tropics. Shifts in the probability of extreme drought seasons generally tracked median precipitation changes; however, some regions skewed toward drought conditions even where median precipitation changes were small. Together, these findings highlight unique seasonal and agricultural region changes in the +1.5 °C and +2.0 °C worlds for adaptation planning in these climate stabilization targets.

Alexander C. Ruane↗

STRANAL-PMC Version 2.0

Version 2.0 of the Strain Rate Dependent Analysis of Polymer Matrix Composites (STRANAL-PMC) software has been released. A prior version was reported in Analyzing Loads and Strains in Polymer- Matrix Composites (LEW-17227), NASA Tech Briefs, Vol. 26, No. 11 (November 2002), page 36. To recapitulate: Modified versions of constitutive equations of viscoplasticity of metals are used to represent deformation of a polymeric matrix. The equations are applied in a micromechanical approach, proceeding upward from slices of unit cells, through the ply level, to the laminate level. The constitutive equations are integrated in time by a Runge- Kutta technique. To predict the ultimate strength of each composite ply, failure criteria are implemented within the micromechanics. The inputs to STRANAL-PMC are the laminate geometry, properties of the fiber and matrix materials, and applied stress or strain versus time. The outputs are time-dependent stresses and strains at the slice, ply, and laminate levels. The improvements in version 2.0 include more rigorous representation of hydrostatic- stress effects in the matrix, refinement and extension of ply failure models, and capabilities to analyze transverse shear stresses. Version 2.0 can be implemented as a material-model code within transient dynamic finite-element codes.

Goldberg, Robert↗

Laboratory Scale Testing of Ignition Overpressure for Space Vehicle Launch Pad Environments

A scale model of a NASA representative space vehicle is used to develop a refined estimate of the transient pressure loads that are expected to form at the base of the vehicle in the event of a vapor cloud explosion. Flammable vapor clouds are known to form prior to engine startup due to the significant amount of unburned hydrogen that is ejected from the combustion chamber. In the event of a vapor cloud explosion, the vehicle and payload must be able to withstand the resulting overpressure waves. The study comprises an array of pressure sensors located along the base heat shield of the scale model space vehicle as well as the interior wall and throat plug plane of the solid rocket booster. A spark source generator is used to simulate the overpressure wave produced by a vapor cloud explosion while measurements are acquired with and without the effect of a mobile launcher. Time- resolved schlieren images of the simulated vapor cloud explosion reveal the path and impact of both the initial wave and several reflected waves on the various components at the base of the space vehicle. In some instances, the reflected waves superpose to create waves that are higher in amplitude than the initial overpressure wave. A time frequency analysis of the pressure waveforms measured inside the solid rocket booster reveal a ring down tone corresponding to a standing wave that is four times the length of the nozzle.

Tinney, Charles E.↗

Integrated Modeling for the Next Generation Space Telescope "Yardstick" Concept

The so-called NASA "Yardstick" design concept for the Next Generation Space Telescope presents unique challenges for systems-level analysis. Simulations that integrate controls, optics, thermal, and structural models are required to evaluate baseline performance, study design sensitivities, and perform design optimization. An integrated modeling approach was chosen using a combination of commercial off-the-shelf and "in-house" developed codes. The resulting capability provides a foundation for linear and non-linear analysis, using both the time and frequency-domain methods. It readily allows various combinations of design parameters and environmental loads to be evaluated directly in terms of key science-related metrics, in this case the scalar RMS (root mean square) line-of-sight and RMS wavefront errors. This presentation first addresses the development of the component, or discipline, models for the Yardstick design. It will then proceed to present the integration of the component models, using linear-systems approaches, in order to support two of the most critical baseline performance analyses: jitter and thermal-elastic stability of the optical telescope assembly (OTA). The results of the jitter analysis indicate that disturbances from the reaction wheels coupled with the lightly-damped and highly-flexible structure present significant challenges to the baseline line-of-sight control architecture. Vibration isolation will be required to meet jitter error requirements. The results of the thermal-elastic analysis indicate that the mirror segment displacements due to ground-to-orbit cool-down of the telescope are within the expected capture range of the segment rigid-body control actuators. This means we will be able to align and phase the primary mirror. However, the results for the analysis of the thermal transient response following an attitude maneuver (slew) show that this telescope design is not sufficiently stable, passively, to meet the wavefront error requirements. Structural re-design is one possibility; alternatively, active thermal control of the OTA may be considered. The Yardstick integrated models were successfully used to demonstrate the feasibility of two thermal control strategies.

Mosier, Gary E.↗

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids

Increasing penetration levels of inverter based distributed energy resources (DERs) impact the legacy distribution system protection. Inverter based DERs provide approximately 1.2-2 pu fault current. In systems with high penetration of inverter based DERs, it is difficult for over-current based protection schemes to differentiate between normal loading conditions and a fault. Directional, distance, and adaptive forms of protection schemes are also affected by low fault currents. This paper analyzes fault generated traveling wave (TW) based high-frequency signatures in the distribution system. In order to simulate such signatures, frequency-dependent distributed parameter line modeling approach is used in this research work to represent distribution lines and underground cables. Modified IEEE 13-bus medium voltage test system is modeled in electromagnetic transient simulation tool and multiple transient scenarios are simulated in this test system. The results are analyzed to understand the high-frequency signatures that can be used to detect and locate faults under high penetration of DERs.

41 EE - Solar Energy Technologies Office (EE-4S)↗

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids

Increasing penetration levels of inverter based distributed energy resources (DERs) impact the legacy distribution system protection. Inverter based DERs provide approximately 1.2-2 pu fault current. In systems with high penetration of inverter based DERs, it is difficult for over-current based protection schemes to differentiate between normal loading conditions and a fault. Directional, distance, and adaptive forms of protection schemes are also affected by low fault currents. This paper analyzes fault generated traveling wave (TW) based high-frequency signatures in the distribution system. In order to simulate such signatures, frequency-dependent distributed parameter line modeling approach is used in this research work to represent distribution lines and underground cables. Modified IEEE 13-bus medium voltage test system is modeled in electromagnetic transient simulation tool and multiple transient scenarios are simulated in this test system. The results are analyzed to understand the high-frequency signatures that can be used to detect and locate faults under high penetration of DERs.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Control of High-Temperature Static and Transient Thermomechanical Behavior of SiMo Ductile Iron by Al Alloying

Abstract Silicon and molybdenum ( SiMo ) ductile iron is commonly used for exhaust manifolds because these components experience thermal cycling in oxidizing environment, which requires resistance to fatigue during transient thermomechanical loads. Previous studies have demonstrated that alloying elements, such as Al , to SiMo ductile iron reduces the amount of surface degradation during static high-temperature exposure. However, deterioration of sphericity of the graphite nodules and a decrease in ductility could affect the tendency of cracking during thermal cycling. In this article, the effect of Al alloying on static and transient thermomechanical behavior of SiMo ductile iron was investigated to optimize the amount of Al alloying. A thermodynamic approach was used to confirm the effect of the Al alloying on the phase transformations in two SiMo cast irons, alloyed by 1.8% Al and 3% Al . These two alloys were cast in a laboratory along with the baseline SiMo ductile iron. Several experimental methods were used to evaluate the dimensional stability, physical properties, static oxidation, and failure resistance during constrained thermal cycling testing to compare their high-temperature capability. Experimental results verified that Al alloying increases the temperature range and decreases volume change during eutectoid transformation, which together with enhancement of oxidation protection improved the dimensional stability. Thermocycling tests showed that the number of cycles to failure depends on the amount of Al alloying and the applied high-temperature exposure during each cycle. SEM/EDX, high-resolution TEM and µCT analysis were used to verify the mechanism resulting from the Al alloying protection. It was shown that an optimal level of Al alloying for balancing oxidation and thermal cracking resistance depends on thermomechanical conditions of application.

36 MATERIALS SCIENCE↗

Scoping Thermal Response Calculations of RNS Waste During Transport to and Disposal at the WIPP

Sandia National Laboratories (SNL) was contracted by the United States Department of Energy Environmental Management (DOE-EM), Los Alamos Field Office to perform mechanical and thermal scoping calculations as part of a study seeking to understand the ignitability risk of the Remediated Nitrate Salts (RNS) waste drums during transportation from the Waste Control Specialists (WCS) facility to Waste Isolation Pilot Plant (WIPP) and permanent disposal of the waste at WIPP. The scoping thermal simulations described in this report pertain to thermal calculations performed with a packaging system consisting of one Standard Waste Box (SWB) loaded with drums placed inside a Standard Large Box 2 (SLB2). During transportation, the SLB2 is inside Transuranic Package Transporter Model III (TRUPACT-III), which provides the third layer of the packaging. Once at the WIPP, it is assumed the SLB2 is extracted from the TRUPACT-III and maintained above ground, and then subsequently placed underground for permanent disposal. In these proposed configurations, the space between the SLB2 and the SWB is always filled by a layer of insulation consisting of air-filled glass microbubbles except for the bottom which rests directly on the SLB2. The thermal scoping calculations described in this report specifically address whether the introduction of external heat inputs, combined with the contributions from the internally generated radiolytic decay heat and chemical reactions, lead to an unstable thermal state during the time of its movement and placement in the permanent disposal location. The external heat inputs are of two forms: 1) ambient thermal irradiation (e.g., solar and ambient storage/disposal temperatures) and 2) accident-induced fire. Three scoping calculation scenarios were derived as representative, conservative scenarios: 1A) TRUPACT-III transient transportation, 1B) SLB2 48-hour outdoor storage with solar radiation, and 2) fully-engulfing fire during SLB2 handling or emplacement following a steady-state analysis in a 38 °C environment. All the simulated scenarios are conservative relative to the operational conditions expected for handling the waste package during transportation and placement in the WIPP underground disposal unit. The predictions obtained from simulating the three exposure scenarios revealed that adding the SLB2 and the air-filled glass microbubbles to the transport and storage/disposal configurations provides additional thermal protection of the drums beyond what the SWB provides alone, both during long-term above ground insolation and underground during a fire accident. Under the current transportation/storage/disposal concepts, the degree of protection provided by the packaging concept is sufficient to prevent the waste from being ignitable. The simulation results demonstrate that there is adequate margin to safely transport and place the RNS waste from WCS to the WIPP under the current operational concept.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-Burnup BWR LOCA Burst Analysis Using High-Fidelity Multiphysics Simulations

The US nuclear industry is looking to improve on the operating economics of the current fleet of light-water reactors (LWRs). One way of achieving this is by operating fuel to higher burnup. In pressurized water reactors (PWRs), relaxing the current burnup limit will allow for cycle length extensions and power uprates; in boiling water reactors (BWRs) it may allow for improved fuel utilization and reduced feed assemblies, as well as more efficient power uprates and increased capacity factors that will support the Administration’s Executive Order to facilitate 5 GW of power uprates at existing nuclear facilities. However, one of the key limitations to operating fuel to higher burnup is the risk of fuel fragmentation, relocation, and dispersal (FFRD). Recognizing the high interest in extending burnup limits, the US Nuclear Regulatory Commission (NRC) has issued Draft Regulatory Guide DG-1434, which defines an approach that would be acceptable to the NRC for addressing FFRD risk. The approach defined will require better understanding of the phenomena leading to FFRD as well as best-estimate simulation methods to understand FFRD risk in high-burnup cores. The Nuclear Energy Advanced Modeling and Simulation program is supporting the FFRD industry challenge problem through development of state-of-the-art, high-fidelity modeling and simulation LWR analysis capabilities; namely, the BISON fuel performance code and the VERA core simulator software. These tools, along with the US NRC TRACE system analysis code, have been utilized for analysis of FFRD risk in both PWR and BWR cores in recent years. The work documented in this report addresses the lack of high-fidelity research for BWRs and builds on a previous activity where the framework has been applied to Cycles 16 through 18 of Limerick Unit 1, a BWR/4, with introduction of 8 high-burnup lead use assemblies (HBLUAs) that were representative of the 8 HBLUAs loaded into Limerick Unit 2 in 2021. VERA was used in this previous activity to model rod-by-rod depletion in these cycles, and its solution was used to initialize a TRACE simulation of a large-break loss-of-coolant accident (LBLOCA) at the end of Cycle 18. In the work documented in this report, the TRACE model was improved by refining the core mesh and utilizing a new feature that allows for capturing the full 3D VERA power distribution in the model. This allows for a more detailed solution for setting BISON boundary conditions. Furthermore, the solutions from VERA and TRACE were used to set up and perform BISON simulations of about 1,000 rods sampled from the core, including all burnup levels. Utilizing two cladding burst models, it was shown that no fuel rods were predicted to burst during the postulated LBLOCA transient. Additionally, a sensitivity study was performed by artificially increasing linear heat rate during the postulated LBLOCA to identify parameters that correlate with rod burst susceptibility. Burnup, fission gas release, and hoop strain were all found to be positively correlated with rod burst susceptibility. Small-break loss-of-coolant accident (SBLOCA) analyses were also performed; these analyses predicted cladding temperature increases that were bounded by the LBLOCA cladding temperatures for all small break sizes studied for this plant. However, future refinements to the plant response assumptions during the SBLOCA could impact the predicted cladding response. Finally, a benchmark study was performed between CTF and TRACE for LOCA conditions to better qualify CTF for BWR LOCA modeling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

Heat Integration Optimization and Dynamic Modeling Investigation for Advancing the Coal-Direct Chemical Looping Process

The purpose of the project is to address the optimization and startup operation of a modular coal direct chemical looping (CDCL) combustion system integrated with a steam cycle for power generation to reduce the risks involved in further scale-up of the technology. The modular reactor design of the CDCL process provides flexibility in the fabrication of the reactor and in its operating capacity (i.e. turndown ratio) at the cost of a more complex heat exchange network (HEN) design and integration. To address the technology gaps and advance the efficiency and economic feasibility of the CDCL technology, the project will perform a detailed and comprehensive analysis of the integration of a modular CDCL reactor system and a steam cycle system under both static and transient conditions via HEN process performance simulations and system dynamic modeling, respectively. The scope of work consists of 1) Experimental and computational studies of the CDCL combustor reactor 2) Comprehensive static (i.e. steady-state) system HEN design analysis in CDCL 550 MWe commercial unit for power generation and 3) Dynamic modeling of site specific design of 10MWe CDCL large pilot plant. The project team has successfully developed and validated a kinetic model for the oxidation of oxygen carriers in the combustor using the unreacted shrinking core model (UCSM). The model is capable of capturing the oxidation kinetics of fully or partially reduced oxygen carrier particles. A computational fluid dynamics (CFD) model is developed to simulate the hydrodynamics, heat transfer, and chemical reaction occurring in the CDCL combustor. The model is developed in MFIX and ANSYS Fluent. Key aspects of CDCL combustor operation, including heat transfer, oxygen carrier oxidation, and the transport of oxygen carrier particles, are simulated using this CFD model. The HEN for a commercial scale 550 MWe CDCL power plant is simulated and optimized using ASPEN Plus. Practical design considerations are incorporated based on industrial experiences. The performance and cost for the commercial CDCL plant is updated based on these analyses. A dynamic model for the 10 MWe CDCL pilot plant is developed in ProTRAX simulation software. The model is based on the pilot plant design developed in project DE-FE0027654 “10 MWe CDCL Large Pilot Plang – Pre-FEED Study” and the steam cycle data obtained from Dover Light & Power plant. The transient behaviors during pilot plant load variation are simulated using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗

Leveraging the checkpoint-restart technique for optimizing CPU efficiency of ATLAS production applications on opportunistic platforms

Data processing applications of the ATLAS experiment, such as event simulation and reconstruction, spend considerable amount of time in the initialization phase. This phase includes loading a large number of shared libraries, reading detector geometry and condition data from external databases, building a transient representation of the detector geometry and initializing various algorithms and services. In some cases the initialization step can take as long as 10-15 minutes. Such slow initialization has a significant negative impact on overall CPU efficiency of the production job, especially when the job is executed on opportunistic, often short-lived, resources such as commercial clouds or volunteer computing. In order to improve this situation, we can take advantage of the fact that ATLAS runs large numbers of production jobs with similar configuration parameters (e.g. jobs within the same production task). This allows us to checkpoint one job at the end of its configuration step and then use the generated checkpoint image for rapid startup of thousands of production jobs. By applying this technique we can bring the initialization time of a job from tens of minutes down to just a few seconds. In addition to that we can leverage container technology for restarting checkpointed applications on the variety of computing platforms, in particular of platforms different from the one on which the checkpoint image was created. We will describe the mechanism of creating checkpoint images of Geant4 simulation jobs with AthenaMP (the multi-process version of the ATLAS data simulation, reconstruction and analysis framework Athena) and the usage of these images for running ATLAS Simulation production jobs on volunteer computing resources (ATLAS@Home) and on Supercomputers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Characterizing Output Bottlenecks of a Production Supercomputer: Analysis and Implications

This article studies the I/O write behaviors of the Titan supercomputer and its Lustre parallel file stores under production load. The results can inform the design, deployment, and configuration of file systems along with the design of I/O software in the application, operating system, and adaptive I/O libraries.We propose a statistical benchmarking methodology to measure write performance across I/O configurations, hardware settings, and system conditions. Moreover, we introduce two relative measures to quantify the write-performance behaviors of hardware components under production load. In addition to designing experiments and benchmarking on Titan, we verify the experimental results on one real application and one real application I/O kernel, XGC and HACC IO, respectively. These two are representative and widely used to address the typical I/O behaviors of applications.In summary, we find that Titan’s I/O system is variable across the machine at fine time scales. This variability has two major implications. First, stragglers lessen the benefit of coupled I/O parallelism (striping). Peak median output bandwidths are obtained with parallel writes to many independent files, with no striping or write sharing of files across clients (compute nodes). I/O parallelism is most effective when the application—or its I/O libraries—distributes the I/O load so that each target stores files for multiple clients and each client writes files on multiple targets in a balanced way with minimal contention. Second, our results suggest that the potential benefit of dynamic adaptation is limited. In particular, it is not fruitful to attempt to identify “good locations” in the machine or in the file system: component performance is driven by transient load conditions and past performance is not a useful predictor of future performance. For example, we do not observe diurnal load patterns that are predictable.

97 MATHEMATICS AND COMPUTING↗