Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Relating Radar Backscatter to Boreal Forest Canopy Parameters

During the first intensive field campaign for the Boreal Ecosystem Atmospheric Study (BOREAS) in August of 1993, AIRSAR data were acquired over the entire BOREAS study area. On of the objectives of the AIRSAR deployment was to examine the sensitivity of the radar signal to vegetation type and biomass distribution in the boreal forest and to develop algorithms for inferring vegetation parameters. During the experiment a set of ground measurements were also made to support the AIRSAR data analysis. The dominant stands in the study area consist of black spruce, young jack pine, old jack pine, and aspen. These stands represent a wide range of biomass and canopy architectural variations which can be distinguished in SAR images.

radar backscatter↗
NASA NTRS2012

The Physics of Imaging with Remote Sensors : Photon State Space & Radiative Transfer

Standard (mono-pixel/steady-source) retrieval methodology is reaching its fundamental limit with access to multi-angle/multi-spectral photo- polarimetry. Next... Two emerging new classes of retrieval algorithm worth nurturing: multi-pixel time-domain Wave-radiometry transition regimes, and more... Cross-fertilization with bio-medical imaging. Physics-based remote sensing: - What is "photon state space?" - What is "radiative transfer?" - Is "the end" in sight? Two wide-open frontiers! center dot Examples (with variations.

clouds↗
NASA NTRS2018

NASA's Black Marble Nighttime Lights Product Suite

NASA's Black Marble nighttime lights product suite (VNP46) is available at 500 meters resolution since January 2012 with data from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band (DNB) onboard the Suomi National Polar-orbiting Platform (SNPP). The retrieval algorithm, developed and implemented for routine global processing at NASA's Land Science Investigator-led Processing System (SIPS), utilizes all high-quality, cloud-free, atmospheric-, terrain-, vegetation-, snow-, lunar-, and stray light-corrected radiances to estimate daily nighttime lights (NTL) and other intrinsic surface optical properties. Key algorithm enhancements include: (1) lunar irradiance modeling to resolve non-linear changes in phase and libration; (2) vector radiative transfer and lunar bidirectional surface anisotropic reflectance modeling to correct for atmospheric and BRDF (Bidirectional Reflectance Distribution Function) effects; (3) geometric-optical and canopy radiative transfer modeling to account for seasonal variations in NTL; and (4) temporal gap-filling to reduce persistent data gaps. Extensive benchmark tests at representative spatial and temporal scales were conducted on the VNP46 time series record to characterize the uncertainties stemming from upstream data sources. Initial validation results are presented together with example case studies illustrating the scientific utility of the products. This includes an evaluation of temporal patterns of NTL dynamics associated with urbanization, socioeconomic variability, cultural characteristics, and displaced populations affected by conflict. Current and planned activities under the Group on Earth Observations (GEO) Human Planet Initiative are aimed at evaluating the products at different geographic locations and time periods representing the full range of retrieval conditions.

Suomi-NPP↗
NASA NTRS2014

Fine Particulate Matter Predictions Using High Resolution Aerosol Optical Depth (AOD) Retrievals

To date, spatial-temporal patterns of particulate matter (PM) within urban areas have primarily been examined using models. On the other hand, satellites extend spatial coverage but their spatial resolution is too coarse. In order to address this issue, here we report on spatial variability in PM levels derived from high 1 km resolution AOD product of Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm developed for MODIS satellite. We apply day-specific calibrations of AOD data to predict PM(sub 2.5) concentrations within the New England area of the United States. To improve the accuracy of our model, land use and meteorological variables were incorporated. We used inverse probability weighting (IPW) to account for nonrandom missingness of AOD and nested regions within days to capture spatial variation. With this approach we can control for the inherent day-to-day variability in the AOD-PM(sub 2.5) relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance among others. Out-of-sample "ten-fold" cross-validation was used to quantify the accuracy of model predictions. Our results show that the model-predicted PM(sub 2.5) mass concentrations are highly correlated with the actual observations, with out-of- sample R(sub 2) of 0.89. Furthermore, our study shows that the model captures the pollution levels along highways and many urban locations thereby extending our ability to investigate the spatial patterns of urban air quality, such as examining exposures in areas with high traffic. Our results also show high accuracy within the cities of Boston and New Haven thereby indicating that MAIAC data can be used to examine intra-urban exposure contrasts in PM(sub 2.5) levels.

aerosol optical depth↗
DOE OSTI2022

U-splines: Splines over unstructured meshes

U-splines are a novel approach to the construction of a spline basis for representing smooth objects in Computer-Aided Design (CAD) and Computer-Aided Engineering (CAE). A spline is a piecewise-defined function that satisfies continuity constraints between adjacent cells in a mesh. U-splines differ from existing spline constructions, such as Non-Uniform Rational B-splines (NURBS), subdivision surfaces, T-splines, and hierarchical B-splines, in that they can accommodate local variation in cell size, polynomial degree, and smoothness simultaneously over more varied mesh configurations. Mixed cell types (e.g., triangle and quadrilateral cells in the same mesh) and T-junctions are also supported, although the continuity of interfaces with triangle and tetrahedral cells is limited in the present work. The U-spline algorithm introduces a new technique for using local null space solutions to construct basis functions for the global spline null space problem. The U-spline construction is presented for curves, surfaces, and volumes with higher dimensional generalizations possible. Lastly, a set of requirements are given to ensure that the U-spline basis is positive, forms a partition of unity, is complete, and is locally linearly independent.

42 ENGINEERING↗
DOE OSTI2020

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Resource-Optimized Fermionic Local-Hamiltonian Simulation on a Quantum Computer for Quantum Chemistry

The ability to simulate a fermionic system on a quantum computer is expected to revolutionize chemical engineering, materials design, nuclear physics, to name a few. Thus, optimizing the simulation circuits is of significance in harnessing the power of quantum computers. Here, we address this problem in two aspects. In the fault-tolerant regime, we optimize the R z and T gate counts along with the ancilla qubit counts required, assuming the use of a product-formula algorithm for implementation. We obtain a savings ratio of two in the gate counts and a savings ratio of eleven in the number of ancilla qubits required over the state of the art. In the pre-fault tolerant regime, we optimize the two-qubit gate counts, assuming the use of the variational quantum eigensolver (VQE) approach. Specific to the latter, we present a framework that enables bootstrapping the VQE progression towards the convergence of the ground-state energy of the fermionic system. This framework, based on perturbation theory, is capable of improving the energy estimate at each cycle of the VQE progression, by about a factor of three closer to the known ground-state energy compared to the standard VQE approach in the test-bed, classically-accessible system of the water molecule. The improved energy estimate in turn results in a commensurate level of savings of quantum resources, such as the number of qubits and quantum gates, required to be within a pre-specified tolerance from the known ground-state energy. We also explore a suite of generalized transformations of fermion to qubit operators and show that resource-requirement savings of up to more than 20 % , in small instances, is possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

Numerical comparison of Kalman filter algorithms - Orbit determination case study

Numerical characteristics of various Kalman filter algorithms are illustrated with a realistic orbit determination study. The case study of this paper highlights the numerical deficiencies of the conventional and stabilized Kalman algorithms. Computational errors associated with these algorithms are found to be so large as to obscure important mismodeling effects and thus cause misleading estimates of filter accuracy. The positive result of this study is that the U-D covariance factorization algorithm has excellent numerical properties and is computationally efficient, having CPU costs that differ negligibly from the conventional Kalman costs. Accuracies of the U-D filter using single precision arithmetic consistently match the double precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity to variations in the a priori statistics.

Bierman, G. J.↗
NASA NTRS1976

Numerical comparison of discrete Kalman filter algorithms - Orbit determination case study

Numerical characteristics of various Kalman filter algorithms are illustrated with a realistic orbit determination study. The case study of this paper highlights the numerical deficiencies of the conventional and stabilized Kalman algorithms. Computational errors associated with these algorithms are found to be so large as to obscure important mismodeling effects and thus cause misleading estimates of filter accuracy. The positive result of this study is that the U-D covariance factorization algorithm has excellent numerical properties and is computationally efficient, having CPU costs that differ negligibly from the conventional Kalman costs. Accuracies of the U-D filter using single precision arithmetic consistently match the double precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity to variations in the a priori statistics.

Bierman, G. J.↗
NASA NTRS1983

A stochastic atmospheric model for remote sensing applications

There are many factors which reduce the accuracy of classification of objects in the satellite remote sensing of Earth's surface. One important factor is the variability in the scattering and absorptive properties of the atmospheric components such as particulates and the variable gases. For multispectral remote sensing of the Earth's surface in the visible and infrared parts of the spectrum the atmospheric particulates are a major source of variability in the received signal. It is difficult to design a sensor which will determine the unknown atmospheric components by remote sensing methods, at least to the accuracy needed for multispectral classification. The problem of spatial and temporal variations in the atmospheric quantities which can affect the measured radiances are examined. A method based upon the stochastic nature of the atmospheric components was developed, and, using actual data the statistical parameters needed for inclusion into a radiometric model was generated. Methods are then described for an improved correction of radiances. These algorithms will then result in a more accurate and consistent classification procedure.

Turner, R. E.↗
NASA NTRS1984

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗
NASA NTRS1993

Convective and stratiform rain: Multichannel microwave sensing over oceans

Measurements made by the Special Sensor Microwave/Imager (SSM/I) radiometer over the oceans, at 19, 37, and 85 GHz in dual polarization, are used to develop a model to classify rain into light-stratiform, moderately convective, and heavy convective types in the mesoscale convective systems (MCS). It is observed that the bulk of the 19- and 37-GHz data are linearly correlated with respect to one another, and generally increase together in brightness as the mean rain rate in the field of view (FOV) of the radiometer increases. However, a significant fraction of the data from these channels departs from this linear relationship, reflecting the nonuniform rain that is convective vs. the relatively light stratiform rain. It is inferred from the SSM/I data, in a MCS, when the slope dT sub 3/dT sub 19 is greater than unity there are optically thin clouds which produce light uniform rain. On the other hand, when dT sub 3/dT sub 19 is close to unity, the rain cells have an open structure and correspond to the convective type of rain. The openings between the cells are apparently a result of the downdrafts and/or entrainment. Relatively low values of 85-GHz brightness temperatures that are present when dT sub 37/dT sub 19 is close to unity support these views and, in addition, leads us to conclude that when the convection is heavy this brightness temperature decreases due to scattering by hydrometeors. On the basis of this explanation of the SSM/I data, an empirical rain retrieval algorithm is developed. Radar backscatter observations over the Atlantic Ocean next to Florida are used to demonstrate the applicability of this method. Three monthly mean maps of rainfall over the oceans from 50 degrees N to 50 degrees S, are presented to illustrate the ability of this method to sense seasonal and interannual variations of rain.

Prabhakara, C.↗
NASA NTRS2001

Precipitation Processes Derived from TRMM Satellite Data, Cloud Resolving Model and Field Campaigns

Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent-heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid. and solid water. Present large-scale weather and climate models can simulate cloud latent heat release only crudely thus reducing their confidence in predictions on both global and regional scales. In this paper, NASA Tropical Rainfall Measuring (TRMM) precipitation radar (PR) derived rainfall information and the Goddard Convective and Stratiform Heating (CSH) algorithm used to estimate the four-dimensional structure of global monthly latent heating and rainfall profiles over the global tropics from December 1997 to October 2000. Rainfall latent heating and radar reflectively structure between ENSO (1997-1998 winter) and non-ENSO (1998-1999 winter) periods are examined and compared. The seasonal variation of heating over various geographic locations (i.e. Indian ocean vs west Pacific; Africa vs S. America) are also analyzed. In addition, the relationship between rainfall latent heating maximum heating level), radar reflectively and SST are examined.

Tao, W.-K.↗
NASA NTRS1973

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗
DOE OSTI2026

Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: a benchmark Quantum Monte Carlo study

The accurate description of non-covalent interactions is critical for understanding the structure, dynamics, and eventual function of biomolecules. The adenine dimer serves as a benchmark system for computational methods due to its role in nucleic acid structures and its rich conformational landscape. In this study, we employ benchmark diffusion quantum Monte Carlo (DMC) methods to investigate the relative energies and role of electron correlation on a set of adenine dimer conformations generated via a search of the potential energy landscape using the global optimizer algorithm. Relative DMC energies are compared against a wide range of density functional theory (DFT) approximation results. We find that although most of the DFT functionals perform well for low-energy structures, their accuracy varies significantly for higher-energy conformations, including stacked and T-shaped structures. A large fraction of the variation is due to the treatment of the van der Waals interaction. BLYP, B3LYP, and PBE0 significantly improve with added D4 dispersion, while the recent r2SCAN-D4 and ωB97M-V functionals show the least scatter and closest agreement with the DMC. These findings highlight the delicate nature of these interactions in biomolecular systems and provide guidance for simulations of their structure and dynamics and for the development of machine learned interatomic potentials.

Washburn, Laurel [ORNL] (ORCID:0000000324179335)↗
NASA NTRS1975

A Monte Carlo analysis of the effects of instrumentation errors on aircraft parameter identification

An output error estimation algorithm was used to evaluate the effects of both static and dynamic instrumentation errors on the estimation of aircraft stability and control parameters. A Monte Carlo analysis, using simulated cruise flight data, was performed for a high performance military aircraft, a large commercial transport, and a small general-aviation aircraft. The effects of variations in the information content of the flight data, resulting from two different choices of control input maneuvers, were also determined. The results indicate that unmodeled instrumentation errors can cause inaccuracies in the estimated parameters which are comparable to their nominal values. Control input errors and angular accelerometer lags were found to be most significant of the instrumentation errors evaluated, and the perturbations they produce are much larger than those arising from the combined effects of static errors and white noise in the output response measurements.

Bryant, W. H.↗
NASA NTRS1982

Numerical solution of transonic wing flow fields

A fast, fully implicit, approximate factorization algorithm designed to solve the conservative full-potential equation is used to compute lifting-wing flow fields with embedded transonic flow. The computer program (TWING - transonic wing analysis) uses an elliptic-solver numerical grid-generation routine and is capable of analyzing arbitrary wing geometries with sweep, twist, taper, and section variation. Computed flowfield results for several wing geometries are presented; they demonstrate substantial improvement in convergence speed relative to the FL028 computer code.

Holst, T. L.↗
NASA NTRS1993

Estimation of solar backscatter ultraviolet albedo using ground-based Umkehr measurements

A retrieval method was developed to estimate the solar backscatter ultraviolet (SBUV) satellite albedo for the ozone profiler wavelengths using ground-based ultraviolet measurements. For the present investigation the Umkehr was used as the ground-based ultraviolet measurement. Simulated SBUV data and Umkehr data theoretically computed from a priori ozone profiles observed by the SAGE II satellite were used to develop the retrieval algorithm and to test its capability. The test indicated that albedos for the SBUV ozone profiler wavelengths should allow estimates to a precision of +/- 5 percent or better, depending on the accuracy of the ultraviolet measurement. Retrievals using actual Umkehr observations were also performed to provide a preliminary look at the magnitude and annual variation of retrieved albedos. A case study was performed, comparing retrieved albedos with SBUV-measured albedos. The SBUV albedo change was seen to be approximately twice as large as the albedo changes estimated by the Umkehr method. Results of the investigation suggest that the method of estimation may be useful for determining the drift rate of the SBUV calibration.

Deluisi, John J.↗