Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 793 records · Page 44

Advances in Westinghouse fuel rod design technology (PAD5)

The Westinghouse Performance Analysis and Design (PAD) code for light-water reactors has been recently updated to version PAD5 and submitted to NRC for review. PAD5 incorporates extensive updates in fuel performance and material property models, including fuel thermal conductivity, pellet fragment relocation, fission gas release, fission gas swelling, and cladding creep and growth. Advances in thermal models are mainly in pellet fragment relocation, radial power distribution within the pellet, and the changes in fuel thermal conductivity due to changes in temperature, burnup and fuel porosity. The PAD5 thermal fission gas release model is a phenomenological model based on the irradiation dependence of the gas release from mechanistic processes for the diffusion of the fission gas generated in the fuel grains to the grain boundaries, formation of fission gas bubbles on the grain boundaries and release of the fission gas when these grain boundary bubbles inter-link. The PAD5 cladding creep model is predominantly a combination of the following mechanisms: thermal creep, irradiation induced creep, and irradiation enhanced creep. The fuel performance database used for the development, calibration and testing of the PAD5 code covers a broad range of both commercial irradiation data and test rod data and requires continuously additional data focusing on the latest development of the PAD5 code. Test programs are also included in the database. These data are checked and stored in a hierarchical database structure that supports the automation of model development and calibration. New data, once incorporated into the database, can be quickly evaluated with respect to key fuel and cladding performance models to confirm expected trends or to evaluate potential model improvements. This automated calibration system has been integrated with state-of-the-art optimization packages in DAKOTA Analysis Kit to provide for advanced statistical assessments and model coefficient optimization. Calibration results are automatically tabulated and plotted to provide easy assessment of trends or new behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification of iterative matrix solutions for multipoint kinetics equations

Multipoint kinetic equation systems have been solved numerically using matrix algebra software and 1st to 4th order implicit schemes based on single-step matrix propagators. These matrix schemes have been validated successfully on demanding point kinetics benchmarks with various prescribed reactivity insertions. Verification tests have also been performed on a simple 3-region fast reactor core model with asymmetric step or ramp reactivity insertions and cross-checked with the multipoint kinetics SACRE code developed at INL and based on existing and validated stiff ODE solver packages. Accurate results can be obtained at a limited time expense. The reason for intriguing results obtained in some transients has been elucidated and linked to the multipoint matrix coefficient interpolation method during the transient.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Erroneous predictions of auxotrophies by CarveMe

We report in their paper, D. Machado et al. claim to identify ‘cooperative’ communities with many amino acid auxotrophies1. A related paper, with many of the same authors, claims that 26% of bacteria from free-living communities are amino acid auxotrophs2. The putative amino acid auxotrophies were derived from genome sequences via the metabolic models of CarveMe3 but were not validated in any way. I checked the predictions of CarveMe for prototrophic bacteria and found a high rate of spurious auxotrophies.

59 BASIC BIOLOGICAL SCIENCES↗

A UNIVERSAL POWER-LAW PROFILE OF PSEUDO-PHASE-SPACE DENSITY-LIKE QUANTITIES IN ELLIPTICAL GALAXIES

We study profiles of mass density, velocity dispersion (VD), and a combination of both using ∼2000 nearly spherical and rotation-free Sloan Digital Sky Survey galaxies. For observational stellar mass density ρ{sub *}(r), we consider a range of dark matter (DM) distribution ρ{sub DM}(r) and VD anisotropy β(r) to investigate radial stellar VD σ{sub *r}(r) using the spherical Jeans equation. While mass and VD profiles vary appreciably depending on DM distribution and anisotropy, the pseudo-phase-space density-like combination ρ(r)/σ{sub ⋆r}{sup 3}(r) with total density ρ(r) = ρ{sub *}(r) + ρ{sub DM}(r) is nearly universal. In the optical region, the negative logarithmic slope has a mean value of (χ) ≈ 1.86-1.90 with a galaxy-to-galaxy rms scatter of ≈0.04-0.06, which is a few times smaller than that of ρ(r) profiles. The scatter of χ can be increased by invoking wildly varying anisotropies that are, however, less likely because they would produce too large a scatter of line of sight VD profiles. As an independent check of this universality, we analyze stellar orbit-based dynamical models of 15 early-type galaxies (ETGs) of the Coma cluster provided by J. Thomas. Coma ETGs, with σ{sub *r}(r) replaced by the rms velocity of stars v {sub *rms}(r) including net rotation, exhibit a similar universality with a slope of χ = 1.93 ± 0.06. Remarkably, the inferred values of χ for ETGs match well the slope ≈1.9 predicted by N-body simulations of DM halos. We argue that the inferred universal nature of ρ(r)/σ{sub ⋆r}{sup 3}(r) cannot be fully explained by equilibrium alone, implying that some astrophysical factors conspire and/or it reflects a fundamental principle in collisionless formation processes.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey year 1 results: Cosmological constraints from galaxy clustering and weak lensing

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg2 of griz imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric-redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while “blind” to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat ΛCDM and wCDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for ΛCDM) or 7 (for wCDM) cosmological parameters including the neutrino mass density and including the 457×457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions and from their combination obtain S8≡σ8(Ωm/0.3)0.5=0.773-0.020+0.026 and Ωm=0.267-0.017+0.030 for ΛCDM; for wCDM, we find S8=0.782-0.024+0.036, Ωm=0.284-0.030+0.033, and w=-0.82-0.20+0.21 at 68% C.L. The precision of these DES Y1 constraints rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for S8 and Ωm are lower than the central values from Planck for both ΛCDM and wCDM, the Bayes factor indicates that the DES Y1 and Planck data sets are consistent with each other in the context of ΛCDM. Combining DES Y1 with Planck, baryonic acoustic oscillation measurements from SDSS, 6dF, and BOSS and type Ia supernovae from the Joint Lightcurve Analysis data set, we derive very tight constraints on cosmological parameters: S8=0.802±0.012 and Ωm=0.298±0.007 in ΛCDM and w=-1.00-0.04+0.05 in wCDM. Upcoming Dark Energy Survey analyses will provide more stringent tests of the ΛCDM model and extensions such as a time-varying equation of state of dark energy or modified gravity.

79 ASTRONOMY AND ASTROPHYSICS↗

What does fault tolerant Deep Learning need from MPI?

Deep Learning (DL) algorithms have become the {\em de facto} Machine Learning (ML) algorithm for large scale data analysis. DL algorithms are computationally expensive -- even distributed DL implementations which use MPI require days of training (model learning) time on commonly studied datasets. Long running DL applications become susceptible to faults -- requiring development of a fault tolerant system infrastructure, in addition to fault tolerant DL algorithms. This raises an important question: {\em What is needed from MPI for designing fault tolerant DL implementations?} In this paper, we address this problem for permanent faults. We motivate the need for a fault tolerant MPI specification by an in-depth consideration of recent innovations in DL algorithms and their properties, which drive the need for specific fault tolerance features. We present an in-depth discussion on the suitability of different parallelism types (model, data and hybrid); a need (or lack thereof) for check-pointing of any critical data structures; and most importantly, consideration for several fault tolerance proposals (user-level fault mitigation (ULFM), Reinit) in MPI and their applicability to fault tolerant DL implementations. We leverage a distributed memory implementation of Caffe, currently available under the Machine Learning Toolkit for Extreme Scale (MaTEx). We implement our approaches by extending MaTEx-Caffe for using ULFM-based implementation. Our evaluation using the ImageNet dataset and AlexNet neural network topology demonstrates the effectiveness of the proposed fault tolerant DL implementation using OpenMPI based ULFM.

Amatya, Vinay C.↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identification and characterization of structural nonlinearities in the Space Shuttle Main Engine

Structural nonlinearities should be considered in the analysis and modeling of mechanical systems for accuracy in predicting dynamic behavior. Moreover, checks for nonlinearity should address the type and class of nonlinearity in the structure and identify the probable cause. A series of experimental modal survey tests were carried out to determine and characterize the structural nonlinearities of the Space Shuttle Main Engine main injector. Different excitation regimes at various amplitudes were utilized to induce vibration and strain, and acceleration measurements were taken at judiciously selected locations. The test results indicated that both stiffness and damping nonlinearities were present.

Panossian, H. V.↗

Spacecraft command verification: The AI solution

Recently, a knowledge-based approach was used to develop a system called the Command Constraint Checker (CCC) for TRW. CCC was created to automate the process of verifying spacecraft command sequences. To check command files by hand for timing and sequencing errors is a time-consuming and error-prone task. Conventional software solutions were rejected when it was estimated that it would require 36 man-months to build an automated tool to check constraints by conventional methods. Using rule-based representation to model the various timing and sequencing constraints of the spacecraft, CCC was developed and tested in only three months. By applying artificial intelligence techniques, CCC designers were able to demonstrate the viability of AI as a tool to transform difficult problems into easily managed tasks. The design considerations used in developing CCC are discussed and the potential impact of this system on future satellite programs is examined.

Fesq, Lorraine M.↗

Rationalizing Burned Carbon with Carbon Monoxide Exported from South America

We present several estimates cross-checking the fluxes of carbon to the atmosphere from burning, comparing models that are based on simple land-surface parameterizations and atmospheric transport dynamics. Both estimates made by NASA Ames and USP modeling techniques are quite high compared to some detailed satellite/land-use studies of emissions. The flux of carbon liberated to the atmosphere via biomass burning is important for several reasons. This flux is a fundamental statistic for the parameterization of the large-scale flux of gases controlling the reactive greenhouse gases methane and ozone. Similarly, it is central to the estimation of the translocation of nitrogen and pyrodenitrification in the tropics. Thirdly, CO2 emitted from rainforest clearing contributes directly to carbon lost from the rainforest system as it contributes to greenhouse gas forcing. While CO2 from pasturage, agriculture, etc, is considered to be reabsorbed seasonally, and so "off budget" for the carbon cycle, it must also be accounted. CO2 anomalies related to daily weather and interannual climatic variation are strong enough to perturb our scientific perception of long-term carbon storage trends. We compare fluxes deduced from land-use statistics (originally, W.M. Hao) and from satellite hot pixels (A. Setzer) with atmospheric fluxes determined by the mesoscale/continental scale models RAMS and MM5, and point to some new work with highly resolved global models (the NASA Data Assimilation Office's GEOS4). Our simulations are tied to events, so that measured tracers like CO tie the models directly to the burning and meteorology of a specific period. We point out a particular sensitivity in estimates based on CO, and indicate how analysis of CO2 along with other biomass-burning tracers may lead to an improved multi-species estimator of carbon burned.

Chatfield, R.↗

The ANMLite Language and Logic for Specifying Planning Problems

We present the basic concepts of the ANMLite planning language. We discuss various aspects of specifying a plan in terms of constraints and checking the existence of a solution with the help of a model checker. The constructs of the ANMLite language have been kept as simple as possible in order to reduce complexity and simplify the verification problem. We illustrate the language with a specification of the space shuttle crew activity model that was constructed under the Spacecraft Autonomy for Vehicles and Habitats (SAVH) project. The main purpose of this study was to explore the implications of choosing a robust logic behind the specification of constraints, rather than simply proposing a new planning language.

Butler, Ricky W.↗

SU-E-T-51: Bayesian Network Models for Radiotherapy Error Detection

Purpose: To develop a probabilistic model of radiotherapy plans using Bayesian networks that will detect potential errors in radiation delivery. Methods: Semi-structured interviews with medical physicists and other domain experts were employed to generate a set of layered nodes and arcs forming a Bayesian Network (BN) which encapsulates relevant radiotherapy concepts and their associated interdependencies. Concepts in the final network were limited to those whose parameters are represented in the institutional database at a level significant enough to develop mathematical distributions. The concept-relation knowledge base was constructed using the Web Ontology Language (OWL) and translated into Hugin Expert Bayes Network files via the the RHugin package in the R statistical programming language. A subset of de-identified data derived from a Mosaiq relational database representing 1937 unique prescription cases was processed and pre-screened for errors and then used by the Hugin implementation of the Estimation-Maximization (EM) algorithm for machine learning all parameter distributions. Individual networks were generated for each of several commonly treated anatomic regions identified by ICD-9 neoplasm categories including lung, brain, lymphoma, and female breast. Results: The resulting Bayesian networks represent a large part of the probabilistic knowledge inherent in treatment planning. By populating the networks entirely with data captured from a clinical oncology information management system over the course of several years of normal practice, we were able to create accurate probability tables with no additional time spent by experts or clinicians. These probabilistic descriptions of the treatment planning allow one to check if a treatment plan is within the normal scope of practice, given some initial set of clinical evidence and thereby detect for potential outliers to be flagged for further investigation. Conclusion: The networks developed here support the use of probabilistic models into clinical chart checking for improved detection of potential errors in RT plans.

60 APPLIED LIFE SCIENCES↗

Decision-directed adaptive recursive estimators - Divergence prevention.

A method is proposed for dealing with the divergence phenomenon problem encountered in applications of minimum-variance recursive estimators when the error covariance calculated by the estimator becomes inconsistent with the actual error covariance. The proposed method differs from previous ones in that the state model and statistics are accepted as given. The form of the optimal estimator is used, but a constant check is made on the consistency of the calculated and actual error covariances. The method is independent of the source of error, whether it be inaccuracies in the system model, incorrect values of the a priori and random process statistics, approximations required in the case of nonlinear systems, or computational roundoff. Several simulated examples, in which inconsistencies in the calculated and actual error covariances exist, show a significant improvement in the performance of the estimator when the given procedure is applied.

Nahi, N. E.↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

SeaWiFS technical report series. Volume 10: Modeling of the SeaWiFS solar and lunar observations

Post-launch stability monitoring of the Sea-viewing Wide Field-of-view Sensor (SeaWifs) will include periodic sweeps of both an onboard solar diffuser plate and the moon. The diffuser views will provide short-term checks and the lunar views will monitor long-term trends in the instrument's radiometric stability. Models of the expected sensor response to these observations were created on the SeaWiFS computer at the National Aeronautics and Space Administration's (NASA) Goddard Space Flight Center (GSFC) using the Interactive Data Language (IDL) utility with a graphical user interface (GUI). The solar model uses the area of intersecting circles to simulate the ramping of sensor response while viewing the diffuser. This model is compared with preflight laboratory scans of the solar diffuser. The lunar model reads a high-resolution lunar image as input. The observations of the moon are simulated with a bright target recovery algorithm that includes ramping and ringing functions. Tests using the lunar model indicate that the integrated radiance of the entire lunar surface provides a more stable quantity than the mean of radiances from centralized pixels. The lunar model is compared to ground-based scans by the SeaWiFS instrument of a full moon in December 1992. Quality assurance and trend analyses routines for calibration and for telemetry data are also discussed.

Woodward, Robert H.↗

The self streamlining wind tunnel

A two dimensional test section in a low speed wind tunnel capable of producing flow conditions free from wall interference is presented. Flexible top and bottom walls, and rigid sidewalls from which models were mounted spanning the tunnel are shown. All walls were unperforated, and the flexible walls were positioned by screw jacks. To eliminate wall interference, the wind tunnel itself supplied the information required in the streamlining process, when run with the model present. Measurements taken at the flexible walls were used by the tunnels computer check wall contours. Suitable adjustments based on streamlining criteria were then suggested by the computer. The streamlining criterion adopted when generating infinite flowfield conditions was a matching of static pressures in the test section at a wall with pressures computed for an imaginary inviscid flowfield passing over the outside of the same wall. Aerodynamic data taken on a cylindrical model operating under high blockage conditions are presented to illustrate the operation of the tunnel in its various modes.

Goodyer, M. J.↗

The DaveMLTranslator: An Interface for DAVE-ML Aerodynamic Models

It can take weeks or months to incorporate a new aerodynamic model into a vehicle simulation and validate the performance of the model. The Dynamic Aerospace Vehicle Exchange Markup Language (DAVE-ML) has been proposed as a means to reduce the time required to accomplish this task by defining a standard format for typical components of a flight dynamic model. The purpose of this paper is to describe an object-oriented C++ implementation of a class that interfaces a vehicle subsystem model specified in DAVE-ML and a vehicle simulation. Using the DaveMLTranslator class, aerodynamic or other subsystem models can be automatically imported and verified at run-time, significantly reducing the elapsed time between receipt of a DAVE-ML model and its integration into a simulation environment. The translator performs variable initializations, data table lookups, and mathematical calculations for the aerodynamic build-up, and executes any embedded static check-cases for verification. The implementation is efficient, enabling real-time execution. Simple interface code for the model inputs and outputs is the only requirement to integrate the DaveMLTranslator as a vehicle aerodynamic model. The translator makes use of existing table-lookup utilities from the Langley Standard Real-Time Simulation in C++ (LaSRS++). The design and operation of the translator class is described and comparisons with existing, conventional, C++ aerodynamic models of the same vehicle are given.

Hill, Melissa A.↗

Investigating decoherence in a simple system

The results of some simple calculations designed to study quantum decoherence are presented. The physics of quantum decoherence are briefly reviewed, and a very simple 'toy' model is analyzed. Exact solutions are found using numerical techniques. The type of incoherence exhibited by the model can be changed by varying a coupling strength. The author explains why the conventional approach to studying decoherence by checking the diagonality of the density matrix is not always adequate. Two other approaches, the decoherence functional and the Schmidt paths approach, are applied to the toy model and contrasted to each other. Possible problems with each are discussed.

Albrecht, Andreas↗