Search NASA⌕ Search

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.

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

At least 307 records · Page 17

Machine learning to accelerate screening for Marcus reorganization energies

Understanding and predicting the charge transport properties of π-conjugated materials is an important challenge for designing new organic electronic devices, such as solar cells, plastic transistors, light-emitting devices, and chemical sensors. A key component of the hopping mechanism of charge transfer in these materials is the Marcus reorganization energy which serves as an activation barrier to hole or electron transfer. While modern density functional methods have proven to accurately predict trends in intramolecular reorganization energy, such calculations are computationally expensive. In this work, we outline active machine learning methods to predict computed intramolecular reorganization energies of a wide range of polythiophenes and their use toward screening new compounds with low internal reorganization energies. Our models have an overall root mean square error (RMSE) of ±0.113 eV, but a much smaller RMSE of only ±0.036 eV on the new screening set. Since the larger error derives from high-reorganization energy compounds, the new method is highly effective to screen for compounds with potentially efficient charge transport parameters.

14 SOLAR ENERGY↗

Mechanisms of Water-Stimulated Mg 2+ Intercalation in Vanadium Oxide: Toward the Development of Hydrated Vanadium Oxide Cathodes for Mg Batteries

As lithium-ion batteries approach their theoretical limits for energy density, magnesium-ion batteries are emerging as a promising next-generation energy storage technology. However, progress in magnesium-ion battery research has been stymied by a lack of available high capacity cathode materials that can reversibly insert magnesium ions. Vanadium Oxide (V 2 O 5 ) has emerged as one of the more promising candidate cathode materials, owing to its high theoretical capacity, facile synthesis methods, and relatively high operating voltage. This review focuses on the outlook of hydrated V 2 O 5 structures as a high capacity cathode material for magnesium-ion batteries. In general, V 2 O 5 structures exhibit poor experimental capacity for magnesium-ion insertion due to sluggish magnesium-ion insertion kinetics and poor electronic conductivity. However, several decades ago, it was discovered that the addition of water to organic electrolytes significantly improves magnesium-ion insertion into V 2 O 5 . This review clarifies the various mechanisms that have been used to explain this observation, from charge shielding to proton insertion, and offers an alternative explanation that examines the possible role of structural hydroxyl groups on the V 2 O 5 surface. While the mechanism still needs to be further studied, this discovery fueled new research into V 2 O 5 electrodes that incorporate water directly as a structural element. The most promising of these hydrated V 2 O 5 materials, many of which incorporate conductive additives, nanostructured architectures, and thin film morphologies, are discussed. Ultimately, however, these hydrated V 2 O 5 structures still face a significant barrier to potential applications in magnesium-ion batteries. During full cell electrochemical cycling, these hydrated structures tend to leach water into the electrolyte and passivate the surface of the magnesium anode, leading to poor cycle life and low capacity retention. Recently, some promising strides have been made to remedy this problem, including the use of artificial solid electrolyte interphase layers as an anode protection scheme, but a call to action for more anode protection strategies that are compatible with trace water and magnesium metal is required.

25 ENERGY STORAGE↗

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of the ring tension test (RTT) for robust determination of material strengths

The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted to fill this gap. The two most used test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Finally, based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Porting WarpX to GPU-accelerated platforms

WarpX is a general purpose electromagnetic particle-in-cell code that was originally designed to run on many-core CPU architectures. We describe the strategy, based on the AMReX library, followed to allow WarpX to use the GPU-accelerated nodes on OLCF’s Summit supercomputer, a strategy we believe will extend to the upcoming machines Frontier and Aurora. We summarize the challenges encountered, lessons learned, and give current performance results on a series of relevant benchmark problems.

97 MATHEMATICS AND COMPUTING↗

Light-trapping structures fabricated in situ for ultrathin III-V solar cells

Here, we describe a fully in situ method of fabricating light-scattering structures on III-V materials that generates a rough morphology via vapor phase etching and redeposition. Fully in situ methods support higher industrial throughput by utilizing the growth reactor to generate the light-trapping structures after device growth without removal from the reactor. We use HCl and PH 3 to etch and redeposit scattering morphologies on Ga 0.5 In 0.5 P in a dynamic hydride vapor phase epitaxy (D-HVPE) reactor. We show that the addition of PH 3 leads to redeposition during the vapor phase HCl etching of Ga 0.5 In 0.5 P and that HCl flow rate and time exposed to HCl-PH 3 each independently cause a linear increase in the redeposited feature size, indicating that redeposition proceeds by island growth in a III-Cl-limited, hydride-enhanced HVPE regime. Auger electron spectroscopy and scanning transmission electron microscopy with energy dispersive spectroscopy (STEM-EDS) reveal redeposition to be highly Ga-rich GaInP, i.e., Ga(In)P. The Ga-rich nature of the redeposition results from the higher thermodynamic driving force for Ga incorporation than for In during HVPE growth and the difference in the volatility of the III-Cl etch products. The resulting morphologies have high broadband scattering, as determined by normal specular reflectance and integrating sphere measurements, indicating effectiveness as light-scattering structures. In a 270-nm-thick GaAs photovoltaic device with a textured back surface, we achieve a 4.9% increase in short circuit current density (J SC ) without any loss in open-circuit voltage (V OC ) relative to a planar control using only a 60 s in situ texturing treatment.

14 SOLAR ENERGY↗

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗

Hybrid gyrokinetic ion/fluid electron simulation of toroidal tearing modes

The effects of toroidicity and kinetic ions on the resistive tearing mode are systematically studied with the gyrokinetic particle-in-cell simulation code GEM and compared with analytic theory. A new field solver in toroidal geometry has been developed for the simulation of low-n (n = 1, 2) modes in tokamaks. It is found that the toroidal effect significantly reduces the growth rate of the tearing mode. The toroidal effect can also increase the radial width of the tearing mode and change the scaling between the radial mode width and resistivity due to the toroidal pressure term in the electron continuity equation. The kinetic effects of ions can decrease the growth rate of the tearing mode. The plasma flux-surface shaping is found to have significant effect on the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗

Standards development for modules in high temperature micro-environments

Photovoltaic (PV) module qualification standards, IEC 61215 and IEC 61730, were designed to apply to “general open-air climates” and IEC 61730 specifically indicated applicability of ambient air temperature of 40°C. Additionally, these standards provided allowances for so-called “open rack mounted PV modules” without a clear definition of “open-rack.” These implied restrictions and allowances meant that hotter climates or thermally restrictive installation methods may not be covered by these often customer-mandated certification standards. This is particularly salient for the significant growth regions of the Middle East and India that would be expected to operate at significantly higher temperatures. The applicability of these documents raised issues over the definition of “open rack” and the fact that the geographic location is just as important as the mounting configuration in assessing the impact of the micro-environment of a PV module. This work summarizes the scientific background for IEC Technical Specification 63126:2020 ED1, titled “Guidelines for qualifying PV modules components and materials for operation at high temperatures.” This standard was recently published by the IEC and is the first step in a systematic effort to rework these standards to address the question of temperature more directly. Instead of specifying a mounting condition, we specify different suites of tests suitable for a system (PV module, mounting style, and location) defined by the 98th percentile cell temperature. With a defined temperature regime to work from, this allowed us to use existing literature research combined with additional modeling work to determine, which tests would need to be modified. This resulted in suggested changes to material thermal indices, thermal cycling temperatures, hot spot testing, ultraviolet testing, and bypass diode testing among other tests and characteristics described in this study.

14 SOLAR ENERGY↗

Specifications of EBR-II Neutron Radiography Method Description and Digitization Approach

All neutron radiography (NRAD) images of fuel pins in Argonne’s collection were originally generated using the NRAD imaging facility established in the Hot Fuel Examination Facility (HFEF) at Idaho National Laboratory (INL). The NRAD reactor facility was built in 1977 and has been operating since. The reactor is a TRIGA-type reactor operating at a power level of 250 kWth to provide a neutron source for radiography imaging. The reactor is equipped with two beam tubes (i.e., east beam tube and north beam tube) to guide the neutron beams to two radiography stations. The east radiography station is directly under the HFEF main cell and is dedicated for specimens already in the HFEF hot cell. The north radiography station is outside of the main HFEF hot cell and allows NRAD imaging of non-irradiated items. The NRAD images of EBR-II irradiated metallic fuel pins were taken in the east radiography station. Thermal neutrons have the capability to transmit through most materials and are ideal for NRAD imaging. However, because of their high thermal neutron absorption cross-section, fissile materials (e.g., highly-enriched nuclear fuels) may not be as transmissible to thermal neutrons. This is also the case for oversize specimens with extraneous thickness. Epithermal neutron imaging is therefore used as a complement to thermal neutron NRAD imaging. At HFEF’s NRAD facility, both thermal and epithermal neutrons can be used for NRAD imaging. Irradiated nuclear fuels emit high levels of γ radiation that can easily darken X-ray films, so direct exposure NRAD cannot be used to image them. Instead, an indirect NRAD imaging method was developed at HFEF’s NRAD facility. In this method, foils made of materials that can be activated by neutrons (i.e., with large neutron absorption cross section) are used to collect transmitted neutron signals. Then the activated foils are then placed against X-ray films and enclosed in a vacuum cassette so that the γ decay from the activated foils can produce images on the X-ray films. Then, general X-ray film processing procedures are used to digitize and store the images. By using different foil materials, different energy neutrons can be used for NRAD imaging. At the HFEF NRAD station, two types of films are commonly used: dysprosium (Dy) foils with thickness of 130 microns are used to capture thermal neutron signal, while indium (In) foils with thickness of 130 microns are used to capture epithermal neutron signal. A cadmium or gadolinium foil is put before the indium foil to work as a thermal neutron filter. The thermal and epithermal NRAD images can be taken simultaneously by using a Dy/Cd/In sandwiched foil combination. The typical NRAD exposure time is approximately 20 minutes. Then the exposed foils are transferred to film vacuum cassettes. The vacuum ensures that there is no gap between the foil and the film. The foil-to-film exposure time is at least three half-lives of the corresponding radioisotopes, which are 3 hours for In and 7.5 hours for Dy, respectively. Exposed films are processed using an automatic film processor to produce completed NRAD images.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗

Extracting Optical Bandgaps from Luminescence Images of Perovskite Solar Cells

We report a camera-based method to capture optical bandgap images of perovskite solar cells (PSCs). It is fast and non-destructive. It can also provide micron-scale spatial resolution. This imaging technique utilizes well-defined and relatively symmetrical band-to-band luminescence spectra emitted from perovskite materials. Their spectra can be fitted with Gaussian functions whose peak locations are similar to absorption thresholds and thus represent the material optical bandgaps. We apply the technique to resolve relative variations in optical bandgaps across various PSCs, and to show optical bandgap inhomogeneity within the same device due to material degradation and impurities. Degradation and impurities were found to both cause optical bandgap shifts inside the materials. Our results are independently confirmed with photoluminescence and absorption spectroscopy.

14 SOLAR ENERGY↗

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interpreting accelerated tests on perovskite modules using photooxidation of MAPbI 3 as an example

Solar panels (modules) based on metal halide perovskites are following a fast track to commercialization. Unlike more established solar cell materials, there are not yet decades-long field observations to increase consumer confidence. The "physics and chemistry of failure" approach is used in other industries and estimates product degradation based on laboratory accelerated tests integrated with an understanding of degradation mechanisms. This work uses that approach to quantify the relationship between accelerated tests and projected product behavior for metal halide perovskite modules. Degradation involving photooxidation of methylammonium lead iodide is used to illustrate the method. Acceleration factors in common accelerated tests are found to be low. Conclusions emphasize that the accelerated tests on photovoltaics should not be interpreted as equivalent across module types or as a green light for commercialization unless supported by the appropriate field data or physics and chemistry of failure analysis.

14 SOLAR ENERGY↗