Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial area”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Concentration-Driven Slowdown of Chain-Exchange Kinetics in Block Copolymer Micelle Solutions

Here, we investigate the dependence of chain-exchange kinetics on polymer concentration for spherical micelles prepared from a polystyrene-b-poly(ethylene-alt-propylene) (SEP) diblock copolymer ( M n = 75 kDa, Ð = 1.05, and f PS = 0.18) in the PEP-selective solvent squalane. Solutions with concentrations of 15 vol% or more adopt a body-centered-cubic (BCC) packing of micelles, or a dense liquid-like packing (LLP), as confirmed by small-angle X-ray scattering. Concurrently, the mean aggregation number of the micelles grows steadily from about 60 in dilute solution to nearly 200 by 50 vol%. Time-resolved small-angle neutron scattering (TR-SANS) analysis reveals a monotonic and substantial (i.e., more than four orders of magnitude) decrease in the rate of chain exchange as the polymer concentration increases from 1 to 50 vol%. The TR-SANS relaxation functions are analyzed, taking into account both changes in contrast (due to chain exchange) and evolution in the structure factor (due to refinement of micelle packing). The activation energy for chain pullout is extracted using an established model and can be expressed as a linear function of polymer concentration, indicating that, in addition to the enthalpic barrier for extracting the core block, the evolution of micelle characteristics with increasing concentration imposes further constraints on chain exchange. Consistent with a theory of Halperin, these results suggest that the dominant factor reflects steric penalties during core block pullout arising from a reduced interfacial area per chain and associated increased corona crowding, due to the increase in aggregation number. The rate of exchange is independent of whether the micelles adopt a BCC or LLP arrangement, confirming that the rate-limiting step is extraction of the core block.

chemical structure↗
DOE OSTI2024

Solution-Processable Polymer Donor-Small Molecule Acceptor Bulk Heterojunction Organocatalysts for Enhancing CO 2 -to-CO Photoconversion

Metal-free organic semiconductors have emerged as a kind of promising star material in photocatalysis. However, their photocatalytic efficacy is impeded by the great recombination rate of the photogenerated charge carriers. Herein, we demonstrate the devising of solution-processable organic semiconductor bulk heterojunction (BHJ) systems with varying ratios by using BP3DT and PC 61 BM as the donor (D) and acceptor (A) candidates, respectively, for enhanced selective CO production from CO 2 photoreduction. The prepared BP3DT/PC 61 BM BHJs are solution-processed on both glass substrates (GS) and molecular sieves (MS) via a drop-casting technique to form an active thin film and porous photocatalysts to catalyze CO production. All BP3DT/PC 61 BM BHJs efficiently produce CO as the primary product in the existence of trace-level TEA/H 2 O (0.5 mL). Our results show that the photocatalytic efficiency depended on the D/A w% ratios and its nanoscale phase separation morphology within the BHJ, contributing to the promotion of exciton dissociation. Additionally, the BHJ undergoes an ostensible morphological transition from large-scale phase-separation to a fine mixed homogeneous nanoscale by varying the BP3DT/PC 61 BM blend ratio from 30/70 w% to 70/30 w%, which increases the interfacial area and shortens the distance of the exciton to the D/A surface, and promotes charge carrier separation, corroborated by photo-optical and photoelectrochemical experimental results. Accordingly, the resulting BP3DT/PC 61 BM BHJ@MS with a D/A ratio of 70/30 w% achieves a CO yield of 819 μmol·g cat -1 h -1 , surpassing those of pure PC 61 BM (126 μmol·g cat -1 h -1 ) and BP3DT (72 μmol·g cat -1 h -1 ) by 6.5- to 11.3-fold and also exceeding that of metal-free linear conjugated polymer-based organocatalysts reported to date. Specifically, this is the first work to develop an organic-semiconductor-based BHJ photocatalyst, which may offer a nascent way to design more efficient organocatalysts for CO 2 reduction.

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

Topological control of liquid-metal-dealloyed structures

Abstract The past few years have witnessed the rapid development of liquid metal dealloying to fabricate nano-/meso-scale porous and composite structures with ultra-high interfacial area for diverse materials applications. However, this method currently has two important limitations. First, it produces bicontinuous structures with high-genus topologies for a limited range of alloy compositions. Second, structures have a large ligament size due to substantial coarsening during dealloying at high temperature. Here we demonstrate computationally and experimentally that those limitations can be overcome by adding to the metallic melt an element that promotes high-genus topologies by limiting the leakage of the immiscible element during dealloying. We further interpret this finding by showing that bulk diffusive transport of the immiscible element in the liquid melt strongly influences the evolution of the solid fraction and topology of the structure during dealloying. The results shed light on fundamental differences in liquid metal and electrochemical dealloying and establish a new approach to produce liquid-metal-dealloyed structures with desired size and topologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hierarchical assembly of ZnO nanowire trunks decorated with ZnO nanosheets for lithium ion battery anodes

Hierarchical architectures composed of nanomaterials in different forms are essential to improve the performance of lithium-ion battery (LIB) anodes. Here, we systematically studied the effects of hierarchical ZnO nanostructures on the electrochemical performance of LIBs. ZnO nanowire (NW) trunks were decorated with ZnO NWs or ZnO nanosheets (NSs) by successive hydrothermal synthesis to create hierarchical three-dimensional nanostructures. The branched ZnO NSs on the ZnO NW trunk exhibited a two-fold higher specific gravimetric capacity compared to ZnO NWs and branched ZnO NWs on ZnO NW trunks after 100 cycles of charging–discharging at 0.2C (197.4 mA g –1 ). The improvement in battery anode performance is attributable to the increased interfacial area between the electrodes and electrolyte, and the void space of the branched NSs that facilitates lithium ion transport and volume changes during cycling.

25 ENERGY STORAGE↗
DOE OSTI2025

Interface instabilities in hafnium hydride entrained iron metal matrix composites

The chemical interactions in Fe–HfH 2 metal matrix composites (MMCs) are studied across multiple length scales to elucidate the decomposition of the parent phases and corresponding reaction zone physics during direct current sintering. Fe–HfH 2 composites were synthesized with increasing as-mixed hydride contents of Fe–25% HfH 2 , Fe–40% HfH 2 , Fe–55% HfH 2 , and Fe–70% HfH 2 (all in vol. %) to demonstrate the ability to achieve sintered MMCs with target hydride contents. Samples were probed across multiple length scales through a multi-modal workflow employing x-ray diffraction, scanning electron microscopy and segmentation analysis, and synchrotron techniques including hard x-ray fluorescence mapping and nanoprobe x-ray absorption near-edge structure measurements. Under the selected sintering temperature and pressure conditions, hydrogen evolution is seen to evolve through parallel paths: thermal decomposition from during the transformation of HfH 2 to HfH x<2 and through subsequent reaction with the Fe matrix leading to intermetallic phase formation. Specifically, HfFe and HfFe 2 intermetallic formation accelerates the release of hydrogen with a subsequent HfO 2 phase forming at grain boundaries. For this MMC, the consumption or loss of hydrogen can be considerable in compacts with initial hydride loading of 25%–40% HfH 2 approaching 83% hydrogen loss for the lower volume fraction composites. Increasing the volume fraction of HfH 2 to 70% enhanced the retained hydrogen content to 53% and attributed to the reduced interfacial area intrinsic to the increased HfH 2 loading in this MMC.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Implementation of a Drift Flux Model into SAM with Development of a Verification and Validation Test Suite for Modeling of Noncondensable Gas Mixtures

The advanced thermal-hydraulic system code, System Analysis Module (SAM), was originally developed for the modeling of single-phase flow in advanced reactors. It has since been expanded to include a four-equation drift flux model for the modeling of two-phase flows containing a noncondensable gas. The model was expanded to support the modeling of molten salt reactor (MSR) designs in which the fuel is directly dissolved in the circulating coolant. These designs have shown that circulating gas bubbles can play an important role in the management of fission products and the operational behavior of the reactor. A drift flux model was implemented to more accurately capture the localized behavior of the void in the core and its impact on the mass transfer of fission products. A thorough assessment of the new model was performed by developing a verification and validation test suite. Verification problems were designed to test all major terms in the new governing equations. The new model converged to the correct solution at the expected order of accuracy for all verification cases. The validation cases included a wide range of flow and void conditions in different pipe geometries. Although higher void experiments show a slight underprediction of void by the drift flux model, experiments that aim to reproduce Molten Salt Reactor Experiment (MSRE) experimental conditions show good agreement with the model. The gas transport model was activated for a SAM model of the MSRE to demonstrate that it can be used in a more complex model. Finally, this gas transport model will be used along with an interfacial area transport equation being implemented in SAM for the prediction of mass transport behavior in MSR conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Facet-dependent dispersion and aggregation of aqueous hematite nanoparticles

Nanoparticle aggregates in solution controls surface reactivity and function. Complete dispersion often requires additive sorbents to impart a net repulsive interaction between particles. Facet engineering of nanocrystals offers an alternative approach to produce monodisperse suspensions simply based on facet-specific interaction with solvent molecules. Here, we measure the dispersion/aggregation of three morphologies of hematite (α-Fe 2 O 3 ) nanoparticles in varied aqueous solutions using ex situ electron microscopy and in situ small-angle x-ray scattering. We demonstrate a unique tendency of (104) hematite nanoparticles to maintain a monodisperse state across a wide range of solution conditions not observed with (001)- and (116)-dominated particles. Density functional theory calculations reveal an inert, densely hydrogen-bonded first water layer on the (104) facet that favors interparticle dispersion. Results validate the notion that nanoparticle dispersions can be controlled through morphology for specific solvents, which may help in the development of various nanoparticle applications that rely on their interfacial area to be highly accessible in stable suspensions.

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

Latent Space Simulation for Carbon Capture Design Optimization

The CO2 capture efficiency in solvent-based carbon capture systems (CCSs) critically depends on the gas-solvent interfacial area (IA), making maximization of IA a foundational challenge in CCS design. While the IA associated with a particular CCS design can be estimated via a computational fluid dynamics (CFD) simulation, using CFD to derive the IAs associated with numerous CCS designs is prohibitively costly. However, previous works such as Deep Fluids (Kim et al., 2019) show that large simulation speed-ups and low error are achievable by replacing CFD simulators with neural-network (NN) surrogates that mimic the CFD simulation process. This raises the possibility for a fast, accurate replacement for a CFD simulator, and thus computationally feasible IA-based CCS-design optimization. As such, here, we explore whether an existing NN-surrogate approach (and variants we develop) can successfully be applied to our complex carbon-capture CFD simulations, with the ultimate goal of obtaining a fast and accurate simulator for our CCS-design application. Our experiments build on the Deep Fluids approach and find that resulting surrogates can produce large speed ups (4000x) while maintaining IA relative errors as low as 4% on unseen CCS configurations (interpolating between configurations seen during training). Thus, despite less faithfulness to the underlying physics of the problem, NN surrogates may be a promising tool for our CCS design optimization problem. Notably, though, the Deep Fluids approach has limitations for our application (such as model non-transferability to CCS-packing changes), which we discuss. We conclude with potential directions for future work that may, like innovations we introduced here (e.g., transformer-based dynamics prediction), improve performance on our complicated dataset of CCS CFD simulations.

AI Surrogate Simulation, Carbon Capture, design op↗
DOE OSTI2021

SOC Synthetic Microstructure Bank

QUICK START: Start with property_library.html (can be found by typing the filename into the query box) and use the interactive table to filter, sort, and select a microstructure with the desired properties. Search for the alphabetic code to obtain the corresponding dataset. Full description: This is a bank of 1,970 unique 3-phase electrode microstructure files. When you account for reassigning phase IDs (e.g. declare that 1=Ni and 2=pore, instead of 1=pore and 2=Ni), it actually represents 5,910 unique electrode microstructures, each of which could be considered to be either an air or a fuel electrode (e.g. declare that the phase IDs correspond to pore, Ni, and YSZ; or that they correspond to pore, LSCF, and GDC; or whatever electron-conductor and ion-conductor combination is being studied). The voxel size is 50 nm and each electrode file contains a 4x4 grid of (12.5 micron)^3 sub-volumes. If placed together in a grid, they comprise a 50x50x12.5 micron electrode (note that the interfaces between sub-volumes will be sharp; this can be mitigated via simulating annealing/relaxation). The sub-volumes can also be used individually for a reasonably sized 12.5 micron cubic region-of-interest. A user can simply consider the voxel size to be a different value to rescale the volumes (and all of their morphological features, including particle size) as desired. These microstructures were generated using DREAM3D. The general procedure is outlined in https://doi.org/10.1016/j.jpowsour.2018.03.025 The file names are an alphabetic code having to do with the input parameters used in DREAM3D when they were generated. Most users would be best served by starting with the file property_library.html or property_library_subvols.html (which lists properties for each individual subvolume). These files contain a catalogue of the actual, measured properties of every microstructure in the database. Any combination of property values can be filtered and sorted until a desired electrode is found, at which point the user can find the file corresponding to that alphabetic code. The properties in the catalogue include connected TPB density, and for each phase: phase fraction, average particle size, polydispersity of particle size, tortuosity, and connected pair-wise interfacial area. They also include what fraction of each property is connected through to the interfaces of the volume. If the desired combination of properties is not found at first, remember that the phase IDs can be re-assigned arbitrarily, e.g. swapping 1s and 2s. In fact, the database was generated with this in mind so as not to generate redundant microstructures. If the database does not contain the desired property combinations, try to search for the other possible permutations of those properties with re-assigned phase IDs. Please cite https://doi.org/10.1149/10301.0909ecst for use. Please contact the maintainer, William K. Epting, for additional information or assistance.

3D microstructure↗
DOE OSTI2025

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗
DOE OSTI2025

SOC Microstructural Analyzer

This program was designed to analyze the 3-phase microstructure of the electrodes of a solid oxide fuel cell (SOFC) or electrolysis cell (SOEC), both referred to in combination as a solid oxide cell (SOC). It is agnostic to the exact system, so it could be repurposed to analyze any 3-phase microstructure. This tool directly analyzes segmented voxel-based data that has been segmented into phase IDs (1,2,3). The voxels will be analyzed directly for: - tortuosity factors - triple phase boundaries - 2-phase interfacial areas, using a meshed isosurface - mean diameters of each phase, using an inscribed sphere method - standard deviation of the diameters of each phase, from the same inscribed sphere data - connectivity information Comprehensive information is available in the readme file (within the zipped repository in Markdown language, and also available here as a rendered PDF). Please cite this page / DOI, as well as https://doi.org/10.1111/jace.14775, for usage.

3D microstructure↗
DOE OSTI2019

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗
DOE OSTI2020

Safe, Low-Cost, High Energy-Density, Solid-State Li-Ion Batteries

We developed intrinsically safe, robust, low-cost, high-energy-density solid-state Li-ion batteries, by integrating high conductivity garnet-type solid Li ion electrolytes and several different cathodes in tailored microstructures, fabricated by low-cost supported thin-film ceramic techniques, to overcome the primary technological barriers to next generation battery technology. In this ARPA-E project we employed Li-garnet solid-state electrolytes pioneered by one of the co-investigators (Thangadurai) which has room temperature (RT) conductivity of ~10 -3 S cm -1 . The highly stable garnet SSE allowed the use of S, NMC, and other cathodes with Li metal anodes without stability or flammability concerns. We studied the garnet’s stability when exposed to CO 2 and humidity (with or without an ALD coating) and showed that the material can be processed in ambient air for most steps prior to battery assembly, significantly reducing production costs. We extensively studied and scaled up the synthesis of the powder and explored external suppliers of garnet to determine the best pathway for commercial success. We used low-cost and scalable tapecasting fabrication methods to create porous-dense-porous “trilayer” solid electrolyte microstructures with a ~20 μm separating dense layer, and >70% porosity porous layers to host the electrodes. These scalable multilayer ceramic fabrication techniques, without need for dry rooms or vacuum equipment, also provide for dramatically reduced manufacturing costs. The tailored microstructured electrode supports (scaffold) increased interfacial area ~50x, overcoming the high impedance typical of planar geometry SSLiBs, resulting in record-setting 10 mA/cm 2 lithium cycling in a symmetric cell. We expanded our tapecasting ability and showed tapes >25m in length by the end of the project. We demonstrated the ability to sinter garnet without the use of a garnet powderbed, which could have been a major obstacle and waste of material at scale. At the same time, we expanded the size of our trilayers from coin cells (0.7 cm 2 ) to the full format design (30 cm 2 ). To develop full cells, we demonstrated 250 Wh/kg at room temperature with commercial NMC electrode films in a bilayer configuration. Similarly, we demonstrated 280 Wh/kg trilayer Li-S cells at room temperature with >90% capacity retention over 300 cycles. At 90°C we demonstrated 350 Wh/kg with a Li-S trilayer. The high thermal stability of the garnet electrolyte also enabled extended cycling (400 cycles) of a Li-TiS 2 cell at 150 °C, showcasing the exceptional safety and nonflammability of the solid-state trilayer batteries. Ion Storage Systems has taken this technology and knowledge to produce trilayer garnet electrolytes and is building out facilities to commercialize high performance solid state batteries. ISS has hired a strong staff, identified and communicated with a large number of potential partners, and successfully secured follow-on funding.

25 ENERGY STORAGE↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2021

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Simulation of Transport Phenomena in Bubble Column Reactors

The production of carbon and hydrogen from methane is a viable pathway for monetizing the vast reserves of US natural gas. Liquid (metals and/or salts) bubble column reactors operating at temperatures greater than 1000° C allow the conversion of methane into hydrogen without CO 2 emission. The low-density carbon floats on top of the liquid surface, allowing an easier removal, compared to other reactor concepts for methane pyrolysis. In this current report, we present results from simulations of transport phenomena in a liquid bubble column reactor. The goal is to investigate and understand parameters that are crucial for bubble column reactor design and scale up. The report presents the formulation and numerical methods used in the multiphase direct numerical simulation software, Quilt. This software was validated for the bubbly flow using numerical and physics based test problems. Quilt was then used to perform simulations of the bubble column reactor as a parametric study, varying the bubble injection parameters. The simulations were visualized for qualitative analysis and followed up with a quantitative analysis. The effect of injection parameters on various quantitative measures of the bubble residence time, interfacial area and motion through the column are presented.

03 NATURAL GAS↗
DOE OSTI2021

High Performance PEFC Electrode Structures

Raytheon Technologies Research Center (RTRC) in collaboration with Ion Power Inc. and the University of Arkansas at Little Rock and the FC-PAD consortium executed a project that: (1) developed comprehensive models of oxygen and proton transport in cathode catalyst layers that provided new insights about losses associated with reducing platinum loading to 0.1 mg/cm 2 , (2) designed, fabricated, characterized, and tested cells featuring carbon-supported catalyst layers with low platinum loadings that met efficiency and power density targets set forth by the Hydrogen and Fuel Cells Technology Office, and (3) designed, fabricated, characterized, and tested novel catalyst structures. The modeling efforts produced eight scientific publications that shared new descriptions of oxygen and proton transport in catalyst layers with the broader technical community. The key physical insights from the models are expressed in straightforward algebraic expressions, facilitating easy adoption by other researchers. The models developed during the project explain why resistance to oxygen transport is inversely proportional to interfacial area of platinum – an observation that eluded quantitative explanation for a decade. The models quantitatively predict polarization when measured morphological parameters and transport properties are input. Membrane-electrode assemblies made during the project met stated goals for current at high voltage and power density at rated voltage for transportation fuel cells at the end of budget period 2. High performance was achieved by developing a diffusion layer with low resistance to oxygen transport and a PtCo/C cathode catalyst layer capable of supporting rapid proton and oxygen transport. The project was subsequently directed to curtail experimental work and focus on modeling during budget period 3. Ion Power Inc. developed new expertise in fabricating catalyst layers containing platinum-cobalt catalysts, and processing very thin membranes. The University of Arkansas at Little Rock developed novel nanocolumnar self-supported thin film electrocatalysts that can be grown on carbon supports using a simple high pressure sputter deposition method. The nanocolumnar microstructure provides adequate surface-to-volume ratio for efficient platinum utilization, and the conformal platinum shell with larger crystal grain sizes covering the carbon support surface may eliminate durability issues associated with catalyst dissolution, agglomeration, and carbon corrosion.

08 HYDROGEN↗