Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse problem”

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 523 records · Page 29

Combining Data with Physical Knowledge for Uncertainty Quantification in Certification and Reliability Analysis

Unifying empirical data with predictive models can enable engineering cost-savings through certification by analysis and reliability-based design. Both concepts require rigorous uncertainty quantification (UQ) and robust understanding and treatment of relevant physics. Combining sampling-based UQ algorithms with high-fidelity simulations creates a computational bottleneck that is often alleviated through the use of machine learning (ML). ML can be used to create computationally efficient surrogates for simulations of complex or high-dimensional physical interactions (e.g., multi-phase interactions associated with melt pools in laser powder bed fusion or spatially-dependent material properties in functionally graded materials). However, negative side effects of ML may include a lack of interpretability and negative correlation between event rarity and simulation accuracy due to a lack of training data. As such, it is important to infuse ML algorithms with physics-based guardrails to provide confidence in their predictions. This talk will provide a brief review of recent NASA research at this intersection of physics-based simulation, ML, and UQ with a focus on certification and reliability analysis.

uncertainty quantification↗

Series solution for the planar asymmetric blunt-body problem

Accurate semianalytic solutions to the inverse blunt-body problem have been obtained using a method of series expansion. Rational fractions are employed for series summation and analytic continuation. Angles of incidence up to 30 deg and Mach numbers as low as 2 have been considered. The maximum-entropy streamline will not wet the body surface in asymmetric flow. It may pass either above or below the stagnation streamline. Limit lines appear in the supersonic portion of the flow field, both in the shock layer and in its upstream analytic continuation.

Schwartz, L. W.↗

One-shot omnidirectional pressure integration through matrix inversion

In this work, we present a method to perform 2D and 3D omnidirectional pressure integration from velocity measurements with a single-iteration matrix inversion approach. This work builds upon our previous work, where the rotating parallel ray approach was extended to the limit of infinite rays by taking continuous projection integrals of the ray paths and recasting the problem as an iterative matrix inversion problem. This iterative matrix equation is now 'fast-forwarded' to the 'infinity' iteration, leading to a different matrix equation that can be solved in a single step, thereby presenting the same computational complexity as the Poisson equation. We observe computational speedups of ~10 6 when compared to brute-force omnidirectional integration methods, enabling the treatment of grids of ~10 9 points and potentially even larger in a desktop setup at the time of publication. Further examination of the boundary conditions of our one-shot method shows that omnidirectional pressure integration implements a boundary condition where the boundary points are treated as interior points to the extent that information is available. Finally, we show how the method can be extended from the regular grids typical of particle image velocimetry to the unstructured meshes characteristic of particle tracking velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-Time Bayesian Inference at Extreme Scale: A Digital Twin for Tsunami Early Warning Applied to the Cascadia Subduction Zone

We present a Bayesian inversion-based digital twin that employs acoustic pressure data from seafloor sensors, along with 3D coupled acoustic–gravity wave equations, to infer earthquake-induced spatiotemporal seafloor motion in real time and forecast tsunami propagation toward coastlines for early warning with quantified uncertainties. Our target is the Cascadia subduction zone, with one billion parameters. Computing the posterior mean alone would require 50 years on a 512 GPU machine. Instead, exploiting the shift invariance of the parameter-to-observable map and devising novel parallel algorithms, we induce a fast offline–online decomposition. The offline component requires just one adjoint wave propagation per sensor; using MFEM, we scale this part of the computation to the full El Capitan system (43,520 GPUs) with 92% weak parallel efficiency. Moreover, given real-time data, the online component exactly solves the Bayesian inverse and forecasting problems in 0.2 seconds on a modest GPU system, a ten-billion-fold speedup.

97 MATHEMATICS AND COMPUTING↗

Inverse scattering for an exterior Dirichlet problem

Scattering caused by a metallic cylinder in the field of a wire carrying a periodic current is studied, with a view to determining the location and shape of the cylinder in light of far field measurements between the cylinder and the wire. The associated direct problem is the exterior Dirichlet problem for the Helmholtz equation in two dimensions, and an improved low frequency estimate for its solution by integral equation methods is shown by inverse scattering calculations to be accurate to this estimate. The far field measurements are related to low frequency boundary integral equations whose solutions may be expressed in terms of a mapping function for the exterior of the unknown curve onto the exterior of a unit disk. The conformal transformation's Laurent expansion coefficients can be related to those of the far field, the first of which leads to the calculation of the distance between the source and the cylinder, while the other coefficients are determined by placing the source in a different location.

Hariharan, S. I.↗

Inversion methods in temperature and aerosol remote sounding: Their commonality and differences, and some unexplored approaches

The two remote sensing problems of temperature profiling and aerosol characterization (complex refractive index, size distribution) are considered. These problems differ only in the explicit form of the source function which, for aerosols, includes contributions from both single and multiple scattering processes. When the observables are the spectral extinction or the single scattering of the source radiation, the associated problem is completely analogous to the linearized temperature inversion problem. Methods for obtaining the solution of the linear problem are classified following three main categories: (1) derivation of properties that all solutions satisfy, which must then be properties of the actual solution; (2) regularization of the ill-posed problem; and (3) data changes within their domain of uncertainty in order to avoid the basic instability. A number of unexplored methods are indicated.

Fymat, A. L.↗

Inversion for the driving forces of plate tectonics

Inverse modeling techniques have been applied to the problem of determining the roles of various forces that may drive and resist plate tectonic motions. Separate linear inverse problems have been solved to find the best fitting pole of rotation for finite element grid point velocities and to find the best combination of force models to fit the observed relative plate velocities for the earth's twelve major plates using the generalized inverse operator. Variance-covariance data on plate motion have also been included. Results emphasize the relative importance of ridge push forces in the driving mechanism. Convergent margin forces are smaller by at least a factor of two, and perhaps by as much as a factor of twenty. Slab pull, apparently, is poorly transmitted to the surface plate as a driving force. Drag forces at the base of the plate are smaller than ridge push forces, although the sign of the force remains in question.

Richardson, R. M.↗

Minimal inversion, command matching and disturbance decoupling in multivariable systems

The present treatment of the related problems of minimal inversion and perfect output control in linear multivariable systems uses a simple analytical expression for the inverse of a square multivariate system's transfer-function matrix to construct a minimal-order inverse of the system. Because the poles of the minimal-order inverse are the transmission zeros of the system, necessary and sufficient conditions for the inverse system's stability are simply stated in terms of the zero polynomial of the original system. A necessary and sufficient condition for the existence of the required controllers is that the plant zero polynomial be neither identical to zero nor unstable.

Seraji, H.↗

Inversion of vegetation canopy reflectance models for estimating agronomic variables. I - Problem definition and initial results using the Suits model

An important but relatively uninvestigated problem in remote sensing is the inversion of vegetative canopy reflectance models to obtain agrophysical parameters, given measured reflectances. The problem is here formally defined and its solution outlined. Numerical nonlinear optimization techniques are used to implement this inversion to obtain the leaf area index using Suits' canopy reflectance model. The results for a variety of cases indicate that this can be done successfully using infrared reflectances at different views or azimuth angles or a combination thereof. The other parameters of the model must be known, although reasonable measurement errors can be tolerated without seriously degrading the accuracy of the inversion. The application of the technique to ground based remote-sensing experiments is potentially useful, but is limited to the degree to which the canopy reflectance model can accurately predict observed reflectances.

Goel, N.↗

Inversion of multiwavelength radiometer measurements by three-dimensional filtering

Remote sensing data from satellites typically have three dimensions: scan position, spacecraft position, and wavelength. Inversion of the radiometric data to infer geophysical parameters is a filtering problem in which the dimension of wavelength (or channel number) is transformed into a dimension of geophysical parameters, and the most general solution is a three-dimensional filter. Linear filters have the advantages of computational speed and easily described transfer functions; but often the measurements are nonlinear functions of the parameters to be inferred. To the extent that the nonlinear inversion problem is overdetermined, it can be modeled by a critically determined linear problem. As an example, inversion of Scanning Multichannel Microwave Radiometer (SMMR) data by means of a three-dimensional Wiener Filter is described. Atmospheric water vapor content, rain liquid water content, surface wind speed and surface temperature are the parameters inferred from the measurements. Nonprecipitating liquid water and water vapor scale height are also modeled but not retrieved. The a priori statistics on which the filter is trained have the effect of governing the selection of a trade-off point of noise as a function of resolution (in all three retrieval dimensions).

Rosenkranz, P. W.↗

Estimation of extreme temperatures in direct solar methane pyrolysis within a porous medium

Porous media have wide application in renewable energy conversion processes, such as solar-thermal fuels production and decarbonization. Heat transport mechanisms within porous media can be highly complex, particularly under extreme conditions encountered in concentrated solar thermal reactors in which direct measurement of temperature is challenging. Here, we implement and report an inverse heat conduction model to estimate the temperature distribution throughout a porous substrate domain in a direct solar methane pyrolysis process. By solving a two-dimensional heat transfer problem and applying an inverse optimization algorithm, we estimate the quasi-steady state spatial temperature distribution in a fibrous porous carbon substrate. The results are validated indirectly by experimentally measured graphite deposition and a simplified reaction kinetic model.

finite difference method↗

Low Velocity Airdrop Tests of an X-38 Backup Parachute Design

The NASA Johnson Space Center's X-38 program designed a new backup parachute system to recover the 25,000 lb X-38 prototype for the Crew Return Vehicle spacecraft. Due to weight and cost constraints, the main backup parachute design incorporated rapid and low cost fabrication techniques using off-the-shelf materials. Near the vent, the canopy was constructed of continuous ribbons, to provide more damage tolerance. The remainder of the canopy was a constructed with a continuous ringslot design. After cancellation of the X-38 program, the parachute design was resized, built, and drop tested for Natick Soldiers Center's Low Velocity Air Drop (LVAD) program to deliver cargo loads up to 22,000 lbs from altitudes as low as 500 feet above the ground. Drop tests results showed that the 500-foot LVAD parachute deployment conditions cause severe skirt inversion and inflation problems for large parachutes. The bag strip occurred at a high angle of attack, causing skirt inversion before the parachute could inflate. The addition of a short reefing line prevented the skirt inversion. Using a lower porosity in the vent area, than is normally used in large parachutes, improved inflation. The drop testing demonstrated that the parachute design could be refined to meet the requirements for the 500-foot LVAD mission.

Stein, Jenny M.↗

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗

On dynamics and control of multi-link flexible space manipulators

In this paper dynamics, inverse dynamics, and control problems for multi-link flexible space manipulators are presented. In deriving the flexible manipulator dynamics the following are assumed: flexible deformations are relatively small; angular rates of the links are much smaller than their fundamental frequencies; nonlinear terms (centrifugal and Coriolis forces) in the flexible manipulator model are the same as those in the rigid body model. These assumptions are reasonable for large space manipulators, such as the space crane. Flexible displacements are measured with respect to the rigid body configuration, for which a linear time-varying system is obtained. The inverse dynamics problem consists of determination of joint torques, given tip trajectory, such that joint angles in flexible configuration are equal to the angles in the rigid body configuration. The manipulator control system consists of the feedforward compensation and feedback control loops. Simulation results of a two-link space crane with large payload show that the performance of this linearized dynamics and control approach is reasonable and robust subject to parameter variations during slew operations.

Gawronski, W.↗

Inverse design for waveguide dispersion with a differentiable mode solver

Inverse design of optical components based on adjoint sensitivity analysis has the potential to address the most challenging photonic engineering problems. However, existing inverse design tools based on finite-difference-time-domain (FDTD) models are poorly suited for optimizing waveguide modes for adiabatic transformation or perturbative coupling, which lies at the heart of many important photonic devices. Among these, dispersion engineering of optical waveguides is especially challenging in ultrafast and nonlinear optical applications involving broad optical bandwidths and frequency-dependent anisotropic dielectric material response. In this work, we develop gradient back-propagation through a general-purpose electromagnetic eigenmode solver and use it to demonstrate waveguide dispersion optimization for second harmonic generation with maximized phase-matching bandwidth. This optimization of three design parameters converges in eight steps, reducing the computational cost of optimization by ∼100x compared to exhaustive search and identifying new designs for broadband optical frequency doubling of laser sources in the 1.3–1.4 µm wavelength range. Furthermore, we demonstrate that the computational cost of gradient back-propagation is independent of the number of parameters, as required for optimization of complex geometries. This technique enables practical inverse design for a broad range of previously intractable photonic devices.

Gray, Dodd (ORCID:000000030469599X)↗

The inversion of VLF-LF sounder data to obtain D-region electron-density distributions

The direct problem is to compute reflection coefficients for an assumed electron-density profile, using a full-wave solution. The inverse problem is treated numerically, using derivatives of the reflection matrix with respect to model parameters. A technique for the rapid computation of the derivatives is described. An algorithm for determining resolution in the profile is discussed.

Shellman, C. H.↗

Inverse kinematic-based robot control

A fundamental problem which must be resolved in virtually all non-trivial robotic operations is the well-known inverse kinematic question. More specifically, most of the tasks which robots are called upon to perform are specified in Cartesian (x,y,z) space, such as simple tracking along one or more straight line paths or following a specified surfacer with compliant force sensors and/or visual feedback. In all cases, control is actually implemented through coordinated motion of the various links which comprise the manipulator; i.e., in link space. As a consequence, the control computer of every sophisticated anthropomorphic robot must contain provisions for solving the inverse kinematic problem which, in the case of simple, non-redundant position control, involves the determination of the first three link angles, theta sub 1, theta sub 2, and theta sub 3, which produce a desired wrist origin position P sub xw, P sub yw, and P sub zw at the end of link 3 relative to some fixed base frame. Researchers outline a new inverse kinematic solution and demonstrate its potential via some recent computer simulations. They also compare it to current inverse kinematic methods and outline some of the remaining problems which will be addressed in order to render it fully operational. Also discussed are a number of practical consequences of this technique beyond its obvious use in solving the inverse kinematic question.

Wolovich, W. A.↗

Analysis of space telescope data collection system

An analysis of the expected performance for the Multiple Access (MA) system is provided. The analysis covers the expected bit error rate performance, the effects of synchronization loss, the problem of self-interference, and the problem of phase ambiguity. The problem of false acceptance of a command word due to data inversion is discussed. A mathematical determination of the probability of accepting an erroneous command word due to a data inversion is presented. The problem is examined for three cases: (1) a data inversion only, (2) a data inversion and a random error within the same command word, and a block (up to 256 48-bit words) containing both a data inversion and a random error.

Ingels, F. M.↗