Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

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,279 records · Page 71

Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods

Quantum algorithms for simulating electronic ground states are slower than popular classical mean-field algorithms such as Hartree–Fock and density functional theory but offer higher accuracy. Accordingly, quantum computers have been predominantly regarded as competitors to only the most accurate and costly classical methods for treating electron correlation. However, here we tighten bounds showing that certain first-quantized quantum algorithms enable exact time evolution of electronic systems with exponentially less space and polynomially fewer operations in basis set size than conventional real-time time-dependent Hartree–Fock and density functional theory. Although the need to sample observables in the quantum algorithm reduces the speedup, we show that one can estimate all elements of the k-particle reduced density matrix with a number of samples scaling only polylogarithmically in basis set size. We also introduce a more efficient quantum algorithm for first-quantized mean-field state preparation that is likely cheaper than the cost of time evolution. We conclude that quantum speedup is most pronounced for finite-temperature simulations and suggest several practically important electron dynamics problems with potential quantum advantage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Geminga neutron star: Evidence for nucleon superfluidity at very high density

A comparison of the recent age and temperature estimates of the Geminga neutron star with cooling models is presented. This star is already in the photon cooling era and it is shown that its temperature can be understood within both the slow and fast neutrino emission scenarios and consequently will not allow discrimination between these two scenarios. However in both cases agreement of the theoretical cooling curves with the observed temperature depends crucially on the presence of nucleon pairing in most, if not all, of the core.

Page, Dany↗

Analysis of Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed ~ 1 ppm random errors for 8-10 km altitudes and ~30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2lidar on the NASA DC-8 and added an 02lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected ~linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. We demonstrated measurements over the California Central Valley, to stratus cloud tops over the Pacific Ocean, over mountain regions with snow, and over several areas with broken clouds. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption and Range During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar only on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear C02 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed in detail and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed -1 ppm random errors for 8-10 km altitudes and -30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2 lidar on the NASA DC-8 and added an O2 lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected -linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. The seven flights in the 2011 Ascends campaign were flown over a wide variety of surface and cloud conditions in the US, which produced a wide variety of lidar signal conditions. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

Offshore lightning more active and intense in U.S. East Coast wind areas than in North Sea

Fast winds over open areas off the U.S. East Coast resulted in many offshore areas being opened for wind power deployment. Yet lightning is active in this same general area. Lightning damage has been noted as the leading cause of unplanned downtime of wind farms. Here we use lightning observations from 2020–2022 to quantify lightning activity—U.S. East Coast Wind Lease and Wind Planning areas received about 14-times the lightning strokes and about 18-times the energy transfer compared to the operational wind farms in the North Sea. A strong north-south U.S. gradient of lightning exists : low off the coast of Maine, growing higher off Chesapeake Bay, and high south of Virginia. Very high lightning activity also occurs about 250 km from the U.S. coastline, which is 50–100 km outside currently designated Wind Areas but likely relevant to planning floating wind farms further from shore. International lightning protection standards, such as IEC 614200-24, use lightning stroke density and turbine top-height to estimate damage to wind turbines, but it remains unclear if such standards could protect U.S. offshore wind turbines in active lightning areas where the turbine top-height can be 100 m taller than the next tallest object for tens to hundreds of kilometers. This analysis suggests U.S. East Coast offshore wind farms will be subjected to much higher lightning stroke and lightning energy densities than the North Sea counterparts. Preparing the turbine designs in anticipation of these more active U.S. lightning risks will help foster wind’s growing contribution to a reliable and resilient low-carbon energy future.

lightning↗

Prediction of Solar Irradiance and Photovoltaic Solar Energy Product Based on Cloud Coverage Estimation Using Machine Learning Methods

Cloud cover estimation from images taken by sky-facing cameras can be an important input for analyzing current weather conditions and estimating photovoltaic power generation. The constant change in position, shape, and density of clouds, however, makes the development of a robust computational method for cloud cover estimation challenging. Accurately determining the edge of clouds and hence the separation between clouds and clear sky is difficult and often impossible. Toward determining cloud cover for estimating photovoltaic output, we propose using machine learning methods for cloud segmentation. We compare several methods including a classical regression model, deep learning methods, and boosting methods that combine results from the other machine learning models. To train each of the machine learning models with various sky conditions, we supplemented the existing Singapore whole sky imaging segmentation database with hazy and overcast images collected by a camera-equipped Waggle sensor node. We found that the U-Net architecture, one of the deep neural networks we utilized, segmented cloud pixels most accurately. However, the accuracy of segmenting cloud pixels did not guarantee high accuracy of estimating solar irradiance. We confirmed that the cloud cover ratio is directly related to solar irradiance. Additionally, we confirmed that solar irradiance and solar power output are closely related; hence, by predicting solar irradiance, we can estimate solar power output. This study demonstrates that sky-facing cameras with machine learning methods can be used to estimate solar power output. This ground-based approach provides an inexpensive way to understand solar irradiance and estimate production from photovoltaic solar facilities.

14 SOLAR ENERGY↗

Vertical neutral wind in the equatorial F-region deduced from electric field and ion density measurements

Direct current (DC) electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 Universal Time (UT) on 25 October 1988, can be explained if a downward neutral wind of 15-30 m/s exists. Simultaneously, the F-region plasma was moving downward at a speed of 30-50 m/s. These events appear in the local time sector of 23.00-23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25-50 m/s are possible at time near midnight in the equatorial F-region.

Laakso, Harri↗

The mass and density of the dwarf planet (225088) 2007 OR10

The satellite of (225088) 2007 OR10 was discovered on archival Hubble Space Telescope images and along with new observations with the WFC3 camera in late 2017 we have been able to determine the orbit. The orbit's notable eccentricity, e ≈ 0.3, may be a consequence of an intrinsically eccentric orbit and slow tidal evolution, but may also be caused by the Kozai mechanism. Dynamical considerations also suggest that the moon is small, D(eff) < 100 km. Based on the newly determined system mass of 1.75 ·10^21 kg, 2007 OR10 is the fifth most massive dwarf planet after Eris, Pluto, Haumea and Makemake. The newly determined orbit has also been considered as an additional option in our radiometric analysis, provided that the moon orbits in the equatorial plane of the primary. Assuming a spherical shape for the primary this approach provides a size of 1230 ± 50 km, with a slight dependence on the satellite orbit orientation and primary rotation rate chosen, and a bulk density of 1.75 ± 0.07 g/cu.cm for the primary. A previous size estimate that assumed an equator-on configuration (1535(−225,+75) km) would provide a density of 0.92(−0.14,+0.46) g/cu.cm, unexpectedly low for a 1000 km-sized dwarf planet.

Csaba Kiss↗

Estimation of proportions in mixed pixels through their region characterization

A region of mixed pixels can be characterized through the probability density function of proportions of classes in the pixels. Using information from the spectral vectors of a given set of pixels from the mixed pixel region, expressions are developed for obtaining the maximum likelihood estimates of the parameters of probability density functions of proportions. The proportions of classes in the mixed pixels can then be estimated. If the mixed pixels contain objects of two classes, the computation can be reduced by transforming the spectral vectors using a transformation matrix that simultaneously diagonalizes the covariance matrices of the two classes. If the proportions of the classes of a set of mixed pixels from the region are given, then expressions are developed for obtaining the estmates of the parameters of the probability density function of the proportions of mixed pixels. Development of these expressions is based on the criterion of the minimum sum of squares of errors. Experimental results from the processing of remotely sensed agricultural multispectral imagery data are presented.

Chittineni, C. B.↗

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗

Extended adiabatic blast waves and a model of the soft X-ray background

An analytical approximation is generated which follows the development of an adiabatic spherical blast wave in a homogeneous ambient medium of finite pressure. An analytical approximation is also presented for the electron temperature distribution resulting from coulomb collisional heating. The dynamical, thermal, ionization, and spectral structures are calculated for blast waves of energy E sub 0 = 5 x 10 to the 50th power ergs in a hot low-density interstellar environment. A formula is presented for estimating the luminosity evolution of such explosions. The B and C bands of the soft X-ray background, it is shown, are reproduced by such a model explosion if the ambient density is about .000004 cm, the blast radius is roughly 100 pc, and the solar system is located inside the shocked region. Evolution in a pre-existing cavity with a strong density gradient may, it is suggested, remove both the M band and OVI discrepancies.

Cox, D. P.↗

Using partially labeled data for normal mixture identification with application to class definition

The problem of estimating the parameters of a normal mixture density when, in addition to the unlabeled samples, sets of partially labeled samples are available is addressed. The density of the multidimensional feature space is modeled with a normal mixture. It is assumed that the set of components of the mixture can be partitioned into several classes and that training samples are available from each class. Since for any training sample the class of origin is known but the exact component of origin within the corresponding class is unknown, the training samples as considered to be partially labeled. The EM iterative equations are derived for estimating the parameters of the normal mixture in the presence of partially labeled samples. These equations can be used to combine the supervised and nonsupervised learning processes.

Shahshahani, Behzad M.↗

Estimation of plasma properties using an extended Kalman filter with plasma global models

A physically-constrained extended Kalman filter (EKF) is applied to various zero-dimensional global models for the estimation of plasma properties using time-dependent experimental data such as the plasma density or ion flux. The capability of the EKF is demonstrated to estimate unknown system states simultaneously, such as reaction rate coefficients and the absorbed electron input power, which can be difficult, if not impossible, to measure experimentally. Global models accounting for pure argon reactions and argon-oxygen reactions are used in this work to demonstrate the ability of the filter to estimate dynamic and complex systems. Furthermore, the results obtained from the EKF plasma global model illustrate that model-data fusion techniques can be used to estimate plasma properties and processes for time-varying systems, such as pulsed discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of matter density profiles on neutrino oscillations for T2HK and T2HKK

This paper explores the effects of changes in matter density profiles on neutrino oscillation probabilities, and whether these could potentially be seen by the future Hyper-Kamiokande long-baseline oscillation experiment (T2HK). The analysis is extended to include the possibility of having an additional detector in Korea (T2HKK). In both cases, we find that these effects will be immeasurable, as the magnitudes of the changes in the oscillation probabilities induced in all density profile scenarios considered here remain smaller than the estimated experimental sensitivity to the oscillation probabilities of each experiment, for both appearance and disappearance channels. Therefore, we conclude that using a constant density profile is sufficient for both the T2HK and T2HKK experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spectral density reconstruction with Chebyshev polynomials

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. Here, in the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

ICESat Observations of Seasonal and Interannual Variations of Sea-Ice Freeboard and Estimated Thickness in the Weddell Sea, Antarctica (2003-2009)

Sea-ice freeboard heights for 17 ICESat campaign periods from 2003 to 2009 are derived from ICESat data. Freeboard is combined with snow depth from Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) data and nominal densities of snow, water and sea ice, to estimate sea-ice thickness. Sea-ice freeboard and thickness distributions show clear seasonal variations that reflect the yearly cycle of growth and decay of the Weddell Sea (Antarctica) pack ice. During October-November, sea ice grows to its seasonal maximum both in area and thickness; the mean freeboards are 0.33-0.41 m and the mean thicknesses are 2.10-2.59 m. During February-March, thinner sea ice melts away and the sea-ice pack is mainly distributed in the west Weddell Sea; the mean freeboards are 0.35-0.46 m and the mean thicknesses are 1.48-1.94 m. During May-June, the mean freeboards and thicknesses are 0.26-0.29 m and 1.32-1.37 m, respectively. The 6 year trends in sea-ice extent and volume are (0.023+/-0.051) x 10(exp 6)sq km/a (0.45%/a) and (0.007+/-1.0.092) x 10(exp 3)cu km/a (0.08%/a); however, the large standard deviations indicate that these positive trends are not statistically significant.

Yi, Donghui↗

A Model for Properties of Basin Ejecta Deposits and Secondary Crater Densities

Ejecta scaling relationships, ballistic sedimentation, and ejecta fragment size distributions are used to give estimates of thicknesses of basin ejecta deposits, proportions of primary ejecta, and secondary crater diameters and surface densities. Additional information is contained in the original extended abstract.

Haskin, L. A.↗