Search NASA⌕ Search

SEARCH · Search NASA

Results for “methods - statistical”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

New method for analysis of X-ray computed tomography scans of TRISO fuel forms

X-ray computed tomography (XCT) is a powerful characterization tool that has been used for examination of both tristructural-isotropic (TRISO) coated particles and fuel forms. Rapid, nondestructive imaging of complete specimens by XCT provides significant characterization advantages in both speed and comprehensiveness over traditional materialographic methods. This paper discusses a method that has been developed to analyze XCT images of pebble fuel forms containing ~18,000 TRISO particles each. The XCT images are three-dimensional with contrast based on relative x-ray transparency. In the case of TRISO fuel forms, XCT image intensity distinguishes between the high-density kernel, the medium-density SiC, and the low-density matrix and carbon layers. Image analysis software has been developed to process XCT data sets to extract spatial distribution information, including particle nearest-neighbor distances, local packing fraction, and fuel-free zone thickness in TRISO fuel forms. Finally, although this type of high-resolution XCT analysis is likely not feasible as a quality control measurement for commercial TRISO fuel, it is highly useful in generating feedback for the development of fuel form fabrication processes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High-dimensional methods for quantum homodyne tomography

We provide optimized recursion relations for homodyne tomography. We improve previous methods by mitigating the divergences intrinsic in the calculation of the pattern functions used previously, and detail how to implement the data analysis through Monte Carlo simulations. Our refinements are necessary for the reconstruction of excited quantum states which populate a high-dimensional subspace of the electromagnetic field Hilbert space. •Stabilization of numerics for quantum homodyne tomography reconstructions.•High dimensional quantum homodyne tomography reconstruction.•Monte Carlo simulations of quantum homodyne tomography.•Quantum Homodyne tomography simulations.

Mosco, Nicola↗
DOE OSTI2015

Method for measurement of diffusivity: Calorimetric studies of Fe/Ni multilayer thin films

A calorimetric method for the measurement of diffusivity in thin film multilayers is introduced and applied to the Fe Ni system. Using this method, the diffusivity in [Fe (25 nm)/Ni (25 nm)](20) multilayer thin films is measured as 4 x 10(-3)exp(-1.6 +/- 0.1 eV/ k(B)T) cm(2)/s, respectively. The diffusion mechanism in the multilayers and its relevance to laboratory synthesis of L1(0) ordered FeNi are discussed. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Liu, JX↗
DOE OSTI2019

Spin-State Ordering in Metal-Based Compounds Using the Localized Active Space Self-Consistent Field Method

Quantitatively accurate calculations for spin-state ordering in transition-metal complexes typically demand a robust multiconfigurational treatment. The poor scaling of such methods with increasing size makes them impractical for large, strongly correlated systems. Density matrix embedding theory (DMET) is a fragmentation approach that can be used to specifically address this challenge. The single-determinantal bath framework of DMET is applicable in many situations, but it has been shown to perform poorly for molecules characterized by strong correlation when a multiconfigurational self-consistent field solver is used. To ameliorate this problem, the localized active space self-consistent field (LASSCF) method was recently described. Here in this work, LASSCF is applied to predict spin-state energetics in mono- and di-iron systems, and we show that the model offers an accuracy equivalent to that of CASSCF but at a substantially lower computational cost. Performance as a function of basis set and active space is also examined.

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

Electrochemical Methods for Assessing Kinetic Factors in the Reduction of CO 2 to Formate: Implications for Improving Electrocatalyst Design

hermochemical insights are often employed in the rationalization of reactivity and in the design of electrocatalysts for CO 2 reduction reactions targeting C–H bond-containing products. This work identifies experimental methods for assessing kinetic aspects of reactivity. These methods are illustrated using [Fe 4 N(CO) 12 ] - , which produces formate from CO 2 at -1.2 V versus SCE in either a MeCN/H 2 O solvent (95:5) or pH 6.5 buffered water. Elementary rates for each reaction step are identified along with the rate-determining step (RDS) as C–H bond formation. Transition state kinetics were determined from an Eyring analysis for the rate-determining C–H bond formation step using temperature-dependent electrochemical measurements. A lower measured Δ G ‡ (298 K, 12.3 ± 0.1 kcal mol –1 ) in a pH 6.5 aqueous solution, compared with a Δ G ‡ (298 K) of 15.0 ± 0.1 kcal mol –1 in a MeCN/H 2 O solvent (95:5), correlates with faster observed reaction rates and provides a kinetic rationalization for the solvent-dependent chemistry. Taken together, the experimentally determined kinetic insights highlight that enhancement of the local concentration of CO 2 at catalyst–hydride sites should be a primary focus of ongoing catalyst design. This will both enhance reaction rates and increase selectivity for C–H bond formation over competing H–H bond formation, because that step is fast in H 2 evolution reactions.

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

Exciton Absorption Spectra by Linear Response Methods: Application to Conjugated Polymers

The theoretical description of the time-evolution of excitons requires, as an initial step, the calculation of their spectra, which has been inaccessible to most users due to the high computational scaling of conventional algorithms and accuracy issues caused by common density functionals. Previously (J. Chem. Phys.2016, 144, 204105), we developed a simple method that resolves these issues. Our scheme is based on a two-step calculation in which a linear-response TDDFT calculation is used to generate orbitals perturbed by the excitonic state, and then a second linear-response TDDFT calculation is used to determine the spectrum of excitations relative to the excitonic state. Herein, we apply this theory to study near-infrared absorption spectra of excitons in oligomers of the ubiquitous conjugated polymers poly(3-hexylthiophene) (P3HT), poly(2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene) (MEH-PPV), and poly(benzodithiophene-thieno[3,4-b]thiophene) (PTB7). For P3HT and MEH-PPV oligomers, the calculated intense absorption bands converge at the longest wavelengths for 10 monomer units, and show strong consistency with experimental measurements. The calculations confirm that the exciton spectral features in MEH-PPV overlap with those of the bipolaron formation. In addition, our calculations identify the exciton absorption bands in transient absorption spectra measured by our group for oligomers (1, 2, and 3 units) of PTB7. For all of the cases studied, we report the dominant orbital excitations contributing to the optically active excited state–excited state transitions, and suggest a simple rule to identify absorption peaks at the longest wavelengths. We suggest our methodology could be considered for further developments in theoretical transient spectroscopy to include nonadiabatic effects, coherences, and to describe the formation of species such as charge-transfer states and polaron pairs.

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

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗
DOE OSTI2023

An energy-optimization method to study gel-swelling in confinement

Using our open-source programming environment Morpho, an energy-optimization method is developed to compute equilibrium shapes of an arbitrarily confined swollen hydrogel. The internal strains and contact forces are measured, probing its mechanics.

Joshi, Chaitanya↗
DOE OSTI2013

On a Riesz basis of exponentials related to the eigenvalues of an analytic operator and application to a non-selfadjoint problem deduced from a perturbation method for sound radiation

In the present paper, we prove that the family of exponentials associated to the eigenvalues of the perturbed operator T(ε) ≔ T{sub 0} + εT{sub 1} + ε{sup 2}T{sub 2} + … + ε{sup k}T{sub k} + … forms a Riesz basis in L{sup 2}(0, T), T > 0, where ε∈C, T{sub 0} is a closed densely defined linear operator on a separable Hilbert space H with domain D(T{sub 0}) having isolated eigenvalues with multiplicity one, while T{sub 1}, T{sub 2}, … are linear operators on H having the same domain D⊃D(T{sub 0}) and satisfying a specific growing inequality. After that, we generalize this result using a H-Lipschitz function. As application, we consider a non-selfadjoint problem deduced from a perturbation method for sound radiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Full dimensional Franck-Condon factors for the acetylene A{sup ~} {sup 1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. I. Method for calculating polyatomic linear—bent vibrational intensity factors and evaluation of calculated intensities for the gerade vibrational modes in acetylene

Franck-Condon vibrational overlap integrals for the A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition in acetylene have been calculated in full dimension in the harmonic normal mode basis. The calculation uses the method of generating functions first developed for polyatomic Franck-Condon factors by Sharp and Rosenstock [J. Chem. Phys. 41(11), 3453–3463 (1964)], and previously applied to acetylene by Watson [J. Mol. Spectrosc. 207(2), 276–284 (2001)] in a reduced-dimension calculation. Because the transition involves a large change in the equilibrium geometry of the electronic states, two different types of corrections to the coordinate transformation are considered to first order: corrections for axis-switching between the Cartesian molecular frames and corrections for the curvilinear nature of the normal modes at large amplitude. The angular factor in the wave function for the out-of-plane component of the trans bending mode, ν{sub 4}{sup ″}, is treated as a rotation, which results in an Eckart constraint on the polar coordinates of the bending modes. To simplify the calculation, the other degenerate bending mode, ν{sub 5}{sup ″}, is integrated in the Cartesian basis and later transformed to the constrained polar coordinate basis, restoring the conventional v and l quantum numbers. An updated A{sup ~}-state harmonic force field obtained recently in the R. W. Field research group is evaluated. The results for transitions involving the gerade vibrational modes are in qualitative agreement with experiment. Calculated results for transitions involving ungerade modes are presented in Paper II of this series [G. B. Park, J. H. Baraban, and R. W. Field, “Full dimensional Franck–Condon factors for the acetylene A{sup ~1}A{sub u}—X{sup ~1}Σ{sup +}{sub g} transition. II. Vibrational overlap factors for levels involving excitation in ungerade modes,” J. Chem. Phys. 141, 134305 (2014)].

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

A complex guided spectral transform Lanczos method for studying quantum resonance states

A complex guided spectral transform Lanczos (cGSTL) algorithm is proposed to compute both bound and resonance states including energies, widths, and wavefunctions. The algorithm comprises of two layers of complex-symmetric Lanczos iterations. A short inner layer iteration produces a set of complex formally orthogonal Lanczos polynomials. They are used to span the guided spectral transform function determined by a retarded Green operator. An outer layer iteration is then carried out with the transform function to compute the eigen-pairs of the system. The guided spectral transform function is designed to have the same wavefunctions as the eigenstates of the original Hamiltonian in the spectral range of interest. Therefore, the energies and/or widths of bound or resonance states can be easily computed with their wavefunctions or by using a root-searching method from the guided spectral transform surface. The new cGSTL algorithm is applied to bound and resonance states of HO{sub 2}, and compared to previous calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Enhanced power factor of higher manganese silicide via melt spin synthesis method

We report on the thermoelectric properties of the Higher Manganese Silicide MnSi₁.₇₅ (HMS) synthesized by means of a one-step non-equilibrium method. The ultrahigh cooling rate generated from the melt-spin technique is found to be effective in reducing second phases, which are inevitable during the traditional solid state diffusion processes. Aside from being detrimental to thermoelectric properties, second phases skew the revealing of the intrinsic properties of this class of materials, for example the optimal level of carrier concentration. With this melt-spin sample, we are able to formulate a simple model based on a single parabolic band that can well describe the carrier concentration dependence of the Seebeck coefficient and power factor of the data reported in the literature. An optimal carrier concentration around 5x10²⁰ cm⁻³ at 300 K is predicted according to this model. The phase-pure melt-spin sample shows the largest power factor at high temperature, resulting in the highest zT value among the three samples in this paper; the maximum value is superior to those reported in the literatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Analytical solutions of the planar cyclic voltammetry process for two soluble species with equal diffusivities and fast electron transfer using the method of eigenfunction expansions

Cyclic voltammetry is a powerful tool that is used for characterizing electrochemical processes. Models of cyclic voltammetry take into account the mass transport of species and the kinetics at the electrode surface. Analytical solutions of these models are not well-known due to the complexity of the boundary conditions. In this study we present closed form analytical solutions of the planar voltammetry model for two soluble species with fast electron transfer and equal diffusivities using the eigenfunction expansion method. Our solution methodology does not incorporate Laplace transforms and yields good agreement with the numerical solution. This solution method can be extended to cases that are more general and may be useful for benchmarking purposes.

36 MATERIALS SCIENCE↗