• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Top 20 Most Cited Articles

The 20 most cited articles over time based on CrossRef data.


Direct numerical simulation of turbulent channel flow up to Reτ = 590

Robert D. Moser, John Kim, and Nagi N. Mansour

Phys. Fluids 11, 943 (1999); http://dx.doi.org/10.1063/1.869966 (3 pages) | Cited 484 times

Full Text: | Download PDF

Show Abstract
Numerical simulations of fully developed turbulent channel flow at three Reynolds numbers up to Reτ = 590 are reported. It is noted that the higher Reynolds number simulations exhibit fewer low Reynolds number effects than previous simulations at Reτ = 180. A comprehensive set of statistics gathered from the simulations is available on the web at http://www.tam.uiuc.edu/Faculty/Moser/channel. © 1999 American Institute of Physics.
Show PACS
47.27.-i Turbulent flows
47.60.-i Flow phenomena in quasi-one-dimensional systems
47.11.-j Computational methods in fluid dynamics

On pressure and velocity boundary conditions for the lattice Boltzmann BGK model

Qisu Zou and Xiaoyi He

Phys. Fluids 9, 1591 (1997); http://dx.doi.org/10.1063/1.869307 (8 pages) | Cited 290 times

Full Text: | Download PDF

Show Abstract
Pressure (density) and velocity boundary conditions are studied for 2-D and 3-D lattice Boltzmann BGK models (LBGK) and a new method to specify these conditions is proposed. These conditions are constructed in consistency with the wall boundary condition, based on the idea of bounceback of the non-equilibrium distribution. When these conditions are used together with the incompressible LBGK model [J. Stat. Phys. 81, 35 (1995)] the simulation results recover the analytical solution of the plane Poiseuille flow driven by a pressure (density) difference. The half-way wall bounceback boundary condition is also used with the pressure (density) inlet/outlet conditions proposed in this paper and in Phys. Fluids 8, 2527 (1996) to study 2-D Poiseuille flow and 3-D square duct flow. The numerical results are approximately second-order accurate. The magnitude of the error of the half-way wall bounceback boundary condition is comparable with that of other published boundary conditions and it has better stability behavior. © 1997 American Institute of Physics.
Show PACS
47.60.-i Flow phenomena in quasi-one-dimensional systems
47.11.-j Computational methods in fluid dynamics
47.70.Nd Nonequilibrium gas dynamics

Apparent fluid slip at hydrophobic microchannel walls

Derek C. Tretheway and Carl D. Meinhart

Phys. Fluids 14, L9 (2002); http://dx.doi.org/10.1063/1.1432696 (4 pages) | Cited 224 times

Full Text: | Download PDF

Show Abstract
Micron-resolution particle image velocimetry is used to measure the velocity profiles of water flowing through 30×300 μm channels. The velocity profiles are measured to within 450 nm of the microchannel surface. When the surface is hydrophilic (uncoated glass), the measured velocity profiles are consistent with solutions of Stokes’ equation and the well-accepted no-slip boundary condition. However, when the microchannel surface is coated with a 2.3 nm thick monolayer of hydrophobic octadecyltrichlorosilane, an apparent velocity slip is measured just above the solid surface. This velocity is approximately 10% of the free-stream velocity and yields a slip length of approximately 1 μm. For this slip length, slip flow is negligible for length scales greater than 1 mm, but must be considered at the micro- and nano scales. © 2002 American Institute of Physics.
Show PACS
47.45.Gx Slip flows and accommodation
47.60.-i Flow phenomena in quasi-one-dimensional systems

Scaling of temperature‐ and stress‐dependent viscosity convection

V. S. Solomatov

Phys. Fluids 7, 266 (1995); http://dx.doi.org/10.1063/1.868624 (9 pages) | Cited 213 times

Full Text: | Download PDF

Show Abstract
Simple scaling analysis of temperature‐ and stress‐dependent viscosity convection with free‐slip boundaries suggests three convective regimes: the small viscosity contrast regime which is similar to convection in a fluid whose viscosity does not depend on temperature, the transitional regime characterized by self‐controlled dynamics of the cold boundary layer and the asymptotic regime in which the cold boundary becomes stagnant and convection involves only the hottest part of the lid determined by a rheological temperature scale. The first two regimes are usually observed in numerical experiments. The last regime is similar to strongly temperature‐dependent viscosity convection with rigid boundaries studied in laboratory experiments. © 1995 American Institute of Physics.
Show PACS
47.27.T- Turbulent transport processes
47.50.-d Non-Newtonian fluid flows
91.35.Dc Heat flow; geothermy

Electrospinning and electrically forced jets. II. Applications

Moses M. Hohman, Michael Shin, Gregory Rutledge, and Michael P. Brenner

Phys. Fluids 13, 2221 (2001); http://dx.doi.org/10.1063/1.1384013 (16 pages) | Cited 201 times

Full Text: | Download PDF

Show Abstract
Electrospinning is a process in which solid fibers are produced from a polymeric fluid stream (solution or melt) delivered through a millimeter-scale nozzle. This article uses the stability theory described in the previous article to develop a quantitative method for predicting when electrospinning occurs. First a method for calculating the shape and charge density of a steady jet as it thins from the nozzle is presented and is shown to capture quantitative features of the experiments. Then, this information is combined with the stability analysis to predict scaling laws for the jet behavior and to produce operating diagrams for when electrospinning occurs, both as a function of experimental parameters. Predictions for how the regime of electrospinning changes as a function of the fluid conductivity and viscosity are presented. © 2001 American Institute of Physics.
Show PACS
47.65.-d Magnetohydrodynamics and electrohydrodynamics
83.50.-v Deformation and flow
47.50.-d Non-Newtonian fluid flows
83.80.Rs Polymer solutions
83.80.Sg Polymer melts
47.60.-i Flow phenomena in quasi-one-dimensional systems
66.20.-d Viscosity of liquids; diffusive momentum transport
47.20.-k Flow instabilities
83.60.Wc Flow instabilities

Scaling laws in granular flows down rough inclined planes

O. Pouliquen

Phys. Fluids 11, 542 (1999); http://dx.doi.org/10.1063/1.869928 (7 pages) | Cited 196 times

Full Text: | Download PDF

Show Abstract
In this paper, new scaling properties for granular flows down rough inclined planes are presented. In the dense steady uniform flow regime, we have systematically measured the mean velocity of the flow as a function of the inclination of the surface θ and of the thickness h of the layer. The results obtained for different systems of beads corresponding to different surface roughness conditions are shown to collapse into a single curve when properly scaled. The scaling is based on the measurement of the minimum thickness hstop(θ) necessary to observe a steady uniform flow at inclination θ. From this experimental observation an empirical description for granular flows down inclined planes is proposed in terms of a dynamic friction coefficient. © 1999 American Institute of Physics.
Show PACS
45.70.Mg Granular flow: mixing, segregation and stratification
47.55.Kf Particle-laden flows
47.10.-g General theory in fluid dynamics

Numerical investigation of 2D convection with extremely large viscosity variations

L.‐N. Moresi and V. S. Solomatov

Phys. Fluids 7, 2154 (1995); http://dx.doi.org/10.1063/1.868465 (9 pages) | Cited 191 times

Full Text: | Download PDF

Show Abstract
Previous experimental studies of convection in fluids with temperature‐dependent viscosity reached viscosity contrasts of the order of 105. Although this value seems large, it still might not be large enough for understanding convection in the interiors of Earth and other planets whose viscosity is a much stronger function of temperature. The reason is that, according to theory, above 104–105 viscosity contrasts, convection must undergo a major transition—to stagnant lid convection. This is an asymptotic regime in which a stagnant lid is formed on the top of the layer and convection is driven by the intrinsic, rheological, temperature scale, rather than by the entire temperature drop in the layer. A finite element multigrid scheme appropriate for large viscosity variations is employed and convection with up to 1014 viscosity contrasts has been systematically investigated in a 2D square cell with free‐slip boundaries. We reached the asymptotic regime in the limit of large viscosity contrasts and obtained scaling relations which are found to be in good agreement with theoretical predictions. © 1995 American Institute of Physics.
Show PACS
44.25.+f Natural convection
47.20.Bp Buoyancy-driven instabilities (e.g., Rayleigh-Benard)

The formation and evolution of synthetic jets

Barton L. Smith and Ari Glezer

Phys. Fluids 10, 2281 (1998); http://dx.doi.org/10.1063/1.869828 (17 pages) | Cited 183 times

Full Text: | Download PDF

Show Abstract
A nominally plane turbulent jet is synthesized by the interactions of a train of counter-rotating vortex pairs that are formed at the edge of an orifice by the time-periodic motion of a flexible diaphragm in a sealed cavity. Even though the jet is formed without net mass injection, the hydrodynamic impulse of the ejected fluid and thus the momentum of the ensuing jet are nonzero. Successive vortex pairs are not subjected to pairing or other subharmonic interactions. Each vortex of the pair develops a spanwise instability and ultimately undergoes transition to turbulence, slows down, loses its coherence and becomes indistinguishable from the mean jet flow. The trajectories of vortex pairs at a given formation frequency scale with the length of the ejected fluid slug regardless of the magnitude of the formation impulse and, near the jet exit plane, their celerity decreases monotonically with streamwise distance while the local mean velocity of the ensuing jet increases. In the far field, the synthetic jet is similar to conventional 2D jets in that cross-stream distributions of the time-averaged velocity and the corresponding rms fluctuations appear to collapse when plotted in the usual similarity coordinates. However, compared to conventional 2D jets, the streamwise decrease of the mean centerline velocity of the synthetic jet is somewhat higher ( ∼ x−0.58), and the streamwise increase of its width and volume flow rate is lower ( ∼ x0.88 and  ∼ x0.33, respectively). This departure from conventional self-similarity is consistent with the streamwise decrease in the jet’s momentum flux as a result of an adverse streamwise pressure gradient near its orifice. © 1998 American Institute of Physics.
Show PACS
47.27.wg Turbulent jets
47.32.C- Vortex dynamics
47.15.Fe Stability of laminar flows
47.27.Cn Transition to turbulence

On the universality of the Kolmogorov constant

Katepalli R. Sreenivasan

Phys. Fluids 7, 2778 (1995); http://dx.doi.org/10.1063/1.868656 (7 pages) | Cited 175 times

Full Text: | Download PDF

Show Abstract
All known data are collected on the Kolmogorov constant in one‐dimensional spectral formula for the inertial range. For large enough microscale Reynolds numbers, the data (despite much scatter) support the notion of a ‘‘universal’’ constant that is independent of the flow as well as the Reynolds number, with a numerical value of about 0.5. In particular, it is difficult to discern support for a recent claim that the constant is Reynolds number dependent even at high Reynolds numbers. © 1995 American Institute of Physics.
Show PACS
47.27.Ak Fundamentals
47.27.Gs Isotropic turbulence; homogeneous turbulence
47.27.Jv High-Reynolds-number turbulence

On a self-sustaining process in shear flows

Fabian Waleffe

Phys. Fluids 9, 883 (1997); http://dx.doi.org/10.1063/1.869185 (18 pages) | Cited 167 times

Full Text: | Download PDF

Show Abstract
A self-sustaining process conjectured to be generic for wall-bounded shear flows is investigated. The self-sustaining process consists of streamwise rolls that redistribute the mean shear to create streaks that wiggle to maintain the rolls. The process is analyzed and shown to be remarkably insensitive to whether there is no-slip or free-slip at the walls. A low-order model of the process is derived from the Navier–Stokes equations for a sinusoidal shear flow. The model has two unstable steady solutions above a critical Reynolds number, in addition to the stable laminar flow. For some parameter values, there is a second critical Reynolds number at which a homoclinic bifurcation gives rise to a stable periodic solution. This suggests a direct link between unstable steady solutions and almost periodic solutions that have been computed in plane Couette flow. It is argued that this self-sustaining process is responsible for the bifurcation of shear flows at low Reynolds numbers and perhaps also for controlling the near-wall region of turbulent shear flows at higher Reynolds numbers. © 1997 American Institute of Physics.
Show PACS
47.15.Fe Stability of laminar flows
47.20.Ft Instability of shear flows (e.g., Kelvin-Helmholtz)
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking
47.27.nb Boundary layer turbulence
47.60.-i Flow phenomena in quasi-one-dimensional systems
47.45.Gx Slip flows and accommodation
47.10.-g General theory in fluid dynamics

A thickened flame model for large eddy simulations of turbulent premixed combustion

O. Colin, F. Ducros, D. Veynante, and T. Poinsot

Phys. Fluids 12, 1843 (2000); http://dx.doi.org/10.1063/1.870436 (21 pages) | Cited 165 times

Full Text: | Download PDF

Show Abstract
A subgrid scale model for large eddy simulations of turbulent premixed combustion is developed and validated. The approach is based on the concept of artificially thickened flames, keeping constant the laminar flame speed sl0. This thickening is simply achieved by decreasing the pre-exponential factor of the chemical Arrhenius law whereas the molecular diffusion is enhanced. When the flame is thickened, the combustion–turbulence interaction is affected and must be modeled. This point is investigated here using direct numerical simulations of flame–vortex interactions and an efficiency function E is introduced to incorporate thickening effects in the subgrid scale model. The input parameters in E are related to the subgrid scale turbulence (velocity and length scales). An efficient approach, based on similarity assumptions, is developed to extract these quantities from the resolved velocity field. A specific operator is developed to exclude the dilatational part of the velocity field from the estimation of turbulent fluctuations. The combustion model is then implemented in a compressible parallel finite volume–element solver able to handle hybrid grids to simulate a lateral injections combustor (LIC). Results are in agreement with the available experimental data. © 2000 American Institute of Physics.
Show PACS
52.35.Bj Magnetohydrodynamic waves (e.g., Alfven waves)
52.25.Fi Transport properties
52.30.-q Plasma dynamics and flow
52.55.Fa Tokamaks, spherical tokamaks
52.75.-d Plasma devices

An approximate deconvolution procedure for large-eddy simulation

S. Stolz and N. A. Adams

Phys. Fluids 11, 1699 (1999); http://dx.doi.org/10.1063/1.869867 (3 pages) | Cited 164 times

Full Text: | Download PDF

Show Abstract
An alternative approach to large-eddy simulation based on approximate deconvolution (ADM) is developed. The main ingredient is an approximation of the nonfiltered field by truncated series expansion of the inverse filter operator. A posteriori tests for decaying compressible isotropic turbulence show excellent agreement with direct numerical simulation. The computational overhead of ADM is similar to that of a scale-similarity model and considerably less than for dynamic models. © 1999 American Institute of Physics.
Show PACS
47.11.-j Computational methods in fluid dynamics
47.27.E- Turbulence simulation and modeling
47.40.-x Compressible flows; shock waves

Momentum transfer of a Boltzmann-lattice fluid with boundaries

M’hamed Bouzidi, Mouaouia Firdaouss, and Pierre Lallemand

Phys. Fluids 13, 3452 (2001); http://dx.doi.org/10.1063/1.1399290 (8 pages) | Cited 159 times

Full Text: | Download PDF

Show Abstract
We study the velocity boundary condition for curved boundaries in the lattice Boltzmann equation (LBE). We propose a LBE boundary condition for moving boundaries by combination of the “bounce-back” scheme and spatial interpolations of first or second order. The proposed boundary condition is a simple, robust, efficient, and accurate scheme. Second-order accuracy of the boundary condition is demonstrated for two cases: (1) time-dependent two-dimensional circular Couette flow and (2) two-dimensional steady flow past a periodic array of circular cylinders (flow through the porous media of cylinders). For the former case, the lattice Boltzmann solution is compared with the analytic solution of the Navier–Stokes equation. For the latter case, the lattice Boltzmann solution is compared with a finite-element solution of the Navier–Stokes equation. The lattice Boltzmann solutions for both flows agree very well with the solutions of the Navier–Stokes equations. We also analyze the torque due to the momentum transfer between the fluid and the boundary for two initial conditions: (a) impulsively started cylinder and the fluid at rest, and (b) uniformly rotating fluid and the cylinder at rest. © 2001 American Institute of Physics.
Show PACS
47.11.-j Computational methods in fluid dynamics
05.60.-k Transport processes
05.50.+q Lattice theory and statistics (Ising, Potts, etc.)
02.60.Ed Interpolation; curve fitting
47.15.-x Laminar flows
47.56.+r Flows through porous media

The development of chaotic advection

Hassan Aref

Phys. Fluids 14, 1315 (2002); http://dx.doi.org/10.1063/1.1458932 (11 pages) | Cited 158 times

Online Publication Date: 28 February 2002

Full Text: | Download PDF

Show Abstract
The concept of chaotic advection was developed some twenty years ago as an outgrowth of work on advection by interacting point vortices. The term “chaotic advection” was first introduced in the title of an abstract for the 35th annual meeting of the American Physical Society (APS) Division of Fluid Dynamics (DFD) in 1982. The main reference, a Journal of Fluid Mechanics paper published in 1984, may be the true “birthdate” of the term. Earlier work from the 1960s by Arnol’d and Hénon on advection by steady three-dimensional flows already contained closely related ideas and results but was not widely appreciated. The present paper, based on the 2000 Otto Laporte Memorial Lecture delivered at the 53rd APS/DFD annual meeting, traces these and other precursors and the development of chaotic advection over the past two decades. Some exciting recent developments, such as application to fluid mixing in micro-electromechanical systems (MEMS) and to materials processing, and the introduction of topological methods of analysis, are highlighted. On balance, chaotic advection is now established as a subtopic of fluid mechanics with wide ramifications and continued promise for theory, experiment and applications. © 2002 American Institute of Physics.
Show PACS
47.52.+j Chaos in fluid dynamics
47.27.T- Turbulent transport processes
01.30.Rr Surveys and tutorial papers; resource letters
05.45.-a Nonlinear dynamics and chaos

Capillary effects during droplet impact on a solid surface

M. Pasandideh‐Fard, Y. M. Qiao, S. Chandra, and J. Mostaghimi

Phys. Fluids 8, 650 (1996); http://dx.doi.org/10.1063/1.868850 (10 pages) | Cited 158 times

Full Text: | Download PDF

Show Abstract
Impact of water droplets on a flat, solid surface was studied using both experiments and numerical simulation. Liquid–solid contact angle was varied in experiments by adding traces of a surfactant to water. Impacting droplets were photographed and liquid–solid contact diameters and contact angles were measured from photographs. A numerical solution of the Navier–Stokes equation using a modified SOLA‐VOF method was used to model droplet deformation. Measured values of dynamic contact angles were used as a boundary condition for the numerical model. Impacting droplets spread on the surface until liquid surface tension and viscosity overcame inertial forces, after which they recoiled off the surface. Adding a surfactant did not affect droplet shape during the initial stages of impact, but did increase maximum spread diameter and reduce recoil height. Comparison of computer generated images of impacting droplets with photographs showed that the numerical model modeled droplet shape evolution correctly. Accurate predictions were obtained for droplet contact diameter during spreading and at equilibrium. The model overpredicted droplet contact diameters during recoil. Assuming that dynamic surface tension of surfactant solutions is constant, equaling that of pure water, gave predicted droplet shapes that best agreed with experimental observations. When the contact angle was assumed constant in the model, equal to the measured equilibrium value, predictions were less accurate. A simple analytical model was developed to predict maximum droplet diameter after impact. Model predictions agreed well with experimental measurements reported in the literature. Capillary effects were shown to be negligible during droplet impact when We≫Re1/2. © 1996 American Institute of Physics.
Show PACS
47.55.D- Drops and bubbles

Measurements of the collision properties of small spheres

Samuel F. Foerster, Michel Y. Louge, Hongder Chang, and Khédidja Allia

Phys. Fluids 6, 1108 (1994); http://dx.doi.org/10.1063/1.868282 (8 pages) | Cited 150 times

Full Text: | Download PDF

Show Abstract
An experiment to measure the properties of the collisions between two small spheres or between a small sphere and a semi‐infinite flat wall are described. The apparatus releases the particles in a free‐fall without initial spin. The impacts are modeled in terms of three coefficients. The first is the coefficient of normal restitution. The second represents the frictional properties of the contact surfaces. The last characterizes the restitution of the tangential components of the velocity of the contact point for impacts that do not involve sliding. The coefficients are calculated from stroboscopic photographs of the ballistics of the particles near the collision. The results establish that the collision model provides an accurate description of the dynamics of the impacts.
Show PACS
81.05.Rm Porous materials; granular materials
45.05.+x General theory of classical mechanics of discrete systems
83.80.Fg Granular solids

Effects of heat release on triple flames

G. R. Ruetsch, L. Vervisch, and A. Liñán

Phys. Fluids 7, 1447 (1995); http://dx.doi.org/10.1063/1.868531 (8 pages) | Cited 149 times

Full Text: | Download PDF

Show Abstract
Heat release effects on laminar flame propagation in partially premixed flows are studied. Data for analysis are obtained from direct numerical simulations of a laminar mixing layer with a uniformly approaching velocity field. The structure that evolves under such conditions is a triple flame, which consists of two premixed wings and a trailing diffusion flame. Heat release increases the flame speed over that of the corresponding planar premixed flame. In agreement with previous analytical work, reductions in the mixture fraction gradient also increase the flame speed. The effects of heat release and mixture fraction gradients on flame speed are not independent, however; heat release modifies the effective mixture fraction gradient in front of the flame. For very small mixture fraction gradients, scaling laws that determine the flame speed in terms of the density change are presented. © 1995 American Institute of Physics.
Show PACS
47.70.Fw Chemically reactive flows

Energy dissipation rate and energy spectrum in high resolution direct numerical simulations of turbulence in a periodic box

Yukio Kaneda, Takashi Ishihara, Mitsuo Yokokawa, Ken’ichi Itakura, and Atsuya Uno

Phys. Fluids 15, L21 (2003); http://dx.doi.org/10.1063/1.1539855 (4 pages) | Cited 146 times

Online Publication Date: 8 January 2003

Full Text: | Download PDF

Show Abstract
High-resolution direct numerical simulations (DNSs) of incompressible homogeneous turbulence in a periodic box with up to 40963 grid points were performed on the Earth Simulator computing system. DNS databases, including the present results, suggest that the normalized mean energy dissipation rate per unit mass tends to a constant, independent of the fluid kinematic viscosity ν as ν→0. The DNS results also suggest that the energy spectrum in the inertial subrange almost follows the Kolmogorov k−5/3 scaling law, where k is the wavenumber, but the exponent is steeper than −5/3 by about 0.1. © 2003 American Institute of Physics.
Show PACS
47.27.Gs Isotropic turbulence; homogeneous turbulence

Filtered density function for large eddy simulation of turbulent reacting flows

P. J. Colucci, F. A. Jaberi, P. Givi, and S. B. Pope

Phys. Fluids 10, 499 (1998); http://dx.doi.org/10.1063/1.869537 (17 pages) | Cited 143 times

Full Text: | Download PDF

Show Abstract
A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.
Show PACS
47.27.-i Turbulent flows
47.70.Fw Chemically reactive flows
47.11.-j Computational methods in fluid dynamics
02.50.Ng Distribution theory and Monte Carlo studies
Close
Google Calendar
ADVERTISEMENT

close