Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Comparison of the Tangent Linear Properties of Tracer Transport Schemes Applied to Geophysical Problems.

A number of geophysical applications require the use of the linearized version of the full model. One such example is in numerical weather prediction, where the tangent linear and adjoint versions of the atmospheric model are required for the 4DVAR inverse problem. The part of the model that represents the resolved scale processes of the atmosphere is known as the dynamical core. Advection, or transport, is performed by the dynamical core. It is a central process in many geophysical applications and is a process that often has a quasi-linear underlying behavior. However, over the decades since the advent of numerical modelling, significant effort has gone into developing many flavors of high-order, shape preserving, nonoscillatory, positive definite advection schemes. These schemes are excellent in terms of transporting the quantities of interest in the dynamical core, but they introduce nonlinearity through the use of nonlinear limiters. The linearity of the transport schemes used in Goddard Earth Observing System version 5 (GEOS-5), as well as a number of other schemes, is analyzed using a simple 1D setup. The linearized version of GEOS-5 is then tested using a linear third order scheme in the tangent linear version.

4DVAR↗
NASA NTRS1981

An analytic solution to the classical two-body problem with drag

The two body problem was analyzed with a specific drag model. The model treats drag as a force proportional to the vector velocity and inversely proportional to the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time of flight is expressed as a quadrature in the true anomaly. The results are: (1) development of a vector differential equation which allows analysis of an infinite number of gravitational and drag models; and (2) obtaining the solution of a linear differential equation using the inverse method of laplace transforms.

Mittleman, D.↗
NASA NTRS1974

Numerical modelling of instantaneous plate tectonics

Assuming lithospheric plates to be rigid, 68 spreading rates, 62 fracture zones trends, and 106 earthquake slip vectors are systematically inverted to obtain a self-consistent model of instantaneous relative motions for eleven major plates. The inverse problem is linearized and solved iteratively by a maximum-likelihood procedure. Because the uncertainties in the data are small, Gaussian statistics are shown to be adequate. The use of a linear theory permits (1) the calculation of the uncertainties in the various angular velocity vectors caused by uncertainties in the data, and (2) quantitative examination of the distribution of information within the data set. The existence of a self-consistent model satisfying all the data is strong justification of the rigid plate assumption. Slow movement between North and South America is shown to be resolvable.

Minster, J. B.↗
NASA NTRS1994

Robust stability of second-order systems

Nonlinear control using feedback linearization or inverse dynamics for robotic manipulators yields good results in the absence of modeling uncertainty. However, modeling uncertainties due to unknown joint friction coefficients and payload variations can give rise to undesirable characteristics when these control systems are implemented. It is shown how passivity concepts can be used to supplement the feedback linearization control design technique, in order to make it robust with respect to the uncertain effects mentioned above. Results are obtained for space manipulators with freely floating base; however, they are applicable to fixed base manipulators as well. The controller guarantees asymptotic tracking of the joint variables. Closed-loop simulation results are illustrated for planar space manipulators for cases where uncertainty exists in friction modeling and payload inertial parameters.

Chuang, C.-H.↗
NASA NTRS2010

A Nonlinear Dynamic Inversion Predictor-Based Model Reference Adaptive Controller for a Generic Transport Model

Presented here is a Predictor-Based Model Reference Adaptive Control (PMRAC) architecture for a generic transport aircraft. At its core, this architecture features a three-axis, non-linear, dynamic-inversion controller. Command inputs for this baseline controller are provided by pilot roll-rate, pitch-rate, and sideslip commands. This paper will first thoroughly present the baseline controller followed by a description of the PMRAC adaptive augmentation to this control system. Results are presented via a full-scale, nonlinear simulation of NASA s Generic Transport Model (GTM).

Campbell, Stefan F.↗
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.↗
NASA NTRS1995

Feedback control laws for highly maneuverable aircraft

During this year, we concentrated our efforts on the design of controllers for lateral/directional control using mu synthesis. This proved to be a more difficult task than we anticipated and we are still working on the designs. In the lateral-directional control problem, the inputs are pilot lateral stick and pedal commands and the outputs are roll rate about the velocity vector and side slip angle. The control effectors are ailerons, rudder deflection, and directional thrust vectoring vane deflection which produces a yawing moment about the body axis. Our math model does not contain any provision for thrust vectoring of rolling moment. This has resulted in limitations of performance at high angles of attack. During 1994-95, the following tasks for the lateral-directional controllers were accomplished: (1) Designed both inner and outer loop dynamic inversion controllers. These controllers are implemented using accelerometer outputs rather than an a priori model of the vehicle aerodynamics; (2) Used classical techniques to design controllers for the system linearized by dynamics inversion. These controllers acted to control roll rate and Dutch roll response; (3) Implemented the inner loop dynamic inversion and classical controllers on the six DOF simulation; (4) Developed a lateral-directional control allocation scheme based on minimizing required control effort among the ailerons, rudder, and directional thrust vectoring; and (5) Developed mu outer loop controllers combined with classical inner loop controllers.

Garrard, William L.↗
NASA NTRS1994

Comparative transduction mechanisms of hair cells in the bullfrog utriculus. II. Sensitivity and response dynamics to hair bundle displacement

1. Hair cells in whole-mount in vitro preparations of the utricular macula of the bullfrog (Rana catesbeiana) were selected according to their macular location and hair bundle morphology. The sensitivity and response dynamics of selected hair cells to natural stimulation were examined by recording their voltage responses to step and sinusoidal hair bundle displacements applied to their longest stereocilia. 2. The voltage responses of 31 hair cells to sinusoidal hair bundle displacements were characterized by their gains and phases, taken with respect to peak hair bundle displacement. The gains of Type B and Type C cells at both 0.5 and 5.0 Hz were markedly lower than those of Type F and Type E cells. Phases, with the exception of Type C cells, lagged hair bundle displacement at 0.5 Hz. Type C cells had phase leads of 25-40 degrees. At 5.0 Hz, response phases in all cells were phase lagged with respect to those at 0.5 Hz. Type C cells had larger gains and smaller phase leads at 5.0 Hz than at 0.5 Hz, suggesting the presence of low-frequency adaptation. 3. Displacement-response curves, derived from the voltage responses to 5.0-Hz sinusoids, were sigmoidal in shape and asymmetrical, with the depolarizing response having a greater magnitude and saturating less abruptly than the hyperpolarizing response. When normalized to their largest displacement the linear ranges of these curves varied from < 0.5 to 1.25 microns and were largest in Type B and smallest in Type F and Type E cells. Sensitivity, defined as the slope of the normalized displacement-response curve, was inversely correlated with linear range. 4. The contribution of geometric factors associated with the hair bundle to linear range and sensitivity were predicted from realistic models of utricular hair bundles created using morphological data obtained from light and electron microscopy. Three factors, including 1) the inverse ratio of the lengths of the kinocilium and longest stereocilia, representing the lever arm between kinociliary and stereociliary displacement; 2) tip link extension/linear displacement, largely a function of stereociliary height and separation; and 3) stereociliary number, an estimate of the number of transduction channels, were considered in this analysis. The first of these factors was quantitatively more important than the latter two factors and their total contribution was largest in Type B and Type C cells. Theoretical models were also used to calculate the relation between rotary and linear displacement.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
NASA NTRS1984

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗
DOE OSTI2023

hIPPYlib-MUQ: A Bayesian Inference Software Framework for Integration of Data with Complex Predictive Models under Uncertainty

Bayesian inference provides a systematic framework for integration of data with mathematical models to quantify the uncertainty in the solution of the inverse problem. However, the solution of Bayesian inverse problems governed by complex forward models described by partial differential equations (PDEs) remains prohibitive with black-box Markov chain Monte Carlo (MCMC) methods. We present hIPPYlib-MUQ, an extensible and scalable software framework that contains implementations of state-of-the art algorithms aimed to overcome the challenges of high-dimensional, PDE-constrained Bayesian inverse problems. These algorithms accelerate MCMC sampling by exploiting the geometry and intrinsic low-dimensionality of parameter space via derivative information and low rank approximation. The software integrates two complementary open-source software packages, hIPPYlib and MUQ. hIPPYlib solves PDE-constrained inverse problems using automatically-generated adjoint-based derivatives, but it lacks full Bayesian capabilities. MUQ provides a spectrum of powerful Bayesian inversion models and algorithms, but expects forward models to come equipped with gradients and Hessians to permit large-scale solution. By combining these two complementary libraries, we created a robust, scalable, and efficient software framework that realizes the benefits of each and allows us to tackle complex large-scale Bayesian inverse problems across a broad spectrum of scientific and engineering disciplines. To illustrate the capabilities of hIPPYlib-MUQ, we present a comparison of a number of MCMC methods available in the integrated software on several high-dimensional Bayesian inverse problems. These include problems characterized by both linear and nonlinear PDEs, various noise models, and different parameter dimensions. The results demonstrate that large (~ 50×) speedups over conventional black box and gradient-based MCMC algorithms can be obtained by exploiting Hessian information (from the log-posterior), underscoring the power of the integrated hIPPYlib-MUQ framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

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

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

58 GEOSCIENCES↗
NASA NTRS1980

Remote sensing of ozone in the middle ultraviolet

Earlier studies of the remote sensing of the ozone distribution from satellites in the middle ultraviolet are extended, giving particular attention to the analytic model method. After a brief introduction, the historical background of this subject is presented. This is followed by a summary and review of early work based upon the analytic model. Then a linearized approach to the problem and a linearized equation for a density correction are discussed. These indicate some of the general limitations of the inversion method. Finally the parametrized model approach is reexamined by using both simulated data and some actual satellite data.

Green, A. E. S.↗
NASA NTRS2011

An Inverse Modeling Approach to Estimating Phytoplankton Pigment Concentrations from Phytoplankton Absorption Spectra

Phytoplankton absorption spectra and High-Performance Liquid Chromatography (HPLC) pigment observations from the Eastern U.S. and global observations from NASA's SeaBASS archive are used in a linear inverse calculation to extract pigment-specific absorption spectra. Using these pigment-specific absorption spectra to reconstruct the phytoplankton absorption spectra results in high correlations at all visible wavelengths (r(sup 2) from 0.83 to 0.98), and linear regressions (slopes ranging from 0.8 to 1.1). Higher correlations (r(sup 2) from 0.75 to 1.00) are obtained in the visible portion of the spectra when the total phytoplankton absorption spectra are unpackaged by multiplying the entire spectra by a factor that sets the total absorption at 675 nm to that expected from absorption spectra reconstruction using measured pigment concentrations and laboratory-derived pigment-specific absorption spectra. The derived pigment-specific absorption spectra were further used with the total phytoplankton absorption spectra in a second linear inverse calculation to estimate the various phytoplankton HPLC pigments. A comparison between the estimated and measured pigment concentrations for the 18 pigment fields showed good correlations (r(sup 2) greater than 0.5) for 7 pigments and very good correlations (r(sup 2) greater than 0.7) for chlorophyll a and fucoxanthin. Higher correlations result when the analysis is carried out at more local geographic scales. The ability to estimate phytoplankton pigments using pigment-specific absorption spectra is critical for using hyperspectral inverse models to retrieve phytoplankton pigment concentrations and other Inherent Optical Properties (IOPs) from passive remote sensing observations.

Moisan, John R.↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

The Prospect for Remote Sensing of Cirrus Clouds with a Submillimeter-Wave Spectrometer

Given the substantial radiative effects of cirrus clouds and the need to validate cirrus cloud mass in climate models, it is important to measure the global distribution of cirrus properties with satellite remote sensing. Existing cirrus remote sensing techniques, such as solar reflectance methods, measure cirrus ice water path (IWP) rather indirectly and with limited accuracy. Submillimeter/wave radiometry is an independent method of cirrus remote sensing based on ice particles scattering the upwelling radiance emitted by the lower atmosphere. A new aircraft instrument, the Far Infrared Sensor for Cirrus (FIRSC), is described. The FIRSC employs a Fourier Transform Spectrometer (FTS). which measures the upwelling radiance across the whole submillimeter region (0.1 1.0-mm wavelength). This wide spectral coverage gives high sensitivity to most cirrus particle sizes and allows accurate determination of the characteristic particle size. Radiative transfer modeling is performed to analyze the capabilities of the submillimeter FTS technique. A linear inversion analysis is done to show that cirrus IWP, particle size, and upper-tropospheric temperature and water vapor may be accurately measured, A nonlinear statistical algorithm is developed using a database of 20000 spectra simulated by randomly varying most relevant cirrus and atmospheric parameters. An empirical orthogonal function analysis reduces the 500-point spectrum (20 - 70/cm) to 15 "pseudo-channels" that are then input to a neural network to retrieve cirrus IWP and median particle diameter. A Monte Carlo accuracy study is performed with simulated spectra having realistic noise. The retrieval errors are low for IWP (rms less than a factor of 1.5) and for particle sizes (rins less than 30%) for IWP greater than 5 g/sq m and a wide range of median particle sizes. This detailed modeling indicates that there is good potential to accurately measure cirrus properties with a submillimeter FTS.

Evans, K. Franklin↗
NASA NTRS1989

Least squares restoration of multi-channel images

In this paper, a least squares filter for the restoration of multichannel imagery is presented. The restoration filter is based on a linear, space-invariant imaging model and makes use of an iterative matrix inversion algorithm. The restoration utilizes both within-channel (spatial) and cross-channel information as constraints. Experiments using color images (three-channel imagery with red, green, and blue components) were performed to evaluate the filter's performance and to compare it with other monochrome and multichannel filters.

Chin, Roland T.↗