Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7

Comparison of removal and spatial mark‐resight models for estimating wild pig density

Density estimation is critical to effectively manage invasive species and elucidate areas of highest concern. For wild pigs (Sus scrofa), the ability to estimate density is complicated because of their variable home range sizes and social structure. Common methods for estimating density (e.g., mark-recapture) may be unsuitable in management applications because additional data needs to be collected before and after management. Removal models offer a suitable alternative to estimate density changes following management and can be applied broadly across areas where management of wild pigs is ongoing. We collected wild pig removal and camera trap data from 25 private properties ranging in size from approximately 0.5 km 2 to 95 km 2 across 3 ecoregions in South Carolina, USA, from 2020–2023. We compared factors affecting consistency and precision of property-level density estimates between removal and spatial mark-resight (SMR) models. In general, excluding 1 large outlier, density estimates from removal models were between 0.60 and 15.85 wild pigs/km 2 (median = 5.34) with a median coefficient of variation (CV) of 0.76 and 95% confidence intervals for the CV between 0.70 and 0.94. Similarly, excluding 1 large outlier, density estimates from SMR were between 0.22 and 30.97 wild pigs/km 2 (median = 5.48) with a median CV of 0.39 and 95% confidence intervals for the CV between 0.38 and 1.20. We found the precision of removal models was affected primarily by the number of wild pigs dispatched in the removal period (3 months) and the ecoregion in which they were removed. None of the covariates, including the number of recaptures (a corresponding measure of sample size), influenced precision of the SMR models, although recaptures did influence the density estimates. At the individual property level, density estimates from our 2 estimators were dissimilar from each other in approximately 80% of instances, although none of the covariates we examined influenced dissimilarity. Our results provide unique insight into how sample size affects density estimates using 2 common methods and into novel SMR models that incorporate both marked and unmarked detections. In addition, the density estimates in this study can be used as a reference for wild pig densities in common land cover types throughout the southeastern United States.

60 APPLIED LIFE SCIENCES↗

Can we estimate plasma density in ICP driver through electrical parameters in RF circuit?

To avoid regular maintenance, invasive plasma diagnostics with probes are not included in the inductively coupled plasma (ICP) based ITER Neutral Beam (NB) source design. Even non-invasive probes like optical emission spectroscopic diagnostics are also not included in the present ITER NB design due to overall system design and interface issues. As a result, negative ion beam current through the extraction system in the ITER NB negative ion source is the only measurement which indicates plasma condition inside the ion source. However, beam current not only depends on the plasma condition near the extraction region but also on the perveance condition of the ion extractor system and negative ion stripping. Nevertheless, inductively coupled plasma production region (RF driver region) is placed at distance (∼ 30cm) from the extraction region. Due to that, some uncertainties are expected to be involved if one tries to link beam current with plasma properties inside the RF driver. Plasma characterization in source RF driver region is utmost necessary to maintain the optimum condition for source operation. In this paper, a method of plasma density estimation is described, based on density dependent plasma load calculation.

43 PARTICLE ACCELERATORS↗

Effects of scale of movement, detection probability, and true population density on common methods of estimating population density

Knowledge of population density is necessary for effective management and conservation of wildlife, yet rarely are estimators compared in their robustness to effects of ecological and observational processes, which can greatly influence accuracy and precision of density estimates. For this study, we simulate biological and observational processes using empirical data to assess effects of animal scale of movement, true population density, and probability of detection on common density estimators. We also apply common data collection and analytical techniques in the field and evaluate their ability to estimate density of a globally widespread species. We find that animal scale of movement had the greatest impact on accuracy of estimators, although all estimators suffered reduced performance when detection probability was low, and we provide recommendations as to when each field and analytical technique is most appropriately employed. The large influence of scale of movement on estimator accuracy emphasizes the importance of effective post-hoc calculation of area sampled or use of methods that implicitly account for spatial variation. In particular, scale of movement impacted estimators substantially, such that area covered and spacing of detectors (e.g. cameras, traps, etc.) must reflect movement characteristics of the focal species to reduce bias in estimates of movement and thus density.

59 BASIC BIOLOGICAL SCIENCES↗

Raccoon densities across four land cover types in the southeastern United States

Raccoons (Procyon lotor) are the primary reservoir for rabies virus in eastern North America. Management of rabies in raccoons is achieved primarily with the use of oral rabies vaccination (ORV) and effective ORV bait densities are determined in part by the densities of raccoons. Decisions regarding ORV bait densities, however, are limited by an incomplete understanding of raccoon densities across the spectrum of landscapes they occupy. We carried out a mark-recapture study of raccoons on the Savannah River Site in South Carolina, USA, from 2017–2019, to develop sex- and landscape-specific raccoon density estimates across 4 rural land cover types in the southeastern United States: bottomland hardwood, riparian forest, isolated wetland, and upland pine (Pinus spp.). We captured 404 unique raccoons 773 times over the 3-year trapping period. Estimated densities were 5.44 ± 0.37 (SE) animals/km 2 in bottomland hardwood forest, 2.62 ± 0.32 animals/km 2 in riparian forest, 2.19 ± 0.29 animals/km 2 in isolated wetlands, and 2.14 ± 0.23 animals/km 2 in upland pine. Densities were significantly higher in bottomland hardwood than all other land cover types, whereas densities among the remaining cover types were similar. These patterns are likely the result of landscape fragmentation and configuration, with riparian forests typically embedded in a matrix of less suitable cover types, leading to low densities despite presumably high resource availability. There were higher densities of males than females in every cover type except upland pine, where the sex ratio was balanced. Densities on our site were low compared to other rural areas, which likely results from the lack of human influence in terms of agriculture or development. The financial cost of baiting for ORV distribution may be reduced by considering the comparatively low densities of raccoons in these rural landscapes in the southeastern United States.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of estimators of probability density functions

Four research projects are summarized which include: (1) the generation of random numbers on the IBM 360/44, (2) statistical tests used to check out random number generators, (3) Specht density estimators, and (4) use of estimators of probability density functions in analyzing large amounts of data.

Speed, F. M.↗

SPITZER OBSERVATIONS OF SUPERNOVA REMNANTS. II. PHYSICAL CONDITIONS AND COMPARISON WITH HH7 AND HH54

We have studied the shock-excited molecular regions associated with four supernova remnants (SNRs)-IC443C, W28, W44, and 3C391-and two Herbig-Haro objects-HH7 and HH54-using Spitzer's Infrared Spectrograph (IRS). The physical conditions within the observed areas (roughly {approx}1' x 1' in size) are inferred from spectroscopic data obtained from IRS and from the Short and Long Wavelength Spectrometers on board the Infrared Space Observatory (ISO), together with photometric data from Spitzer's Infrared Array Camera. Adopting a power-law distribution for the gas temperature in the observed region, with the mass of gas at temperature T to T+dT assumed proportional to T{sup -b} dT, the H{sub 2} S(0) to S(7) spectral line maps obtained with IRS were used to constrain the gas density, yielding estimated densities n(H{sub 2}) in the range of {approx} (2-4)x10{sup 3} cm{sup -3}. The excitation of H{sub 2} S(9) to S(12) and high-J CO pure rotational lines, however, require environments several times denser. The inconsistency among the best-fit densities estimated from different species can be explained by density fluctuations within the observed regions. The best-fit power-law index b is smaller than the value 3.8 predicted for a paraboloidal C-type bow shock, suggesting that the shock front has a 'flatter' shape than that of a paraboloid. The best-fit parameters for SNRs and Herbig-Haro objects do not differ significantly between the two classes of sources, except that for the SNRs the ortho-to-para ratio (OPR) of hot gas (T> 1000 K) is close to the local thermal equilibrium value 3, while for HH7 and HH54 even the hottest gas exhibits an OPR smaller than 3; we interpret this difference as resulting from environmental differences between these classes of sources, molecular material near SNRs being subject to stronger photodissociation that results in faster para-to-ortho conversion. Finally, we mapped the physical parameters within the regions observed with IRS and found that the mid-lying H{sub 2} emissions-S(3) to S(5)-tend to trace the hot component of the gas, while the intensities of S(6) and S(7) are more sensitive to the density of the gas compared to S(3) to S(5).

79 ASTRONOMY AND ASTROPHYSICS↗

Method for estimating the density of high-level nuclear waste glass

A database of over 1100 silicate glass compositions and densities was compiled and used to evaluate the efficacy of an algorithm for estimating the density of silicate glass compositions. We sought to develop a parsimonious algorithm based on the additivity of partial molar volumes of individual oxide components weighted by their mole fraction in a glass composition. Bound molar volumes were used for oxides in which the density of the oxide bound in a glass matrix was previously determined. The bound molar volumes were known for oxides covering 97.5 mole percent of the database compositional space. The measured glass densities were plotted against the estimated glass densities and a linear regression yielded an R 2 adj. = 0.95 and a slope and intercept of approximately one and zero, respectively. This regression suggests that glass densities estimated by the algorithm, within analysis uncertainty, are equal to the measured densities of the glasses. In addition to the development of the density estimation, we corroborated many of the referenced bound molar volume data used in the parameterization of the estimation algorithm via linear regression of the individual partial molar volumes versus the inverse measured densities (specific volumes) of the glasses in the database.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-F-J-213: Feasibility Study of Using a Dual-Energy Cone Beam CT (DECBCT) in Proton Therapy Treatment Planning

Purpose: The aim of this study is to evaluate the feasibility of using a dual-energy CBCT (DECBCT) in proton therapy treatment planning to allow for accurate electron density estimation. Methods: For direct comparison, two scenarios were selected: a dual-energy fan-beam CT (high: 140 kVp, low: 80 kVp) and a DECBCT (high: 125 kVp, low: 80 kVp). A Gammex 467 tissue characterization phantom was used, including the rods of air, water, bone (B2–30% mineral), cortical bone (SB3), lung (LN-300), brain, liver and adipose. For the CBCT, Hounsfield Unit (HU) numbers were first obtained from the reconstructed images after a calibration was made based on water (=0) and air materials (=−1000). For each tissue surrogate, region-of-interest (ROI) analyses were made to derive high-energy and low-energy HU values (HUhigh and HUlow), which were subsequently used to estimate electron density based on the algorithm as previously described by Hunemohr N., et al. Parameters k1 and k2 are energy dependent and can be derived from calibration materials. Results: While for the dual-energy FBCT, the electron density is found be within +/−3% error relative to the values provided by the phantom vendor: −1.8% (water), 0.03% (lung), 1.1% (brain), −2.82% (adipose), −0.49% (liver) and −1.89% (cortical bones). While for the DECBCT, the estimation of electron density exhibits a relatively larger variation: −1.76% (water), −36.7% (lung), −1.92% (brain), −3.43% (adipose), 8.1% (liver) and 9.5% (cortical bones). Conclusion: For DECBCT, the accuracy of electron density estimation is inferior to that of a FBCT, especially for materials of either low-density (lung) or high density (cortical bone) compared to water. Such limitation arises from inaccurate HU number derivation in a CBCT. Advanced scatter-correction and HU calibration routines, as well as the deployment of photon counting CT detectors need be investigated to minimize the difference between FBCT and CBCT.

60 APPLIED LIFE SCIENCES↗

Physics-Based SNOWPACK Model Improves Representation of Near-Surface Antarctic Snow and Firn Density

Estimates of snow and firn density are required for satellite-altimetry-based retrievals of ice sheet mass balance that rely on volume-to-mass conversions. Therefore, biases and errors in presently used density models confound assessments of ice sheet mass balance and by extension ice sheet contribution to sea level rise. Despite this importance, most contemporary firn densification models rely on simplified semi-empirical methods, which are partially reflected by significant modeled density errors when compared to observations. In this study, we present a new drifting-snow compaction scheme that we have implemented into SNOWPACK, a physics-based land surface snow model. We show that our new scheme improves existing versions of SNOWPACK by increasing simulated near-surface (defined as the top 10 m) density to be more in line with observations (near-surface bias reduction from −44.9 to −5.4 kg m−3). Furthermore, we demonstrate high-quality simulation of near-surface Antarctic snow and firn density at 122 observed density profiles across the Antarctic ice sheet, as indicated by reduced model biases throughout most of the near-surface firn column when compared to two semi-empirical firn densification models (SNOWPACK mean bias=−9.7 kg m−3, IMAU-FDM mean bias=−32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the near surface where firn density is most variable due to accumulation and compaction variability driven by synoptic weather and seasonal climate variability. Additionally, the GSFC-FDM exhibits lower mean density bias from 7–10 m (SNOWPACK bias=−22.5 kg m−3, GSFC-FDM bias=10.6 kg m−3) and throughout the entire near surface at high-accumulation sites (SNOWPACK bias=−31.4 kg m−3, GSFC-FDM bias=−4.7 kg m−3). However, we found that the performance of SNOWPACK did not degrade when applied to sites that were not included in the calibration of semi-empirical models. This suggests that SNOWPACK may possibly better represent firn properties in locations without extensive observations and under future climate scenarios, when firn properties are expected to diverge from their present state.

Eric Keenan↗

A COMPARATIVE STUDY OF DENSITY FIELD ESTIMATION FOR GALAXIES: NEW INSIGHTS INTO THE EVOLUTION OF GALAXIES WITH ENVIRONMENT IN COSMOS OUT TO z ∼ 3

It is well-known that a galaxy’s environment has a fundamental influence in shaping its properties. We study the environmental effects on galaxy evolution, with an emphasis on the environment defined as the local number density of galaxies. The density field is estimated with different estimators (weighted adaptive kernel smoothing, 10th and 5th nearest neighbors, Voronoi and Delaunay tessellation) for a K{sub s} < 24 sample of ∼190,000 galaxies in the COSMOS field at 0.1 < z < 3.1. The performance of each estimator is evaluated with extensive simulations. We show that overall there is a good agreement between the estimated density fields using different methods over ∼2 dex in overdensity values. However, our simulations show that adaptive kernel and Voronoi tessellation outperform other methods. Using the Voronoi tessellation method, we assign surface densities to a mass complete sample of quiescent and star-forming galaxies out to z ∼ 3. We show that at a fixed stellar mass, the median color of quiescent galaxies does not depend on their host environment out to z ∼ 3. We find that the number and stellar mass density of massive (>10{sup 11} M{sub ⊙}) star-forming galaxies have not significantly changed since z ∼ 3, regardless of their environment. However, for massive quiescent systems at lower redshifts (z ≲ 1.3), we find a significant evolution in the number and stellar mass densities in denser environments compared to lower density regions. Our results suggest that the relation between stellar mass and local density is more fundamental than the color–density relation and that environment plays a significant role in quenching star-formation activity in galaxies at z ≲ 1.

79 ASTRONOMY AND ASTROPHYSICS↗

Control algorithms for aerobraking in the Martian atmosphere

The Analytic Predictor Corrector (APC) and Energy Controller (EC) atmospheric guidance concepts were adapted to control an interplanetary vehicle aerobraking in the Martian atmosphere. Changes are made to the APC to improve its robustness to density variations. These changes include adaptation of a new exit phase algorithm, an adaptive transition velocity to initiate the exit phase, refinement of the reference dynamic pressure calculation and two improved density estimation techniques. The modified controller with the hybrid density estimation technique is called the Mars Hybrid Predictor Corrector (MHPC), while the modified controller with a polynomial density estimator is called the Mars Predictor Corrector (MPC). A Lyapunov Steepest Descent Controller (LSDC) is adapted to control the vehicle. The LSDC lacked robustness, so a Lyapunov tracking exit phase algorithm is developed to guide the vehicle along a reference trajectory. This algorithm, when using the hybrid density estimation technique to define the reference path, is called the Lyapunov Hybrid Tracking Controller (LHTC). With the polynomial density estimator used to define the reference trajectory, the algorithm is called the Lyapunov Tracking Controller (LTC). These four new controllers are tested using a six degree of freedom computer simulation to evaluate their robustness. The MHPC, MPC, LHTC, and LTC show dramatic improvements in robustness over the APC and EC.

Ward, Donald T.↗

Picturing Data With Uncertainty

NASA is in the business of creating maps for scientific purposes to represent important biophysical or geophysical quantities over space and time. For example, maps of surface temperature over the globe tell scientists where and when the Earth is heating up; regional maps of the greenness of vegetation tell scientists where and when plants are photosynthesizing. There is always uncertainty associated with each value in any such map due to various factors. When uncertainty is fully modeled, instead of a single value at each map location, there is a distribution expressing a set of possible outcomes at each location. We consider such distribution data as multi-valued data since it consists of a collection of values about a single variable. Thus, a multi-valued data represents both the map and its uncertainty. We have been working on ways to visualize spatial multi-valued data sets effectively for fields with regularly spaced units or grid cells such as those in NASA's Earth science applications. A new way to display distributions at multiple grid locations is to project the distributions from an individual row, column or other user-selectable straight transect from the 2D domain. First at each grid cell in a given slice (row, column or transect), we compute a smooth density estimate from the underlying data. Such a density estimate for the probability density function (PDF) is generally more useful than a histogram, which is a classic density estimate. Then, the collection of PDFs along a given slice are presented vertically above the slice and form a wall. To minimize occlusion of intersecting slices, the corresponding walls are positioned at the far edges of the boundary. The PDF wall depicts the shapes of the distributions very dearly since peaks represent the modes (or bumps) in the PDFs. We've defined roughness as the number of peaks in the distribution. Roughness is another useful summary information for multimodal distributions. The uncertainty of the multi-valued data can also be interpreted by the number of peaks and the widths of the peaks as shown by the PDF walls.

Kao, David↗

Estimating the density of states of frustrated spin systems

Abstract Estimating the density of states (DOS) of systems with rugged free energy landscapes is a notoriously difficult task of the utmost importance in many areas of physics ranging from spin glasses to biopolymers. DOS estimation has also recently become an indispensable tool for the benchmarking of quantum annealers when these function as samplers. Some of the standard approaches suffer from a spurious convergence of the estimates to metastable minima, and these cases are particularly hard to detect. Here, we introduce a sampling technique based on population annealing enhanced with a multi-histogram analysis and report on its performance for spin glasses. We demonstrate its ability to overcome the pitfalls of other entropic samplers, resulting in some cases in large scaling advantages that can lead to the uncovering of new physics. The new technique avoids some inherent difficulties in established approaches and can be applied to a wide range of systems without relevant tailoring requirements. Benchmarking of the studied techniques is facilitated by the introduction of several schemes that allow us to achieve exact counts of the degeneracies of the tested instances.

Barash, Lev (ORCID:000000022298785X)↗

Detecting Anomalies in Time Series Using Kernel Density Approaches

This paper introduces a novel anomaly detection approach tailored for time series data with exclusive reliance on normal events during training. Our key innovation lies in the application of kernel-density estimation (KDE) to scrutinize reconstruction errors, providing an empirically derived probability distribution for normal events post-reconstruction. This non-parametric density estimation technique offers a nuanced understanding of anomaly detection, differentiating it from prevalent threshold-based mechanisms in existing methodologies. In post-training, events are encoded, decoded, and evaluated against the estimated density, providing a comprehensive notion of normality. In addition, we propose a data augmentation strategy involving variational autoencoder-generated events and a smoothing step for enhanced model robustness. The significance of our autoencoder-based approach is evident in its capacity to learn normal representation without prior anomaly knowledge. Through the KDE step on reconstruction errors, our method addresses the versatility of anomalies, departing from assumptions tied to larger reconstruction errors for anomalous events. Our proposed likelihood measure then distinguishes normal from anomalous events, providing a concise yet comprehensive anomaly detection solution. The extensive experimental results support the feasibility of our proposed method, yielding significantly improved classification performance by nearly 10% on the UCR benchmark data.

Frehner, Robin↗

Estimated probability density functions for the times between flashes in the storms of 12 September 1975, 26 August 1975, and 13 July 1976

A report is given to supplement the progress report of June 17, 1977. In that progress report gamma, lognormal, and Rayleigh probability density functions were fitted to the times between lightning flashes in the storms of 9/12/75, 8/26/75, and 7/13/76 by the maximum likelihood method. The goodness of fit is checked by the Kolmogoroff-Smirnoff test. Plots of the estimated densities along with normalized histograms are included to provide a visual check on the goodness of fit. The lognormal densities are the most peaked and have the highest tails. This results in the best fit to the normalized histogram in most cases. The Rayleigh densities have too broad and rounded peaks to give good fits. In addition, they have the lowest tails. The gamma densities fall inbetween and give the best fit in a few cases.

Tretter, S. A.↗

Directional Statistics of Preferential Orientations of Two Shapes in Their Aggregate and Its Application to Nanoparticle Aggregation

Nanoscientists have long conjectured that adjacent nanoparticles aggregate with one another in certain preferential directions during chemical syntheses of nanoparticles, which is so called the oriented attachment. For the study of the oriented attachment phenomenon, the microscopy and nanoscience community have used dynamic electron microscopy for direct observations of nanoparticle aggregation and have been so far relying on the manual and qualitative analysis of the observations. In this paper, we present a statistical approach for studying the oriented attachment phenomenon with multiple aggregation examples observed by dynamic microscope data. We abstract a nanoparticle aggregation by an event of two primary geometric objects merging into a secondary geometric object. We define a metric space for geometrical bodies, and the endowed metric in the space allows us to define the likelihood function of the orientation angles of two primary objects within the corresponding secondary object. We maximize the likelihood to estimate the orientation angles. We fit a mixture of bivariate von Mises distributions to the estimated orientation angles for the density estimation. By analyzing the density estimate, we are able to demonstrate that nanoparticles of certain geometries have indeed preferential orientations in their aggregates.

Mixture of bivariate von Mises distributions, Bi-a↗

Foundations of Feature Selection and Classification for Non-Gaussian Distributed Targets

This report describes the author’s research toward creating algorithms for selecting useful statistical features in pattern/target classification problems in which the features are non-Gaussian distributed. In engineering practice, it is common to either (1) Not perform any feature selection procedure, or (2) Use a feature selection algorithm that assumes the features are Gaussian distributed. This process can be far from optimal if the data are actually non-Gaussian, as they frequently are in practice. This research has the goal of mitigating that problem by creating a useful algorithm that can be used in practice. The approach is to focus on the performance indices used in common feature selection algorithms, such as Branch and Bound algorithms and the Sequential Forward Selection Algorithm. Ordinarily, the performance index used to measure the class separation in feature space involves assuming the data are Gaussian and deriving tractable performance indices that use closed form equations that represent the probability density functions of the class data. The advantage of this approach is that it produces feature selection algorithms that have low computational complexity and do not require knowledge of the actual data densities. The disadvantage is that these algorithms may not perform reasonably when the data are non-Gaussian, as they commonly are in practice. This research examines the use of information-theoretic class separability measures that can deal with the non-Gaussian case. In particular, this work shows that the Hellinger Distance (a divergence measure), has very desirable mathematical properties and can be useful for feature selection. The work in this paper employs a multivariate kernel density estimator along with the Hellinger Distance and an optimal subset selection algorithm to do feature selection. This report (1) Defines the classification and feature selection problems, (2) Describes the continuous variable Hellinger distance in one dimension and discusses its mathematical properties in comparison with other distance measures, (3) Develops the discrete variable Hellinger distance expression for one and two dimensions, and shows that the discrete Hellinger distance for multiple dimensions is a straightforward extension, and (4) Presents experimental density estimation results using a new MATLAB implementation of a multivariate kernel density estimator. The draft of this report served as a starting point for Ensign Matthew Wilder’s MSEE thesis work at the Naval Postgraduate School (NPS) in 2010 and 2011. Ensign Wilder completed his thesis and graduated from NPS in June 2011 (M. J. Wilder, “Automatic Target Recognition: Statistical Feature Selection of Non-Gaussian Distributed Target Classes,” MSEE Thesis, Naval Postgraduate School, Thesis Advisor: Grace A. Clark, Second Reader: Monique P. Fargues, June 2011). In the thesis, Ensign Wilder created implementations of the Branch and Bound and Sequential Feature Selection algorithms, extended the multivariate kernel density estimator code, and evaluated the performance of Hellinger Distance-based feature selection algorithms using both simulated data and data from a real-world benchmark data set. The new algorithm was shown to be very effective for both Gaussian and non-Gaussian data sets. The advantages gained for non-Gaussian data come at the cost of increased computational complexity. Using a laptop computer for practical problems, we generally use initial feature vectors containing no more than about ten or twelve features. Proposed future research lies in creating faster kernel density estimator algorithms that would allow for the use of higher dimensional data.

42 ENGINEERING↗

An analysis of lead-free (Bi{sub 0.5}Na{sub 0.5}){sub 0.915}-(Bi{sub 0.5}K{sub 0.5}){sub 0.05}Ba{sub 0.02}Sr{sub 0.015}TiO{sub 3} ceramic for efficient refrigeration and thermal energy harvesting

This article demonstrates the colossal energy harvesting capability of a lead-free (Bi{sub 0.5}Na{sub 0.5}){sub 0.915}-(Bi{sub 0.5}K{sub 0.5}){sub 0.05}Ba{sub 0.02}Sr{sub 0.015}TiO{sub 3} ceramic using the Olsen cycle. The maximum harvestable energy density estimated for this system is found to be 1523 J/L (1523 kJ/m{sup 3}) where the results are presented for extreme ambient conditions of 20–160 °C and electric fields of 0.1–4 MV/m. This estimated energy density is 1.7 times higher than the maximum reported to date for the lanthanum-doped lead zirconate titanate (thin film) system. Moreover, this study introduces a generalized and effective solid state refrigeration cycle in contrast to the ferroelectric Ericson refrigeration cycle. The cycle is based on a temperature induced polarization change on application of an unipolar electric field to ferroelectric ceramics.

36 MATERIALS SCIENCE↗