Search NASA⌕ Search

SEARCH · Search NASA

Results for “synthetic model”

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 217 records · Page 12

Investigating laser ablated plume dynamics of carbon and aluminum targets

Recently acquired high-resolution images of nanosecond laser ablation plumes suggest a strong correlation between the internal plume structure and the type of material being ablated. However, the details of this relation are currently not well understood. In this work, we attempt to explore this correlation using a 2D radiation hydrodynamics model to study the dependence of internal plume structure formation on the ablation material. Here spatio-temporal emission maps and plume expansion velocities from experimental measurements are compared with the model predictions, including synthetic emission maps. The shape and expansion rate of an outer air plume region are found to be in good agreement for both carbon and aluminum, as are the inner material plume dynamics for carbon ablation. The largest disagreement is observed in the case of a polished aluminum target, where the chaotic inner plume features seen in the experimental images are not observed in the model. The possible physical mechanisms responsible for this discrepancy are discussed. This effort constitutes a continued development toward a predictive model of ablation plume dynamics and chemistry for various materials in extreme environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A hybrid deep neural operator/finite element method for ice-sheet modeling

One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work we develop a hybrid approach to approximate existing ice sheet models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while we retain a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model, compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier and can therefore be used for generating samples for uncertainty quantification. Further, we study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗

Maven: a multimodal foundation model for supernova science

Abstract A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. At the same time, no data-driven models exist to understand these photometric and spectroscopic observables in a common context. Contrastive learning objectives, which have grown in popularity for aligning distinct data modalities in a shared embedding space, provide a potential solution to extract information from these modalities. We present Maven, the first foundation model for supernova science. To construct Maven, we first pre-train our model to align photometry and spectroscopy from 0.5 M synthetic supernovae using a contrastive objective. We then fine-tune the model on 4702 observed supernovae from the Zwicky transient facility. Maven reaches state-of-the-art performance on both classification and redshift estimation, despite the embeddings not being explicitly optimized for these tasks. Through ablation studies, we show that pre-training with synthetic data improves overall performance. In the upcoming era of the Vera C. Rubin observatory, Maven will serve as a valuable tool for leveraging large, unlabeled and multimodal time-domain datasets.

Zhang, Gemma (ORCID:0000000280198082)↗

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗

Virus-Like Particles (VLPs) as a Platform for Hierarchical Compartmentalization

Hierarchically self-assembled structures are common in biology, but it is often challenging to design and fabricate synthetic analogs. The archetypal cell is defined by hierarchically organized multicompartmentalized structures with boundaries that delineate the interior from exterior environments and is an inspiration for complex functional materials. Here, we have demonstrated an approach to the design and construction of a nested protein cage system that can additionally incorporate the packing of other functional macromolecules and exhibit some of the features of a minimal synthetic cell-like material. We have demonstrated a strategy for controlled co-packaging of subcompartments, ferritin (Fn) cages, together with active cellobiose-hydrolyzing β-glycosidase enzyme macromolecules, CelB, inside the sequestered volume of the bacteriophage P22 capsid. Using controlled in vitro assembly, we were able to modulate the stoichiometry of Fn cages and CelB encapsulated inside the P22 to control the degree of compartmentalization. The co-encapsulated enzyme CelB showed catalytic activity even when packaged at high total macromolecular concentrations comparable to an intracellular environment. Here, this approach could be used as a model to create synthetic protein-based protocells that can confine smaller functionalized proto-organelles and additional macromolecules to support a range of biochemical reactions.

59 BASIC BIOLOGICAL SCIENCES↗

Using deep machine learning to interpret proton radiography data from a pulsed power experiment

Deep machine learning is used to analyze a proton radiograph from a tin pulsed power experiment and determine density values for each pixel in the image. Two promising convolutional neural network architectures that have proven to be effective for image analysis in other applications are applied to analyze a proton radiograph and find density values. The process of creating a suitable training dataset is described, involving the Lagrangian hydrodynamic model used for simulations of the experiment, the proton radiography forward model to make synthetic images for training, and the manner in which data augmentation is used to expand the resulting image dataset. It is shown that machine learning not only produces a reasonable density field but is also able to predict features in the density field that are suggested by the proton radiograph but not captured by simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗

Simulated JWST Data Sets for Multispectral and Hyperspectral Image Fusion

The James Webb Space Telescope (JWST) will provide multispectral and hyperspectral infrared images of a large number of astrophysical scenes. Multispectral images will have the highest angular resolution, while hyperspectral images (e.g., with integral field unit spectrometers) will provide the best spectral resolution. This paper aims at providing a comprehensive framework to generate an astrophysical scene and to simulate realistic hyperspectral and multispectral data acquired by two JWST instruments, namely, NIRCam Imager and NIRSpec IFU. We want to show that this simulation framework can be resorted to assess the benefits of fusing these images to recover an image of high spatial and spectral resolutions. To do so, we make a synthetic scene associated with a canonical infrared source, the Orion Bar. We develop forward models including corresponding noises for the two JWST instruments based on their physical features. JWST observations are then simulated by applying the forward models to the aforementioned synthetic scene. We test a dedicated fusion algorithm we developed on these simulated observations. We show that the fusion process reconstructs the high spatio-spectral resolution scene with a good accuracy on most areas, and we identify some limitations of the method to be tackled in future works. The synthetic scene and observations presented in the paper can be used, for instance, to evaluate instrument models, pipelines, or more sophisticated algorithms dedicated to JWST data analysis. Besides, fusion methods such as the one presented in this paper are shown to be promising tools to fully exploit the unprecedented capabilities of the JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

DisasterAWARE – A Global Alerting Platform for Flood Events

The rising number of flooding events combined with increased urbanization is contributing to significant economic losses due to damages to structures and infrastructures. From a risk reduction and resilience perspective, it is not only essential to forecast flood risk and potential impacts, but also to disseminate the information to stakeholders on the ground for rapid implementation of mitigation and response measures. This paper provides (i) an introduction to DisasterAWARE®, a global alerting system, that is used to disseminate flood risk information to stakeholders across the globe, and (ii) a discussion of the models implemented using earth observation data (Synthetic Aperture Radar and optical imagery) for near real-time assessment of flood severity and potential flood impacts to infrastructures. While the models are still in their nascent stage, a case study implementation of the models for the 2020 flooding event in Africa is presented to showcase the model integration with DisasterAWARE®.

54 ENVIRONMENTAL SCIENCES↗

Multilabel proportion prediction and out-of-distribution detection on gamma spectra of short-lived fission products

In the machine learning problem of multilabel classification, the objective is to determine for each test instance which classes the instance belongs to. In this work, we consider an extension of multilabel classification, called multilabel proportion prediction, in the context of radioisotope identification (RIID) using gamma spectra data. We aim to not only predict radioisotope proportions, but also identify out-of-distribution (OOD) spectra. We achieve this goal by viewing gamma spectra as discrete probability distributions, and based on this perspective, we develop a custom semi-supervised loss function that combines a traditional supervised loss with an unsupervised reconstruction error function. Our approach was motivated by its application to the analysis of short-lived fission products from spent nuclear fuel. In particular, we demonstrate that a neural network model trained with our loss function can successfully predict the relative proportions of 37 radioisotopes simultaneously. The model trained with synthetic data was then applied to measurements taken by Pacific Northwest National Laboratory (PNNL) to conduct analysis typically done by subject-matter experts. Here, we also extend our approach to successfully identify when measurements are OOD, and thus should not be trusted, whether due to the presence of a novel source or novel proportions.

Anomaly detection↗

Unsupervised domain adaptation for radioisotope identification in gamma spectroscopy

Training machine learning models for radioisotope identification using gamma spectroscopy remains an elusive challenge for many practical applications, largely stemming from the difficulty of acquiring and labeling large, diverse experimental datasets. Simulations can mitigate this challenge, but the accuracy of models trained on simulated data can deteriorate substantially when deployed to an out-of-distribution operational environment. In this study, we demonstrate that unsupervised domain adaptation (UDA) can improve the ability of a model trained on synthetic data to generalize to a new testing domain, provided unlabeled data from the target domain are available. Conventional supervised techniques are unable to utilize this data because the absence of isotope labels precludes defining a supervised classification loss. Instead, we first pretrain a spectral classifier using labeled synthetic data and subsequently leverage unlabeled target data to align the learned feature representations between the source and target domains. We compare a range of different UDA techniques, finding that minimizing the maximum mean discrepancy (MMD) between source and target feature vectors yields the most consistent improvement to testing scores. For instance, using a custom transformer-based neural network, we achieved a testing accuracy of $0.904 \pm 0.022$ on an experimental LaBr test set after performing unsupervised feature alignment via MMD minimization, compared to $0.754 \pm 0.014$ before alignment. Overall, our results highlight the potential of using UDA to adapt a radioisotope classifier trained on synthetic data for real-world deployment.

Lalor, Peter W.↗

Enhancing power grid resilience to winter storms via generator winterization with equity considerations

Here we develop two-stage stochastic programming models for generator winterization that enhance power grid resilience while incorporating social equity. The first stage in our models captures the investment decisions for generator winterization, and the second stage captures the operation of a degraded power grid, with the objective of minimizing load shed and social inequity. To incorporate equity into our models, we propose a concept called adverse effect probability that captures the disproportionate effects of power outages on communities with varying vulnerability levels. Grid operations are modeled using DC power flow, and equity is captured through mean or maximum adverse effects experienced by communities. We apply our models to a synthetic Texas power grid, using winter storm scenarios created from the generator outage data from the 2021 Texas winter storm. Our extensive numerical experiments show that more equitable outcomes, in the sense of reducing adverse effects experienced by vulnerable communities during power outages, are achievable with no impact on total load shed through investing in winterization of generators in different locations and capacities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

The Imprint of Southern Ocean Storms on Modeled Surface Chlorophyll, Their Drivers and Satellite Biases

Southern Ocean (SO) phytoplankton chlorophyll is highly variable on sub-seasonal time scales. Although the SO is the windiest ocean basin globally, it is not conclusively understood how storms impact SO phytoplankton dynamics. Much of our existing knowledge stems from satellites, but biases due to data gaps from cloud cover and low solar angles remain unquantified. Here, we use ocean–sea-ice simulations with the Community Earth System Model to quantify the climatological 1997–2018 imprint of storms on chlorophyll and phytoplankton dynamics in the ice-free SO. Additionally, by comparing the full-field model output to synthetic satellite observations, we quantify sampling biases in satellite-derived estimates. We find that both the sign and the magnitude of the average surface chlorophyll imprint vary substantially across storms but last for at least 4 days after the storm passing. Based on our analysis, more than one third of the storms explain the majority of local non-seasonal chlorophyll variability, but satellite-derived storm imprints are often too large in magnitude. On the day of the storm passing, changes in vertical mixing predominantly cause surface chlorophyll anomalies, and reduced light availability due to enhanced cloud cover outweighs the enhanced nutrient availability due to entrainment. Interestingly, storms imprint differently on total net primary production than on surface chlorophyll, demonstrating the difficulty to derive carbon-cycle impacts from a surface-chlorophyll assessment. With SO future storm activity projected to increase, complementing satellite observations with other observing technologies, for example, profiling floats, is necessary to better constrain how storms impact biological carbon cycling in the SO.

Nissen, Cara [University of Colorado, Boulder, CO ↗

Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials

There is a need to characterize complex materials and their dynamics under reaction conditions to accelerate materials design. Adsorbate vibrational excitations are selective to adsorbate/surface interactions and infrared (IR) spectra associated with activating adsorbate vibrational modes are accurate, capture details of most modes, and can be obtained operando. Current interpretation depends on heuristic peak assignments for simple spectra, precluding the possibility of obtaining detailed structural information. Here, we combine data-based approaches with chemistry-dependent problem formulation to develop physics-driven surrogate models that generate synthetic IR spectra from first-principles calculations. Using synthetic IR spectra of carbon monoxide on platinum, we implement multinomial regression via neural network ensembles to learn probability distributions functions (pdfs) that describe adsorption sites and quantify uncertainty. We use these pdfs to infer detailed surface microstructure from experimental spectra and extend this methodology to other systems as a first step towards characterizing complex interfaces and closing the materials gap.

Science & Technology - Other Topics↗

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-dimensional talc as a natural abundant ultra-broadband hyperbolic material

Here, we demonstrate that two-dimensional talc, a naturally abundant phyllosilicate mineral, supports hyperbolic phonon-polaritons (HPhPs) across the mid- and far-infrared wavelengths. Using scattering scanning near-field optical microscopy (s-SNOM) and synchrotron infrared nano-spectroscopy (SINS), we reveal tunable HPhP modes in talc flakes with long lifetimes, high confinement, and quality factors of up to 5. We further observe Fabry–Pérot cavity modes in tapered flakes, confirmed by simulations and analytical modeling. Compared to synthetic crystals, talc offers an ultra-broadband, low-cost, and sustainable platform for infrared nanophotonics and optoelectronics.

36 MATERIALS SCIENCE↗