SEARCH · Search NASA

Results for “Generalized Method of Cells”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Cathodes and electrolytes for rechargeable magnesium batteries and methods of manufacture

The invention relates to Chevrel-phase materials and methods of preparing these materials utilizing a precursor approach. The Chevrel-phase materials are useful in assembling electrodes, e.g., cathodes, for use in electrochemical cells, such as rechargeable batteries. The Chevrel-phase materials have a general formula of Mo 6 Z 8 (Z=sulfur) or Mo 6 Z 1 8-yZ 2 y (Z 1 =sulfur; Z 2 =selenium), and partially cuprated Cu 1 Mo 6 S 8 as well as partially de-cuprated Cu 1-x Mg x Mo 6 S 8 and the precursors have a general formula of M x Mo 6 Z 8 or M x Mo 6 Z 1 8-yZ 2 y , M=Cu. The cathode containing the Chevrel-phase material in accordance with the invention can be combined with a magnesium-containing anode and an electrolyte.

Kumta, Prashant N.↗

DOE OSTI2020

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Investigation on capacity loss mechanisms of lithium-ion pouch cells under mechanical indentation conditions

Capacity loss was observed in Li-ion cells after mechanical deformation approaching the onset of internal short circuit (ISCr). In this paper, a series of indentation tests were carried out on commercial Li-ion cells of three capacities (500, 1500 and 2000 mAh). Both in-situ and ex-situ methods were used to investigate the mechanisms of indentation-induced capacity loss. We note after indentation test, the cell capacity reduced by 0.5%–6% of its original value. The incremental capacity (IC) analysis results showed that IC curves generally shifted to lower voltage region, indicating the increase in cell internal resistance. In addition, the fitting results of electrochemical impedance spectroscopy (EIS) indicated that mechanical indentation can result in a reduction in ohmic resistance and the increase in polarization resistance. Scanning electron microscopy (SEM) and X-ray computed tomography (XCT) results showed crushing of graphite, mud cracking of copper current collectors and enlarged pores in separator, which is proposed to be the main reasons for the increase in polarization resistance and permanent capacity loss. The rapid capacity loss due to mechanical abuse was compared with the long-term capacity fading.

25 ENERGY STORAGE↗

DOE OSTI2020

Semi-transparent p-type barium copper sulfide as a back contact interface layer for cadmium telluride solar cells

Optically transparent p-type materials play a critical role in transparent electronics including photovoltaic (PV) devices. P-type sulfide materials offer an alternative to oxides for PV application due to improved hole transport properties. Here, we report the solution-based synthesis of earth-abundant p-type transparent conducting barium copper sulfide (α-BaCu 4 S 3 , BCS) thin films. These films were characterized using scanning electron microscopy, X-ray diffraction, UV–Vis–NIR spectrophotometry, Raman spectroscopy, and spectroscopic ellipsometry. BCS films of ~100 nm thickness transmit >70% of visible light. We report on tests of the hole transport properties of these BCS films for cadmium telluride (CdTe) photovoltaics, finding that the BCS deposition process forms a beneficial tellurium (Te) rich surface on CdTe by selectively removing Cd from the surface. Based on our study, the BCS interface layer plays dual functions for CdTe PV devices as a hole transport material and as an etchant, enhancing the resulting device performance. We observed a significant increase in open-circuit voltage of CdTe solar cells with the BCS buffer layer. Furthermore, we discuss semitransparent CdTe solar cells with BCS as a hole transport layer and indium tin oxide as a finishing electrode. Semitransparent CdTe solar cells shows 13.3% conversion efficiency for the front side illumination and 1.2% efficiency for back side illumination, indicating high recombination of charge carriers generated close to the rear CdTe/BCS/ITO contact.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Similarity of capacitive radio-frequency discharges in nonlocal regimes

Similarity transformations are essential for correlating discharges at different scales, which are mostly utilized with local field or local energy approximations. In this work, we report the fully kinetic results from particle-in-cell/Monte Carlo collision simulations that unambiguously demonstrate the similarity of radio frequency (rf) discharges in nonlocal regimes where the electron energy relaxation length is much larger than the gap dimension. It is found that at a constant rf voltage amplitude, discharges will be similar if the gas pressure, inverse of gap distance, and rf driving frequency are all changed by the same scaling factor. The scaling relations of fundamental parameters are illustrated for rf discharges in the alpha-mode with secondary electron emission ignored, and the temporal electron kinetics are shown to have invariance in similar discharges. The results explicitly validate the scaling laws in nonlocal kinetic regimes, indicating promising application potentials of the similarity transformations across a wide range of kinetic regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Adaptive Angular Quadrature Scheme for a backwards-in-time Method of Characteristics Solution to the Radiative Transfer Equation [Slides]

Radiative transfer/radiation transport are important problems to solve in astrophysics and high energy density physics. Various methods exist to solve radiation transport, such as Monte Carlo (MC), Discrete Ordinates (S N ), Method of Characteristics (MOC), and the spherical harmonics (P N ) method. Method of Characteristics requires “launching” of rays in discrete directions. Unresolved details of angular mesh create ray effects and can miss sources in the domain. Ray effects can lead to unphysical “stepping” in solution and incorrect energy deposition. Adaptive quadrature schemes can be used to detect and mitigate these effects. The Method of Characteristics (MOC) is a common method for solving hyperbolic PDEs in radiation transport and supersonic flow problems. Generally in MOC for radiation transport, virtual particles are tracked from birth to the end of a timestep. This requires interpolation to go from final location to cell averaged or corner values of angular intensity. Backwards-in-Time (BIT) particle tracking avoids this by prescribing the final position of the virtual particle at the cell nodes/corners. Angular intensities are computed at time k + 1 by launching ray back to previous timestep(s), or t = 0. Scheme allows solution to be computed as the characteristic ray is traced backwards in time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Phonons and phase symmetries in bulk CrCl 3 from scattering measurements and theory

Phonon-derived behaviors are important indicators of novel phenomena in transition metal trihalides, including spin liquid behavior, two-dimensional magnetism, and spin-lattice coupling. However, phonons and their dependence on spin structure and excitations have not been adequately explored. In this work, we probe and critically examine the vibrational properties of the prototype ferromagnetic honeycomb lattice material CrCl 3 using inelastic neutron scattering and density functional theory. We demonstrate that magnetic and van der Waals interactions are essential to describing the structure and phonons in CrCl 3 ; however, the specific spin configuration is unimportant. Further, this provides context for understanding thermal transport measurements as governed by dynamical spin-lattice couplings. More importantly, we introduce an efficient dynamic method that exploits translational symmetries in large conventional unit cells that generates insights into phonon dispersions, interactions, and measured spectra in terms of quantum phase interference conditions. This work opens new avenues for understanding phonons in layered magnets and more generally in conventional cell geometries of a variety of materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Deep material network via a quilting strategy: visualization for explainability and recursive training for improved accuracy

Recent developments integrating micromechanics and neural networks offer promising paths for rapid predictions of the response of heterogeneous materials with similar accuracy as direct numerical simulations. The deep material network is one such approaches, featuring a multi-layer network and micromechanics building blocks trained on anisotropic linear elastic properties. Once trained, the network acts as a reduced-order model, which can extrapolate the material’s behavior to more general constitutive laws, including nonlinear behaviors, without the need to be retrained. However, current training methods initialize network parameters randomly, incurring inevitable training and calibration errors. Here, we introduce a way to visualize the network parameters as an analogous unit cell and use this visualization to “quilt” patches of shallower networks to initialize deeper networks for a recursive training strategy. The result is an improvement in the accuracy and calibration performance of the network and an intuitive visual representation of the network for better explainability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Evaluating the Accuracy of Various Irradiance Models in Detecting Soiling of Irradiance Sensors

We evaluate the feasibility of using various clear-sky models or purchased satellite data for estimating the soiling of a reference cell irradiance sensor. We find results to be more accurate for models that consider local meteorological conditions. We conclude that given the data sets considered, and depending on the requirements of the data analyst, choosing to use purchased satellite irradiance data from Solargis to estimate the soiling of a reference cell sensor tends to yield more accurate results, although there are instances where a clear-sky model performs better. The SOLIS clear-sky model in PVLIB with variable P wat provided useful soiling results, implying that the general method of using a clear-sky model with local meteorological data may provide a low-cost tool for detecting soiling of irradiance sensors.

14 SOLAR ENERGY↗

DOE OSTI2022

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Sucrose‐Based Dense, Pure, and Highly‐Crystalline Graphitic Materials for Lithium‐Ion Batteries

At present, most synthetic graphite materials commonly used as anode active ingredients in lithium-ion cells are produced by graphitization of petroleum cokes. The carbon footprint associated with synthetic graphite production is significant. Thus, bio-derived and cheap precursors, such as saccharides, would be an attractive alternative for the sustainable production of graphitic carbons. However, they are non-graphitizing at temperatures as high as 3000 °C, preserving the curved, fullerene-like structure of graphene layers and microporosity. Consequently, many lithium ions are consumed during the formation of solid electrolyte interphase films and passivated in the nanovoids. Here, a method for the production of pure, crystalline, graphitic materials based on sucrose disposed of microporosity is presented, which also works with a variety of saccharides and other organic precursors of hard carbons—generally considered incapable of such transformation. This process employs catalytic graphitization by Si particles at high temperatures. Finally, the electrochemical response of such derived sucrose-based graphite in Li-ion half-cells demonstrated its feasibility to serve as an anode active material for rechargeable Li-ion batteries.

Si particles↗

DOE OSTI2020

Finding wombling boundaries in LHC data with Voronoi and Delaunay tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Advanced diagnostics to evaluate heterogeneity in lithium-ion battery modules

Battery packs for electric and stationary applications experience varied operating conditions, including abuse—e.g., fast charging, overcharging, thermal, vibration, shock, etc.—throughout their lifetimes. Innovative diagnostic tools and algorithms that go beyond single cells and deal with modules and packs are essential for early detection of off-normal issues. High-resolution tools with known detection limits are key to developing appropriate mitigation strategies. With the advent of rapid impedance spectroscopy that can generate a broadband impedance spectrum in ~10 sec, the case for impedance-based diagnostics that can be readily aligned with other methods, such as incremental capacity or dQ.dV -1 , has become promising. This study used the aforementioned diagnostic methods to identify realistic in-vehicle battery abnormalities (e.g., localized self-discharge and non-uniform aging), in series (up to 10S) and parallel (4P) strings, using 16 Ah graphite/NMC cells. The impedance-based diagnostic is found to be sensitive to the string size and state. Depending on the type of abnormality, detection frequency varied. The dQ.dV -1 method showed the potential to detect long-term aging-related heterogeneity in modules. In general, both the impedance and dQ.dV -1 methods were able to detect series strings’ abnormalities, but struggled to find those issues within parallel modules.

25 ENERGY STORAGE↗

DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗