Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Synthetic code”

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 235 records · Page 13

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA's GMAO Atmospheric Motion Vectors Simulator: Description and Application to the MISTiC Winds Concept

An atmospheric wind vectors (AMVs) simulator was developed by NASA's GMAO to simulate observations from future satellite constellation concepts. The synthetic AMVs can then be used in OSSEs to estimate and quantify the potential added value of new observations to the present Earth observing system and, ultimately, the expected impact on the current weather forecasting skill. The GMAO AMV simulator is a tunable and flexible computer code that is able to simulate AMVs expected to be derived from different instruments and satellite orbit configurations. As a case study and example of the usefulness of this tool, the GMAO AMV simulator was used to simulate AMVs envisioned to be provided by the MISTiC Winds, a NASA mission concept consisting of a constellation of satellites equipped with infrared spectral midwave spectrometers, expected to provide high spatial and temporal resolution temperature and humidity soundings of the troposphere that can be used to derive AMVs from the tracking of clouds and water vapor features. The GMAO AMV simulator identifies trackable clouds and water vapor features in the G5NR and employs a probabilistic function to draw a subset of the identified trackable features. Before the simulator is applied to the MISTiC Winds concept, the simulator was calibrated to yield realistic observations counts and spatial distributions and validated considering as a proxy instrument to the MISTiC Winds the Himawari-8 Advanced Imager (AHI). The simulated AHI AMVs showed a close match with the real AHI AMVs in terms of observation counts and spatial distributions, showing that the GMAO AMVs simulator synthesizes AMVs observations with enough quality and realism to produce a response from the DAS equivalent to the one produced with real observations. When applied to the MISTiC Winds scanning points, it can be expected that the MISTiC Winds will be able to collect approximately 60,000 wind observations every 6 hours, if considering a constellation composed of 12 satellites (4 orbital planes). In addition, one of the main expected impacts of the MISTiC Winds concept is the ability to derive water vapor feature tracking AMVs below 500-400 hPa, an unique feature among the water vapor AMVs derived from the current Earth observing system.

Carvalho, David↗

NASA's GMAO Atmospheric Motion Vectors Simulator: Description and Application to the MISTiC Winds Concept

An atmospheric wind vectors (AMVs) simulator was developed by NASA's GMAO to simulate observations from future satellite constellation concepts. The synthetic AMVs can then be used in OSSEs to estimate and quantify the potential added value of new observations to the present Earth observing system and, ultimately, the expected impact on the current weather forecasting skill. The GMAO AMV simulator is a tunable and flexible computer code that is able to simulate AMVs expected to be derived from different instruments and satellite orbit configurations. As a case study and example of the usefulness of this tool, the GMAO AMV simulator was used to simulate AMVs envisioned to be provided by the MISTiC Winds, a NASA mission concept consisting of a constellation of satellites equipped with infrared spectral midwave spectrometers, expected to provide high spatial and temporal resolution temperature and humidity soundings of the troposphere that can be used to derive AMVs from the tracking of clouds and water vapor features. The GMAO AMV simulator identifies trackable clouds and water vapor features in the G5NR and employs a probabilistic function to draw a subset of the identified trackable features. Before the simulator is applied to the MISTiC Winds concept, the simulator was calibrated to yield realistic observations counts and spatial distributions and validated considering as a proxy instrument to the MISTiC Winds the Himawari-8 Advanced Imager (AHI). The simulated AHI AMVs showed a close match with the real AHI AMVs in terms of observation counts and spatial distributions, showing that the GMAO AMVs simulator synthesizes AMVs observations with enough quality and realism to produce a response from the DAS equivalent to the one produced with real observations. When applied to the MISTiC Winds scanning points, it can be expected that the MISTiC Winds will be able to collect approximately 60,000 wind observations every 6 hours, if considering a constellation composed of 12 satellites (4 orbital planes). In addition, one of the main expected impacts of the MISTiC Winds concept is the ability to derive water vapor feature tracking AMVs below 500-400 hPa, an unique feature among the water vapor AMVs derived from the current Earth observing system.

Carvalho, David↗

Review of Recent Results in Global MHD Modeling: ISTP Project Scientist for Theory and Ground-Based Observations

Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.

Curtis, Steven↗

MODIS Science Team Member Semi-annual Report

This paper presents a semi-annual report of the MODerate resolution imaging Spectroradiometer (MODIS) Science Team Members. The most important activities undertaken during this reporting period are the following: 1) Versions 2.1 and 2.2 surface reflectance L2/L3 DAAC/SDST delivery; 2) Version 2.0 1km and 250m VI product delivery (assist Arizona); 3) Version 2.1 surface reflectance L2 testing; 4) Land Synthetic data set generator improvements; 5) QA; 6) Surface reflectance error budget generation (SWAMP request); 7) SCF Hardware; 8) Aerosol transport modeling; 9) Aerosol optical depth retrieval from AVHRR data; 10) Aerosol characteristics retrieval from SeaWIFS/AVHRR fusioned data; 11) Validation activities; 12) Aerosol climatology; and 13) 6S code. The report includes summaries of the topics above.

Vermote, Eric↗

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Game-Based Learning Theory

Persistent Immersive Synthetic Environments (PISE) are not just connection points, they are meeting places. They are the new public squares, village centers, malt shops, malls and pubs all rolled into one. They come with a sense of 'thereness" that engages the mind like a real place does. Learning starts as a real code. The code defines "objects." The objects exist in computer space, known as the "grid." The objects and space combine to create a "place." A "world" is created, Before long, the grid and code becomes obscure, and the "world maintains focus.

Laughlin, Daniel↗

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗

Assessing Edge Localized Mode (ELM) Mitigation by Pellet Triggering in Smallest-Possible Affordable Robust Compact (SPARC) Low-Collisionality Discharges

Simulations of pellet edge-localized mode (ELM) triggering with the M3D-C1 code confirmed experimentally observed pellet mass thresholds in low-collisionality DIII-D discharges and predicted extremely large thresholds in ITER [A. Wingen et al, Nucl. Fusion 64, 066015(2024)]. The linear (static) and nonlinear (dynamic) workflows developed in this previous work are now applied to the SPARC tokamak. Multiple synthetic equilibria for the high-performance H-mode scenario are tested for MHD stability without pellets using the linear workflow and ELITE simulations; a stable one is identified. Pellets of various mass are injected along an inboard, near-midplane trajectory in the M3D-C1 simulation. No pellet mass threshold for ELM triggering can be found for realistic pellet masses, concluding that fueling from this injection location would not trigger ELMs. A poloidal scan of possible injection locations shows very high thresholds on the high-field side and reasonable thresholds on the low-field-side (LFS). Based on results of this work the location for a LFS injection was finalized in the SPARC design. For this LFS injection the pellet mass threshold is determined using the linear workflow, aligning well with the poloidal scan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiative transfer in the dynamic atmospheres of long period variable stars

An iterative procedure is presented for determining the thermal structure and dynamics of Mira-type stellar atmospheres, where the non-LTE radiative transfer code PANDORA is used in conjunction with the Bowen hydrodynamics code of Iowa State University. Preliminary results are reported for an atmospheric model of a pulsating AGB star of 1 solar mass, 240 solar radii, Teff = 3000 K, and a period of 320 days. At the present time, H, H(-), Mg I, and Mg II radiative transfer calculations have been completed and synthetic spectra are shown for H-alpha. The radiative transfer calculations demonstrate that cooling in the innermost shock of the original Bowen model is underestimated due to the omission of various hydrogen transitions. These initial results suggest that the main shock of the Bowen models are too hot and/or too deep.

Luttermoser, Donald G.↗

Active Sonic Boom Control

A theory and simulation code are developed to study nonsteady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of flight path. An aircraft in nonsteady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender nonsteady sources but requires no far field approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two nonsteady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and convert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on the average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical

Crow, Steven C.↗

Robust reconstruction of single-cell RNA-seq data with iterative gene weight updates

Single-cell RNA-sequencing technologies have greatly enhanced our understanding of heterogeneous cell populations and underlying regulatory processes. However, structural (spatial or temporal) relations between cells are lost during cell dissociation. These relations are crucial for identifying associated biological processes. Many existing tissue-reconstruction algorithms use prior information about subsets of genes that are informative with respect to the structure or process to be reconstructed. When such information is not available, and in the general case when the input genes code for multiple processes, including being susceptible to noise, biological reconstruction is often computationally challenging. We propose an algorithm that iteratively identifies manifold-informative genes using existing reconstruction algorithms for single-cell RNA-seq data as subroutine. We show that our algorithm improves the quality of tissue reconstruction for diverse synthetic and real scRNA-seq data, including data from the mammalian intestinal epithelium and liver lobules.

59 BASIC BIOLOGICAL SCIENCES↗

Central America Disasters: Using Earth Observations to Map Flooding for Disaster Monitoring, Inform Potential Risk, and Prepare for Possible Response

In November 2020, Hurricanes Eta and Iota hit Central America within weeks of each other, causing severe flooding, landslides, and widespread damage. NASA DEVELOP partnered with Comité Regional de Recursos Hidráulicos (CRRH), Centro de Coordinación para la Prevención de los Desastres en América Central y República Dominicana (CEPREDENAC), and Sistema de la Integración Centroamericana (SICA) to better understand how flooding throughout Central America has impacted and will continue to affect communities, focusing on sites in Guatemala, Honduras, El Salvador, Belize, Nicaragua, western Panama, and eastern Costa Rica from January 2015 to October 2021. The team utilized surface reflectance data from Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI), Sentinel-2 MultiSpectral Instrument (MSI), Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS), and Terra Moderate Resolution Imaging Spectroradiometer (MODIS). This project also utilized backscatter data from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and elevation data from the Shuttle Radar Topography Mission (SRTM). Incorporating these Earth observations in NASA SERVIR’s Hydrologic Remote Sensing Analysis for Floods (HYDRAFloods) tool run on Google Earth Engine (GEE), the team produced historical surface water maps, a case study analysis of the two hurricanes, and a code tutorial. These results indicated that surface water increased in priority sites from 2015 to 2021, optical and SAR imagery detected similar flood patterns and extent after the hurricanes, and rainfall was concentrated on the east coast of the region. These products allow partners to make informed decisions around flooding preparation and disaster mitigation.

Caroline Williams↗

Machine Learning-based Prediction of Departure from Nucleate Boiling Power for the PSBT Benchmark

Machine Learning (ML) has seen an exponential growth in its applications due to its advanced data driven prediction capabilities. The study presents a data-driven approach as a preliminary attempt to predict the power at which departure from nucleate boiling (DNB) occurs in pressurized water reactors (PWRs) by constructing an advanced ML algorithm that takes outlet pressure, inlet temperature and inlet mass flux as the input features. DNB is a critical heat flux (CHF) phenomenon seen in PWRs. The experimental data from the PWR subchannel and bundle tests (PSBT) benchmark is first used to train an artificial neural network (ANN) to predict the DNB power, which produces a root mean square error (RMSE) of 6.89 kW/m when tested on a blind subset of the PSBT data. Since the PSBT dataset is relatively small to train an accurate ANN, a data augmentation methodology based on generative adversarial networks (GANs) is used to expand the training dataset. By assuming that the real data follows a certain distribution, GANs try to learn that underlying distribution to generate similar synthetic data to augment the database and to improve the predictive capabilities of the ANN. The data generated from GANs are validated using 1-nearest neighbor and kernel maximum mean discrepancy. To further ensure data from GAN is similar to PSBT, the data is tested and filtered out using the sub-channel thermal-hydraulic code CTF. The results indicate that with the addition of 120 data points from GAN the RMSE reduces to 4.84 kW/m showing promising results for future developments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

GOOML Big Kahuna Forecast Modeling and Genetic Optimization Files

This submission includes example files associated with the Geothermal Operational Optimization using Machine Learning (GOOML) Big Kahuna fictional power plant, which uses synthetic data to model a fictional power plant. A forecast was produced using the GOOML data model framework and fictional input data, and a genetic optimization is included which determines optimal flash plant parameters. The inputs and outputs associated with the forecast and genetic optimization are included. The input and output files consist of data, configuration files, and plots. A link to the Physics-Guided Neural Networks (phygnn) GitHub repository is also included, which augments a traditional neural network loss function with a generic loss term that can be used to guide the neural network to learn physical or theoretical constraints. phygnn is used by the GOOML framework to help integrate its machine learning models into the relevant physics and engineering applications. Note that the data included in this submission are intended to provide a demonstration of GOOML's capabilities. Additional files that have not been released to the public are needed for users to run these models and reproduce these results. Units can be found in the readme data resource.

15 GEOTHERMAL ENERGY↗

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗

Analysis of the MODIS Above-Cloud Aerosol Retrieval Algorithm Using MCARS

The Multi-sensor Cloud and Aerosol Retrieval Simulator (MCARS) presently produces synthetic radiance data from Goddard Earth Observing System version 5 (GEOS-5) model output as if the Moderate Resolution Imaging Spectroradiometer (MODIS) was viewing a combination of atmospheric column inclusive of clouds, aerosols and a variety of gases and land/ocean surface at a specific location. In this paper we use MCARS to study the MODIS Above-Cloud AEROsol retrieval algorithm (MOD06ACAERO). MOD06ACAERO is presently a regional research algorithm able to retrieve aerosol optical thickness over clouds, in particular absorbing biomass burning aerosols overlying marine boundary layer clouds in the Southeastern Atlantic Ocean. The algorithm’s ability to provide aerosol information in cloudy conditions makes it a valuable source of information for modeling and climate studies n an area where current clear sky-only operational MODIS aerosol retrievals effectively have a data gap between the months of June and October. We use MCARS for a verification and closure study of the MOD06ACAERO algorithm. The purpose of this study is to develop a set of constraints a model developer might use during assimilation of MOD06ACAERO data. Our simulations indicate that the MOD06ACAERO algorithm performs well for marine boundary layer clouds in the SE Atlantic provided some specific screening rules are observed. For the present study, a combination of five simulated MODIS data granules was used for a dataset of 13.5 million samples with known input conditions. When pixel retrieval uncertainty was less than 30%, optical thickness of the underlying cloud layer was greater than 4 and scattering angle range within the cloud bow was excluded, MOD06ACAERO retrievals agreed with the underlying ground truth (GEOS-5 cloud and aerosol profiles used to generate the synthetic radiances) with a slope of 0.913, offset of 0.06, and RMSE=0.107. When only near-nadir pixels were considered (view zenith angle within +/-20 degrees) the agreement with source data further improved (0.977, 0.051 and 0.096 respectively). Algorithm closure was examined using a single case out of the five 38 used for verification. For closure, the MOD06ACAERO code was modified to use GEOS-5 temperature and moisture profiles as ancillary. Agreement of MOD06ACAERO retrievals with source data for the closure study had a slope of 0.996 with offset -0.007 and RMSE of 0.097 at pixel uncertainty level of less than 40%, illustrating the benefits of high-quality ancillary atmospheric data for such retrievals.

MODIS↗

SDSS-IV MaStar: [α/Fe] for the MaNGA Stellar Library from synthetic model spectra

ABSTRACT We calculate the α-enhancement ratio [α/Fe] for the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar) while also fitting for the fundamental atmospheric parameters effective temperature, surface gravity, and metallicity – Teff, log g, [Fe/H]. This approach builds upon a previous catalogue of stellar parameters, whereby only the fundamental atmospheric parameters are fit with solar-scaled models. Here, we use the same Markov Chain Monte Carlo method with the additional free parameter [α/Fe]. Using the full spectral fitting code pPXF, we are able to fit multiple lines sensitive to [α/Fe] for a more robust measurement. Quality flags based on the convergence of the sampler, errors in [α/Fe] and a cut in the χ2 of the model fit are used to clean the final catalogue, returning 17 214 spectra and values in the range of −0.25 < [α/Fe] < 0.48. Comparing our calculated [α/Fe] with literature values reveals a degeneracy in cool stars with log g ≥ ∼4; this comparison is then used to create an alternative and calibrated parameter set. We also plot the final catalogue in an [Fe/H] versus [α/Fe] diagram and recover the expected result of increasing [α/Fe] with decreasing [Fe/H] for Milky Way disc-halo stars. We apply our method to a subsample of spectra of uniform resolution and higher signal to noise that finds that our results are independent of this higher signal to noise. In the context of stellar population models, we are able to cover a parameter space for the creation of intermediate to old age models at solar-scaled [α/Fe], high [Fe/H] and enhanced [α/Fe], low [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗