Results for “inverse problems” — Page 55

NASA NTRS1989

An analytic development of orbit determination for a distant, planetary orbiter

With the advent of the Mariner '71 Mission, NASA has been sending spacecraft to orbit various distant bodies within the solar system. At present, there is still no adequate theory describing the inherent state estimation accuracy, based on two-way, coherent range-rate data. It is the purpose of this article to lay the groundwork for a general elliptic theory, and in addition to provide an analytic solution for the special case of circular orbits. It is shown that circular orbits about distant planets may suffer singularities in over-all position error estimation. These singularities are due to orbit inclination, placement of the line-of-nodes, and insignificant cross-velocity at the start and end of retrograde motion when orbiting a superior planet. Even though these conclusions appear to yield poor state estimation, one should not be unduly alarmed inasmuch as the stated conditions for singularity are not maintained for extended periods during typical mission scenarios. However, mission analysts should be aware of these potential pitfalls and realize that spuriously large results for circular orbiters can be obtained and are not the result of incorrect assumptions or faulty software. The general elliptic problem appears so involved that analytic inversion at this time is just not feasible, and in any case the resulting expression for the position error would likely be so lengthy that any understanding would be lost in the maze.

Russell, R. K.↗

NASA NTRS1995

Reverberation mapping by regularized linear inversion

Reverberation mapping of active galactic nucleus (AGN) emission-line regions requires the numerical deconvolution of two time series. We suggest the application of a new method, regularized linear inversion, to the solution of this problem. This method possesses many good features; it imposes no restrictions on the sign of the response function; it can provide clearly defined uncertainty estimates; it involves no guesswork about unmeasured data; it can give a clear indication of when the underlying convolution model is inadequate; and it is computationally very efficient. Using simulated data, we find the minimum S/N and length of the time series in order for this method to work satisfactorily. We also define guidelines for choosing the principal tunable parameter of the method and for interpreting the results. Finally, we reanalyze published data from the 1989 NGC 5548 campaign using this new method and compare the results to those previously obtained by maximum entropy analysis. For some lines we find good agreement, but for others, especially C III lambda(1909) and Si IV lambda(1400), we find significant differences. These can be attributed to the inability of the maximum entropy method to find negative values of the response function, but also illustrate the nonuniqueness of any deconvolution technique. We also find evidence that certain line light curves (e.g., C IV lambda(1549)) cannot be fully described by the simple linear convolution model.

Krolik, Julian H.↗

NASA NTRSDate not supplied

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

A mixed-fidelity low-boom multidisciplinary optimization (MDO) problem is formulated for integration of low-fidelity MDO and computational fluid dynamics (CFD) based low-boom inverse design optimization. The mixed-fidelity low-boom MDO problem aims to enforce the weight consistency for CFD-based low-boom inverse design optimization: the optimum low-boom configuration is designed for the weight at the start of cruise of the same configuration for the low-boom overland mission. Moreover, it also seeks the optimum trades among the maximum takeoff gross weight (MTOGW), cruise Mach, and range for the low-boom overland mission while meeting the requirement of flying a transatlantic overwater mission. A block coordinate optimization (BCO) method is developed to find an approximate solution of the mixed-fidelity low-boom MDO problem. The BCO method is successfully applied to generate two CFD-based low-boom configurations that closely match two reversed equivalent area targets with ground noise levels below 70 PLdB, respectively. Moreover, these CFD-based low-boom configurations can be obtained with minor wing modifications from the solutions of the low-fidelity MDO problem. The low-fidelity MDO solutions are Pareto points for constrained multiobjective optimization of MTOGW, low-boom cruise Mach, and low-boom range. The best low-boom concept has a predetermined fuselage shape tailored for passengers and main gear storage, carries 40 passengers at seat pitch of 48 in, flies a low-boom overland mission with cruise Mach of 1.8 and range of 2,950 nm, cruises overwater at Mach 1.8 with range of 3,600 nm, satisfies the specified constraints for landing/cruise/takeoff, has MTOGW of 145,164 lb, trims the low-boom cruise flight with fuel redistributions instead of control surface deflections, and has a reversed equivalent area distribution closely matching a target with ground noise level below 70 PLdB.

multidisciplinary optimization↗

NASA NTRS2021

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

A mixed-fidelity low-boom multidisciplinary optimization (MDO) problem is formulated for integration of low-fidelity MDO and computational fluid dynamics (CFD) based low-boom inverse design optimization. The mixed-fidelity low-boom MDO problem aims to enforce the weight consistency for CFD-based low-boom inverse design optimization: the optimum low-boom configuration is designed for the weight at the start of cruise of the same configuration for the low-boom overland mission. Moreover, it also seeks the optimum trades among the maximum takeoff gross weight (MTOGW), cruise Mach, and range for the low-boom overland mission while meeting the requirement of flying a transatlantic overwater mission. A block coordinate optimization (BCO) method is developed to find an approximate solution of the mixed-fidelity low-boom MDO problem. The BCO method is successfully applied to generate two CFD-based low-boom configurations that closely match two reversed equivalent area targets with ground noise levels below 70 PLdB, respectively. Moreover, these CFD-based low-boom configurations can be obtained with minor wing modifications from the solutions of the low-fidelity MDO problem. The low-fidelity MDO solutions are Pareto points for constrained multiobjective optimization of MTOGW, low-boom cruise Mach, and low-boom range. The best low-boom concept has a predetermined fuselage shape tailored for passengers and main gear storage, carries 40 passengers at seat pitch of 48 in, flies a low-boom overland mission with cruise Mach of 1.8 and range of 2,950 nm, cruises overwater at Mach 1.8 with range of 3,600 nm, satisfies the specified constraints for landing/cruise/takeoff, has MTOGW of 145,164 lb, trims the low-boom cruise flight with fuel redistributions instead of control surface deflections, and has a reversed equivalent area distribution closely matching a target with ground noise level below 70 PLdB.

low-boom supersonic transports↗

NASA NTRSDate not supplied

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

A mixed-fidelity low-boom multidisciplinary optimization (MDO) problem is formulated for integration of low-fidelity MDO and computational fluid dynamics (CFD) based low-boom inverse design optimization. The mixed-fidelity low-boom MDO problem aims to enforce the weight consistency for CFD-based low-boom inverse design optimization: the optimum low-boom configuration is designed for the weight at the start of cruise of the same configuration for the low-boom overland mission. Moreover, it also seeks the optimum trades among the maximum takeoff gross weight (MTOGW), cruise Mach, and range for the low-boom overland mission while meeting the requirement of flying a transatlantic overwater mission. A block coordinate optimization (BCO) method is developed to find an approximate solution of the mixed-fidelity low-boom MDO problem. The BCO method is successfully applied to generate two CFD-based low-boom configurations that closely match two reversed equivalent area targets with ground noise levels below 70 PLdB, respectively. Moreover, these CFD-based low-boom configurations can be obtained with minor wing modifications from the solutions of the low-fidelity MDO problem. The low-fidelity MDO solutions are Pareto points for constrained multiobjective optimization of MTOGW, low-boom cruise Mach, and low-boom range. The best low-boom concept has a predetermined fuselage shape tailored for passengers and main gear storage, carries 40 passengers at seat pitch of 48 in, flies a low-boom overland mission with cruise Mach of 1.8 and range of 2,950 nm, cruises overwater at Mach 1.8 with range of 3,600 nm, satisfies the specified constraints for landing/cruise/takeoff, has MTOGW of 145,164 lb, trims the low-boom cruise flight with fuel redistributions instead of control surface deflections, and has a reversed equivalent area distribution closely matching a target with ground noise level below 70 PLdB.

low-boom supersonic transports↗

NASA NTRS1975

Nozzle Free Jet Flows Within the Strong Curved Shock Regime

A study based on inviscid analysis was conducted to examine the flow field produced from a convergent-divergent nozzle when a strong curved shock occurs. It was found that a certain constraint is imposed on the flow solution of the problem which is the unique feature of the flow within this flow regime, and provides the reason why the inverse method of calculation cannot be employed for these problems. An approximate method was developed to calculate the flow field, and results were obtained for two-dimensional flows. Analysis and calculations were performed for flows with axial symmetry. It is shown that under certain conditions, the vorticity generated at the jet boundary may become infinite and the viscous effect becomes important. Under other conditions, the asymptotic free jet height as well as the corresponding shock geometry were determined.

Shih, Tso-Shin↗

NASA NTRS1986

Solution of the symmetric eigenproblem AX=lambda BX by delayed division

Delayed division is an iterative method for solving the linear eigenvalue problem AX = lambda BX for a limited number of small eigenvalues and their corresponding eigenvectors. The distinctive feature of the method is the reduction of the problem to an approximate triangular form by systematically dropping quadratic terms in the eigenvalue lambda. The report describes the pivoting strategy in the reduction and the method for preserving symmetry in submatrices at each reduction step. Along with the approximate triangular reduction, the report extends some techniques used in the method of inverse subspace iteration. Examples are included for problems of varying complexity.

Thurston, G. A.↗

DOE OSTI2022

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗

NASA NTRS1992

Robust dynamic inversion control laws for aircraft control

Dynamic inversion is a technique for control law design in which feedback is used to simultaneously cancel system dynamics and achieve desired dynamic response characteristics. However, dynamic inversion control laws lack robustness to modeling errors if improperly designed. This paper examines a simple linear example, control of roll rate about the body axis of high performance aircraft, to illustrate some robustness problems which may occur with a simple dynamic inversion control law. The paper demonstrates how structured singular value synthesis techniques can be used to enhance the robustness properties of the dynamic inversion controller.

Balas, Gary J.↗

DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Quantum Time-Space Tradeoffs for Matrix Problems

We consider the time and space required for quantum computers to solve a wide variety of problems involving matrices, many of which have only been analyzed classically in prior work. Our main results show that for a range of linear algebra problems—including matrix-vector product, matrix inversion, matrix multiplication and powering—existing classical time-space tradeoffs, several of which are tight for every space bound, also apply to quantum algorithms with at most a constant factor loss. For example, for almost all fixed matrices 𝐴, including the discrete Fourier transform matrix, we prove that quantum circuits with at most 𝑇 input queries and 𝑆 qubits of memory require 𝑇 = Ω⁢(𝑛 2 /𝑆) to compute matrix-vector product 𝐴⁢𝑥 for 𝑥 ∈{0,1 𝑛 . We similarly prove that matrix multiplication for 𝑛 ×𝑛 binary matrices requires 𝑇 = Ω⁢(𝑛 3 /$\sqrt{𝑆}$). Because many of our lower bounds are matched by deterministic algorithms with the same time and space complexity, our results show that quantum computers cannot provide any asymptotic advantage for these problems with any space bound. We obtain matching lower bounds for the stronger notion of quantum cumulative memory complexity—the sum of the space per layer of a circuit. We also consider Boolean (i.e., AND-OR) matrix multiplication and matrix-vector products, improving the previous quantum time-space tradeoff lower bounds for 𝑛 × 𝑛 Boolean matrix multiplication to 𝑇 = Ω⁢(𝑛 2.5 /𝑆 1/4 ) from 𝑇 = Ω⁢(𝑛 2.5 /𝑆 1/2 ). Our improved lower bound for Boolean matrix multiplication is based on a new coloring argument that extracts more from the strong direct product theorem that was the basis for prior work. To obtain our tight lower bounds for linear algebra problems, we require much stronger bounds than strong direct product theorems. We obtain these bounds by adding a new bucketing method to the quantum recording-query technique of Zhandry that lets us apply classical arguments to upper bound the success probability of quantum circuits.

lower bounds↗

DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2010

Inversion Method for Early Detection of ARES-1 Case Breach Failure

A document describes research into the problem of detecting a case breach formation at an early stage of a rocket flight. An inversion algorithm for case breach allocation is proposed and analyzed. It is shown how the case breach can be allocated at an early stage of its development by using the rocket sensor data and the output data from the control block of the rocket navigation system. The results are simulated with MATLAB/Simulink software. The efficiency of an inversion algorithm for a case breach location is discussed. The research was devoted to the analysis of the ARES-l flight during the first 120 seconds after the launch and early prediction of case breach failure. During this time, the rocket is propelled by its first-stage Solid Rocket Booster (SRB). If a breach appears in SRB case, the gases escaping through it will produce the (side) thrust directed perpendicular to the rocket axis. The side thrust creates torque influencing the rocket attitude. The ARES-l control system will compensate for the side thrust until it reaches some critical value, after which the flight will be uncontrollable. The objective of this work was to obtain the start time of case breach development and its location using the rocket inertial navigation sensors and GNC data. The algorithm was effective for the detection and location of a breach in an SRB field joint at an early stage of its development.

Mackey, Ryan M.↗

NASA NTRS1978

The impact of chlorofluoromethane and NO/x/ injections on stratospheric ozone

The chlorofluoromethane problem and the SST problem have been subject to continuing evaluation by the use of one-dimensional stratospheric photochemistry models. The problems to some extent appear to be inversely related to one another through hydroxyl radical reactions; in other words, processes which tend to increase the amount of OH in a model tend to increase the sensitivity of the ozone column to ClX perturbations, while at the same time decreasing the sensitivity to NO(x) perturbations. Evaluations of sensitivity to ClX and NO(x) have varied as new input information has become available and modeling techniques have improved. It is clear that either NO(x) or ClX, or both, will affect the distribution of ozone with altitude.

Stolarski, R. S.↗

NASA NTRS1983

A noncanonical analytic solution to the J2 perturbed two-body problem

The motion of a satellite subject to an inverse-square gravitational force of attraction and a perturbation due to the earth's oblateness as the J2 term is analyzed, and a uniform, analytic solution correct to first-order in J2, is obtained using a noncanonical approach. The basis for the solution is the transformation and uncoupling of the differential equations for the model. The resulting solution is expressed in terms of elementary functions of the independent variable (the 'true anomaly'), and is of a compact and simple form. Numerical results are comparable to existing solutions.

Jezewski, D. J.↗

NASA NTRS1985

Effects of solar magnetic field on cosmic rays

Aspects of the problem of galactic cosmic ray propagation, including inversion of the solar total magnetic field and an analysis of data related to the heliomagnetic cycle are discussed. It is noted that the global structure of the solar magnetic field results in an additional flux of galactic cosmic rays generated by curvature and gradient drifts. An analysis of heliomagnetic cycle data shows that the latitudinal gradient results in a N-S asymmetry, with the amplitude of the effect growing with depth in the atmosphere. The inversion of the solar total magnetic field, drift effects, and other space distributions are found to contribute to a 22-year cycle of solar activity.

Goncher, G. A.↗

NASA NTRS1994

Coupled 2-dimensional cascade theory for noise and unsteady aerodynamics of blade row interaction in turbofans. Volume 1: Theory development and parametric studies

Typical analytical models for interaction between rotor and stator in a turbofan analyze the effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a separate calculation. In those models, there is a one-to-one relationship between wake harmonics and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the BPF noise harmonic, etc. This report presents a more complete model in which flow tangency boundary conditions are satisfied on two cascades in relative motion for several harmonics simultaneously. By an extension of S.N. Smith's code for two dimensional flat plate cascades, the noise generation/frequency scattering/blade row reflection problem is solved in a single matrix inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy in the higher BPF harmonics due to relative motion between the blade rows. Furthermore, when swirl between the rotor and stator is modeled, a 'mode trapping' effect occurs which explains observations on fans operating at rotational speeds below BFP cuton: the BPF mode amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. However, energy scattered into higher harmonics does propagate and dominates the spectrum at two and three times BPF. This report presents the complete derivation of the theory, comparison with a previous (more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either upstream or downstream of the stator and to counter rotation stages. The theory has been coded in a FORTRAN program called CUP2D, documented in Volume 2 of this report. It is concluded that the new features of this analysis - unsteady coupling, frequency scattering, and flow turning between rotor and stator - have a profound effect on noise generation caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate for noise analysis and prediction.

Hanson, Donald B.↗

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