Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

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 847 records · Page 47

Non‐Linear Dimensionality Reduction With a Variational Encoder Decoder to Understand Convective Processes in Climate Models

Abstract Deep learning can accurately represent sub‐grid‐scale convective processes in climate models, learning from high resolution simulations. However, deep learning methods usually lack interpretability due to large internal dimensionality, resulting in reduced trustworthiness in these methods. Here, we use Variational Encoder Decoder structures (VED), a non‐linear dimensionality reduction technique, to learn and understand convective processes in an aquaplanet superparameterized climate model simulation, where deep convective processes are simulated explicitly. We show that similar to previous deep learning studies based on feed‐forward neural nets, the VED is capable of learning and accurately reproducing convective processes. In contrast to past work, we show this can be achieved by compressing the original information into only five latent nodes. As a result, the VED can be used to understand convective processes and delineate modes of convection through the exploration of its latent dimensions. A close investigation of the latent space enables the identification of different convective regimes: (a) stable conditions are clearly distinguished from deep convection with low outgoing longwave radiation and strong precipitation; (b) high optically thin cirrus‐like clouds are separated from low optically thick cumulus clouds; and (c) shallow convective processes are associated with large‐scale moisture content and surface diabatic heating. Our results demonstrate that VEDs can accurately represent convective processes in climate models, while enabling interpretability and better understanding of sub‐grid‐scale physical processes, paving the way to increasingly interpretable machine learning parameterizations with promising generative properties.

54 ENVIRONMENTAL SCIENCES↗

Factors That Influence Variability in Stress-Drop Measurements Using Spectral Decomposition and Spectral-Ratio Methods for the 2019 Ridgecrest Earthquake Sequence

Stress drop is a fundamental parameter related to earthquake source physics, but is hard to measure accurately. To better understand how different factors influence stress-drop measurements, we compare two different methods using the Ridgecrest stress-drop validation data set: spectral decomposition (SD) and spectral ratio (SR), each with different processing options. Here, we also examine the influence of spectral complexity on source parameter measurement. Applying the SD method, we find that frequency bandwidth and time-window length could influence spectral magnitude calibration, while depth-dependent attenuation is important to correctly map stress-drop variations. For the SR method, we find that the selected source model has limited influence on the measurements; however, the Boatwright model tends to produce smaller standard deviation and larger magnitude dependence than the Brune model. Variance reduction threshold, frequency bandwidth, and time-window length, if chosen within an appropriate parameter range, have limited influence on source parameter measurement. For both methods, wave type, attenuation correction, and spectral complexity strongly influence the result. The scale factor that quantifies the magnitude dependence of stress drop show large variations with different processing options, and earthquakes with complex source spectra deviating from the Brune-type source models tend to have larger scale factor than earthquakes without complexity. Based on these detailed comparisons, we make a few specific suggestions for data processing workflows that could help future studies of source parameters and interpretations.

58 GEOSCIENCES↗

Variational inference of effective range parameters for 3 He− 4 He scattering

We use two different methods, Monte Carlo sampling and variational inference (VI), to perform a Bayesian calibration of the effective-range parameters in 3 He– 4 He elastic scattering. The parameters are calibrated to data from a recent set of 3 He– 4 He elastic scattering differential cross section measurements. Analysis of these data for E lab ≤ 4.3 MeV yields a unimodal posterior for which both methods obtain the same structure. However, the effective-range expansion amplitude does not account for the 7/2 − state of 7 Be so, even after calibration, the description of data at the upper end of this energy range is poor. The data up to E lab = 2.6 MeV can be well described, but calibration to this lower-energy subset of the data yields a bimodal posterior. After adapting VI to treat such a multi-modal posterior we find good agreement between the VI results and those obtained with parallel-tempered Monte Carlo sampling.

effective field theory↗

Variational principle for the parallel-symplectic representation of electromagnetic gyrokinetic theory

The nonlinear (full-f) electromagnetic gyrokinetic Vlasov-Maxwell equations are derived in the parallel-symplectic representation from an Eulerian gyrokinetic variational principle. The gyrokinetic Vlasov-Maxwell equations are shown to possess an exact energy conservation law, which is derived by the Noether method from the gyrokinetic variational principle. Here, the gyrocenter Poisson bracket and the gyrocenter Jacobian contain contributions from the perturbed magnetic field. Lastly, in the full-f formulation of the gyrokinetic Vlasov-Maxwell theory presented here, the gyrocenter parallel-Ampere equation contains a second-order contribution to the gyrocenter current density that is derived from the second-order gyrocenter ponderomotive Hamiltonian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microwave and Millimeter Wave Imaging Using Synthetic Aperture Focusing and Holographical Techniques

Microwave and millimeter wave nondestructive testing and evaluation (NDT&E) methods have shown great potential for determining material composition in composite structures, determining material thickness or debond thickness between two layers, and determining the location and size of flaws, defects, and anomalies. The same testing methods have also shown great potential to produce relatively high-resolution images of voids inside Spray On Foam Insulation (SOFI) test panels using real focused methods employing lens antennas. An alternative to real focusing methods are synthetic focusing methods. The essence of synthetic focusing is to match the phase of the scattered signal to measured points spaced regularly on a plane. Many variations of synthetic focusing methods have already been developed for radars, ultrasonic testing applications, and microwave concealed weapon detection. Two synthetic focusing methods were investigated; namely, a) frequency-domain synthetic aperture focusing technique (FDSAFT), and b) wide-band microwave holography. These methods were applied towards materials whose defects were of low dielectric contrast like air void in SOFI. It is important to note that this investigation used relatively low frequencies from 8.2 GHz to 26.5 GHz that are not conducive for direct imaging of the SOFI. The ultimate goal of this work has been to demonstrate the capability of these methods before they are applied to much higher frequencies such as the millimeter wave frequency spectrum (e.g., 30-300 GHz).

Case, Joseph Tobias↗

Description of ground and excited electronic states by ensemble density functional method with extended active space

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units.We demonstrate that the newmethod correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New operando method for quantifying the relative half-cycle rates of the NO SCR redox cycle over Cu-exchanged zeolites

An operando experimental method for monitoring relative variations in half-cycle rates of the Cu redox-cycle model of selective-catalytic-reduction (SCR) is demonstrated, along with a kinetic model for predicting the same. Conversion inflections (CI) can occur at SCR onset within Cu/SCR catalysts; transient CI involves initial conversion growth to an intermediate value greater than exists at steady state (SS) before degrading to the SS value. While similar CI has been reported for Fe/SCR catalysts as being due to NH 3 inhibition, Cu/SCR CI is inconsistent with a similar origin. Rather, Cu/SCR CI is due to imbalances in the redox-half-cycle rates, and will occur when reduction is faster than oxidation. The temporal CI shape varies with the redox-cycle kinetic parameters and model formulation, and can be used to study the same. Here, we demonstrate the first measurements of dynamic half-cycle rate balancing at SCR onset, use CI and onset transients to study variations in half-cycle rates under Standard, Fast and NO 2 SCR, and present the first redox-cycle model showing CI. The new operando method and model provides a basis for better understanding the SCR redox process, and determining the individual half-cycles kinetic model formulation and parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimum propellant usage for reaction jet systems of space vehicles

The on-off type control for reaction jet systems is proven to be the optimal fuel scheme. However, due to the nonlinear characteristics of this type control, no direct method of solution is known for this optimal process being applied to the attitude control of space vehicles. This paper will discuss the optimization of fuel usage for an attitude control of space vehicles. A computation technique is developed for the calculation of optimal control law. The method of calculus of variations is applied to the estimation of the changes of performance index as well as terminal constraints. Thus an algorithm is obtained by the steepest descent method. A numerical example is given in the paper.

Liu, T. C.↗

Theoretical and experimental investigation of additive drag

The significance of additive drag is discussed and equations for determining its approximate value are derived. Charts are presented giving values of additive drag for open-nose inlets and for annular-nose inlets with conical flow at the inlet. The effects of variable inlet total-pressure recovery and static pressures on the center body are investigated, and an analytical method of predicting the variation of pressure on the center body with mass-flow ratio is given. Experimental values of additive drag are compared with values predicted by the methods.

AIR INLETS↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

An adaptive synchronous extraction (ASE) method for estimating intensity and footprint of surface urban heat islands: A case study of 254 North American cities

The urban heat island (UHI) effect has attracted great attention due to its potential impacts on rapidly growing urban areas. Using remotely sensed estimates of land surface temperature (LST), a large number of studies have focused on the surface UHI (SUHI) effect, which can be characterized by its two fundamental properties: intensity and footprint. The SUHI intensity reflects the LST difference between the urban area and the background reference area (BRA), and the SUHI footprint indicates the spatial extent influenced by the heat island. Currently, numerous methods have been developed to estimate the SUHI intensity and footprint, but are still greatly challenged by three main issues. Namely, the discrepancy in BRA selection criterion brings great uncertainty to the estimated SUHI intensity, the estimation of SUHI footprint is largely constrained by the predefined models, and the quantification of SUHI effect is potentially influenced by several confounding factors. Here, we proposed an adaptive synchronous extraction (ASE) method, which is capable of adaptively selecting the most optimal BRA while removing the influence of confounding factors, and achieving synchronous estimation of SUHI intensity and footprint. We applied the ASE method to 254 North American cities and conducted an in-depth comparative analysis to discuss its applicability and benefits. The main results include: (1) The ASE method avoids the limitations of existing methods in BRA selection and model presetting, and shows resilience to parameter variations. This makes the ASE method highly applicable to quantify the SUHI intensity and footprint in cities with various thermal characteristics. (2) The ASE method can better highlight the spatial, seasonal and day-night contrasts in the estimated SUHI intensity. This superiority is particularly evident when comparing it to methods based on the equal-area buffer or the simplified urban-extent algorithm. (3) Confounding factors pose non-negligible impacts on the quantification of the SUHI effect. Typically, ignoring the influence of topographic relief or missing LST data can lead to an overall overestimation of the SUHI intensity, while not removing surrounding urban areas will cause some underestimation of the SUHI intensity. In conclusion, overall, the proposed ASE method provides a new generalizable tool for quantifying the SUHI effect, which has great potentials for future studies and urban climate assessments.

54 ENVIRONMENTAL SCIENCES↗

Total variation vs L1 regularization: a comparison of compressive sensing optimization methods for chemical detection

One of the fundamental assumptions of compressive sensing (CS) is that a signal can be reconstructed from a small number of samples by solving an optimization problem with the appropriate regularization term. Two standard regularization terms are the L1 norm and the total variation (TV) norm. We present a comparison of CS reconstruction results based on these two approaches in the context of chemical detection, and we demonstrate that optimization based on the L1 norm outperforms optimization based on the TV norm. Our comparison is driven by CS sampling, reconstruction, and chemical detection in two real-world datasets: the Physical Sciences Inc. Fabry-Perot interferometer sensor multispectral dataset and the Johns Hopkins Applied Physics Lab FTIR-based longwave infrared sensor hyperspectral dataset. Both datasets contain the release of a chemical simulant such as glacial acetic acid, triethyl phosphate, and sulfur hexafluoride. For chemical detection we use the adaptive coherence estimator (ACE) and bulk coherence, and we propose algorithmic ACE thresholds to define the presence or absence of a chemical of interest in both un-compressed data cubes and reconstructed data cubes. The un-compressed data cubes provide an approximate ground truth. We demonstrate that optimization based on either the L1 norm or TV norm results in successful chemical detection at a compression rate of 90%, but we show that L1 optimization is preferable. We present quantitative comparisons of chemical detection on reconstructions from the two methods, with an emphasis on the number of pixels with an ACE value above the threshold.

Hyperspectral imagery, L1-norm, total variation, A↗

Methods of removing a ceramic coating from a substrate

A method for removing a ceramic coating from a substrate is presented. The method includes contacting the ceramic coating with a composition including a fluoride source and nitric acid. A method of forming a component having a variation in saturation magnetization is presented. The method includes masking selected portions of a surface of a metallic component using a ceramic coating to form a masked metallic component; selectively diffusing nitrogen into the metallic component by exposing the masked metallic component to a nitrogen-rich atmosphere; and removing the ceramic coating from the surface of the metallic component by contacting the ceramic coating with a composition including the fluoride source and nitric acid.

36 MATERIALS SCIENCE↗

Methods of making a component with variable magnetization and related components

A method of forming a component having a variation in saturation magnetization is presented. The method includes selectively diffusing nitrogen into a metallic component of a masked metallic component by exposing the masked metallic component to a nitrogen-rich atmosphere. The masked metallic component includes a patterned oxide layer formed on a surface of the metallic component, and the patterned oxide layer includes an oxide of a metal present in the metallic component. A related component is also presented.

36 MATERIALS SCIENCE↗

Spaceborne Simulations of Two Direct-Detection Doppler Lidar Techniques

Direct-detection (or incoherent) lidar is now a proven technique for measuring winds in the atmosphere. Over the last few years, several types of direct-detection lidar have evolved. These methods rely on Fabry-Perot interferometers(also termed etalons) or other narrow-passband filters to provide the required spectral resolution. One method, now called the edge (EDG) technique, uses a sharply-sloping filter and measures changes in the filter transmission caused by Doppler shifting of the laser wavelength. A variation of the EDG method, called the double-edge (DEDG) technique, uses two filters. The molecular DEDG method was first demonstrated by Chanin et al. for stratospheric measurements and more recently Korb et al. successfully demonstrated the aerosol DEDG through the troposphere. A second method, here termed the multi-channel (MC) technique, measures Doppler shifts by observing angular displacement of a Fabry-Perot fringe in a spatially resolving detector. The EDG technique thus employs the Fabry-Perot to convert the frequency shift into an amplitude signal, while the MC technique uses the Fabry-Perot to resolve the spectral signature which is then fitted to determine the centroid. The focus of this presentation is on the DEDG and MC methods because these are viewed as the current state of the art in direct-detection lidar. Successful ground-based demonstrations of direct-detection wind measurements have resulted in proposals for spaceborne systems. With this new emphasis on spaceborne systems comes the need for accurate prediction of spaceborne direct-detection Doppler lidar performance. Previously, the EDG and MC methods have been compared although only for aerosol Doppler systems. A recent paper by McGill and Spinhirne compares the DEDG and MC methods in a non-system specific manner for both the aerosol and molecular Doppler systems. The purpose of this presentation is to extend the previous work of McGill and Spinhirne to examine the performance of spaceborne profiling systems. Particular emphasis will be placed on the molecular systems, as these are viewed as the strength of direct-detection Doppler lidar.

McGill, Matthew J.↗

A simple method for correcting spatially resolved solar intensity oscillation observations for variations in scattered light

A measurement of the intensity distribution in an image of the solar disk will be corrupted by a spatial redistribution of the light that is caused by the earth's atmosphere and the observing instrument. A simple correction method is introduced here that is applicable for solar p-mode intensity observations obtained over a period of time in which there is a significant change in the scattering component of the point spread function. The method circumvents the problems incurred with an accurate determination of the spatial point spread function and its subsequent deconvolution from the observations. The method only corrects the spherical harmonic coefficients that represent the spatial frequencies present in the image and does not correct the image itself.

Jefferies, S. M.↗