Search NASA⌕ Search

SEARCH · Search NASA

Results for “constrained variational analysis”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

TPSAS-NF1676L-17943-DND

SAGE II used the solar occultation technique to measure atmospheric ozone for 21 years between 1984 and 2005 and has continued to play a key role in numerous international assessments. Traditionally, the determination of trends in ozone using SAGE II data has been done via time series analysis of monthly mean data within specific latitude bands. However, given the sparse sampling of SAGE II data, this will introduce temporal and/or spatial biases, particularly in the period after late 2000 when the instrument operated at half-duty cycle. Herein we present a new method for performing time series analysis of unevenly sampled data (in particular SAGE II); namely, we utilize the daily mean data and perform a simultaneous temporal and spatial analysis, making use of a Fourier series to constrain latitudinal variation. This methodology is applied to both the previous (v6.2) and current (v7.0) versions of the SAGE II data. Results are discussed that show the improvements to the data between versions 6.2 and 7.0 that make time series analysis produce reduced residuals as well as the preliminary results of the time series itself, which reveal some of the limitations in currently used proxies for the regression.

Robert P Damadeo↗
NASA NTRS2003

Cloud Modeling Using Field Project Data for the Study of Precipitation Processes

The use of cloud-resolving models (CRMs) in the study of precipitation process and their relation to the large-scale environment can be generally categorized into two approaches. The first approach is so called "cloud ensemble modeling". In this approach, many clouds of different size in various stages of their lifecycles can be present at any model simulation time. Large-scale effects are derived from observations and imposed into the model as the main forcing. The advantage of this approach is that the modeled convection will be forced to have the same intensity, thermodynamic budget and organization as the obserations.This approach will also allow CRMs to perform multi-day or multi-week time integrations. The second approach usually requires initial temperature and water vapor profiles that have a medium to large CAPE, and open lateral boundary conditions are used. The modeled clouds could be termed "self-forced convection". Model improvements, such as in the microphysics, are achieved using the second approach. In cloud ensemble modeling, accurate large-scale advective tendencies for temperature and water vapor are the main forcing for the CRMs. We found that the large-scale advective terms for temperature and water vapor are not always consistent, For example, large-scale forcing could indicate strong drying which would produce cooling in the model through evaporation but not contain large-scale advective heating to compensate. This discrepancy in forcing would cause differences between the observed and modeled latent heating profiles. Good measurements of other quantities (i.e., surface fluxes and radiation) are also required to perform variational objective analysis that computes and minimizes a "cost function" that constrains the difference between the large-scale advective forcing in temperature and water vapor. With self-forced convection, accurate vertical distributions of temperature, moisture (water vapor), and horizontal winds are required. The timing of the measurements relative to cloud development is crucial (i.e., prior to cloud triggering). Microphysical measurements (i.e., the cloud number concentration and size distribution) can also be used in this second approach but are of secondary importance with cloud ensemble modeling. In this paper, data collected during TRMM field campaigns (FCs; i.e., SCSMEX, LBA and KWAJEX) which were aimed at validating TRMM products (i.e., rainfall and the vertical distribution of latent heating) will be used to examine the impact of errors in the initial conditions (e.g.., soundings an large-scale forcing) on simulated rainfall distributions and brightness. Rainfall and precipitations simulated from a CRM will also be compared with those estimated by a schocastic model.

Tao, W.-K.↗
NASA NTRS1983

Analysis of regolith electromagnetic scattering as constrained by high resolution Earth-based measurements of the lunar microwave emission

When high resolution measurements of the phase variation of the lunar disk center brightness temperature revealed that in situ regolith electrical losses were larger than those measured on returned samples by a factor of 1.5 to 2.0 at centimeter wavelengths, the need for a refinement of the regolith model to include realistic treatment of scattering effects was identified. Two distinct scattering regimes are considered: vertial variations in dielectric constant and volume scattering due to subsurface rock fragments. Models of lunar regolith energy transport processes are now at the state for which a maximum scientific return could be realized from a lunar orbiter microwave mapping experiment. A detailed analysis, including the effects of scattering produced a set of nominal brightness temperature spectra for lunar equatorial regions, which can be used for mapping as a calibration reference for mapping variations in mineralogy and heat flow.

Keihm, S. J.↗
DOE OSTI2025

Jensen–Shannon divergence based novel loss functions for Bayesian neural networks

Bayesian neural networks (BNNs) are state-of-the-art machine learning methods that can naturally regularize and systematically quantify uncertainties using their stochastic parameters. Kullback–Leibler (KL) divergence-based variational inference used in BNNs suffer from unstable optimization and challenges in approximating light-tailed posteriors due to the unbounded nature of the KL divergence. To resolve these issues, we formulate a novel loss function for BNNs based on a new modification to the generalized Jensen–Shannon (JS) divergence, which is bounded. In addition, we propose a Geometric JS divergence-based loss, which is computationally efficient since it can be evaluated analytically. We found that the JS divergence-based variational inference is intractable, and hence employed a constrained optimization framework to formulate these losses. Our theoretical analysis and empirical experiments on multiple regression and classification data sets suggest that the proposed losses perform better than the KL divergence-based loss, especially when the data sets are noisy or biased. Specifically, there are approximately 5% and 8% improvements in accuracy for a noise-added CIFAR-10 dataset and a regression dataset, respectively. There is about 13% reduction in false negative predictions of a biased histopathology dataset. Additionally, we quantify and compare the uncertainty metrics for the regression and classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Unsupervised machine learning for unbiased chemical classification in X-ray absorption spectroscopy and X-ray emission spectroscopy

Here we report a comprehensive computational study of unsupervised machine learning for extraction of chemically relevant information in X-ray absorption near edge structure (XANES) and in valence-to-core X-ray emission spectra (VtC-XES) for classification of a broad ensemble of sulphorganic molecules. By progressively decreasing the constraining assumptions of the unsupervised machine learning algorithm, moving from principal component analysis (PCA) to a variational autoencoder (VAE) to t-distributed stochastic neighbour embedding (t-SNE), we find improved sensitivity to steadily more refined chemical information. Surprisingly, when embedding the ensemble of spectra in merely two dimensions, t-SNE distinguishes not just oxidation state and general sulphur bonding environment but also the aromaticity of the bonding radical group with 87% accuracy as well as identifying even finer details in electronic structure within aromatic or aliphatic sub-classes. We find that the chemical information in XANES and VtC-XES is very similar in character and content, although they unexpectedly have different sensitivity within a given molecular class. We also discuss likely benefits from further effort with unsupervised machine learning and from the interplay between supervised and unsupervised machine learning for X-ray spectroscopies. Our overall results, i.e., the ability to reliably classify without user bias and to discover unexpected chemical signatures for XANES and VtC-XES, likely generalize to other systems as well as to other one-dimensional chemical spectroscopies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Unified Theory of Fractional Nonlocal and Weighted Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green's identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green's identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. The proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Saturn meteorology - A diagnostic assessment of thin-layer configurations for the zonal flow

Voyager imaging, infrared, and radio observations for Saturn have been recently interpreted by Smith et al. (1982) as an indication that the jet streams observed at the cloud tops extend to depths greater than the 10,000-bar level. This analysis assumes a maximum latitudinal temperature contrast of a few percent, a mean atmospheric rotation rate at depth given by Saturn's ratio period, and no variation with latitude of the bottom pressure level for the zonal flow system. These assumptions are not, however, firmly constrained by observation. The diagnostic analysis of plausible alternative configurations for Saturn's atmospheric structure demonstrates that a thin weather layer system (confined at mid to high latitudes to levels above 200 bar) cannot be excluded by any of the available observations. A quantitative estimate of the effects of moisture condensation (including the differentiation of mean molecular weight) suggests that these might provide the buoyancy contrasts necessary to support a thin-layer flow provided that Saturn's outer envelope is enriched approximately 10 times in water abundance relative to a solar composition atmosphere and strongly differentiated with latitude at the condensation level.

Allison, M.↗
NASA NTRS1989

A decoupled recursive approach for constrained flexible multibody system dynamics

A variational-vector calculus approach is employed to derive a recursive formulation for dynamic analysis of flexible multibody systems. Kinematic relationships for adjacent flexible bodies are derived in a companion paper, using a state vector notation that represents translational and rotational components simultaneously. Cartesian generalized coordinates are assigned for all body and joint reference frames, to explicitly formulate deformation kinematics under small deformation kinematics and an efficient flexible dynamics recursive algorithm is developed. Dynamic analysis of a closed loop robot is performed to illustrate efficiency of the algorithm.

Lai, Hao-Jan↗
DOE OSTI2021

Towards a Unified theory of Fractional and Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green’s identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green’s identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. Here, the proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

A Variational Analysis of Divergence Profiles Based upon Column-Integrated Mass, Moisture and Energetic Constraints with Satellite-Derived Boundary Fluxes

A diagnostic study is made of the mean global divergent circulation based upon a constrained least action principle that minimizes column-integrated divergent kinetic energy subject to constraints on mass, moisture, available potential energy (ape) and total kinetic energy. The concept of gross moist stability was incorporated in the prescription of Lagrange weight function associated with the ape constraint in order to simulate the net effects of cumulus convective heating in the tropics. The variational analyses were validated satisfactorily with the original NCEP/NCAR-reanalysis divergence fields for the Septembers of 1987 and 1988. Further analyses show that in the tropical ascending regions, the analyzed divergences are dominated by the mass and ape constraints; the moisture constraint is implicitly satisfied while the kinetic energy constraint is highly dependent on the ape constraint. In the subtropical descending regions, the analyzed divergences are dominated by the mass, moisture and kinetic energy constraints; the ape constraint is implicitly satisfied. When the constraint integrals were blended with the satellite-derived boundary flux data from GPCP precipitation and ERBE/SRB radiation estimates, the newly analyzed divergences are significantly stronger than the reanalysis divergences in the areas where the estimates of precipitation rates are higher. With few exceptions, the increases in upper-layer divergences are coupled with nearly equal increases in lower-layer convergences.

Lu, Huei-Lin↗
DOE OSTI2026

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗
DOE OSTI2021

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

Temperature Humidity and Sea Level Pressure Increments Induced by 1DVAR Analysis of GPS Refractivity

The Global Positioning System (GPS) transmitted signals are affected by the atmosphere. Using the radio occultation technique, where a receiver is placed on a low-Earth-orbiting platform. it is possible to perform soundings. by measuring the bending angles of the rays. The information can be converted into atmospheric refractivity. We have developed a one dimensional variational (1DVAR) analysis that uses GPS/MET 1995 refractivity and 6-hour FVDAS (Finite Volume Data Assimilation System) forecasts as background information to constrain the retrievals. The analysis increments are defined as 1DVAR minus background temperature, humidity and sea level pressure. Before assimilating the 1DVAR profiles into the FVDAS. the increments need to be understood. First, some bias could be induced in the retrievals when confronted with actual biased data: second. bias in the back-round could create undesired bias in the retrievals. Anv bias in the analyses will ultimately change the climatology of the model the retrievals will be assimilated into. We relate the increments to the reduction of the difference between observed minus computed refractivity profiles. We also point out the difference in the mean increments using backgrounds which have assimilated either NESDIS TIROS Operational Vertical Sounder (TOVS) operational retrievals or Data Assimilation Office (DAO) TOVS interactive retrievals. The climatology of the model in terms of refractivity is significantly different and this impacts the GPS 1DVAR increments. This testifies that changing the basic load of assimilated data has an influence on the impact the GPS data may have in a DAS.

Poli, Paul↗
NASA NTRS2017

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗
NASA NTRSDate not supplied

Contribution to Solar Brightness of Small-Size Magnetic Elements

Irradiance variability is mostly driven by surface magnetism, each magnetic feature contributing in a different manner. The contribution of small-size magnetic elements observed ubiquitously at high-spatial resolution on the solar photosphere, is still debated, as such features are mostly unresolved on full-disk images employed to model irradiance variability. Understanding the contribution to solar brightness of small-size magnetic elements, especially in quiet regions, will help understanding irradiance variability on the decadal and longer temporal scales, which, in turn, are fundamental to understand the role of the Sun on the Earth’s climate. We present a preliminary study of the brightness of small-size magnetic elements using high-spatial resolution spectro-polarimetric observations acquired with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) during Commissioning Phase 1. We focus on the Intensity contrast vs magnetic field relation, which is a fundamental observable in irradiance studies. Previous studies mostly focused on the FeI 630.1/630.2 nm spectral range and so to contextualize our results we focus here on this spectral region as well. By comparing DKIST results with MURaM magneto-hydrodynamic simulations, and with previous results published in the literature, we conclude that our results are in agreement with studies conducted on HINODE observations. We plan to extend the analysis to UV and IR ViSP observations and VBI blue continuum in order to better constrain the quiet sun contribution to spectral irradiance variations.

SMD↗
DOE OSTI2020

Enriched East Asian oxygen isotope of precipitation indicates reduced summer seasonality in regional climate and westerlies

Speleothem oxygen isotope records over East Asia reveal apparently large and rapid paleoclimate changes over the last several hundred thousand years. However, what the isotopic variation actually represent in terms of the regional climate and circulation is debated. We present an answer that emerges from an analysis of the interannual variation in amount-weighted annual δ 18 O of precipitation over East Asia as simulated by an isotope-enabled model constrained by large-scale atmospheric reanalysis fields. 18 O-enriched years have reduced summer seasonality both in terms of precipitation isotopes and in the large-scale circulation. Changes occur between June and October, where the δ 18 O of precipitation (δ 18 O p ) transitions from the isotopically heavier winter to the lighter summer regime. For 18 O-enriched years, this transition is less pronounced. Variations in precipitation amount alone are insufficient to explain the amount-weighted annual δ 18 O p between 18 O-enriched and 18 O-depleted years. Reduced summer seasonality is also expressed in the low-level monsoonal southerlies and upper-level westerlies; for the latter, the northward migration across the Tibetan Plateau in the summer is less pronounced. Our result thus implicates the westerlies across the plateau as the proximate cause of East Asian paleomonsoon changes, manifested as a modulation of its summer peak.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Seasonal variations of snow depth on Mars

Using topography collected over one martian year from the Mars Orbiter Laser Altimeter on the Mars Global Surveyor (MGS) spacecraft, we have measured temporal changes in the elevation of the martian surface that correlate with the seasonal cycle of carbon dioxide exchange between the surface and atmosphere. The greatest elevation change (1.5 to 2 meters) occurs at high latitudes ( above 80 degrees ), whereas the bulk of the mass exchange occurs at lower latitudes (below 75 degrees N and below 73 degrees S). An unexpected period of sublimation was observed during northern hemisphere autumn, coincident with dust storms in the southern hemisphere. Analysis of MGS Doppler tracking residuals revealed temporal variations in the flattening of Mars that correlate with elevation changes. The combined changes in gravity and elevation constrain the average density of seasonally deposited carbon dioxide to be 910 +/- 230 kilograms per cubic meter, which is considerably denser than terrestrial snow.

Carbon Dioxide↗
DOE OSTI2021

A Constrained Data Assimilation Algorithm Based on GSI Hybrid 3D-EnVar and Its Application

Abstract Data assimilation (DA) at mesoscales is important for severe weather forecasts, yet the techniques of data assimilation at this scale remain a challenge. This study introduces dynamical constraints in the Gridpoint Statistical Interpolation (GSI) three-dimensional ensemble variational (3D-EnVar) data assimilation algorithm to enable the use of high-resolution surface observations of precipitation to improve atmospheric analysis at mesoscales. The constraints use the conservations of mass and moisture. Mass constraint suppresses the unphysical high-frequency oscillation, while moisture conservation constrains the atmospheric states to conform with the observed high-resolution precipitation. We show that the constrained data assimilation (CDA) algorithm significantly reduced the spurious residuals of the mass and moisture budgets compared to the original data assimilation (ODA). A case study is presented for a squall line over the Southern Great Plains on 20 May 2011 during Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program by using ODA or CDA analysis as initial condition of forecasts. The state variables, and the location and intensity of the squall line are better simulated in the CDA experiment. Results show how surface observation of precipitation can be used to improve atmospheric analysis through data assimilation by using the dynamical constraints of mass and moisture conservations.

Wang, Jia↗