Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load uncertainty”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Ultrasonic Studies of Composites Undergoing Thermal and Fatigue Loading

New composite materials possess attractive properties for use in advanced aircraft. A necessary requirement for their introduction into aeronautic use is an accurate understanding of their long term aging processes so that proper design criteria can be established. In order to understand those properties, these composites must be exposed to thermal and load cycles that are characteristic of flight conditions. Additionally, airline companies will require nondestructive evaluation (NDE) methods that can be used in the field to assess the condition of these new materials as they age. As part of an effort to obtain the required information about new composites for aviation use, we are performing ultrasonic measurements both in the NDE laboratory and in the materials testing laboratory at NASA. The materials testing laboratory is equipped with environmental chambers mounted on load frames so that composite samples can be exposed to thermal and loading cycles representative of flight protocols. Applying both temperature and load simultaneously will help to highlight temperature and load interactions during the aging of these composite materials. This study reports on our initial ultrasonic attenuation results from thermoset and thermoplastic composite samples. Ultrasonic attenuation measurements have been used reliably to assess the effects of material degradation. For example, recently, researchers have shown that by using frequencies of ultrasound on the order of 24 MHz, they could obtain adequate contrast in the evaluation of thermal degradation in these composites. This paper will present data that shows results at a lower frequency range. In addition, we report results on the frequency dependence of attenuation as the slope of attenuation with respect to frequency, beta = delta alpha (f) / delta f. The slope of attenuation is an attractive parameter since it is quantitative, yet does not require interface corrections like conventional quantitative attenuation measurements. This is a consequence of the assumption that interface correction terms are frequently independent. Uncertainty in those corrections terms compromises the value of conventional quantitative attenuation data.

Madaras, Eric I.↗
DOE OSTI2022

Effect of Thermal Expansion Coefficient, Viscosity and Melting Range in Simulation of PCM Embedded Heat Exchangers With and Without Fins

Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.

CFD↗
DOE OSTI2020

ZPPR-15 Small Reactivity Worth Experiments

In addition to measuring the worth of sodium voiding and the worth of simulated control rods, the ZPPR staff measured the Doppler worth of four standard samples, the worth of axial expansion and the worth of radial bowing in ZPPR-15. Doppler sample worths were measured in ZPPR-15A, ZPPR-15B and ZPPR-15D. The Doppler samples were natural UO 2 , depleted uranium metal, depleted U-10Zr alloy and 33% enriched UO 2 . These samples were placed in a sample capsule in the inner core of ZPPR-15 and heated to various temperatures between 300 K and 1100 K. The resulting changes in reactivity relative to the reference configuration were measured to determine the Doppler worths of the samples. Special segmented drawers were constructed to measure the worth of axial expansion in ZPPR15. The segments were connected to each other, and the back segment was connected to a spring-loaded cable. Pulling the cable introduced small gaps between drawer segments to simulate axial expansion, and releasing the cable eliminated those gaps. The worth of this simulated axial expansion was measured in ZPPR-15A, ZPPR-15B and ZPPR-15D. The worth of radial bowing was measured in ZPPR-15, but the experimenters regarded the result as unsatisfactory. Radial bowing was measured by two methods in ZPPR-17A, so the ZPPR-17A bowing measurements were analyzed in place of the unpublished ZPPR-15 radial bowing measurement. The first ZPPR-17A measurement consisted of rearranging the plates in 96 core drawers at the core-radial blanket boundary to simulate outward motion of the fuel during bowing. The second ZPPR-17A measurement used a new bowing oscillator to move the core plates a small distance vertically in a special drawer. The published uncertainties for the Doppler worth, axial expansion worth and radial bowing worth measurements are the statistical uncertainties in the measurements. In reality, the published uncertainty for each of these measurements is just one component of the total uncertainty. There are additional uncertainties specific to each measurement and a common uncertainty related to conversion from the natural measurement units, cents, to pcm for calculations. A full uncertainty analysis was performed for each measurement. The significant uncertainties were quantified, and a total uncertainty was determined for each measurement. The ZPPR-15 and ZPPR-17A experimental records were used to create detailed as-built Monte Carlo models of the Doppler worth, axial expansion and radial bowing measurement configurations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Regulatory Testing of RPP-WTP HLW Glasses to Support Delisting Compliance, VSL-04R4780-1, Rev. 0 (Sep 2004)

The primary goal of the testing described in this report was to collect data to demonstrate compliance of the immobilized high-level waste (IHLW) glasses with delisting requirements. The collected data will be used to support a petition to delist the IHLW glasses destined for the national disposal facility. The Delisting Data Quality Objectives (DQO) (Cook and Blumenkranz 2003) identified a list of (16) inorganic constituents of potential concern (COPCs) and their associated limits for delisting. These COPCs can be divided into three groups, Cases 1, 2, and 3, based on their Toxicity Characteristic Leaching Procedure (TCLP) responses versus their respective delisting limits. To briefly summarize, Case 1 COPCs are those that, when loaded at their highest expected concentration in Waste Treatment Plant (WTP) glasses, are not expected to leach at their respective delisting limits when the glasses are exposed to the TCLP. Case 2 COPCs may reach the delisting limits in TCLP leachates of WTP glasses if loaded to concentrations near their maximum expected concentrations in glass. Finally, Case 3 COPCs are components that are likely to be present in concentrations sufficient to exceed their respective delisting limits in TCLP leachates of some possible glasses. The test objective was to show that, for (i) the expected range of inorganic contents in the waste feed to the Hanford WTP high level waste (HLW) vitrification facility, (ii) the expected range of glass product compositions, and (iii) the Case 1 and Case 2 COPCs identified by the DQO, the IHLW glasses meet all the relevant requirements for delisting. For Case 3 COPCs (i.e., Cd), the testing was to demonstrate the relationship between glass composition and TCLP cadmium (Cd) release, and then employ the results to develop TCLP-composition response models. During WTP operations, TCLP-composition models can be used to predict, within the required statistical uncertainties, TCLP responses of IHLW production glasses that are within the compositional region used to develop the model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Statistical Correlation of Heliostat Pointing Deviation With Wind

This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.

17 WIND ENERGY↗
DOE OSTI2023

Off-policy deep reinforcement learning with automatic entropy adjustment for adaptive online grid emergency control

Electric overloading conditions and contingencies put modern power systems at risk of voltage collapse and blackouts. Load shedding is crucial to maintain voltage stability for grid emergency control. However, the rule- or model-based schemes rely on accurate dynamic system models and face considerable challenges in adapting to various operating conditions and uncertain event occurrences. Here, to address these issues, this paper proposes a novel deep reinforcement learning (DRL)-based voltage stability control algorithm with automatic entropy adjustment (AEA) for grid emergency control. Various dynamic network components for complex system operations are modeled to construct the DRL environment. An off-policy soft actor-critic architecture is developed to maximize the expected reward and policy entropy simultaneously. The AEA mechanism is proposed to facilitate the policy maximum entropy procedure, and the proposed method can automatically provide effective discrete and continuous actions against various fault scenarios. Our approach accomplishes high sampling efficiency, scalability, and auto-adaptivity of the control policies under high uncertainties. Comparative studies with the existing DRL-based control methods in IEEE benchmarks indicate salient performance improvement of the proposed method for dynamic system emergency control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2008

Assessment of the LC-2 Prelaunch Fatigue Spectra of the CM-to-SM Flange Weld

The pad stay and rollout components of the Ares I-X life cycle can generate cyclic stress oscillations to the vehicle that could initiate and grow fatigue cracks from weld defects. The Ares I-X Project requested that a study be performed to determine if stabilization of the vehicle is required to reduce the stresses that could initiate and grow fatigue cracks at the flange-to-skin weld of the Command Module (CM) and Service Module (SM) interface. A fatigue crack growth analysis was conducted that used loads (LC-2) and stress analyses developed by the Ares I-X Project and utilized material data and analysis methods developed by a critical initial flaw size (CIFS) analysis conducted by NASA Engineering and Safety Center (NESC) for the Upper Stage Simulator (USS) of the Ares I-X vehicle. A full CIFS analysis for the CM-to-SM flange-to-skin weld was not performed because the full flight spectrum was not provided and was not necessary to answer the question posed by the Ares I-X Project. Instead, an approach was developed to determine if the crack growth due to the pad stay and rollout components of the flight spectrum would adversely influence the CIFS. The approach taken used a number of conservative assumptions that eliminated the need for high-fidelity analyses and additional material testing, but still provided a bounding solution for the uncertainties of the problem. The results from this analysis indicate that the LC-2 pad stay and rollout spectrum components would not produce significant fatigue crack growth on the CM-to-SM flange-to-skin weld. Thus, from a fatigue crack growth standpoint, no stabilization is required to reduce the LC-2 pad stay and rollout cyclic stresses on the CM-to-SM flange-to-skin weld.

Dawicke, David S.↗
NASA NTRS2011

Measurement of CPAS Main Parachute Rate of Descent

The Crew Exploration Vehicle Parachute Assembly System (CPAS) is being designed to land the Orion Crew Module (CM) at a safe rate of descent at splashdown. Flight test performance must be measured to a high degree of accuracy to ensure this requirement is met with the most efficient design possible. Although the design includes three CPAS Main parachutes, the requirement is that the system must not exceed 33 ft/s under two Main parachutes, should one of the Main parachutes fail. Therefore, several tests were conducted with clusters of two Mains. All of the steady-state rate of descent data are normalized to standard sea level conditions and checked against the limit. As the Orion design gains weight, the system is approaching this limit to within measurement precision. Parachute "breathing," cluster interactions, and atmospheric anomalies can cause the rate of descent to vary widely and lead to challenges in characterizing parachute terminal performance. An early test had contradictory rate of descent results from optical trajectory and Differential Global Positioning Systems (DGPS). A thorough analysis of the data sources and error propagation was conducted to determine the uncertainty in the trajectory. It was discovered that the Time Space Position Information (TSPI) from the optical tracking provided accurate position data. However, the velocity from TPSI must be computed via numerical differentiation, which is prone to large error. DGPS obtains position through pseudo-range calculations from multiple satellites and velocity through Doppler shift of the carrier frequency. Because the velocity from DGPS is a direct measurement, it is more accurate than TSPI velocity. To remedy the situation, a commercial off-the-shelf product that combines GPS and an Inertial Measurement Unit (IMU) was purchased to significantly improve rate of descent measurements. This had the added benefit of solving GPS dropouts during aircraft extraction. Statistical probability distributions for CPAS Main parachute rate of descent and drag coefficient were computed and plotted. Using test data, a terminal rate of descent at splashdown can be estimated as a function of canopy loading.

Ray, Eric S.↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

High-Fidelity Measurements for Flattop-HEU Benchmark Reevaluation

Flattop was first built in the 1950’s at Los Alamos National Laboratory. Flattop-HEU is composed of a sphere of highly enriched uranium (HEU) surrounded by a thick spherical natural uranium (NU) reflector. The reflector is composed of three parts: a stationary hemisphere and two movable quarter spheres. For fine control of the reactivity of the system, there are three control rods of natural uranium located in voids in the stationary hemisphere. The final components that make Flattop a useful critical assembly are the glory hole and mass adjustment pieces. These pieces can be loaded in various configurations into the glory hole and the core pedestal to control the known worth of the system. The glory hole and mass adjustment pieces are mostly small pieces of HEU with some mass adjustment pieces fabricated from NU. This allows for the irradiation of samples to a specified level. To better document the system, Flattop was evaluated and included in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook. The original benchmark evaluation of Flattop-HEU was written in 1999 based on an experiment completed in the 1960’s. This original evaluation was written to provide a single diameter that defined critical mass; however, as computational capabilities have increased, the focus for benchmark evaluations has shifted to include detailed modelswith all physical dimensions. Thus, as Flattop is a lynchpin in critical experiment work, the benchmark is being reevaluated at current standards. This summary discusses some of the largest known uncertainties from the evaluation and the high-fidelity measurements taken to reduce these uncertainties.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Voltage positioning using co-optimization of controllable grid assets in radial networks

With increasing penetration of solar PV, some distribution feeders are experiencing highly variable net-load flows and even reverse flows. To optimize distribution systems under such conditions, the scheduling of mechanical devices, such as OLTCs and capacitor banks, needs to take into account forecasted solar PV and actual grid conditions. However, these legacy switching assets are operated on a daily or hourly timescale, due to the wear and tear associated with mechanical switching, which makes them unsuitable for real-time control. Therefore, there is a natural timescale- separation between these slower mechanical assets and the responsive nature of inverter-based resources. In this paper, we present a network admissible convex formulation for holistically scheduling controllable grid assets to position voltage optimally against solar PV. An optimal hourly schedule is presented that utilizes mechanical resources to position the predicted voltages close to nominal values, while minimizing the use of inverter-based resources (i.e., DERs), making them available for control at a faster time-scale (after the uncertainty reveals itself). A convex, inner approximation of the OPF problem is adapted to a mixed-integer linear program that minimizes voltage deviations from nominal (i.e., maximizes voltage margins). Here, the resulting OPF solution respects all the network constraints and is, hence, robust against modeling simplifications. Simulation based analysis on IEEE distribution feeders validates the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2000

Novel High Gas-Temperature Calibration System Demonstrated

Accurate measurement of high gas temperatures, typically above 1300 K, has always presented challenges to researchers. Thermocouples often perturb the local gas flow and temperature field; they provide indirect information; and at high temperatures, they require large corrections so that actual gas temperatures can be determined. The physical and chemical stability of thermocouples to withstand the thermal loads and reactive environments prevailing at high temperatures naturally limits their life and maximum use temperature. Optical systems have their own drawbacks since accurate results depend on well-characterized emissivity, optical thickness, and gas composition information. These properties are rarely well known, especially in high-temperature, chemically reacting environments. In addition, optical systems usually require independent calibrations, which often involve the use of thermocouples, and hence, suffer from their aforementioned limitations. A new technique developed by researchers at the NASA Glenn Research Center at Lewis Field exploits an abrupt increase in the emittance of optically thin materials at their unique melting temperatures for a direct determination of gas temperature. Pure metallic-oxide fibers, varying in diameter from 60 to 400 mm, have been used in measurements over a temperature range of 2050 to 2700 K. The accuracy and reproducibility of the technique is estimated to be 15 K: that is, within the uncertainty in the melting points of the materials. Other fiber materials with different, but unique, melting points could be used to extend the technique over a larger temperature range.

Gokoglu, Suleyman A.↗
NASA NTRS2015

Mars Ascent Vehicle-First Stage Motor

This project is development effort of a first stage solid motor based on a two-stage solid motor Mars Ascent Vehicle (MAV) design for the robotic Mars Sample Return (MSR) mission (fig. 1). The MSR MAV has been studied for decades and multiple concepts have been shown to meet the mission objectives as posed.1 However, there remains significant uncertainty with the MAV requirements. The sample container and sample cache itself is immature. Additionally, MAV-specific requirements ranging from full three-axis controlled and strict communication requirements to minimal capability concepts are still under consideration. Given the maturity of the overall mission requirements, the MAV has been limited to a large number of parametric analyses and paper studies. Recently, a Jet Propulsion Laboratory study highlighted the flexibility of a two-stage solid motor concept. The MAV itself is driven by the constraints of the Entry, Decent, and Landing (EDL) system. Within the EDL constraints, there is a range of MAV options ranging in complexity from simple spun upper stage options to higher capability three-axis controlled solutions. There are also options to trade the ratio of mission (Delta)V between the first and second stage. Finally, sensitivity studies also indicated that solid motors with a high percentage of off-load flexability only had minor impact on the total system mass over a single point design optimized motor. This flexibility in the first stage motor has allowed NASA to mature the design of the motor beyond parametric analyses and start to address known design challenges of the motor.

Dankanich, John↗
DOE OSTI2024

Probing the atmospheric boundary layer with integrated remote-sensing platforms during the American WAKE ExperimeNt (AWAKEN) campaign

The American WAKE ExperimeNt (AWAKEN) collaboration is an observational-based field campaign in northern Oklahoma intended to analyze the potential influence of onshore wind farms and their collective wakes on wind power production, turbine structural loads, and on the atmospheric boundary layer (ABL). Focusing on the ABL effects, the University of Oklahoma and the Lawrence Livermore National Laboratory collected continuous high-resolution kinematic and thermodynamic profile measurements during 2022 and Summer 2023. The deployment strategy for these campaigns is detailed first, followed by an initial comparison of data from two sites in the AWAKEN domain: a near-farm site to examine collective wake impacts on the ABL, and a far-field site remaining outside the wind farm-waked region. Here, we summarize the datasets available and demonstrate the benefits of these observations and multiple value-added products (VAPs) for investigation of ABL features observed during AWAKEN. We also highlight examples of preliminary analyses, including ABL height detection and nocturnal low-level jet examination, which are produced using novel VAPs based on optimal estimation to retrieve deeper Doppler lidar wind profiles than previously resolved, along with their uncertainty. By including the near-farm and far-field site in these analyses, we identified a pattern of stronger lower-atmospheric mixing at the near-farm site than the far-field site, motivating deeper investigation into the relationship between wind farms and general ABL characteristics. Future analysis will delve deeper into this relationship by examining other ABL characteristics, such as atmospheric stability and convection.

17 WIND ENERGY↗
DOE OSTI2025

Sensitivity-based voltage constraints for optimal power flow in low-voltage distribution feeders

The optimal power flow (OPF) problem for distribution systems can include network details down to the low-voltage (LV) points of interconnection of individual customers. This paper addresses the implementation of voltage magnitude constraints, and sets forth a practicable approach for capturing the effects on voltage from the switching behavior of loads (e.g., heat pumps, air conditioners, water heaters, or pool pumps) and from the variability of renewable generation (e.g., rooftop solar). The proposed method adjusts the OPF voltage constraints based on forecasts of load and generation upper and lower bounds, in conjunction with sensitivity factors derived from the power flow equations. An illustrative OPF formulation is also provided, which incorporates transformer models that include core loss. We demonstrate that accurate modeling of these LV network components is critical to avoid voltage violations at customer points of interconnection. Furthermore, the ideas are validated through numerical case studies on a realistic distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

MECHANISTIC MULTIPHYSICS MODELING OF CLADDING RUPTURE IN NUCLEAR FUEL RODS DURING LOSS-OF-COOLANT ACCIDENT CONDITIONS

The Loss of Coolant Accident (LOCA) is a design basis accident that is included as part of the safety analysis of nuclear power plants. As the nuclear industry desires to increase the cycle length and discharge burnup of existing nuclear power plants they must demonstrate safe operation during a LOCA on high burnup fuel. During a LOCA transient on high burnup fuel rods, the rods may undergo a process known as fuel fragmentation, relocation, and dispersal (FFRD). To permit dispersal, the cladding encapsulating the fuel must rupture with an opening size large enough to allow the fragmented fuel particles to release. Current licensing tools used by industry and the United States Nuclear Regulatory Commission are limited in geometric fidelity and materials that can be analyzed. These simulation tools generally employ a quasi-two-dimensional (1.5D or Layered1D) or 2D-RZ axisymmetric geometric representations exclusively. While a valid approach under some instances, there are times when important physics have an asymmetric behavior in the fuel rod. Examples include fuel fragmentation, thermal-hydraulic boundary conditions, and cladding rupture, all of which are important for LOCA analysis. As industry pursues burnup extensions it must be demonstrated that fuel dispersal can be mitigated or eliminated. To do this, an understanding of the rupture opening after cladding failure is required. This work presents the development of a model for predicting the size and location of the rupture opening in failed fuel rods during LOCA conditions using advanced modeling and simulation tools. In order to supply the rupture model with appropriate boundary conditions, improvements to fuel fragmentation, axial relocation and oxidation modeling were required. First, the eXtended Finite Element Method (XFEM) is used to mechanistically predict the number of fuel fragments that form due to material strength randomization, criteria for strength randomization, mesh density, power ramping rates and irradiation effects. These predictions with associated uncertainty were compared to empirical correlations developed for UO2 verifying that they can be used with increased confidence in subsequent axial relocation analyses. Secondly, a new first-of-its-kind Layered2D computational framework was developed that provides the ability to apply azimuthally varying boundary conditions while providing discrete layers to track fuel movement during the LOCA. An existing fuel axial relocation model developed for Layered1D was extended to work within the Layered2D framework. A large sensitivity study was performed on the initial version of the model to identify modeling parameters of particular importance, with the emissivity used for radiation after blowdown being the primary source of uncertainty. Then, a simplistic approach to incorporate mechanical degradation of the cladding due to oxidation was also developed to investigate the impact of reduced cladding thickness on predictions of the time to failure of cladding tubes. It was found that the cladding will typically rupture prior to a reduction in thickness that is sufficient to impact the rupture behavior. A model was then developed for predicting cladding rupture that transfers the cladding surface temperatures, rod internal and external pressures, fast neutron fluence, and fast neutron flux from a more detailed Layered1D, Layered2D, or 2D-RZ analysis to a 3D cladding only analysis. Comparisons of the rupture model to a few experiments indicate reasonable predictions. The rupture model was then applied to two accident tolerant fuel concepts (FeCrAl and Cr-coated Zircaloy) where it predicted that both ATF concepts would have smaller rupture openings and delayed rupture times than the standard Zircaloy-4 cladding material under identical loading conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Three-dimensional microstructure-explicit and void-explicit mesoscale simulations of detonation of HMX at millimeter sample size scale

Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.

Miller, Christopher (ORCID:0000000206844381)↗
NASA NTRS2011

Ares I-X Flight Test Validation of Control Design Tools in the Frequency-Domain

A major motivation of the Ares I-X flight test program was to Design for Data, in order to maximize the usefulness of the data recorded in support of Ares I modeling and validation of design and analysis tools. The Design for Data effort was intended to enable good post-flight characterizations of the flight control system, the vehicle structural dynamics, and also the aerodynamic characteristics of the vehicle. To extract the necessary data from the system during flight, a set of small predetermined Programmed Test Inputs (PTIs) was injected directly into the TVC signal. These PTIs were designed to excite the necessary vehicle dynamics while exhibiting a minimal impact on loads. The method is similar to common approaches in aircraft flight test programs, but with unique launch vehicle challenges due to rapidly changing states, short duration of flight, a tight flight envelope, and an inability to repeat any test. This paper documents the validation effort of the stability analysis tools to the flight data which was performed by comparing the post-flight calculated frequency response of the vehicle to the frequency response calculated by the stability analysis tools used to design and analyze the preflight models during the control design effort. The comparison between flight day frequency response and stability tool analysis for flight of the simulated vehicle shows good agreement and provides a high level of confidence in the stability analysis tools for use in any future program. This is true for both a nominal model as well as for dispersed analysis, which shows that the flight day frequency response is enveloped by the vehicle s preflight uncertainty models.

Johnson, Matthew↗