Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tikhonov regularization”

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.

28 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Extinction and Backscatter retrievals for High Spectral Resolution Lidar for high noise situations

High Spectral Resolution Lidars are uniquely valuable for characterizing atmospheric aerosol and clouds, both alone and in combination with passive instruments such as polarimeters. The benefits of HSRL are the ability to characterize the vertical variability of the atmosphere, to a greater degree than any passive instruments, and – compared to other lidar—additional channels of measurements that permit the direct retrieval of particulate extinction. HSRLs have been successfully deployed from both ground-based and airborne platforms, and will soon be included on satellites as well. From space, new challenges await, due to the need for fully au-tonomous cost-effective instrumentation and lower signal-to-noise ratios. Accordingly, algorithms for the retrieval of HSRL backscatter and extinction must also be reexamined and optimized for best performance in these more challenging circumstances. Motivated by the need for algorithms for the Clio lidar on the future Atmosphere Observing System polar platform, we explore and compare three methods of aerosol extinction retrieval for HSRL data: the direct differentiation method used in standard processing of LaRC airborne HSRL, optimal estimation, and Twomey-Tikhonov regularization. Pros and cons of each methodology will be dis-cussed and explored.

Sharon P. Burton↗
NASA NTRS2024

Decade-Long Ozone Profile Record from Suomi NPP OMPS Limb Profiler: Assessment of Version 2.6 Data

We evaluate a decadal ozone profile record derived from the Suomi National Polar-orbiting Partnership (SNPP) Ozone Mapping and Profiler Suite (OMPS) Limb Profiler (LP) satellite instrument. In 2023, the OMPS LP data were re-processed with the new version 2.6 retrieval algorithm that combines measurements from the ultraviolet (UV) and visible (VIS) parts of the spectra and employs the second order Tikhonov regularization to retrieve a single vertical ozone profile between 12.5 km (or cloud tops) and 57.5 km with the vertical resolution of about 1.9 - 2.5 km between 20-55 km. The algorithm uses radiances measured at six UV ozone-sensitive wavelengths (295, 302, 306, 312, 317 and 322 nm) paired with 353 nm, and one VIS wavelength at 606 nm combined with 510 nm and 675 nm to form a triplet. Each wavelength pair or triplet is used over a limited range of tangent altitudes where the sensitivity to ozone changes are strongest. A new implemented aerosol correction scheme is based on a gamma-function particle size distribution. Numerous calibration changes that affected ozone retrievals were also applied to measured LP radiances, including updates in altitude registration, radiometric calibration, stray light, and spectral registration. The key version 2.6 improvement is the reduction in relative drifts between LP ozone and correlative measurements, linked previously to a drift in the version 2.5 LP altitude registration. We compare LP ozone profiles with those from Aura Microwave Limb Sensor (MLS) to quantify ozone changes in version 2.6.

Limb scattering↗
NASA NTRS1992

Inverse problems: Fuzzy representation of uncertainty generates a regularization

In many applied problems (geophysics, medicine, and astronomy) we cannot directly measure the values x(t) of the desired physical quantity x in different moments of time, so we measure some related quantity y(t), and then we try to reconstruct the desired values x(t). This problem is often ill-posed in the sense that two essentially different functions x(t) are consistent with the same measurement results. So, in order to get a reasonable reconstruction, we must have some additional prior information about the desired function x(t). Methods that use this information to choose x(t) from the set of all possible solutions are called regularization methods. In some cases, we know the statistical characteristics both of x(t) and of the measurement errors, so we can apply statistical filtering methods (well-developed since the invention of a Wiener filter). In some situations, we know the properties of the desired process, e.g., we know that the derivative of x(t) is limited by some number delta, etc. In this case, we can apply standard regularization techniques (e.g., Tikhonov's regularization). In many cases, however, we have only uncertain knowledge about the values of x(t), about the rate with which the values of x(t) can change, and about the measurement errors. In these cases, usually one of the existing regularization methods is applied. There exist several heuristics that choose such a method. The problem with these heuristics is that they often lead to choosing different methods, and these methods lead to different functions x(t). Therefore, the results x(t) of applying these heuristic methods are often unreliable. We show that if we use fuzzy logic to describe this uncertainty, then we automatically arrive at a unique regularization method, whose parameters are uniquely determined by the experts knowledge. Although we start with the fuzzy description, but the resulting regularization turns out to be quite crisp.

Kreinovich, V.↗
DOE OSTI2021

Bayesian Spatially Varying Multi-Regularization Image Deblurring

Many scientific experiments such as those found in astronomy, geology, microbiology, and X-ray radiography require the use of high-energy instruments to capture images. Since blur and noise are inevitably present in any imaging system, the images must be \deblurred" to extract the full information content. Mathematically, image deblurring is an ill-posed inverse problem that requires regularization. The regularization, in turn, has a large effect on the deblurred image: different regularization strengths, and types, lead to drastically different reconstructions. Moreover, many images contain a mixture of smooth and sharp features which suggests the use of multi-regularization, i.e., varying the type of regularization (e.g. Tikhonov or total variation) across the image. We address these issues by formulating the image deblurring problem within a hierarchical Bayesian framework in which we spatially adapt the strength of the regularization and also vary the regularization type across the image. In this way, the image itself, along with corresponding regularization strength at each pixel, are described jointly by a posterior distribution which we can sample by Markov chain Monte Carlo (MCMC) methods. We illustrate our techniques on simplified test problems and apply them to high-energy X-ray images taken at the Nevada National Security Site. Numerical tests show that our new method is robustly applicable and increases the quality of the image reconstruction when compared to other (Bayesian) methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Regularized Variance-Reduced Modified Extragradient Method for Stochastic Hierarchical Games

We consider an N -player hierarchical game in which the i th player’s objective comprises of an expectation-valued term, parametrized by rival decisions, and a hierarchical term. Such a framework allows for capturing a broad range of stochastic hierarchical optimization problems, Stackelberg equilibrium problems, and leader-follower games. We develop an iteratively regularized and smoothed variance-reduced modified extragradient framework for iteratively approaching hierarchical equilibria in a stochastic setting. We equip our analysis with rate statements, complexity guarantees, and almost-sure convergence results. We then extend these statements to settings where the lower-level problem is solved inexactly and provide the corresponding rate and complexity statements. Our model framework encompasses many game theoretic equilibrium problems studied in the context of power markets. We present a realistic application to the study of virtual power plants, emphasizing the role of hierarchical decision making and regularization. Preliminary numerics suggest that empirical behavior compares well with theoretical guarantees.

Tikhonov regularization↗
DOE OSTI2022

Bayesian operator inference for data-driven reduced-order modeling

This work proposes a Bayesian inference method for the reduced-order modeling of time-dependent systems. Informed by the structure of the governing equations, the task of learning a reduced-order model from data is posed as a Bayesian inverse problem with Gaussian prior and likelihood. The resulting posterior distribution characterizes the operators defining the reduced-order model, hence the predictions subsequently issued by the reduced-order model are endowed with uncertainty. The statistical moments of these predictions are estimated via a Monte Carlo sampling of the posterior distribution. Since the reduced models are fast to solve, this sampling is computationally efficient. Furthermore, the proposed Bayesian framework provides a statistical interpretation of the regularization term that is present in the deterministic operator inference problem, and the empirical Bayes approach of maximum marginal likelihood suggests a selection algorithm for the regularization hyperparameters. The proposed method is demonstrated on two examples: the compressible Euler equations with noise-corrupted observations, and a single-injector combustion process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Automated shaker placement and regularized input estimation for MIMO testing.

Multi-input, multi-output (MIMO) testing is used in component qualification to reproduce operational responses in the laboratory. It is often preferred to single-input and base-shake testing because of the potential for equivalent or better tests using smaller actuators and shorter test suites. Given a target response, two key steps in MIMO test design are selecting actuator locations and solving for input loads. Actuator locations are often manually selected using expert judgment. If an automatic method is used, locations are usually determined by simulating the vibration control problem and minimizing a combination of the input energy and control residuals. To select a configuration, the relative importance of input energy and residuals must be specified. Specifying relative weights is, in general, a manual and subjective process. This paper develops an objective function that compares actuator configurations based on control accuracy and required input energy without any manual parameter tuning. The objective function uses an optimally selected tradeoff parameter for each candidate configuration. To choose actuator locations using the new objective function, a pivoting algorithm for integer programming problems is developed. Starting with an initial configuration (such as the one generated by a greedy algorithm), the pivoting algorithm guarantees an objective function decrease in each iteration until convergence is reached. In a simulation featuring a structure excited by a diffuse acoustic field, electrodynamic shaker locations and regularized inputs are solved for without any analyst-specified parameters. Simulations are performed in MIMO configurations where the number of target responses is less than, equal to, and greater than the number of actuators.

Multi-input multi-output↗
NASA NTRS1992

Estimation of local planetary gravity fields using line of sight gravity data and an integral operator

A novel method is presented for mapping line-of-sight gravity data (LOSGD) joining planetary probes and observers during Doppler tracking operations, with a view to geodetic and geophysical applications. LOSGD are in this case mapped as gravity anomalies along a radial direction, at constant altitude, using an inversion procedure in conjunction with a Tikhonov-Arsenine regularization method. The application of different regularization-parameter choices to a synthetic case is followed by application to the real case of Pioneer-Venus orbiter data for Venus' Gula Mons.

Barriot, J. P.↗
DOE OSTI2022

DIAS: A Data-Informed Active Subspace Regularization Framework for Inverse Problems

This paper presents a regularization framework that aims to improve the fidelity of Tikhonov inverse solutions. At the heart of the framework is the data-informed regularization idea that only data-uninformed parameters need to be regularized, while the data-informed parameters, on which data and forward model are integrated, should remain untouched. We propose to employ the active subspace method to determine the data-informativeness of a parameter. The resulting framework is thus called a data-informed (DI) active subspace (DIAS) regularization. Four proposed DIAS variants are rigorously analyzed, shown to be robust with the regularization parameter and capable of avoiding polluting solution features informed by the data. They are thus well suited for problems with small or reasonably small noise corruptions in the data. Furthermore, the DIAS approaches can effectively reuse any Tikhonov regularization codes/libraries. Though they are readily applicable for nonlinear inverse problems, we focus on linear problems in this paper in order to gain insights into the framework. Various numerical results for linear inverse problems are presented to verify theoretical findings and to demonstrate advantages of the DIAS framework over the Tikhonov, truncated SVD, and the TSVD-based DI approaches.

Nguyen, Hai↗
NASA NTRS2017

Optimizing the Regularization in Broadband Wavefront Control Algorithm for WFIRST Coronagraph

The WFIRST/AFTA 2.4 m space telescope currently under study includes a stellar coronagraph for the imaging and the spectral characterization of extrasolar planets. The coronagraph employs sequential deformable mirrors to compensate for phase and amplitude errors. Using the optical model of an Occulting Mask Coronagraph (OMC) testbed at the Jet Propulsion Laboratory (JPL), we have investigated and compared through modeling and simulations the performance of several actuator regularization-schemes in broadband wavefront control algorithm used to generate dark holes in an OMC, such as a Hybrid Lyot Coronagraph (HLC). Using the concept of a Tikhonov filter constituting the G-matrix, we have explained what the different regularization schemes do to singular-modes during a wavefront control process called Electric Field Conjugation (EFC). In some cases we confirmed the numerical predictions with the testbed measured results. We present our findings in this paper.

Sidick, Erkin↗