Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion constraints”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Visual display aid for orbital maneuvering - Design considerations

This paper describes the development of an interactive proximity operations planning system that allows on-site planning of fuel-efficient multiburn maneuvers in a potential multispacecraft environment. Although this display system most directly assists planning by providing visual feedback to aid visualization of the trajectories and constraints, its most significant features include: (1) the use of an 'inverse dynamics' algorithm that removes control nonlinearities facing the operator, and (2) a trajectory planning technique that separates, through a 'geometric spreadsheet', the normally coupled complex problems of planning orbital maneuvers and allows solution by an iterative sequence of simple independent actions. The visual feedback of trajectory shapes and operational constraints, provided by user-transparent and continuously active background computations, allows the operator to make fast, iterative design changes that rapidly converge to fuel-efficient solutions. The planning tool provides an example of operator-assisted optimization of nonlinear cost functions.

Grunwald, Arthur J.↗
DOE OSTI2021

NuSTAR Observations of Abell 2163: Constraints on Non-thermal Emission

Since the first non-thermal reports of inverse Compton (IC) emission from the intracluster medium (ICM) of galaxy clusters at hard X-ray energies, we have yet to unambiguously confirm IC emission in observations with newer facilities. RXTE detected IC emission in one of the hottest known clusters, Abell 2163 (A2163), a massive merging cluster with a giant radio halo—the presumed source of relativistic electrons IC scattering CMB photons to X-ray energies. The cluster's redshift (z ~ 0.2) allows its thermal and non-thermal radio emission to fit the field of view of the Nuclear Spectroscopic Telescope Array (NuSTAR), permitting a deep observation capable of confirming or ruling out the RXTE report. The IC flux provides constraints on the average magnetic field strength in a cluster. To determine the global diffuse IC emission in A2163, we fit its global NuSTAR spectrum with four models: single- (1T) and two-temperature (2T), 1T+power-law component (T+IC), and multi-temperature+power law (9T+IC). Each represents different characterizations of the thermal ICM emission, with power-law components added to represent IC emission. We find that the 3–30 keV spectrum can be described by purely thermal emission, with a global average temperature of kT = (11.8 ± 0.2) keV. The IC flux is constrained to <4.0 × 10 –12 erg s –1 cm –2 using the 1T+IC model and <1.6 × 10 –12 erg s –1 cm –2 with the more physical 9T+IC model, both to 90% confidence levels. Here, combining these limits with 1.4 GHz diffuse radio data from the VLA, we find the average magnetic field strength to be >0.22 μG and >0.35 μG, respectively, providing the strongest constraints on these values in A2163 to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Visualization and efficient generation of constrained high-dimensional theoretical parameter spaces

We describe a set of novel methods for efficiently sampling high-dimensional parameter spaces of physical theories defined at high energies, but constrained by experimental measurements made at lower energies. Often, theoretical models such as supersymmetry are defined by many parameters, $\mathcal{O}$(10 – 100), expressed at high energies, while relevant experimental constraints are often defined at much lower energies, preventing them from directly ruling out portions of the space. Instead, the low-energy constraints define a complex, potentially non-contiguous subspace of the theory parameters. Naive scanning of the theory space for points which satisfy the low-energy constraints is hopelessly inefficient due to the high dimensionality, and the inverse problem is considered intractable. As a result, many theoretical spaces remain under-explored. We introduce a class of modified generative autoencoders, which attack this problem by mapping the high-dimensional parameter space to a structured low-dimensional latent space, allowing for easy visualization and efficient generation of theory points which satisfy experimental constraints. An extension without dimensional compression, which focuses on limiting potential information loss, is also introduced.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2022

Employing Relaxed Smoothness Constraints on Imaginary Part of Refractive Index in AERONET Aerosol Retrieval Algorithm

In the Aerosol Robotic Network (AERONET) retrieval algorithm, smoothness constraints on the imaginary part of the refractive index provide control of retrieved spectral dependence of aerosol absorption by preventing the inversion code from fitting the noise in optical measurements and thus avoiding unrealistic oscillations of retrievals with wavelength. The history of implementation of the smoothness constraints in the AERONET retrieval algorithm is discussed. It is shown that the latest version of the smoothness constraints on the imaginary part of refractive index, termed standard and employed by Version 3 of the retrieval algorithm, should be modified to account for strong variability of light absorption by brown-carbon-containing aerosols in UV through mid-visible parts of the solar spectrum. In Version 3 strong spectral constraints were imposed at high values of the Ångström exponent (440–870 nm) since black carbon was assumed to be the primary absorber, while the constraints became increasingly relaxed as aerosol exponent deceased to allow for wavelength dependence of absorption for dust aerosols. The new version of the smoothness constraints on the imaginary part of the refractive index assigns different weights to different pairs of wavelengths, which are the same for all values of the Ångström exponent. For example, in the case of four-wavelength input, the weights assigned to short-wavelength pairs (440–675, 675–870 nm) are small so that smoothness constraints do not suppress natural spectral variability of the imaginary part of the refractive index. At longer wavelengths (870–1020 nm), however, the weight is 10 times higher to provide additional constraints on the imaginary part of refractive index retrievals of aerosols with a high Ångström exponent due to low sensitivity to aerosol absorption for longer channels at relatively low aerosol optical depths. The effect of applying the new version of smoothness constraints, termed relaxed, on retrievals of single-scattering albedo is analyzed for case studies of different aerosol types: black- and brown-carbon-containing fine mode aerosols, mineral dust coarse mode aerosols, and urban industrial fine mode aerosol. It is shown that for brown-carbon-containing aerosols employing the relaxed smoothness constraints resulted in significant reduction in retrieved single-scattering albedo and spectral residual errors (compared to standard) at the short wavelengths. For example, biomass burning smoke cases showed a reduction in single-scattering albedo and spectral residual error at 380 nm of ∼ 0.033 and ∼ 17 %, respectively, for the Rexburg site and ∼ 0.04 and ∼ 12.7 % for the Rimrock site, both AERONET sites in Idaho, USA. For a site with very high levels of black-carbon-containing aerosols (Mongu, Zambia), the effect of modification in the smoothness constraints was minor. For mineral dust aerosols at small Ångström exponent values (Mezaira site, UAE), the spectral constraint on the imaginary part of the refractive index was already relaxed in Version 3; therefore the new relaxed constraint results in minimal change. In the case of weakly absorbing urban industrial aerosols at the GSFC site, there are significant changes in retrieved single-scattering albedo using relaxed assumption, especially reductions at longer wavelengths: ∼ 0.016 and ∼ 0.02 at 875 and 1020 nm, respectively, for 440 nm aerosol optical depth (AOD) ∼ 0.3. The modification of smoothness constraints on the imaginary part of the refractive index has a minor effect on retrievals of other aerosol parameters such as the real part of the refractive index and parameters of the aerosol size distribution. The implementation of the relaxed smoothness constraints on the imaginary part of the refractive index in the next version of the AERONET inversion algorithm will produce significant impacts at some sites in short wavelength channels (380 and 440 nm) for some biomass burning smoke cases with significant brown carbon content and possibly in mid-visible channels (500 and 675 nm) to near-infrared channels (870 to 1020 nm) for some urban industrial aerosol types. However, most differences in single-scattering albedo retrievals between those applying the new relaxed constraint and the standard constraint will be within the uncertainty of the single-scattering albedo retrievals, depending on the level of aerosol optical depth, Ångström exponent, brown carbon content and wavelength.

Aliaksandr Siniuk↗
NASA NTRS2002

Site distribution and aliasing effects in the inversion for load coefficients and geocenter motion from GPS data

Precise GPS measurements of elastic relative site displacements due to surface mass loading offer important constraints on global surface mass transport. We investigate effects of site distribution and aliasing by higher- degree (n ≥ 2) loading terms on inversion of GPS data for n = 1 load coefficients and geocenter motion. Covariance and simulation analyses are conducted to assess the sensitivity of the inversion to aliasing and mismodeling errors and possible uncertainties in the n = 1 load coefficient determination. We found that the use of center-of-figure approximation in the inverse formulation could cause 10– 15% errors in the inverted load coefficients. n = 1 load estimates may be contaminated significantly by unknown higher-degree terms, depending on the load scenario and the GPS site distribution. The uncertainty in n = 1 zonal load estimate is at the level of 80–95% for two load scenarios.

Webb, Frank H.↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

DIGS: deep inference of galaxy spectra with neural posterior estimation

Abstract With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of simulation-based inference (SBI) and amortized neural posterior estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is (a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; and (b) amortized calculations of posterior distributions of said galaxy parameters at the modest cost of a few galaxy fits with Markov chain Monte Carlo (MCMC) methods. We utilise the SED generator and inference framework Prospector to generate simulated spectra, and train a dataset of 2 × 10 6 spectra (corresponding to a five-parameter SED model) with NPE. We show that SBI—with its combination of fast and amortized posterior estimations—is capable of inferring accurate galaxy stellar masses and metallicities. Our uncertainty constraints are comparable to or moderately weaker than traditional inverse-modelling with Bayesian MCMC methods (e.g. 0.17 and 0.26 dex in stellar mass and metallicity for a given galaxy, respectively). We also find that our inference framework conducts rapid SED inference (0.9–1.2 × 10 5 galaxy spectra via SBI/NPE at the cost of 1 MCMC-based fit). With this work, we set the stage for further work that focuses of SED fitting of galaxy spectra with SBI, in the era of JWST galaxy survey programs and the wide-field Roman Space Telescope spectroscopic surveys.

spectroscopy↗
NASA NTRS2012

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools.

Anderson, George R.↗
DOE OSTI2022

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Site Distribution and Aliasing Effects in the Inversion for Load Coefficients and Geocenter Motion from GPS Data

Precise GPS measurements of elastic relative site displacements due to surface mass loading offer important constraints on global surface mass transport. We investigate effects of site distribution and aliasing by higher-degree (n greater than or equal 2) loading terms on inversion of GPS data for n = 1 load coefficients and geocenter motion. Covariance and simulation analyses are conducted to assess the sensitivity of the inversion to aliasing and mismodeling errors and possible uncertainties in the n = 1 load coefficient determination. We found that the use of center-of-figure approximation in the inverse formulation could cause 10- 15% errors in the inverted load coefficients. n = 1 load estimates may be contaminated significantly by unknown higher-degree terms, depending on the load scenario and the GPS site distribution. The uncertainty in n = 1 zonal load estimate is at the level of 80 - 95% for two load scenarios.

global positioning system↗
DOE OSTI2021

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗
NASA NTRS1985

Algorithm Solves Constrained and Unconstrained Optimization Problems

Is quasi-Newton iteration utilizing Broyden/Fletcher/Goldfarb/Shanno update on inverse Hessian matrix. Capable of solving constrained optimization unconstrained optimization and constraints only problems with one to five independent variables from one to five constraint functions and one dependent function optimized.

Denson, M. A.↗
NASA NTRS1990

The shape of Eros

Monte Carlo simulations are presently used to optimize estimation, ascertain associated errors, and guide bias-correction procedures, for the Eros polar silhouette convex hull that has been estimated from radar echo spectra. This hull is trapezoidal; this nonaxisymmetric shape may account for odd harmonics in Eros' echo spectral signature as a function of rotation phase. Additional constraints have been obtained for the figure of Eros through the inversion of the optical lightcurve to estimate the asteroid's two-dimensional average of the three-dimensional shape. This 'mean cross-section' and the polar silhouette exhibit similar elongations.

Ostro, S. J.↗
NASA NTRS2002

Turbomachinery Airfoil Design Optimization Using Differential Evolution

An aerodynamic design optimization procedure that is based on a evolutionary algorithm known at Differential Evolution is described. Differential Evolution is a simple, fast, and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems, including highly nonlinear systems with discontinuities and multiple local optima. The method is combined with a Navier-Stokes solver that evaluates the various intermediate designs and provides inputs to the optimization procedure. An efficient constraint handling mechanism is also incorporated. Results are presented for the inverse design of a turbine airfoil from a modern jet engine and compared to earlier methods. The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated. Substantial reductions in the overall computing time requirements are achieved by using the algorithm in conjunction with neural networks.

Madavan, Nateri K.↗
NASA NTRS2002

Turbomachinery Airfoil Design Optimization Using Differential Evolution

An aerodynamic design optimization procedure that is based on a evolutionary algorithm known at Differential Evolution is described. Differential Evolution is a simple, fast, and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems, including highly nonlinear systems with discontinuities and multiple local optima. The method is combined with a Navier-Stokes solver that evaluates the various intermediate designs and provides inputs to the optimization procedure. An efficient constraint handling mechanism is also incorporated. Results are presented for the inverse design of a turbine airfoil from a modern jet engine. The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated. Substantial reductions in the overall computing time requirements are achieved by using the algorithm in conjunction with neural networks.

Madavan, Nateri K.↗
DOE OSTI2026

Search for New Physics via Low-Energy Electron Recoils with a 4.2 Tonne-Year Exposure from the LZ Experiment

We report results from searches for new physics models through electron recoils using data collected by the LUX-ZEPLIN experiment during its first two science runs, with a total exposure of 4.2 tonne−years. The observed data are consistent with a background-only hypothesis. Constraints are derived for electromagnetic interactions of solar neutrinos, solar axionlike particles (ALPs), mirror dark matter, and the absorption of bosonic dark matter candidates. The inverse Primakoff process for 57 Fe deexcitation solar ALPs is considered for the first time. These results represent the most stringent constraints to date on keV-scale Primakoff and 57 Fe solar ALPs, bosonic dark matter, mirror dark matter, and neutrino millicharge, while remaining competitive for the other signal models investigated.

Axion-like particles↗
DOE OSTI2023

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗
DOE OSTI2024

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

Here, we analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call “complexity reduction”. In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature T g of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2–6 nm at T g . Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

cancer↗