Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse analysis”

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 181 records · Page 10

Why matter effects matter for JUNO

In this paper we focus on the Earth matter effects for the solar parameter determination by a medium baseline reactor experiment such as JUNO. We derive perturbative expansions for the mixing angles θ 12 and θ 13 as well as the Δm$^{2}_{21}$ and Δm$^{2}_{31}$ in terms of the matter potential relevant for JUNO. These expansions, up to second order in the matter potential, while simple, allow one to calculate the electron antineutrino survival probability to a precision much better than needed for the JUNO experiment. We use these perturbative expansions to semi-analytically explain and confirm the shift caused by the matter effects on the solar neutrino mixing parameters θ 12 and Δm$^{2}_{21}$ which were previously obtained by a purely numerical χ 2 analysis. Since these shifts do not satisfy the naive expectations and are significant given the precision that can be achieved by the JUNO experiment, a totally independent cross check using a completely different method is of particular importance. We find that these matter effect shifts do not depend on any of the details of the detector characteristics apart from the baseline and earth mass density between reactor(s) and detector, but do depend on the normalized product of reactor neutrino spectrum times the inverse-beta decay cross-section. The results of this manuscript suggests an alternative analysis method for measuring sin 2 θ 12 and Δm$^{2}_{21}$ in JUNO which would be a useful cross check of the standard analysis and for the understanding of the Wolfenstein matter effect. The explanation of these shifts together with a quantitative understanding, using a semi-analytical method, is the principal purpose of this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗

Search for a Sub-eV Sterile Neutrino Using Daya Bay’s Full Dataset

This Letter presents results of a search for the mixing of a sub-eV sterile neutrino with three active neutrinos based on the full data sample of the Daya Bay Reactor Neutrino Experiment, collected during 3158 days of detector operation, which contains 5.55 × 106 reactor $\overline{v}$ e candidates identified as inverse beta-decay interactions followed by neutron capture on gadolinium. The analysis benefits from a doubling of the statistics of our previous result and from improvements of several important systematic uncertainties. No significant oscillation due to mixing of a sub-eV sterile neutrino with active neutrinos was found. Exclusion limits are set by both Feldman-Cousins and CLs methods. Light sterile neutrino mixing with sin 2⁡ 2⁢θ 14 ≳ 0.01 can be excluded at 95% confidence level in the region of 0.01 eV 2 ≲ |Δ⁢$m$$^{2}_{41}$| ≲ 0.1 eV 2 . This result represents the world-leading constraints in the region of 2 × 10 –4 eV 2 ≲ |Δ⁢$m$$^{2}_{41}$| ≲ 0.2 eV 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uranium Mass Attribute System (UMAS) Study: Evaluation of Options against the Arms Control Evaluation Criteria

In this FY19 study, a range of possible measurement systems for confirming a U-235 mass attribute were identified. The systems were reviewed for their applicability in a potential future arms control agreement, and a subset of the eight most plausible systems were selected for further evaluation against the Arms Control Evaluation Criteria (ACEC). The evaluation results were then used to identify a primary and secondary recommended measurement system for confirming a U-235 mass attribute on warheads, nuclear components of warheads, and/or materials as part of a hypothetical arms control agreement.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Dissolution of Mg-enriched defects in implanted GaN and increased p-type dopant activation

Annealing Mg-implanted homoepitaxial GaN at temperatures above 1400 °C eliminates the formation of inversion domains and leads to improved dopant activation efficiency. Extended defects, in the form of inversion domains, contain electrically inactive Mg after post-implantation annealing at temperatures as high as 1300 °C (one GPa N 2 overpressure), which results in a low dopant activation efficiency. Triple-axis x-ray data reveal that implant-induced strain is fully relieved after annealing at 1300 °C for 10 min, indicating that strain-inducing point defects formed during implantation have reconfigured and inversion domains are formed. However, annealing at temperatures of 1400–1500 °C (one GPa N 2 overpressure) eliminates the presence of the inversion domains. While residual defects, such as dislocation loops, still exist after annealing at and above 1400 °C, chemical analysis at multiple dislocation loops shows no sign of Mg segregation. Meanwhile, an overall decreasing trend in the dislocation loop density is observed after annealing at the higher temperatures and longer times. Additionally, once inversion domains are formed and the samples are cooled to room temperature, they are shown to dissolve with subsequent annealing above 1400 °C. While such defects have been observed before, the important finding that such defects can be dissolved with a short, higher temperature step is key. Earlier work [Breckenridge et al., J. Appl. Phys. Lett. 118, 022101 (2021)] addressing electrical measurements of these types of samples showed that annealing at 1400 °C leads to a dopant activation efficiency that is an order of magnitude higher than that observed at 1300 °C. In conclusion, this work complements earlier work by identifying the inversion domains, which incorporate Mg, and points to the benefits, in terms of defect density and p-type dopant activation, of using higher temperature (>1400 °C) annealing cycles to activate Mg in GaN, even if the Mg-containing inversion domains had been formed during lower temperature annealing.

36 MATERIALS SCIENCE↗

Towards Interactive, Reproducible Analytics at Scale on HPC Systems

The growth in scientific data volumes has resulted in a need to scale up processing and analysis pipelines using High Performance Computing (HPC) systems. These workflows need interactive, reproducible analytics at scale. The Jupyter platform provides core capabilities for interactivity but was not designed for HPC systems. In this paper, we outline our efforts that bring together core technologies based on the Jupyter Platform to create interactive, reproducible analytics at scale on HPC systems. Our work is grounded in a real world science use case-applying geophysical simulations and inversions for imaging the subsurface. Our core platform addresses three key areas of the scientific analysis workflow-reproducibility, scalability, and interactivity. We describe our implemention of a system, using Binder, Science Capsule, and Dask software. We demonstrate the use of this software to run our use case and interactively visualize real-Time streams of HDF5 data.

containers↗

Automated CT registration, segmentation, and quantification (AutoCT) v1.0

Processing and analyzing brain imaging is crucial in both scientific development and clinical field. In this software package, we build a pipeline that integrates automatic registration, segmentation, and quantitative analysis for subjects' CT scans. Leveraging diffeomorphic transformations, we enable optimized forward and inverse mappings between an image and the reference. Furthermore, we extract localized features from deformation field based on an online template process, which advances statistical learning downstream. The created templates, atlas as well as our methods provide the brain imaging community tools for AI implementations.

Essiari, Abdelilah↗

Automated CT registration, segmentation, and quantification (AutoCT) v1.1

Processing and analyzing brain imaging is crucial in both scientific development and clinical field. In this software package, we build a pipeline that integrates automatic registration, segmentation, and quantitative analysis for subjects' CT scans. Leveraging diffeomorphic transofrmations, we enable optimized forward and inverse mappings between an image and the reference. Furthermore, we extract localized features from deformation field based on an online template process, which advances statistical learning downstream. The created templates, atlas as well as our methods provide the brain imaging community tools for AI implementations

Bai, Zhe↗

pnnl/MCRASTA

McRasta (Markov Chain Rate and State Analysis) was developed to estimate parameter uncertainty in constitutive friction models via Bayesian inverse and Markov Chain Monte Carlo (MCMC) methods.

Fichera, Marissa [Pacific Northwest National Labor↗

Tethered balloon measurements reveal enhanced aerosol occurrence aloft interacting with Arctic low-level clouds

Low-level clouds in the Arctic affect the surface energy budget and vertical transport of heat and moisture. The limited availability of cloud-droplet-forming aerosol particles strongly impacts cloud properties and lifetime. Vertical particle distributions are required to study aerosol–cloud interaction over sea ice comprehensively. This article presents vertically resolved measurements of aerosol particle number concentrations and sizes using tethered balloons. The data were collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition in the summer of 2020. Thirty-four profiles of aerosol particle number concentration were observed in 2 particle size ranges: 12–150 nm (N 12−150 ) and above 150 nm (N >150 ). Concurrent balloon-borne meteorological measurements provided context for the continuous profiles through the cloudy atmospheric boundary layer. Radiosoundings, cloud remote sensing data, and 5-day back trajectories supplemented the analysis. The majority of aerosol profiles showed more particles above the lowest temperature inversion, on average, double the number concentration compared to below. Increased N 12−150 up to 3,000 cm −3 were observed in the free troposphere above low-level clouds related to secondary particle formation. Long-range transport of pollution increased N >150 to 310 cm −3 in a warm, moist air mass. Droplet activation inside clouds caused reductions of N >150 by up to 100%, while the decrease in N 12−150 was less than 50%. When low-level clouds were thermodynamically coupled with the surface, profiles showed 5 times higher values of N 12−150 in the free troposphere than below the cloud-capping temperature inversion. Enhanced N 12−150 and N >150 interacting with clouds were advected above the lowest inversion from beyond the sea ice edge when clouds were decoupled from the surface. Vertically discontinuous aerosol profiles below decoupled clouds suggest that particles emitted at the surface are not transported to clouds in these conditions. It is concluded that the cloud-surface coupling state and free tropospheric particle abundance are crucial when assessing the aerosol budget for Arctic low-level clouds over sea ice.

54 ENVIRONMENTAL SCIENCES↗

Single-shot Quantum Signal Processing Interferometry

Quantum systems of infinite dimension, such as bosonic oscillators, provide vast resources for quantum sensing. Yet, a general theory on how to manipulate such bosonic modes for sensing beyond parameter estimation is unknown. We present a general algorithmic framework, quantum signal processing interferometry (QSPI), for quantum sensing at the fundamental limits of quantum mechanics by generalizing Ramsey-type interferometry. Our QSPI sensing protocol relies on performing nonlinear polynomial transformations on the oscillator's quadrature operators by generalizing quantum signal processing (QSP) from qubits to hybrid qubit-oscillator systems. We use our QSPI sensing framework to make efficient binary decisions on a displacement channel in the single-shot limit. Theoretical analysis suggests the sensing accuracy, given a single-shot qubit measurement, scales inversely with the sensing time or circuit depth of the algorithm. We further concatenate a series of such binary decisions to perform parameter estimation in a bit-by-bit fashion. Numerical simulations are performed to support these statements. Our QSPI protocol offers a unified framework for quantum sensing using continuous-variable bosonic systems beyond parameter estimation and establishes a promising avenue toward efficient and scalable quantum control and quantum sensing schemes beyond the NISQ era.

Physics↗

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗

Multiscale Data-Driven Seismic Full-Waveform Inversion With Field Data Study

Seismic full-waveform inversion (FWI), which uses iterative methods to estimate high-resolution subsurface models from seismograms, is a powerful imaging technique in exploration geophysics. In recent years, the computational cost of FWI has grown exponentially due to the increasing size and resolution of seismic data. Moreover, it is a nonconvex problem and can encounter local minima due to the limited accuracy of the initial velocity models or the absence of low frequencies in the measurements. To overcome these computational issues, we develop a multiscale data-driven FWI method based on fully convolutional networks (FCNs). In preparing the training data, we first develop a real-time style transform method to create a large set of synthetic subsurface velocity models from natural images. We then develop two convolutional neural networks with encoder-decoder structures to reconstruct the low- and high-frequency components of the subsurface velocity models, separately. To validate the performance of our data-driven inversion method and the effectiveness of the synthesized training set, we compare it with conventional physics-based waveform inversion approaches using both synthetic and field data. Finally, these numerical results demonstrate that, once our model is fully trained, it can significantly reduce the computation time and yield more accurate subsurface velocity models in comparison with conventional FWI.

58 GEOSCIENCES↗

DEEPEN: Newberry Volcano MT and Gravity Data - 2022 and 2023 Acquisition and Processing

As part of DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments), a 3D play fairway analysis (PFA) was conducted at Newberry Volcano in Central Oregon for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). For use in this PFA, combined full tensor broadband magnetotelluric (MT) and gravity data were acquired, processed, and inverted by Enthalpion Energy LLC (Enthalpion) with support from NREL staff. The data collection efforts took place from June 19th to July 24. Data were collected with the goal of gaining an improved understanding of the South Flank and the extent of the magma chamber. This GDR submission includes the raw data, single inversions, joint inversions, resolution matrices, and a report on these processes. More detailed information about the folder structure of the datasets included in this submission may be found on pages 94-97 of Magnetotelluric and Gravity Survey Report.pdf below.

15 GEOTHERMAL ENERGY↗

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part II

X-ray Bragg coherent diffraction imaging (BCDI) has been demonstrated as a powerful 3D microscopy approach for the investigation of sub-micrometre-scale crystalline particles. The approach is based on the measurement of a series of coherent Bragg diffraction intensity patterns that are numerically inverted to retrieve an image of the spatial distribution of the relative phase and amplitude of the Bragg structure factor of the diffracting sample. This 3D information, which is collected through an angular rotation of the sample, is necessarily obtained in a non-orthogonal frame in Fourier space that must be eventually reconciled. To deal with this, the approach currently favored by practitioners (detailed in Part I) is to perform the entire inversion in conjugate non-orthogonal real- and Fourier-space frames, and to transform the 3D sample image into an orthogonal frame as a post-processing step for result analysis. In this article, which is a direct follow-up of Part I, two different transformation strategies are demonstrated, which enable the entire inversion procedure of the measured data set to be performed in an orthogonal frame. The new approaches described here build mathematical and numerical frameworks that apply to the cases of evenly and non-evenly sampled data along the direction of sample rotation ( i.e. the rocking curve). The value of these methods is that they rely on the experimental geometry, and they incorporate significantly more information about that geometry into the design of the phase-retrieval Fourier transformation than the strategy presented in Part I. Finally, two important outcomes are (1) that the resulting sample image is correctly interpreted in a shear-free frame and (2) physically realistic constraints of BCDI phase retrieval that are difficult to implement with current methods are easily incorporated. Computing scripts are also given to aid readers in the implementation of the proposed formalisms.

36 MATERIALS SCIENCE↗

Search for Non-resonant Higgs Boson Pair Production in the Four Bottom Quark Decay Channel With the CMS Experiment

This dissertation presents a search for non-resonant Higgs boson pair production, focusing on the four bottom quark decay channel. It explores the gluon fusion and vector boson fusion production mechanisms. The analysis is performed with a dataset of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS detector at the LHC, corresponding to an integrated luminosity of around 138 inverse femtobarns. Innovative techniques in the areas of object identification, event categorization, signal identification, and background modeling are used to maximize the analysis sensitivity. No excess of signal events are observed relative to the background-only expectation, and 95% CL upper limits on the production cross section are set. The observed upper limit on the standard model production cross section is set at 3.6 times the theoretical expectation. At the time of writing this dissertation, it is the most stringent constraint at the LHC from an individual channel. Furthermore, the observed constraint on the coupling modifier of the Higgs boson self-interaction is set between -2.3 and 9.4. The observed constraint on the coupling modifier of the di-vector-boson-di-Higgs-boson interaction is set between -0.1 and 2.2.

Guerrero Ibarra, Guerrero Fernando↗

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗