skip to main content

Explore Research Products in the NSF-PAR

Title: Clustering in Euler–Euler and Euler–Lagrange simulations of unbounded homogeneous particle-laden shear
Particle-laden flows of sedimenting solid particles or droplets in a carrier gas have strong inter-phase coupling. Even at low particle volume fractions, the two-way coupling can be significant due to the large particle to gas density ratio. In this semi-dilute regime, the slip velocity between phases leads to sustained clustering that strongly modulates the overall flow. The analysis of perturbations in homogeneous shear reveals the process by which clusters form: (i) the preferential concentration of inertial particles in the stretching regions of the flow leads to the formation of highly concentrated particle sheets, (ii) the thickness of the latter is controlled by particle-trajectory crossing, which causes a local dispersion of particles, (iii) a transverse Rayleigh–Taylor instability, aided by the shear-induced rotation of the particle sheets towards the gravity normal direction, breaks the planar structure into smaller clusters. Simulations in the Euler–Lagrange formalism are compared to Euler–Euler simulations with the two-fluid and anisotropic-Gaussian methods. It is found that the two-fluid method is unable to capture the particle dispersion due to particle-trajectory crossing and leads instead to the formation of discontinuities. These are removed with the anisotropic-Gaussian method which derives from a kinetic approach with particle-trajectory crossing in mind.  more » « less
Award ID(s):
1505795
NSF-PAR ID:
10095329
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
859
ISSN:
0022-1120
Page Range / eLocation ID:
174 to 203
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Journal of Fluid Mechanics)
    null (Ed.)
    In this work, model closures of the multiphase Reynolds-averaged Navier–Stokes (RANS) equations are developed for homogeneous, fully developed gas–particle flows. To date, the majority of RANS closures are based on extensions of single-phase turbulence models, which fail to capture complex two-phase flow dynamics across dilute and dense regimes, especially when two-way coupling between the phases is important. In the present study, particles settle under gravity in an unbounded viscous fluid. At sufficient mass loadings, interphase momentum exchange between the phases results in the spontaneous generation of particle clusters that sustain velocity fluctuations in the fluid. Data generated from Eulerian–Lagrangian simulations are used in a sparse regression method for model closure that ensures form invariance. Particular attention is paid to modelling the unclosed terms unique to the multiphase RANS equations (drag production, drag exchange, pressure strain and viscous dissipation). A minimal set of tensors is presented that serve as the basis for modelling. It is found that sparse regression identifies compact, algebraic models that are accurate across flow conditions and robust to sparse training data. 
    more » « less
  2. ; ; ; ; ; ; ( , Atmospheric Measurement Techniques)
    Flow tube reactors are widely employed to study gas-phase atmospheric chemistry and secondary organic aerosol (SOA) formation. The development of a new laminar-flow tube reactor, the Caltech Photooxidation Flow Tube (CPOT), intended for the study of gas-phase atmospheric chemistry and SOA formation, is reported here. The present work addresses the reactor design based on fluid dynamical characterization and the fundamental behavior of vapor molecules and particles in the reactor. The design of the inlet to the reactor, based on computational fluid dynamics (CFD) simulations, comprises a static mixer and a conical diffuser to facilitate development of a characteristic laminar flow profile. To assess the extent to which the actual performance adheres to the theoretical CFD model, residence time distribution (RTD) experiments are reported with vapor molecules (O3) and submicrometer ammonium sulfate particles. As confirmed by the CFD prediction, the presence of a slight deviation from strictly isothermal conditions leads to secondary flows in the reactor that produce deviations from the ideal parabolic laminar flow. The characterization experiments, in conjunction with theory, provide a basis for interpretation of atmospheric chemistry and SOA studies to follow. A 1-D photochemical model within an axially dispersed plug flow reactor (AD-PFR) framework is formulated to evaluate the oxidation level in the reactor. The simulation indicates that the OH concentration is uniform along the reactor, and an OH exposure (OHexp) ranging from ∼ 109 to ∼ 1012 molecules cm−3 s can be achieved from photolysis of H2O2. A method to calculate OHexp with a consideration for the axial dispersion in the present photochemical system is developed. 
    more » « less
  3. ; ; ( , Journal of Fluid Mechanics)
    Using a numerical model, we analyse the effects of shape on both the orientation and transport of anisotropic particles in wavy flows. The particles are idealized as prolate and oblate spheroids, and we consider the regime of small Stokes and particle Reynolds numbers. We find that the particles preferentially align into the shear plane with a mean orientation that is solely a function of their aspect ratio. This alignment, however, differs from the Jeffery orbits that occur in the residual shear flow (that is, the Stokes drift velocity field) in the absence of waves. Since the drag on an anisotropic particle depends on its alignment with the flow, this preferred orientation determines the effective drag on the particles, which in turn impacts their net downstream transport. We also find that the rate of alignment of the particles is not constant and depends strongly on their initial orientation; thus, variations in initial particle orientation result in dispersion of anisotropic-particle plumes. We show that this dispersion is a function of the particle’s eccentricity and the ratio of the settling and wave time scales. Due to this preferential alignment, we find that a plume of anisotropic particles in waves is on average transported farther but dispersed less than it would be if the particles were randomly oriented. Our results demonstrate that accurate prediction of the transport of anisotropic particles in wavy environments, such as microplastic particles in the ocean, requires the consideration of these preferential alignment effects. 
    more » « less
  4. ; ; ( , Journal of Fluid Mechanics)
    A computational study was performed both of a single agglomerate and of the collision of two agglomerates in a shear flow. The agglomerates were extracted from a direct numerical simulation of a turbulent agglomeration process, and had the loosely packed fractal structure typical of agglomerate structures formed in turbulent agglomeration processes. The computation was performed using a discrete-element method for adhesive particles with four-way coupling, accounting both for forces between the fluid and the particles (and vice versa ) as well as force transmission directly between particles via particle collisions. In addition to understanding and characterizing the particle dynamics, the study focused on illuminating the fluid flow field induced by the agglomerate in the presence of a background shear and the effect of collisions on this particle-induced flow. Perhaps the most interesting result of the current work was the observation that the flow field induced by a particle agglomerate rotating in a shear flow has the form of two tilted vortex rings with opposite-sign circulation. These rings are surrounded by a sea of stretched vorticity from the background shear flow. The agglomerate rotates in the shear flow, but at a slower rate than the ambient fluid elements. In the computations with two colliding agglomerates, we observed cases resulting in agglomerate merger, bouncing and fragmentation. However, the bouncing cases were all observed to also result in an exchange of particles between the two colliding agglomerates, so that they were influenced both by elastic rebound of the agglomerate structures as well as by tearing away of particulate matter between the agglomerates. Overall, the problems of agglomerate–flow interaction and of the collision of two agglomerates in a shear flow are considerably richer in physical phenomena and more complex than can be described by the common approximation that represents each agglomerate by an ‘equivalent sphere’. 
    more » « less
  5. ; ; ; ( , Journal of Fluid Mechanics)
    We examine the linear stability of a homogeneous gas–solid suspension of small Stokes number particles, with a moderate mass loading, subject to a simple shear flow. The modulation of the gravitational force exerted on the suspension, due to preferential concentration of particles in regions of low vorticity, in response to an imposed velocity perturbation, can lead to an algebraic instability. Since the fastest growing modes have wavelengths small compared with the characteristic length scale ( $U_{g}/{\it\Gamma}$ ) and oscillate with frequencies large compared with ${\it\Gamma}$ , $U_{g}$ being the settling velocity and ${\it\Gamma}$ the shear rate, we apply the WKB method, a multiple scale technique. This analysis reveals the existence of a number density mode which travels due to the settling of the particles and a momentum mode which travels due to the cross-streamline momentum transport caused by settling. These modes are coupled at a turning point which occurs when the wavevector is nearly horizontal and the most amplified perturbations are those in which a momentum wave upstream of the turning point creates a downstream number density wave. The particle number density perturbations reach a finite, but large amplitude that persists after the wave becomes aligned with the velocity gradient. The growth of the amplitude of particle concentration and fluid velocity disturbances is characterised as a function of the wavenumber and Reynolds number ( $\mathit{Re}=U_{g}^{2}/{\it\Gamma}{\it\nu}$ ) using both asymptotic theory and a numerical solution of the linearised equations. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email