Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Dynamic Programming”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Mixed Integer Programming and Heuristic Scheduling for Space Communication

Optimal planning and scheduling for a communication network was created where the nodes within the network are communicating at the highest possible rates while meeting the mission requirements and operational constraints. The planning and scheduling problem was formulated in the framework of Mixed Integer Programming (MIP) to introduce a special penalty function to convert the MIP problem into a continuous optimization problem, and to solve the constrained optimization problem using heuristic optimization. The communication network consists of space and ground assets with the link dynamics between any two assets varying with respect to time, distance, and telecom configurations. One asset could be communicating with another at very high data rates at one time, and at other times, communication is impossible, as the asset could be inaccessible from the network due to planetary occultation. Based on the network's geometric dynamics and link capabilities, the start time, end time, and link configuration of each view period are selected to maximize the communication efficiency within the network. Mathematical formulations for the constrained mixed integer optimization problem were derived, and efficient analytical and numerical techniques were developed to find the optimal solution. By setting up the problem using MIP, the search space for the optimization problem is reduced significantly, thereby speeding up the solution process. The ratio of the dimension of the traditional method over the proposed formulation is approximately an order N (single) to 2*N (arraying), where N is the number of receiving antennas of a node. By introducing a special penalty function, the MIP problem with non-differentiable cost function and nonlinear constraints can be converted into a continuous variable problem, whose solution is possible.

Lee, Charles H.↗
NASA NTRS2013

NASA Tech Briefs, July 2013

Dielectrophoresis-Based Particle Sensor Using Nanoelectrode Arrays; Multi-Dimensional Damage Detection for Surfaces and Structures; ULTRA: Underwater Localization for Transit and Reconnaissance Autonomy; Autonomous Cryogenic Leak Detector for Improving Launch Site Operations; Submillimeter Planetary Atmospheric Chemistry Exploration Sounder; Method for Reduction of Silver Biocide Plating on Metal Surfaces; Silicon Micromachined Microlens Array for THz Antennas; Forward-Looking IED Detector Ground Penetrating Radar; Fully Printed, Flexible, Phased Array Antenna for Lunar Surface Communication, Battery Charge Equalizer with Transformer Array; An Efficient, Highly Flexible Multi-Channel Digital Downconverter Architecture; Dimmable Electronic Ballast for a Gas Discharge Lamp; Conductive Carbon Nanotube Inks for Use with Desktop Inkjet Printing Technology; Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates; High-Performance, Low-Temperature-Operating, Long-Lifetime Aerospace Lubricants; Carbon Nanotube Microarrays Grown on Nanoflake Substrates; Differential Muon Tomography to Continuously Monitor Changes in the Composition of Subsurface Fluids; Microgravity Drill and Anchor System; 20 Granular Media-Based Tunable Passive Vibration Suppressor; 21 Miga Aero Actuator and 2D Machined Mechanical Binary Latch; Micro-XRF for In Situ Geological Exploration of Other Planets; Hydrogen-Enhanced Lunar Oxygen Extraction and Storage Using Only Solar Power; Uplift of Ionospheric Oxygen Ions During Extreme Magnetic Storms; Miniaturized, High-Speed, Modulated X-Ray Source; Hollow-Fiber Spacesuit Water Membrane Evaporator 25 High-Power Single-Mode 2.65-micrometers InGaAsSb/AlInGaAsSb Diode Lasers; Optical Device for Converting a Laser Beam Into Two Co-aligned but Oppositely Directed Beams; A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging; X-Ray Diffractive Optics; SynGenics Optimization System (SynOptSys); 29 CFD Script for Rapid TPS Damage Assessment; radEq Add-On Module for CFD Solver Loci-CHEM; Science Opportunity Analyzer (SOA) Version 8; 30 Autonomous Byte Stream Randomizer; Distributed Engine Control Empirical/Analytical Verification Tools; Dynamic Server-Based KML Code Generator Method for Level-of-Detail Traversal of Geospatial Data; Automated Planning of Science Products Based on Nadir Overflights and Alerts for Onboard and Ground Processing; Linked Autonomous Interplanetary Satellite Orbit Navigation; Risk-Constrained Dynamic Programming for Optimal Mars Entry, Descent, and Landing; Scheduling Operations for Massive Heterogeneous Clusters; Deepak Condenser Model (DeCoM); Flight Software Math Library; Recirculating 1-K-Pot for Pulse-Tube Cryostats; 35 Method for Processing Lunar Regolith Using Microwaves; Wells for In Situ Extraction of Volatiles from Regolith (WIEVR); and Estimating the Backup Reaction Wheel Orientation Using Reaction Wheel Spin Rates Flight Telemetry from a Spacecraft.

Source record↗
DOE OSTI2018

Final Report on CFD and Thermal-Mechanical Stress Analysis of PWR Surge Line under Transient Condition Thermal Stratification and an Evolutionary Cyclic Plasticity Based Transformative Fatigue Evaluation Approach without Using S~N Curve: Rev. 1

This report presents an evolutionary cyclic plasticity based transformative approach for fatigue evaluation of safety critical components. This approach is generic which can be used not only for the present generation light water reactor components (which is our focus) but also for advanced reactor components. The approach is based on the fundamental concept of time-evolution of progressive fatigue damage rather than the end-of-life data based conventional S~N curve approaches. Series of fatigue tests on 316 stainless steel uniaxial specimens were conducted under strain-controlled and stress-controlled cyclic loading in two different environment such as in-air at 300 °C and PWR coolant water condition at 300 °C. Material models are developed to capture the uniaxial test data and are subsequently programmed into commercially available ABAQUS code to transform the uniaxial material behavior to multiaxial loading domain. The developed evolutionary cyclic-plasticity approach is verified by both analytical and 3D-FE modeling of uniaxial fatigue test specimens. Results show that the developed material models can capture the time-dependence of material hardening and softening with great accuracy. In addition, results show that the material models not only can capture the material behavior under constant amplitude loading but also can capture the load-sequence effect under variable and random amplitude loading. The ANL developed approach is one of its kind, which shows first in the fatigue research community that fatigue evaluation of a component (in our case verified for laboratory scale fatigue test specimens with great accuracy) can be performed for its entire fatigue life, including the hardening and softening behavior, without using a conventional S~N curve based approach. Additionally results are presented related to fatigue evaluation of real reactor components such as of a pressurized water reactor (PWR) surge line (SL), under upset transient conditions. SLs of PWR experience complex thermal stratification due to the mixing of differential temperature flow from hot leg and pressurizer. A detailed parametric study of thermal stratification in SL is performed through both standalone computational fluid dynamics (CFD) analysis and fluid-structure-interaction based coupled CFD and heat transfer analysis. Results show that the mass flow rate in SL substantially affects the extent of its thermal stratification. The resulting thermal boundaries conditioned are used as inputs for thermal-mechanical stress analysis and component level fatigue evaluation. Fatigue evaluation of SL pipe under example loading conditions was performed using the above-mentioned evolutionary-plasticity based approach and ASME code (for in-air life) and NUREG-6909 (for PWR environment life) based approaches. Results show that the evolutionary plasticity based approach estimate substantially different life compared to the conventional ASME based approach. To show the entirety of the proposed evolutionary fatigue-modeling framework, this report includes the results presented in semi-annual March report (ANL/LWRS-17/01). Sections 2-5 of March report are directly included as Sections 2-5 in this report, without any significant changes. However, Sections 6-9 of this report present the work performed during March-September of 2017. Whereas, Sections 1 and 10 present the organization of this report and the summary and future work, respectively. This report (ANL/LWRS-17/03 Rev. 1) is a slightly revised version of the earlier published report (ANL/LWRS-17/03). In this revised version, some of the results and associated discussion have been corrected/updated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1996

Workshop on evolution of igneous asteroids: Focus on Vesta and the HED meteorites

Recently, the geology of the surface of Vesta has been coming to light. In 1983 Gaffey first began showing maps of the surface geology of Vesta constructed from numerous spectra obtained at different times as the asteroid rotated. By noting the details of spectral variation with rotation, he was able to develop two possible gross-scale geologic maps of Vesta showing the distributions of mafic and ultramafic materials. These maps were published in 1997. Finally, the capabilities of the Hubble Space Telescope were brought to bear on Vesta and images with a resolution of about 50 km were obtained using four different filters by Binzel and co-workers. Maps produced by this team published in 1997 began to reveal the geology of Vesta in sufficient detail that crude interpretations of the geologic history of the asteroid could be attempted. Additionally, in 1993 Binzel and Xu published a study of small asteroids in the region near Vesta in orbital-element space. In this study, they showed that there are a number of asteroids a few kilometers in size with reflectance spectra like that of Vesta that form a trail in orbital-element space from near Vesta to near resonances that can more easily supply material to near-Earth space. Binzel and Xu thus concluded that these small asteroids were spells of Vesta ejected by impact and that some of their brethren had been perturbed to Earth-approaching orbits. They suggested that these latter were the immediate parents of HED meteorites. This seemed to remove a long-standing dynamical objection to Vesta as the HED parent body, as discussed by Wasson and Wetherill in 1979. Within the last few years, NASA has initiated the Discovery program of low-cost, rapid-timescale development, exploration missions. Vesta has been proposed as an object worthy of study by a Discovery mission, although a Vesta mission has not yet been selected. With all the recent activity aimed at studying Vesta and the HED meteorites, and the possibility of a space mission to Vesta, we felt that time was ripe to convene a workshop bringing together astronomers, meteoriticists, and planetary geologists to focus on what could be learned about the geologic evolution of Vesta through integrating astronomical and HED meteorite studies. This, of course, assumes that the HED meteorites are from Vesta, and this issue was specifically addressed (but not resolved) in the workshop. Indeed, it seems likely that this issue can only be resolved by returning samples from Vesta for detailed study on Earth. The workshop was held at the LPI on October 16-18,1996, and was attended by some 70 scientists. Sessions included a set of talks on Earth- and space-based astronomical observations of Vesta plus the evidence pro and con for Vesta being the HED parent body, talks on the petrology and geochemistry of HED meteorites, talks on the formation and dynamics of ejecta from Vesta, talks on the thermal history of asteroids and HED meteorites, volcanic processes and differentiation history, and a short session devoted to possible missions to Vestal By all accounts, the workshop was considered a great success, although this is the opinion of a biased set of observers.

Mittlefehldt, D. W.↗
NASA NTRS2017

Effects of 2G on Gene Expression of Stress-Related Hormones in Rat Placenta

Understanding the effects of spaceflight on mammalian reproductive and developmental physiology is important to future human space exploration and permanent settlement beyond Earth orbit. Fetal developmental programming, including modulation of the HPA axis, is thought to originate at the placental-uterine interface, where both transfer of maternal hormones to the fetus and synthesis of endogenous hormones occurs. In healthy rats, fetal corticosterone levels are kept significantly lower by 11BetaHSD-2, which inactivates corticosterone by conversion into cortisone. Placental tissues express endogenous HPA axis-associated hormones including corticotropin-releasing hormone (CRH), pre-opiomelanocortin (POMC), and vasopressin, which may contribute to fetal programming alongside maternal hormones. DNA methylase 3A, 11BetaHSD-2, and 11BetaHSD-1, which are involved in the regulation of maternal cortisol transfer and modulation of the HPA axis, are also expressed in placental tissues along with glucocorticoid receptor and may be affected by differential gravity exposure during pregnancy. Fetuses may respond differently to maternal glucocorticoid exposure during gestation through sexually dimorphic expression of corticosterone-modulating hormones. To elucidate effects of altered gravity on placental gene expression, here we present a ground-based analogue study involving continuous centrifugation to produce 2g hypergravity. We hypothesized that exposure to 2g would induce a decrease in 11BetaHSD-2 expression through the downregulation of DNA methylase 3a and GC receptor, along with concurrent upregulation in endogenous CRH, POMC, and vasopressin expression. Timed pregnant female rats were exposed to 2G from Gestational day 6 to Gestational day 20, and comparisons made with Stationary Control (SC) and Vivarium Control (VC) dams at 1G. Dams were euthanized and placentas harvested on G20. We homogenized placental tissues, extracted and purified RNA, synthesized cDNA, and quantified the expression levels of the genes of interest relative to the GAPDH housekeeping gene, using RT-qPCR and gene-specific cDNA probes. Elucidation of glucocorticoid transfer and synthesis in the placenta can provide new insights into the unique dynamics of mammalian development in microgravity and guide future multi-generational studies in space.

Benson, S.↗
DOE OSTI2016

FY 2016 Status Report: Documentation of All CIRFT Data including Hydride Reorientation Tests (Draft M2)

The first portion of this report provides a detailed description of fiscal year (FY) 2015 test result corrections and analysis updates based on FY 2016 updates to the Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) program methodology, which is used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal conditions of transport (NCT). The CIRFT consists of a U-frame test setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages connecting to a universal testing machine. The curvature SNF rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are clamped to the side connecting plates of the U-frame and used to capture deformation of the rod. The second portion of this report provides the latest CIRFT data, including data for the hydride reorientation test. The variations in fatigue life are provided in terms of moment, equivalent stress, curvature, and equivalent strain for the tested SNFs. The equivalent stress plot collapsed the data points from all of the SNF samples into a single zone. A detailed examination revealed that, at the same stress level, fatigue lives display a descending order as follows: H. B. Robinson Nuclear Power Station (HBR), LMK, and mixed uranium-plutonium oxide (MOX). Just looking at the strain, LMK fuel has a slightly longer fatigue life than HBR fuel, but the difference is subtle. The third portion of this report provides finite element analysis (FEA) dynamic deformation simulation of SNF assemblies . In a horizontal layout under NCT, the fuel assembly’s skeleton, which is formed by guide tubes and spacer grids, is the primary load bearing apparatus carrying and transferring vibration loads within an SNF assembly. These vibration loads include interaction forces between the SNF assembly and the canister basket walls. Therefore, the integrity of the guide tubes and spacer grids critically affects the vibration intensity of the fuel assembly during transport and must be considered when developing the multipurpose purpose canister (MPC) design for safe SNF transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2019

Modal Analysis of the Orion Capsule Two Parachute System

As discussed in Ref [1], it is apparent from flight tests that the system made up of two main parachutes and a capsule can undergo several distinct dynamical behaviors. The most significant and problematic of these is the pendulum mode in which the system develops a pronounced swinging motion with an amplitude of up to 24 deg. Large excursions away from vertical by the capsule could cause it to strike the ground at a large horizontal or vertical speed and jeopardize the safety of the astronauts during a crewed mission. In reference [1], Ali et al. summarized a series of efforts taken by the Capsule Parachute Assembly System (CPAS) Program to understand and mitigate the pendulum issue. The period of oscillation and location of the system's pivot point are determined from post-flight analysis. Other noticeable but benign modes include: 1) flyout (scissors) mode, where the parachutes move back and forth symmetrically with respect to the vertical axis similar to the motion of a pair of scissors; 2) maypole mode, where the two parachutes circle around the vertical axis at a nearly constant radius and period; and 3) breathing mode, in which deformation of the non-rigid canopies affects the axial acceleration of the system in an oscillatory manner. Because these modes are relatively harm- less, little effort has been devoted to analyzing them in comparison with the pendulum motion. Motions of the actual system made up of two parachutes and a capsule are extremely complicated due to nonlinearities and flexibility effects. Often it is difficult to obtain insight into the fundamental dynamics of the system by examining results from a multi-body simulation based on nonlinear equations of motion (EOMs). As a part of this study, the dynamics of each mode observed during flight is derived from first principles on an individual basis by making numerous simplifications along the way. The intent is to gain a better understanding into the behavior of the complex multi-body system by studying the reduced set of differential equations associated with each mode. This approach is analogous to the traditional modal analysis technique used to study airplane flight dynamics, in which the full nonlinear behavior of the airframe is decomposed into the phugoid and short period modes for the longitudinal dynamics and the spiral, roll-subsidence, and dutch-roll modes for the lateral dynamics. It is important to note that the study does not address the mechanisms that cause the system to transition from one mode to another, nor does it discuss motions during which two or more modes occur simultaneously.

Pei, J.↗
DOE OSTI2024

Investigating the Influence of River Geomorphology on Human Presence Using Night Light Data: A Case Study in the Indus Basin

Human settlements have historically thrived near rivers due to enhanced navigation and trade, and the availability of water supply and resources. The use of night light data, representing economic activities, provides a novel approach to studying the interactions between human activity and rivers over time. Here, we use the Defense Meteorological Satellite Program (DMSP) stable night light data from 2000 to 2013 as a proxy for human presence and activities to quantify the statistical relationships between night light presence and intensity in the Indus Basin, Asia. We test how these data are affected by proximity to trunk channels and by channel type (single/multi-thread) in the study area. We find that night light presence is enhanced by 26% within a 0 to 5 km proximity range of the Indus River and its tributaries, relative to the basin as a whole. We interpret this to represent increased human presence and activity within this zone. However, the mean intensity is lower near the river and higher away from the river, signifying denser settlements, such as towns and cities, which are preferentially located away from the Indus and its tributaries. Moreover, the enhancement of lit pixels signifying human presence and activities is increased by 18% near single-thread sections of the Indus River, compared to segments of the Indus displaying multi-thread morphologies. We suggest that this is due to the enhanced stability of single-threaded channels, relative to mobile multi-threaded channel reaches. This study demonstrates how night lights are an important tool in studying the relationship between human presence and river dynamics in large catchments such as the Indus, and we suggest that these data will have an important role in assessing differential flood spatial and social vulnerability at a regional scale.

Environmental Sciences & Ecology↗
NASA NTRS1998

X-33 Integrated Test facility, Extended Range Simulation

In support of the X-33 Single Stage To Orbit program, NASA Dryden Flight Research Center was selected to provide continuous communications coverage of the X-33 vehicle from launch, through landing at Malmstrom Air Force Base, Montana and Michaels Army Air Field, Utah. An extensive real-time range simulation capability is being developed to ensure successful communications with the autonomous X-33 vehicle. This paper will provide an overview of the various levels of simulation being developed to support the X-33 extended range subsystems. These subsystems include the Flight Termination System, L-Band command uplink subsystem and the S-Band telemetry downlink subsystem. In addition, the radar model developed provides continuous azimuth, elevation and range information based on the flight trajectory. The Dynamic Ground Station Analysis model developed by NASA Goddard Space Flight Center, calculate the received signal strength at each ground station. This model takes into consideration Radio Frequency (RF) link parameters such as frequency, antenna gain, space loss, plasma effects and the vehicle's position and attitude at any point in time during the flight path. All three RF links are then attenuated based on this calculated level and the RF signals are sent into telemetry receivers to emulate remote sites, or the power incident on the vehicle from uplinked signals. The best source received telemetry data is then passed back to the Launch and Mission Control Monitoring System (LMCMS) resident in the Operations Control Center. The LMCMS also provides the range simulation system the uplink command combined with differential GPS corrections. Later stages will require the progressive integration of actual range hardware with this simulation effort, leading to communication between telemetry, uplink and FTS antennas at NASA Dryden Flight Research Center, with vehicle antennas mounted on the Walter C. Williams Research Aircraft Integration Facility (RAIF). Decommutated Pulse Code Modulated (PCM) data is displayed on one of the four monitors that comprise the Range Safety Officer's (RSO) station. Also displayed are instantaneous impact prediction models, and Federal Aviation Administration (FAA) data for notification of other traffic in the area. Aside from initiating the flight termination command and validating communication links, the RSO station with the range simulation will be used to provide both range control and range safety officers training. The training is necessary to perform their respective functions with greater levels of confidence prior to first flight.

Sharma, Ashley↗
NASA NTRS2001

Venus Multiprobe Mission

A summary of the proposal for the Discovery Program Venus Multiprobe Mission (VMPM) which was submitted in 1994 will be presented. The major goal of the mission was to provide sufficient observations of winds and thermal structure in the Venus atmosphere that the mechanism producing the atmospheric superrotation could be identified. The superrotation has been a major unsolved problem of planetary science for over 30 years. In spite of many theoretical efforts and many spacecraft missions to Venus, the principal driving mechanism(s) remain undetermined. Sixteen small probes, modeled after the Pioneer Venus small probes, would enter the Venus atmosphere, and measure winds, temperatures, and pressures from 60 km to the surface. Towards the bottom of the atmosphere, measurements would be made at 10 meter intervals. Temperature and pressure measurements are similar to those used on Pioneer Venus, but are four times more sensitive. Winds would be obtained by dual-frequency, differential. long-baseline interferometry, using S- and X-band transmitters on the probes, and four receiving telescopes round the Pacific rim. The precision of wind measurements is 0.05 ms- 1 or better, at least twenty times more precise than for Pioneer Venus. Ground based observations occurring in conjunction with the multiprobe entries would provide a temporal context for the upper boundary to the probe observations. UV and near IR features would be tracked for periods of hours to weeks to obtain cloud level zonal mean and transient eddy circulations with accuracies of 10 /ms or better. The data would be processed into wind and opacity maps. A hierarchy of atmospheric numerical models would be used to distinguish different dynamical elements existing in the data. Primary among the models was the development of a Venus general circulation model containing all known physical processes.

Young, Richard E.↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
NASA NTRS1991

System monitoring and diagnosis with qualitative models

A substantial foundation of tools for model-based reasoning with incomplete knowledge was developed: QSIM (a qualitative simulation program) and its extensions for qualitative simulation; Q2, Q3 and their successors for quantitative reasoning on a qualitative framework; and the CC (component-connection) and QPC (Qualitative Process Theory) model compilers for building QSIM QDE (qualitative differential equation) models starting from different ontological assumptions. Other model-compilers for QDE's, e.g., using bond graphs or compartmental models, have been developed elsewhere. These model-building tools will support automatic construction of qualitative models from physical specifications, and further research into selection of appropriate modeling viewpoints. For monitoring and diagnosis, plausible hypotheses are unified against observations to strengthen or refute the predicted behaviors. In MIMIC (Model Integration via Mesh Interpolation Coefficients), multiple hypothesized models of the system are tracked in parallel in order to reduce the 'missing model' problem. Each model begins as a qualitative model, and is unified with a priori quantitative knowledge and with the stream of incoming observational data. When the model/data unification yields a contradiction, the model is refuted. When there is no contradiction, the predictions of the model are progressively strengthened, for use in procedure planning and differential diagnosis. Only under a qualitative level of description can a finite set of models guarantee the complete coverage necessary for this performance. The results of this research are presented in several publications. Abstracts of these published papers are presented along with abtracts of papers representing work that was synergistic with the NASA grant but funded otherwise. These 28 papers include but are not limited to: 'Combined qualitative and numerical simulation with Q3'; 'Comparative analysis and qualitative integral representations'; 'Model-based monitoring of dynamic systems'; 'Numerical behavior envelopes for qualitative models'; 'Higher-order derivative constraints in qualitative simulation'; and 'Non-intersection of trajectories in qualitative phase space: a global constraint for qualitative simulation.'

Kuipers, Benjamin↗
DOE OSTI2024

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

FY 2016 Status Report: CIRFT Testing Data Analyses and Updated Curvature Measurements

This report provides a detailed description of FY15 test result corrections/analysis based on the FY16 Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) test program methodology update used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal transportation conditions. The CIRFT consists of a U-frame testing setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages to a universal testing machine. The curvature of rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are used and clamped to the side connecting plates of the U-frame to capture the deformation of the rod. The contact-based measurement, or three-LVDT-based curvature measurement system, on SNF rods has been proven to be quite reliable in CIRFT testing. However, how the LVDT head contacts the SNF rod may have a significant effect on the curvature measurement, depending on the magnitude and direction of rod curvature. It has been demonstrated that the contact/curvature issues can be corrected by using a correction on the sensor spacing. The sensor spacing defines the separation of the three LVDT probes and is a critical quantity in calculating the rod curvature once the deflections are obtained. The sensor spacing correction can be determined by using chisel-type probes. The method has been critically examined this year and has been shown to be difficult to implement in a hot cell environment, and thus cannot be implemented effectively. A correction based on the proposed equivalent gauge-length has the required flexibility and accuracy and can be appropriately used as a correction factor. The correction method based on the equivalent gauge length has been successfully demonstrated in CIRFT data analysis for the dynamic tests conducted on Limerick (LMK) (17 tests), North Anna (NA) (6 tests), and Catawba mixed oxide (MOX) (10 tests) SNF samples. These CIRFT tests were completed in FY14 and FY15. Specifically, the data sets obtained from measurement and monitoring were processed and analyzed. The fatigue life of rods has been characterized in terms of moment, curvature, and equivalent stress and strain..

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Simulator for an Accelerator-Driven Subcritical Fissile Solution System

LANL has developed a process to generate a progressive family of system models for a fissile solution system. This family includes a dynamic system simulation comprised of coupled nonlinear differential equations describing the time evolution of the system. Neutron kinetics, radiolytic gas generation and transport, and core thermal hydraulics are included in the DSS. Extensions to explicit operation of cooling loops and radiolytic gas handling are embedded in these systems as is a stability model. The DSS may then be converted to an implementation in Visual Studio to provide a design team the ability to rapidly estimate system performance impacts from a variety of design decisions. This provides a method to assist in optimization of the system design. Once design has been generated in some detail the C++ version of the system model may then be implemented in a LabVIEW user interface to evaluate operator controls and instrumentation and operator recognition and response to off-normal events. Taken as a set of system models the DSS, Visual Studio, and LabVIEW progression provides a comprehensive set of design support tools.

42 ENGINEERING↗
DOE OSTI2015

Laboratory Experiments and Instrument Intercomparison Studies of Carbonaceous Aerosol Particles

Aerosols containing black carbon (and some specific types of organic particulate matter) directly absorb incoming light, heating the atmosphere. In addition, all aerosol particles backscatter solar light, leading to a net-cooling effect. Indirect effects involve hydrophilic aerosols, which serve as cloud condensation nuclei (CCN) that affect cloud cover and cloud stability, impacting both atmospheric radiation balance and precipitation patterns. At night, all clouds produce local warming, but overall clouds exert a net-cooling effect on the Earth. The effect of aerosol radiative forcing on climate may be as large as that of the greenhouse gases, but predominantly opposite in sign and much more uncertain. The uncertainties in the representation of aerosol interactions in climate models makes it problematic to use model projections to guide energy policy. The objective of our program is to reduce the uncertainties in the aerosol radiative forcing in the two areas highlighted in the ASR Science and Program Plan. That is, (1) addressing the direct effect by correlating particle chemistry and morphology with particle optical properties (i.e. absorption, scattering, extinction), and (2) addressing the indirect effect by correlating particle hygroscopicity and CCN activity with particle size, chemistry, and morphology. In this connection we are systematically studying particle formation, oxidation, and the effects of particle coating. The work is specifically focused on carbonaceous particles where the uncertainties in the climate relevant properties are the highest. The ongoing work consists of laboratory experiments and related instrument inter-comparison studies both coordinated with field and modeling studies, with the aim of providing reliable data to represent aerosol processes in climate models. The work is performed in the aerosol laboratory at Boston College. At the center of our laboratory setup are two main sources for the production of aerosol particles: (a) two well-characterized source of soot particles and (b) a flow reactor for controlled OH and/or O3 oxidation of relevant gas phase species to produce well-characterized SOA particles. After formation, the aerosol particles are subjected to physical and chemical processes that simulate aerosol growth and aging. A suite of instruments in our laboratory is used to characterize the physical and chemical properties of aerosol particles before and after processing. The Time of Flight Aerosol Mass Spectrometer (ToF-AMS) together with a Scanning Mobility Particle Sizer (SMPS) measures particle mass, volume, density, composition (including black carbon content), dynamic shape factor, and fractal dimension. The–ToF-AMS was developed at ARI with Boston College participation. About 120 AMS instruments are now in service (including 5 built for DOE laboratories) performing field and laboratory studies world-wide. Other major instruments include a thermal denuder, two Differential Mobility Analyzers (DMA), a Cloud Condensation Nuclei Counter (CCN), a Thermal desorption Aerosol GC/MS (TAG) and the new Soot Particle Aerosol Mass Spectrometer (SP-AMS). Optical instrumentation required for the studies have been brought to our laboratory as part of ongoing and planned collaborative projects with colleagues from DOE, NOAA and university laboratories. Optical instruments that will be utilized include a Photoacoustic Spectrometer (PAS), a Cavity Ring Down Aerosol Extinction Spectrometer (CRD-AES), a Photo Thermal Interferometer (PTI), a new 7-wavelength Aethalometer and a Cavity Attenuated Phase Shift Extinction Monitor (CAPS). These instruments are providing aerosol absorption, extinction and scattering coefficients at a range of atmospherically relevant wavelengths. During the past two years our work has continued along the lines of our original proposal. We report on 12 completed and/or continuing projects conducted during the period 08/14 to 0814/2015. These projects are described in 17 manuscripts published in refereed journals.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Observations of Fast-Moving Features in the Debris Disk of AU Mic on a Three-Year Timescale: Confirmation and New Discoveries

Context. The nearby and young M star AU Mic is surrounded by a debris disk in which we previously identified a series of large-scale arch-like structures that have never been seen before in any other debris disk and that move outward at high velocities. Aims. We initiated a monitoring program with the following objectives: (1) track the location of the structures and better constrain their projected speeds, (2) search for new features emerging closer in, and ultimately (3) understand the mechanism responsible for the motion and production of the disk features. Methods. AU Mic was observed at 11 different epochs between August 2014 and October 2017 with the IR camera and spectrograph of SPHERE. These high-contrast imaging data were processed with a variety of angular, spectral, and polarimetric differential imaging techniques to reveal the faintest structures in the disk. We measured the projected separations of the features in a systematic way for all epochs. We also applied the very same measurements to older observations from the Hubble Space Telescope (HST) with the visible cameras STIS and ACS. Results. The main outcomes of this work are (1) the recovery of the five southeastern broad arch-like structures we identified in our first study, and confirmation of their fast motion (projected speed in the range 4–12 km/s); (2) the confirmation that the very first structures observed in 2004 with ACS are indeed connected to those observed later with STIS and now SPHERE; (3) the discovery of two new very compact structures at the northwest side of the disk (at 0.40′′ and 0.55′′ in May 2015) that move to the southeast at low speed; and (4) the identification of a new arch-like structure that might be emerging at the southeast side at about 0.4′′ from the star (as of May 2016). Conclusions. Although the exquisite sensitivity of SPHERE allows one to follow the evolution not only of the projected separation, but also of the specific morphology of each individual feature, it remains difficult to distinguish between possible dynamical scenarios that may explain the observations. Understanding the exact origin of these features, the way they are generated, and their evolution over time is certainly a significant challenge in the context of planetary system formation around M stars.

Boccaletti, A.↗