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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Putting all the X in one basket: Updated X-ray constraints on sub-GeV Dark Matter

Sub-GeV dark matter particles can annihilate or decay producing e ± pairs which upscatter the low-energy photon fields in the Galaxy and generate an X-ray emission (via the Inverse Compton effect). Here, using X-ray data from X MM -N EWTON , I NTEGRAL , N U S TAR and S UZAKU , we derive new constraints on this class of dark matter (DM). In the annihilation case, our new bounds are the strongest available for DM masses above 180 MeV, reaching $\langle$σv$\rangle$ ≲ 10 -28 cm 3 /s for m DM ≃ 1 GeV. In the decay case, our bounds are the strongest to date over a large fraction of the considered mass range, constraining τ ≳ 10 28 s for m DM ≃ 1 GeV and improving by up to 3 orders of magnitude upon existing limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Simplified, inverse, ejector design tool

A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and entrainment estimates subject to operational, acoustic, and design constraints. These constraints are manifested through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity, and static pressure matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with direct application to High Speed Civil Transport (HSCT) design studies.

Dechant, Lawrence J.↗
NASA NTRS1991

Operator-assisted planning and execution of proximity operations subject to operational constraints

Future multi-vehicle operations will involve multiple scenarios that will require a planning tool for the rapid, interactive creation of fuel-efficient trajectories. The planning process must deal with higher-order, non-linear processes involving dynamics that are often counter-intuitive. The optimization of resulting trajectories can be difficult to envision. An interaction proximity operations planning system is being developed to provide the operator with easily interpreted visual feedback of trajectories and constraints. This system is hosted on an IRIS 4D graphics platform and utilizes the Clohessy-Wiltshire equations. An inverse dynamics algorithm is used to remove non-linearities while the trajectory maneuvers are decoupled and separated in a geometric spreadsheet. The operator has direct control of the position and time of trajectory waypoints to achieve the desired end conditions. Graphics provide the operator with visualization of satisfying operational constraints such as structural clearance, plume impingement, approach velocity limits, and arrival or departure corridors. Primer vector theory is combined with graphical presentation to improve operator understanding of suggested automated system solutions and to allow the operator to review, edit, or provide corrective action to the trajectory plan.

Grunwald, Arthur J.↗
NASA NTRS1991

Reliable automatic plan selection for visual robotic positioning

Reliability analysis techniques developed for robotics are applied to a six-DOF (degree-of-freedom) visual positioning system. The task is to measure an object's pose and then grip the object reliably subject to time and accuracy constraints. Errors due to pixel truncation in the stadimetric vision system, forward kinematic and inverse kinematic uncertainties, and joint position and velocity noise are stochastically modeled. Simulation results for PUMA 500 kinematic and dynamic characteristics are presented.

Mcinroy, John E.↗
NASA NTRS2008

Estimating Small-Body Gravity Field from Shape Model and Navigation Data

This paper presents a method to model the external gravity field and to estimate the internal density variation of a small-body. We first discuss the modeling problem, where we assume the polyhedral shape and internal density distribution are given, and model the body interior using finite elements definitions, such as cubes and spheres. The gravitational attractions computed from these approaches are compared with the true uniform-density polyhedral attraction and the level of accuracies are presented. We then discuss the inverse problem where we assume the body shape, radiometric measurements, and a priori density constraints are given, and estimate the internal density variation by estimating the density of each finite element. The result shows that the accuracy of the estimated density variation can be significantly improved depending on the orbit altitude, finite-element resolution, and measurement accuracy.

Covariance Analysis↗
DOE OSTI2023

Physics‐Informed Machine Learning for Inverse Design of Optical Metamaterials

Optical metamaterials manipulate light through various confinement and scattering processes, offering unique advantages like high performance, small form factor and easy integration with semiconductor devices. However, designing metasurfaces with suitable optical responses for complex metamaterial systems remains challenging due to the exponentially growing computation cost and the ill‐posed nature of inverse problems. To expedite the computation for the inverse design of metasurfaces, a physics‐informed deep learning (DL) framework is used. A tandem DL architecture with physics‐based learning is used to select designs that are scientifically consistent, have low error in design prediction, and accurate reconstruction of optical responses. The authors focus on the inverse design of a representative plasmonic device and consider the prediction of design for the optical response of a single wavelength incident or a spectrum of wavelength in the visible light range. The physics‐based constraint is derived from solving the electromagnetic wave equations for a simplified homogenized model. The model converges with an accuracy up to 97% for inverse design prediction with the optical response for the visible light spectrum as input, and up to 96% for optical response of single wavelength of light as input, with optical response reconstruction accuracy of 99%.

Sarkar, Sulagna↗
NASA NTRSDate not supplied

OCO-2 Satellite-imposed Constraints on Terrestrial Biospheric CO2 Flux over South Asia

The spatiotemporal variability of terrestrial biospheric carbon dioxide (CO2) flux over South Asia is poorly understood. The inverse model “top-down” CO2 flux estimates which rely on atmospheric CO2 observational data are impeded by sparse in situ measurements over this region. The Orbiting Carbon Observatory 2 (OCO-2) provides much needed global retrievals of column-average CO2 dry-air mole fraction (XCO2) with the finest spatial and temporal resolution and highest sensitivity to surface fluxes available to-date. This study conducted global inverse model simulations as part of the second OCO-2 Multi-model Intercomparison Project assimilating version 9 of the OCO-2 XCO2 retrievals in land nadir (LN) and land glint (LG) observing modes and global in situ (IS) measurements. The four-dimensional variational assimilation system with the GEOS-Chem global chemical transport model was used to estimate CO2 fluxes from 2015 to 2018. We then assessed the spatiotemporal variability of optimized CO2 Net ecosystem Exchange (NEE) fluxes aggregated over the South Asia region. The most robust result found by assimilating OCO-2 observations was the constraints imposed on the seasonal cycle of terrestrial biospheric fluxes over South Asia. The seasonality of South Asian NEE estimated by assimilating OCO-2 or global IS data showed a larger seasonal cycle compared to the current understanding of NEE in this region (represented by the prior NEE used in the model). The satellite-in situ joint inversion (IS + LN + LG) led to land flux seasonal amplitude (absolute magnitude of the difference between peak and trough of monthly mean values over a year) of 4.0 PgC yr-1, compared to the prior model (1.5 PgC yr-1). Moreover, OCO-2 data imposed a phase shift in the seasonal cycle, resulting in a large CO2 source in April and a large uptake in September. Most of the features in the seasonal cycle imposed by OCO-2 data are in agreement with previous “top-down” studies assimilating regional aircraft observations. Nevertheless, we recommend a closer examination of the novel features highlighted by the OCO-2 column data over South Asia in future studies by conducting regional inverse modelling with denser in situ and vertically-resolved regional observations along with satellite data.

OCO-2↗
DOE OSTI2022

Probabilistic constrained Bayesian inversion for transpiration cooling

Abstract To enable safe operations in applications such as rocket combustion chambers, the materials require cooling to avoid material damage. Here, transpiration cooling is a promising cooling technique. Numerous studies investigate possibilities to simulate and evaluate the complex cooling mechanism. One naturally arising question is the amount of coolant required to ensure a safe operation. To study this, we introduce an approach that determines the posterior probability distribution of the Reynolds number using an inverse problem and constraining the maximum temperature of the system under parameter uncertainties. Mathematically, this chance inequality constraint is dealt with by a generalized polynomial chaos expansion of the system. The posterior distribution will be evaluated by different Markov chain Monte Carlo based methods. A novel method for the constrained case is proposed and tested among others on two‐dimensional transpiration cooling models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Scalable Space-Time Domain Decomposition Approach for Solving Large Scale Nonlinear Regularized Inverse Ill Posed Problems in 4D Variational Data Assimilation

We address the development of innovative algorithms designed to solve the strong-constraint Four Dimensional Variational Data Assimilation (4DVar DA) problems in large scale applications. We present a space-time decomposition approach which employs the whole domain decomposition, i.e. both along the spacial and temporal direction in the overlapping case, and the partitioning of both the solution and the operator. Starting from the global functional defined on the entire domain, we get to a sort of regularized local functionals on the set of sub domains providing the order reduction of both the predictive and the Data Assimilation models. The algorithm convergence is developed. Performance in terms of reduction of time complexity and algorithmic scalability is discussed on the Shallow Water Equations on the sphere. The number of state variables in the model, the number of observations in an assimilation cycle, as well as numerical parameters as the discretization step in time and in space domain are defined on the basis of discretization grid used by data available at repository Ocean Synthesis/Reanalysis Directory of Hamburg University.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Low-speed single-element airfoil synthesis

The use of recently developed airfoil analysis/design computational tools to clarify, enrich and extend the existing experimental data base on low-speed, single element airfoils is demonstrated. A discussion of the problem of tailoring an airfoil for a specific application at its appropriate Reynolds number is presented. This problem is approached by use of inverse (or synthesis) techniques, wherein a desirable set of boundary layer characteristics, performance objectives, and constraints are specified, which then leads to derivation of a corresponding viscous flow pressure distribution. Examples are presented which demonstrate the synthesis approach, following presentation of some historical information and background data which motivate the basic synthesis process.

Mcmasters, J. H.↗
NASA NTRS1997

Mars Surveyor '98 Landers MVACS Robotic Arm Control System Design Concepts

This paper describes the control system design concepts for the Mars Volatiles and Climate Surveyor (MVACS) Robotic Arm which supports the scientific investigations to be conducted as part of the Mars Surveyor '98 Lander project. Solutions are presented to some of the problems encountered in this demanding space application with its tight constraints on mass, power, volume, and computing resources. Problems addressed include 4-DOF forward and inverse kinematics, trajectory planning to minimize potential impact damage, joint drive train protection, Lander tilt prevention, hardware fault monitoring, and collision avoidance.

Bonitz, Robert G.↗
DOE OSTI2026

Stable and tunable MeV $$\gamma$$-ray generation via dual-laser inverse Thomson scattering from a laser-plasma accelerator

Abstract Inverse Thomson scattering from laser-plasma accelerators offers a pathway to compact, tunable MeV $$\gamma$$ -ray sources for reduced-dose radiography and enhanced performance in nuclear resonance fluorescence (NRF)-based isotope identification. However, photon yield and spectral quality are often limited by constraints on interaction geometry and scatter-laser tunability. Here we demonstrate a MeV $$\gamma$$ -ray source based on a dual-laser inverse Thomson scattering configuration driven by a 100-TW laser-plasma accelerator. Electron beams tunable from 122 to 204 MeV with $$<5$$ mrad divergence and $$<1$$ mrad pointing stability generate $$\gamma$$ rays with peak energies from 276 keV to 1.2 MeV and yields up to $$2\times 10^{7}$$ photons per shot. By independently controlling the interaction position and the scatter-pulse duration, we experimentally match the scatter pulse to the walk-off-limited interaction length. Extending the scatter pulse to 200 fs increases photon production by approximately $$15\%$$ while maintaining operation in the linear Thomson regime, thereby preserving narrow spectral bandwidth and controlled radiation divergence. Radiographic characterization demonstrates MeV-level penetration and $$\approx 0.1$$ mm spatial resolution, while stable operation is sustained over multi-hour timescales across multiple days. These results show that interaction-length optimization provides a scalable strategy for improving photon yield, spectral control, and operational stability in compact laser-plasma-accelerator-driven $$\gamma$$ -ray sources.

Tsai, Hai-En↗
DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗
DOE OSTI2021

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗
DOE OSTI2023

Investigation of the Role of the Role of Nuclear Physics in Heavy Element Nucleosynthesis, through the Study of Key Reactions, and the Improvement of Theoretical Reaction Rates (Final Report)

Neutrino-driven winds in core-collapse supernovae have been identified as a possible site for the production of elements heavier than iron. Traditionally, these neutrino-driven winds have been proposed as the site of the main r-process. Recent simulations fail to reproduce the conditions required for the main r-process. while they remain a promising site for producing the lightest elements beyond iron, e.g., Sr, Y, and Zr through the νp process. The efficiency of the νp process depends on the hydrodynamical conditions, the electron fraction (which is related to the neutrino properties), and the nuclear reactions on many short-lived nuclei with limited (if at all) experimental information. The reaction rates on these nuclei are based on theoretical predictions using the Hauser-Feshbach model. Recent sensitivity studies have highlighted the importance of neutron-induced reactions on these nuclei along the νp process path. This work aimed to experimentally constrain reaction rates that are known to play a key role in the neutrino-p process nucleosynthesis. A secondary subsequently-added objective was to start the implementation of techniques that improve the description of nuclear properties in the Hauser-Feshbach model by extending the microscopic nuclear level density description offered via the shell model to high excitation energies without using experiment-based renormalizations. The main objective of this work was the experimental constraint of the 56 Ni(n,p) 56 Co reaction rate via a measurement of the inverse reaction 56 Co(p,n) 56 Ni at the National Superconducting Cyclotron Laboratory (NSCL) and later the Facility for Rare Isotope Beams (FRIB). This reaction is considered the key one for determining the yields possible by the neutrino-p process. A technique for this type of measurement in inverse kinematics at low energies did not exist before this work. The work also had two secondary objectives. First, to contribute to efforts to measure the same reaction in direct kinematics using a radioactive target at Los Alamos National Laboratory (LANL), and second, to advance work to implement shell-model-deduced microscopic level densities in Hauser-Feshbach calculations. The project has resulted in the development of the first technique to perform (p,n) cross-section measurements in relevant-for-astrophysics low energies in inverse kinematics using a magnetic spectrometer or separator, and neutron detectors for neutron tagging. It has also resulted in the precise measurement of the cross-section of the 40 Ar(p,n) 40 K reaction in a proof-of-principle experiment realized by using a beam-line quadrupole of the ReA3 accelerator of NSCL/FRIB. As part of this project the technique was successfully adapted to make use of the superior acceptance of the Separator for Capture Reactions (SECAR) at FRIB. In this project, the required experimental setup simulations and beam optics were developed and tested with the measurement of the 58 Fe(p,n) 58 Cu reaction cross-section. Additionally, this project contributed with simulation work to the development of a technique to measure (n,p) reactions with radioactive targets at LANL, and which resulted in the measurement of the key 56 Ni(n,p) 56 Co reaction cross-section at neutron energies above ≈1 MeV. Last, the project initiated work in the development of shell model based level densities using the moments method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Gulf stream velocity structure through combined inversion of hydrographic and acoustic Doppler data

Near-surface velocities from an acoustic Doppler instrument are used in conjunction with CTD/O2 data to produce estimates of the absolute flow field off Cape Hatteras. The data set consists of two transects across the Gulf Stream made by the R/V Endeavor cruise EN88 in August 1982. An inverse procedure is applied which makes use of both the acoustic Doppler data and property conservation constraints. Velocity sections at approximately 73 deg. W and 71 deg. W are presented with formal errors of 1-2 cm/s. The net Gulf Stream transports are estimated to be 116 + or - 2 Sv across the south leg and 161 + or - 4 Sv across the north. A Deep Western Boundary Current transport of 4 + or - 1 Sv is also estimated. While these values do not necessarily represent the mean, they are accurate estimates of the synoptic flow field in the region.

Pierce, S. D.↗
NASA NTRS1989

Least squares restoration of multi-channel images

In this paper, a least squares filter for the restoration of multichannel imagery is presented. The restoration filter is based on a linear, space-invariant imaging model and makes use of an iterative matrix inversion algorithm. The restoration utilizes both within-channel (spatial) and cross-channel information as constraints. Experiments using color images (three-channel imagery with red, green, and blue components) were performed to evaluate the filter's performance and to compare it with other monochrome and multichannel filters.

Chin, Roland T.↗
NASA NTRS2001

Estimates of M2 Tidal Energy Dissipation from TOPEX/Poseidon Altimeter Data

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

Egbert, Gary D.↗