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 559 records · Page 31

Analysis of the thermal performance of heat pipe radiators

A comprehensive mathematical model and computational methodology are presented to obtain numerical solutions for the transient behavior of a heat pipe radiator in a space environment. The modeling is focused on a typical radiator panel having a long heat pipe at the center and two extended surfaces attached to opposing sides of the heat pipe shell in the condenser section. In the set of governing equations developed for the model, each region of the heat pipe - shell, liquid, and vapor - is thermally lumped to the extent possible, while the fin is lumped only in the direction normal to its surface. Convection is considered to be the only significant heat transfer mode in the vapor, and the evaporation and condensation velocity at the liquid-vapor interface is calculated from kinetic theory. A finite-difference numerical technique is used to predict the transient behavior of the entire radiator in response to changing loads.

Boo, J. H.↗

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Computer program for analysis of high speed, single row, angular contact, spherical roller bearing, SASHBEAN. Volume 2: Mathematical formulation and analysis

Spherical roller bearings have typically been used in applications with speeds limited to about 5000 rpm and loads limited for operation at less than about 0.25 million DN. However, spherical roller bearings are now being designed for high load and high speed applications including aerospace applications. A computer program, SASHBEAN, was developed to provide an analytical tool to design, analyze, and predict the performance of high speed, single row, angular contact (including zero contact angle), spherical roller bearings. The material presented is the mathematical formulation and analytical methods used to develop computer program SASHBEAN. For a given set of operating conditions, the program calculates the bearings ring deflections (axial and radial), roller deflections, contact areas stresses, depth and magnitude of maximum shear stresses, axial thrust, rolling element and cage rotational speeds, lubrication parameters, fatigue lives, and rates of heat generation. Centrifugal forces and gyroscopic moments are fully considered. The program is also capable of performing steady-state and time-transient thermal analyses of the bearing system.

Aggarwal, Arun K.↗

HVDC Macrogrid Modeling for Power-flow and Transient Stability Studies in North American Continental-Level Interconnections

High Voltage Direct Current (HVDC) transmission lines are being constructed throughout the world, aided by advancements in power electronics and the potential value to transfer power between distant areas and off-shore locations. Multiple HVDC lines within and across large AC interconnections could bring about economic benefits such as interregional capacity exchange, and transfer of low-cost, distant electric energy directly to load centers. In addition, network configuration of HVDC lines could result in additional benefits that have not been deeply studied. This paper describes the modeling process for continental-level power system interconnections with the addition of multiple HVDC lines configured as a macrogrid. The models used for study are based on industry-accepted power-flow and dynamic system models for the North American Eastern and Western Interconnections. The model provides insight on feasibility and initial steady-state and stability tests of the HVDC macrogrid and its interactions with the existing electricity infrastructure, opening the door to analysis of the technical value of such a macrogrid.

HVDC transmission, North American Power Grids, AC ↗

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗

Dynamic Finite Element Predictions for Mars Sample Return Cellular Impact Test #4

The nonlinear, transient dynamic finite element code, MSC.Dytran, was used to simulate an impact test of an energy absorbing Earth Entry Vehicle (EEV) that will impact without a parachute. EEVOs are designed to return materials from asteroids, comets, or planets for laboratory analysis on Earth. The EEV concept uses an energy absorbing cellular structure designed to contain and limit the acceleration of space exploration samples during Earth impact. The spherical shaped cellular structure is composed of solid hexagonal and pentagonal foam-filled cells with hybrid graphite-epoxy/Kevlar cell walls. Space samples fit inside a smaller sphere at the center of the EEVOs cellular structure. Pre-test analytical predictions were compared with the test results from a bungee accelerator. The model used to represent the foam and the proper failure criteria for the cell walls were critical in predicting the impact loads of the cellular structure. It was determined that a FOAM1 model for the foam and a 20% failure strain criteria for the cell walls gave an accurate prediction of the acceleration pulse for cellular impact.

Fasanella, Edwin L.↗

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗

Transient response for interaction of two dynamic bodies

During the launch sequence of any space vehicle complicated boundary interactions occur between the vehicle and the launch stand. At the start of the sequence large forces exist between the two; contact is then broken in a short but finite time which depends on the release mechanism. The resulting vehicle response produces loads which are very high and often form the design case. It is known that the treatment of the launch pad as a second dynamic body is significant for an accurate prediction of launch response. A technique was developed for obtaining loads generated by the launch transient with the effect of pad dynamics included. The method solves uncoupled vehicle and pad equations of motion. The use of uncoupled models allows the simulation of vehicle launch in a single computer run. Modal formulation allows a closed-form solution to be written, eliminating any need for a numerical integration algorithm. When the vehicle is on the pad the uncoupled pad and vehicle equations have to be modified to account for the constraints they impose on each other. This necessitates the use of an iterative procedure to converge to a solution, using Lagrange multipliers to apply the required constraints. As the vehicle lifts off the pad the coupling between the vehicle and the pad is eliminated point by point until the vehicle flies free. Results obtained by this method were shown to be in good agreement with observed loads and other analysis methods. The resulting computer program is general, and was used without modification to solve a variety of contact problems.

Prabhakar, A.↗

Gateway Power Quality Lessons Learned

Power Quality is a physical description of the electrical characteristics that allow the system to function properly without significant loss of performance or life. This physical description includes steady state voltage limits, transient voltage limits in normal/abnormal conditions, ripple voltage, stability, fault conditions, and more. All which are vital for improving reliability, ensuring stable operation, defining proper fault recovery, and ensuring a ‘plug and play’ approach to design and integration. Typically, a specification for Power Quality is created based on expected system performance or an existing standard, such as the International Space Power System Interoperability Standards (ISPSIS). One example of such a specification is the Gateway Power Quality Specification. It defines the requirements and characteristics of the 120 Volt direct current electrical power system for the Gateway Electrical Power System (EPS) and the Gateway Electrical Power Consuming Equipment (EPCE). This specification also maintains a separate requirement verification section that defines test methods for requirement verification. The test methods include analysis, test, inspection, and demonstration. These test methodologies and requirements are used to ensure that the loads operate when connected to the specified power quality and performance as defined by this specification. The challenge with developing a specification is that desired system characteristics are not always fully matured before the specification is needed and many performance requirements may be application specific. This drives the need to utilize lessons learned through extensive analysis and testing as well as historical knowledge to finalize requirements. Some example requirements where this is important are Small and Large Signal Stability, Ripple Voltage, Inrush/Surge Currents, and Fault Containment. Lessons learned are also important in the testing, analysis, and verification to ensure consistent and accurate results to verify performance. This presentation will cover the lessons learned for power quality relative to ripple, inrush/surge, fault containment, testing/verification, and more.

Power Quality↗

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Impact-induced twinning in a magnesium alloy under different stress conditions

We investigate deformation twinning in a highly textured magnesium alloy plate of Mg-3.1%Al-0.9%Zn-0.4%Mn under edge-on impact with in situ, synchrotron-based, ultrafast X-ray diffraction measurements, and corresponding stress states are simulated with the finite element method. Deformation twinning and its anisotropy under the triaxial stress condition are explained by a statistical analysis of resolved shear stress. The critical resolved shear stress criterion is applicable under complicated stress conditions induced by high strain rate impact loading. Three typical simple stress conditions are further explored as verification and application cases: uniaxial-stress, uniaxial-strain, and plane-stress. Extension twinning in the magnesium alloy is prone to occur for impact loading applied perpendicular to the crystallographic c-axis, regardless of the exact stress conditions.

36 MATERIALS SCIENCE↗

Numerical Modeling of the Chilldown of Cryogenic Transfer Lines Using a Sinda/GFSSP Integrated Solver

An important first step in cryogenic propellant loading is the chilldown of transfer lines. During the chilldown of the transfer line, the flow is two-phase and unsteady, with solid to fluid heat transfer and therefore a coupled thermo-fluid analysis is necessary to model the system. This paper describes a numerical model of pipe chilldown that utilizes the Sinda/GFSSP Conjugate Integrator (SGCI). SGCI is a new analysis tool developed at NASA's Marshall Space Flight Center (MSFC). SGCI facilitates the solution of thermofluid problems in interconnected solid-fluid systems. The solid component of the system is modeled in MSC Patran and translated into an MSC Sinda thermal network model. The fluid component is modeled in GFSSP, the Generalized Fluid System Simulation Program. GFSSP is a general network flow solver developed at NASA/MSFC. GFSSP uses a finite-volume approach to model fluid systems that can include phase change, multiple species, fluid transients, and heat transfer to simple solid networks. SGCI combines the GFSSP Fortran code with the Sinda input file and compiles the integrated model. Sinda solves for the temperatures of the solid network, while GFSSP simultaneously solves the fluid network for pressure, temperature, and flow rate. The two networks are coupled by convection heat transfer from the solid wall to the cryogenic fluid. The model presented here is based on a series of experiments conducted in 1966 by the National Bureau of Standards (NBS). A vacuum-jacketed, 200 ft copper transfer line was chilled by liquid nitrogen and liquid hydrogen. The predictions of transient temperature profiles and chilldown time of the integrated Sinda/GFSSP model will be compared to the experimental measurements.

LeClair, Andre↗

Modular Controllable Transformers (MCT)

Large Power Transformers (LPTs) are at the heart of the electricity infrastructure and are critical pieces in today’s grid. Aging assets, long turn-around times, transportation delays and dependence on foreign manufacturing for these components have created conditions for points of failure across the system. Moreover, there is no available dynamic control on these passive sections of the grid. For all the automation embedded on the low voltage side of the low voltage side of the system there is absolutely no flexibility on the high voltage side. Changes in power flow patterns are brought about by altering dispatch patterns for generators which requires a coordinated effort and is often time consuming. Traditional FACTS based approaches such as UPFC or HVDC light involve massive modifications to the existing infrastructure. Recent instances of physical attacks on LPTs and the slow recovery process associated with fixing or replacing these damaged units has created an increased awareness for this issue. LPT designs are highly customized owing to specific field requirements. Long lead times coupled with transportation delays make procuring and installing LPTs a lengthy process. In an attempt to minimize the recovery time following an LPT failure, transformer manufacturers have developed mobile transformer units. These units typically consist of either a single-phase or three-phase transformer mounted on a truck and characterized by compact form factors. By reducing the assembly time and transport time these units offer a fast and temporary solution to LPT failures. Several efforts led by DOE, Edison Electric Company (EEI) and North American Electric Reliability Corporation (NERC) have led to the RecX program. RecX transformers consist of single phase transformers transported on specialized assemblies to enable fast restoration and minimal installation time. NERC and EEI have instituted transformer sharing programs like Spare Equipment Database (SED) as well as Spare Transformer Equipment Database (STEP) which encourage utilities to maintain and share a fleet of LPTs. Mobile transformers are typically limited to less than 100 MVA to enable transportation on a trailer. Also, the efficiency of mobile transformer units and RecX transformers is lower since they are temporary solutions. The ideal replacement solution would be one that could utilize a set of standardized transformers while incorporating flexibility to influence metrics like power flow, apparent impedance and voltage profiles. In view of these issues the team at Georgia Tech proposed a concept called Modular Controllable Transformers (MCTs). This project was commissioned as a part of DOE’s Transformer Resiliency and Advance Component (TRAC) program. The approach consists of splitting one large unit rated upwards of 100 MVA with smaller modular, standardized units. Each of these smaller units are augmented with a converter to add power flow, voltage and impedance control. Splitting the power level into multiple units has two important implications. By making smaller units, the transportation times can be reduced drastically. Moreover, even though these units have very specific impedance and other parameters, the augmented converters can make it emulate the unit being replaced. Thus, faster replacement times and added controls are achieved making this an extremely resilient approach. In the event of one modular units failure, the other two controllable units can operate in the most resilient fashion keeping the net system overloads to a minimum. Moreover, for the same probability of failure only a fraction of the system capacity is lost. The converter is augmented on to the tertiary winding of the transformer making it retrofittable in the field without massive modifications to the system. The material proposed in this report highlights design considerations to minimize recovery time as well as results displaying overall improvement in resiliency. The cost analysis of this approach shows significant savings compared to UPFC and HVDC lite based approaches. Extensive analysis was conducted on numerous aspects of this project over the course of this project. Simulation studies were conducted to prove the 5MVA converter design. Challenging BIL constraints and transient issues were addressed over the course of this project. Simulations on EMTP platform presented in consequent sections were used to finely tune the control for the same. A final design for a back-to-back 5MVA neutral-point-clamped (NPC) converter was completed in this project. Furthermore, switching components like high power IGBTs from IXYS Westcode Press-Pack series and ABB Stackpack series were identified for this application. These devices exhibit a blocking voltage of 6.5 kV and a nominal collector current of 900 A. Potential housing options as well as existing commercial option were identified for this converter design. A 56 MVA modular transformer unit was designed with Delta Star Inc. with consultation and inputs from all project members. The transformer cooling systems were chosen to be Oil Natural Air Forced (ONAF) to optimize ease of transport. A full Product Requirement Document (PRD) for the converter and the transformer was generated. The team also investigated key aspects of resiliency and how to quantify the same. A metric was developed to analyze resiliency and numerous simulation cases were presented. Increased (N-1) resiliency and outage resistance was showcased through simulation on the Texas Large scale power system. Dynamic control capabilities were also showcased on numerous system level simulations. An IEEE 30 bus system was used to show the ability of the MCT approach to reduce line congestion. Power flow control; a part of the MCTs abilities, showcased a 29% reduction in transformer loading on the same system. Full control over tie line flows was demonstrated on a modified IEEE 13 bus system. Dynamic LTC capabilities were shown through simulation and were compared to traditional LTCs. Furthermore, the ability of the MCTs to alter apparent transformer impedance was simulated to promote standardization of modular transformer units while retaining the ability to emulate a varied range of impedance characteristics. Thus granular control over power flow, voltages, impedances and was presented. Regulatory analysis conducted by Oak Ridge National Laboratory (ORNL) further bolstered the efficacy of the MCT approach and verified the added resiliency achieved due to this approach. As a conclusion, it was proved that the resilience improvement obtained from using the MCT approach is justified by just one failure in a 30-year period. This means that the benefits obtained from using the MCT approach far outweigh the cost of setting this approach up. As a part of the future efforts, the team is proposing building a (12.47-24) kV/400 kVA scaled down unit to control 5MVA of power flows. Emphasis would be laid on using commercial components to speed up the build of these units. Further the actual transport of the modular transformer unit would be a part of this operation and would provide insight into transportation issues. The approach shows much promise and adds key dynamic control on passive sub transmission and transmission sectors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

Weight function procedure for reduced order fracture analysis of arbitrary flaws in cylindrical pressure vessels

Fracture mechanics calculations using the finite element method can be computationally expensive, which makes them challenging to use in engineering evaluations of flaws in pressure vessels. The weight function (WF) technique is a reduced order fracture modeling method that is widely used to greatly reduce these computational costs. Its computational efficiency is essential for use in probabilistic fracture mechanics evaluations of embrittled nuclear reactor pressure vessels (RPVs), due to the large number of sampled flaws that must be evaluated. Although the WF technique is general, it is typically used only for axis-aligned flaws. Recent discoveries of off-axis flaws in operating nuclear reactors have necessitated detailed simulations of such flaws and the interactions between neighboring flaws. This paper presents a generalized WF approach applicable for analyzing arbitrary flaw geometries in thick-walled cylindrical vessels, including surface-breaking and subsurface flaws, which can either axis-aligned or off-axis, and can account for interactions with other flaws. Herein, this approach is demonstrated on representative simulations of multiple flaw geometries in a RPV subjected to transient loading conditions. In all cases, the WF technique gives good comparison with benchmark results from direct simulation with greatly reduced computational effort.

42 ENGINEERING↗

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Establishing a neutronics design and equilibrium cycle analysis for the I 2 S-LWR reactor with UO 2 and U 3 Si 2 Fuel

The Integral Inherently Safe (I 2 S-LWR) reactor was developed not only to enhance passive reactor safety performance, but also to improve reactor economics. This paper describes the steady-state and depletion analysis of the I 2 S-LWR core using both UO 2 and U 3 Si 2 fuel designs. The development and validation of the core models are described, to include the generation of HELIOS cross sections and comparison with the WEC lattice code PARAGON and the Monte Carlo code SERPENT. The few group homogenized cross sections were then prepared for the U.S. NRC PARCS core simulator using the GenPMAXS code. The multicycle capability in PARCS was used to determine the equilibrium core condition for three batch cores loaded with the UO 2 and then with the U 3 Si 2 fuel designs. Comparing the results from the WEC code ANC at beginning, middle, and end of cycle showed good agreement. Core conditions at beginning and end of cycle were then used to determine the converged steady-state conditions of a coupled RELAP5/PARCS model, which provided the initial condition for various transient accident simulations which are reported in a companion paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗