Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse problems”

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 649 records · Page 36

Weighted current sheets supported in normal and inverse configurations - A model for prominence observations

A technique which incorporates both photospheric and prominence magnetic field observations is used to analyze the magnetic support of solar prominences in two dimensions. The prominence is modeled by a mass-loaded current sheet which is supported against gravity by magnetic fields from a bipolar source in the photosphere and a massless line current in the corona. It is found that prominence support can be achieved in three different kinds of configurations: an arcade topology with a normal polarity; a helical topology with a normal polarity; and a helical topology with an inverse polarity. In all cases the important parameter is the variation of the horizontal component of the prominence field with height. Adding a line current external to the prominence eliminates the nonsupport problem which plagues virtually all previous prominence models with inverse polarity.

Demoulin, P.↗

The analysis of a rocket tomography measurement of the N2+3914A emission and N2 ionization rates in an auroral arc

Techniques were developed for recovering two-dimensional distributions of auroral volume emission rates from rocket photometer measurements made in a tomographic spin scan mode. These tomographic inversion procedures are based upon an algebraic reconstruction technique (ART) and utilize two different iterative relaxation techniques for solving the problems associated with noise in the observational data. One of the inversion algorithms is based upon a least squares method and the other on a maximum probability approach. The performance of the inversion algorithms, and the limitations of the rocket tomography technique, were critically assessed using various factors such as (1) statistical and non-statistical noise in the observational data, (2) rocket penetration of the auroral form, (3) background sources of emission, (4) smearing due to the photometer field of view, and (5) temporal variations in the auroral form. These tests show that the inversion procedures may be successfully applied to rocket observations made in medium intensity aurora with standard rocket photometer instruments. The inversion procedures have been used to recover two-dimensional distributions of auroral emission rates and ionization rates from an existing set of N2+3914A rocket photometer measurements which were made in a tomographic spin scan mode during the ARIES auroral campaign. The two-dimensional distributions of the 3914A volume emission rates recoverd from the inversion of the rocket data compare very well with the distributions that were inferred from ground-based measurements using triangulation-tomography techniques and the N2 ionization rates derived from the rocket tomography results are in very good agreement with the in situ particle measurements that were made during the flight. Three pre-prints describing the tomographic inversion techniques and the tomographic analysis of the ARIES rocket data are included as appendices.

Mcdade, Ian C.↗

Airfoil design method using the Navier-Stokes equations

An airfoil design procedure is described that was incorporated into an existing 2-D Navier-Stokes airfoil analysis method. The resulting design method, an iterative procedure based on a residual-correction algorithm, permits the automated design of airfoil sections with prescribed surface pressure distributions. The inverse design method and the technique used to specify target pressure distributions are described. It presents several example problems to demonstrate application of the design procedure. It shows that this inverse design method develops useful airfoil configurations with a reasonable expenditure of computer resources.

Malone, J. B.↗

Experimental Flight Validation of the Prandtl 1933 Bell Spanload

This report describes the validation of the 1933 Prandtl bell spanload. This spanload is the minimum induced drag of a wing for a given structural weight with properties that eliminate adverse yaw. Aircraft using the Prandtl bell spanload were flown and investigated. The results of this research show that many previously held assumptions should be rethought, and the creation of aircraft using the Prandtl bell spanload will require considerable new techniques. Part of this work centered on the use of inverse methods. The usual first-step approach to a computational fluids problem is to create a geometry of the aircraft. Once this geometry exists, the computational fluids solution has been solved; however, the problem of creating the geometry still exists. A very small segment of the computational fluids world has concentrated on inverse solutions. This design approach begins with an end result - the computational fluids solution; from this end result, the geometry is sought. A more generic inverse tool was desired that would allow for the design of wings - specifically Prandtl bell spanload wings. Such a tool has been developed and allows for variable taper, aspect ratio, sweep, airfoils, and design-lift coefficients; output from the tool results in twist distribution of wings.

Albion H. Bowers↗

Inversion of solar extinction data from the Apollo-Soyuz Test Project Stratospheric Aerosol Measurement (ASTP/SAM) experiment

The inversion methods are reported that have been used to determine the vertical profile of the extinction coefficient due to the stratospheric aerosols from data measured during the ASTP/SAM solar occultation experiment. Inversion methods include the onion skin peel technique and methods of solving the Fredholm equation for the problem subject to smoothing constraints. The latter of these approaches involves a double inversion scheme. Comparisons are made between the inverted results from the SAM experiment and near simultaneous measurements made by lidar and balloon born dustsonde. The results are used to demonstrate the assumptions required to perform the inversions for aerosols.

Pepin, T. J.↗

Estimation of the two-dimensional ocean current shear field with a synthetic aperture radar

A double inverse method to convert SAR data to estimates of current shear was developed theoretically. Problems of applying the technique to real SAR data were examined. The first step of the inversion process involved elimination of noise in the (x-vector, k-vector) wave spectra deduced from SAR images. The second step, estimation of current shear, required the existence of well defined swell and a favorable single of interaction between swell and ocean current. Application of the method in the Antarctic circumpolar current encountered sea conditions that both satisfied and violated these requirements. Where the wave-current geometry was favorable, realistic surface current velocity fields were obtained, suggesting the presence of a substantial mesoscale eddy. Another attempt to apply the method in the Gulf Stream produced wave spectral estimates that showed no clear peak and hence produced no current estimates.

Barnett, Tim P.↗

Robust, nonlinear, high angle-of-attack control design for a supermaneuverable vehicle

High angle-of-attack flight control laws are developed for a supermaneuverable fighter aircraft. The methods of dynamic inversion and structured singular value synthesis are combined into an approach which addresses both the nonlinearity and robustness problems of flight at extreme operating conditions. The primary purpose of the dynamic inversion control elements is to linearize the vehicle response across the flight envelope. Structured singular value synthesis is used to design a dynamic controller which provides robust tracking to pilot commands. The resulting control system achieves desired flying qualities and guarantees a large margin of robustness to uncertainties for high angle-of-attack flight conditions. The results of linear simulation and structured singular value stability analysis are presented to demonstrate satisfaction of the design criteria. High fidelity nonlinear simulation results show that the combined dynamics inversion/structured singular value synthesis control law achieves a high level of performance in a realistic environment.

Adams, Richard J.↗

Fusion of Visible and Thermal-Infrared Imagery for SLAM for Landing on Icy Moons

This paper addresses the problem of localization for landing on the surface of icy moons, like Europa or Enceladus. Due to the possibility of specular reflection as well as high bulk albedo, icy surfaces present new challenges that make traditional vision-based navigation systems relying on visible imagery unreliable. We propose augmenting visible light cameras with a thermal-infrared camera using inverse-depth parameterized monocular EKF-SLAM to address problems arising from the appearance of icy moons. Results were obtained from a novel procedural Europa surface simulation which models the appearance and the thermal properties simultaneously from physically-based methods. In this framework, we show that thermal features improve localization by 23% on average when compared to a visible camera. Moreover, fusing both sensing modalities increases the improvement in localization to 31% on average, compared to using a visible light camera alone.

thermal vision↗

Development of an unsteady wake theory appropriate for aeroelastic analyses of rotors in hover and forward flight

The purpose of this research is the development of an unsteady aerodynamic model for rotors such that it can be used in conventional aeroelastic analysis (e.g., eigenvalue determination and control system design). For this to happen, the model must be in a state-space formulation such that the states of the flow can be defined, calculated and identified as part of the analysis. The fluid mechanics of the problem is given by a closed-form inversion of an acceleration potential. The result is a set of first-order differential equations in time for the unknown flow coefficients. These equations are hierarchical in the sense that they may be truncated at any number of radial or azimuthal terms.

Peters, David A.↗

Computing Optimal Multiarc Trajectories

Optimal Multi-Arc Trajectories (OMAT) computer program designed to calculate solution to optimal-trajectory problem in cases of low thrust-to-weight ratios. Developed for two-body exoatmospheric problem with three degrees of freedom in inverse-square force field. Offers two different options: impulsive changes in velocity and finite burns with low thrust-to-weight ratios. Two distinct solutions available from program: optimal multi-impulse (OMI) solution and optimal multiburn (OMB) solution. Two solutions obtained separately, or results of OMI solution used to guess unknown parameters of OMB solution. Written for DEC VAX-series computer. Written completely in FORTRAN 77.

Jezewski, Donald J.↗

Recovery of excitation intensity dependence in pulsed, focused beams - Nonsaturated case

It is shown how the dependence of a nonsaturated excitation process on local, instantaneous intensity can be obtained from temporally and spatially unresolved data by inversion, even for a nonuniform intensity distribution in the probe volume. This treatment is in contrast with the usual approach in which effects of nonuniform excitation are either disregarded or simulated. For exponential intensity profiles, the solution is obtained in the form of an Abel inversion. The ill-conditioned nature of the problem is demonstrated, and extension to the study of intensity-dependent line shapes is made.

Ruyten, Wilhelmus M.↗

Doppler radar imaging of spherical planetary surfaces

This paper describes a technique for using echo power-spectra for making radar images of spherical planetary targets. In developing the Doppler-radar imaging system, the target's reflectivity distribution is expanded in a truncated spherical harmonic series; the distribution of echo power in rotational phase and the Doppler frequency are obtained in form of a system of linear equations which express the echo spectra as functions of the series coefficients. To estimate the coefficients from an observed phase-Doppler distribution, the inversion is cast as a least-squares problem, and is solved using singular value decompositions, yielding a linear imaging system which can be fully characterized by its impulse response as a function of latitude. Simulations were designed and used to explore the capabilities of Doppler radar imaging, and the sensitivity of imaging to such factors as subradar altitude coverage and signal-to-noise ratio was investigated.

Hudson, R. Scott↗

Rayleigh-Taylor Gravity Waves and Quasiperiodic Oscillation Phenomenon in X-ray Binaries

Accretion onto compact objects in X-ray binaries (black hole, neutron star (NS), white dwarf) is characterized by non-uniform flow density profiles. Such an effect of heterogeneity in presence of gravitational forces and pressure gradients exhibits Rayleigh-Taylor gravity waves (RTGW). They should be seen as quasiperiodic wave oscillations (QPO) of the accretion flow in the transition (boundary) layer between the Keplerian disk and the central object. In this paper the author shows that the main QPO frequency, which is very close to the Keplerian frequency, is split into separate frequencies (hybrid and low branch) under the influence of the gravitational forces in the rotational frame of reference. The RTGWs must be present and the related QPOs should be detected in any system where the gravity, buoyancy and Coriolis force effects cannot be excluded (even in the Earth and solar environments). The observed low and high QPO frequencies are an intrinsic signature of the RTGW. The author elaborates the conditions for the density profile when the RTGW oscillations are stable. A comparison of the inferred QPO frequencies with QPO observations is presented. The author finds that hectohertz frequencies detected from NS binaries can be identified as the RTGW low branch frequencies. The author also predicts that an observer can see the double NS spin frequency during the NS long (super) burst events when the pressure gradients and buoyant forces are suppressed. The Coriolis force is the only force which acts in the rotational frame of reference and its presence causes perfect coherent pulsations with a frequency twice of the NS spin. The QPO observations of neutron binaries have established that the high QPO frequencies do not go beyond of the certain upper limit. The author explains this observational effect as a result of the density profile inversions. Also the author demonstrates that a particular problem of the gravity waves in the rotational frame of reference in the approximation of very small pressure gradients is reduced to the problem of the classical oscillator in the rotational frame of reference which was previously introduced and applied for the interpretation of kHZ QPO observation by Osherovich & Titarchuk.

Titarchuk, Lev↗

Nanorod-Based Fast-Response Pressure-Sensitive Paints

A proposed program of research and development would be devoted to exploitation of nanomaterials in pressuresensitive paints (PSPs), which are used on wind-tunnel models for mapping surface pressures associated with flow fields. Heretofore, some success has been achieved in measuring steady-state pressures by use of PSPs, but success in measuring temporally varying pressures has been elusive because of the inherent slowness of the optical responses of these materials. A PSP contains a dye that luminesces in a suitable wavelength range in response to photoexcitation in a shorter wavelength range. The luminescence is quenched by oxygen at a rate proportional to the partial pressure of oxygen and thus proportional to the pressure of air. As a result, the intensity of luminescence varies inversely with the pressure of air. The major problem in developing a PSP that could be easily applied to a wind-tunnel model and could be useful for measuring rapidly varying pressure is to provide very high gas diffusivity for rapid, easy transport of oxygen to and from active dye molecules. Most PSPs include polymer-base binders, which limit the penetration of oxygen to dye molecules, thereby reducing responses to pressure fluctuations. The proposed incorporation of nanomaterials (somewhat more specifically, nanorods) would result in paints having nanostructured surfaces that, relative to conventional PSP surfaces, would afford easier and more nearly complete access of oxygen molecules to dye molecules. One measure of greater access is effective surface area: For a typical PSP as proposed applied to a given solid surface, the nanometer-scale structural features would result in an exposed surface area more than 100 times that of a conventional PSP, and the mass of proposed PSP needed to cover the surface would be less than tenth of the mass of the conventional PSP. One aspect of the proposed development would be to synthesize nanorods of Si/SiO2, in both tangle-mat and regular- array forms, by use of chemical vapor deposition (CVD) and wet chemical processes, respectively. The rods would be coated with a PSP dye, and the resulting PSP signals would be compared with those obtained from PSP dye coats on conventional support materials. Another aspect of the proposed development would be to seek to exploit the quantum properties of nanorods of a suitable semiconductor (possibly GaN), which would be synthesized by CVD. These quantum properties of semiconductor nanorods include narrow-wavelength-band optical absorption and emission characteristics that vary with temperature. The temperature sensitivity might enable simultaneous measurement of fluctuating temperature and pressure and to provide a temperature correction for the PSP response.

Bencic, Timothy↗

Sign Problem in Tensor-Network Contraction

We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025

Chen, Jielun (ORCID:0000000178411545)↗

First order impulsive solutions

A mathematically rigorous derivation is given of first order corrections to multi-impulse approximations to the solutions to space flight optimization problems with bang-bang control. The rocket was subjected to an inverse square gravitational force and to a thrust force with constant magnitude. The mass decreased linearly with time. An optimal impulsive solution was obtained for a problem with given initial and final conditions. The method was then used to obtain first-order corrections to the initial values of the costate variables. Indications are given on how the theory may be extended to higher order corrections. The theory was applied to intercept and rendezvous problems.

Andrus, J. F.↗

ForceFinder

SAND2025-11750O ForceFinder extends the Structural Dynamics Python Libraries (SDynPy) with comprehensive tools for inverse source estimation (ISE) tasks via frequency response function (FRF) matrix inversion. The software is designed for transfer path analysis and multiple-input/multiple-output (MIMO) vibration control problems. It allows users to estimate sources through various algorithms, from the basic Moore-Penrose pseudo-inverse to statistical learning methods such as Tikhonov regularization via an L-curve and elastic net regularization via an information criterion. ForceFinder uses an object-oriented framework, where all components of the ISE problem—such as FRFs, responses, and transformations—are stored in a "SourcePathReceiver" object. This software can be applied to any noise and vibration problem. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Carter, Steven [Sandia National Lab. (SNL-CA), Liv↗

A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations

A direct, implicit, numerical solution algorithm, with second-order accuracy in space and time, is constructed for the three-dimensional unsteady incompressible Navier-Stokes equations formulated in terms of velocity and vorticity, using generalized orthogonal coordinates to achieve the accurate solution of complex viscous flow configurations. A numerically stable, efficient, direct inversion procedure is developed for the computationally intensive divergence-curl elliptic velocity problem. This overdetermined partial differential operator is first formulated as a uniquely determined, nonsingular matrix-vector problem; this aspect of the procedure is a unique feature of the present analysis. The three-dimensional vorticity-transport equation is solved by a modified factorization technique which completely eliminates the need for any block-matrix inversions and only scalar tridiagonal matrices need to be inverted. The method is applied to the test problem of the three-dimensional flow within a shear-driven cubical box. Coherent streamwise vortex structures are observed within the steady-state flow at Re = 100.

Osswald, G. A.↗