Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle heterogeneities”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental Study of Microgel Conformance-Control Treatment for a Polymer-Flooding Reservoir Containing Superpermeable Channels

Summary Polymer flooding has been widely used to improve oil recovery. However, its effectiveness would be diminished when channels (e.g., fractures, fracture-like channels, void-space conduits) are present in a reservoir. In this study, we designed a series of particular sandwich-like channel models and tested the effectiveness and applicable conditions of micrometer-sized preformed particle gels (PPGs, or microgels) in improving the polymer-flooding efficiency. We studied the selective penetration and placement of the microgel particles, and their abilities for fluid diversion and oil-recovery improvement. The results suggest that polymer flooding alone would be inefficient to achieve a satisfactory oil recovery as the heterogeneity of the reservoir becomes more serious (e.g., permeability contrast kc/km > 50). The polymer solution would vainly flow through the channels and leave the majority of oil in the matrices behind. Additional conformance-treatment efforts are required. We tried to inject microgels in an attempt to shut off the channels. After the microgel treatment, impressive improvement of the polymer-flooding performance was observed in some of our experiments. The water cut could be reduced significantly by as high as nearly 40%, and the sweep efficiency and overall oil recovery of the polymer flood were improved. The conditions under which the microgel-treatment strategy was effective were further explored. We observed that the microgels form an external impermeable cake at the very beginning of microgel injection and prevent the gel particles from entering the matrices. Instead, the microgel particles could selectively penetrate and shut off the superpermeable channels under proper conditions. Our results suggest that the 260-µm microgel particles tested in this study are effective to attack the excessive-water-production problem and improve the oil recovery when the channel has a high permeability (>50 darcies). The gels are unlikely to be effective for channels that are less than 30 darcies because of the penetration/transport difficulties. After the gels effectively penetrate and shut off the superpermeable channel, the subsequent polymer solution is diverted to the matrices (i.e., the unswept oil zones) to displace the bypassed oil. Overall, this study provides important insights to help achieve successful polymer-flooding applications in reservoirs with superpermeable channels.

Engineering↗
DOE OSTI2024

Graph neural network for neutrino physics event reconstruction

Liquid argon time projection chamber (LArTPC) detector technology offers a wealth of high-resolution information on particle interactions, and leveraging that information to its full potential requires sophisticated automated reconstruction techniques. Here, this article describes NUGRAPH 2, a graph neural network for low-level reconstruction of simulated neutrino interactions in a LArTPC detector. Simulated neutrino interactions in the MicroBooNE detector geometry are described as heterogeneous graphs, with energy depositions on each detector plane forming nodes on planar subgraphs. The network utilizes a multihead attention message-passing mechanism to perform background filtering and semantic labeling on these graph nodes, identifying those associated with the primary physics interaction with 98.0% efficiency and labeling them according to particle type with 94.9% efficiency. The network operates directly on detector observables across multiple two-dimensional representations but utilizes a three-dimensional-context-aware mechanism to encourage consistency between these representations. Model inference takes 0.12 s / event on a CPU and 0.005 s / event batched on a GPU. This architecture is designed to be a general-purpose solution for particle reconstruction in neutrino physics, with the potential for deployment across a broad range of detector technologies, and offers a core convolution engine that can be leveraged for a variety of tasks beyond the two described in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗
NASA NTRS1982

The Halley dust model

The properties of dust ejecta from Comet Halley are studied on the basis of (a) evidence from the comet's past apparitions and (b) analogy with recent, physically similar comets. Specifically discussed are the light curve and spectrum, discrete phenomena in the head, the physical properties of the nucleus (size, albedo, rotation, surface temperature, and morphology), and an interaction between the nucleus and dust atmosphere. Also reviewed are constraints on the size and mass distributions of dust particles, information on submicron-size and submillimeter-size grains from the comet's dust tail and antitail, and the apparent existence of more than one particle type. Similarities between the jet patterns of Halley and the parent comet of the Perseid meteor stream are depicted, and effects of the surface heterogeneity (discrete active regions) on the dust flow are assessed. Current dust models for Halley are summarized and the existence of short-term variations in the dust content in the comet's atmosphere is suggested.

Sekanina, Z.↗
NASA NTRS2016

HST-STIS Spectra of Saturn's Rings and Implications for Their Reddening Agent

We obtained HST-STIS spectra of Saturn's main rings in May 2011, using the G230L (and G430L) gratings, with final averaged radial resolution of 160 (and 330) km/pixel. The dataset filled a previous 200-330nm "spectral gap" between Cassini and ground-based spectra. The data provide radial profiles as a function of wavelength, but our most basic product at this point is a set of very low-noise spectra, radially averaged over broad regions of the rings (A, B, C, and Cassini Division). The raw spectra required special processing to remove artifacts due to extended-source grating scatter. We have modeled the spectra using a new particle surface model, which corrects for on-surface shadowing due to the likely very rough ring particle surfaces, and avoids overestimation of intra-mixed "neutral absorber". We correct for non-classical layer effects and finite ring optical depth, and relate our observed reflectivities to the spherical albedos of individual smooth particles. We model these smooth particle albedos using standard Hapke theory for regolith grain mixtures that are either homogeneous and "intramixed" (nonicy absorbers dispersed in water ice regolith grains) or heterogeneous "intimate" mixtures. As candidates for the nonicy contaminants we have considered amorphous carbon, aromatic-rich and aliphatic-rich organic tholins, silicates, hematite and iron metal. For the A and B rings, we find that iron metal (including a new theoretical estimate of the refractive indices of nanometer-sized grains of iron) is not spectrally steep enough in the 200-300nm range, and that aliphatic-rich tholins are either too steep at short wavelengths or too flat at long wavelengths. However, less than 1% by mass of aromatic-rich tholins provides a very good fit across the entire spectral range with no gratuitous "neutral absorber" needed, and a minimum of additional free parameters. The best fits require forward-scattering regolith grains. For the C Ring and Cassini Division, additional absorbers are needed (updated results will be given).

Saturn's Rings↗
NASA NTRSDate not supplied

Airborne High Spectral Resolution Lidar Measurements of Aerosol Distributions and Properties During the NASA CAMP2Ex Mission

The NASA Clouds, Aerosols and Monsoon Processes-Philippines Experiment (CAMP2Ex) occurred during August-October, 2019,primarily over the ocean near the Philippines, including the ocean west of Luzon. CAMP2Ex attempts to characterize the role of anthropogenic and natural aerosol particles in modulating warm and mixed-phase precipitation in the maritime continent during the Southwest Monsoon. This mission is also examining relationships between spatial heterogeneity in the aerosol and cloud fields and radiative heating rates and studying the effects of aerosol concentration and composition on the microphysical properties of shallow cumulus and cumulus congestus clouds. Critical to the successful investigation of these objectives is profiling of aerosol, cloud and meteorological properties throughout the atmospheric column. For the CAMP2Ex experiment, we deployed the second generation NASA Langley airborne High Spectral Resolution Lidar (HSRL-2) on the NASA P-3 aircraft to obtain profile measurements of aerosol and thin cloud extinction and optical thickness via the HSRL technique at 355 and 532 nm, and of backscatter and depolarization at 355, 532, and 1064 nm. The HSRL technique allows for the independent measurement of extinction and backscatter without the need to assume modeled lidar ratios. Additional HSRL-2 data products include aerosol type, mixed layer depth, and range-resolved aerosol microphysical parameters. The HSRL measurements reveal the temporal, spatial, and vertical variability of aerosol optical properties over these locations. The HSRL-2 measured these backscatter, extinction, and depolarization profiles during all nineteen science flights and the transit flights to and from the aircraft operations base at Clark International Airport.In this presentation, we use these HSRL-2 measurements to show the horizontal and vertical distributions of aerosols in the CAMP2Ex operational area as well as provide insight into the optical and physical characteristics of these aerosols. We also discuss how these aerosols are distributed within and above the mixed layer.

aerosol↗
DOE OSTI2022

Current‐Dependent Lithium Metal Growth Modes in “Anode‐Free” Solid‐State Batteries at the Cu|LLZO Interface

Abstract Controlling the lithium growth morphology in lithium reservoir‐free cells (RFCs), so‐called “anode‐free” solid‐state batteries, is of key interest to ensure stable battery operation. Despite several benefits of RFCs like improved energy density and easier fabrication, issues during the charging of the cell hinder the transition from lithium metal batteries with a lithium reservoir layer to RFCs. In RFCs, the lithium metal anode is plated during the first charging step at the interface between a metal current collector and the solid electrolyte, which is prone to highly heterogeneous growth instead of the desired homogeneous film‐like growth. Herein, the lithium morphology during the first charging step in RFCs is explored as a function of current density and current collector thickness. Using operando scanning electron microscopy, an increase in the lithium particle density is observed with increasing current density at the Cu|Li 6.25 Al 0.25 La 3 Zr 2 O 12 interface. This observation is then applied to improve the area coverage of lithium by pulsed plating. It is also shown that thin current collectors ( d = 100 nm) are unsuited for RFCs, as lithium whiskers penetrate them, resulting in highly heterogeneous interfaces. This suggests the use of thicker metal layers (several µm) to mitigate whisker penetration and facilitate homogeneous lithium plating.

25 ENERGY STORAGE↗
DOE OSTI2021

Impacts of long-range-transported mineral dust on summertime convective cloud and precipitation: a case study over the Taiwan region

As one of the most abundant atmospheric aerosols and effective ice nuclei, mineral dust affects clouds and precipitation in the Earth system. Here numerical experiments are carried out to investigate the impacts of dust aerosols on summertime convective clouds and precipitation over the mountainous region of Taiwan by acting as ice-nucleating particles. We run the Weather Research and Forecasting model (WRF) with the Morrison two-moment and spectral-bin microphysics (SBM) schemes at 3 km resolution, using dust number concentrations from a global chemical transport model (GEOS-Chem-APM). The case study indicates that the long-range-transported mineral dust, with relatively low number concentrations, can notably affect the properties of convective clouds (ice and liquid water contents, cloud top height, and cloud coverage) and precipitation (spatial pattern and intensity). The effects of dust are evident during strong convective periods, with significantly increased ice water contents in the mixed-phase regime via the enhanced heterogeneous freezing. With both the Morrison and SBM schemes, we see the invigoration effects of dust aerosols on the convective intensity through enhanced condensation and deposition latent heating. The low-altitude dust particles are uplifted to the freezing level by updrafts, which, in turn, enhance the convective cloud development through immersion freezing and convective invigoration. Compared to the Morrison scheme, the SBM scheme predicts more realistic precipitation and different invigoration effects of dust. The differences are partially attributed to the saturation adjustment approach utilized in the bulk scheme, which leads to a stronger enhancement of condensation at midlatitudes to low altitudes and a weaker deposition increase at the upper level.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Experimental and theoretical analysis of carbon driven detonation waves in a heterogeneously premixed Rotating Detonation Engine

Coal dust explosions can be hazardous; however, they can also generate a significant rise in stagnation pressure if adequately harnessed. Rotating detonation combustors seek to take advantage of the stagnation pressure rise phenomenon in a more sustained and controlled manner via confinement to a physical annulus, leading to increased overall thermodynamic efficiency. Here this investigation presents an analysis of detonations fueled by Carbon Black, a solid particulate consisting of virtually pure carbon molecules and lean Hydrogen-Air mixtures. It is realized that with the addition of Carbon Black, an increase of lean mixture detonability and detonation velocities extending the operating limit over that of a pure hydrogen-air mixture is experienced. For all testing conditions, the total equivalence ratio is held at φ = 1, while the fuel mixture's carbon mass fraction is increased from 0 to 0.7 while the hydrogen is decreased. Detonation wave velocities are extracted from high-speed imaging through applying a Discrete Fourier Transform algorithm to determine changes to the wave speed as Carbon Black particles are introduced. As a result, due to the addition of Carbon Black as an auxiliary fuel source, detonations were formed instead of deflagrations in operating conditions where one would expect deflagrations at the same hydrogen-air equivalence ratios without Carbon Black addition. The detonation formation provides evidence that the coal particles are reacting within the detonation wave in a large enough capacity to support a detonation wave within the annulus. Furthermore, the wave speed is shown to increase with the additional of carbon particles. At a constant global equivalence ratio, the detonation wave velocities were found to decrease with hydrogen's incremental replacement with coal particles. Whereby, through a theoretical comparison of the heat of combustion as computed from the experimentally derived detonation wave velocities, a linear relationship of the two was shown to exist. Therefore, the heat of combustion has the potential to describe an operational limit to sustaining a detonation wave.

42 ENGINEERING↗
DOE OSTI2022

Particulate Fuel Modeling of MC 2 -3 using Iterative Local Spatial Self-shielding Method

We report a new spatial self-shielding method for particulate fuels has been developed based on disadvantage factors and implemented in the MC 2 -3 code. This method named the iterative local spatial self-shielding (ILSS) method considers the shadowing effect of randomly distributed particles on spatial self-shielding in particles through a homogenized composition region added outside the particle of interest at the center. The self-shielded cross sections of the central particle are determined iteratively since they are used in determining the cross sections of the homogenized composition region. The ILSS method was verified for infinite stochastic medium problems of single and multiple types of particles, VHTR unit cell problems, and HTTR assembly problems. The verification test results show that the ILSS method accurately predicts the stochastic particle shadowing effect and reaction rates in particles, whereas the regular array model and the stochastic collision probability method underpredict the particle shadowing effect and overestimate reaction rates in particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Dependent scattering and fractal microstructure determine the transparency of aerogel monoliths

This study reveals how dependent scattering and microstructure significantly affect electromagnetic wave propagation through aerogel monoliths, contributing to their transparency. Light scattering by particle ensembles is considered “dependent” when the scattering properties rely not only on particle size and optical constants but also on their spatial distribution, typically occurring when the average interparticle distance is small in comparison with the wavelength of incident radiation. Addressing dependent scattering requires solving Maxwell’s equations for complex heterogeneous structures, which is computationally demanding and usually limited to sample thicknesses on the same scale as the wavelength. This study combines computer-generated ambigel microstructures of fractal aggregates of polydisperse nanoparticles and the radiative transfer with reciprocal transaction method to predict the transmittance of thick ambigel slabs. Transmittance measurements of ambiently dried aerogel monoliths (ambigels) with porosities from about 50% to 90% closely matched the predicted values for their digital twins. However, ignoring dependent scattering or particle aggregation led to inaccurate predictions. This study validated the computational framework, and its findings offer insights for designing photonic metamaterials and analyzing their interactions with electromagnetic waves.

Yalcin, Refet A. (ORCID:0000000339973494)↗
DOE OSTI2023

Reaction Current Heterogeneity at the Interface between a Lithium Electrode and Polymer/Ceramic Composite Electrolytes

Although lithium metal anodes are expected to increase the energy density of next-generation batteries, dendrite growth during charge remains a major bottleneck preventing widespread implementation. Composite solid electrolytes with ceramic particles embedded in a polymer matrix have the potential to prevent dendrites owing to the higher mechanical stiffness while also possessing the flexibility to maintain contact with the electrode. However, microscopically, the different mechanical and electrochemical properties of the polymer and ceramic domains can cause inhomogeneous charge transfer at the Li/electrolyte interface, which can lead to nonuniform Li deposition and propagation of dendrites. Here, we computationally examine the coupled electrochemical, transport, and mechanical processes at the interface to determine the propensity for dendrite formation and possible approaches to mitigate this issue. Predictions of two possible microstructures at the interface, namely, (i) where both the polymer and ceramic come in contact with Li metal and (ii) when only the polymer comes in contact with Li metal, suggest that the former has a greater tendency for nonuniform plating. In addition, predictions suggest that minimizing the interfacial resistance between polymers and ceramics and incorporating interlayers between the electrode and electrolyte help mitigate current heterogeneity. In conclusion, these predictions provide guidance for experimental approaches to prevent dendrites in composite electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2021

Supermolecular structures of recrystallized starches with amylopectin side chains modified by amylosucrase to different chain lengths

Understanding the supermolecular structures of recrystallized starches is imperative to manufacture novel starches with controllable physicochemical properties enabling novel applications. In this study, native amylopectin (AP) from maize was modified to different chain lengths using amylosucrase from Neisseria polysaccharea. Native AP granules showed well-organized lamellar structures with a periodicity of ca. 9.0 nm whereas modified starches formed heterogeneous structures with the average crystallite size ranging from 19.77 to 29.88 nm. Acidic treatments eroded amorphous domains in both native AP granules and modified starch particles, producing nanocrystals composed of double helices. A-type nanocrystals (~40 nm) from native AP granules showed quadrangular shape with an acute angle of ca. 60°, while irregular B-type nanocrystals were observed for the modified starches. Molecular characterizations suggested that the nanocrystals had relatively uniform chain lengths (DP 16.5–25.2), and the helix length (the thickness of nanocrystals) was positively correlated to the thermal stability of starches.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2004

Determination of the Three-Dimensional Morphology of ALH84001 and Biogenic MV-1 Magnetite: Comparison of Results from Electron Tomography and Classical Transmission Electron Microscopy

Dated at approximately 3.9 billion years of age, carbonate disks, found within fractures of the host rock of Martian meteorite ALH84001, have been interpreted as secondary minerals that formed at low temperature in an aqueous medium. Heterogeneously distributed within these disks are magnetite nanocrystals that are of Martian origin. Approximately one quarter of these magnetites have morphological and chemical similarities to magnetite particles produced by magnetotactic bacteria strain MV-1, which are ubiquitous in aquatic habitats on Earth. Moreover, these types of magnetite particles are not known or expected to be produced by abiotic means either through geological processes or synthetically in the laboratory. The remaining three quarters of the ALH84001 magnetites are likely products of multiple processes including, but not limited to, precipitation from a hydrothermal fluid, thermal decomposition of the carbonate matrix in which they are embedded, and extracellular formation by dissimilatory Fe-reducing bacteria. We have proposed that the origins of magnetites in ALH84001 can be best explained as the products of multiple processes, one of which is biological. Recently the three-dimensional (3-D) external morphology of the purported biogenic fraction of the ALH84001 magnetites has been the subject of considerable debate. We report here the 3-D geometry of biogenic magnetite crystals extracted from MV-1 and of those extracted from ALH84001 carbonate disks using a combination of high resolution classical and tomographic transmission electron microscopy (TEM). We focus on answering the following questions: (1) which technique provides adequate information to deduce the 3-D external crystal morphology?; and, (2) what is the precise 3-D geometry of the ALH84001 and MV-1 magnetites?

Thomas-Keprta, Kathie L.↗
DOE OSTI2023

A Lattice Boltzmann Method for Predicting Porous Transport Layer Performance During Electrolysis

Electrolysis, the splitting of water into oxygen and hydrogen using electricity, is a sustainable way to produce green hydrogen for energy storage. In polymer electrolyte membrane (PEM) water electrolysis, water is brought into contact with charged catalyst layers and electrochemically separated into oxygen and hydrogen. The hydrogen product formed at the cathode is carried through the catalyst layer for eventual collection, while the oxygen by-product formed at the anode is removed from the surface via a multiphase interaction with circulating water and a solid porous transport layer (PTL). The design of this PTL aids in the detachment and advection of the oxygen by-product and thereby plays a role in the overall efficacy of the catalyst. In this presentation, we present our initial results modeling this multiphase system using a single-component, multiphase lattice Boltzmann method. We use the Shan-Chen model describing inter-particle forces to capture both the cohesion of the water (liquid) and oxygen (gas) phases and their interaction with the PTL (solid) (Shan and Chen, 1993). We use a Carnahan-Starling equation of state to model the effective density governing these inter-particle interactions which allows us to model this relatively high density ratio system (Carnahan and Starling, 1969). With these simulations, we show that the geometry and heterogeneity of the PTL geometry plays a large role in its ability to move oxygen away from the catalyst layer and the resulting bubble structures that are formed within the PTL. The current work demonstrates these effects using synthesized PTL geometries and 2D physics, which will be extended to experimentally-imaged PTL sections and 3D algorithms in the near future.

Boltzmann↗
DOE OSTI2024

Image-Based Failure Assessment of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a heterogeneous Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. IM-RK is additionally used to inform how crack evolution in turn affects the coupled electro-chemo-mechanical behavior of the Li-ion battery cathode.

image-based modeling↗
NASA NTRS2001

Mechanistic Studies of Combustion and Structure Formation During Synthesis of Advanced Materials

Combustion in a variety of heterogeneous systems, leading to the synthesis of advanced materials, is characterized by high temperatures (2000-3500 K) and heating rates (up to 10(exp 6) K/s) at and ahead of the reaction front. These high temperatures generate liquids and gases which are subject to gravity-driven flow. The removal of such gravitational effects is likely to provide increased control of the reaction front, with a consequent improvement in control of the microstructure of the synthesized products. Thus, microgravity (mu-g) experiments lead to major advances in the understanding of fundamental aspects of combustion and structure formation under the extreme conditions of the combustion synthesis (CS) wave. In addition, the specific features of microgravity environment allow one to produce unique materials, which cannot be obtained under terrestrial conditions. The current research is a logic continuation of our previous work on investigations of the fundamental phenomena of combustion and structure formation that occur at the high temperatures achieved in a CS wave. Our research is being conducted in three main directions: 1) Microstructural Transformations during Combustion Synthesis of Metal-Ceramic Composites. The studies are devoted to the investigation of particle growth during CS of intermetallic-ceramic composites, synthesized from nickel, aluminum, titanium, and boron metal reactants. To determine the mechanisms of particle growth, the investigation varies the relative amount of components in the initial mixture to yield combustion wave products with different ratios of solid and liquid phases, under 1g and mu-g conditions; 2) Mechanisms of Heat Transfer during Reactions in Heterogeneous Media. Specifically, new phenomena of gasless combustion wave propagation in heterogeneous media with porosity higher than that achievable in normal gravity conditions, are being studied. Two types of mixtures are investigated: clad powders, where contact between reactants occurs within each particle, and mixtures of elemental powders, where interparticle contacts are important for the reaction; and 3) Mechanistic Studies of Phase Separation in Combustion of Thermite Systems. Studies are devoted to experiments on thermite systems (metal oxide-reducing metal) where phase separation processes occur to produce alloys with tailored compositions and properties. The separation may be either gravity-driven or due to surface forces, and systematic studies to elucidate the true mechanism are being conducted. The knowledge obtained will be used to find the most promising ways of controlling the microstructure and properties of combustion-synthesized materials. Low-gravity experiments are essential to create idealized an environment for insights into the physics and chemistry of advanced material synthesis processes.

Varma, A.↗