Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle-laden flows”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Deep learning for drag force modelling in dilute, poly-dispersed particle-laden flows with irregular-shaped particles

Here, this study applies machine learning-based approaches to develop a drag force model for irregular-shaped particles in incompressible flows. The particle-laden flows are studied through an in-house particle-resolved direct numerical simulation (PR-DNS) at low-intermediate Reynolds numbers (Re). We utilize the PR-DNS to obtain drag force coefficients and flow fields of single particles. A variational auto-encoder model is used to obtain latent vectors to represent the geometrical features of the particles, and artificial neural networks (ANN) are developed to predict drag force coefficients and flow fields of the single particles. This study applies a pairwise interaction extended point-particle (PIEP) model to obtain the coefficients assuming the flow fields of neighboring particles can be linearly superposed. The PIEP method shows significant improvement on prediction for a few neighbored particles. In addition, the results reveal R2 scores of 0.56-0.62 and errors of 9.1-10.0 % for the dilute, polydispersed systems with a volume fraction of 0.5 %.

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

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Debris-induced consequences on turbulence and vorticity in solar photovoltaic module-generated array wakes

Particle-laden flows in solar photovoltaic (PV) systems are inevitable, where wind-swept debris in open environments are carried by high winds and turbulence, coating panel surfaces or damaging structures. Particle deposition, or soiling, is a well-known issue for large-scale plants which rely on uninhibited solar rays for optimal production. But understanding the mechanisms leading to soiling requires a physical and fluid dynamics-centered focus, since turbulence dominates PV panel wakes and is also known to alter particle concentration and trajectories. This study presents an experimental campaign toward consequences of particle-laden flow between two model PV panels using time-resolved particle image velocimetry. The model array was subjected to varied particle volume fractions, including a tracer particle case and a water droplet case. Characterization of mean velocity, turbulence statistics, and mean kinetic energy within the single phase and, separately, particle phase flows showed modified features due to particle inertia. Images captured at a frequency of 1 kHz in the near wake of the upstream panel allow for a first experimental look at vorticity and convective velocity of vortex structures for single-phase and particle-phase flows which are crucial to debris transport and soiling in PV environments.

Energy & Fuels↗
DOE OSTI2024

Interactions of laser-driven tin ejecta microjets over phase transition boundaries

Ejecta microjets offer an experimental methodology to study high-speed particle laden-flow interactions, as microjets consist of millions of particulates traveling at velocities of several kilometers per second and are easily generated by most common shock drives. Previous experiments on the OMEGA Extended Performance laser found that collisions between two counter-propagating laser-driven tin ejecta microjets varied as a function of drive pressure; jets generated near shock pressures of 10 GPa passed through each other without interacting, whereas jets generated at shock pressures of over 100 GPa interacted strongly, forming a cloud around the center interaction point. In this paper, we present a more systematic scan of tin ejecta microjet collisions over intermediate pressure regimes to identify how and at what shock pressure interaction behavior onsets. Radiographs of interacting microjets at five different laser drive energies qualitatively demonstrate that interaction behavior onsets slowly as a function of laser drive energy. Quantitative mass and density metrics from each radiograph provide trends on jet characteristics and collisional mass dispersion. It is observed that jetting mass, jet densities, and mass dispersion increase with increasing drive pressures and that the increased jet density at the higher drive energies may account for the increased mass dispersion. This work provides an important step in the understanding of tin ejecta microjet collisions and paves the way for future studies on the physics dominating high-speed particle-laden flow interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

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↗
DOE OSTI2023

Particle transport-driven flow dynamics and heat transfer modulation in solar photovoltaic modules: Implications on soiling

Wind-swept debris and excessive heating reduce solar photovoltaic efficiency and accelerate panel degradation. For many solar plants in open environments, particle deposition is well-understood to exacerbate efficiency losses by blocking useable light and insulating panel surfaces. However, little is known about particle behavior in multi-panel systems, where complex wakes dominate and alter trajectories. This study examines debris-laden flow within a scaled array of two model solar panels, observing convective cooling and focusing on turbulent inter-panel flow dynamics. Convection is compared between the first and second panels (P 1 and P 2 ) and for a single inflow velocity at four particle concentrations (Φ v ), including a neutrally buoyant particle (i.e. tracer) case and three water droplet cases. Tracers and Φ v > 0 concentrations show increased convective heat transfer for P 2 compared to P 1 , while inertial particle presence and increased concentration reduces this improvement from +170% to +90%. Particle image velocimetry surrounding P 2 reveals dampened inter-panel wake features for particle-laden flow compared to the tracer case, where diminished mean flow and mean kinetic energy breakdown are observed for the particles due to their inertia. Together, the reduction in comparative convection and dampened wake effects accompanying inertial particle concentration motivates further studies on coupled array physics in debris-filled flow.

14 SOLAR ENERGY↗
DOE OSTI2024

Multi-frame x-ray radiography and image tracking for quantification of expansion in laser-driven tin ejecta microjets

One regime of experimental particle-laden flow study involves ejecta microjets—often defined as a stream of micrometer-scale particles generated through shock interaction with a non-uniform surface and generally travel above 1 km/s. In order to capture the change in characteristics as a function of propagation time, we apply a multi-frame x-ray radiography platform to observe and track the jet transport dynamics. A synchrotron x-ray source allows us to perform quantitative analyses and comparisons between the eight images captured by the imaging system. Observation of a single jet through time allows the use of a cross correlation algorithm to independently track various regions within the jet and quantify the jet expansion over time using normalized area and normalized areal density values. Through a comparison with the calculated values of ballistic transport, these findings show less expansion than expected for ballistically transporting particles. In conclusion, this work combines multi-frame synchrotron radiography with image tracking to establish a foundation for future studies on jet transport and particle interaction dynamics.

Sun, Yuchen [Lawrence Livermore National Laborator↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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.↗
DOE OSTI2024

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↗
DOE OSTI2024

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↗
DOE OSTI2026

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↗
DOE OSTI2025

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↗