Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lognormal distribution”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A statistical representation of bond coating oxidation under environmental barrier coatings

Environmental barrier coatings (EBCs) protect SiC-based ceramic matrix composites (CMCs) in turbine hot sections from high-temperature volatilization in combustion gases. The formation of a SiO 2 thermally grown oxide (TGO) is expected under the EBC after long-term operation. The oxidation resistance of the EBC is understood as a life-limiting factor for the CMC, and this work predicts long-term oxidation behavior under EBCs through a simple statistical approach. Specimens were exposed to 1350°C isothermal conditions for 100-h thermal cycles in flowing steam for up to 1000 h. The EBC morphology, SiO 2 thickness, and SiO 2 cracking behavior were assessed. Using thousands of SiO 2 thickness measurements across many millimeters of the interface, a realistic representation of the entire TGO was captured via a lognormal distribution. The lognormal fit parameters were extrapolated out to 25 000 h to assess the degree of SiO 2 growth, the spread of SiO 2 thicknesses related to the rough oxidizing interface, and percentages of the intermediate bond coating consumed. In conclusion, local interfacial defects from the coating deposition process are identified as local failure points for EBC—CMC systems.

SiO 2↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗
DOE OSTI2020

PV O&M Cost Model [SWR-21-14]

This tool presents a method for calculating costs associated with the operation and maintenance (O&M) of photovoltaic (PV) systems. The tool compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($/year) for each year of an analysis period, the net present value ($) of life cycle costs accumulated over the analysis period, and the reserve account amount ($) that might be required to fund unexpected repairs. The reserve account includes parts inventory and is important for providing a source of funds to make repairs quickly and avoid lost production. This method allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string, or central inverter), and site and environmental conditions (e.g., pollen, bird populations) which is an improvement over simple per unit valuations of O&M costs ($/kW/year). This model also distinguishes costs that vary from year to year and increase at different rates over time as modeled by heuristic failure distributions (e.g., Weibull or Lognormal distribution) based on actuarial data for many of the services.

Walker, Andy↗
DOE OSTI2020

Model of Operation-and-Maintenance Costs for Photovoltaic Systems

This article presents a method for calculating costs associated with operation and maintenance (O&M) of photovoltaic (PV) systems. It compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($\$$/year) for each year of an analysis period, the net present value ($\$$) of life cycle costs accumulated over the analysis period, and the reserve account amount ($\$$). Here we show that this method is an improvement over the previous averaged or levelized per-unit ($\$$/kW/year) valuations for estimating PV O&M costs, because it allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string or central inverter) and site and environmental conditions (e.g., snow, pollen, bird populations). This model also distinguishes costs that vary from year to year and increase at different rates over time because of heuristic failure distributions (e.g. Weibull or Lognormal distribution) based on actuarial data for many of the services. This cost model was created by the PV O&M Working Group of researchers and industry, sponsored by DOE Solar Energy Technologies Office, and has been published in an on-line version hosted by SunSpec Alliance at apsuite.sunspec.org. A spreadsheet version is included with this paper as supplementary material.

14 SOLAR ENERGY↗
DOE OSTI2023

Non-parametric projections of national income distribution consistent with the Shared Socioeconomic Pathways

Understanding and projecting income distributions within countries and regions is important to understanding consumption trends and the distributional consequences of climate impacts and responses. Several global, country-level projections of income distribution are available but most project only the Gini coefficient (a summary statistic of the distribution) or utilize the Gini along with the assumption of a lognormal distribution. We test the lognormal assumption and find that it typically underestimates income in the highest deciles and over-estimates it in others. We find that a new model based on two principal components of national time series data for income distribution provides a better fit to the data for all deciles, especially for the highest and lowest. We also construct a projection model in which the first principal component is driven by the Gini coefficient and the second captures deviations from this relationship. We use the model to project income distribution by decile for all countries for the five shared socioeconomic pathways. We find that inequality is consistently higher than projections based on the Gini and the lognormal functional form, with some countries reaching ratios of the highest to lowest income deciles that are almost three times their value using the lognormal assumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Thermophysical modeling of raw glaze liquidus temperature and viscosity

Raw glazes, which do not contain pre-melted glass frits, are widely used due to their low cost. The natural variability of the mineral compositions used in raw glazes affects their viscosity during firing. Molten glaze viscosity plays an important role in determining the final surface quality; therefore, it is important to determine how raw glaze composition and firing conditions affect viscosity. Thermophysical modeling provides a way to evaluate the liquidus temperature and viscosity of molten glazes as a function of composition and oxygen activity. Equilibrium phase volume fraction and composition of 242 raw glaze recipes were evaluated under oxidizing and reducing conditions. Undissolved or precipitation of solid phases explain the secondary flux transition from anti-fluxing to fluxing near the liquidus temperature. The liquidus temperature and viscosity of iron glazes decrease as a function of increasing iron content under a reducing atmosphere. The empirical probability distribution of molten glaze viscosity follows a lognormal distribution, with lead glazes having significantly lower viscosity compared with glazes without lead. Here, the peak firing temperature viscosity of lead-free glazes is near the working point of glass. Multiple linear regression analysis shows that the peak firing temperature viscosity and liquidus temperature significantly predict the viscosity factor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Climate Signatures on Lake And Wetland Size Distributions in Arctic Deltas

Abstract Understanding how thermokarst lakes on arctic river deltas will respond to rapid warming is critical for projecting how carbon storage and fluxes will change in those vulnerable environments. Yet, this understanding is currently limited partly due to the complexity of disentangling significant interannual variability from the longer‐term surface water signatures on the landscape, using the short summertime window of optical spaceborne observations. Here, we rigorously separate perennial lakes from ephemeral wetlands on 12 arctic deltas and report distinct size distributions and climate trends for the two waterbodies. Namely, we find a lognormal distribution for lakes and a power‐law distribution for wetlands, consistent with a simple proportionate growth model and inundated topography, respectively. Furthermore, while no trend with temperature is found for wetlands, a statistically significant decreasing trend of mean lake size with warmer temperatures is found, attributed to colder deltas having deeper and thicker permafrost preserving larger lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets

Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.

02 PETROLEUM↗
DOE OSTI2020

Dust Growth by Accretion of Molecules in Supersonic Interstellar Turbulence

We show that the growth rate of dust grains in cold molecular clouds is enhanced by the high degree of compressibility of a turbulent, dilute gas. By means of high-resolution (1024 3 ) numerical simulations, we confirm the theory that the spatial mean growth rate is proportional to the gas-density variance. This also results in broadening of the grain-size distribution (GSD) due to turbulence-induced variation of the grain-growth rate. We show, for the first time in a detailed numerical simulation of hydrodynamic turbulence, that the GSD evolves toward a shape that is a reflection of the gas-density distribution, regardless of the initial distribution. That is, in case of isothermal, rotationally forced turbulence, the GSD tends to be a lognormal distribution. We also show that in hypersonic turbulence, decoupling of gas and dust becomes important and that this leads to an even further accelerated grain growth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients: Impact of Particle Deposition in ICRP 66 Human Respiratory Tract Model

Inhaled radioactive materials can pose a long-term health concern, as the material can be incorporated into the body’s metabolic pathways and remain in organs and tissues for extended durations. During the retention period, the radioactive material may localize in a source organ and irradiate adjacent target organs and tissues. Distribution of these materials changes over time, requiring biokinetic modeling to evaluate their movement through various tissues and organs. The evolving distribution depends on multiple inputs characterizing the inhaled material, such as particle size and size distribution, particle density, aspect ratio, specific radionuclide, the chemical form, and solubility. In addition, biological parameters such as breathing rate, breathing type (nasal or nasal/oral), respiratory system morphometry, tidal volume, functional residual capacity, and anatomical dead space all influence material transport. These aerosol properties and physiological characteristics of the respiratory tract jointly define a range of initial conditions that influence the time-dependent distribution of radioactive material. To evaluate both uncertainty in the initial conditions of inhalation exposure and the final output (committed effective dose) from biokinetic models, a Python-based software tool, Radiological Exposure Dose Calculator (REDCAL), was developed to propagate uncertainty within the human respiratory tract model. Focusing on deposition fraction uncertainty, the primary objective was to characterize the initial activity distribution across respiratory regions as a function of anticipated particle sizes and distributions. The impact of the deposition fraction uncertainty was propagated to committed effective dose coefficients for selected radionuclides in a companion publication. For each particle size, a lognormal distribution, characterized by its geometric mean as defined within ICRP Publication 66, serves as the basis for introducing uncertainty into the physical processes governing deposition in various lung regions. Finally, this study addresses the deposition process and examines how uncertainty in deposition mechanisms affects activity distribution in the airways, ultimately presenting the expected range and standard deviation of deposited activity as a function of particle size.

International Commission on Radiological Protectio↗
DOE OSTI2020

Characterization of Radiation Induced Defects Across Scales

The Whole Powder Pattern Modelling (WPPM) approach of diffraction line profile analysis (LPA) was applied for the first time on PuO 2 in order to understand the evolution of the microstructural properties as a function of process conditions. Crystallite size distributions and defect concentrations were determined by WPPM on three different PuO 2 samples - including one with a bimodal distribution and compared to the result of more common LPA analyses used in PuO 2 literature studies (e.g. Scherrer equation, double Voigt method). The conclusion was that a lognormal distribution of spherical crystallites with dislocations most closely represent the microstructure of PuO 2 produced by heating Pu 2 (C 2 O 4 ) 3 (H 2 O) 6 •3H 2 O between 650°C and 850°C. The volume-weighted mean column length (Lvol) parameter obtained from the double Voigt method for the monomodal samples were the same as the Lvol obtained from the WPPM within uncertainty. The two LPA methods did not agree on the Lvol parameters of the bimodal distribution sample with only the WPPM getting values that were consistent with the individual components. The WPPM can potentially provide more insight into future studies of the influences on the microstructure properties of PuO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Constraining the scatter in the galaxy–halo connection at Milky Way masses

In this work, we develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter haloes. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive haloes. In both of our methods, we use a galaxy group finder to locate central galaxies in the Sloan Digital Sky Survey main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between haloes and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighbouring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, Vpeak, as our property to connect galaxies to haloes, and all comparisons between theory and observation are made after first passing the model through the group-finding algorithm. We parametrize scatter as a lognormal distribution in M* at fixed Vpeak, σ[M*|Vpeak]. The cross-correlation technique yields a constraint of σ[M*|Vpeak] = 0.27 ± 0.05 dex at a mean Vpeak of 168 km s-1, corresponding to a scatter in log M* at fixed Mh of σ[M*|Mh] = 0.38 ± 0.06 dex at $M_\text{h}=10^{11.8}\, \text{M}_\odot$. The kurtosis technique yields σ[M*|Vpeak] = 0.30 ± 0.03 at Vpeak = 209 km s-1, corresponding to σ[M*|Mh] = 0.34 ± 0.04 at $M_\text{h}=10^{12.2}\, \text{M}_\odot$. The values of σ[M*|Mh] are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between Vpeak and Mh, and it represents an increase of nearly a factor of 2 relative to the values inferred from clustering and group studies at high masses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Estimation of Aerosol Columnar Size Distribution from Spectral Extinction Data in Coastal and Maritime Environment

Aerosol columnar size distributions (SDs) are commonly provided by aerosol inversions based on measurements of both spectral extinction and sky radiance. These inversions developed for a fully clear sky offer few SDs for areas with abundant clouds. Here, we estimate SDs from spectral extinction data alone for cloudy coastal and maritime regions using aerosol refractive index (RI) obtained from chemical composition data. Our estimation involves finding volume and mean radius of lognormally distributed modes of an assumed bimodal size distribution through fitting of the spectral extinction data. We demonstrate that vertically integrated SDs obtained from aircraft measurements over a coastal site have distinct seasonal changes, and these changes are captured reasonably well by the estimated columnar SDs. We also demonstrate that similar seasonal changes occur at a maritime site, and columnar SDs retrieved from the combined extinction and sky radiance measurements are approximated quite well by their extinction only counterparts (correlation exceeds 0.9) during a 7-year period (2013–2019). The level of agreement between the estimated and retrieved SDs depends weakly on wavelength selection within a given spectral interval (roughly 0.4–1 µm). Since the extinction-based estimations can be performed frequently for partly cloudy skies, the number of periods where SDs can be found is greatly increased.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

SootLib: A soot model library for combustion simulation

Soot formation in combustion is an important process that affects radiative heat transfer, flame temperatures, and emissions with health and environmental impacts. Soot formation involves complex chemistry for nucleation, growth, oxidation, and coagulation processes. The soot particles vary widely in size and accurate modeling requires representation of the particle size distribution (PSD). Modeling soot is not trivial, and is only one of several physical processes active in combustion systems. This paper presents a software package called SootLib, which is an open-source library for modeling soot formation and other aerosol systems. SootLib is written in C++, is documented with Doxygen, and is available on GitHub. The library includes several models for soot chemistry and coagulation, and it represents the PSD using either a sectional model or the method of moments (MOM). Four closure approaches for the MOM are implemented allowing up to eight moments: monodispersed, an assumed-shape lognormal distribution, the quadrature method of moments, and the method of moments with interpolative closure. SootLib provides an interface for inclusion in other combustion packages including CFD or reacting flow solvers. The range of models allows comparisons and sensitivity studies, and the modularity facilitates extension to other soot models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Corrosion of rebar in concrete. Part II: Literature survey and statistical analysis of existing data on chloride threshold

The literature is reviewed with respect to defining, measuring, and interpreting the chloride threshold (CT). CT was defined with a strong fundamental basis in the Point Defect Model. which also dictates the nature of the statistical analysis of surveyed CT values. A database of CT for various primary and secondary independent variables has been established. Statistical analyses reveal that CT is lognormally distributed, in accordance with the PDM, with most probable values being 0.85, 0.69, and 1.19 for %total Cl/cem, %free Cl/cem, and [Cl - ]/[OH - ], respectively. The breakdown and corrosion potentials follow normal distributions with means of 0.056 and -0.322 V SCE , respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Parasitic Modulation of Microwave Signals by a Hypersonic Plasma Layer

During hypersonic flight, compressional and viscous heating of the air can form a plasma layer which encases the aircraft. If the boundary layer becomes turbulent, then the electron density fluctuations can effect a parasitic modulation in microwave signals transmitted through the plasma. We developed an approach for studying the interaction of microwave signals with a turbulent, hypersonic plasma layer. The approach affords a great deal of flexibility in both the plasma layer model and the antenna configuration. We then analyzed a situation in which microwaves, transmitted from a rectangular aperture antenna, propagate through a turbulent plasma layer to a distant receiver. We characterized the first-and second-order statistics of the computed parasitic modulation and quantified the depolarization of the signal. The amplitude fluctuations are lognormally distributed at low frequencies and Rice-distributed at high frequencies. Fluctuations in the copolarized phase and amplitude of the far-field signal are strongly anticorrelated. Here, we used a multioutput Gaussian process (MOGP) to model these quantities. The efficacy of the MOGP model is demonstrated by recovering the time evolution of the copolarized phase given the copolarized amplitude and occasional measurements of the phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probabilistic Assessment of Damage from Same Shock Response Spectra due to Variations in Damping

Interpretation of field data from shock tests and subsequent assessment of product safety margins via laboratory testing are based on the shock response spectra (SRS). The SRS capture how a single degree of freedom (SDOF) structure responds to the shock at differing frequencies and, therefore, no longer contain the duration or other temporal parameters pertaining to the shock. A single duration can often be included in the technical specification or in the recreation of acceleration vs. time history from the specified SRS; however, there is little basis for that beyond technical judgment. The loss of such temporal information can result in the recreated SRS being the same while its effect on a system or component can be different. This paper attempts to quantify this deficiency as well as propose a simple method of capturing damping from shock waves that can allow the original waveform to be more accurately reconstructed from the SRS. In this study the decay rate associated with various frequencies that comprise the overall shock was varied. This variation in the decay rate leads to a variation in the acceleration vs. time history, which can be correlated to a “Damage Index” that captures the fatigue damage imparted to the object under shock. Several waveforms that have the same SRS but varying rates of decay for either high- or low-frequency components of the shock were investigated. The resulting variation in stress cycles and Damage Index is discussed in the context of the lognormal distribution of fatigue failure data. It is proposed that, along with the SRS, the decay rate is also captured to minimize the discrepancy between field data and representative laboratory tests.

36 MATERIALS SCIENCE↗