Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “segregation”

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 1,189 records · Page 66

Microstructure Changes Observed in the Near-Surface Region of SRF Nb Cavities Cutouts upon Cooling/Heating Cycles Using GI-Synchrotron XRD

We have mapped microstructural changes in the near-surface region of Nb from SRF cavity-cutouts upon thermal cycles in the range from 300 to 30 K using grazing incidence synchrotron X-Ray diffraction (GIXRD). Segregation of secondary phases was observed after the thermal cycle, and their nature has been clarified and discussed in view of previous studies on hydrides formation in SRF bulk Nb cavities. The temperature dependence of the relative population of these formed phases was obtained from GIXRD patterns profile fitting. Both, Nb bulk matrix and the new phases formed after cool-down show specific structural features as thermal contraction/expansion, structural transitions, and Nb lattice variation due to the induced strain by precipitates formation. The information derived from this structural study can explain some phenomena as the dissipation at high accelerating field (i.e. High Field Q Slope, HFQS) in the Nb SRF performance as well as new mechanisms never addressed in previous studies.

43 PARTICLE ACCELERATORS↗

Novel Proppant Logging Technique for Infill Drilling of Unconventional Shale Wells

Summary During the development of an unconventional play, wells are drilled and completed in batches, and depending on the development plans, current and expected energy market trends, as well as other developmental considerations, new wells are drilled and hydraulically fractured later near existing producing laterals. This creates challenges in terms of optimizing resource recovery and reducing interwell communication. A novel approach is proposed that utilizes systematic composite sampling and analysis of drilling mud returns to look for and quantitatively identify sand particles. The workflow involves cleaning, drying, and segregation of samples into sizes of interest to us (size distribution of pumped proppant in offset parent wells). These samples are imaged at a very high resolution and analyzed for grains using characteristic optical imaging properties to classify proppant sand particles using computer vision algorithms. Further analysis, such as elemental compositional analysis, is used to validate the results from the imaging workflow. We present a case study from the Permian Basin, where a new child well was used as a test case to prove this technology at the Hydraulic Fracturing Test Site (HFTS-2) in Delaware Basin. We introduce new proppant parameters that help identify sustained proppant zones vs. localized propped fractures. We have used additional diagnostics and data collected at the test site to validate observations from the proppant log and have successfully interpreted significantly propped vs. unpropped zones. A key finding from this test has been the significant proppant transport distances observed away from parent wells. Observable proppant was found at a lateral distance of approximately 425 m for one set of parent wells and more than 915 m for another set of parent wells. While a major limitation of this technique is the sampling rate, given adequate sampling, the proposed technology represents a systematic and one-of-a-kind interpretation of spatial proppant distribution while drilling infill wells. It provides us with unique opportunities to better understand the current state of the reservoir being targeted, including zones that are likely highly drained relative to others, and how the planned hydraulic fracturing of child wells can be improved.

Energy & Fuels↗

Genomic dissection of anthracnose resistance response in sorghum [Sorghum bicolor (L.) Moench]

Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important grain crop behind maize, wheat, rice, and barley. Today, it is of interest as a source of fermentable sugars for the production of renewable fuels and chemicals, and as a source of biomass for co-firing. The productivity and profitability of sorghum are limited by several biotic constraints, most notably anthracnose caused by the fungal pathogen Colletotrichum sublineolum. The most cost-effective and environmentally benign strategy to control anthracnose is through the incorporation of resistance genes. Over the last three years, our research efforts have been directed to identify new sources of resistance in temperate adapted and tropical germplasm, and to delimited genomic regions associated with the observe anthracnose resistant response. Three biparental mapping populations derived from the resistant lines SC112-14, QL3 and IS18760 were evaluated for anthracnose resistance response in Texas, Georgia, Florida and Puerto Rico. In parallel, three high density recombination maps were constructed and used to identify resistant loci. Anthracnose resistant response in line SC112-14 is controlled by a major locus on chromosome 5. Segregation analysis of 1,500 progenies delimited the resistance locus on chromosome 5 to a 23-kb region harboring three candidate genes, including Sobic.005G17230 identified by GWAS of the sorghum association panel (SAP). The latter gene belongs to a family of genes encoding F-box proteins indicating that this resistance response involved in signaling cascades and transcriptional reprograming, rather than recognition of pathotype-associated molecular patterns. In contrast, anthracnose resistant response in lines QL3 and IS18760 is controlled by multiple small-effect genes. Greenhouse evaluation of a representative subset of the three mapping populations against nine pathotypes found that lines susceptible in the field could be resistant to a single pathotype in the greenhouse. Thus, the activation of a resistance response system by a single pathotype could not provide a broader resistance response against multiple pathotypes. The screening of 1,801 sweet sorghum accessions from the National Plant Germplasm System identified 654 accessions with Brix value larger than 10, which in turn was used to select a subset of 233 accessions for evaluation of anthracnose resistant response. Even though most of the accessions were not completely infected by anthracnose, 28 accessions were completely resistant against pathotypes from Texas, Georgia, Florida and Puerto Rico. Genotyping-by-sequencing analysis of this subset identified 157,843 single nucleotide polymorphisms. Population structure analysis of the subset based on a subset of 2,345 unlinked SNPs found that the genetic diversity could be divided into four populations. The genetic relatedness among accessions within populations suggests most of the resistant germplasm may contain few different resistance sources. These resistance sources present in sweet sorghum germplasm could expedite the development of new resistant sweet sorghum cultivars and hybrids by avoiding time-consuming introgression breeding approaches with non-sweet sorghums serving as donor of the resistance alleles.

59 BASIC BIOLOGICAL SCIENCES↗

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

Reflectin: Protein Driver of Dynamically Tunable Biophotonics; New Paradigm for Tunably Reconfigurable Materials

We discovered that the reflectins, with architecture and sequences unlike those of any other proteins known, exhibit a novel interplay between signal-controlled nucleation and growth of dynamically arrested liquid-liquid phase-segregated assemblies that uniquely enables the continuous and precise fine-tuning of physical properties governing the color and brightness of light reflected from the membrane-enclosed nanostructures that contain them. Our research included high-resolution structural analyses conducted at the DOE-BES user facilities of the Stanford Synchrotron Light Source and LBL National Laboratory’s Molecular Foundry, leveraging their capabilities for time-resolved SAXS and cryo-TEM with our ongoing genetic engineering, biophysical and computational analyses to identify the structural determinants and mechanisms governing the dynamic behaviors of the reflectins. We augmented these methods with EPR analyses that revealed for the first time the existence of a specific pathway of reflectin assembly. Using dynamic light scattering and continuously monitored circular dichroism and Raman microscopy, we elucidated the progression of structural changes underlying the mechanisms of signal-induced condensation, folding, assembly, liquefaction, dynamic arrest, regulation and reversal.

59 BASIC BIOLOGICAL SCIENCES↗

Degradation & Performance Studies of ALD-Stabilized Nano-Composite SOFC Cathodes

This work showed, definitively, for the first time, that zirconia Atomic Layer Deposited (ALD) overcoats improve the stability of Solid Oxide Fuel Cell (SOFC) cathodes by chemically gettering interfacial Sr that has segregated out of lanthanum strontium cobalt iron oxide (LSCF). The optimum ALD ZrO 2 overcoat thickness was found to be ≥5 nm and <10 nm. Gadolinium doped ceria (GDC) pre-infiltration was also shown, for the first time, to improve SOFC cathode stability through the chemical getting of Sr at the LSCF-backbone and LSCF-gas interfaces. Due to the fact that it can be performed easily in air on large-area electrodes, GDC pre-infiltration is likely more commercially viable than ALD ZrO 2 overcoating as an interfacial Sr gettering technique. Combining ZrO 2 overcoating with GDC pre-infiltration was not successful because the GDC nano-particles gettered the ZrO 2 overcoat instead of the interfacial Sr.

30 DIRECT ENERGY CONVERSION↗

Tuning Surface Stoichiometry of SOFC Electrodes at the Molecular and Nano-scale for Enhanced Performance and Durability

This project achieved the following objectives. Different cation segregation behaviors of different common SOFC cathodes, including La 0.6 Sr 0.4 Co0.2Fe 0.8 O 3-δ (LSCF), Sr 0.5 Sm 0.5 O 3-δ (SSC), and PrBa 0.5 Sr 0.5 Co 1.5 Fe 0.5 O5 +δ (PBSCF), were determined. We showed how oxygen partial pressure and gas impurities impact the stability of SOFC cathodes. We developed atomic layer deposition (ALD) coating techniques for electrodes and showed that by introducing different elements and different ALD oxidizers, electrode surface chemistry can be altered, resulting in enhanced oxygen reduction kinetics. In addition, we developed solution infiltration technique to enhance performance and durability of electrodes. Different infiltrates as well as different thermal treatment processes were screened to identify the optimal surface modification process that yields both low impedance and high durability. The modified cathode shows excellent stability and has a low area specific resistance (ASR) of only 0.2 Ωcm 2 at 600 °C after over 2000 hours of operation. Further, we developed ceramic anodes, SrFe(Co,Mo)O 3 (SFCM) and SrFe(Ni,Mo)O 3 (SFNM), and enhanced anode oxidation kinetics by solution infiltration or in situ catalyst exsolution from the ceramic anode surface. The modified anodes showed improved performance in full SOFCs, and the optimized anode modification shows high stability for over 400 hours at 550 °C. Moreover, the modified anode also demonstrates high durability in methane. This work provides fundamental understanding of electrode surface chemistry and demonstrates a simple, facile, cost-effective approach to enhance catalytic activity and durability of SOFC electrodes.

20 FOSSIL-FUELED POWER PLANTS↗

Intermediate-Term Thermal Aging Effect Evaluation for Grade 92 and 316L at The LWR Relevant Temperature

Life extension of the existing nuclear reactors imposes accumulated damages, such as higher fluences and longer period of corrosion, to structural materials, which would result in significant challenges to the traditional reactor materials such as type 304 and 316 stainless steels. Advanced alloys with superior radiation resistance will increase safety margins, design flexibility, and economics for not only the life extension of the existing fleet but also new builds with advanced reactor designs. The Electric Power Research Institute (EPRI) teamed up with Department of Energy (DOE) on the Advanced Radiation Resistant Materials (ARRM) program, aiming to develop and test degradation resistant alloys from current commercial alloy specifications by 2021 to a new advanced alloy with superior degradation resistance in light water reactor (LWR)-relevant environments by 2024. Based on a comprehensive microstructure and property screening performed in Phase-1 of the ARRM program, a total of five alloys, together with 316L and X-750 as references, were down-selected for Phase-2 neutron irradiation studies. Because thermal aging could exert a synergistic effect on neutron irradiation due to the low neutron damage rate on the order of 10 –7 displacements per atom per second (dpa/s), Grade 92 (one of the five down-selected alloys) and two heats of 316L were selected in this task to study the effect of aging at 350°C for 12.5–12.7 kh on microstructure and mechanical properties. The aging time is approximately corresponding to 5 dpa neutron irradiation, which can be served as a reference for the 5-dpa-irradiated samples of the alloys to help understand the independent neutron damage influence on microstructure and mechanical properties. Optical microscopy and scanning electron microscopy were used for microstructural characterization. Hardness, tensile, Charpy impact toughness, and fracture toughness in the ductile regime of the aged samples were examined. The aging of Grade 92 led to the formation of many Laves phase in sizes of ~100–200 nm and resulted in some reduction in hardness and yield/ultimate tensile strength with some increases in uniform and total plastic elongations, which may have helped the enhancement of Charpy impact toughness, e.g., ~4 J increase in upper-shelf energy and 20.5°C reduction in ductile-brittle transition temperature compared with the unaged condition. The fracture toughness of the aged Grade 92 showed decent toughness of ~303 MPa√m (K Jq ) with ~87 tearing modulus at 22°C, which reduced to ~241 MPa√m (~20% reduction) with ~75 tearing modulus at 300°C. The two heats of 316L are differentiated by their amounts of δ-ferrite, i.e., ~1 vol% in heat T1103 and ~4 vol% in heat N5B8. The two heats of 316L were subjected to 15% cold work (CW) as this condition is often used in nuclear reactors to help trapping radiation-induced defects and thus postpone the steady state swelling stage. The aging did not result in any noticeable microstructure changes, except for possible segregations at grain boundaries that need to be further investigated. The aging resulted in slight reduction in hardness of 316L-T1103 and slight increase in hardness of 316L-N5B8, which is consistent with their tensile testing results. The 15%CW might have introduced some inhomogeneity, resulting in large standard deviations in hardness, which were reduced after the aging. The aged 316L followed the same trend as the unaged condition for the strength and elongation results from the tensile testing. Minimum elongations appeared within ~300–600°C, above which the elongations seem to increase. Unlike the negligible or minor changes in microstructure, hardness and tensile results of the two heats of 316L, their fracture toughness showed noticeable difference. The aged 316L-T1103 showed good toughness of ~350 MPa√m with ~80 tearing modulus at 22°C, which reduced to ~271 MPa√m (~23% reduction) with ~83 tearing modulus at 300°C. In contrast, the aged 316L-N5B8 showed lower toughness of ~244 MPa√m with ~53 tearing modulus at 22°C, which reduced to ~180 MPa√m (~26% reduction) with ~31 tearing modulus at 300°C. The preliminary results indicate that the presence of high volume of δ-ferrite would noticeable impair fracture toughness although it may not influence hardness and tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Understanding and Overcoming Water-induced Interfacial Degradation in Si Modules (Final Technical Report)

Moisture ingress is an established issue for photovoltaic module durability. Durability studies probing moisture effects typically evaluate performance losses at the module level, attributing global power losses to the overall humidity condition of the test environment while leaving local module behavior unknown. This project successfully develops an in situ short-wave infrared probe of moisture content in PV modules with demonstrated detection limits better than 100 ppm in EVA-encapsulated Al-BSF and PERC architectures (66 µg/cm3) and applicability across the range of terrestrial conditions. Combining this optical moisture quantification method (water reflectometry detection, WaRD) with peel testing of interfaces, biased photoluminescence imaging of performance, and first principles computation, we correlate module moisture content and cell performance over the course of accelerated damp heat tests. By performing multi-level humidity and thermal accelerated testing over thousands of hours (0, 65 and 85% RH and 25, 65, and 85°C), we are able to attribute effects of individual environmental stresses and the dependence on module architecture. Using first-principles computational chemistry, we find that the decomposition of EVA to acetic acid (Norrish II) is thermodynamically preferred over a pathway to acetaldehyde (Norrish I), elucidating the underlying chemistry of encapsulant degradation. We also find it is likely that water segregates to the stable interfaces of the cell, especially the SiNx/EVA interface. We find that moisture is strongly correlated with reduction in adhesion at the interfaces of the cells, particularly at Ag fingers. We show uniquely that the peel strength required to delaminate the encapsulant from the front of the cell depends not only on the exposure to a high heat and humid environment, but also to the moisture content during the mechanical testing, indicating an important interaction between when mechanical stresses are faced in terms of the moisture content of the module and its durability. We find that the background sheet resistance (that is, the sheet resistance rise not attributable to finger interruptions or cracks) is substantially higher in the backsheet mini-modules vs glass-glass packages when humidity and temperature are faced. This increase in power loss due to the background resistance increase over time in damp heat in glass-backsheet modules resulted in ~3-4% larger decrease in relative PCE compared to glass-glass modules. In glass-backsheet modules, the effect of the moisture dose alone is comparable or greater than that of the combined temperature and water term, suggesting that glass-backsheet modules are overall more susceptible to moisture induced performance loss. Finally, we show that thin film PV modules are amenable to WaRD measurement, based on the ability to detect moisture in POE encapsulants or in the cell stack itself. The ease of WaRD measurement and its few requirements on the bill of materials should enable broad applicability of this technique for studying the effects of moisture in PV.

14 SOLAR ENERGY↗

Informative Background on the Interoperability Requirements in IEEE Std 1547-2018

The intent of this document is to support local authorities governing interconnection requirements in adopting the most current requirements for interconnecting distributed energy resources (DERs) to the electric distribution system, as specified in the Institute of Electrical and Electronics Engineers Standard for Interconnection and Interoperability of Distributed Energy Resources with Associated Electric Power Systems Interfaces (IEEE Std 1547-2018 ). To create interconnection rules that successfully address policy goals, market trends, and technical requirements, electric service regulators need to be well informed of DER integration considerations. Interconnection processes and technical issues often overlap; the two are not easily segregated, and they are often considered and evaluated concurrently within the context of interconnection rulemaking. Accordingly, this document addresses both the technical issues and related process considerations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Guide to Updating Interconnection Rules and Incorporating IEEE Standard 1547

The intent of this paper is to support local authorities in adopting the most current technical requirements for interconnecting distributed energy resources (DERs) to the electric distribution system as specified in the Institute of Electrical and Electronics Engineers Standard for Interconnection and Interoperability of Distributed Energy Resources with Associated Electric Power Systems Interfaces (IEEE Std 1547-2018). To create interconnection rules that successfully address policy goals, market trends, and technical requirements, state regulators need to be well informed of DER integration considerations. Interconnection processes and technical issues often overlap; the two are not easily segregated, and they are often considered and evaluated concurrently within the context of interconnection rulemaking. Accordingly, this paper addresses both the technical issues and related process considerations. A lack of uniform and transparent procedures for addressing interconnection rule changes can result in implementation issues or inefficiencies, such as unclear, lengthy, and complicated interconnection rules that can increase distributed generation "soft costs" (i.e., non-hardware costs). This can further delay the deployment of DERs and jeopardize time-constrained national and state policy goals, such as renewable energy targets, or limit participation in emerging techno-economic trends. Additionally, misconstructed interconnection rules are at risk of becoming mired in stakeholder conflicts. This guide addresses the concerns of electric service regulators from both the process and technical standpoints by presenting a structured, step-by-step approach to developing and updating existing interconnection rules. In this document, the process of developing and updating interconnection rules is subdivided into three steps: (1) determining the context (stakeholders and major drivers); (2) developing the rule, including updating technical requirements; and (3) maintaining and revising the rule over time.

14 SOLAR ENERGY↗

Streaming Generalized Canonical Polyadic Tensor Decompositions

In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.

97 MATHEMATICS AND COMPUTING↗

Developing Accelerated Test Protocols and Tuning Microstructures of the Common Materials to Improve Robustness, Reliability, and Endurance of SOFC Cells

This work is built on our prior work where we found that phase transformation in praseodymium nickelates, e.g. Pr 2 NiO 4 (PNO) and (Pr 1-x Nd x ) 2 NiO 4 (PNNO), can be electrochemically driven, and is substantially faster when compared to thermal annealing studies. The first task aims at an attempt to further accelerate the phase transformation in the oxygen electrode by alternating the current input in the cells, which lead to the development of accelerated test protocols (ATPs). ATPs showed up to 60x faster phase transformation and up to 10x faster performance degradation in (Pr 0.50 Nd 0.50 ) 2 NiO 4 electrodes, when compared to long-term operation under constant current density. Furthermore, the phase stable Nd 2 NiO 4 and (La 0.6 Sr 0.4 )(Co 0.8 Fe 0.2 )O 3 (LSCF6482) electrodes were tested in full cells under ATPs, and showed up to 10x faster performance degradation within 1,100 hours in a comparison with long-term thermal annealing studies and electrochemical operation under constant current density. The second task aims at the quantification of the contributions of cell components to the total impedance of a solid oxide fuel cell (SOFC) using electrochemical impedance spectroscopy (EIS) and distribution of relaxation times (DRT). Specifically, the role of gas composition at both anode and cathode was systematically studied, aiming at deconvoluting, identifying and quantifying the contributions of different electrode processes. This was achieved by first tuning the partial pressure of H 2 at the anode and subsequently varying the partial pressure of O 2 at the cathode. The results suggest that, while DRT offers a viable way of deconvoluting different times distributions, additional attention is needed before assigning a peak to a specific electrode process due to the significant overlap of the contributions from the cathode and the anode. Density function theory studies show that both Pr-vacancies and O-defects play a key role on the activity and stability for nickelates towards oxygen reduction reaction. The resident O-interstitials and oxygen ions in the PrO layer form peroxide (O 2 2- ) nearby Pr vacancies. The O 2 2- limits oxygen-ion transport due to the required additional energy to break its O-O bond. We further calculated the formation and segregation energies for different Ln ions (La, Pr, Nd, Pm, Sm, Gd, Tb, Dy, and Ho) in PNO and CeO 2 (111) surfaces. In addition to Nd, Pm and La are suggested as potential dopants in PNO to enhance it stability without decomposition due to their more negative formation energies, lower diffusion energies, and positive separation energies.

01 COAL, LIGNITE, AND PEAT↗

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗

Predictive Tools for Customizing Heat Treatment of Additively Manufactured Aerospace Components

Laser-bed powder fusion (LBPF) additive manufacturing is increasingly being used to produce components of complex geometries using the Ni-base superalloy Inconel 718. The composition and the microstructure of the alloy are currently well optimized for wrought components made using conventional manufacturing processes such as rolling, forging, extrusion, etc. The attractive mechanical properties of the alloy result from the underlying austenitic matrix with fine equiaxed grains, and a high density and uniform distribution of the precipitation hardening phase, γ". Heat treatment steps such as homogenization, solutioning and aging are well documented for the wrought alloy. However, when the same wrought alloy compositions are used for the additive manufacturing (AM) processes, the asprocessed microstructure is significantly different, because of the different thermal history associated with LBPF, including rapid solidification and multiple temperature excursions that lead to multiple re-melting and reheating in the solid state. Rapid solidification introduces potential non-equilibrium effects at the moving solid-liquid interfaces that impact the extent of solute segregation, as well as the morphology of the dendritic grains that form. In order to recover the target mechanical properties, AM components have to undergo post-process heat treatments. However, such heat treatments have to be custom designed for the AM process and the component geometry because of the expected vast differences in the microstructure at various locations of a component with complex geometry. The homogenization and precipitation steps should be optimized for the component so that target mechanical properties can be obtained throughout the part. The objective of this research is to utilize High Performance Computing in phase field simulations of microstructure evolution during post-processing of AM components. The physics-based modeling will be beneficial in reducing the experimental effort required for heat treatment process selection, optimization, and certification, thus leading to a significant reduction in energy consumption for AM and post-processing heat treatment. The optimization study will help identify heat treatments steps that are critical for development of a final desired microstructure with the minimum energy input. This combined with shortening of the production cycle (time-to-market) by reducing the number of failed parts (property targets), and reduction in the number of iterations for process optimization, will enable 30-40% savings in the energy costs. Phase field simulations of the degree of homogenization and the effect of local matrix composition on the nucleation and growth of competing precipitating phases were performed using the Microstructure Evolution Using Massively Parallel Phase Field Simulations code developed in-house at the Oak Ridge National Laboratory. The simulations were able to successfully capture the kinetics of nucleation and growth, and morphologies of various precipitating phases as a function of local matrix compositions and composition gradients characteristic of local microstructures arising from location-dependent variations in the thermal conditions. Future work will involve extending the simulations to a length scale consisting of multiple dendrites, so that the effect of homogenization on the coarsening of the dendrites can be simulated and used as an additional input to the optimization of the heat treatment process.

36 MATERIALS SCIENCE↗

Cost of Capturing CO 2 from Industrial Sources

This systems analysis by the National Energy Technology Laboratory's Strategic Systems Analysis and Engineering directorate) evaluates the cost and performance impacts of capturing CO 2 emissions from nine industrial sources (ammonia, ethylene oxide, and ethanol production, natural gas process, coal- and gas-to-liquids, refinery hydrogen production, iron and steel, and cement manufacturing). The industrial sectors examined are segregated according to the CO 2 purity level of the flue gas stream, prior to treatment. Certain sectors naturally produce a gas stream that is inherently high in CO 2 purity, and these sectors can achieve 99-100% removal. Other sectors produce a lower purity CO 2 flue gas stream and achieving 90-99% removal requires deeper levels of treatment, adding cost. In addition to the report that documents the analysis, a Carbon Capture Retrofit Database tool was also created that allows users to apply CO 2 capture to selected industries, to evaluate the cost of capture and compare across multiple plants, as well as across different industries. Users have the ability to change select input parameters (such as fuel price, capture rate, and financing assumptions) to evaluate the impact on industrial CO 2 capture economics.

20 FOSSIL-FUELED POWER PLANTS↗

Development of Integrated Thermal Fluids Modeling Capability for MSRs

The DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad class of light-water and advanced reactor concepts. The research and development approach under thermal fluids technical area synergistically combines three length and time scales in a hierarchal multi-scale approach. To demonstrate the feasibility and capabilities of a multi-scale thermal fluids capability using these codes, a key joint effort has been pursued to develop an integrated system-and engineering-scale thermal fluids analysis capability with the MOOSE-based codes, through integration of SAM and Pronghorn, both based on the MOOSE framework. This report summarizes the progress in developing an integrated system- and engineering-scale thermal fluids analysis capability based on SAM and Pronghorn for molten salt reactors (MSRs), which gained significant interest in recently years. Two coupling approaches were studied, i.e. separate domain or domain-segregated coupling approach and the domain-overlapping approach. A series of coupled multi-physics models have been developed for a common reference molten salt fast reactor (MSFR) concept, ranging from standalone SAM system model to integrated SAM-Pronghorn-Griffin models. Both the steady state and the transient simulations are performed to the state-of-the-art simulation capabilities of NEAMS software suite in MSFR system applications. This report also covers further development and testing of the gas transport model in SAM for MSR modeling support. The presence of noncondensable gas in MSR systems would have strong impacts on fission gas removal and transport of noble metals throughout the system. Fission products removed through the off-gas system can also impact reactivity and can act as an additional point of heat removal. To ensure that the gas transport model is adequately tested, the supported modeling capabilities and features of SAM were identified, and a suite of tests were developed to test the model for each identified feature. Validation and UQ testing were also performed for the model and demonstrated that the buoyancy term, which was originally developed using non-salt/helium gas experimental data, can capture the experimental gas bubble velocity and diameter to within experimental and code uncertainty. An existing MSRE model was also modified in this work to include the gas transport model and to demonstrate the gas model behavior for realistic conditions of interest.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Influence of Environment on Post-Detonation Chemistry and Debris Formation (Abbreviated Final Report: 20-SI-006)

Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.

07 ISOTOPE AND RADIATION SOURCES↗