SEARCH · Engineering Papers

Results for “theoretical and computational techniques”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2024

3D Cloud Masking Across A Broad Swath Using Multi-Angle Polarimetry and Deep Learning

Understanding the 3-dimensional structure of clouds is of crucial importance to modeling our changing climate. Active sensors, such as radar and lidar, provide accurate vertical cloud profiles, but are mostly restricted to along-track sampling. Passive sensors can capture a wide swath, but struggle to see beneath cloud tops. In essence, both types of products are restricted to two dimensions: as a cross-section in the active case, and an image in the passive case. However, multi-angle sensor configurations contain implicit information about 3D structure, due to parallax and atmospheric path differences. Extracting that implicit information can be challenging, requiring computationally expensive radiative transfer techniques. Machine learning, as an alternative, may be able to capture some of the complexity of a full 3D radiative transfer solution with significantly less computational expense. In this work, we make three contributions towards understanding 3D cloud structure from multi-angle polarimetry. First, we introduce a large-scale, publicly available dataset that fuses existing cloud products into a format more amenable to machine learning. This dataset treats multi-angle polarimetry as an input, and radar-based vertical cloud profiles as an output. Second, we describe and evaluate strong baseline machine learning models based that predict these profiles from the passive imagery. Notably, these models are trained only on center-swath labels, but can predict cloud profiles over the entire passive imagery swath. Third, we leverage the information-theoretic nature of machine learning to draw conclusions about the relative utility of various sensor configurations, including spectral channels, viewing angles, and polarimetry. These findings have implications for Earth-observing missions such as NASA’s Plankton, Aerosol, Cloud-ocean Ecosystem (PACE) and Atmosphere Observing System (AOS) missions, as well as in informing future applications of computer vision to atmospheric remote sensing.

clouds↗

DOE OSTI2024

Emerging Jets Search, Triton Server Deployment, and Track Quality Development: Machine Learning Applications in High Energy Physics

Machine learning is becoming prevalent in high energy physics, with numerous applications in physics analyses and event reconstruction showing great improvements compared to traditional computing methods. This thesis studies three projects which each propose new avenues for machine learning applications within the high energy physics CMS experiment located at CERN. In the first project, a search for a dark matter signal called “emerging jets” is performed, using graph neural networks to greatly increase sensitivity to the signal’s signature within the data. The result of this dark matter search sets the most stringent exclusion limits to date on theoretical emerging jet models. Motivated by inefficiencies encountered when processing the emerging jet graph neural network at Fermi National Accelerator Laboratory’s computing centers, the second project re-optimizes the computing centers for machine learning inference. This re-optimization uses NVIDIA Triton Inference Servers to process users’ analysis code heterogeneously, therefore achieving high processing throughput and decreasing user time-to-insight. The last project focuses on an upgrade to the CMS experiment’s real-time event selection system which improves physics object reconstruction under harsh processing conditions. A boosted decision tree is used to quickly and efficiently quantify a reconstructed particle’s “track quality” in order to remove particle tracks reconstructed erroneously. In summary, this thesis will not only present examples of how high energy physics can greatly benefit by leveraging machine learning techniques for physics analysis and reconstruction, but will also provide guidance on how the field can prepare for the inevitable increase in machine learning applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1985

Probabilistic Finite Element: Variational Theory

The goal of this research is to provide techniques which are cost-effective and enable the engineer to evaluate the effect of uncertainties in complex finite element models. Embedding the probabilistic aspects in a variational formulation is a natural approach. In addition, a variational approach to probabilistic finite elements enables it to be incorporated within standard finite element methodologies. Therefore, once the procedures are developed, they can easily be adapted to existing general purpose programs. Furthermore, the variational basis for these methods enables them to be adapted to a wide variety of structural elements and to provide a consistent basis for incorporating probabilistic features in many aspects of the structural problem. Tasks concluded include the theoretical development of probabilistic variational equations for structural dynamics, the development of efficient numerical algorithms for probabilistic sensitivity displacement and stress analysis, and integration of methodologies into a pilot computer code.

Belytschko, T.↗

NASA NTRS1999

Estimation of Non-Seasonal Oceanic Mass Redistribution and its Geodynamic Effects from the Topex/Poseidon Observation

Mass redistribution in terms of angular momentum will change the Earth's rotation under the conservation of angular momentum. It also changes the external gravitational field according to Newton's gravitational law. These small geodynamic variations can be measured very accurately by modem space geodetic techniques. We report on new estimates of the oceanic tidal angular momentum and hence the tidal variations in the Earth's rotation rate and polar motion. These are deduced from the new tidal solutions GOT99 and TPXO.3, plus M2 and Mf estimates (with error bars) from TPXO.4a. Comparisons with observed diurnal/semidiurnal tidal variations are generally similar to previous model estimates. Comparisons of models and observations of Mf polar motion are confused, with neither models nor observations agreeing among themselves or with each other. Part of this is caused by differences in Mf currents; model TPXO.4 shows current-intensive Rossby modes along western boundaries of major basins, which are indeed expected on theoretical grounds. Unfortunately, there are no reliable in situ observations of Mf currents useful for comparing models. On non-tidal oceanic mass redistribution we compute the principal components of the multi-year T/P sea surface height (SSH) field using EOF (empirical orthogonal functions) taking into consideration of the latitude-dependent area weighting. Emphasis is placed upon non-seasonal and interannual oscillations on regional scales, including ENSO and the North Atlantic Oscillation. Steric corrections to the observed SSH due to sea surface temperature changes (which has no geodynamic effects) is now under study using POCM ocean model output.

Chao, Benjamin F.↗

DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗

NASA NTRS1977

Prelaunch testing of the laser geodynamic satellite (LAGEOS)

The LAGEOS was extensively tested optically prior to launch. The measurement techniques used are described and resulting data is presented. Principal emphasis was placed on pulse spreading characteristics, range correction for center of mass tracking, and pulse distortion due to coherent effects. A mode-locked freqeuncy doubled Nd:YAG laser with a pulse width of about 60 ps was used as the ranging transmitter and a crossfield photo-multiplier was used in the receiver. High speed sampling electronics were employed to increase receiver bandwidth. LAGEOS reflected pulses typically had a width of 250 ps with a variability in the range correction of less than 2 mm rms. Pulse distortion due to coherent effects was inferred from average waveforms and appears to introduce less than + or - 50 ps jitter in the location of the pulse peak. Analytic results on this effect based on computer simulations are also presented. Theoretical and experimental data on the lidar cross section were developed in order to predict the strength of lidar echoes from the satellite. Cross section was measured using a large aperture laser collimating system to illuminate the LAGEOS. Reflected radiation far-field patterns were measured using the collimator in an autocollimating mode. Data were collected with an optical data digitzer and displayed as a three-dimensional plot of intensity versus the two far-field coordinates. Measurements were made at several wavelengths, for several types of polarizations, and as a function of satellite orientation.

Fitzmaurice, M. W.↗

NASA NTRS2015

Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, most of them are either unable to capture interference effects from bodies, or require an extremely large computational budget. The objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference using a novel Computational Fluid Dynamics (CFD) solver for rotorcraft, RotCFD. In RotCFD, a mid-fidelity URANS solver is used with an incompressible flow model and a realizable k-ε turbulence model. The rotor is, however, not modeled using a computationally expensive, unsteady viscous body-fitted grid, but is instead modeled using a blade element model with a momentum source approach. Various flight modes of the XV-15 isolated rotor, including hover, tilt and airplane mode, have been simulated and correlated to existing experimental and theoretical data. The rotor model is subsequently used for wind tunnel wall interference simulations in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames Research Center in California. The results from the validation of the isolated rotor performance showed good correlation with experimental and theoretical data. The results were on par with known theoretical analyses. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes, which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. For both test sections of the NFAC wall interference was examined by simulating the XV-15 rotor in the test section of the wind tunnel and with an identical grid but extended boundaries in free field. Both cases were also examined with an isolated rotor or with the rotor mounted on the modeled geometry of the Tiltrotor Test Rig (TTR). A 'quasi linear trim' was used to trim the thrust for the rotor to compare the power as a unique variable. Power differences between free field and wind tunnel cases were found from -7 % to 0 % in the 80- by 120-Foot Wind Tunnel test section and -1.6 % to 4.8 % in the 40- by 80-Foot Wind Tunnel, depending on the TTR orientation, tunnel velocity and blade setting. The TTR will be used in 2016 to test the Bell 609 rotor in a similar fashion to the research in this report.

CFD↗

DOE OSTI2022

Theoretical Perspective on Operando Spectroscopy of Fluxional Nanocatalysts

Improvements in operando spectroscopy have enabled the catalysis community to investigate the dynamic nature of catalysts under operating conditions with increasing detail. Still, the highly dynamic nature of some catalysts, such as fluxional supported subnano clusters, presents a formidable challenge even for the most state-of-the-art techniques. Here, the reason is that such fluxional catalytic interfaces contain a variety of thermally accessible states. Operando spectroscopies used in catalysis generally fall into two categories: ensemble-based techniques, which provide spectra containing the signals of the entire ensemble of states of the catalyst and are not necessarily dominated by the most active species, and localized techniques, which provide atomistic-level information about the dynamics of active sites in a very small area, which might not include the most active species. Combining many different kinds of techniques can provide detailed insight; however, we propose that effective utilization of specific computational techniques and approaches within the fluxionality paradigm can fill the gap and enable atomistic characterization of the most relevant catalytic sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1975

Solid-propellant rocket motor ballistic performance variation analyses

Results are presented of research aimed at improving the assessment of off-nominal internal ballistic performance including tailoff and thrust imbalance of two large solid-rocket motors (SRMs) firing in parallel. Previous analyses using the Monte Carlo technique were refined to permit evaluation of the effects of radial and circumferential propellant temperature gradients. Sample evaluations of the effect of the temperature gradients are presented. A separate theoretical investigation of the effect of strain rate on the burning rate of propellant indicates that the thermoelastic coupling may cause substantial variations in burning rate during highly transient operating conditions. The Monte Carlo approach was also modified to permit the effects on performance of variation in the characteristics between lots of propellants and other materials to be evaluated. This permits the variabilities for the total SRM population to be determined. A sample case shows, however, that the effect of these between-lot variations on thrust imbalances within pairs of SRMs is minor in compariosn to the effect of the within-lot variations. The revised Monte Carlo and design analysis computer programs along with instructions including format requirements for preparation of input data and illustrative examples are presented.

Sforzini, R. H.↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

DOE OSTI2024

Identifiability and characterization of transmon qutrits through Bayesian experimental design

Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A Miniaturized, High-Bandwidth Optical Fiber Fabry–Perot Cavity Vibration Sensor Demonstrated up to 800 °C

A typical structural health monitoring technique involves measuring the vibrational characteristics of components or systems to detect signs of degradation or damage. Many industrial applications require engineered systems to safely operate under extreme, high-temperature environments that pose challenges not only to materials but also to sensors that would be used for structural health monitoring. Here, in this study, miniaturized optical Fabry-Perot cavities (FPCs) were developed and tested as a means of measuring the resonant frequencies of metal components that are most relevant to extreme-environment applications. Two of the three candidate FPC designs tested up to 800 ° C provided accurate measurements (validated by theoretical models and laser Doppler vibrometry) of the fundamental vibrational mode of the specimen to which each was bonded, although both sensors failed during thermal cycling. An analysis of the reflected optical spectrum from the FPC and X-ray computed tomography revealed two opportunities to improve the sensor reliability. First, the Cu optical fiber coating that was used could either be replaced with a more oxidation-resistant material or protected with commercially available films. Second, the adhesives used to bond the fibers to metal capillaries and establish the FPC could be replaced with a more robust solution, although the Resbond 907TS adhesive appeared to outperform Resbond 907.

Birri, Anthony [Oak Ridge National Laboratory (ORN↗

NASA NTRS2004

Combining Satellite Data, Trajectory Modeling and Surface Insolation Measurements to Deduce the Direct Radiative Effect of Smoke Aerosol

In this paper, we have introduced a method of inferring the radiative effect of smoke aerosols using a technique that combines satellite remote sensing with trajectory modeling. The results shown here clearly show large flux biases between theoretical and measured radiative fluxes correlate with the arrival of smoke aerosol to the area. Further analysis is required to convincingly demonstrate that the reason for these differences is the radiative effect of the smoke aerosols. To do this, the estimated fluxes taken from the ERA-15 will be recomputed every 3 hours using International Satellite Cloud Climatology Project (ISCCP) data set entitled DX gridded to a 1o equal angle resolution (see paper 7B.2 for details). Surface radiometric and ancillary data for several more Canadian surface sites are being obtained at minute temporal resolution. The ultimate purpose of this research is to derive aerosol smoke maps for fire events such as this to be included in an aerosol climatology and be incorporated in the computation of the earth's surface radiation budget to better understand the radiative effect of aerosols.

Stackhouse, Paul W., Jr.↗

NASA NTRS2008

A Multi-Scale Modeling System with Unified Physics

Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 km, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. In recent years, exponentially increasing computer power has extended cloud-resolving-mode1 integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 km2 in three-dimensions. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique.

Tao, Wei-Kuo↗

NASA NTRS2011

Multiple-Beam Detection of Fast Transient Radio Sources

A method has been designed for using multiple independent stations to discriminate fast transient radio sources from local anomalies, such as antenna noise or radio frequency interference (RFI). This can improve the sensitivity of incoherent detection for geographically separated stations such as the very long baseline array (VLBA), the future square kilometer array (SKA), or any other coincident observations by multiple separated receivers. The transients are short, broadband pulses of radio energy, often just a few milliseconds long, emitted by a variety of exotic astronomical phenomena. They generally represent rare, high-energy events making them of great scientific value. For RFI-robust adaptive detection of transients, using multiple stations, a family of algorithms has been developed. The technique exploits the fact that the separated stations constitute statistically independent samples of the target. This can be used to adaptively ignore RFI events for superior sensitivity. If the antenna signals are independent and identically distributed (IID), then RFI events are simply outlier data points that can be removed through robust estimation such as a trimmed or Winsorized estimator. The alternative "trimmed" estimator is considered, which excises the strongest n signals from the list of short-beamed intensities. Because local RFI is independent at each antenna, this interference is unlikely to occur at many antennas on the same step. Trimming the strongest signals provides robustness to RFI that can theoretically outperform even the detection performance of the same number of antennas at a single site. This algorithm requires sorting the signals at each time step and dispersion measure, an operation that is computationally tractable for existing array sizes. An alternative uses the various stations to form an ensemble estimate of the conditional density function (CDF) evaluated at each time step. Both methods outperform standard detection strategies on a test sequence of VLBA data, and both are efficient enough for deployment in real-time, online transient detection applications.

Thompson, David R.↗

DOE OSTI2025

Raman Digital Twin of Monolayer Janus Transition Metal Dichalcogenides

Monolayer transition metal dichalcogenides (TMDs) are a key class of two-dimensional (2D) materials with broad technological potential. Their Janus counterparts exhibit unique properties due to broken out-of-plane symmetry and further enrich the functionalities of TMDs. However, experimental synthesis and identification of Janus TMDs remain challenging. It is thus highly desirable to have a rapid, simple, and in situ characterization technique to monitor, in real time, the conversion process from the parent to Janus structure. Raman spectroscopy stands out for such a task as it is a powerful, nondestructive, and very commonly used tool to characterize 2D materials both in situ and ex situ. To realize the full potential of Raman spectroscopy on rapid characterization of Janus TMDs, we present a computational “Raman digital twin” library for various monolayer Janus TMDs in both 2H and Td phases. We focus on group-6 TMDs: MoS 2 , WS 2 , MoSe 2 , WSe 2 , MoTe 2 , WTe 2 and their Janus variants: MoSSe, MoSTe, MoSeTe, WSSe, WSTe, and WSeTe. Using first-principles density functional theory (DFT), we calculate their vibrational properties and predict distinct Raman fingerprints. These phonon and Raman signatures reflect each material’s structural symmetry and atomic composition, enabling clear identification via Raman spectroscopy. Our theoretical work supports experimental efforts by providing benchmarks for material identification, structural analysis, and quality control. In conclusion, the computational library expedites the discovery and development of Janus 2D materials, facilitating tighter integration between theoretical predictions and experimental validation.

Chalcogenides↗

DOE OSTI2021

Laplacian Smoothing Stochastic Gradient Markov Chain Monte Carlo

As an important Markov chain Monte Carlo (MCMC) method, the stochastic gradient Langevin dynamics (SGLD) algorithm has achieved great success in Bayesian learning and posterior sampling. Furthermore, SGLD typically suffers from a slow convergence rate due to its large variance caused by the stochastic gradient. In order to alleviate these drawbacks, we leverage the recently developed Laplacian smoothing technique and propose a Laplacian smoothing stochastic gradient Langevin dynamics (LS-SGLD) algorithm. We prove that for sampling from both log-concave and non-log-concave densities, LS-SGLD achieves strictly smaller discretization error in 2-Wasserstein distance, although its mixing rate can be slightly slower. Experiments on both synthetic and real datasets verify our theoretical results and demonstrate the superior performance of LS-SGLD on different machine learning tasks including posterior sampling, Bayesian logistic regression, and training Bayesian convolutional neural networks.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

A comparison of ISCCP and FIRE satellite cloud parameters

One of the goals of the First ISCCP Regional Experiment (FIRE) is the quantification of the uncertainties in the cloud parameter products derived by the International Satellite Cloud Climatology Project (ISCCP). This validation effort has many facets including sensitivity analyses and comparisons to similar data or theoretical results with known accuracies. The FIRE provides cloud-truth data at particular points or along particular lines from surface and aircraft measurement systems. Relating these data to the larger, area-averaged ISCCP results requires intermediate steps using higher resolution satellite data analyses. Errors in the cloud products derived with a particular method can be determined by performing analyses of high resolution satellite data over the area surrounding the point or line measurement. This same analysis technique may then be used to derive cloud parameters over a larger area containing similar cloud fields. It is assumed that the uncertainties found for the small scale analyses are the same for the large scale so that the method has been calibrated for the particular cloud type; i.e., its accuracy is known. Differences between the large scale results using the ISCCP technique and the calibrated method can be computed and used to determine if any significant biases or rms errors occur in the ISCCP results. Selected ISCCP results are compared to cloud parameters derived using the hybrid bispectral threshold method over the FIRE IFO and extended observation areas.

Gibson, Gary G.↗