SEARCH · Engineering Papers
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.
Slides discussing general INL capabilities in a variety of research areas including LWRS, Plant Modernization, Risk-Informed Systems Analysis, the INL Digital Innovation Center of Excellence DICE, integrated online dynamic plant models, online deterministic failure-mode and reliability models of plant systems, equipment and components, and plant information models.
Recent works exploring deep learning application to dynamical systems modeling have demonstrated that embedding physical priors into neural networks can yield more effective, physically-realistic, and data-efficient models. However, in the absence of complete prior knowledge of a dynamical system's physical characteristics, determining the optimal structure and optimization strategy for these models can be difficult. In this work, we explore methods for discovering neural state space dynamics models for system identification. Starting with a design space of block-oriented state space models and structured linear maps with strong physical priors, we encode these components into a model genome alongside network structure, penalty constraints, and optimization hyperparameters. Demonstrating the overall utility of the design space, we employ an asynchronous genetic search algorithm that alternates between model selection and optimization and obtains accurate physically consistent models of three physical systems: an aerodynamics body, a continuous stirred tank reactor, and a two tank interacting system.
Laminated continuous fiber polymeric composites are capable of sustaining substantial load induced microstructural damage prior to component failure. Because this damage eventually leads to catastrophic failure, it is essential to capture the mechanics of progressive damage in any cogent life prediction model. For the past several years the authors have been developing one solution approach to this problem. In this approach the mechanics of matrix cracking and delamination are accounted for via locally averaged internal variables which account for the kinematics of microcracking. Damage progression is predicted by using phenomenologically based damage evolution laws which depend on the load history. The result is a nonlinear and path dependent constitutive model which has previously been implemented to a finite element computer code for analysis of structural components. Using an appropriate failure model, this algorithm can be used to predict component life. In this paper the model will be utilized to demonstrate the ability to predict the load path dependence of the damage and stresses in plates subjected to fatigue loading.
We propose a novel Coarse-Grained Reconfigurable Architecture-based (CGRA) QNN acceleration framework, CQNN. CQNN has a large number of basic components for binary functions. By programming CQNN at runtime according to the target QNN models, these basic components are integrated efficiently to support QNN functions with any data-width and hyper-parameter requirements and CQNN is reconfigured to have the optimal architecture for the target models. The framework includes compiler, architecture design, simulator and RTL generator. Experimental results show CQNNs can complete the inference of AlexNet and VGG-16 within 0.13ms and 2.63ms.
Optical measurements on an axisymmetrical quartz component engine research model were made to evaluate the flow field encountered during induction. The measurement technique is LIPA (Laser Induced Photochemical Anemometry), a non-intrusive velocimetry concept that provides an investigator of fluid flow with a tool to attain planar information about three-dimensional velocity and vorticity vectors in a single measurement step. The goal of this investigation is to further develop this measurement technique and apply it to study the induction stroke of a water analog model of a four-stroke internal combustion engine. The research conducted in the water analog model is a fundamental scientific inquiry into the flow fields that develop in the induction stroke of an engine at idling engine speeds. As this is the first investigation of its kind using LIPA technique, our goal has been to quantify, in a preliminary manner, the flow field features that develop during the intake stroke. In the process a more comprehensive understanding of the flow field features was developed, and tied to the quantification. The study evaluated the flow field of the intake stroke by estimating fields of velocity and vorticity. On the basis of these data, information about fluid dynamics during induction at engine speeds of 10, 20, and 30 RPM (corresponding to 170, 340, and 510 RPM respectively, when air is the flowing medium) for three different valve lifts was obtained. The overall development of the flow field, its energy content (kinetic, fluctuation) for the different settings of the engine parameters, vorticity information, and cyclic variations have been quantified. These have been discussed in terms of mixing performance.
The beginning of the 21st century is marked by a rapid growth of land surface satellite data and model sophistication. This offers new opportunities to estimate multiple components of the water cycle via satellite-based land data assimilation (DA) across multiple scales. By resolving more processes in land surface models and by coupling the land, the atmosphere, and other Earth system compartments, the observed information can be propagated to constrain additional unobserved variables. Furthermore, access to more satellite observations enables the direct constraint of more and more components of the water cycle that are of interest to end users. However, the finer level of detail in models and data is also often accompanied by an increase in dimensions, with more state variables, parameters, or boundary conditions to estimate, and more observations to assimilate. This requires advanced DA methods and efficient solutions. One solution is to target specific observations for assimilation based on a sensitivity study or coupling strength analysis, because not all observations are equally effective in improving subsequent forecasts of hydrological variables, weather, agricultural production, or hazards through DA. This paper offers a perspective on current and future land DA development, and suggestions to optimally exploit advances in observing and modeling systems.
This study shows that liquid mixtures of methanol and water are attractive candidates for microgravity droplet combustion experiments and associated numerical modeling. The gas phase chemistry for these droplet mixtures is conceptually simple, well understood and substantially validated. In addition, the thermodynamic and transport properties of the liquid mixture have also been well characterized. Furthermore, the results obtained in this study predict that the extinction of these droplets may be observable in ground-based drop to tower experiments. Such experiments will be conducted shortly followed by space-based experiments utilizing the NASA FSDC and DCE experiments.
The topics for the three lectures for the Canadian Summer School are Lagrangian Models, numerical transport schemes, and chemical and transport models. In the first lecture I will explain the basic components of the Lagrangian model (a trajectory code and a photochemical code), the difficulties in using such a model (initialization) and show some applications in interpretation of aircraft and satellite data. If time permits I will show some results concerning inverse modeling which is being used to evaluate sources of tropospheric pollutants. In the second lecture I will discuss one of the core components of any grid point model, the numerical transport scheme. I will explain the basics of shock capturing schemes, and performance criteria. I will include an example of the importance of horizontal resolution to polar processes. We have learned from NASA's global modeling initiative that horizontal resolution matters for predictions of the future evolution of the ozone hole. The numerical scheme will be evaluated using performance metrics based on satellite observations of long-lived tracers. The final lecture will discuss the evolution of chemical transport models over the last decade. Some of the problems with assimilated winds will be demonstrated, using satellite data to evaluate the simulations.
Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.
DNA nanotechnology provides a route towards the synthesis of custom nano-structured materials and circuits through hierarchical assembly processes. Although predictive kinetic models are being developed for the assembly of DNA nanostructures from small monomeric components, a general model for the hierarchical assembly of DNA nanostructures remains elusive. DNA tile nanotubes provide an ideal model system for the study of hierarchical assembly via end-to-end joining. In this study, we experimentally characterize the length-dependence of the end-to-end joining rate of DNA tile nanotubes. We then test the ability of three different models of polymer end-to-end joining to reproduce experimentally measured changes in nanotube lengths during a joining reaction using an ODE model for nanotube joining. All three models predict physically realistic joining rates that are consistent with prior measurements, with a length-independent end-to-end joining rate model providing the best fit to the experimental data. A length-independent constant joining rate is consistent with other DNA self-assembly processes across a broad range of length scales and also suggests how tractable models for hierarchical DNA nanostructure could be developed.
The existing theory of subset regression is examined, taking into account optimality criteria, small experiments, nonlinear models, colinearities, and special techniques. Approaches based on chain pooling coupled with principal components regression are discussed, giving attention to a comparison of half-normal plotting with chain pooling, a procedure based on prior ordering, deletion under the F-test, the largest of a set of chi-square variates, and principal components regression and model deletion. The choice of a true (population) model for simulations is considered along with the evaluation of the decision procedure and suitable computer programs.
A coal industry model was developed to support the Jet Propulsion Laboratory in its investigation of advanced underground coal extraction systems. The model documentation includes the programming for the coal mining cost models and an accompanying users' manual, and a guide to reading model output. The methodology used in assembling the transportation, demand, and coal reserve components of the model are also described. Results presented for 1986 and 2000, include projections of coal production patterns and marginal prices, differentiated by coal sulfur content.
A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.
The impact of numerical weather prediction (NWP) and satellite meteorology on operational weather forecasting has become overwhelming in the past few years. The paper looks toward the merger of these technologies in making short range 6-18 h forecasts through the use of mesoscale NWP models. A short-range (2-18 h) mesoscale forecast system envisioned for the near future is described that includes four components: hydrodynamic numerical models, large-scale and mesoscale; satellites, polar orbiting for high latitudes and geostationary for low latitudes; mesoclimatology, derived in large part from satellite data; and special-purpose simple models and empirical relations. It is important that the components of the forecast system be developed in parallel rather than in series if the system is to be completed within five years. There is enough evidence to substantiate the revolution in the mesoscale weather prediction in the next five years.
The space exposure of LDEF for 5.75 years, forming a host target in low earth orbit (LEO) orbit to a wide distribution of hypervelocity particulates of varying dimensions and different impact velocities, has yielded a multiplicity of impact features. Although the projectile parameters are generally unknown and, in fact not identical for any two impacts on a target, the great number of impacts provides statistically meaningful basis for the valid comparison of the response of different targets. Given sufficient impacts for example, a comparison of impact features (even without knowledge of the project parameters) is possible between: (1) differing material types (for the same incident projectile distribution); (2) differing target configurations (e.g., thick and thin targets for the same material projectiles; and (3) different velocities (using LDEF's different faces). A comparison between different materials is presented for infinite targets of aluminum, Teflon, and brass in the same pointing direction; the maximum finite-target penetration (ballistic limit) is also compared to that of the penetration of similar materials comprising of a semi-infinite target. For comparison of impacts on similar materials at different velocities, use is made of the pointing direction relative to LDEF's orbital motion. First, however, care must be exercised to separate the effect of spatial flux anisotropies from those resulting from the spacecraft velocity through a geocentrically referenced dust distribution. Data comprising thick and thin target impacts, impacts on different materials, and in different pointing directions is presented; hypervelocity impact parameters are derived. Results are also shown for flux modeling codes developed to decode the relative fluxes of Earth orbital and unbound interplanetary components intercepting LDEF. Modeling shows the west and space pointing faces are dominated by interplanetary particles and yields a mean velocity of 23.5 km/s at LDEF, corresponding to a V(infinity) Earth approach velocity = 20.9 km/s. Normally resolved average impact velocities on LDEF's cardinal point faces are shown. As 'excess' flux on the east, north, and south faces is observed, compatible with an Earth orbital component below some 5 microns in particle diameter.
Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.
To design failure management expert systems, engineers mentally analyze the effects of failures and procedures as they propagate through device configurations. CONFIG is a generic device modeling tool for use in discrete event simulation, to support such analyses. CONFIG permits graphical modeling of device configurations and qualitative specification of local operating modes of device components. Computation requirements are reduced by focussing the level of component description on operating modes and failure modes, and specifying qualitative ranges of variables relative to mode transition boundaries. Simulation processing occurs only when modes change or variables cross qualitative boundaries. Device models are built graphically, using components from libraries. Components are connected at ports by graphical relations that define data flow. The core of a component model is its state transition diagram, which specifies modes of operation and transitions among them.
Structured grid solvers such as NASA's OVERFLOW compressible Navier-Stokes flow solver can generate large data files that contain convergence histories for flow equation residuals, turbulence model equation residuals, component forces and moments, and component relative motion dynamics variables. Most of today's large-scale problems can extend to hundreds of grids, and over 100 million grid points. However, due to the lack of efficient tools, only a small fraction of information contained in these files is analyzed. OVERSMART (OVERFLOW Solution Monitoring And Reporting Tool) provides a comprehensive report of solution convergence of flow computations over large, complex grid systems. It produces a one-page executive summary of the behavior of flow equation residuals, turbulence model equation residuals, and component forces and moments. Under the automatic option, a matrix of commonly viewed plots such as residual histograms, composite residuals, sub-iteration bar graphs, and component forces and moments is automatically generated. Specific plots required by the user can also be prescribed via a command file or a graphical user interface. Output is directed to the user s computer screen and/or to an html file for archival purposes. The current implementation has been targeted for the OVERFLOW flow solver, which is used to obtain a flow solution on structured overset grids. The OVERSMART framework allows easy extension to other flow solvers.