Search NASA⌕ Search

SEARCH · Search NASA

Results for “Normal 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Revisiting short-term earthquake triggered volcanism

It has been noted for centuries that earthquakes appear to trigger the eruption of volcanoes. For example, analyses of global volcanic and seismic records since 1500 AD have shown that explosive eruptions with Volcanic Explosivity Index (VEI) values ≥2 are preceded within days by nearby major earthquakes (magnitude M8 or larger) about 4 times more often than expected due to coincidence, suggesting that large earthquakes can trigger eruptions. Here we expand the definition of a triggered eruption to include the possibility of M6 or greater earthquakes within 5 days and 800 km of a VEI 2 or greater eruption. Removing pre-1964 records, to ensure complete and accurate catalogs, we find 30 volcanoes that at some point experienced a potentially triggered eruption and define these volcanoes as “sensitive” volcanoes. Within this group of sensitive volcanoes, normalized distributions of volcano-centric factors such as tectonic setting, dominant rock type, and type of volcano are practically indistinguishable from those of sensitive volcanoes in which the time of eruption is randomized. Comparisons of sensitive volcanoes and insensitive volcanoes (i.e., volcanoes that have never experienced a triggered eruption) reveal that sensitive volcanoes are simply more active than insensitive volcanoes: They erupt more frequently, are located solely in subduction zones, and erupt primarily andesites and basaltic-andesites. The potentially triggered eruptions do not show the magnitude-distance relationship expected for seismically induced responses (e.g., hydrologic responses), and eruptions that do occur within days of nearby earthquakes do so within rates expected by random chance. There is, however, a 5–12% increase in the number of explosive eruptions in the 2 months to 2 years following major earthquakes. Finally, we conclude that short-term seismically triggered explosive eruptions occur less frequently than previously inferred, an important conclusion when considering volcanic hazards and for understanding the nature of earthquake-volcano interactions.

58 GEOSCIENCES↗
DOE OSTI2018

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.

58 GEOSCIENCES↗
DOE OSTI2018

On the Five-Moment Hamburger Maximum Entropy Reconstruction

We consider the Maximum Entropy Reconstruction (MER) as a solution to the five-moment truncated Hamburger moment problem in one dimension. In the case of five monomial moment constraints, the probability density function (PDF) of the MER takes the form of the exponential of a quartic polynomial. This implies a possible bimodal structure in regions of moment space. An analytical model is developed for the MER PDF applicable near a known singular line in a centered, two-component, third- and fourth-order moment (μ{sub 3}, μ{sub 4}) space, consistent with the general problem of five moments. The model consists of the superposition of a perturbed, centered Gaussian PDF and a small-amplitude packet of PDF-density, called the outlying moment packet (OMP), sitting far from the mean. Asymptotic solutions are obtained which predict the shape of the perturbed Gaussian and both the amplitude and position on the real line of the OMP. The asymptotic solutions show that the presence of the OMP gives rise to an MER solution that is singular along a line in (μ{sub 3}, μ{sub 4}) space emanating from, but not including, the point representing a standard normal distribution, or thermodynamic equilibrium. We use this analysis of the OMP to develop a numerical regularization of the MER, creating a procedure we call the Hybrid MER (HMER). Compared with the MER, the HMER is a significant improvement in terms of robustness and efficiency while preserving accuracy in its prediction of other important distribution features, such as higher order moments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Study of the Thermomechanical Performance of FR-4 Laminates During the Reflow Process

The behavior of FR4 laminates for printed circuit boards (PCBs) during a double reflow process and the main factors associated with warpage are discussed. The microstructure, glass transition temperature (T{sub g}), decomposition temperature (T{sub d}), coefficient of thermal expansion (CTE), time to delamination and % water absorption were characterized, and the temperature profiles in the board during its pass through the reflow oven were recorded. Bow and twist measurements were taken on a set of 30 PCBs. Both T{sub d} and T{sub g} were found to be below their nominal values. The CTE is highly anisotropic and vary at different locations, and it is more than three times higher in the through thickness middle section with respect to the bottom and top locations. Temporary temperature gradients of more than 50 degrees between different locations were detected at the start of the reflow step as a result, the bow values were found to be strongly anisotropic, particularly after the first reflow, in which different sides can present bow values up to more than six times higher than others, with maxima values of 0.47 mm. The anisotropy of the CTE appears to affect only slightly the twist values. Compared with the first reflow, the bow and twist values are lower, exhibit less dispersion and are more homogeneous for the second reflow. The histograms of both bow and twist do not seem to follow a normal distribution, possibly as a result of the combined effect of spatial variations and anisotropy of properties as well as the temperature gradients during reflow.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Tannic Acid: a Novel Calibrator for Facile and Accurate Mass Measurement of Electrospray Ionization Mass Spectrometry

Accurate mass calibration is beneficial to the identification of the unknown compounds quickly and accurately. The ESI mass spectrum of tannic acid (TA) tends to a normal distribution of the cluster ion peaks in m/z range from 371.0368 to 1739.1169. Based on the interesting result, we reported the use of TA, a natural plant polyphenol, as a novel calibrator for electrospray ionization mass spectrometry (ESI MS), which has the following three advantages, including (1) easy preparation, (2) the calibration range of m/z 200~2000, and (3) the calibration error is around 3.00 ppm in positive ion mode, which is less than the use of sodium formate (SF) and Prod #88323 calibrators. This TA calibrator has great potential for the wide applications in biological, chemical, and pharmacal analysis. .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Efficient computation of basic sums for random polydispersed composites

The main goal of the paper was to develop algorithms and methods for computation of basic sums, the mathematical objects of great importance in computational materials science having applications to description of the representative volume element (RVE) and to the effective properties of 2D composites. The previously used algorithm had the exponential complexity. We propose a linearly complex algorithm. All the presented algorithms can be easily implemented in modern scientific computing packages, while maintaining both efficient calculations and a high level of abstraction. The proposed approach is applied to derivation of a polynomial approximation of the effective conductivity formula for 2D random material with non-overlapping circular inclusions with normally distributed radii. The obtained formulas are applied to the optimal packing problem of disks on the plane.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Disordered crystals from first principles I: Quantifying the configuration space

This work represents the first chapter of a project on the foundations of first-principle calculations of the electron transport in crystals at finite temperatures. We are interested in the range of temperatures, where most electronic components operate, that is, room temperature and above. The aim is a predictive first-principle formalism that combines ab-initio molecular dynamics and a finite-temperature Kubo-formula for homogeneous thermodynamic phases. The input for this formula is the ergodic dynamical system (Ω,G,dP) defining the thermodynamic crystalline phase, where Ω is the configuration space for the atomic degrees of freedom, G is the space group acting on Ω and dP is the ergodic Gibbs measure relative to the G-action. The present work develops an algorithmic method for quantifying (Ω,G,dP) from first principles. Using the silicon crystal as a working example, we find the Gibbs measure to be extremely well characterized by a multivariate normal distribution, which can be quantified using a small number of parameters. The latter are computed at various temperatures and communicated in the form of a table. Using this table, one can generate large and accurate thermally-disordered atomic configurations to serve, for example, as input for subsequent simulations of the electronic degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2019

A probabilistic perspective on thermodynamic parameter uncertainties: Understanding aqueous speciation of mercury

Speciation plays an important role in determining the fate and transport of metals in terrestrial surface and subsurface systems. Equilibrium speciation modeling in aqueous systems relies on thermodynamic constants (log K values) of complexes, which are subject to uncertainties. In this study, using Monte Carlo (MC) simulations with Latin hypercube sampling (LHS) we systematically analyze the propagation of thermodynamic constant uncertainty through speciation modeling of an inorganic mercury-sulfide-chloride-water system. We find that seemingly small variances of the input log K normal distributions can lead to output species concentrations spanning multiple orders of magnitude, with highly skewed probability distributions. When equilibrium with mineral metacinnabar (β-HgS(s)) is neglected, the relative uncertainty of each output species is strongly positively correlated with the skewness of its concentration probability distribution, i.e., the lowest uncertainty occurs when the species concentration probability distribution is the most negatively skewed as its concentration approaches the total element concentration limit. The highest uncertainty in the identity of dominant species is located around species equivalence points. For cases where the mineral equilibrium is included, we derive analytical probability density functions for the log concentrations of all major species in the system. The mineral log K uncertainty is found to be an important contributor to all output concentration uncertainties. The analysis of combined effects of both pH and total sulfide concentration on output concentration uncertainties shows that high concentration uncertainties occur under highly sulfidic alkaline conditions and low uncertainties at low pH and sulfide concentrations. Finally, an analysis as presented here can distinguish between conditions that require a full uncertainty analysis and those for which the classical deterministic speciation modeling suffices.

58 GEOSCIENCES↗
DOE OSTI2023

Determining the tensile strength of fuel surrogate TRISO-coated particle buffer, IPyC, and buffer-IPyC interlayer regions

A novel micro tensile sample fabrication technique for determining the tensile strength of the buffer, IPyC, and buffer -IPyC interlayer regions of surrogate (ZrO2) TRISO fuel particle layers was refined and implemented. Copper micro tensile samples served as baseline materials to verify the methods used. Tensile tests performed in this study, while limited in number, were analyzed using standard and Weibull statistics. As expected, the buffer layer was weakest, with an average ultimate tensile strength of 138.70 MPa, and the IPyC layer samples, were strongest, with an average ultimate tensile strength of 189.74 MPa. In the buffer -IPyC interface samples, all breaks occurred in the buffer region, though the average ultimate tensile strength of the samples, 159.80 MPa, was between the pure buffer and IPyC samples. These results suggest the interlayer region has unique properties, perhaps associated with pyrocarbon infiltration into the buffer layer during particle coating. All interlayer samples fractured within the buffer side; however, the stress strain behavior of some of these samples resembled the behavior of the IPyC layer samples. Here, the buffer and IPyC layer strengths had a normal distribution under Weibull analysis, while the interlayer region had a Rayleigh distribution. Further testing is needed to clarify both the standard and Weibull statistical results

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Improving Enzyme Optimum Temperature Prediction with Resampling Strategies and Ensemble Learning

Accurate prediction of the optimal catalytic temperature ( T opt ) of enzymes is vital in biotechnology, as enzymes with high T opt values are desired for enhanced reaction rates. Recently, a machine learning method (temperature optima for microorganisms and enzymes, TOME) for predicting T opt was developed. TOME was trained on a normally distributed data set with a median T opt of 37 °C and less than 5% of T opt values above 85 °C, limiting the method’s predictive capabilities for thermostable enzymes. Due to the distribution of the training data, the mean squared error on T opt values greater than 85 °C is nearly an order of magnitude higher than the error on values between 30 and 50 °C. Here, we apply ensemble learning and resampling strategies that tackle the data imbalance to significantly decrease the error on high T opt values (>85 °C) by 60% and increase the overall R 2 value from 0.527 to 0.632. The revised method, temperature optima for enzymes with resampling (TOMER), and the resampling strategies applied in this work are freely available to other researchers as Python packages on GitHub.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Surface Morphology and Electrical Properties of Cu 3 BTC 2 Thin Films Before and After Reaction with TCNQ

HKUST-1 or Cu 3 BTC 2 (BTC = 1,3,5-benzenetricarboxylate) is a prototypical metal–organic framework (MOF) that holds a privileged position among MOFs for device applications, as it can be deposited as thin films on various substrates and surfaces. Recently, new potential applications in electronics have emerged for this material when HKUST-1 was demonstrated to become electrically conductive upon infiltration with 7,7,8,8-tetracyanoquinodimethane (TCNQ). However, the factors that control the morphology and reactivity of the thin films are unknown. Here, we present a study of the thin-film growth process on indium tin oxide and amorphous Si prior to infiltration. From the unusual bimodal, non-log-normal distribution of crystal domain sizes, we conclude that the nucleation of new layers of Cu 3 BTC 2 is greatly enhanced by surface defects and thus difficult to control. We then show that these films can react with methanolic TCNQ solutions to form dense films of the coordination polymer Cu(TCNQ). This chemical conversion is accompanied by dramatic changes in surface morphology, from a surface dominated by truncated octahedra to randomly oriented thin platelets. The change in morphology suggests that the chemical reaction occurs in the liquid phase and is independent of the starting surface morphology. The chemical transformation is accompanied by 10 orders of magnitude change in electrical conductivity, from <10 –11 S/cm for the parent Cu 3 BTC 2 material to 10 –1 S/cm for the resulting Cu(TCNQ) film. In conclusion, the conversion of Cu 3 BTC 2 films, which can be grown and patterned on a variety of (nonplanar) substrates, to Cu(TCNQ) opens the door for the facile fabrication of more complex electronic devices.

2D nucleation↗
DOE OSTI2020

Parameterizing Subgrid Variations of Land Surface Heat Fluxes to the Atmosphere Improves Boreal Summer Land Precipitation Simulation With the NCAR CESM1.2

Subgrid horizontal variations of land surface heat fluxes to the atmosphere resulting from subgrid land cover heterogeneity are important in land-atmosphere interaction in global climate models (GCMs). To incorporate it, a parameterization using stochastic sampling based on truncated normal distributions diagnosed from the land model and internal ensemble mean of multiple calls to the planetary boundary layer and deep convection schemes is developed. After its implementation in the National Center for Atmospheric Research (NCAR) Community Earth System Model version 1.2 (CESM1.2), large changes of simulated land precipitation are found to occur in regions where large subgrid variations of surface heat fluxes exist. The simulated precipitation is improved in boreal summer, especially over eastern China and the coastal areas of the Bay of Bengal. The improved precipitation is mainly a result of improved large-scale moisture convergence and advection by altered vertical diffusion and convection.

58 GEOSCIENCES↗
DOE OSTI2018

Detection of vapor nanobubbles by small angle neutron scattering (SANS)

Experiments using boiling water on untreated (roughness 100–300 nm) metal surfaces using small-angle neutron scattering (SANS) show the appearance of structures that are 50–70 nm in size when boiling is present. The scattering signal disappears when the boiling ceases, and no change in the signal is detected at any surface temperature condition below saturation. This confirms that the signal is caused by vapor nanobubbles. Two boiling regimes are evaluated herein that differ by the degree of subcooling (3–10 °C). A polydisperse spherical model with a log-normal distribution fits the SANS data well. In conclusion, the size distribution indicates that a large number of nanobubbles exist on the surface during boiling, and some of them grow into large bubbles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗
DOE OSTI2025

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗
DOE OSTI2019

Stochastic particle production in a de Sitter background

Here, we explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass m 2 (t) of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass M of the spectator field: the conformal case with M 2 =2H 2 , and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.

79 ASTRONOMY AND ASTROPHYSICS↗