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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1996

An Inversion of Gravity and Topography for Mantle and Crustal Structure on Mars

Analysis of the gravity and topography of Mars presently provides our primary quantitative constraints on the internal structure of Mars. We present an inversion of the long-wavelength (harmonic degree less than or equal to 10) gravity and topography of Mars for lateral variations of mantle temperature and crustal thickness. Our formulation incorporates both viscous mantle flow (which most prior studies have neglected) and isostatically compensated density anomalies in the crust and lithosphere. Our nominal model has a 150-km-thick high-viscosity surface layer over an isoviscous mantle, with a core radius of 1840 km. It predicts lateral temperature variations of up to a few hundred degrees Kelvin relative to the mean mantle temperature, with high temperature under Tharsis and to a lesser extent under Elysium and cool temperatures elsewhere. Surprisingly, the model predicts crustal thinning beneath Tharsis. If correct, this implies that thinning of the crust by mantle shear stresses dominates over thickening of the crust by volcanism. The major impact basins (Hellas, Argyre, Isidis, Chryse, and Utopia) are regions of crustal thinning, as expected. Utopia is also predicted to be a region of hot mantle, which is hard to reconcile with the surface geology. An alternative model for Utopia treats it as a mascon basin. The Utopia gravity anomaly is consistent with the presence of a 1.2 to 1.6 km thick layer of uncompensated basalt, in good agreement with geologic arguments about the amount of volcanic fill in this area. The mantle thermal structure is the dominant contributor to the observed geoid in our inversion. The mantle also dominates the topography at the longest wavelengths, but shorter wavelengths (harmonic degrees greater than or equal to 4) are dominated by the crustal structure. Because of the uncertainty about the appropriate numerical values for some of the model's input parameters, we have examined the sensitivity of the model results to the planetary structural model (core radius and core and mantle densities), the mantle's viscosity stratification, and the mean crustal thickness. The model results are insensitive to the specific thickness or viscosity contrast of the high-viscosity surface layer and to the mean crustal thickness in the range 25 to 100 km. Models with a large core radius or with an upper mantle low-viscosity zone require implausibly large lateral variations in mantle temperature.

Kiefer, Walter S.↗

NASA NTRS2021

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗

NASA NTRS2021

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗

NASA NTRS1991

An inverse method with regularity condition for transonic airfoil design

It is known from Lighthill's exact solution of the incompressible inverse problem that in the inverse design problem, the surface pressure distribution and the free stream speed cannot both be prescribed independently. This implies the existence of a constraint on the prescribed pressure distribution. The same constraint exists at compressible speeds. Presented here is an inverse design method for transonic airfoils. In this method, the target pressure distribution contains a free parameter that is adjusted during the computation to satisfy the regularity condition. Some design results are presented in order to demonstrate the capabilities of the method.

Zhu, Ziqiang↗

NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗

DOE OSTI2022

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

NASA NTRS2009

A Relationship Between Constraint and the Critical Crack Tip Opening Angle

Of the various approaches used to model and predict fracture, the Crack Tip Opening Angle (CTOA) fracture criterion has been successfully used for a wide range of two-dimensional thin-sheet and thin plate applications. As thicker structure is considered, modeling the full three-dimensional fracture process will become essential. This paper investigates relationships between the local CTOA evaluated along a three-dimensional crack front and the corresponding local constraint. Previously reported tunneling crack front shapes were measured during fracture by pausing each test and fatigue cycling the specimens to mark the crack surface. Finite element analyses were run to model the tunneling shape during fracture, with the analysis loading conditions duplicating those tests. The results show an inverse relationship between the critical fracture value and constraint which is valid both before maximum load and after maximum load.

Johnston, William M.↗

DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Exploring theoretical uncertainties in the hydrodynamic description of relativistic heavy-ion collisions

We explore theoretical uncertainties in the hydrodynamic description of relativistic heavy-ion collisions by examining the full nonlinear causality conditions and quantifying the second-order transport coefficients' role on flow observables. The causality conditions impose physical constraints on the maximum allowed values of inverse Reynolds numbers during the hydrodynamic evolution. Including additional second-order gradient terms in the Denicol-Niemi-Molnár-Rischke (DNMR) theory significantly shrinks the casual regions compared to those in the Israel-Stewart hydrodynamics. For Au + Au collisions, we find the variations of flow observables are small with and without imposing the necessary causality conditions, suggesting a robust extraction of the quark-gluon plasma's transport coefficients in previous model-to-data comparisons. However, sizable sensitivity is present in small p + Au collisions, which poses challenges to study the small systems' collectivity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1990

The dynamics of orbital maneuvering: Design and evaluation of a visual display aid for human controllers

An interactive proximity operations planning system, which allows on-site planning of fuel-efficient, multi-burn maneuvers in a potential multi-spacecraft environment was developed and tested. Though this display system most directly assists planning by providing visual feedback to aid visualization of the trajectories and constraints, its most significant features include an inverse dynamics algorithm that removes control nonlinearities facing the operator and a trajectory planning technique that reduces the order of control and creates, through a geometric spread-sheet the illusion of an inertially stable environment. This synthetic environment provides the user with control of relevant static and dynamic properties of way-points during small orbital changes allowing independent solutions to the normally coupled problems of orbital maneuvering. An experiment was carried out in which experienced operators were required to plan a trajectory to retrieve an object accidently separated from a dual-keel space station. The time required to plan these maneuvers was found to be predicted by the direction of the insertion thrust and did not depend on the point of separation from the space station.

Ellis, Stephen R.↗

NASA NTRS1989

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗

DOE OSTI2024

Constraining solar electron number density via neutrino flavor data at Borexino

Understanding the physics of the deep solar interior, and the more exotic environs of core-collapse supernovae (CCSN) and binary neutron-star (NS) mergers, is of keen interest in many avenues of research. To date, this physics is based largely on simulations via forward integration. While these simulations provide valuable constraints, it could be insightful to adopt the "inverse approach" as a point of comparison. Within this paradigm, parameters of the solar interior are not output based on an assumed model, but rather are inferred based on real data. Here, we take the specific case of solar electron number density, which historically is taken as output from the standard solar model. We show how one may arrive at an independent constraint on that density profile based on available neutrino flavor data from the Earth-based Borexino experiment. The inference technique's ability to offer a unique lens on physics can be extended to other datasets, and to analogous questions for CCSN and NS mergers, albeit with simulated data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1994

The control of satellites with microgravity constraints: The COMET Control System

The COMET attitude determination and control system, using inverse dynamics and a novel torque distribution/momentum management technique, has shown great flexibility, performance, and robustness. Three-axis control with two wheels is an inherent consequence of inverse dynamics control which allows for reduction in spacecraft weight and cost, or alternatively, provides a simple means of failure-redundancy for three-wheel spacecraft. The control system, without modification, has continued to perform well in spite of large changes in spacecraft mass properties and mission orbit altitude that have occurred during development. This flexibility has obviated imposition of early stringent ADACS design constraints and has greatly reduced commonly incurred ADACS modification costs and delay associated with program maturation.

Grossman, Walter↗

NASA NTRS1978

Organic chemistry on Titan

Observations of nonequilibrium phenomena on the Saturn satellite Titan indicate the occurrence of organic chemical evolution. Greenhouse and thermal inversion models of Titan's atmosphere provide environmental constraints within which various pathways for organic chemical synthesis are assessed. Experimental results and theoretical modeling studies suggest that the organic chemistry of the satellite may be dominated by two atmospheric processes: energetic-particle bombardment and photochemistry. Reactions initiated in various levels of the atmosphere by cosmic ray, Saturn wind, and solar wind particle bombardment of a CH4 - N2 atmospheric mixture can account for the C2-hydrocarbons, the UV-visible-absorbing stratospheric haze, and the reddish color of the satellite. Photochemical reactions of CH4 can also account for the presence of C2-hydrocarbons. In the lower Titan atmosphere, photochemical processes will be important if surface temperatures are sufficiently high for gaseous NH3 to exist. Hot H-atom reactions initiated by photo-dissociation of NH3 can couple the chemical reactions of NH3 and CH4 and produce organic matter.

Chang, S.↗

NASA NTRS1985

A review of problems and progress in studies of satellite magnetic anomalies

A review is conducted of studies performed during the Magsat project. The obtained data are considered, taking into account questions of data availability, aspects of orbit attitude determination, ionospheric noise, a field model, and an anomaly field presentation. Models for interpretation are discussed, giving attention to forward modeling, and equivalent layer inverse modeling. In an evaluation of rock property constraints, the magnetic bottom is discussed along with Curie points, metamorphism and magnetization, and the direction of magnetization.

Mayhew, M. A.↗

NASA NTRS1991

A synthetic environment for visualization and planning of orbital maneuvers

An interactive proximity operations planning system, which allows on-site planning of fuel-efficient, multi-burn maneuvers in a potential multi-space-craft environment has been developed and tested. This display system most directly assists planning by providing visual feedback in a synthetic virtual space that aids visualization of trajectories and their constraints. Its most significant features include (1) an 'inverse dynamics' algorithm that removes control nonlinearities facing the operator and (2) a stack-oriented action-editor that reduces the order of control and creates, through a 'geometric spreadsheet,' the illusion of an inertially stable environment. This synthetic environment provides the user with control of relevant static and dynamic properties of way-points during small orbital changes allowing independent solutions to otherwise coupled problems of orbital maneuvering.

Ellis, Stephen R.↗

NASA NTRS1996

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗