Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation 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 685 records · Page 38

Evaluating Clouds in Long-Term Cloud-Resolving Model Simulations with Observational Data

Two 20-day, continental midlatitude cases are simulated with a three-dimensional (3D) cloud-resolving model (CRM) and compared to Atmospheric Radiation Measurement (ARM) data. This evaluation of long-term cloud-resolving model simulations focuses on the evaluation of clouds and surface fluxes. All numerical experiments, as compared to observations, simulate surface precipitation well but over-predict clouds, especially in the upper troposphere. The sensitivity of cloud properties to dimensionality and other factors is studied to isolate the origins of the over prediction of clouds. Due to the difference in buoyancy damping between 2D and 3D models, surface precipitation fluctuates rapidly with time, and spurious dehumidification occurs near the tropopause in the 2D CRM. Surface fluxes from a land data assimilation system are compared with ARM observations. They are used in place of the ARM surface fluxes to test the sensitivity of simulated clouds to surface fluxes. Summertime simulations show that surface fluxes from the assimilation system bring about a better simulation of diurnal cloud variation in the lower troposphere.

Zeng, Xiping↗

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗

Cloud model-based simulations of satellite microwave data and their application to an SSM/I rain algorithm

Satellite microwave observations of convective rain systems are simulated using a three-dimensional cloud model and a microwave radiative-transfer model. The simulated data sets are then used as a basis to derive rain estimation algorithms which are tested on actual satellite data. A tropical simulation is discussed, and the algorithms are applied to SSM/I data.

Adler, Robert F.↗

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

Cooldown Strategy for Cryogenic Wind Tunnels

Criteria substantiated by experimental and real-time simulation data ensure optimal utlization of liquid nitrogen injected for cooling a 0.3-m transonic cryogenic tunnel (TCT). Research includes simulated cooldowns using procedures normally employed by each of several TCT operators, simulated cooldowns at constant mach number, simulated cooldowns at various constant fan speeds and cooldowns at constant metal-to-gas temperature differences.

Thibodeaux, J. J.↗

Atmospheric effects on remote sensing of surface reflectance

This paper reviews the atmospheric effects on remote sensing of surface reflectance. The scattering and absorption of sunlight by atmospheric molecules and aerosols affects the quality of images of the surface remotely sensed from satellites and aircrafts. The concentration and characteristics of the atmospheric aerosols vary from place to place and vary with time. The effect of atmospheric aerosols on the upward radiance depends on their optical thickness, scattering phase function and absorption. These parameters result from the aerosol concentration, composition, and the relative humidity. For high resolution images the aerosol scale height is also of importance. The radiative transfer theory that predicts the atmospheric radiances for a given surface and atmosphere is a well established theory for the case of uniform surfaces (or low resolution data). Some radiative transfer models exist for nonuniform surfaces and others are being developed. Recent field experiment and laboratory simulation data confirm the need for these models and can be used for their testing. It is shown that the atmospheric effect reduces the apparent resolution of satellite imagery and causes errors in the classification of surface fields. Suggestions for correction procedures are given. Such corrections can be based on ground observations, on satellite radiances above dark areas, or on climatologic information, depending on the accuracy of the corrections needed. The chosen correction algorithm depends also on the image resolution and the specific remote sensing application.

Kaufman, Y. J.↗

Isotropy of the small scales of turbulence at low Reynolds number

Spectral local isotropy tests are applied to direct numerical simulation data, mainly at the centerline of a fully developed turbulent channel flow. Despite the small Reynolds number of the simulation, the high-wavenumber behavior of velocity and vorticity spectra is consistent with local isotropy. This consistency is verified by the relationship between streamwise wavenumber spectra and spanwise wavenumber spectra. The high-wavenumber behavior of the pressure spectrum is also consistent with local isotropy and compares favorably with the calculation of Batchelor (1951), which assumes isotropy and joint normality of the velocity field at two points in space. The latter assumption is validated by the shape but not the magnitude of the quadruple correlation of the streamwise velocity fluctuation at small separations. There is only partial support for local spectral isotropy away from the centerline as the magnitude of the mean strain rate increases.

Kim, J.↗

Atmospheric infrared sounder

This report summarizes the activities of two Atmospheric Infrared Sounder (AIRS) team members during the first half of 1995. Changes to the microwave first-guess algorithm have separated processing of Advanced Microwave Sounding Unit A (AMSU-A) from AMSU-B data so that the different spatial resolutions of the two instruments may eventually be considered. Two-layer cloud simulation data was processed with this algorithm. The retrieved water vapor column densities and liquid water are compared. The information content of AIRS data was applied to AMSU temperature profile retrievals in clear and cloudy atmospheres. The significance of this study for AIRS/AMSU processing lies in the improvement attributable to spatial averaging and in the good results obtained with a very simple algorithm when all of the channels are used. Uncertainty about the availability of either a Microwave Humidity Sensor (MHS) or AMSU-B for EOS has motivated consideration of possible low-cost alternative designs for a microwave humidity sensor. One possible configuration would have two local oscillators (compared to three for MHS) at 118.75 and 183.31 GHz. Retrieval performances of the two instruments were compared in a memorandum titled 'Comparative Analysis of Alternative MHS Configurations', which is attached.

Rosenkranz, Philip, W.↗

Comparing Aircraft Observations of Snowfall to Forecasts Using Single or Two Moment Bulk Water Microphysics Schemes

High resolution weather forecast models with explicit prediction of hydrometeor type, size distribution, and fall speed may be useful in the development of precipitation retrievals, by providing representative characteristics of frozen hydrometeors. Several single or double-moment microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, allowing for the prediction of up to three ice species. Each scheme incorporates different assumptions regarding the characteristics of their ice classes, particularly in terms of size distribution, density, and fall speed. In addition to the prediction of hydrometeor content, these schemes must accurately represent the vertical profile of water vapor to account for possible attenuation, along with the size distribution, density, and shape characteristics of ice crystals that are relevant to microwave scattering. An evaluation of a particular scheme requires the availability of field campaign measurements. The Canadian CloudSat/CALIPSO Validation Project (C3VP) obtained measurements of ice crystal shapes, size distributions, fall speeds, and precipitation during several intensive observation periods. In this study, C3VP observations obtained during the 22 January 2007 synoptic-scale snowfall event are compared against WRF model output, based upon forecasts using four single-moment and two double-moment schemes available as of version 3.1. Schemes are compared against aircraft observations by examining differences in size distribution, density, and content. In addition to direct measurements from aircraft probes, simulated precipitation can also be converted to equivalent, remotely sensed characteristics through the use of the NASA Goddard Satellite Data Simulator Unit. Outputs from high resolution forecasts are compared against radar and satellite observations emphasizing differences in assumed crystal shape and size distribution characteristics.

Molthan, Andrew L.↗

Software Modules for the Proximity-1 Space Link Interleaved Time Synchronization (PITS) Protocol

The Proximity-1 Space Link Interleaved Time Synchronization (PITS) protocol provides time distribution and synchronization services for space systems. A software prototype implementation of the PITS algorithm has been developed that also provides the test harness to evaluate the key functionalities of PITS with simulated data source and sink. PITS integrates time synchronization functionality into the link layer of the CCSDS Proximity-1 Space Link Protocol. The software prototype implements the network packet format, data structures, and transmit- and receive-timestamp function for a time server and a client. The software also simulates the transmit and receive-time stamp exchanges via UDP (User Datagram Protocol) socket between a time server and a time client, and produces relative time offsets and delay estimates.

Woo, Simon S.↗

Developing a Spectral Correlation based Method for mitigating Angular Mismatch Effects in Intercalibration Using a Benchmark Hyperspectral Sensor

The CLARREO Pathfinder (CPF) mission will implement a state-of-the-art intercalibration method for transferring CPF’s in-orbit Système Internationale (SI)-traceable reference to the shortwave channel of the Clouds and the Earth’s Radiance Energy System (CERES) and the reflective solar bands of the Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the NOAA-20 satellite platform with a targeted intercalibration methodology uncertainty of 0.3% (k=1). In order to achieve such a high intercalibration accuracy, the CPF intercalibration measurements to be scheduled need to closely match those from CERES and VIIRS in time, space, angles, and wavelength. Despite’s CPF’s two-axis pointing capability to match its boresight line of sight to that of a target sensor, there will be finite residual differences in angular samplings of the two instruments. The potential angular mismatch between CPF and the target sensors can introduce systematic errors in the intercalibration results, and therefore needs to be corrected using a rigorously designed algorithm. The CPF intercalibration team has been developing a correction method for mitigating the impact of these angular anisotropic effects in the CPF-VIIRS and CPF-CERES intercalibration. The method explores the spectral correlation relationship between the reflected solar radiances from the same surface target that would be measured by CPF at two adjacent angles. Our studies have shown that the spectral information based on CPF measurements can be used to accurately predict the spectral radiance or reflectance difference due to a given mismatch in the viewing and solar geometry. The hyper-spectral information can also be used to provide scene stratification without using any auxiliary data, which is critical to reduce the angular correction uncertainty. The angular correction relationship can be well established using simulated CPF-like top-of-atmosphere spectral radiances observed at all sorts of viewing geometry and solar angles and for different scenes. Intensive radiative transfer simulations have been completed using a state-of-art radiative transfer model developed by the CPF team members. The implementation of the algorithm on high-fidelity event simulation data and the characterization for the angular adjustment uncertainty will be presented.

Wan Wu↗

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

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) 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. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

A Gaussian process based surrogate approach for the optimization of cylindrical targets

Simulating direct-drive inertial confinement experiments presents significant computational challenges, both due to the complexity of the codes required for such simulations and the substantial computational expense associated with target design studies. Machine learning models, and in particular, surrogate models, offer a solution by replacing simulation results with a simplified approximation. In this study, we apply surrogate modeling and optimization techniques that are well established in the existing literature to one-dimensional simulation data of a new cylindrical target design containing deuterium–tritium fuel. These models predict yields without the need for expensive simulations. We find that Bayesian optimization with Gaussian process surrogates enhances sampling efficiency in low-dimensional design spaces but becomes less efficient as dimensionality increases. Nonetheless, optimization routines within two-dimensional and five-dimensional design spaces can identify designs that maximize yield, while also aligning with established physical intuition. Optimization routines, which ignore constraints on hydrodynamic instability growth, are shown to lead to unstable designs in 2D, resulting in yield loss. However, routines that utilize 1D simulations and impose constraints on the in-flight aspect ratio converge on novel cylindrical target designs that are stable against hydrodynamic instability growth in 2D and achieve high yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

September 2023 - Mini Blast Tube Experimental Results [Slides]

Ryan Rettinger (M-6): provided the python script for extracting the pressure data. Joseph Kerwin (E-14): provided the countless pictures taken during the flow shots. Chris Long (W-13): provided the simulation data shown in the upcoming slides.

42 ENGINEERING↗

Technique to determine location of radio sources from measurements taken on spinning spacecraft

The procedure developed to extract average source direction and average source size from spin-modulated radio astronomy data measured on the IMP-6 spacecraft is described. Because all measurements are used, rather than just finding maxima or minima in the data, the method is very sensitive, even in the presence of large amounts of noise. The technique is applicable to all experiments with directivity characteristics. It is suitable for onboard processing on satellites to reduce the data flow to Earth. The application to spin-modulated nonpolarized radio astronomy data is made and includes the effects of noise, background, and second source interference. The analysis was tested with computer simulated data and the results agree with analytic predictions. Applications of this method with IMP-6 radio data have led to: (1) determination of source positions of traveling solar radio bursts at large distances from the Sun; (2) mapping of magnetospheric radio emissions by radio triangulation; and (3) detection of low frequency radio emissions from Jupiter and Saturn.

Fainberg, J.↗

Novel Turbulence Model Tailored for Complex Rocket Motor Exhaust Jets

This paper describes a new RANS turbulence model intended to give improved accuracy in turbulence that is affected by compressibility, rotation and curvature, and variable density, making it well-suited to complex rocket motor exhaust jets. Four key elements of the model are: the use of simple-averaged (not mass-averaged) turbulence decomposition; a second-order Taylor series for the equilibrated turbulent fluxes; empirical modeling that accounts for the combined effects of aforementioned turbulence modifiers; and model tuning on large-eddy simulation data for mixing layers for a wide range of these turbulence modifiers. The model is implemented as an extension of the k-w SST model. Results include comparisons to detached eddy simulation results for hot and cold co-flowing exhaust jets, and to experimental data for supersonic retropropulsion on the CobraMRV Mars-lander concept vehicle. The model is named TMRC for its constituents elements, the Taylor series with Mach number, rotation and curvature effects.

turbulence modeling↗

Inference finds consistency between a neutrino flavor evolution model and Earth-based solar neutrino measurements

We continue examining statistical data assimilation (SDA), an inference methodology, to infer solutions to neutrino flavor evolution, for the first time using real - rather than simulated - data. The model represents neutrinos streaming from the Sun's center and undergoing a Mikheyev-Smirnov-Wolfenstein (MSW) resonance in flavor space, due to the radially-varying electron number density. The model neutrino energies are chosen to correspond to experimental bins in the Sudbury Neutrino Observatory (SNO) and Borexino experiments, which measure electron-flavor survival probability at Earth. In conclusion, the procedure successfully finds consistency between the observed fluxes and the model, if the MSW resonance - that is, flavor evolution due to solar electrons - is included in the dynamical equations representing the model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Data Analysis and Simulation Study of Urban Air Mobility

For the Urban Air Mobility (UAM) industry, NASA has defined a series of UAM Maturity Levels (UML) corresponding to increasingly more complex and operationally dense UAM operations. In support of the gradual progression towards higher UML levels, NASA is currently conducting a set of UAM air traffic simulations—collectively referred to as X4. This paper describes a set of system effectiveness measures, and their associated metrics, for data analysis of X4 simulations. The descriptions, rationales, and calculation procedures for two metrics to be used in data analysis of simulation results, the number of predicted demand-capacity imbalances and the pre-departure delays, are described. Results from data analysis of one set of simulation runs are presented to demonstrate how these metrics support the assessment of performance of the system architecture for X4 simulations and the verification of experiment requirements.

Urban Air Mobility↗