Search NASA⌕ Search

SEARCH · Search NASA

Results for “large-scale inverse problems”

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.

64 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Isospin symmetry breaking in the mirror pair 73 Sr- 73 Br

The recent experimental observation of isospin symmetry breaking (ISB) in the ground states of the T=3/2 mirror pair 73 Sr- 73 Br is theoretically studied using large-scale shell-model calculations. The large valence space and the successful PFSDG-U effective interaction used for the nuclear part of the problem capture possible structural changes and provide a robust basis to treat the ISB effects of both electromagnetic and nonelectromagnetic origin. The calculated shifts and mirror-energy differences are consistent with the inversion of the I π =1/2 - ,5/2 - states between 73 Sr- 73 Br and suggest that the role played by the Coulomb interaction is dominant. Finally, an isospin breaking contribution of nuclear origin is estimated to be ≈25keV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The Artificial Scientist: in-Transit Machine Learning of Plasma Simulations

Large-scale simulations or scientific experiments produce petabytes of data per run. This poses massive challenges for I/O and storage when scientific analysis workflows are run manually offline. Unsupervised deep learning-based techniques to extract patterns and non-linear relations from these large amounts of data provide a way to build scientific understanding from raw data, reducing the need for manual pre-selection of analysis steps, but require exascale compute and memory to process the full dataset available. In this paper, we demonstrate a heterogeneous streaming workflow in which plasma simulation data is streamed directly to a Machine Learning (ML) application training a model on the simulation data in-transit, completely circumventing the capacity-constrained filesystem bottleneck. This workflow employs openPMD to provide a high level interface to describe scientific data and also uses ADIOS2, to transfer volumes of data that exceed the capabilities of the filesystem. We employ experience replay to avoid catastrophic forgetting in learning from this non-steady state process in a continual manner and adapt it to improve model convergence while learning in-transit. As a proof-of-concept, we approach the ill-posed inverse problem of predicting particle dynamics from radiation in a particle-incell (PIConGPU) simulation of the Kelvin-Helmholtz instability (KHI). We detail hardware-software co-design challenges as we scale PIConGPU to full Frontier, the Top-1 system as of June 2024 Top500 list.

Kelling, Jeffrey [Helmholtz-Zentrum Dresden Rossen↗
DOE OSTI2022

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗
DOE OSTI2023

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗
DOE OSTI2020

Turbulent Model of Crab Nebula Radiation

We construct a turbulent model of the Crab Nebula's nonthermal emission. The present model resolves a number of long-standing problems of the Kennel–Coroniti model: (i) the sigma problem, (ii) the hard spectrum of radio electrons, (iii) the high peak energy of gamma-ray flares, (iv) and the spatial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I, accelerated at the wind termination shock via the Fermi-I mechanism; and Component-II, accelerated in reconnecting turbulence in highly magnetized (σ >> 1) plasma in the central part of the Crab Nebula. The reconnecting turbulence in Component-II extends from radio to gamma-rays: it accelerates radio electrons with a hard spectrum, destroys the large-scale magnetic flux (and thus resolves the sigma problem), and occasionally produces gamma-ray flares (from the largest-scale reconnection events). The model reproduces the broadband spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as the spatially resolved evolution of the spectral indices in the IR and optical bands.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Large-scale kinetic energy spectra from Eulerian analysis of EOLE wind data

A data set of 56,000 winds determined from the horizontal displacements of EOLE balloons at the 200 mb level in the Southern Hemisphere during the period October 1971-February 1972 is utilized for the computation of planetary- and synoptic-scale kinetic energy space spectra. However, the random distribution of measurements in space and time presents some problems for the spectral analysis. Two different approaches are used, i.e., a harmonic analysis of daily wind values at equi-distant points obtained by space-time interpolation of the data, and a correlation method using the direct measurements. Both methods give similar results for small wavenumbers, but the second is more accurate for higher wavenumbers (k above or equal to 10). The spectra show a maximum at wavenumbers 5 and 6 due to baroclinic instability and then decrease for high wavenumbers up to wavenumber 35 (which is the limit of the analysis), according to the inverse power law k to the negative p, with p close to 3.

Desbois, M.↗
NASA NTRS1977

Mass drivers. 2: Structural dynamics

Various structural and dynamical problems related to both small-scale forces between the drive coils and within the bucket structure as well as the overall combined large-scale dynamical interaction of the bucket stream and MDRE (Mass Drive Reaction Engine) structure are examined. The large-scale dynamics appear weakly stable. Finally, MDRE operation in an inverse-square-law gravitational field is discussed and the required curved shape of the guideway is computed.

William H. Arnold↗
NASA NTRS1992

Large-scale sparse singular value computations

Four numerical methods for computing the singular value decomposition (SVD) of large sparse matrices on a multiprocessor architecture are presented. Lanczos and subspace iteration-based methods for determining several of the largest singular triplets (singular values and corresponding left and right-singular vectors) for sparse matrices arising from two practical applications: information retrieval and seismic reflection tomography are emphasized. The target architectures for implementations are the CRAY-2S/4-128 and Alliant FX/80. The sparse SVD problem is well motivated by recent information-retrieval techniques in which dominant singular values and their corresponding singular vectors of large sparse term-document matrices are desired, and by nonlinear inverse problems from seismic tomography applications which require approximate pseudo-inverses of large sparse Jacobian matrices.

Berry, Michael W.↗
NASA NTRS1996

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗
NASA NTRS2021

Modeling a Coronal Mass Ejection from an Extended Filament Channel. I. Eruption and Early Evolution

We present observations and modeling of the magnetic field configuration, morphology, and dynamics of a large-scale, high-latitude filament eruption observed by the Solar Dynamics Observatory. We analyze the 2015 July 9–10 filament eruption and the evolution of the resulting coronal mass ejection (CME) through the solar corona. The slow streamer-blowout CME leaves behind an elongated post-eruption arcade above the extended polarity inversion line that is only poorly visible in extreme ultraviolet (EUV) disk observations and does not resemble a typical bright flare-loop system. Magnetohydrodynamic (MHD) simulation results from our data-inspired modeling of this eruption compare favorably with the EUV and white-light coronagraph observations. We estimate the reconnection flux from the simulation’s flare-arcade growth and examine the magnetic-field orientation and evolution of the erupting prominence, highlighting the transition from an erupting sheared-arcade filament channel into a streamer-blowout flux-rope CME. Our results represent the first numerical modeling of a global-scale filament eruption where multiple ambiguous and complex observational signatures in EUV and white light can be fully understood and explained with the MHD simulation. In this context, our findings also suggest that the so-called “stealth CME” classification, as a driver of unexpected or “problem” geomagnetic storms, belongs more to a continuum of observable/non-observable signatures than to separate or distinct eruption processes.

Benjamin J Lynch↗