SEARCH · Engineering Papers

Results for “fluid bubbles”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) using COBRA-TF

With the resurgence of interest in Molten Salt Reactor (MSR) technology among commercial and governmental agencies, there is a need for new experiments and new modelling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory. It is designed to be a dynamic fuel-bearing salt loop placed in the Advanced Test Reactor. One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this can have on fuel-salt properties, fission product transport and interactions, source term generation, corrosion and redox potential control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system and noble metals (Mo, Tc, Te etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in the general species transport capability of COBRA-TF (CTF) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Visualization of R410A Flow Boiling Inside Horizontal Smooth Tubes under Diabatic Conditions

This article presents a novel approach to visualize the flow boiling in round tubes under “diabatic conditions”. In the visualization test section, a 3D printed round tube is placed inside a glass tube with secondary fluid at a higher temperature flowing in the annulus between the two tubes thus providing diabatic conditions. Experimental work was conducted for evaporating R410A flow in a horizontal smooth round tube with the inner diameter of 6.3 mm. Compared to the traditional adiabatic visualization technique, the more realistic evaporative condition is simulated and the developing evaporation process is captured. New insights for the flow behavior at low vapor quality are brought by this diabatic method. The vapor plug occurs at very low vapor quality, and this is attributed to the coalescence of bubbles generated on the heated tube wall. For the intermittent flow, liquid slug triggered by Kelvin–Helmholtz instability affects the nucleation on the tube wall. It brings the fresh liquid to rewet the dry surface and activate the nucleation sites. Finally, heat flux condition affects the heat transfer mechanisms during the flow boiling, and the heat transfer through the nucleate boiling mechanism increases as the heat flux increases.

3D printed↗

DOE OSTI2020

Theory of ignition and burn propagation in inertial fusion implosions

A detailed analytic model is presented here to investigate the physics of burn propagation in inertially confined plasmas. The onset of ignition and burn propagation occurs when alpha heating of the hot spot causes rapid ablation of shell mass into the hot spot. This allows large energy gains to be achieved since most of the fuel mass is located in the shell. Here, we first present a comprehensive review of previous analytic models that have been used to describe the physics of hot-spot evolution and ignition; we then show that a proper description of a propagating burn wave requires a comprehensive model of hot spot and shell evolution that includes proper mass conservation in the shell, fusion reactivity, and fuel depletion. The analytic theory is in good agreement with detailed radiation-hydrodynamic simulations that predict the onset of burn propagation as occurring when the yield enhancement caused by alpha heating is between 15- and 25-fold, $f_α$ ~ 1.4, where $f_α$ = alpha energy deposited/hot-spot energy at bang time, and the hot-spot burnup fraction is approximately 2%. We show that the definition of ignition is not sensitive to the alpha-particle stopping power nor asymmetries provided that the absorbed fraction of alpha particles $θ_α$ is correctly accounted for. Finally, we use the results of 2-D simulations to show that even when $θ_α$ is small and unknown (as is true in hot spots with mid modes that have significant leakage of alpha particles into the surrounding cold bubbles), one can still relate the experimentally measureable parameter $\chi^{53}_α$ to the yield amplification and the burning-plasma parameter $Q^{\text{hs}}_α$ = alpha energy deposited/total input work delivered to the hot spot

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Amphiphilic Block Copolymers for Flocculation and Hydrophobization of Legacy Waste Suspensions in Flotation Driven Dewatering Operations

There exists a degrading legacy fuel pond at the Sellafield site (UK), known as the First-Generation Magnox Storage Pond (FGMSP). After being placed into a passive care and maintenance regime resulting in a long storage period, the fuel rods, primarily their cladding, corroded in the pond forming magnesium hydroxide based sludges and suspended material. These ponds are a grave concern to the British government in terms of hazard and risk reduction. Therefore, decommissioning of these ponds is a top priority beginning with the contaminant retrievals process [1]. The pond has accumulated significant quantities of waste materials amongst the skips of fuel, including but not limited to: large inventories of corroded Magnox sludge, fuel rod fragments, metal fragments (from fuel skips), concrete degradation products (from the pond infrastructure), wind-blown sand, and other materials such as bird guano and animal remains. This challenge requires a chemically robust technology to complete sludge retrievals [2]. As nuclear is different, a more stringent process operation criteria is required. From the criteria, flotation was selected as a viable technology. Flotation involves the application of collector molecules which modify the hydrophobicity of suspended particles allowing them to adsorb to rising bubbles. As some particles lack inertia for flotation due to their size, dual flocculant/collector agents can be deployed for greater particle recovery, in this case, amphiphilic diblock copolymers. Two copolymers of Poly(acrylic acid)-b-poly(n-butyl acrylate) (or PAA-b-PnBA), of different hydrophobic chain lengths (PBA) were synthesized for flotation campaigns. Flotation has shown promise to be a valuable rapid dewatering strategy for decommissioning of legacy waste ponds whilst upholding the required operational criteria. This research has shown: PAA-b-PnBA copolymers promote flocculation of Mg(OH){sub 2} particulates. Diblock copolymers retain more fluid than traditional surfactant based collectors wrt. relative particulate recovery performance- where longer PnBA chain length polymers performed best. Performance appears to be hydrodynamically limited due to over-flocculation of Mg(OH){sub 2}. Potential to combine with sedimentation for high particle recovery with low fluid carryover.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Investigating biomass composition and size effects on fast pyrolysis using global sensitivity analysis and CFD simulations

It is notoriously difficult to build an accurate universal model for biomass pyrolysis due to its sensitivity to a wide number of critical material attributes such as chemical species and physical sizes. In this work, a biomass pyrolysis kinetics with 32 heterogeneous reactions and 59 species was implemented in an open-source multiphase computational fluid dynamics (CFD) software MFiX and validated against two different experimental pyrolysis data sets that provided detailed data describing chemical component yields. The reaction scheme was then used to build a surrogate model and assess the sensitivity of pyrolysis yields to feedstock compositions. The sensitivity analysis determined that the yield of bio-char showed a strong positive sensitivity to the carbon-rich lignin and tannin pseudo-species in the reaction scheme while the bio-oil and bio-gas were correlated to oxygen-rich lignin pseudo-species. The reaction scheme was then integrated into a coarse-grained discrete element model to simulate fast pyrolysis in a bubbling fluidized bed over a range of feedstock particle sizes. The reactor simulations showed further sensitivity to particle size and hydrodynamics. Notably, particles under 0.5 mm have small heat transfer limitations but left the reactor before completely converting and thus reduced the bio-oil yield. Results from this study can be used to guide future development of highly accurate models for fast pyrolysis reactors with a variety of feedstock properties and operating conditions.

09 BIOMASS FUELS↗

DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

DOE OSTI2022

Effects of variable deceleration periods on Rayleigh-Taylor instability with acceleration reversals

The dynamics of an interfacial flow that is initially Rayleigh-Taylor unstable but becomes statically stable for some intermediate period due to the reversal of the externally imposed acceleration field is studied. We discuss scenarios that consider both single and double-acceleration reversals. The accel-decel (AD) case consists of a single reversal imposed at an instant after the constant acceleration instability has entered a self-similar regime. The layer of mixed fluid ceases to grow upon acceleration reversal, and the dominant mechanics are due to internal wave oscillations. Variation of mass flux and the Reynolds stress anisotropy is observed due to the action of the internal waves. Here, a second reversal of the AD case that is termed as accel-decel-accel, ADA is then explored; the response of the mixing layer is shown to depend strongly on the duration and the periodicity of the Reynolds stress anisotropy of the mixing layer during the deceleration period. We explore the effect of this variable deceleration period after the second acceleration reversal where the flow once again becomes Rayleigh-Taylor unstable based on metrics that include the integral mixing-layer width, bubble and spike amplitudes, mass flux, Reynolds stress anisotropy tensor, and the molecular mixing parameter.

42 ENGINEERING↗

DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Experimental and Numerical Investigation of Flow Boiling in Additive Manufactured Foam Structures With Vapor Pathways

Abstract The unique properties of metal foams make them potential candidates for a range of applications, including microsystem thermal management. Using additive manufacturing to create foam-type structures can improve upon prior thermal solutions by eliminating thermal interface materials and allowing for customization/local control of parameters. In the present investigation, flow boiling in additive-manufactured metal foams is investigated both experimentally and numerically. Two test samples, one with uniform structure and the other with pathways for vapor removal, are compared both experimentally and numerically. A conjugate computational fluid dynamics and heat transfer (CFD-HT) model utilizing a three-dimensional volume of fluid (VOF) model with accompanying evaporation/condensation model provided in-depth visualization of the boiling flow phenomena. The experiments generated the thermohydraulic performance over a range of heat fluxes, demonstrating that the sample incorporating dedicated vapor pathways performed better in both pressure and heat transfer performance metrics compared to the uniform foam. Additionally, negative system-level effects (i.e., hydraulic oscillations) were shown to be abated using the vapor removal structures. The numerical model yielded further insight into the factors contributing to the improved performance. Results indicated the pathways functioned as vapor removal channels, allowing the generated vapor to vent from the foam structure into the lanes. Further computational investigations demonstrated changes in flow regimes, where the addition of vapor channels caused the flow to change from churn to annular. Bubble behavior unique to the vapor pathway structure was studied, showing stagnant regions that eject vapor into the channel.

Engineering↗

DOE OSTI2024

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Mitigation of Boiling-Induced Thermal Degradation Using Microporous Nickel Inverse Opals Structures

Engineered microporous structures have received much attention in high-heat-flux electronics cooling due to their high thermal conductivity and permeability, and large surface area for heat transfer, but are susceptible to boiling-induced thermal degradation. Here, this study investigates the efficacy of nickel inverse opals (NiIOs) in mitigating structural degradation caused by corrosion-assisted erosion during pool boiling with water as the working fluid. First, we compared the reliability of NiIOs to copper inverse opals (CuIOs) for a 3-day pool boiling test at constant heat flux. The NiIOs demonstrated superior resistance to thermal degradation due to their inherent corrosion resistance and mechanical strength. Subsequently, we conducted a more controlled experiment to show the effect of heat flux on the degradation of the NiIOs while excluding the effect of temperature variations. Pool boiling tests of 20-μm-thickness NiIOs covering an area of ∼11 × 11 mm 2 with a 2.5 × 2.5 mm 2 heater at the center were conducted at heat flux levels of 20%, 40%, and 60% of the critical heat flux (CHF) for 3 days. The NiIOs subjected to heat flux levels of 20% and 40% CHF showed minimal degradation, while the sample subjected to 60% CHF showed erosion on the top surface due to higher bubble formation and departure rate. These results show the potential of NiIOs as a promising solution for long-term thermal management in high-power electronic devices, although design considerations for maximum allowable heat flux are necessary for reliable operation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Demonstration of a multi-channel fluidized bed particle–supercritical carbon dioxide heat exchanger for concentrating solar applications

High-temperature thermal energy storage in oxide particles at temperatures above 600°C can couple concentrated solar energy with high-efficiency thermal power cycles to provide dispatchable solar-driven electricity. Challenges remain in developing cost-effective primary heat exchangers, which require expensive alloys, to extract the high-temperature thermal energy from the particles to power cycle fluids, such as supercritical CO 2 (sCO 2 ) in recuperated Brayton cycles. To explore one pathway for cost-effective, high-temperature particle heat exchangers, the current study demonstrates a shell-and-plate, particle–sCO 2 heat exchanger with narrow- channel fluidized beds coupled with micro-channel sCO 2 flows in the heat exchanger walls. This study evaluates the feasibility of multiple parallel, narrow-channel fluidized beds in shell-and-plate particle–sCO 2 HXs, to achieve high bed-wall heat fluxes at elevated temperatures. A reduced-order model simulates the narrow- channel, fluidized-bed particle–sCO 2 heat exchanger to design the fluidized bed geometry, in terms of depth, height, and number of channels,for a nominal 40-kWth heat exchanger at particle and sCO 2 inlet temperatures up to 600 °C and 400 °C respectively. The resulting shell-and-plate heat exchanger design operates with bubbling fluidization of the downward-flowing oxide particles to enhance bed-wall heat transfer. The heat exchanger core is fabricated with etched sCO 2 micro-channels in thin wall plates that are diffusion bonded to spacer frames to form the shell-and-plate structure with 12 parallel, fluidized bed channels, 10.4 mm deep. The heat exchanger is tested at the National Solar Thermal Test Facility at Sandia National Laboratories with CARBOBEAD HSP particles at design particle flow rates of 0.20 kg s –1 and inlet temperatures up to 525 °C. Results show that fluidization across multiple parallel channel beds can maintain uniform particle inventory with a common freeboard zone above the heat exchanger core. Bubbling fluidization improves particle–wall heat transfer coefficients but also increases axial dispersion of particle thermal energy, which lowers the log- mean temperature difference such that total heat transfer remains relatively constant to within ±10% over a broad range of fluidization gas velocities. The axial dispersion required particle and sCO 2 flow rates to be increased by 25% over model-designed conditions to achieve the targeted 40 kWth, which indicates the importance of incorporating axial dispersion into heat exchanger design models and of deploying bed structures to suppress it. Furthermore, this study demonstrates the feasibility and preferred fluidizing gas conditions for particle heat exchangers for releasing high-temperature thermal energy storage systems.

14 SOLAR ENERGY↗

DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗