skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Improved Discrete Random Walk Stochastic Model for Simulating Particle Dispersion and Deposition in Inhomogeneous Turbulent Flows
Abstract The performance of different versions of the discrete random walk models in turbulent flows with nonuniform normal root-mean-square (RMS) velocity fluctuations and turbulence time scales were carefully investigated. The OpenFOAM v2−f low Reynolds number turbulence model was used for evaluating the fully developed streamwise velocity and the wall-normal RMS velocity fluctuations profiles in a turbulent channel flow. The results were then used in an in-house matlab particle tracking code, including the drag and Brownian forces, and the trajectories of randomly injected point-particles with diameters ranging from 10 nm to 30 μm were evaluated under the one-way coupling assumption. The distributions and deposition velocities of fluid-tracer and finite-size particles were evaluated using the conventional-discrete random walk (DRW) model, the modified-DRW model including the velocity gradient drift correction, and the new improved-DRW model including the velocity and time gradient drift terms. It was shown that the conventional-DRW model leads to superfluous migration of fluid-point particles toward the wall and erroneous particle deposition rate. The concentration profiles of tracer particles obtained by using the modified-DRW model still are not uniform. However, it was shown that the new improved-DRW model with the velocity and time scale drift corrections leads to uniform distributions for fluid-point particles and reasonable concentration profiles for finite-size heavy particles. In addition, good agreement was found between the estimated deposition velocities of different size particles by the new improved-DRW model with the available data.  more » « less
Award ID(s):
1925596
PAR ID:
10287206
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluids Engineering
Volume:
142
Issue:
10
ISSN:
0098-2202
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The lattice Boltzmann method is employed to conduct direct numerical simulations of turbulent open channel flows with the presence of finite-size spherical sediment particles. The uniform particles have a diameter of approximately 18 wall units and a density of ρp=2.65ρf, where ρp and ρf are the particle and fluid densities, respectively. Three low particle volume fractions ϕ=0.11%, 0.22%, and 0.44% are used to investigate the particle-turbulence interactions. Simulation results indicate that particles are found to result in a more isotropic distribution of fluid turbulent kinetic energy (TKE) among different velocity components, and a more homogeneous distribution of the fluid TKE in the wall-normal direction. Particles tend to accumulate in the near-wall region due to the settling effect and they preferentially reside in low-speed streaks. The vertical particle volume fraction profiles are self-similar when normalized by the total particle volume fractions. Moreover, several typical transport modes of the sediment particles, such as resuspension, saltation, and rolling, are captured by tracking the trajectories of particles. Finally, the vertical profiles of particle concentration are shown to be consistent with a kinetic model. 
    more » « less
  2. A new Anechoic Wall Jet Wind Tunnel was built at Virginia Tech. A detailed design based on the old wall jet tunnel was done to improve the quality of the resultant flow. Aerodynamic and acoustic calibrations were performed in order to understand properties and characteristics of the flow generated by this new facility which can be used for various aeroacoustic studies. Far-field acoustics were measured using half-inch B&K microphones in a streamwise array to characterize and reduce the background noise. Sound pressure levels were lower by 10 dB for frequencies up to 700Hz in comparison to the old facility. The turbulent surface pressure fluctuations of the wall-jet flow were studied using Sennheiser microphones placed along streamwise and spanwise locations to record surface pressure fluctuations. Comparison of the autocorrelation plotted for microphones along the same span indicate uniform flow features. A decay in the turbulence levels is observed along the downstream direction as expected. Aerodynamic calibrations included mean velocity measurements along different spanwise locations, wall-jet boundary layer profiles and streamwise cross-sections. Spanwise and cross-sectional velocity profiles show good uniformity of the flow. Detailed boundary layer analyses were performed with the parameters obtained from the experiments. 
    more » « less
  3. We use a novel experimental setup to obtain the vertical velocity and acceleration statistics of snowflakes settling in atmospheric surface-layer turbulence, for Taylor microscale Reynolds numbers (Reλ) between 400 and 67 000, Stokes numbers (St) between 0.12 and 3.50, and a broad range of snowflake habits. Despite the complexity of snowflake structures and the non-uniform nature of the turbulence, we find that mean snowflake acceleration distributions can be uniquely determined from the value of St. Ensemble-averaged snowflake root mean square (rms) accelerations scale nearly linearly with St. Normalized by the rms value, the acceleration distribution is nearly exponential, with a scaling factor for the (exponent) of −3/2 that is independent of Reλ and St; kurtosis scales with Reλ, albeit weakly compared to fluid tracers in turbulence; gravitational drift with sweeping is observed for St < 1. Surprisingly, the same exponential distribution describes a pseudo-acceleration calculated from fluctuations of snowflake terminal fall speed in still air. This equivalence suggests an underlying connection between how turbulence determines the trajectories of particles and the microphysics determining the evolution of their shapes and sizes. 
    more » « less
  4. null (Ed.)
    ABSTRACT Recently, Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of a uniform external force (such as radiation pressure or gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvén or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma β, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvénic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme overdensities (up to ≫109 times mean), and exhibit substantial substructure even in nearly incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions. 
    more » « less
  5. Abstract Because of their buoyancy, rigidity, and finite size, inertial particles do not obey the same dynamics as fluid parcels. The motion of small spherical particles in a fluid flow is described by the Maxey–Riley equations and depends nonlinearly on the velocity of the fluid in which the particles are immersed. Fluid velocities in the ocean often have a strong small-scale turbulent component which is difficult to observe or model, presenting a challenge to predicting the evolution of distributions of inertial particles in the ocean. To overcome this challenge, we assume that the turbulent velocity imposes a random force on particles and consider a stochastic analog of the Maxey–Riley equations. By performing a perturbation analysis of the stochastic Maxey–Riley equations, we obtain a simple and accurate partial differential equation for the spatial distribution of particles. The equation is of the advection–diffusion type and handles the uncertainty introduced by unresolved turbulent flow features. In several numerical test cases, distributions of particles obtained by solving the newly derived equation compare favorably with distributions obtained from Monte Carlo simulations of individual particle trajectories and with theoretical predictions. The advection–diffusion form of our newly derived equation is amenable to inclusion within many existing ocean circulation models. Significance StatementWe introduce a new model for describing spatial distributions of small rigid objects, such as plastic debris, in the ocean. The model takes into account the effects of finite particle size and particle buoyancy, which cause particle trajectories to differ from fluid parcel trajectories. Our model also represents small-scale turbulence stochastically. 
    more » « less