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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1989

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.↗

NASA NTRS1988

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.↗

NASA NTRS1990

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↗

NASA NTRS2002

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↗

NASA NTRS2007

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↗

DOE OSTI2025

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)↗

NASA NTRS1976

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.↗

DOE OSTI2025

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↗

NASA NTRS1987

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.↗

DOE OSTI2026

Dual Representations and H ∞ -Optimal Control of Partial Differential Equations

We consider H ∞ -optimal state-feedback control of the class of linear Partial Differential Equations (PDEs) which admit a Partial Integral Equation (PIE) representation. While linear matrix inequalities are commonly used for optimal control of Ordinary Differential Equations (ODEs), the absence of a universal state-space representation and suitable dual form prevents such methods from being applied to optimal control of PDEs. Specifically, for ODEs, the controller synthesis problem is defined in state-space, and duality is used to resolve the bilinearity of that synthesis problem. Recently, the PIE representation was proposed as a universal state-space representation for linear PDE systems. In this paper, we show that any PDE system represented by a PIE admits a dual PIE with identical stability and I/O properties. This result allows us to reformulate the stabilizing and optimal state-feedback control problems as convex optimization over the cone of positive Partial Integral (PI) operators. Operator inversion formulae then allow us to construct feedback gains for the original PDE system. The results are verified through application to several canonical problems in optimal control of PDEs and indicate the resulting bounds on H ∞ norm are not conservative.

42 ENGINEERING↗

NASA NTRS1974

A study to determine the usefulness of interval analysis in solving problems in celestial mechanics

This investigation was undertaken to determine the usefulness of interval analysis to numerical integration and matrix inversion techniques and to combine these results to determine the value of interval analysis in bounding computational errors in the two-body problem. Conclusions were that interval analysis may be worthwhile in certain small scale isolated problems, but its usefulness in any large scale problem is doubtful.

Walling, D.↗

NASA NTRS1985

Adaptive Inverse Control for Rotorcraft Vibration Reduction

This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.

Jacklin, Stephen A.↗

DOE OSTI2025

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗

NASA NTRS1976

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗

NASA NTRS1976

Free vibration analysis of spinning flexible space structures

Gupta, K. K.↗

NASA NTRS1979

Ion beam probing of electrostatic fields

The determination of a cylindrically symmetric, time-independent electrostatic potential V in a magnetic field B with the same symmetry by measurements of the deflection of a primary beam of ions is analyzed and substantiated by examples. Special attention is given to the requirements on canonical angular momentum and total energy set by an arbitrary, nonmonotone V, to scaling laws obtained by normalization, and to the analogy with ionospheric sounding. The inversion procedure with the Abel analysis of an equivalent problem with a one-dimensional fictitious potential is used in a numerical experiment with application to the NASA Lewis Modified Penning Discharge. The determination of V from a study of secondary beams of ions with increased charge produced by hot plasma electrons is also analyzed, both from a general point of view and with application to the NASA Lewis SUMMA experiment. Simple formulas and geometrical constructions are given for the minimum energy necessary to reach the axis, the whole plasma, and any point in the magnetic field. The common, simplifying assumption that V is a small perturbation is critically and constructively analyzed; an iteration scheme for successively correcting the orbits and points of ionization for the electrostatic potential is suggested.

Persson, H.↗

NASA NTRS1987

The magnetic field of Uranus

Aspherical harmonic model of the planetary magnetic field of Uranus is obtained from the Voyager 2 encounter observations using generalized inverse techniques which allow partial solutions to complex (underdetermined) problems. The Goddard Space Flight Center 'Q3' model is characterized by a large dipole tilt (58.6 deg) relative to the rotation axis, a dipole moment of 0.228 G R(Uranus radii cubed) and an unusually large quadrupole moment. Characteristics of this complex model magnetic field are illustrated using contour maps of the field on the planet's surface and discussed in the context of possible dynamo generation in the relatively poorly conducting 'ice' mantle.

Connerney, J. E. P.↗

NASA NTRS1996

Quasi-Optimal Schwarz Methods for the Conforming Spectral Element Discretization

Fast methods are proposed for solving the system K(sub N)x = b resulting from the discretization of self-adjoint elliptic equations in three dimensional domains by the spectral element method. The domain is decomposed into hexahedral elements, and in each of these elements the discretization space is formed by polynomials of degree N in each variable. Gauss-Lobatto-Legendre (GLL) quadrature rules replace the integrals in the Galerkin formulation. This system is solved by the preconditioned conjugate gradients method. The conforming finite element space on the GLL mesh consisting of piecewise Q(sub 1) elements produces a stiffness matrix K(sub h) that is spectrally equivalent to the spectral element stiffness matrix K(sub N). The action of the inverse of K(sub h) is expensive for large problems, and is therefore replaced by a Schwarz preconditioner B(sub h) of this finite element stiffness matrix. The preconditioned operator then becomes B(sub h)(exp -l)K(sub N). The technical difficulties stem from the nonregularity of the mesh. Tools to estimate the convergence of a large class of new iterative substructuring and overlapping Schwarz preconditioners are developed. This technique also provides a new analysis for an iterative substructuring method proposed by Pavarino and Widlund for the spectral element discretization.

Casarin, Mario↗