Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 163 records · Page 9

Modified Cascading Generalized Inverse Control Allocation

The current aviation revolution towards electric propulsion aircraft (e.g., electric vertical takeoff-and-landing) brings unique control challenges. These vehicles are typically over-actuated (more effectors than desired control outcomes), may require control strategies for the three phases of flight (hover, transition and cruise), and currently have limited electric power availability. These vehicle challenges bring the need for optimal control allocation to the forefront of research. A leading control allocation algorithm, used in current flight vehicles, is the Cascading Generalized Inverse (CGI). Unfortunately, the Cascading Generalized Inverse algorithm is unable to achieve some desired outcomes, it intermittently provides non-optimal allocations, and it may fail to preserve moment direction near maximal achievable outcomes. In this research, the shortcomings of the Cascading Generalized Inverse algorithm are addressed by augmenting the algorithm with Scalar Difference Quadratic unsaturation identification and location at each iteration. Rigorous theory is shown that the Modified Cascading Generalized Inverse performs better at obtaining optimal allocations for all attainable outcomes. Numerical case studies for over-actuated vehicles demonstrate resolution to the aforementioned deficiencies.

Control Allocation↗

Modified Cascading Generalized Inverse Control Allocation

The current aviation revolution towards electric propulsion aircraft (e.g., electric vertical takeoff-and-landing) brings unique control challenges. These vehicles are typically over-actuated (more effectors than desired control outcomes), may require control strategies for the three phases of flight (hover, transition and cruise), and currently have limited electric power availability. These vehicle challenges bring the need for optimal control allocation to the forefront of research. A leading control allocation algorithm, used in current flight vehicles, is the Cascading Generalized Inverse (CGI). Unfortunately, the Cascading Generalized Inverse algorithm is unable to achieve some desired outcomes, it intermittently provides non-optimal allocations, and it may fail to preserve moment direction near maximal achievable outcomes. In this research, the shortcomings of the Cascading Generalized Inverse algorithm are addressed by augmenting the algorithm with Scalar Difference Quadratic unsaturation identification and location at each iteration. Rigorous theory is shown that the Modified Cascading Generalized Inverse performs better at obtaining optimal allocations for all attainable outcomes. Numerical case studies for over-actuated vehicles demonstrate resolution to the aforementioned deficiencies.

Control Allocation↗

A parameterization of the depth of the entrainment zone

A theory of the parameterization of the entrainment zone depth has been developed based on conservation of energy. This theory suggests that the normalized entrainment zone depth is proportional to the inverse square root of the Richardson number. A comparison of this theory with atmospheric observations indicates excellent agreement. It does not adequately predict the laboratory data, although it improves on parcel theory, which is based on a momentum balance.

Boers, Reinout↗

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

Direct Characterization of Type-I Band Alignment in 2D Ruddlesden–Popper Perovskites

2D Ruddlesden–Popper halide perovskites have attracted considerable attention due to their desirable optoelectronic properties, high chemical and structural tunability, and improved environmental stability. However, the understanding of their structure–properties relationships is still limited. In particular, the energy level positions and band alignments at interfaces involving these materials, which are important features to control in the context of any applications, are still under debate. Here, the electronic structure of high-purity films of BA 2 MA n-1 Pb n I 3n+1 for n = 1–5 (where BA stands for butylammonium and MA for methylammonium) is investigated, using optical absorption, ultraviolet, and inverse photoemission spectroscopies, and density functional theory calculations. Further, this study determines the ionization energy and electron affinity of each compound and demonstrates a type-I band alignment for the BA 2 MA n-1 Pb n I 3n+1 series. This study further describes the evolution of the exciton binding energy as a function of the thickness of the inorganic layers.

2D Ruddlesden–Popper perovskites↗

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗

Phases of dark matter from inverse decays

Inverse decays are an interesting avenue for producing dark matter in the early Universe. We study in detail various phases of dark matter parameter space where inverse decays control its abundance, expanding on our work of inverse decay (INDY) dark matter and going beyond. The role of initial conditions and the impact of departure from kinetic equilibrium are investigated as well. We show how these inverse decay phases can arise in theories of a kinetically mixed dark photon and dark Higgs, with promising prospects for detection at upcoming experiments.

Cosmology↗

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗

The Induced Magnetic Field of the Moon: Conductivity Profiles and Inferred Temperature

Electromagnetic induction in the moon driven by fluctuations of the interplanetary magnetic field is used to determine the lunar bulk electrical conductivity. The present data clearly show the north-south and east-west transfer function difference as well as high frequency rollover. The difference is shown to be compatible over the mid-frequency range with a noise source associated with the compression of the local remanent field by solar wind dynamic pressure fluctuations. Models for two, three, and four layer; current layer, double current layer, and core plus current layer moons are generated by inversion of the data using a theory which incorporates higher order multipoles. Core radii conductivities generally are in the range 1200 to 1300 km and 0.001 to 0.003 mhos/m; and for the conducting shell 1500 to 1700 km with 0.0001 to 0.0007 mhos/m with an outer layer taken as nonconducting. Core temperature based on available olivine data is 700 to 1000 C.

Sonett, C. P.↗

Vertical resolution of temperature profiles obtained from remote radiation measurements

The Backus-Gilbert theory, originally developed for analysis of inversion problems associated with the physics of the solid earth, was applied to the problem of the vertical sounding of the atmosphere by means of remote radiation measurements. An application was made to spectral intervals 2.8/cm wide in the 667/cm band CO2, and tradeoff curves are presented which quantitatively define the relationship between intrinsic vertical resolution and random error in temperature profile estimates. It is found that for a 1-2 K random error with state-of-the-art instrumentation, the intrinsic vertical resolution ranges from approximately 0.5 locale scale height (l.s.h.) in the lower troposphere to greater than 2 l.s.h. in the upper stratosphere with approximately 1 l.s.h. resolution in the vicinity of the tropopause. These values are somewhat smaller than the widths of the radioactive transfer kernels at similar levels. Increasing the number of spectral intervals from 7 to 16 is found to produce only a marginal improvement in vertical resolution.

Conrath, B. J.↗

The induced magnetic field of the moon - Conductivity profiles and inferred temperature.

Electromagnetic induction in the moon driven by fluctuations of the interplanetary magnetic field is used to determine the lunar bulk electrical conductivity. The earlier data are now augmented by an order of magnitude. The present data clearly show the north-south and east-west transfer function difference as well as the high-frequency rollover suggested earlier. The difference is shown to be compatible over the midfrequency range (0.001 to 0.01 Hz) with a noise source associated with the compression of the local remanent field by solar wind dynamic pressure fluctuations. The rollover of the transfer functions is shown to result from higher order magnetic multipole radiation; electric multipoles appear supressed, although a vestigial TM interaction may still be present. Models for two, three, and four layer; current layer, double current layer, and core plus current layer moons are generated by inversion of the data, using a theory that incorporates higher-order multipoles.

Sonett, C. P.↗

Surface roughness studies for wind tunnel models used in high Reynolds number testing

This paper focuses on stylus and optical techniques for the measurement of surface roughness in wind tunnel models. The stylus instruments provide detailed information, such as surface profiles and area maps, that may then be used either to calculate statistical properties (i.e., the rms surface roughness) or to study individual surface peaks or other features. By contrast, certain optical techniques yield area-averaged statistical properties of the surface roughness directly. Two instruments that use the technique of optical angular scattering are compared. One is a research instrument that has been developed to study the basic scattering phenomena by testing the optical theories and surface models used in inverse calculations of statistical roughness parameters. The second instrument is more compact and is under development as a hand held, on-line device to be used during manufacture of wind tunnel models for the National Transonic Facility at NASA Langley Research Center. The scattering geometries for the two instruments are compared and results from these instruments and the stylus technique are shown for roughness specimens that are typical of the surface finish of wind tunnel models.

Vorburger, T. V.↗

Viscous induced drag

A fundamentally new approach to the aircraft minimum induced drag problem is presented. The method, a 'viscous lifting line', is based on the minimum entropy production principle and does not require the planar wake assumption. An approximate, closed form solution is obtained and compared with several classical results. In addition, the problem of optimizing in-plane wing sweep with constant wing root bending moment is considered. Like the classical lifting line theory, this theory predicts that induced drag is proportional to the square of the lift coefficient and inversely proportional to the wing aspect ratio. Unlike the classical theory, it predicts that in-plane wing sweep may significantly reduce induced drag, that induced drag is Reynolds number dependent, and that the optimum spanwise circulation distribution is non-elliptic.

Greene, George C.↗

An entropy method for induced drag minimization

A fundamentally new approach to the aircraft minimum induced drag problem is presented. The method, a 'viscous lifting line', is based on the minimum entropy production principle and does not require the planar wake assumption. An approximate, closed form solution is obtained for several wing configurations including a comparison of wing extension, winglets, and in-plane wing sweep, with and without a constraint on wing-root bending moment. Like the classical lifting-line theory, this theory predicts that induced drag is proportional to the square of the lift coefficient and inversely proportioinal to the wing aspect ratio. Unlike the classical theory, it predicts that induced drag is Reynolds number dependent and that the optimum spanwise circulation distribution is non-elliptic.

Greene, George C.↗

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗

Success Stories in Control: Nonlinear Dynamic Inversion Control

NASA plays an important role in advancing the state of the art in flight control systems. In the case of Nonlinear Dynamic Inversion (NDI) NASA supported initial implementation of the theory in an aircraft and demonstration in a space vehicle. Dr. Dale Enns of Honeywell Aerospace Advanced Technology performed this work in cooperation with NASA and under NASA contract. Honeywell and Lockheed Martin were subsequently contracted by AFRL to create "Design Guidelines for Multivariable Control Theory". This foundational work directly contributed to the advancement of the technology and the credibility of the control law as a design option. As a result Honeywell collaborated with Lockheed Martin to produce a Nonlinear Dynamic Inversion controller for the X-35 and subsequently Lockheed Martin did the same for the production Lockheed Martin F-35 vehicle. The theory behind NDI is to use a systematic generalized approach to controlling a vehicle. Using general aircraft nonlinear equations of motion and onboard aerodynamic, mass properties, and engine models specific to the vehicle, a relationship between control effectors and desired aircraft motion can be formulated. Using this formulation a control combination is used that provides a predictable response to commanded motion. Control loops around this formulation shape the response as desired and provide robustness to modeling errors. Once the control law is designed it can be used on a similar class of vehicle with only an update to the vehicle specific onboard models.

Bosworth, John T.↗

Time domain version of the uniform GTD

The uniform geometrical theory of diffraction (UTD) solutions can be inversely transformed analytically to obtain a time-domain version of the UTD. The time-domain solutions are valid in the early time period where an observation time t is close to the time after the arrival of the first diffracted wavefront. Comparisons with GTD (geometrical theory of diffraction) and also with available rigorous results reveal that the UTD solutions are accurate for substantial early time periods while the GTD results are valid for very early time periods.

Veruttipong, Thavath Watt↗