Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow field”

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

Filamentary surface plasma discharge flow length and time scales

Nanosecond surface dielectric barrier discharges (ns-SDBDs) are a class of plasma actuators that utilize a high-voltage pulse of nanosecond duration between two surface-mounted electrodes to create an electrical breakdown of air, along with rapid heating. These actuators usually produce multiple filaments when operated at high pulse frequencies, and the rapid heating leads to the formation of shock waves and complex flow fields. In this work we replicate a single filament of the ns-SDBDs and characterize the induced flow using velocity measurements from particle image velocimetry and density measurements from background-oriented schlieren. The discharge is produced by a high voltage electrical pulse between two copper electrodes on an acrylic base. A hot gas kernel characterizes the flow field formed close to the electrodes that expands and cools over time and a vortex ring that propagates away from the surface while entraining cold ambient fluid. The gas density deficit inside the kernel displays a power-law decay over time. Based on the observations, we develop a simplified theoretical model based on vortex-driven cooling and perform a scaling analysis to obtain the induced flow length and time scales. The results show that the cooling process's time scales correspond to a circulation-based time scale of the vortex ring, and the length scale of the kernel corresponds to the vortex ring radius. Furthermore, these findings can guide the choice of optimal filament spacing and pulse frequencies in the design, deployment, and operation of ns-SDBDs for flow control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗
DOE OSTI2026

Impact of Porous Transport Layer In-Plane Conduction on Spatially Resolved Current and EIS Measurements in a Proton Exchange Membrane Water Electrolyzer

An XY segmented cell was developed for low temperature PEM water electrolysis (PEMWE). The system can assess the local performance by enabling in situ measurements of spatial currents and impedances. In this work, we show through experiments, as well as through modelling work, that the porous transport layer (PTL) must be segmented to eliminate crosstalk. Accurate measurements are only possible when crosstalk is fully eliminated. The XY segmented cell is applied to a case study characterizing the impact of a PTL platinum coating void on spatial performance. The localized performance impact of the coating void is found to be orientation specific: coating voids facing the catalyst layer reduce performance significantly more than coating voids facing the flow field. The results suggest that the tolerances for PTL coating uniformity can be lower at the side facing the flow field. The work showcases the feasibility of the XY segmented cell for impact assessment studies. The presented XY segmented cell enables the characterization of spatial phenomena in PEMWE devices and is envisioned to support modeling efforts and the investigation of manufacturing related tolerances for mass produced PEMWE devices.

08 HYDROGEN↗
DOE OSTI2025

Particle removal from a flat surface using a translating bounded vortex flow

A bounded vortex flow is a hydrodynamic approach for removal of particles from a surface without scattering the particles onto nearby surfaces. The bounded vortex flow field is generated by a nozzle that combines azimuthally tilted jets arranged in a circular pattern and a central suction port. When the nozzle face is directed toward an ‘impingement surface’, the flow develops a wall-normal intake vortex below the suction outlet, which causes high shear stress on the impingement surface. When particles are present on the impingement surface, the high shear stress causes particles to roll along the surface and to be lifted off the surface and transported up the core of the wall-normal vortex into the suction outlet. In typical applications, the nozzle would be translated along the impingement surface to clean particles from the surface. The current paper reports on an experimental study of the effect of nozzle translation on the effectiveness of the bounded vortex flow field for particle mitigation. The effectiveness of particle mitigation was examined as a function of flow rate through the nozzle, particle size, and nozzle translation velocity relative to the impingement surface. As a result, numerical computations are used to relate the flow rate to the maximum shear stress on the impingement surface, which is then used to theoretically predict onset of particle motion.

42 ENGINEERING↗
DOE OSTI2024

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗
DOE OSTI2022

Motion tomography via occupation kernels

The goal of motion tomography is to recover a description of a vector flow field using measurements along the trajectory of a sensing unit. In this paper, we develop a predictor corrector algorithm designed to recover vector flow fields from trajectory data with the use of occupation kernels developed by Rosenfeld et al. [9,10]. Specifically, we use the occupation kernels as an adaptive basis; that is, the trajectories defining our occupation kernels are iteratively updated to improve the estimation in the next stage. Initial estimates are established, then under mild assumptions, such as relatively straight trajectories, convergence is proven using the Contraction Mapping Theorem. We then compare the developed method with the established method by Chang et al. [5] by defining a set of error metrics. Here, we found that for simulated data, where a ground truth is available, our method offers a marked improvement over [5]. For a real-world example, where ground truth is not available, our results are similar results to the established method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Comparative analysis of characteristics of a supersonic cw chemical HF laser on molecular fluorine and nitrogen trifluoride

We have experimentally studied the characteristics of a supersonic cw chemical HF laser with a flat nozzle array based on a nozzle – nozzle reagent mixing scheme operating with the use of F{sub 2} – D{sub 2} – He and NF{sub 3} – D{sub 2} – He fuel compositions in an atomic fluorine generator. The active medium flow field pattern, its gas-dynamic characteristics, and the laser radiation spectrum are considered. The comparison of these characteristics has made it possible to establish a strong impact on the flow field (lasing zone length) of the chemical composition of the combustion products from the atomic fluorine generator. The gas-dynamic parameters of the active medium formed by using a NF{sub 3} – D{sub 2} – He fuel have turned out to be less acceptable in terms of the gas flow pressure recovery than in the case of F{sub 2} – D{sub 2} – He fuel, while the radiation spectrum was more preferable from the viewpoint of passing through the atmosphere. (paper)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Localized Flow Modification to Increase Power Capture of a Small-Scale Floating Undershot Waterwheel

The goal of the work described in this thesis is to design a flow augmentation device to increase the power capture and efficiency of a small-scale floating Under-Shot Water Wheel (USWW) currently being developed by Florida Atlantic University research funded by the U.S Department of Energy. The flow concentrator subsystem is intended to maximize the kinetic energy extracted by the marine hydrokinetic (MHK) energy collection device through modification of the local flow field across the capture plane. The primary objective is to increase the velocity and/or rate of mass inflow through the turbine through inserting a streamlined body in the region of interest. By utilizing the resulting flow field to increase hydraulic forcing on the waterwheel blades, the torque and/or RPM of the USWW can be increased. Based on experimental testing in the FAU wave tank at 1:5 prototype scale (280 mm wheel diameter) the flow concentrator was shown to produce an increase in device power coefficient of 17-55% measured over a velocity range of 0.16-0.45 m/s.

Hess, Sullivan↗
DOE OSTI2021

Particle image velocimetry analysis with simultaneous uncertainty quantification using Bayesian neural networks

Particle image velocimetry (PIV) is an effective tool in experimental fluid mechanics for extracting flow fields from images. Recently, convolutional neural networks (CNNs) have been used to perform PIV analysis with accuracy on par with classical methods. Here we extend the use of CNNs to analyze PIV data while providing simultaneous uncertainty quantification on the inferred flow field. The method we apply in this paper is a Bayesian convolutional neural network (BCNN) which learns distributions of the CNN weights through variational Bayes. In order to demonstrate the utility of BCNNs for the PIV task, we compare the performance of three distinct BCNN models with simple architectures. The first network estimates flow velocity from image interrogation regions only. Our second model learns to infer velocity from both the image interrogation regions and interrogation region cross-correlation maps. Finally, our best performing network infers velocities from interrogation region cross-correlation maps only. We find that BCNNs using interrogation region cross-correlation maps as inputs perform better than those using interrogation windows only as inputs and discuss reasons why this may be the case. Additionally, we test the best performing BCNN on a full synthetic test image pair and a real image pair from the 1st International PIV Challenge. We show that ~98% of true particle displacements from the full synthetic image pair can be captured within the BCNN's 95% confidence intervals, and that the BCNN's performance on the real image pair is quantitatively similar to that of algorithms tested in the 1st International PIV Challenge. Finally, we show that BCNNs can be generalized to be used with multi-pass PIV algorithms with a moderate loss in accuracy, which may be overcome by future work on finetuning and training schemes. So to our knowledge, this is the first use of Bayesian neural networks to perform PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This document details the application of the Plateau-to-River (P2R) Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the CA. The simulated flow field will support the simulation of fate and transport of contaminants for use in predicting dose as part of the CA. This document describes the changes made to model inputs in order to simulate the predictive flow field for the CA.

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

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Comparison of three DWM-based wake models at above-rated wind speeds

In this study we investigate three mid-fidelity wind turbine wake models based on the dynamic wake meandering (DWM) model principle, and compare their performance with a reference dataset, produced with large-eddy simulations using the actuator line model. The models are compared with respect to flow field, power, and loads on a row of four 5MW reference turbines experiencing above-rated wind conditions. In general, the DWM models show fairly good agreement with large-eddy simulation for the time-averaged flow fields, blade forces and power, with increasing differences along the turbine row. Also when comparing fatigue loads of blade root moments, the differences between the models increase further into the row, with deviations up to 25 % of the reference case. However, while the development in blade root moment fatigue along the turbine row is predominantly driven by the energy content at the frequency corresponding to the turbine's rotational period (1 P ) for the DWM models, the large-eddy simulation results suggest that the key drivers for the blade root and tower loads are the increase in meandering and energy at higher frequencies (> 1 P ) deeper into the turbine row. For the tower loads, the DWM models highly underestimate the fatigue for the waked turbines. From these results, we suggest priorities for future model developments so that robust model implementations can be used in wind farm design and operation.

17 WIND ENERGY↗
DOE OSTI2022

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
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↗