Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-laden”

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 37 records · Page 2

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling Run Data

Repository for the data generated as part of the 2023-2024 ALCC project "Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling." The data was generated with MFIX-Exa's CFD-DEM model. The problem of interest is gravity driven, particle-laden, gas-solid flow in a triply-periodic domain of length 2048 particle diameters with an aspect ratio of 4. The mean particle concentration ranges from 1% to 40% and the Archimedes number ranges from 18 to 90. The particle-to-fluid density ratio, particle-particle restitution and friction coefficients and domain aspect ratio are held constant at values of 1000, 0.9, 0.25 and 4, respectively. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 using NERSC award ALCC-ERCAP0025948.

AMReX↗

Evaluation of the Barracuda Software Package for Simulating Bubble Motion in Vibrating Liquid-Filled Containers

The commercial software package Barracuda, developed by CPFD Software for simulating particle-laden fluid flows, is evaluated as a means to simulate the motion of bubbles in vibrating liquid-filled containers. Demonstration simulations of bubbles rising due to buoyancy forces in a cylinder filled with silicone oil and angled at 0, 30, 45, and 60 degrees from the vertical were performed by CPFD Software. The results of these simulations are discussed, and the capabilities of Barracuda for simulating bubble motion are assessed. It was determined that at present Barracuda does not meet the needs of the desired application. Further developments that would enable its use for this application are highlighted.

42 ENGINEERING↗

Tools for Design and Scale-Up of Solar Thermochemical Reactors: Cooperative Research and Development Final Report, CRADA Number CRD-13-00530

NREL will be collaborating with the Participant on a United States - Australia Solar Energy Collaboration (USASEC) Project Number 1-USO034 "Tools for design and scale-up of solar thermochemical reactors." The grant funds for the Participant's 3.5 year project number 1-US034 commencing on 1 February 2013 have been awarded to the Participant by the Australian Renewable Energy Agency and NREL will be collaborating with the Participant during the final 28 months of this project. This project seeks to provide basic knowledge required to design solar thermochemical reactors able to perform the required energy conversions. In several proposed and demonstrated reactors, concentrated sunlight directly irradiates small solid particles suspended in fluid, enabling very high heat transfer rates to the particles which are the sites of chemical reaction. The reactors, therefore, involve the complex and couple dynamics of turbulent, chemically reacting, particle-laden flows and their interaction with concentrating solar radiation. A strong understanding of these coupled interactions will be crucial important in predicting and optimizing the performance of prototype reactors, but this understanding does not yet exist, since they have never been studied in any fundamental way. The project has a assembled an internationally leading team from The University of New South Wales (UNSW) and the University of Adelaide in Australia and the NREL in the United States to address this key gap in available know-how. The project will use U.S. Dept. of Energy (DOE) supercomputers, among the most powerful available worldwide, with cutting-edge software tools to perform first-principles simulations of the relevant interactions. These studies will be combined with detailed laser-based measurements in Australia to provide the first comprehensive databases concerning the governing phenomena in directly irradiated solar-thermochemical reactors. The outcomes will be the basic scientific knowledge, engineering knowhow and modeling tools necessary to design new reactor concepts and then scale up from the laboratory bench to practical size systems.

14 SOLAR ENERGY↗

Granular Temperature Modeling with Large-Scale CFD-DEM Data

The granular temperature is a measure of solids-phase fluctuating kinetic energy in particle-laden gas-solid flows. This report focuses on reduced-order models of granular temperature in triply-periodic cluster induced turbulent flow. The granular temperature data was generated by post-processing large-scale CFD-DEM simulations that were carried out as part of an ALCC awarded study. The algebraic model of Tang et al, originally fit to DNS data, was found to be amenable to the present data. The results of several regression analyses are presented building on the original expression. This work was performed in part through the National Nuclear Security Administration (NNSA) Minority Serving Institutions Internship Program (MSIIP).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a multiphase PIC model for slurry flow modeling

MFIX-Exa is a recently released multiphase CFD code originally developed for the simulation of particle-laden gas-solid flows. Due to its high-performance computing capabilities, MFIX-Exa is an ideal candidate for scale-up studies of slurry reactors, specifically the coarse-grained particle-in-cell (PIC) model with its statistical treatment of the particle phase. Unfortunately, several physical models that were neglected during original development because they are not relevant for high-density ratio gas-solid flows are important in slurry flows where the particle-to-fluid density ratio is near unity. In this preliminary work we focus on the effective (suspension) viscosity. The models of Brinkman (1952), Krieger and Dougherty (1956), and Cheng and Law (2003) are considered. The impact of the effective viscosity model is studied on horizontal pipe flow. The experimental data of Gillies et al. (2002) is used to assess the pressure drop predictions.

Fullmer, William D.↗

Using AMReX's Embedded Boundaries to Support MFIX-Exa’s Geometry Capabilities

MFIX-Exa is an AMReX-based application used for the simulation of chemically reacting multiphase flows inside reactors. The reactor geometry is mostly a combination of primitive shapes such as cylinders, cubes and spheres of varying dimensions connected together with potentially a lot of dead space between them. Additionally, the reactors may have particle-laden jets introduced at different parts of their geometry to enhance reactant mixing. This talk presents how MFIX-Exa uses OpenSCAD’s csg format to allow users to define the embedded boundary (EB) geometry and how it supports specifying inflow through certain sections of the EB. The talk will describe how MFIX-Exa leverages the EB features of grid-pruning, restart files and redistribution. Finally, the talk goes over some capabilities which are currently not available in MFIX-Exa but planned for the future, such as moving parts of the reactor geometry.

embedded boundaries↗

ExINP NSA Ice Nucleating Particle Concentrations

This data set comprises cumulative ambient ice-nucleating particle (INP) concentrations measured at the National Oceanic and Atmospheric Administration's (NOAA’s) Barrow Atmospheric Baseline Observatory (71.3230° N, 156.6114° W, “BRW” hereafter), next to the Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) site and ~ 6 km northeast of the town of Utqiaġvik. Our INP abundance data were generated using a combination of online instrument, the Portable Ice Nucleation Experiment chamber ver. 3 (PINE-03), and an offline cold stage, the West Texas Cryogenic Refrigerator Applied to Freezing Test system (WT-CRAFT). Our online INP data are all from the Examining the Ice-Nucleating Particles from NSA (ExINP-NSA) campaign conducted from October 19, 2021 to May 24, 2024. The offline INP concentration analysis was performed at West Texas A&M University for aerosol particle samples collected on polycarbonate filters (with 0.2-micron diameter pores). The PINE-03 measurements, as well as sampling activities for offline INP measurements, were conducted using the BRW site. An inset laminar sampling stack was mounted to the instrument platform, allowing PINE-03 to intake particle-laden air. For most of the campaign period, the semi-autonomous PINE-03 chamber was remotely controlled from West Texas A&M University using the LabView interface through the BeyondTrust remote-access console. PINE-03 was set to conduct an immersion freezing expansion experiment (i.e., simulated adiabatic cooling along with RHw at or above 100%). Except during the scheduled maintenance periods, the time resolution of each expansion experiment was approximately 12 minutes. PINE-03 continuously measured INP concentrations during the entire campaign without any substantial breaks. For most of the campaign period, PINE scanned its set-point vessel air temperatures from -14 °C to -31 °C and back to -14 °C about every 120 minutes.

54 ENVIRONMENTAL SCIENCES↗

ExINP NSA Ice Nucleating Particle Concentrations

This data set comprises cumulative ambient ice-nucleating particle (INP) concentrations measured at the National Oceanic and Atmospheric Administration's (NOAA’s) Barrow Atmospheric Baseline Observatory (71.3230° N, 156.6114° W, “BRW” hereafter), next to the Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) site and ~ 6 km northeast of the town of Utqiaġvik. Our INP abundance data were generated using a combination of online instrument, the Portable Ice Nucleation Experiment chamber ver. 3 (PINE-03), and an offline cold stage, the West Texas Cryogenic Refrigerator Applied to Freezing Test system (WT-CRAFT). Our online INP data are all from the Examining the Ice-Nucleating Particles from NSA (ExINP-NSA) campaign conducted from October 19, 2021 to May 24, 2024. The offline INP concentration analysis was performed at West Texas A&M University for aerosol particle samples collected on polycarbonate filters (with 0.2-micron diameter pores). The PINE-03 measurements, as well as sampling activities for offline INP measurements, were conducted using the BRW site. An inset laminar sampling stack was mounted to the instrument platform, allowing PINE-03 to intake particle-laden air. For most of the campaign period, the semi-autonomous PINE-03 chamber was remotely controlled from West Texas A&M University using the LabView interface through the BeyondTrust remote-access console. PINE-03 was set to conduct an immersion freezing expansion experiment (i.e., simulated adiabatic cooling along with RHw at or above 100%). Except during the scheduled maintenance periods, the time resolution of each expansion experiment was approximately 12 minutes. PINE-03 continuously measured INP concentrations during the entire campaign without any substantial breaks. For most of the campaign period, PINE scanned its set-point vessel air temperatures from -14 °C to -31 °C and back to -14 °C about every 120 minutes.

54 ENVIRONMENTAL SCIENCES↗

Development of a Deep Learning Model for Predicting the Drag Coefficients of Spherical and Non-Spherical Particles,

There is yet to be a well-established drag model for non-spherical particles required in a particle-laden flow that could cover a wide range of sphericities. This talk will explore the development of a general drag model for non-spherical particles by applying deep learning using available experimental data available in the literature. The integration of several raw experimental measurements from different sources and research directions allows the training of robust Artificial Intelligence and Machine Learning (ML) models. Neural networks are an ML approach inspired by the inner biological workings of the brain. This work aims to develop a Deep Neural Network (DNN) that predicts drag coefficient values with the ability to adapt appropriately to unseen data. Given the limited number of data points available and the variance found within the data collected from various sources, challenges may arise when looking to train the model. Our study tests and implements various model regularization techniques and assesses different loss and activation functions for the proposed DNN. The proposed model considers a broader range of features other than sphericity and Reynold number. These features include density ratio, solid volume fraction, lengthwise and crosswise sphericity, and more. Furthermore, we present the features that play a significant role in predicting different drag coefficients through feature importance. Within the investigated parameter ranges in this study, the following conclusions can be achieved and summarized below: • An improved drag coefficient model can be developed by considering more features such as, aspect ratio, lengthwise sphericity, crosswise sphericity, and density ratio. • DNN model can predict better results compared to traditional methods using MAE metric. • The proposed model addresses data challenges such as limited data and extreme data points through expanded feature-set and regularization. • Three major features that mostly affect the drag coefficient were identified from a feature importance analysis.

Presa-Reyes, Maria↗

Two-dimensional glass transition–like behavior of Janus particle–laden interface

Abstract Understanding the interactive behavior of Janus particles (JPs) is a growing field of research. The enhancement in binding energy, in comparison to homogenous particles, and the dual characteristic of JPs open up new possibilities for novel applications. In many such applications, interfacial materials become subjected to flows that produce dilational and shear stresses. Therefore, it is important to understand the impact that the Janus character brings to interfaces. In this work, we study the microstructure of two-dimensional (2D) JP monolayers formed at the air–water interface and examine the shear viscoelasticity with an interface rheometer that was adapted for in situ surface pressure control via a Langmuir trough. We extend concepts from bulk rheology to data obtained from interfacial rheology as a tool to understand and predict the monolayer’s viscoelastic behavior. Finally, by calculating the time relaxation spectrum from the measured 2D dynamic moduli, we conclude that a phenomenon similar to glass transition is taking place by analogy.

42 ENGINEERING↗

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. Furthermore, the model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗

Shock interacting with a random array of stationary particles underwater

Accurately predicting the force experienced by particles within a random bed during shock traversal is a challenging problem. In this paper, particle-resolved inviscid simulations of shock propagation over randomly distributed beds are performed, with the goal of quantifying the force on individual particles. Water was considered the continuum material, and then results were compared with previous research completed in air and the differences resulting from these different mediums were examined. Simulations were conducted for four different combinations of incident Mach number and bed volume fraction with the particles remaining stationary. Additionally, time-resolved streamwise and transverse force coefficients were calculated for each particle in the bed. It was observed that while the average force coefficient was similar to that of an isolated particle, the force on individual particles was substantially different. An important observation was that the inviscid drag force on some particles remained consistently positive, but others consistently negative even long after the passage of the shock. This persistent force on the particle was attributed to the potential flow fields that result from neighbors' specific locations upstream and downstream. Similar behavior was observed in the transverse force, with particles experiencing sustained force in the same lateral direction. These persistent forces are inviscid quasisteady contributions which with the viscous counterpart will play an essential role in the particle's long-term dispersion if allowed to move in response to the force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗