Comprehensive quasi-steady force correlations for compressible flow through random particle suspensions
Not Available
42 ENGINEERING↗
SEARCH · Search NASA
Results for “Particle Drag”
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.
Estimate of Lorentz force effects on charged particles
Thermonuclear (TN) fusion reaction between deutrons are D + D → T + P D + D → 3 He + n where P is the proton, D is the deuteron, T is the triton, n is the neutron, and 3 He is the helion. They are born with high energies (0.82 MeV for 3 He, 1.01 MeV for T, 3.02 MeV for P, and 2.45 MeV for n). The charged particles slow down due to the drag force (stopping power). While slowing down, the charged particles deposit their energy, which may affect the TN burn characteristics. When magnetic field is present, the velocity of the charged particle would be affected due to the Lorentz force, which can affect the stopping of the charged particle. As a consequence, TN burn characteristics may change, Our purpose is to asses the importance of the Lorentz force in stopping of the charged particle, based on a simple comparison between the Lorentz force and stopping power.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Shock and contact interaction with a simple cubic array of particles
Shock-particle interaction is a fundamental pillar of multiphase compressible flows that has been studied at length for many decades. However, little attention has been paid to the interaction of particles with a contact interface that follows a shock in shock tube experiments and applications relating to blast waves. Presently, the phenomenon is studied at the microscale via particle resolved simulations of shock contact systems interacting with a structured array of particles as well as isolated particles. Simulations are conducted at particle volume fractions of 0%, 5%, 10%, 20%, and 40% at three contact Mach numbers. Additionally, the diaphragm position is varied, which controls the timing of the shock arrival time in relation to the contact arrival time. The modification to the drag on these stationary particles by the contact is analyzed and compared to the compressible Maxey–Riley–Gatignol model, which is adequate for the single particle cases but does not account for fluid mediated particle–particle interactions.
Mechanics↗
Shifting and Splitting of Resonance Lines due to Dynamical Friction in Plasmas
We report a quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Neural-network-based drag force model for polydisperse assemblies of irregular-shaped particles
Not provided.
Engineering↗
Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion
The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.
36 MATERIALS SCIENCE↗
The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation
We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.
42 ENGINEERING↗
Direct numerical simulation of turbulent dispersion of evaporative aerosol clouds produced by an intense expiratory event
Airborne particles are a major route for transmission of COVID-19 and many other infectious diseases. When a person talks, sings, coughs, or sneezes, nasal and throat secretions are spewed into the air. After a short initial fragmentation stage, the expelled material is mostly composed of spherical particles of different sizes. While the dynamics of the largest droplets are dominated by gravitational effects, the smaller aerosol particles, mostly transported by means of hydrodynamic drag, form clouds that can remain afloat for long times. In subsaturated air environments, the dependence of pathogen-laden particle dispersion on their size is complicated due to evaporation of the aqueous fraction. Particle dynamics can significantly change when ambient conditions favor rapid evaporation rates that result in a transition from buoyancyto-drag dominated dispersion regimes. To investigate the effect of particle size and evaporation on pathogen-laden cloud evolution, a direct numerical simulation of a mild cough was coupled with an evaporative Lagrangian particle advection model. The results suggest that while the dispersion of cough particles in the tails of the size distribution are unlikely to be disrupted by evaporative effects, preferential aerosol diameters (30-40 μm) may exhibit significant increases in the residence time and horizontal range under typical ambient conditions. Using estimations of the viral concentration in the spewed fluid and the number of ejected particles in a typical respiratory event, we obtained a map of viral load per volume of air at the end of the cough and the number of virus copies per inhalation in the emitter vicinity.
60 APPLIED LIFE SCIENCES↗
Numerical simulations of liquid jetting with solid inclusions
The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Evolution of the age-included nearest pair distribution in disperse multiphase flows
The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗
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↗
Initial Atmospheric Transport of Particles
A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.
54 ENVIRONMENTAL SCIENCES↗
Slow manifold reduction for plasma science
The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
On the Scaling and Growth Limit of Fluvial Dunes
Abstract Bedforms emerge in a variety of shapes and sizes when a granular erodible surface is subject to a strong enough shear flow, as observed in topographic data from submarine canyons, rivers, deserts, and planetary bodies. The two salient features of bedforms are the ability to collectively transport particles, and to generate form drag thus increasing flow resistance. These two mechanisms are in competition and contribute to force bedforms of increasing size migrating more slowly. In a dedicated large‐scale open channel flow facility, hierarchies of fluvial bedforms were generated and measured in equilibrium conditions. The corresponding scale‐dependent migration velocity and mass flux contributions were quantified in the wave number and frequency domains. Experimental results are paired with a validated set of theoretical models to demonstrate that ripples or dunes reach an equilibrium state when drag partitioning ensures high enough frictional drag to sustain the bedload transport of sediment, and low enough form drag to enable the migration of the largest bedform size. This mechanism is inferred to constrain the growth of bedforms when sediment supply is not the limiting factor.
58 GEOSCIENCES↗
Energy spectrum of lost alpha particles in magnetic mirror confinement
In a magnetic mirror fusion reactor, capturing the energy of fusion-produced alpha particles is essential to sustaining the reaction. However, since alpha particles are born at energies much higher than the confining potential, a substantial fraction are lost due to pitch-angle scattering before they can transfer their energy to the plasma via drag. The energy of lost alpha particles can still be captured through direct conversion, but designing an effective mechanism requires a description of the energies and times at which they become deconfined. Here we present analytical solutions for the loss velocity, energy, and time distributions of alpha particles in a magnetic mirror. After obtaining the Fokker-Planck collision operator, we asymptotically solve for the eigenfunctions of the Legendre operator to reveal a closed-form solution. In conclusion, our framework applies to any high-energy species, for any applied potential and mirror ratio R > 1, making this work broadly applicable to mirror devices.
Confinement↗
Predictive models enhance feedstock quality of corn stover via air classification
Feedstock heterogeneity is a fundamental obstacle to cost-competitive biobased products. Agricultural products like corn stover have anatomical components that vary in their chemical composition, mechanical properties, structure, and response to chemical and biological treatments. A technique that can enrich streams in select anatomical fractions would allow a tailored deconstruction approach to increase overall process efficiency. Air classification can be leveraged for such refining, however, fundamental characterization and understanding of the particle properties that underly the physics of air classification are only modestly documented. Here, we determine fundamental particle properties including mass-to-area ratio, drag coefficient, and partition velocity that describe how anatomical tissues of corn stover behave during air classification. In this work, mass-to-area ratios of anatomical tissues vary by nearly two orders of magnitude from 2.3 mg/mm 2 for cob to 0.04 mg/mm 2 for leaf. Drag coefficients of longer, fibrous materials (i.e., rind, husk, and sheath) are shown to correlate with particle area (p-value < 0.001) whereas granular tissues (i.e., cob, pith, and leaf) correlate better with mass-to-area ratio (p-values < 0.001). When compared to experimental observations, a simulated two-stage air classification and size reduction scenario predicts the overall partitioning of anatomical tissues within 15% for pith, husk, rind, and cob tissues. The model predicts an air-classified fraction preferentially enriched in cob (purity = 20%), rind (purity = 74%), and pith (purity = 4.5%) with a mass yield of 47%. Empirical relations for these properties can be used to predict the partitioning of corn stover during air classification based on anatomical type and size.
09 BIOMASS FUELS↗
A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach
Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.
42 ENGINEERING↗