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 37 records · Page 2

Removing the Effects of Tropical Dynamics from North Pacific Climate Variability

Teleconnections from the Tropics energize variations of the North Pacific climate, but detailed diagnosis of this relationship has proven difficult. Simple univariate methods, such as regression on El Niño-Southern Oscillation (ENSO) indices, may be inadequate since the key dynamical processes involved -- including ENSO diversity in the Tropics, re-emergence of mixed layer thermal anomalies, and oceanic Rossby wave propagation in the North Pacific -- have a variety of overlapping spatial and temporal scales. Here we use a multivariate Linear Inverse Model to quantify tropical and extra-tropical multi-scale dynamical contributions to North Pacific variability, in both observations and CMIP6 models. In observations, we find that the Tropics are responsible for almost half of the seasonal variance, and almost three quarters of the decadal variance, along the North American coast and within the subtropical front region northwest of Hawaii. SST anomalies that are generated by local dynamics within the Northeast Pacific have much shorter time scales, consistent with transient weather forcing by Aleutian low anomalies. Variability within the Kuroshio-Oyashio Extension (KOE) region is considerably less impacted by the Tropics, on all time scales. Consequently, without tropical forcing the dominant pattern of North Pacific variability would be a KOE pattern, rather than the Pacific Decadal Oscillation (PDO). In contrast to observations, most CMIP6 historical simulations produce North Pacific variability that maximizes in the KOE region, with amplitude significantly higher than observed. Correspondingly, the simulated North Pacific in all CMIP6 models is shown to be relatively insensitive to the Tropics, with a dominant spatial pattern generally resembling the KOE pattern, not the PDO.

54 ENVIRONMENTAL SCIENCES↗

The Role of Seasonality and the ENSO Mode in Central and East Pacific ENSO Growth and Evolution

Abstract A cyclostationary linear inverse model (CSLIM) is used to investigate the seasonal growth of tropical Pacific Ocean El Niño–Southern Oscillation (ENSO) events with canonical, central Pacific (CP), or eastern Pacific (EP) sea surface temperature (SST) characteristics. Analysis shows that all types of ENSO events experience maximum growth toward final states occurring in November and December. ENSO events with EP characteristics also experience growth into May and June, but CP events do not. A single dominant “ENSO mode,” growing from an equatorial heat content anomaly into a characteristic ENSO-type SST pattern in about 9 months (consistent with the delayed/recharge oscillator model of ENSO), is essential for the predictable development of all ENSO events. Notably, its seasonality is responsible for the late-calendar-year maximum in ENSO amplification. However, this ENSO mode alone does not capture the observed growth and evolution of diverse ENSO events, which additionally involve the seasonal evolution of other nonorthogonal Floquet modes. EP event growth occurs when the ENSO mode is initially “covered up” in combination with other Floquet modes. The ENSO mode’s slow seasonal evolution allows it to emerge while the other modes rapidly evolve and/or decay, leading to strongly amplifying and more predictable EP events. CP events develop when the initial state has a substantial contribution from Floquet modes with meridional mode–like SST structures. Thus, while nearly all ENSO events involve the seasonally varying ENSO-mode dynamics, the diversity and predictability of ENSO events cannot be understood without identifying contributions from the remaining Floquet modes. Significance Statement The purpose of this study is to identify structures that lead to seasonal growth of diverse types of El Niño–Southern Oscillation (ENSO) events. An important contribution from this study is that it uses an observationally constrained, empirically derived seasonal model. We find that processes affecting the evolution of diverse ENSO events are strongly seasonally dependent. ENSO events with eastern equatorial Pacific sea surface temperature (SST) characteristics are closely related to a single “ENSO mode” that resembles theoretical models of ENSO variability. ENSO events that have central equatorial Pacific SST characteristics include contributions from additional “meridional mode” structures that evolve via different physical processes. These findings are an important step in evaluating the seasonal predictability of ENSO diversity.

Vimont, Daniel J.↗

Investigating the Relative Contribution from Tropical Indo-Pacific SST to Asian Monsoon Precipitation Variability Using LIM

Abstract A critical issue is determining the factors that control the year-to-year variability in precipitation over southern Asia. In this study, we employ a cyclostationary linear inverse model (CS-LIM) to quantify the relative contribution of tropical Pacific and Indian Ocean sea surface temperature anomalies (SSTAs) to the interannual variability of the Asian monsoon, especially Indian summer monsoon rainfall (ISMR). Through a series of CS-LIM experiments, we isolate the impacts of the direct forcing from Pacific SSTAs, Indian Ocean SSTAs, and their interaction on Asian monsoon rainfall variability. Our results reveal distinct patterns of influence with the direct forcing from the Pacific (Indian) Ocean tending to enhance (reduce) the magnitude of precipitation variability, while the Indo-Pacific interaction acts to strongly damp the variability of Asian monsoon precipitation, especially over India. We further investigate these specific impacts on ISMR by analyzing the relationship between tropical Indo-Pacific SSTAs and the leading three empirical orthogonal functions (EOFs) of ISMR. The results from our CS-LIM experiments indicate that the direct forcing from El Niño–Southern Oscillation (ENSO) enhances the variability of the first and third EOFs, while the Indian Ocean SSTA opposes ENSO’s effects, which is consistent with previous studies. Our new results show that the tropical Indo-Pacific interaction strongly damps ISMR variability, which is due to the ENSO-induced Indian Ocean dipole (IOD) opposing the direct impacts from ENSO on ISMR. Additionally, reduced ENSO amplitude and duration associated with the Indo-Pacific interaction may also contribute to the damping effect on ISMR, but this requires further study to understand the relevant mechanisms.

Meteorology & Atmospheric Sciences↗

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗

Simplified Linear Multivariable Control Of Robots

Simplified method developed to design control system that makes joints of robot follow reference trajectories. Generic design includes independent multivariable feedforward and feedback controllers. Feedforward controller based on inverse of linearized model of dynamics of robot and implements control law that contains only proportional and first and second derivatives of reference trajectories with respect to time. Feedback controller, which implements control law of proportional, first-derivative, and integral terms, makes tracking errors converge toward zero as time passes.

Seraji, Homayoun↗

Comparison of some biased estimation methods (including ordinary subset regression) in the linear model

Ridge, Marquardt's generalized inverse, shrunken, and principal components estimators are discussed in terms of the objectives of point estimation of parameters, estimation of the predictive regression function, and hypothesis testing. It is found that as the normal equations approach singularity, more consideration must be given to estimable functions of the parameters as opposed to estimation of the full parameter vector; that biased estimators all introduce constraints on the parameter space; that adoption of mean squared error as a criterion of goodness should be independent of the degree of singularity; and that ordinary least-squares subset regression is the best overall method.

Sidik, S. M.↗

Design of multivariable controllers for robot manipulators

The paper presents a simple method for the design of linear multivariable controllers for multi-link robot manipulators. The control scheme consists of multivariable feedforward and feedback controllers. The feedforward controller is the minimal inverse of the linearized model of robot dynamics and contains only proportional-double-derivative (PD2) terms. This controller ensures that the manipulator joint angles track any reference trajectories. The feedback controller is of proportional-integral-derivative (PID) type and achieves pole placement. This controller reduces any initial tracking error to zero as desired and also ensures that robust steady-state tracking of step-plus-exponential trajectories is achieved by the joint angles. The two controllers are independent of each other and are designed separately based on the linearized robot model and then integrated in the overall control scheme. The proposed scheme is simple and can be implemented for real-time control of robot manipulators.

Seraji, H.↗

Linear multivariable control of robot manipulators

The paper presents a simple scheme for the linear multivariable control of multi-link robot manipulators. The control scheme consists of multivariable feedforward and feedback controllers. The feedforward controller is the minimal inverse of the linearized model of robot dynamics and contains only proportional-double-derivative terms. This controller ensures that the manipulator joint angles track any reference trajectories. The feedback controller is of proportional-derivative type and achieves stability and pole placement. This controller reduces any initial tracking error to zero as desired. The two controllers are independent of each other and are designed separately based on the linearized robot model and then integrated in an overall control scheme.

Seraji, H.↗

An approach to multivariable control of manipulators

The paper presents simple schemes for multivariable control of multiple-joint robot manipulators in joint and Cartesian coordinates. The joint control scheme consists of two independent multivariable feedforward and feedback controllers. The feedforward controller is the minimal inverse of the linearized model of robot dynamics and contains only proportional-double-derivative (PD2) terms - implying feedforward from the desired position, velocity and acceleration. This controller ensures that the manipulator joint angles track any reference trajectories. The feedback controller is of proportional-integral-derivative (PID) type and is designed to achieve pole placement. This controller reduces any initial tracking error to zero as desired and also ensures that robust steady-state tracking of step-plus-exponential trajectories is achieved by the joint angles. Simple and explicit expressions of computation of the feedforward and feedback gains are obtained based on the linearized model of robot dynamics. This leads to computationally efficient schemes for either on-line gain computation or off-line gain scheduling to account for variations in the linearized robot model due to changes in the operating point. The joint control scheme is extended to direct control of the end-effector motion in Cartesian space. Simulation results are given for illustration.

Seraji, H.↗

Method and apparatus for adaptive force and position control of manipulators

The present invention discloses systematic methods and apparatus for the design of real time controllers. Real-time control employs adaptive force/position by use of feedforward and feedback controllers, with the feedforward controller being the inverse of the linearized model of robot dynamics and containing only proportional-double-derivative terms is disclosed. The feedback controller, of the proportional-integral-derivative type, ensures that manipulator joints follow reference trajectories and the feedback controller achieves robust tracking of step-plus-exponential trajectories, all in real time. The adaptive controller includes adaptive force and position control within a hybrid control architecture. The adaptive controller, for force control, achieves tracking of desired force setpoints, and the adaptive position controller accomplishes tracking of desired position trajectories. Circuits in the adaptive feedback and feedforward controllers are varied by adaptation laws.

Seraji, Homayoun↗

An equivalent layer magnetization model for Australia based on Magsat data

An equivalent layer magnetization model for Australia and adjacent oceanic areas is presented. The model is obtained by linear inversion of Magsat anomaly data measured in the altitude range 325-550 km. The anomaly data set has been isolated from the raw data set by use of models of the core field and very long wavelength external fields, and is internally consistent. Certain major structural features of the Australian continent are geographically associated with magnetization anomalies. A first-order difference is seen between the Tasman Zone and the Precambrian cratonic areas: magnetization anomalies are much more subdued in the former, possibly reflecting a shallowing of the Curie isotherm within the crust. A profile of the vertical integral of magnetization is presented for a crustal section extending from the Gawler Block to the southeast coast. It is shown that the magnetization variations are probably due partly, but not wholly, to depth to Curie isotherm variations; gross magnetization variations among at least three distinct crustal units must be involved.

Mayhew, M. A.↗

Steady induction effects in geomagnetism. Part 1C: Geomagnetic estimation of steady surficial core motions: Application to the definitive geomagnetic reference field models

In the source-free mantle/frozen-flux core magnetic earth model, the non-linear inverse steady motional induction problem was solved using the method presented in Part 1B. How that method was applied to estimate steady, broad-scale fluid velocity fields near the top of Earth's core that induce the secular change indicated by the Definitive Geomagnetic Reference Field (DGRF) models from 1945 to 1980 are described. Special attention is given to the derivation of weight matrices for the DGRF models because the weights determine the apparent significance of the residual secular change. The derived weight matrices also enable estimation of the secular change signal-to-noise ratio characterizing the DGRF models. Two types of weights were derived in 1987-88: radial field weights for fitting the evolution of the broad-scale portion of the radial geomagnetic field component at Earth's surface implied by the DGRF's, and general weights for fitting the evolution of the broad-scale portion of the scalar potential specified by these models. The difference is non-trivial because not all the geomagnetic data represented by the DGRF's constrain the radial field component. For radial field weights (or general weights), a quantitatively acceptable explication of broad-scale secular change relative to the 1980 Magsat epoch must account for 99.94271 percent (or 99.98784 percent) of the total weighted variance accumulated therein. Tolerable normalized root-mean-square weighted residuals of 2.394 percent (or 1.103 percent) are less than the 7 percent errors expected in the source-free mantle/frozen-flux core approximation.

Voorhies, Coerte V.↗

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗

Inversion of multiwavelength radiometer measurements by three-dimensional filtering

Remote sensing data from satellites typically have three dimensions: scan position, spacecraft position, and wavelength. Inversion of the radiometric data to infer geophysical parameters is a filtering problem in which the dimension of wavelength (or channel number) is transformed into a dimension of geophysical parameters, and the most general solution is a three-dimensional filter. Linear filters have the advantages of computational speed and easily described transfer functions; but often the measurements are nonlinear functions of the parameters to be inferred. To the extent that the nonlinear inversion problem is overdetermined, it can be modeled by a critically determined linear problem. As an example, inversion of Scanning Multichannel Microwave Radiometer (SMMR) data by means of a three-dimensional Wiener Filter is described. Atmospheric water vapor content, rain liquid water content, surface wind speed and surface temperature are the parameters inferred from the measurements. Nonprecipitating liquid water and water vapor scale height are also modeled but not retrieved. The a priori statistics on which the filter is trained have the effect of governing the selection of a trade-off point of noise as a function of resolution (in all three retrieval dimensions).

Rosenkranz, P. W.↗

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗

Machine learning based inverse modeling of full-field strain distribution for mechanical characterization of a linear elastic and heterogeneous membrane

Heterogeneous membranes or films are thin and soft structures with spatial variations in material property and thickness. Mechanical behavior of heterogeneous membranes is not well understood, mainly due to the difficulty in obtaining accurate and reliable material property data. To understand the mechanical behavior of these materials, accurate and efficient characterization methods for heterogeneous membranes are needed. Here, in this paper, an inverse method based on machine learning is developed to efficiently extract mechanical properties from full-field strain distributions. This approach is demonstrated on a flat heterogeneous membrane with uniform thickness formed by up to four linear elastic synthetic materials in a grid arrangement, and deforming in a moderate strain range (true strain ~10%). The results show that the machine learning method achieves accuracy comparable to the traditional inverse finite element method, and is 6 orders of magnitude faster in the demonstrated case studies.

36 MATERIALS SCIENCE↗

Microbial community structure at the U.S.-Joint Global Ocean Flux Study Station ALOHA: Inverse methods for estimating biochemical indicator ratios

Modeling biogeochemical fluxes in the marine plankton requires the application of factors for extrapolation of biomass indicators measured in the field (chlorophyll a, adenosine triphosphate, bacterial counts) to biomass carbon or nitrogen. These are often inferred from culture studies and are poorly constrained for natural populations. At least squares inverse method with a simple linear model constrains the values of several common indicator ratios, giving self-consistent solutions that provide useful information about the structure of the microbial community at our North Pacific Ocean study site (Station ALOHA (A Long-term Oligotrophic Habitat Assessment)). These results indicate that the fraction of the microbial biomass that is autotrophic (pigmented) is greater in the mixed layer than at the deep chlorophyll maximum layer and that heterotrophic bacteria are a significant but not necessarily predominant component of the microbial community in the euphotic zone.

Christian, James R.↗