Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The Radiative Effects of Martian Water Ice Clouds on the Local Atmospheric Temperature Profile

Mars Pathfinder made numerous discoveries, one of which was a deep temperature inversion that extended from about 15 km down to 8 km above the surface. It has been suggested by Haberle et al. (1999. J. Geophys. Res. 104, 8957-8974.) that radiative cooling by a water ice cloud may generate such an inversion. Clouds can strongly affect the local air temperature due to their ability to radiate efficiently in the infrared and due to the low air mass of the martian atmosphere, which allows the temperature to change during the relatively short lifetime of a cloud. We utilize a time-dependent microphysical aerosol model coupled to a radiative--convective model to explore the effects water ice clouds have on the local martian temperature profile. We constrain the dust and water vapor abundance using data from the Viking Missions and Mars Pathfinder. Water t ice clouds with visible optical depths of r > 0.1 form readily in these simulations. These clouds alter the local air temperature directly, through infrared cooling, and indirectly, by redistributing atmospheric dust. With this model we are able to reproduce the temperature inversions observed by Mars Pathfinder and Mars Global t Surveyor 2000 Academic Press

Colaprete, Anthony↗
NASA NTRS2014

Polarimetric SAR Interferometry Evaluation in Mangroves

TanDEM-X (TDX) enables to generate an interferometric coherence without temporal decorrelation effect that is the most critical factor for a successful Pol-InSAR inversion, as have recently been used for forest parameter retrieval. This paper presents mangrove forest height estimation only using single-pass/single-baseline/dual-polarization TDX data by means of new dual-Pol-InSAR inversion technique. To overcome a lack of one polarization in a conventional Pol- InSAR inversion (i.e. an underdetermined problem), the ground phase in the Pol-InSAR model is directly estimated from TDX interferograms assuming flat underlying topography in mangrove forest. The inversion result is validated against lidar measurement data (NASA's G-LiHT data).

Mangrove↗
NASA NTRS2015

An Inversion Analysis of Recent Variability in Natural CO2 Fluxes Using GOSAT and In Situ Observations

About one-half of the global CO2 emissions from fossil fuel combustion and deforestation accumulates in the atmosphere, where it contributes to global warming. The rest is taken up by vegetation and the ocean. The precise contribution of the two sinks, and their location and year-to-year variability are, however, not well understood. We use two different approaches, batch Bayesian synthesis inversion and variational data assimilation, to deduce the global spatiotemporal distributions of CO2 fluxes during 2009-2010. One of our objectives is to assess different sources of uncertainties in inferred fluxes, including uncertainties in prior flux estimates and observations, and differences in inversion techniques. For prior constraints, we utilize fluxes and uncertainties from the CASA-GFED model of the terrestrial biosphere and biomass burning driven by satellite observations and interannually varying meteorology. We also use measurement-based ocean flux estimates and two sets of fixed fossil CO2 emissions. Here, our inversions incorporate column CO2 measurements from the GOSAT satellite (ACOS retrieval, filtered and bias-corrected) and in situ observations (individual flask and afternoon-average continuous observations) to estimate fluxes in 108 regions over 8-day intervals for the batch inversion and at 3 x 3.75 weekly for the variational system. Relationships between fluxes and atmospheric concentrations are derived consistently for the two inversion systems using the PCTM atmospheric transport model driven by meteorology from the MERRA reanalysis. We compare the posterior fluxes and uncertainties derived using different data sets and the two inversion approaches, and evaluate the posterior atmospheric concentrations against independent data including aircraft measurements. The optimized fluxes generally resemble those from other studies. For example, the results indicate that the terrestrial biosphere is a net CO2 sink, and a GOSAT-only inversion suggests a shift in the global sink from the tropics south to the north relative to the prior and to an in-situ-only inversion. We also find a smaller terrestrial sink in higher-latitude northern regions in boreal summer of 2010 relative to 2009.

Carbon dioxide↗
DOE OSTI2022

Quantum error mitigation by hidden inverses protocol in superconducting quantum devices

We present a method to improve the convergence of variational algorithms based on hidden inverses (HIs) to mitigate coherent errors. In the context of error mitigation, this means replacing the hardware implementation of certain Hermitian gates with their inverses. Doing so results in noise cancellation and a more resilient quantum circuit. This approach improves performance in a variety of two-qubit error models where the noise operator also inverts with the gate inversion. We apply the mitigation scheme on superconducting quantum processors running the variational quantum eigensolver (VQE) algorithm to find the H 2 ground-state energy. When implemented on superconducting hardware we find that the mitigation scheme effectively reduces the energy fluctuations in the parameter learning path in VQE, reducing the number of iterations for a converged value. We also provide a detailed numerical simulation of VQE performance under different noise models and explore how HIs and randomized compiling affect the underlying loss landscape of the learning problem. These simulations help explain our experimental hardware outcomes, helping to connect lower-level gate performance to application-specific behavior in contrast to metrics like fidelity which often do not provide an intuitive insight into observed high level performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

The relationship between wind vector and normalized radar cross section used to derive Seasat-A Satellite Scatterometer winds

The Seasat-A Satellite Scatterometer (SASS) ocean normalized radar cross section (NRCS) dependence on the 19.5-m neutral stability wind vector may be specified as a function of radar incidence angle, the angle between wind direction and radar azimuth, and the neutral stability wind speed expressed in m/sec at a height of 19.5 m. An account is given of the development of models both expressing this relationship and providing the basis of inversion of NRCS to SASS winds, from initially aircraft scatterometer measurement-based forms to three Seasat field-validation experiments which furnish model NRCS versus surface windspeed data for comparison with SASS data.

Schroeder, L. C.↗
NASA NTRS1989

The evolution and discharge of electric fields within a thunderstorm

An analysis of the present three-dimensional thunderstorm electrical model and its finite-difference approximations indicates unconditional stability for the discretization that results from the approximation of the spatial derivatives by a box-schemelike method and of the temporal derivative by either a backward-difference or Crank-Nicholson scheme. Lightning propagation is treated through numerical techniques based on the inverse-matrix modification formula and Cholesky updates. The model is applied to a storm observed at the Kennedy Space Center in 1978, and numerical comparisons are conducted between the model and the theoretical results obtained by Wilson (1920) and Holzer and Saxon (1952).

Hager, William W.↗
NASA NTRS1991

Combined structures-controls optimization of lattice trusses

The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.

Balakrishnan, A. V.↗
DOE OSTI2026

Landau–Zener Transition Enhanced Quantum Sensing in Spin Defects of Hexagonal Boron Nitride

Negatively charged boron vacancies (V $^{–}_{B}$ ) in hexagonal boron nitride (hBN) comprise a promising quantum sensing platform, optically addressable at room temperature and transferable onto samples. However, broad hyperfine-split spin transitions of the ensemble pose challenges for quantum sensing with conventional resonant excitation due to limited spectral coverage. While V $^{–}_{B}$ in isotopically enriched hBN using 10 B and 15 N isotopes (h 10 B 15 N) exhibits sharper spectral features, significant inhomogeneous broadening persists. We show that, implemented via frequency modulation on an FPGA, a frequency-ramped microwave pulse achieves around 4-fold greater |0⟩→|−1⟩ spin-state population transfer and thus contrast than resonant microwave excitation and thus 16-fold shorter measurement time for spin relaxation-based quantum sensing. Quantum dynamics simulations reveal that an effective two-state Landau–Zener model captures the complex relationship between population inversion and pulse length with relaxations incorporated. Our approach is robust and valuable for quantum relaxometry with spin defects in hBN, especially in noisy environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗
DOE OSTI2024

Asymmetric reheating via inverse symmetry breaking

Asymmetric reheating is a generic requirement for models of dark sectors with light species, but its implementation is usually in tension with unique phenomenologies otherwise possible in compelling theories containing dark copies of the Standard Model. We present a simple module to implement asymmetric reheating during a Z 2 -breaking phase above some critical temperature. This reinvigorates the possibility of an exactly degenerate mirror sector and the striking phenomenology of composite particles oscillating into their mirror counterparts. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

Bicircular Light Floquet Engineering of Magnetic Symmetry and Topology and Its Application to the Dirac Semimetal Cd 3 As 2

Here, we show that bicircular light (BCL) is a versatile way to control magnetic symmetries and topology in materials. The electric field of BCL, which is a superposition of two circularly polarized light waves with frequencies that are integer multiples of each other, traces out a rose pattern in the polarization plane that can be chosen to break selective symmetries, including spatial inversion. Using a realistic low-energy model, we theoretically demonstrate that the three-dimensional Dirac semimetal Cd 3 As 2 is a promising platform for BCL Floquet engineering. Without strain, BCL irradiation induces a transition to a noncentrosymmetric magnetic Weyl semimetal phase with tunable energy separation between the Weyl nodes. In the presence of strain, we predict the emergence of a magnetic topological crystalline insulator with exotic unpinned surface Dirac states that are protected by a combination of twofold rotation and time reversal (2') and can be controlled by light.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Test of the Gravitational Force Law on Cosmological Scales Using the Kinematic Sunyaev-Zeldovich Effect

The mean pairwise velocity of massive halos reflects the gravitational force law on cosmic scales. For this work, we combine cosmic microwave background intensity maps from the Atacama Cosmology Telescope and a galaxy catalog from the Sloan Digital Sky Survey to estimate the mean pairwise velocity using the kinematic Sunyaev-Zeldovich (kSZ) effect. On scales from 30 to 230 megaparsecs, we constrain the gravitational acceleration between pairs of halos at separation 𝑟 to be 𝑔 ∝ 1/𝑟 𝑛 with 𝑛 = 2.1 ± 0.3, which is consistent with Newtonian gravity in an expanding spacetime (i.e., the standard Λ⁢ CDM model). This constraint shows agreement with an inverse quadratic radial dependence over the large distances separating galaxy halos, as expected in standard cosmology. Upcoming surveys have the potential to rule out 𝑛 = 1 at 10⁢𝜎 significance. Our results establish the kSZ effect as a powerful tool for testing gravity on cosmological scales.

alternative gravity theories↗
DOE OSTI2020

Utah FORGE: Phase 3 Magnetotelluric (MT) Data

New high-quality tensor MT data at 122 sites, including the vertical magnetic field and utilizing ultra-remote referencing, have been acquired over the Utah FORGE project area. The results will be used to delineate the densities of faults and fractures in crystalline basement rocks, to define the heat sources, and to derive baseline 3D resistivity structure for later MT monitoring of temporal changes in resistivity structure following well stimulation in the EGS reservoir. There are three files here related to Utah FORGE magnetotelluric (MT) data acquisition and processing. The FORGE MT EDIs zip file contains the observed MT responses in industry-standard EDI format. For each site, there is an EDI response file that utilized a local independent reference for noise cancellation, and a file that utilized a distant reference cancelling noise associated with the DC transmission line of the Delta IPP passing down the west side of Milford Valley. These two site files could be merged as appropriate. The FORGE Model Cell Center file contains the model volume of the 3D Forge MT inversion for characterizing the resistivity structure in the project area. It was derived using finite element inversion methodology described in Wannamaker et al, in the attached FORGE Phase3 Geophysics paper, from the MT observation EDI files. It is ASCII format (.dat) and entries are defined at top of the file in a simple x-y-z-Rho listing in UTM coordinates. The element layers drape the topography so the Rho value layers are not purely horizontal slices. This greatly simplifies the listing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Scientific Computational Imaging Code (SCICO)

Scientific Computational Imaging Code (SCICO) is a Python package for solving the inverse problems that arise in scientific imaging applications. Its primary focus is providing methods for solving ill-posed inverse problems by using an appropriate prior model of the reconstruction space. SCICO includes a growing suite of operators, cost functionals, regularizers, and optimization routines that may be combined to solve a wide range of problems, and is designed so that it is easy to add new building blocks. SCICO is built on top of JAX rather than NumPy, enabling GPU/TPU acceleration, just-in-time compilation, and automatic gradient functionality, which is used to automatically compute the adjoints of linear operators. An example of how to solve a multi-channel tomography problem with SCICO is shown in Figure 1. The SCICO source code is available from GitHub, and pre-built packages are available from PyPI. It has extensive online documentation, including API documentation and usage examples, which can be run online at Google Colab and binder.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Characterization of Partially Observed Epidemics - Application to COVID-19

This report documents a statistical method for the "real-time" characterization of partially observed epidemics. Observations consist of daily counts of symptomatic patients, diagnosed with the disease. Characterization, in this context, refers to estimation of epidemiological parameters that can be used to provide short-term forecasts of the ongoing epidemic, as well as to provide gross information for the time-dependent infection rate. The characterization problem is formulated as a Bayesian inverse problem, and is predicated on a model for the distribution of the incubation period. The model parameters are estimated as distributions using a Markov Chain Monte Carlo (MCMC) method, thus quantifying the uncertainty in the estimates. The method is applied to the COVID-19 pandemic of 2020, using data at the country, provincial (e.g., states) and regional (e.g. county) levels. The epidemiological model includes a stochastic component due to uncertainties in the incubation period. This model-form uncertainty is accommodated by a pseudo-marginal Metropolis-Hastings MCMC sampler, which produces posterior distributions that reflect this uncertainty. We approximate the discrepancy between the data and the epidemiological model using Gaussian and negative binomial error models; the latter was motivated by the over-dispersed count data. For small daily counts we find the performance of the calibrated models to be similar for the two error models. For large daily counts the negative-binomial approximation is numerically unstable unlike the Gaussian error model. Application of the model at the country level (for the United States, Germany, Italy, etc.) generally provided accurate forecasts, as the data consisted of large counts which suppressed the day-to-day variations in the observations. Further, the bulk of the data is sourced over the duration before the relaxation of the curbs on population mixing, and is not confounded by any discernible country-wide second wave of infections. At the state-level, where reporting was poor or which evinced few infections (e.g., New Mexico), the variance in the data posed some, though not insurmountable, difficulties, and forecasts were able to capture the data with large uncertainty bounds. The method was found to be sufficiently sensitive to discern the flattening of the infection and epidemic curve due to shelter-in-place orders after around 90% quantile for the incubation distribution (about 10 days for COVID-19). The proposed model was also used at a regional level to compare the forecasts for the central and north-west regions of New Mexico. Modeling the data for these regions illustrated different disease spread dynamics captured by the model. While in the central region the daily counts peaked in the late April, in the north-west region the ramp-up continued for approximately three more weeks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluation of Electrical Resistivity Tomography to Monitor the Transport of Past Releases Beneath Tank Farms

Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗