Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated data”

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 1,081 records · Page 60

The Analysis of Image Segmentation Hierarchies with a Graph-based Knowledge Discovery System

Currently available pixel-based analysis techniques do not effectively extract the information content from the increasingly available high spatial resolution remotely sensed imagery data. A general consensus is that object-based image analysis (OBIA) is required to effectively analyze this type of data. OBIA is usually a two-stage process; image segmentation followed by an analysis of the segmented objects. We are exploring an approach to OBIA in which hierarchical image segmentations provided by the Recursive Hierarchical Segmentation (RHSEG) software developed at NASA GSFC are analyzed by the Subdue graph-based knowledge discovery system developed by a team at Washington State University. In this paper we discuss out initial approach to representing the RHSEG-produced hierarchical image segmentations in a graphical form understandable by Subdue, and provide results on real and simulated data. We also discuss planned improvements designed to more effectively and completely convey the hierarchical segmentation information to Subdue and to improve processing efficiency.

Tilton, James C.↗

The WFIRST CGI Integral Field Spectrograph: Requirements and Performance Predictions

The WFIRST coronagraphic instrument (CGI) will demonstrate exoplanet spectroscopy using an integral field spectrograph (IFS). The CGI IFS, being designed and built at Goddard Space Flight Center, has a spectral resolution of R50 and is designed to accommodate a 20% bandpass spanning 600- 970 nm. The IFS is principally targeting the abundance of Methane features, with the primary coronagraph band being centered around 770nm. Key to the performance estimates are the achievable signal-to-noise (SNR) ratios and the stability of the microspectra over the course of tens and hundreds of hours. As a technology demonstration for CGI, the ability to close a wavefront control loop around the IFS, maintain a stable dark hole, and provide time resolved data that simultaneously spans spatial and spectral dimensions are crucial demonstrations for future observatories. The IFS is optimized both for coronagraphs and science observations with a potential future starshade. We highlight how the long duration observations, and requirements for both starshades and coronagraphs drive the IFS requirements and the calibrations required both onorbit and on the ground. We also provide further detail on the optomechanical design, its stability based on thermal and structural predictions, anticipated performance, and operations concept of the CGI IFS. The impact of these performance metrics are projected into simulated data products, demonstrating cube extraction of noisy images and the subsequent planet spectrum that can be extracted from them. These demonstrations and performance predictions are key to future missions such as LUVOIR and HabEx, whose principal science case relies on efficient spectroscopy of exoplanets.

Tyler Groff↗

Subhourly Clipping Correction Model Comparison

This work will compare the Allen method and Walker method of accounting for subhourly inverter clipping power losses in hourly PV performance models. The Allen method uses a matrix lookup based on DNI clearness and clipping potential to assign a clipping correction loss at each simulation timestep. The Walker method models the PV DC power input to the inverter as a distribution over the hourly timestep and uses integration over the timestep to determine the amount of clipping that occurs within the timestep. Both these models have been recently implemented in the System Advisor Model's (SAM) open-source code, and will be applied to hourly SURFRAD datasets to analyze the subhourly clipping loss predicted by each model for different system designs and inverter loading conditions. Both models will be compared to "true" 1-minute SURFRAD data simulations to see their accuracy against more accurate 1-minute clipping correction loss predictions. This model comparisons will be investigated in more detail at the PVSC conference in Seattle, Washington June 2024.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗

Parallel Implementation of Nonadditive Gaussian Process Potentials for Monte Carlo Simulations

A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical simulation of the geographical sources of water for Continental Scale Experiments (CSEs) Precipitation

There are several important research questions that the Global Energy and Water Cycle Experiment (GEWEX) is actively pursuing, namely: What is the intensity of the water cycle and how does it change? And what is the sustainability of water resources? Much of the research to address these questions is directed at understanding the atmospheric water cycle. In this paper, we have used a new diagnostic tool, called Water Vapor Tracers (WVTs), to quantify the how much precipitation originated as continental or oceanic evaporation. This shows how long water can remain in the atmosphere and how far it can travel. The model-simulated data are analyzed over regions of interest to the GEWEX community, specifically, their Continental Scale Experiments (CSEs) that are in place in the United States, Europe, Asia, Brazil, Africa and Canada. The paper presents quantitative data on how much each continent and ocean on Earth supplies water for each CSE. Furthermore, the analysis also shows the seasonal variation of the water sources. For example, in the United States, summertime precipitation is dominated by continental (land surface) sources of water, while wintertime precipitation is dominated by the Pacific Ocean sources of water. We also analyze the residence time of water in the atmosphere. The new diagnostic shows a longer residence time for water (9.2 days) than more traditional estimates (7.5 days). We emphasize that the results are based on model simulations and they depend on the model s veracity. However, there are many potential uses for the new diagnostic tool in understanding weather processes and large and small scales.

Bosilovich, Michael G.↗

Ambient-Frequency-Data Based System-Level Inertia Estimation Using Physical Equation and its Practice on Hawaii Islands

Here, in this paper, a practical ambient-frequency-data-based inertia estimation method using a physical equation is proposed and validated by real measurement data from Hawaii island grids. With high renewable penetration, accurate inertia estimation is important and urgent. Ambient frequency oscillation always exists in power grids, so the proposed method has advantages of real-time inertia estimation and no need for additional disturbances. This paper first developed the physical equation for inertia estimation to offer a clear mechanism for easy implementation in practice. To apply the proposed method in actual grids, a practical method to extract the ambient frequency oscillation from frequency measurement is further proposed. The inertia estimation using the physical equation is validated by KIUC simulation data and HECO field data, in which error rates are around 2% and 8%, respectively. Practical inertia estimation is challenging due to the large amounts of resources contributing to the power grid's effective inertia, but the method provided in this paper can offer a novel way for practical inertia estimation, which can help renewable penetration to boost carbon-free grid.

14 SOLAR ENERGY↗

Slope-induced errors in radar altimetry over continental ice sheets

Altimeter measurements of ranges to modelled irregular surfaces are simulated and two correction schemes are used to reconstruct the modelled surfaces from the simulated data. When the cross-track slope is negligible, the problem can be treated as two-dimensional and the relocation method is preferred to the slope correction method. In a two-dimensional test case, 85 percent of the slope-induced rms error was removed by the relocation method. The mean error along profiles of about 75 km or longer is usually reduced more than the rms error. An alternative slope correction scheme, which uses the local slope to calculate the expected error, is less effective in the two-dimensional case. Over a simulated three-dimensional surface, where groundtracks are widely spaced and cross-track slopes are significant, the slope-correction method must be used in at least the cross-track direction.

Brenner, A. C.↗

Real-World Testing of LiDAR-Inertial Based Navigation and Mapping for Precision Landing

The fusion of LiDAR and inertial measurements duringspacecraft descent and landing can be used to estimate alander’s navigation state and map the terrain below. Together,these data products can be used to enable safe and preciselanding on celestial bodies for which a priori orbital reconnaissanceis insufficient for hazard detection and avoidance.Unlike camera images used in visual terrain relative navigation,LiDAR scans are insensitive to changes in illumination; as aresult, the technique can be used to land in poorly lit areas,or at times of day when the lighting conditions are incongruentwith existing orbital imagery. In this paper, we extend previouswork in which we introduced a factor graph based smoothingapproach for LiDAR-inertial navigation and mapping. Whereasthe algorithms were previously tested on simulated data, thispaper presents testing on real-world data. Data from theAutonomous Landing Hazard Avoidance Technology (ALHAT)airplane flight tests in the Yucca Flats and Death Valley in2009 (FT3), the Morpheus vertical take off and landing flighttests at Kennedy Space Center in 2014 (FT6), and the landingof Perseverance and Ingenuity on Mars in 2021 (M2020) wereused to evaluate algorithm performance. In this paper, weextend our LiDAR-inertial technique to work with a variety ofranging technologies: single point laser altimetry (FT3), denseflash LiDAR (FT6), and six-beam radar (M2020). A thoroughperformance analysis for all three datasets is presented. Datasetpreparation, improvements in algorithm robustness, and outlierrejection, which were necessitated by the transition to realworlddata, are discussed.

Trawny, Nikolas↗

A new method for recognizing quadric surfaces from range data and its application to telerobotics and automation

Pose and orientation of an object is one of the central issues in 3-D recognition problems. Most of today's available techniques require considerable pre-processing such as detecting edges or joints, fitting curves or surfaces to segment images, and trying to extract higher order features from the input images. We present a method based on analytical geometry, whereby all the rotation parameters of any quadric surface are determined and subsequently eliminated. This procedure is iterative in nature and was found to converge to the desired results in as few as three iterations. The approach enables us to position the quadric surface in a desired coordinate system, and then to utilize the presented shape information to explicitly represent and recognize the 3-D surface. Experiments were conducted with simulated data for objects such as hyperboloid of one and two sheets, elliptic and hyperbolic paraboloid, elliptic and hyperbolic cylinders, ellipsoids, and quadric cones. Real data of quadric cones and cylinders were also utilized. Both of these sets yielded excellent results.

Alvertos, Nicolas↗

Geo Thermal Cloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The project is motivated by the challenges, risks, and costs associated with geothermal exploration and production. Many processes and parameters impacting geothermal conditions are poorly understood. Diverse datasets are available to help characterize subsurface geothermal conditions (public and proprietary; satellite, airborne surveys, vegetation/water sampling, geological, geophysical, etc.). Yet, it is not clear how to properly leverage these datasets for geothermal exploration due to an incomplete understanding of how physical processes impacting subsurface geothermal conditions are represented in these observations. Recent advancements in machine learning (ML) provide great promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. Our goals and work under Phases 1 and 2 (as proposed) of this project address all these needs.

15 GEOTHERMAL ENERGY↗

Deconvoluting experimental decay energy spectra: The O 26 case

In nuclear reaction experiments, the measured decay energy spectra can give insights into the shell structure of decaying systems. However, extracting the underlying physics from the measurements is challenging due to detector resolution and acceptance effects. The Richardson-Lucy (RL) algorithm, a deblurring method that is commonly used in optics and has proven to be a successful technique for restoring images, was applied to our experimental nuclear physics data. The only inputs to the method are the observed energy spectrum and the detector's response matrix also known as the transfer matrix. We demonstrate that the technique can help access information about the shell structure of particle-unbound systems from the measured decay energy spectrum that is not immediately accessible via traditional approaches such as χ-square fitting. For a similar purpose, we developed a machine learning model that uses a deep neural network (DNN) classifier to identify resonance states from the measured decay energy spectrum. We tested the performance of both methods on simulated data and experimental measurements. Then, we applied both algorithms to the decay energy spectrum of 26 O → 24 O + n + n measured via invariant mass spectroscopy. Here, the resonance states restored using the RL algorithm to deblur the measured decay energy spectrum agree with those found by the DNN classifier. Both deblurring and DNN approaches suggest that the raw decay energy spectrum of 26 O exhibits three peaks at approximately 0.15 MeV, 1.50 MeV, and 5.00 MeV, with half-widths of 0.29 MeV, 0.80 MeV, and 1.85 MeV, respectively.

Spectrometers & spectroscopic techniques↗

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Model-based Reconstruction for Single Particle Cryo-Electron Microscopy

Single particle cryo-electron microscopy is a vital tool for 3D characterization of protein structures. A typical workflow involves acquiring projection images of a collection of randomly oriented particles, picking and classifying individual particle projections by orientation, and finally using the individual particle projections to reconstruct a 3D map of the electron density profile. The reconstruction is challenging because of the low signal-to-noise ratio of the data, the unknown orientation of the particles, and the sparsity of data especially when dealing with flexible proteins where there may not be sufficient data corresponding to each class to obtain an accurate reconstruction using standard algorithms. In this paper we present a model-based image reconstruction technique that uses a regularized cost function to reconstruct the 3D density map by assuming known orientations for the particles. Our method casts the reconstruction as minimizing a cost function involving a novel forward model term that accounts for the contrast transfer function of the microscope, the orientation of the particles and the center of rotation offsets. We combine the forward model term with a regularizer that enforces desirable properties in the volume to be reconstructed. Using simulated data, we demonstrate how our method can significantly improve upon the typically used approach.

Venkatakrishnan, Singanallur↗

A review of near-wall Reynolds-stress

The advances made in second-order near-wall turbulence closures are summarized. All closures examined are based on some form of high Reynolds number models for the Reynolds stress and the turbulent kinetic energy dissipation rate equations. Consequently, most near-wall closures proposed to data attempt to modify the high Reynolds number models for the dissipation rate equation so that the resultant models are applicable all the way to the wall. The near-wall closures are examined for their asymptotic behavior so that they can be compared with the proper near-wall behavior of the exact equations. A comparison of the closure's performance in the calculation of a low Reynolds number plane channel flow is carried out. In addition, the closures are evaluated for their ability to predict the turbulence statistics and the limiting behavior of the structure parameters compared to direct simulation data.

So, R. M. C.↗

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗