SEARCH · Engineering Papers

Results for “Inverse Optimization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72

All sources NASA NTRS DOE OSTI

DOE OSTI2024

All-optical image denoising using a diffractive visual processor

Abstract Image denoising, one of the essential inverse problems, targets to remove noise/artifacts from input images. In general, digital image denoising algorithms, executed on computers, present latency due to several iterations implemented in, e.g., graphics processing units (GPUs). While deep learning-enabled methods can operate non-iteratively, they also introduce latency and impose a significant computational burden, leading to increased power consumption. Here, we introduce an analog diffractive image denoiser to all-optically and non-iteratively clean various forms of noise and artifacts from input images – implemented at the speed of light propagation within a thin diffractive visual processor that axially spans <250 × λ, where λ is the wavelength of light. This all-optical image denoiser comprises passive transmissive layers optimized using deep learning to physically scatter the optical modes that represent various noise features, causing them to miss the output image Field-of-View (FoV) while retaining the object features of interest. Our results show that these diffractive denoisers can efficiently remove salt and pepper noise and image rendering-related spatial artifacts from input phase or intensity images while achieving an output power efficiency of ~30–40%. We experimentally demonstrated the effectiveness of this analog denoiser architecture using a 3D-printed diffractive visual processor operating at the terahertz spectrum. Owing to their speed, power-efficiency, and minimal computational overhead, all-optical diffractive denoisers can be transformative for various image display and projection systems, including, e.g., holographic displays.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Two-dimensional topological crystalline quantum spin Hall effect in transition metal intercalated compounds

While most of the two-dimensional (2D) topological crystalline insulators (TCIs) belong to group IV-VI narrow-band-gap semiconductors in a square lattice, in the present work we predict a TCI family based on transition metal intercalated compounds in a hexagonal lattice. First-principles calculations combined with a substrate-fixed globally optimal structural search technique show that a layer of Os prefers a uniform distribution between two graphene sheets. Band dispersion calculations reveal a Dirac point and a Dirac nodal ring near the Fermi level. The Dirac point is ascribed to the hybridization of e 2 and e 2 * orbitals, and the Dirac ring is formed due to dispersion of s and e 1 * orbitals. Upon inclusion of spin-orbit coupling, these Dirac states open topologically nontrivial local band gaps, which are characterized by nonzero mirror Chern numbers. The quantum spin Hall effect is also observed by integrating the spin Berry curvature in the Brillouin zone. In contrast to the 2D group IV-VI TCIs whose band inversions at X and Y points are “locked” by C 4 rotation symmetry, here the relative energy of two local band gaps can be manipulated by in-plane biaxial strains. Some other similar intercalation compounds are also shown to be topologically nontrivial. Our work extends the 2D TCI family into a hexagonal lattice composed of transition metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

DECT-MULTRA: Dual-Energy CT Image Decomposition With Learned Mixed Material Models and Efficient Clustering

Dual-energy computed tomography (DECT) imaging plays an important role in advanced imaging applications due to its material decomposition capability. Image-domain decomposition operates directly on CT images using linear matrix inversion, but the decomposed material images can be severely degraded by noise and artifacts. This paper proposes a new method dubbed DECT-MULTRA for image-domain DECT material decomposition that combines conventional penalized weighted-least squares (PWLS) estimation with regularization based on a mixed union of learned transforms (MULTRA) model. Our proposed approach pre-learns a union of common-material sparsifying transforms from patches extracted from all the basis materials, and a union of cross-material sparsifying transforms from multi-material patches. The common-material transforms capture the common properties among different material images, while the cross-material transforms capture the cross-dependencies. The proposed PWLS formulation is optimized efficiently by alternating between an image update step and a sparse coding and clustering step, with both of these steps having closed-form solutions. The effectiveness of our method is validated with both XCAT phantom and clinical head data. The results demonstrate that our proposed method provides superior material image quality and decomposition accuracy compared to other competing methods.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Imaging of Calibrated Defects in Additively Manufactured Materials: Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components

Additive manufacturing (AM, or 3D printing) for commercial nuclear energy applications is an emerging method for cost-efficient manufacturing aimed at replacing aging nuclear reactor parts and reducing costs for new construction. Known material flaws in AM include low-density regions consisting of non-sintered powder, which have to be detected to ensure the safety of long-term performance nuclear reactor components. Currently, limited options are available for nondestructive evaluation (NDE), either during or post manufacturing. As a solution to NDE of AM, we are developing pulsed thermal imaging which is non-contact, one-sided, and scalable to arbitrary size and shape of the AM parts. Pulsed thermography (PT) system utilizes a high intensity flash lamp to rapidly heat surface of sample, and a high sensitivity fast-frame megapixel infrared (IR) camera to capture data of surface temperature variations. The acquired data cube consists of a stack of surface temperature images taken at different times. Information about material internal defects is extracted by analyzing the data cube. This report provides results of preliminary performance evaluation of pulsed thermal imaging capability in detection of imprinted flaws in AM metallic structures. The flaws were introduced into AM parts as imprinted hemispherical low density regions, consisting of trapped un-sintered metallic powder. Specimens for developed for this study consisted of AM stainless steel 316 and Inconel 718 plates. The diameters of imprinted defects varied from 1mm to 8mm, and their depths below the plate flat surface varied between 1mm and 6mm. Pulsed thermal tomography (PTT) processes the measured data cube to obtain 3D reconstructions of material effusivity using a unique inversion algorithm developed at Argonne. PTT has been previously used in imaging of similar size flat bottom hole (FBH) simulated defects in stainless steel 316 and Inconel 718 specimens. In the study involving AM specimens, PTT imaging results have shown that 1mm-diameter defects located 1mm and 2mm below the surface of specimens were detectable. Larger size defects were detectable at greater depth. We also explored an alternative approach to detection of material flaws in PT data cube, which is using neural learning-based approach to blind source separation. Detection of small material defects requires finding features in the data cube which have signal contrast levels approaching sensitivity limit of IR camera. In this study, an optimized Neural Learning based Blind Source Separation (NLBSS) algorithm, including Principal Component Analysis (PCA), and Independent Component Analysis (ICA) is demonstrated to automatically extract principal temporal and spatial features of thermography frames to enhance flaw detection. By using the NLBSS algorithm, material internal defects can be automatically detected. Furthermore, this processing approach compensates for experimental thermal imaging artifacts, such as noise and uneven heating. By merging artificial intelligence with phtotothermics, the NDE system detects internal calibrated defects of various sizes and depths in AM nuclear-grade metallic alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2026

The NREL Sensor Laboratory Detection of Hydrogen Emissions

The development of a functional hydrogen detection system is a multifaceted process that integrates hardware, deployments strategies, and analytics which can be supported by the NREL Sensor Laboratory: 1. Support of the design, validation and optimization of sensing prototypes; 2. Guide optimized sensing element development, including control electronics; 3. Laboratory testing to validate/optimize metrological performance (measurement range, detection limit, etc.); 4. Provide test sites for field deployments representative of real-world scenarios with controlled hydrogen releases; 5. Develop sensor placement and operation guidance; 6. Provide guidance on electronics to accommodate facility integration; 7. Electrical safety designs to allow for operation within restricted zones; 8. Integration into facility monitoring and control systems; 9. Guide incorporation of cyber security elements to protect facilities from malicious attacks; 10. Modeling and application of advanced analytics to detect and quantify emissions; 11. Higher Order dispersion models to guide sensor placement for reliable detection; 12. Advanced analytics for improved metrological performances, and to inform inverse modeling; 13. Market support and commercialization (national and international markets); 14. Commercial deployments in H2@SCALE markets (e.g., HUBs and other large-scale hydrogen markets); and 15. Leverage off international collaborations/partnerships (e.g., NREL is on the advisory board for the European initiative "pre-Normative Research on Hydrogen Releases Assessment"-NHyRA).

08 HYDROGEN↗

NASA NTRS2023

Using Orbiting Carbon Observatory-2 (OCO-2) column CO2 retrievals to rapidly detect and estimate biospheric surface carbon flux anomalies

The global carbon cycle is experiencing continued perturbations via increases in atmospheric carbon concentrations, which are partly reduced by terrestrial biosphere and ocean carbon uptake. Greenhouse gas satellites have been shown to be useful in retrieving atmospheric carbon concentrations and observing surface and atmospheric CO2 seasonal-to-interannual variations. However, limited attention has been placed on using satellite column CO2 retrievals to evaluate surface CO2 fluxes from the terrestrial biosphere without advanced inversion models at low latency. Such applications could be useful to monitor, in near real time, biosphere carbon fluxes during climatic anomalies like drought, heatwaves, and floods, before more complex terrestrial biosphere model outputs and/or advanced inversion modelling estimates become available. Here, we explore the ability of Orbiting Carbon Observatory-2 (OCO-2) column-averaged dry air CO2 (XCO2) retrievals to directly detect and estimate terrestrial biosphere CO2 flux anomalies using a simple mass-balance approach. An initial global analysis of surface–atmospheric CO2 coupling and transport conditions reveals that the western US, among a handful of other regions, is a feasible candidate for using XCO2 for detecting terrestrial biosphere CO2 flux anomalies. Using the CarbonTracker model reanalysis as a test bed, we first demonstrate that a well-established mass-balance approach can estimate monthly surface CO2 flux anomalies from XCO2 enhancements in the western United States. The method is optimal when the study domain is spatially extensive enough to account for atmospheric mixing and has favorable advection conditions with contributions primarily from one background region. We find that errors in individual soundings reduce the ability of OCO-2 XCO2 to estimate more frequent, smaller surface CO2 flux anomalies. However, we find that OCO-2 XCO2 can often detect and estimate large surface flux anomalies that leave an imprint on the atmospheric CO2 concentration anomalies beyond the retrieval error/uncertainty associated with the observations. OCO-2 can thus be useful for low-latency monitoring of the monthly timing and magnitude of extreme regional terrestrial biosphere carbon anomalies.

Andrew F. Feldman↗

DOE OSTI2016

Single-view phase retrieval of an extended sample by exploiting edge detection and sparsity

We propose a new approach to robustly retrieve the exit wave of an extended sample from its coherent diffraction pattern by exploiting sparsity of the sample's edges. This approach enables imaging of an extended sample with a single view, without ptychography. We introduce nonlinear optimization methods that promote sparsity, and we derive update rules to robustly recover the sample's exit wave. We test these methods on simulated samples by varying the sparsity of the edge-detected representation of the exit wave. Finally, our tests illustrate the strengths and limitations of the proposed method in imaging extended samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1976

Elliptic integral solutions to a class of space flight optimization problems

This paper is initially concerned with the minimum-time, exoatmospheric flight of a rocket with constant thrust acceleration magnitude, as in the cases of nuclear and solar electric propulsion. Gravitational acceleration is assumed to be a constant scalar multiple of the radius vector, plus a correction term which is a given function of time. The solution to the state equations is obtained in terms of elliptic integrals. A method is presented for the solution of the two-point boundary-condition problem associated with orbital transfer. At most, the latter method requires iteration upon final time, angle of injection, and two other parameters which are bounded. An example problem is provided which involves a rocket with very low thrust and a spiraling trajectory of many revolutions, but an altitude change of only several hundred miles above the earth. Finally, the original elliptic integral solution is extended to a larger class of low and intermediate thrust problems with constant thrust magnitude, mass decreasing with time, and an inverse square gravitational force.

Andrus, J. F.↗

DOE OSTI2024

Stomatal Parameters in a Changing Environment

Here, we recommend that stomatal slope parameters (g 1 ) be inferred by inversion so that variations in g 1 may be attributed to variations physiological and environmental conditions. Understanding g1 will advance predictions of plant gas exchange and performance under global climate.

Stomatal conductance modeling↗

NASA NTRS1994

Evaluating the design of an earth radiation budget instrument with system simulations. Part 2: Minimization of instantaneous sampling errors for CERES-I

Much of the new record of broadband earth radiation budget satellite measurements to be obtained during the late 1990s and early twenty-first century will come from the dual-radiometer Clouds and Earth's Radiant Energy System Instrument (CERES-I) flown aboard sun-synchronous polar orbiters. Simulation studies conducted in this work for an early afternoon satellite orbit indicate that spatial root-mean-square (rms) sampling errors of instantaneous CERES-I shortwave flux estimates will range from about 8.5 to 14.0 W/m on a 2.5 deg latitude and longitude grid resolution. Rms errors in longwave flux estimates are only about 20% as large and range from 1.5 to 3.5 W/sq m. These results are based on an optimal cross-track scanner design that includes 50% footprint overlap to eliminate gaps in the top-of-the-atmosphere coverage, and a 'smallest' footprint size to increase the ratio in the number of observations lying within to the number of observations lying on grid area boundaries. Total instantaneous measurement error also depends on the variability of anisotropic reflectance and emission patterns and on retrieval methods used to generate target area fluxes. Three retrieval procedures from both CERES-I scanners (cross-track and rotating azimuth plane) are used. (1) The baseline Earth Radiaton Budget Experiment (ERBE) procedure, which assumes that errors due to the use of mean angular dependence models (ADMs) in the radiance-to-flux inversion process nearly cancel when averaged over grid areas. (2) To estimate N, instantaneous ADMs are estimated from the multiangular, collocated observations of the two scanners. These observed models replace the mean models in computation of satellite flux estimates. (3) The scene flux approach, conducts separate target-area retrievals for each ERBE scene category and combines their results using area weighting by scene type. The ERBE retrieval performs best when the simulated radiance field departs from the ERBE mean models by less than 10%. For larger perturbations, both the scene flux and collocation methods produce less error than the ERBE retrieval. The scene flux technique is preferable, however, because it involves fewer restrictive assumptions.

Stowe, Larry↗

NASA NTRSDate not supplied

Optimization of a Parametric Lobe Mixer for Internally Mixed Nozzles

The results of the first phase of an optimization study on a parametric procedurally-generated lobe mixer are presented. This mixer is part of a dual-stream plug nozzle system, designed under NASA’s CST project for use in near-term commercial supersonic applications. Data was obtained numerically using NASA’s in-house CFD code FUN3D and analyzed to determine mixer performance. Phase I of the study examined five independently-variable geometric parameters describing the lobe mixer. Methods of Design of Experiments (DoE) were used to obtain a parametric model to quickly evaluate mixer performance and determine an optimum configuration. Mixers that induced strong axial vorticity and shear layers without causing flow separation were found to perform best, with exit temperature variance highly inversely correlated to mixing effectiveness.

Lobe Mixer↗

NASA NTRSDate not supplied

Optimization of a Parametric Lobe Mixer for Internally Mixed Nozzles

Lobe mixer↗

DOE OSTI2024

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗

DOE OSTI2014

Chemical kinetics parameters and model validation for the gasification of PCEA nuclear graphite

A series of gasification experiments, using two right cylinder specimens (~ 12.7 x 25.4 mm and 25.4 x 25.4 mm) of PCEA nuclear graphite in ambient airflow, measured the total gasification flux at weight losses up to 41.5% and temperatures (893-1015 K) characteristics of those for in-pores gasification Mode (a) and in-pores diffusion-limited Mode (b). The chemical kinetics parameters for the gasification of PCEA graphite are determined using a multi-parameters optimization algorithm from the measurements of the total gasification rate and transient weight loss in experiments. These parameters are: (i) the pre-exponential rate coefficients and the Gaussian distributions and values of specific activation energies for adsorption of oxygen and desorption of CO gas; (ii) the specific activation energy and pre-exponential rate coefficient for the breakup of stable un-dissociated C(O2) oxygen radicals to form stable (CO) complexes; (iii) the specific activation energy and pre-exponential coefficient for desorption of CO2 gas and; (iv) the initial surface area of reactive free sites per unit mass. This area is consistently 13.5% higher than that for nuclear graphite grades of NBG-25 and IG-110 and decreases inversely proportional with the square root of the initial mass of the graphite specimens in the experiments. Experimental measurements successfully validate the chemical-reactions kinetics model that calculates continuous Arrhenius curves of the total gasification flux and the production rates of CO and CO2 gases. The model results at different total weight losses agree well with measurements and expand beyond the temperatures in the experiments to the diffusion-limited mode of gasification. Also calculated are the production rates of CO and CO2 gases and their relative contributions to the total gasification rate in the experiments as functions of temperature, for total weight losses of 5% and 10%.

kinetics↗

DOE OSTI2014

Risk Analysis Virtual ENvironment

RAVEN has 3 major functionalities: 1. Provides a Graphical User Interface for the pre- and post-processing of the RELAP-7 input and output. 2. Provides the capability to model nuclear power plants control logic for the RELAP-7 code and dynamic control of the accident scenario evolution. This capability is based on a software structure that realizes a direct connection between the RELAP-7 solver engine (MOOSE) and a python environment where the variables describing the plant status are accessible in a scripting environment. RAVEN support the generation of the probabilistic scenario control by supplying a wide range of probability and cumulative distribution functions and their inverse functions. 3. Provides a general environment to perform probability risk analysis for RELAP-7, RELAP-5 and any generic MOOSE based applications. The probabilistic analysis is performed by sampling the input space of the coupled code parameters and it is enhanced by using modern artificial intelligence algorithms that accelerate the identification of the areas of major risk (in the input parameter space). This environment also provides a graphical visualization capability to analyze the outcomes. Among other approaches, the classical Monte Carlo and Latin Hypercube sampling algorithms are available. For the acceleration of the convergence of the sampling methodologies, Support Vector Machines, Bayesian regression, and collocation stochastic polynomials chaos are implemented. The same methodologies here described could be used to solve optimization and uncertainties propagation problems using the RAVEN framework.

Rabiti, Cristian↗

DOE OSTI2019

Setting Upper Limits on the Local Burst Rate Density of Primordial Black Holes Using HAWC

Primordial Black Holes (PBHs) are black holes that may have been created by density fluctuations in the early Universe and could be as massive as supermassive black holes or as small as the Planck scale. It is believed that a black hole has a temperature inversely proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of ~5e14 g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV – TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to particles with energies of 300 GeV – 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field of view of ~2 sr and a duty cycle above 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search on 959 days of gamma-ray data from HAWC by optimizing a previous gamma-ray burst transient search to the PBH burst energy spectrum, and placed the strongest upper limits on the local PBH burst rate density at the 99% confidence level.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1990

Recursive dynamics for flexible multibody systems using spatial operators

Due to their structural flexibility, spacecraft and space manipulators are multibody systems with complex dynamics and possess a large number of degrees of freedom. Here the spatial operator algebra methodology is used to develop a new dynamics formulation and spatially recursive algorithms for such flexible multibody systems. A key feature of the formulation is that the operator description of the flexible system dynamics is identical in form to the corresponding operator description of the dynamics of rigid multibody systems. A significant advantage of this unifying approach is that it allows ideas and techniques for rigid multibody systems to be easily applied to flexible multibody systems. The algorithms use standard finite-element and assumed modes models for the individual body deformation. A Newton-Euler Operator Factorization of the mass matrix of the multibody system is first developed. It forms the basis for recursive algorithms such as for the inverse dynamics, the computation of the mass matrix, and the composite body forward dynamics for the system. Subsequently, an alternative Innovations Operator Factorization of the mass matrix, each of whose factors is invertible, is developed. It leads to an operator expression for the inverse of the mass matrix, and forms the basis for the recursive articulated body forward dynamics algorithm for the flexible multibody system. For simplicity, most of the development here focuses on serial chain multibody systems. However, extensions of the algorithms to general topology flexible multibody systems are described. While the computational cost of the algorithms depends on factors such as the topology and the amount of flexibility in the multibody system, in general, it appears that in contrast to the rigid multibody case, the articulated body forward dynamics algorithm is the more efficient algorithm for flexible multibody systems containing even a small number of flexible bodies. The variety of algorithms described here permits a user to choose the algorithm which is optimal for the multibody system at hand. The availability of a number of algorithms is even more important for real-time applications, where implementation on parallel processors or custom computing hardware is often necessary to maximize speed.

Jain, A.↗

DOE OSTI2016

An Optimal Parameterization Framework for Infrasonic Tomography of the Stratospheric Winds Using Non-Local Sources

A method is developed to apply acoustic tomography methods to a localized network of infrasound arrays with intention of monitoring the atmosphere state in the region around the network using non-local sources without requiring knowledge of the precise source location or non-local atmosphere state. Closely spaced arrays provide a means to estimate phase velocities of signals that can provide limiting bounds on certain characteristics of the atmosphere. Larger spacing between such clusters provide a means to estimate celerity from propagation times along multiple unique stratospherically or thermospherically ducted propagation paths and compute more precise estimates of the atmosphere state. In order to avoid the commonly encountered complex, multimodal distributions for parametric atmosphere descriptions and to maximize the computational efficiency of the method, an optimal parametrization framework is constructed. This framework identifies the ideal combination of parameters for tomography studies in specific regions of the atmosphere and statistical model selection analysis shows that high quality corrections to the middle atmosphere winds can be obtained using as few as three parameters. Lastly, comparison of the resulting estimates for synthetic data sets shows qualitative agreement between the middle atmosphere winds and those estimated from infrasonic traveltime observations.

58 GEOSCIENCES↗