Direct numerical simulations of sedimenting spherical particles at non-zero Reynolds number

RSC Advances ◽  
2014 ◽  
Vol 4 (96) ◽  
pp. 53681-53693 ◽  
Author(s):  
Adnan Hamid ◽  
John J. Molina ◽  
Ryoichi Yamamoto
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Dynamic Properties ◽  
Direct Numerical Simulations ◽  
Spherical Particles ◽  
Inertial Effects ◽  
Profile Method ◽  
Volume Fractions ◽  
Wide Range ◽  
Smoothed Profile Method

We performed direct numerical simulations, using a smoothed profile method to investigate the inertial effects on the static and dynamic properties of a sedimenting suspension over a wide range of volume fractions from 0.01 to 0.4.

scholarly journals Smoothed profile method for direct numerical simulations of hydrodynamically interacting particles

Soft Matter ◽  
2021 ◽  
Author(s):  
Ryoichi Yamamoto ◽  
John Jairo Molina ◽  
Yasuya Nakayama
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Standard Procedure ◽  
Direct Numerical Simulations ◽  
Interacting Particles ◽  
Profile Method ◽  
Smoothed Profile Method ◽  
General Method

A general method is presented for computing the motions of hydrodynamically interacting particles in various kinds of host fluids for arbitrary Reynolds number. The method follows the standard procedure for...

Homogeneity of turbulence generated by static-grid structures

2010 ◽  
Vol 654 ◽  
pp. 473-500 ◽  
Author(s):  
Ö. ERTUNÇ ◽  
N. ÖZYILMAZ ◽  
H. LIENHART ◽  
F. DURST ◽  
K. BERONOV
Keyword(s):  
Reynolds Number ◽  
Velocity Field ◽  
Numerical Simulations ◽  
Direct Numerical Simulations ◽  
Uniform Grid ◽  
Reynolds Stresses ◽  
Hot Wire ◽  
Mean Velocity ◽  
Wide Range ◽  
The Mean

Homogeneity of turbulence generated by static grids is investigated with the help of hot-wire measurements in a wind-tunnel and direct numerical simulations based on the Lattice Bolztmann method. It is shown experimentally that Reynolds stresses and their anisotropy do not become homogeneous downstream of the grid, independent of the mesh Reynolds number for a grid porosity of 64%, which is higher than the lowest porosities suggested in the literature to realize homogeneous turbulence downstream of the grid. In order to validate the experimental observations and elucidate possible reasons for the inhomogeneity, direct numerical simulations have been performed over a wide range of grid porosity at a constant mesh Reynolds number. It is found from the simulations that the turbulence wake behind the symmetric grids is only homogeneous in its mean velocity but is inhomogeneous when turbulence quantities are considered, whereas the mean velocity field becomes inhomogeneous in the wake of a slightly non-uniform grid. The simulations are further analysed by evaluating the terms in the transport equation of the kinetic energy of turbulence to provide an explanation for the persistence of the inhomogeneity of Reynolds stresses far downstream of the grid. It is shown that the early homogenization of the mean velocity field hinders the homogenization of the turbulence field.

scholarly journals Waves and vortices in the inverse cascade regime of stratified turbulence with or without rotation

2016 ◽  
Vol 806 ◽  
pp. 165-204 ◽  
Author(s):  
Corentin Herbert ◽  
Raffaele Marino ◽  
Duane Rosenberg ◽  
Annick Pouquet
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Direct Numerical Simulations ◽  
Vertical Shear ◽  
Energy Spectra ◽  
Stratified Turbulence ◽  
Viscous Effects ◽  
Inverse Cascade ◽  
Main Effects ◽  
The Waves

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behaviour of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number $Re\approx 1000$, with the ratio between the Brunt–Väisälä frequency $N$ and the inertial frequency $f$ varying from $1/4$ to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range $0.02\leqslant Fr\leqslant 0.32$. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio $N/f$ increases, also in the stratified case, although less and less so. However, when $N/f$ increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing $N/f$ can be interpreted in terms of the waves and vortices, and identify the main effects that can explain this transition based on both inviscid invariants arguments and viscous effects due to vertical shear.

scholarly journals Turbulent Coagulation of Aerosol Particles: New Insights From Direct Numerical Simulations and Holographic Imaging Experiments

2003 ◽  
Author(s):  
Lance R. Collins ◽  
Hui Meng ◽  
Aruj Ahluwalia ◽  
Lujie Cao ◽  
Gang Pan
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Search Algorithm ◽  
Concentration Field ◽  
Direct Numerical Simulations ◽  
Three Dimensions ◽  
Collision Kernel ◽  
Coagulation Rate ◽  
Subsequent Work ◽  
Turbulent Fluctuations

Particle collisions driven by turbulent fluctuations play a key role in such diverse problems as cloud formation, aerosol powder manufacturing and inhalation drug therapy to name a few. In all of these examples (and many others) turbulent fluctuations increase the rate of collisions relative to the background collision rate driven by Brownian motion. Furthermore, turbulence can spontaneously generate very large fluctuations in the particle concentration field. This “clustering” is caused by the inertial mismatch between the heavy particles and the lighter surrounding gas; vortices in the flow “centrifuge” the heavier particles out of vortex cores and into the straining regions that lie in between the vortices. Because collision is a binary process, concentration fluctuations further enhance the turbulent coagulation rate by as much as two orders of magnitude. An effect of this size must be accounted for in a rational model of turbulent coagulation. Sundaram & Collins (J. Fluid Mech. 1997) showed that the radial distribution function (RDF) of the particle population, evaluated at contact, precisely corrects the collision kernel for clustering. Subsequent work has explored the dependence of the RDF on the system parameters (e.g., particle size, concentration, response time and Reynolds number) using direct numerical simulations. These results have improved our understanding and ability to predict the effect of the first three parameters; however, owing to the limited range of Reynolds number that can be reached in a numerical simulation, questions remain over the scaling of the RDF with Reynolds number. This is a critical issue for high-Reynolds-number applications such as cloud physics, where values of the Reynolds number can be 1–2 orders of magnitude greater than can be simulated. We will present our highest Reynolds number simulations to date and show our attempts to resolve this issue. Recently, the ability to measure three-dimensional particle positions using holography has been realized (e.g., Meng & Pu, J. Opt. Soc. Am. 2003). With holography, the optical image that is produced contains fringes that, upon inverting the laser, reproduce the original image in three dimensions. The hologram can then be scanned using a digital camera to obtain the particle positions. An important consideration with this study is the need to differentiate individual particles. We developed a search algorithm that locates particle centers, even in the presence of optical aberations and speckle noise. The algorithm has been used to obtain the first experimental RDF measurements to date. Thus far we see good agreement between the experimentally obtained RDF and the simulations. Besides validating the simulations, experiments can span a much broader range of Reynolds numbers, providing critical data that may help resolve the open questions associated with this parameter.

scholarly journals Comparison between experiments and direct numerical simulations in a channel flow with roughness on one wall

2008 ◽  
Vol 600 ◽  
pp. 403-426 ◽  
Author(s):  
P. BURATTINI ◽  
S. LEONARDI ◽  
P. ORLANDI ◽  
R. A. ANTONIA
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Friction Velocity ◽  
Streamwise Velocity ◽  
Negative Energy ◽  
Direct Numerical Simulations ◽  
Rough Wall ◽  
Local Isotropy ◽  
Taylor’S Hypothesis ◽  
Yielding Support

The turbulent flow in a two-dimensional channel with roughness on one wall is investigated using experiments and direct numerical simulations (DNS). The elements have a square cross-section with height k=0.1H (H is the channel half-width) and a streamwise spacing of 4k. The Reynolds number Reτr, based on the friction velocity at the rough wall and H, is in the range 300–1100. Particular attention is given to the rough-wall side. Measured turbulence intensities, length scales, leading terms in the turbulent kinetic energy budget, and velocity spectra are compared with those obtained from the DNS. Close agreement is found, yielding support for the simplifying assumptions in the experiment (notably local isotropy and Taylor's hypothesis) and the adequacy of the spatial resolution in the simulation. Overall, the profiles of the Reynolds normal stresses on the roughness side are almost independent of Reτr, when normalized by outer variables. Energy spectra at different locations above the rough wall collapse well at high wavenumbers, when normalized by Kolmogorov scales. In contrast to previous studies, a region of negative energy production near the location of the maximum streamwise velocity is not observed. Comparison with a smooth-wall channel, at similar values of the friction-velocity Reynolds number, highlights differences only in the streamwise velocity component near the wall.

Lagrangian characteristics of turbulence and scalar transport in direct numerical simulations

2001 ◽  
Vol 427 ◽  
pp. 241-274 ◽  
Author(s):  
P. K. YEUNG
Keyword(s):  
Reynolds Number ◽  
Correlation Function ◽  
Numerical Simulations ◽  
Time Scales ◽  
Time Scale ◽  
Cross Correlation ◽  
Scalar Fields ◽  
Direct Numerical Simulations ◽  
Cross Correlation Function ◽  
Differential Diffusion

A study of the Lagrangian statistical properties of velocity and passive scalar fields using direct numerical simulations is presented, for the case of stationary isotropic turbulence with uniform mean scalar gradients. Data at higher grid resolutions (up to 5123 and Taylor-scale Reynolds number 234) allow an update of previous velocity results at lower Reynolds number, including intermittency and dimensionality effects on vorticity time scales. The emphasis is on Lagrangian scalar time series which are new to the literature and important for stochastic mixing models. The variance of the ‘total’ Lagrangian scalar value (ϕ˜+, combining contributions from both mean and fluctuations) grows with time, with the velocity–scalar cross-correlation function and fluid particle displacements playing major roles. The Lagrangian increment of ϕ˜+ conditioned upon velocity and scalar fluctuations is well represented by a linear regression model whose parameters depend on both Reynolds number and Schmidt number. The Lagrangian scalar fluctuation is non-Markovian and has a longer time scale than the velocity, which is due to the strong role of advective transport, and is in contrast to results in an Eulerian frame where the scalars have shorter time scales. The scalar dissipation is highly intermittent and becomes de-correlated in time more rapidly than the energy dissipation. Differential diffusion for scalars with Schmidt numbers between 1/8 and 1 is characterized by asymmetry in the two-scalar cross-correlation function, a shorter time scale for the difference between two scalars, as well as a systematic decrease in the Lagrangian coherency spectrum up to at least the Kolmogorov frequency. These observations are consistent with recent work suggesting that differential diffusion remains important in the small scales at high Reynolds number.

Effect of Size and Slip Coefficient of a Porous Manifold on Thermal Stratification in Storage Tanks

2009 ◽  
Author(s):  
N. M. Brown ◽  
F. C. Lai
Keyword(s):  
Experimental Data ◽  
Reynolds Number ◽  
Numerical Simulations ◽  
Storage Tank ◽  
Richardson Number ◽  
Thermal Stratification ◽  
Temperature Fields ◽  
Storage Tanks ◽  
Slip Coefficient ◽  
Wide Range

Numerical simulations have been performed to study the effects of size and slip coefficient of a porous manifold on the thermal stratification in a storage tank. The model is used to predict the development of flow and temperature fields during a charging process. Computations have covered a wide range of the Grashof number (1.8 × 105 < Gr < 1.8 × 108) and Reynolds number (10 ≤ Re ≤ 104), or in terms of the Richardson number, 10−2 < Ri < 105. The results obtained compare favorably well with the experimental data. In addition, the present results have confirmed the effectiveness of porous manifold in the promotion of thermal stratification and provide useful information for the design of such system.

scholarly journals Direct numerical simulations of dense suspensions: wave instabilities in liquid-fluidized beds

2007 ◽  
Vol 587 ◽  
pp. 303-336 ◽  
Author(s):  
J. J. DERKSEN ◽  
S. SUNDARESAN
Keyword(s):  
Numerical Simulations ◽  
Bulk Viscosity ◽  
Fluid Model ◽  
Travelling Waves ◽  
System Size ◽  
Direct Numerical Simulations ◽  
Fluid Particle ◽  
Spherical Particles ◽  
Particle Phase ◽  
Viscosity Term

We present results of direct numerical simulations of travelling waves in dense assemblies of monodisperse spherical particles fluidized by a liquid. The cases we study have been derived from the experimental work of others. In these simulations, the flow of interstitial fluid is solved by the lattice-Boltzmann method (LBM) and the particles move under the influence of gravity, hydrodynamic forces stemming from the LBM, subgrid-scale lubrication forces and hard-sphere collisions. We first show that the propagating inhomogeneous structures seen in the simulations are in agreement with those observed experimentally. We then use the detailed information contained in the simulation results to assess aspects of two-fluid model closures, namely, fluid–particle drag, and the various contributions to the effective stresses. We show that the rates of compaction and dilation of the particle phase in the travelling waves are comparable to the rate at which the microstructure relaxes, and that there is a pronounced effect of the rate of compaction on the average collisional normal stress. Although this effect can be expressed as an effective bulk viscosity term, this approach would require the use of a path-dependent bulk viscosity. We also find that the effective fluid–particle drag coefficient can be described well with the often-used closure motivated by the experiments of Richardson & Zaki (Trans. Inst. Chem. Engng vol. 32, 1954, p. 35). In this respect, the effect of the system size for determining the drag requires specific care. The shear viscosity of the particle phase manifests small, but clearly noticeable dependence on the rate of compaction/dilation of the particle phase. Our observations point to the need for higher-order closures that recognize the slow evolution of the microstructure in these flows and account for the effects of non-equilibrium microstructure on the stresses.

Sign in / Sign up

Export Citation Format

Share Document