Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel cell modeling”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A mathematical design framework for membrane pre-concentration in energy-efficient recovery of fermentation products

Due to the dilute nature of products manufactured via fermentation and cell-free bioprocessing, dewatering is a common unit operation in downstream processing (DSP) for bioproduct recovery, but it is typically energy intensive. To improve DSP energy efficiency for bio-based small molecules, integrating high-pressure membrane pre-concentration is a promising process option. However, this approach is typically constrained by a tradeoff between concentration factor (CF) and product recovery (PR), namely increasing the CF typically results in greater product loss, and vice versa. Here we developed a model that enables process design guidelines to: (i) identify scenarios in which the additional energy consumption and product loss from membrane pre-concentration are justified for use in DSP, and (ii) determine the optimal CF that minimizes process specific energy consumption. We compared the energy consumption of high-pressure membrane-integrated processes to evaporation-only processes and applied the model to an experimental case study for the separation and purification of butyric acid from Clostridium tyrobutyricum fermentation using an in situ product recovery (ISPR) process. The model estimated that integrating a tangential-flow reverse osmosis (RO) pre-concentration unit could reduce process energy consumption up to 45%. The use of advanced membrane pre-concentration technologies, such as negative rejection membranes and organic solvent reverse osmosis (OSRO), have the potential to further reduce the overall process specific energy consumption up to 96%, projected based on modeling. Overall, membrane pre-concentration, especially when strategically integrated prior to an evaporation step with optimized process conditions, holds significant potential for improving DSP energy efficiency, particularly in applications requiring substantial solvent removal for product recovery from dilute mixtures.

09 BIOMASS FUELS↗
DOE OSTI2020

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

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

Upgrade of Gamma Spectrometry Systems for ORNL TRISO Fuel PIE

Gamma spectrometry is a key element in much of the post-irradiation examination (PIE) work performed under the Advanced Gas Reactor Fuel Development and Qualification (AGR) Program (Demkowicz et al. 2015; Stempien et al. 2021). Gamma spectrometers are integrated into three major capabilities used at the Oak Ridge National Laboratory (ORNL) Irradiated Fuels Examination Laboratory (IFEL) for PIE of tristructural-isotropic (TRISO) coated particles and fuel compacts: the Core Conduction Cooldown Test Facility (CCCTF), the Vertical Counting System (VCS), and the Irradiated Microsphere Gamma Analyzer (IMGA). The CCCTF includes liquid-nitrogen-cooled traps to extract 85 Kr out of the He sweep gas that passes through the furnace in which the fuel compacts are heated during safety testing. Analysis of the 85 Kr activity in the traps is the primary indicator for TRISO failure during safety testing. The VCS is a system used to accurately measure gamma emission from components placed in a lead-shielded chamber. It is used to count the CCCTF deposition cups after removal from furnace. Each cup resides in the CCCTF furnace for typically 12–24 h and is periodically replaced with a fresh cup throughout the safety test. Metallic fission products collect on the water-cooled cups and several gamma-emitting isotopes ( 110 mAg, 134 Cs, 137 Cs, 154 Eu, and 155 Eu) are often measured and provide indication of the retention performance of the TRISO coatings. The VCS is also used to measure the presence of these isotopes on the CCCTF tantalum liner and sweep gas inlet tube for the determination of cup collection efficiency, as well as support other gamma spectrometry needs related to calibration of the 85 Kr fission gas traps and various other special PIE tasks. The IMGA uses gamma spectrometry to measure the inventory of gamma-emitting isotopes in individual TRISO particles. An automated particle handling system within the IMGA hot cell removes each particle from a source vial and positions it in front of a gamma detector, and output from the gamma spectrometer is used by the IMGA software to determine a destination vial such that particles are sorted according to their inventory and retention characteristics. At the conclusion of the AGR-1 and AGR-2 PIE campaigns, the gamma spectrometer systems used at ORNL to support that PIE had reached the end of its life cycle due to gradual obsolescence of the hardware and software. Upgrade of the Canberra Genie 2000 software used by these systems to a Windows 10 version was not a viable option, because the newest Windows 10 version offered by Mirion (the new owner of the Canberra technology) did not include the dynamic-link libraries (DLLs) needed for integration with the custom PIE software used with the CCCTF and IMGA, and Mirion had no current plans for development and release of Windows 10 versions of these DLLs with the Model S560 Genie 2000 Programming Library. Ultimately a switch was made to ORTEC gamma spectrometry systems, which appeared to be a more sustainable solution due to more proactive vendor support. The ORTEC conversion involved replacing the aging detector preamplifier and multichannel analyzer (MCA) hardware, upgrading the obsolete Windows 7 computers to Windows 10 compatible models, adopting ORTEC GammaVision software, and extensive modification of the ORNL-developed Visual Basic .NET (VB.NET) programs that provide the CCCTF and IMGA user interfaces.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

Xylan Is Critical for Proper Bundling and Alignment of Cellulose Microfibrils in Plant Secondary Cell Walls

Plant biomass represents an abundant and increasingly important natural resource and it mainly consists of a number of cell types that have undergone extensive secondary cell wall (SCW) formation. These cell types are abundant in the stems of Arabidopsis, a well-studied model system for hardwood, the wood of eudicot plants. The main constituents of hardwood include cellulose, lignin, and xylan, the latter in the form of glucuronoxylan (GX). The binding of GX to cellulose in the eudicot SCW represents one of the best-understood molecular interactions within plant cell walls. The evenly spaced acetylation and 4-O-methyl glucuronic acid (MeGlcA) substitutions of the xylan polymer backbone facilitates binding in a linear two-fold screw conformation to the hydrophilic side of cellulose and signifies a high level of molecular specificity. However, the wider implications of GX–cellulose interactions for cellulose network formation and SCW architecture have remained less explored. In this study, we seek to expand our knowledge on this by characterizing the cellulose microfibril organization in three well-characterized GX mutants. The selected mutants display a range of GX deficiency from mild to severe, with findings indicating even the weakest mutant having significant perturbations of the cellulose network, as visualized by both scanning electron microscopy (SEM) and atomic force microscopy (AFM). We show by image analysis that microfibril width is increased by as much as three times in the severe mutants compared to the wild type and that the degree of directional dispersion of the fibrils is approximately doubled in all the three mutants. Further, we find that these changes correlate with both altered nanomechanical properties of the SCW, as observed by AFM, and with increases in enzymatic hydrolysis. Results from this study indicate the critical role that normal GX composition has on cellulose bundle formation and cellulose organization as a whole within the SCWs.

09 BIOMASS FUELS↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2020

Open Source Visualization and Analysis Platform for 3D Reconstructions of Materials by Transmission Electron Microscopy

Three-dimensional characterization of materials at the nano- and meso-scale has become possible with transmission and scanning transmission electron microscopes (S/TEM). Its importance has extended to a wide class of nanomaterials such as hydrogen fuel cells, solar cells, industrial catalysts, new battery materials and semiconductor devices, as well as spanning high-tech industry, universities, and national labs. While capable instrumentation is abundant, this rapidly expanding demand for high-resolution tomography is bottlenecked by software that is instead tailored for lower-dose, biological applications and not optimized for higher-resolution materials applications. Existing tomography fails to utilize the chemical information provided by S/TEM spectrometers. To address this problem, this project delivered a scalable, fully functional, freely-distributable, open Source materials tomography package with a modern user interface that enables automated acquisition, alignment, and real-time reconstruction of raw tomography data, and provides advanced segmentation, three-dimensional chemical visualization and analysis optimized for materials applications. It has established an extendable framework capable of automation for high-throughput for the tomography of materials from data acquisition to visualization. Phase I and II delivered a full-featured, cross-platform, clean and integrated application for S/TEM materials tomography. It can read projection data from the microscope, with graphical tools for alignment of data, tomographic reconstruction, segmentation, and visualization of the reconstructed 3D volume. The entire pipeline can be saved to an XML state file, enabling fully reproducible data processing and analysis with a full Python environment for the development of custom algorithms, processing, and analysis requirements. Phase IIB extended the application to provide real-time tomographic reconstruction, "live updates" as data is processed, analysis of big data, and multi-channel tomography. Here, quantitative assessment of nanomaterials can occur as data is being recorded on a 3D visualization platform that accommodates additional information from multiple channels. We provide unmatched high-throughput tomography, where the complete tomographic pipeline from measurement to 3D visualization can occur rapidly. The project was extended to offer capabilities for micro-CT, X-ray, neutron, focused ion beam, atomic electron tomography and atom probe tomography. With over 600 transmission electron microscopes worldwide and approximately 50 coming online each year, the demand and impact of an open-source tomography tool is large. Significant opportunities exist in high-tech industry, universities, and national labs to enable or enhance three-dimensional imaging at the nanoscale and bring automated high-throughput approaches that will accelerate progress in materials characterization and metrology. The project supports a service based business model by enabling lab-specific acquisition and processing customization and integration-support and development that will be provided into Phase III and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microbial community diversity changes during voltage reversal repair in a 12-unit microbial fuel cell

Microbial fuel cell stacks (MFC-Stack) are often confronted with voltage reversals, likely due to an interplay between microbial community dynamics and insufficient electric circuit balancing. Herein, we provide new insight into voltage reversals by examining the microbiomes of twelve MFC units of a 12-liter Pilot-MFC-Stack during repair. Different biofilm repair methods (self-healing, electrostimulation, and re-acclimatization upon cross-inoculation) were used to evaluate the microbial community response. In addition, MFC-Stack simulation was performed based on Kirchhoff’s Second Law to predict values for source potentials and post-evaluate internal resistances. Analysis of the 16S rRNA amplicon sequencing data suggests that the biofilm repair methods could slowly heal damaged biofilms. Notably, severely voltage reversed MFC units had low electrogen relative abundances (18%) and positive anode potentials, while strong bioanodes and contained more than 50% electrogens and had negative anode potentials. Between-community analyses (beta diversity ordination and multinomial regression) of the voltage reversed MFC units revealed differences among biofilms in contrast to healthy/strong MFC units. Permutational multivariate analysis of variance (PERMANOVA) confirmed that reversed biofilms were, indeed, significantly (p < 0.05) different from stronger ones. Overall, these analyses demonstrated the utility of combining electrotechnical and microbial community analyses, especially beta diversity ordination and multinomial regression, to understand problematic MFC units and the potential success of a biofilm repair method. Finally, thicker biofilms were usually healthier and stronger, although thickness was no guarantee for proper structure and power function as all factors were interdependent. There was an evolutionary trend that strong anodes became stronger/healthier and others weaker. This spontaneous trend has to be considered to avoid irreversible voltage reversals and to repair electrogenic biofilms in an MFC-Stack.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Turboexpander for Direct Cooling in Hydrogen Vehicle Fueling Infrastructure

Hydrogen fuel cell electric vehicles (FCEVs) have been identified as one of a few options for zero carbon emissions transportation. A major advantage of FCEVs is that they can fuel quickly and follow a familiar fueling behavior to hydrocarbon-fueled vehicles. Whether light duty or heavy duty, the goal for a hydrogen dispenser is to fuel a vehicle in the same amount of time as the fossil fuel equivalent. When hydrogen is dispensed into the vehicle storage system, however, the temperature rises due to the Joule-Thomson effect and the heat of compression. Typically, vehicles store the compressed hydrogen in composite overwrapped pressure vessels that have a polymer liner with an operational temperature limit of 85°C. This temperature limit can be exceeded during fast fueling if hydrogen is not precooled. Precooling allows for a dispenser to fuel a vehicle at a faster flow rate by preventing the storage tank on the vehicle from overheating. Fueling protocols and requirements are presented in SAE J2601 Fueling Protocols for Light Duty Gaseous Hydrogen Surface Vehicles [1]. A heavy-duty equivalent is under development with similar requirements for precooling. Currently, conventional precooling for light-duty vehicle refueling uses a heat exchanger and chiller to cool the hydrogen gas to -40°C before entering the vehicle. The precooling system represents a significant part of the station capital and operating costs, so if the cost of the precooling system can be reduced by improving its efficiency, the overall station capital and operating cost can be reduced. In this project, National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) researchers teamed up to investigate the turboexpander precooling application. A turboexpander is a device that places a turbine in a flow path where a pressure differential can be attained. This expansion device will extract work and lower the temperature of the fluid as the pressure reduces. While initial calculations based on established principles showed potential for a turboexpander to generate cooled gas, much work needs to be done to prove the concept. Turboexpanders typically work best under steady state conditions, while the dispenser is a very dynamic flow system. Dynamic turboexpander systems have been proven, such as a turbocharger on a gasoline vehicle. The inlet pressure at a dispenser is also much higher than any other known turboexpander system but should behave similarly to higher density fluids at lower pressures. Having both performed initial calculations, NREL and SNL researchers teamed up to investigate the turboexpander precooling application further. A project was soon built around the idea with SNL performing system modeling using previously proven capabilities and NREL performing hardware characterization with established station capabilities. Creare LLC was contracted as the turbomachinery expert to design and build the concept device. Part way through the project, however, contracting issues with the funding partner caused the project to terminate early before building and characterizing the concept device. While the project could not continue, many key findings were already learned. This paper is a summary of those findings.

08 HYDROGEN↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2022

Bioproducts from high-protein algal biomass: An economic and environmental sustainability review and risk analysis

This report has High-protein algal biomass as an important bio-commodity that has the potential to provide a new source of sustainable protein products. Herein is a critical review that identifies 1) the most relevant sustainability findings related to the processing of proteinaceous algal biomass to higher value protein products and 2) the potential pathways to improve life cycle assessment (LCA) and techno-economic analysis (TEA) metrics, including life-cycle carbon dioxide equivalent (CO 2 eq), life cycle energy, and minimum selling price (MSP) of these products. The critical review of the literature revealed a large variation in model input parameters relating to these metrics. Therefore, a Monte Carlo analysis was conducted to assess the risk associated with these input variations. To understand the uncertainties that propagate into high-protein algae to products’ systems, we reviewed more than 20 state-of-the-art unit operations for algal biomass processing., including cell disruption, protein solubilization, protein precipitation and purification, and protein concentration. We evaluated displacement of proteinaceous products by algal-bioproducts, including ruminant feed, aquaculture feed, protein tablets, and biopolymers and biopolyesters, with prices in the market ranging from 1.9 to 120 $\$$ kg -1 protein. This review realized that the MSP of ruminant and non-ruminant feed ranges from 0.65±0.56 to 2.9±1.1 $\$$ kg -1 protein, and bioplastics’ MSP ranges from 0.97 to 7.0 $\$$ kg -1 protein. Regarding LCA metrics, there is limited research on life cycle energy in proteinaceous biomass concentration and bioproduct systems, reported at 32.7 MJ kg protein -1 , for animal feed displacement. Animal feed emissions in the literature report negative fluxes, representing environmental benefits, as low as -3.7 kgCO 2 eq kg -1 protein and positive fluxes, i.e., global warming potential, as high as 12.8 kgCO 2 eq kg -1 protein. There is limited research on bioplastics life cycle emissions reported at 0.6 kgCO 2 eq kg -1 protein. In general, the studies to date of algae-derived protein bioproducts showed similar life cycle emissions to soybean meals, nylon, polymers, and polystyrenes. Our risk analysis realized that more than 50% of scenarios can result in negative-net life cycle CO 2 eq emissions. This review and risk analysis assess and demonstrate the scenarios that improve economic and environmental sustainability metrics in high-protein algal bioproducts systems.

09 BIOMASS FUELS↗
DOE OSTI2024

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Transforming our understanding of chloroplast-associated genes through comprehensive characterization of protein localizations and protein-protein interactions

Bioenergy crops are a renewable source of fuels and are a critical base for building a carbon-neutral economy. Rational engineering of bioenergy crops has the potential to enhance the yields. However, our ability to engineer plants is limited because the functions of most genes remain unknown. Systematic characterization of gene function in plants thus has the potential to greatly accelerate bioenergy research. Here, we focus on the chloroplast, an underexplored energy-producing organelle that is a hallmark of plants. The chloroplast is one of the promising targets of biofuel crop engineering efforts because of its central role in photosynthesis, metabolism, and intracellular signaling. However, the protein composition of the chloroplast and the functions of most of its proteins remain poorly characterized. At the core of this project, we sought to comprehensively determine the localization of chloroplast-associated proteins and generate a spatially defined protein-protein interaction network for chloroplast. For this purpose, we used the leading model alga Chlamydomonas reinhardtii, which greatly increased experimental speed and throughput. We illustrated the value of our findings to land plants by determining the localization of Arabidopsis thaliana land plant homologs of the Chlamydomonas proteins. Altogether, we were successful in determining the localization of 1,034 chloroplast-associated proteins in Chlamydomonas. The localizations provide numerous insights into the spatial organization of chloroplasts and how they function to support photosynthesis. The localization patterns of distinct proteins revealed new chloroplast structures and revealed new spatial organization inside the chloroplast. We also identified new components of known chloroplast structures, such as the chloroplast envelope, nucleoid, plastoglobuli, and pyrenoid. We identified these new components by investigating the interacting partners of known proteins. Many proteins localized in both the chloroplast and other cellular structures, thereby hinting at new functions and communication between cellular structures. We also applied machine learning on the atlas to generate predictions for the location of all of the proteins in Chlamydomonas. This enabled us to assign putative functions to many uncharacterized proteins based on their cellular location. Altogether, this research establishes a rich resource that opens new avenues of investigation and guides future work in deciphering and manipulating chloroplast function. Next, we developed an extensive protein-protein interaction network for the chloroplast by performing affinity purification-mass spectrometry on ~1,150 tagged chloroplast-associated proteins, the first such large-scale study in any photosynthetic organism. This dataset reveals 4,694 high-confidence protein-protein interactions, offering insights into the functions of thousands of conserved poorly-characterized chloroplast proteins. This systematic identification of protein-protein interactions in the chloroplast also provides multiple exciting new research directions and a detailed blueprint of the chloroplast's operation. This research lays the groundwork to decipher the inner workings of the chloroplast, the cell structure at the heart of photosynthesis. The spatial atlas and protein-protein interactions reveal chloroplast organizational features that would not have been accessible with traditional approaches. The localization mapping, insights into the function, and research materials generated further provide a rich resource for the research community to advance the understanding of how the chloroplast is organized to enable engineering of enhanced photosynthetic organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2023

Technical Evaluation of Accelerated Basin De-Inventory Material Addition to Sludge Batch 11 (Rev. 1)

The Accelerated Basin De-inventory (ABD) program involves discarding spent nuclear fuel that is currently stored in L-Basin to the Defense Waste Processing Facility (DWPF) for vitrification. The first ABD discards will occur during the preparation of Sludge Batch (SB) 11. Savannah River Mission Completion has requested that the Savannah River National Laboratory assess the technical gaps related to the increased gadolinium poisoning requirement and the impacts of performing the Low Temperature Aluminum Dissolution (LTAD) process in Tank 51 with H-Canyon discards present. The following summarizes the evaluation of the impacts of increasing the quantity of gadolinium (and related topics) from what was previously evaluated in the SRNL studies of gadolinium-poisoned ABD material solubility, the overall ABD flowsheet review, and increasing the fissile mass loading in glass: 1) Based on literature surveys, there is no indication that organic interactions with gadolinium will be significant at the high pH (typically >13) conditions of the Concentration, Storage, and Transfer Facilities. Any interactions of gadolinium with organics in DWPF are not expected to adversely impact DWPF or downstream facilities. Thus, there is little-to-no residual risk from organic interactions with gadolinium [Gap closed]; 2) Adding depleted uranium to ABD material, targeting 235 U enrichment of 4.90% within each transfer window, will mitigate potential impacts from an increase in soluble 235 U enrichment during sludge washing and LTAD. The plan to take advantage of previous transfers and allow 235 U enrichment of >5% during the final transfer window carries a risk that Tank 51 supernate will have a 235 U enrichment of >5%, which should be evaluated for acceptance; 3) Increasing the gadolinium mass ratio to 3.0:1 Gd: 235 U(eq SLU ) should lead to the same or higher partitioning of gadolinium into the solid phase within the DWPF Chemical Process Cell, resulting in both liquid and solid phases with expected partitioning of Gd consistent with the prior solubility study [Gap closed for SB11]; 4) There are no expected impacts on DWPF melt temperature and melter operations due to the minimal ~0.2 weight percent (wt%) increase in Gd concentration relative to previous sludge batches [Gap closed for SB11]; 5) As observed previously, Gd is expected to enter the off-gas system via physical entrainment, but at a slightly higher concentration than what was observed for SB9 melter off-gas pluggage deposits (0.07 wt%) [Gap closed for SB11] ; 6) There are no expected impacts on DWPF recycle or the Recycle Collection Tank glycolate destruction process. [Gap closed for SB11]; 7) Gd is projected to be a trace component in the SB11 glass (<0.5 wt%) and can be ignored for process control. Trace components do not significantly impact glass durability, thus the conclusions of the previous Product Consistency Test evaluation at a fissile mass loading of 2,500 g fissile/m3 glass still applies to SB11. The ~0.1 wt% increase in Gd2O3 concentration relative to the previous study will not impact the predictability of SB11 glass with the DWPF Product Composition Control System (PCCS) models for durability or the acceptability of glass according to the Waste Acceptance Product Specifications (WAPS) criterion for product consistency [Gap closed for SB11]; 8) No additional Toxicity Characteristic Leaching Procedure testing is necessary for SB11 and the hazardous waste specification of the SB11 DWPF waste form is unchanged after the addition of the ABD stream [Gap closed for SB11]. The following summarizes the evaluation of the impacts of adding two-thirds of the ABD material to Tank 51 prior to LTAD: 1) The addition of two-thirds of the ABD increases overall aluminum mass from 1.39×10 4 kg to 1.64×10 4 kg (15.5% ABD Al). The form of the insoluble portion of the Al resulting from ABD addition should be the more readily dissolved Al(OH) 3 and amorphous forms. The portion of the ABD aluminum that is processed by LTAD is expected to be completely soluble, thus requiring that less of the boehmite in the sludge be dissolved to reach the same Al target in the SB. [Gap closed for SB11]; The expected LTAD impact on other components, as related primarily to the components in ABD, are discussed. Gd is expected to remain insoluble during LTAD and not impact the solubility of other components. [Gap closed for SB11]; The addition of two-thirds of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 13,100 kg (63% ABD U) and the projected plutonium mass from 86.0 kg to 89.5 kg (3.9% ABD Pu). The addition of all of the ABD increases overall projected SB11 uranium mass from 4,740 kg to 16,100 kg (70% ABD U) and the projected plutonium mass from 86.0 kg to 90.4 kg (5.3% ABD Pu). The 235 U enrichment will be ≤5%. The fissile uranium will be adequately poisoned by Gd and the fissile Pu will be adequately poisoned by Fe from the sludge. [Gap closed for SB11]; There is a low risk that ABD addition will impact the rheology or pumpability of the slurry. There is a low but higher risk of ABD addition prior to LTAD impacting the settling rate; Based on the evaluation of adding two-thirds of the ABD material and all of the ABD material prior to the LTAD process, there is no volume or mass limit that would need to be imposed on ABD additions prior to LTAD. [Gap closed for SB11]. Revision 1 of this report addresses a variation on the ABD additions and LTAD strategy where sodium hydroxide additions for LTAD may be performed intermittently or concurrently with an ABD addition window. The proposed change does not alter the conclusions of this evaluation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Optimized Carbon Fiber Intermediate Development to Enable High-Volume Manufacturing of Lightweight Automotive Composites

The ongoing pursuit of improved fuel economy and reduced greenhouse gas emissions has resulted in sustained interest for lightweight materials technologies. In this context, carbon fiber composites have captured the imagination of automotive engineers due to the potential to achieve substantial mass reduction when compared to traditional steel construction. That stated, the use of carbon fiber composites in automotive has been limited for the most part to premium supercars and other derivative platforms. In these cases, manufacturing costs are less of an obstacle to implementation, and the performance benefits of carbon fiber have enabled production of structures offering more than 50% weight savings. In practice, translating these low-volume demonstrations onto high-volume vehicle platforms has remained challenging. This can be attributed to several factors, with the absence of suitable high throughput production methodologies being a key impediment. To date, structural, crash critical components have relied upon manufacturing techniques born out of the aerospace industry. This has created a disconnect between automotive production systems that are accustomed to manufacturing multiple parts per minute and the aerospace technologies that have cycle times in the order of hours. Consequently, the focus of this project is the development of manufacturing process technology for carbon fiber composites that can support a mainstream vehicle program at an assumed throughput of 100,000 vehicles per year. In practice, this translates to a part-to-part cycle time of less than 3 minutes. Project participation included contributions from a broad range of academic, industrial and national lab partners. The primary scope of work, being the development of new carbon fiber epoxy compounds that are stable at room temperature and suited to high throughput automated processing. For project management, the work streams were divided into six key areas, with the lead organization in parentheses. • Carbon fiber/epoxy materials formulation development and scale up (Dow) • Simulation of discontinuous near isotropic meso-structure intermediates (Purdue) • Simulation of mechanical performance of compression molded components (Purdue) • Meso-Scale morphological analysis and correlation with structural performance (UTK) • Paint and adhesion durability analysis (MSU) • Demonstrator part design, prototype production, and validation testing (Ford). The primary goal at the commencement of the project was development of a chopped carbon fiber sheet molding compound (SMC) that offered a three times improvement in tensile modulus over a comparable glass-based SMC. In addition to meeting mechanical performance targets, the resin kinetics were modified to achieve a processing cycle time of less than 3 minutes. Other critical-to-quality (CTQ) specifications were also stipulated to account for a broad range of materials and processing characteristics. To achieve the above, staff scientists at Dow Chemical created an extensive series of new epoxy blends for testing and validation. Throughout this development, a key challenge was attaining material performance goals without impacting processing behavior and paintability of finished components. The latter required a new internal mold release system being developed by Dow that was designed to complement the kinetics of the rapid cure epoxy. As a complement to work studies at the industrial partners, the teams from academia executed a series of analytical and experimental studies to investigate potential factors influencing CF-SMC performance. Unit cell models were developed to capture the meso-scale representations of the fiber matrix architecture. Results of this analysis and subsequent morphological investigations led to the design of a novel composite derivative comprising carbon fiber platelets embedded in an epoxy matrix; the platelet size and aspect ratio playing significant role in final composite properties. This approach was a departure from previous research in CF-SMC development whereby bulk filamentization or disassembly of the carbon fiber rovings had been considered the most effective means of achieving both fiber wet through and wet-out. As the course of the academia studies progressed, the aspect ratio of the fiber constituents was further optimized before finalizing material attributes and processing conditions. For the purposes of technology validation, the Ford team led a work stream devoted to the design, fabrication and testing of demonstration components. The carbon fiber SMC material has the potential to displace numerous stampings and castings on an automotive structure but ultimately vehicle closure applications were selected to showcase the abilities of the CF-SMC to achieve both mass reduction and business case for large complex structures. Using target properties established by the Dow staff scientists, the complete closure system for a full-size sedan decklid and a mid-size wagon liftgate were engineered. Prototypes for both applications were fabricated using production representative processing methods to allow for physical testing and performance validation of the CF-SMC structures. Following completion of a testing program that concluded with a FMVSS301 55 mph offset rear crash, the CF-SMC formulation was declared by the Ford team to have met all engineering requirements. To summarize, the joint development activities during this project led to significant technical breakthroughs and achievement of all milestones. The result was the development of a novel, tack-free carbon fiber molding compound that is suited to automated processing. This combined room temperature stability, fast cure kinetics, and internal mold release system facilitates cycle times that are conducive to high-volume production. The VORAFUSE M6400 successfully passed technology validation at Ford and is now eligible for consideration on future production commercial vehicle programs.

36 MATERIALS SCIENCE↗