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

Flickr Twitter iResearch App Facebook

Search Issue | RSS Feeds RSS
Previous Issue Next Issue

Jul 2000

Volume 12, Issue 7, pp. 1629-1881

back to top
RSS Feeds

An approach to wall modeling in large-eddy simulations

N. V. Nikitin, F. Nicoud, B. Wasistho, K. D. Squires, and P. R. Spalart

Phys. Fluids 12, 1629 (2000); http://dx.doi.org/10.1063/1.870414 (4 pages) | Cited 97 times

Full Text: | Download PDF

Show Abstract
Channel flow with friction Reynolds number Reτ as high as 80 000 is treated by large-eddy simulation at a moderate cost, using the subgrid-scale model designed for detached-eddy simulations. It includes wall modeling, and was not adjusted for this flow. The grid count scales with the logarithm of the Reynolds number. Three independent codes are in fair agreement with each other. Reynolds-number variations and grid refinement cause trades between viscous, modeled, and resolved shear stresses. The skin-friction coefficient is too low, on the order of 15%. The velocity profiles contain a “modeled” logarithmic layer near the wall and some suggest a “resolved” logarithmic layer farther up, but the two layers have a mismatch of several units in U+. © 2000 American Institute of Physics.
Show PACS
47.60.-i Flow phenomena in quasi-one-dimensional systems
47.11.-j Computational methods in fluid dynamics
back to top
RSS Feeds

Nonlinear dynamics of three-dimensional long-wave Marangoni instability in thin liquid films

Alexander Oron

Phys. Fluids 12, 1633 (2000); http://dx.doi.org/10.1063/1.870415 (13 pages) | Cited 29 times

Full Text: | Download PDF

Show Abstract
The three-dimensional evolution of the long-wave Marangoni instability of thin liquid films is studied. According to earlier theoretical predictions no continuous steady states exist and the film ruptures. As in two dimensions, the mechanism of fingering is found to be a main route to rupture. A four-fold rotational symmetry of the film interface is retained, when a square periodic domain and the harmonic initial disturbance are used. A use of initial random disturbances in general eliminates the square symmetry of the solution. An increase of the domain size results in a growing complexity of the emerging patterns. In contrast with the dynamics in two dimensions the evolution of the interface in three dimensions and in particular the pattern emerging at rupture may strongly depend on the choice of the initial condition. The two-dimensional evolution of the film is found to be unstable to small three-dimensional random disturbances. © 2000 American Institute of Physics.
Show PACS
47.20.Dr Surface-tension-driven instability
47.27.T- Turbulent transport processes

Spreading and imbibition of viscous liquid on a porous base. II

S. H. Davis and L. M. Hocking

Phys. Fluids 12, 1646 (2000); http://dx.doi.org/10.1063/1.870416 (10 pages) | Cited 24 times

Full Text: | Download PDF

Show Abstract
The competition between viscous spreading of a liquid on a substrate and its absorption by a dry porous substrate is studied in two dimensions. The effect of capillary suction into the pores vies with contact-line slip on the substrate to determine the lifetime of the drop. © 2000 American Institute of Physics.
Show PACS
47.56.+r Flows through porous media
47.55.D- Drops and bubbles
68.08.Bc Wetting

Optimization of fluid front dynamics in porous media using rate control. I. Equal mobility fluids

Bagus Sudaryanto and Yannis C. Yortsos

Phys. Fluids 12, 1656 (2000); http://dx.doi.org/10.1063/1.870417 (15 pages) | Cited 9 times

Full Text: | Download PDF

Show Abstract
In applications involving the injection of a fluid in a porous medium to displace another fluid, a main objective is the maximization of the displacement efficiency. For a fixed arrangement of injection and production points (sources and sinks), such optimization is possible by controling the injection rate policy. Despite its practical relevance, however, this aspect has received scant attention in the literature. In this paper, we provide a fundamental approach based on optimal control theory, for the simplified case when the fluids are miscible, of equal viscosity, and in the absence of dispersion and gravity effects. Both homogeneous and heterogeneous porous media are considered. From a fluid dynamics viewpoint, this is a problem in the deformation of material lines in porous media, as a function of time-varying injection rates. It is shown that the optimal injection policy that maximizes the displacement efficiency, at the time of arrival of the injected fluid, is of the “bang–bang” type, in which the rates take their extreme values in the range allowed. This result applies to both homogeneous and heterogeneous media. Examples in simple geometries and for various constraints are shown, illustrating the efficiency improvement over the conventional approach of constant rate injection. In the heterogeneous case, the effect of the permeability heterogeneity, particularly its spatial correlation structure, on diverting the flow paths, is analyzed. It is shown that bang–bang injection remains the optimal approach, compared to constant rate, particularly if they were both designed under the assumption that the medium was homogeneous. Experiments in a homogeneous Hele-Shaw cell are found to be in good agreement with the theory. © 2000 American Institute of Physics.
Show PACS
47.56.+r Flows through porous media
47.85.L- Flow control

Electrohydrodynamic effects on the deformation and orientation of a liquid capsule in a linear flow

Jong-Wook Ha and Seung-Man Yang

Phys. Fluids 12, 1671 (2000); http://dx.doi.org/10.1063/1.870418 (14 pages) | Cited 3 times

Full Text: | Download PDF

Show Abstract
The role of a uniform electric field on the deformation and orientation of a liquid capsule with a viscoelastic membrane is considered analytically in the small deformation limit. The capsule is freely suspended either in a quiescent fluid or in a shear flow. The viscoelasticity of the membrane is taken into account by the Kelvin–Voigt model and the electrohydrodynamic flow is analyzed on the basis of the leaky dielectric model. In this article, we consider three different prototype models of capsules; viz., a neo-Hookean (incompressible isotropic) membrane, a red blood cell-type (area-preserving) membrane, and an interfacial-tension droplet. The deformed capsule shape from its initial sphericity and its orientaion are determined from the linearized governing equations and boundary conditions in the limit of small deformations. The asymptotic theory shows that the degree of capsule deformation induced by a uniform electric field alone is independent of the surface viscosity of the capsule as well as the viscosity ratio between the two fluids inside and outside the capsule. Meanwhile, in the presence of an imposed shear flow, the degree of deformation depends on the surface viscosity with preserving still the independence of the viscosity ratio. For an illustrative purpose, experimental results for the role of a uniform electric field on the orientation of an interfacial-tension droplet in a shear flow are discussed briefly. © 2000 American Institute of Physics.
Show PACS
47.65.-d Magnetohydrodynamics and electrohydrodynamics
83.60.Bc Linear viscoelasticity
87.19.U- Hemodynamics
87.19.Wx Pneumodyamics, respiration

Interfacial instabilities of a fluid annulus in a rotating Hele–Shaw cell

Lluís Carrillo, Jordi Soriano, and Jordi Ortín

Phys. Fluids 12, 1685 (2000); http://dx.doi.org/10.1063/1.870419 (14 pages) | Cited 21 times

Full Text: | Download PDF

Show Abstract
We have studied the interfacial instabilities experienced by a liquid annulus as it moves radially in a circular Hele–Shaw cell rotating with angular velocity Ω. The instability of the leading interface (oil displacing air) is driven by the density difference in the presence of centrifugal forcing, while the instability of the trailing interface (air displacing oil) is driven by the large viscosity contrast. A linear stability analysis shows that the stability of the two interfaces is coupled through the pressure field already at a linear level. We have performed experiments in a dry cell and in a cell coated with a thin fluid layer on each plate, and found that the stability depends substantially on the wetting conditions at the leading interface. Our experimental results of the number of fingers resulting from the instability compare well with the predictions obtained through a numerical integration of the coupled equations derived from a linear stability analysis. Deep in the nonlinear regime we observe the emission of liquid droplets through the formation of thin filaments at the tip of outgrowing fingers. © 2000 American Institute of Physics.
Show PACS
47.20.Gv Viscous and viscoelastic instabilities
81.15.Lm Liquid phase epitaxy; deposition from liquid phases (melts, solutions, and surface layers on liquids)
68.03.-g Gas-liquid and vacuum-liquid interfaces
47.32.-y Vortex dynamics; rotating fluids
47.55.D- Drops and bubbles

On the rapid estimation of permeability for porous media using Brownian motion paths

Chi-Ok Hwang, James A. Given, and Michael Mascagni

Phys. Fluids 12, 1699 (2000); http://dx.doi.org/10.1063/1.870420 (11 pages) | Cited 8 times

Full Text: | Download PDF

Show Abstract
We describe two efficient methods of estimating the fluid permeability of standard models of porous media by using the statistics of continuous Brownian motion paths that initiate outside a sample and terminate on contacting the porous sample. The first method associates the “penetration depth” with a specific property of the Brownian paths, then uses the standard relation between penetration depth and permeability to calculate the latter. The second method uses Brownian paths to calculate an effective capacitance for the sample, then relates the capacitance, via angle-averaging theorems, to the translational hydrodynamic friction of the sample. Finally, a result of Felderhof is used to relate the latter quantity to the permeability of the sample. We find that the penetration depth method is highly accurate in predicting permeability of porous material. © 2000 American Institute of Physics.
Show PACS
47.56.+r Flows through porous media
05.40.Jc Brownian motion

Mechanism of air entrainment by a disturbed liquid jet

C. D. Ohl, H. N. Oguz, and A. Prosperetti

Phys. Fluids 12, 1710 (2000); http://dx.doi.org/10.1063/1.870421 (5 pages) | Cited 17 times

Full Text: | Download PDF

Show Abstract
It was shown in recent work that the crests of surface disturbances on a falling jet are a powerful agent for air entrainment at the free surface of a liquid pool. The paper explores the opposite case in which the jet is disturbed so as to form an axisymmetric trough, rather than a crest. It is found that no air is entrained in this case. The paper concludes with some considerations on the validity of a recently proposed model for air entrainment. © 2000 American Institute of Physics.
Show PACS
47.55.Hd Stratified flows
47.55.Kf Particle-laden flows
47.27.wg Turbulent jets

Weakly nonlinear saturation of short-wave instabilities in a strained Lamb–Oseen vortex

Denis Sipp

Phys. Fluids 12, 1715 (2000); http://dx.doi.org/10.1063/1.870422 (15 pages) | Cited 15 times

Full Text: | Download PDF

Show Abstract
A Lamb–Oseen vortex in a planar straining field is known to be subject to 3D (three-dimensional) short-wave instabilities which are due to the resonance of the straining field and two stationary Kelvin waves characterized by the same axial wave number and by azimuthal wave numbers equal to −1 and +1. The linear regime has been described by Moore and Saffman (1975). In this article, we extend this analysis to the weakly nonlinear regime. The emerging eigenmode is characterized by a complex amplitude A = ∣Aeiϕ, whose behavior is governed by an amplitude equation. It is shown that the unstable perturbation corresponds to an oscillation of the vortex in a plane inclined at an angle ϕ, while the amplitude of these oscillations is proportional to A∣. The vortex centers are defined as the points where the velocity of the vortex is zero, which also corresponds to the points where the pressure is minimum. We show that these instabilities saturate. The saturation amplitudes are evaluated numerically and expressed in terms of oscillation amplitudes of the vortex centers. If a denotes the internal radius of the vortex and if the straining field is due to a counter-rotating vortex of same strength, located at a distance b, then the maximum amplitude Δ of the vortex oscillations is Δ/b = 6.1a2/b2. This result is in agreement with those of the experiments of Leweke and Williamson (1998) for which a/b = 0.2. It also shows that in aeronautical situations, for which a/b is smaller, i.e., a/b<0.1, the considered short-wave instability will saturate at very low amplitude. © 2000 American Institute of Physics.
Show PACS
47.32.C- Vortex dynamics
47.20.Ky Nonlinearity, bifurcation, and symmetry breaking
47.35.-i Hydrodynamic waves

Axisymmetric instabilities of Bödewadt flow

R. Fernandez-Feria

Phys. Fluids 12, 1730 (2000); http://dx.doi.org/10.1063/1.870423 (10 pages) | Cited 3 times

Full Text: | Download PDF

Show Abstract
A spatial linear stability analysis of Bödewadt’s self-similar solution for the rotating flow over a flat plate is performed. In particular, considered is the stability of axisymmetric perturbations propagating towards the axis of rotation, which are the most important ones observed experimentally. Viscous and nonparallel effects on the stability of the perturbations are retained up to the order of the inverse of the local Reynolds number R. The resulting parabolic stability equations are solved numerically using a spectral collocation method varying the nondimensional frequency q and R. The instability region on the (q,R)-plane is discussed and compared with existing experimental data and direct numerical simulation results. The circular waves observed experimentally and in numerical simulations are shown to correspond to an inertial instability mode which becomes stabilized as R decreases below a critical value. © 2000 American Institute of Physics.
Show PACS
47.32.-y Vortex dynamics; rotating fluids
47.20.-k Flow instabilities

Three-dimensional centrifugal-type instabilities of two-dimensional flows in rotating systems

Denis Sipp and Laurent Jacquin

Phys. Fluids 12, 1740 (2000); http://dx.doi.org/10.1063/1.870424 (9 pages) | Cited 25 times

Full Text: | Download PDF

Show Abstract
This paper deals with the stability of incompressible inviscid planar basic flows in a rotating frame. We give a sufficient condition for such flows to undergo three-dimensional shortwave centrifugal-type instabilities. This criterion reduces to the Bradshaw–Richardson (1969) or Pedley (1969) criterion in the specific case of parallel shear flows subject to rotation, to Rayleigh’s centrifugal criterion (1916) in the case of axisymmetric vortices in inertial frames, to the Kloosterziel and van Heijst (1991) criterion in the case of axisymmetric vortices subject to rotation and to Bayly’s criterion (1988) in the case of general two-dimensional flows in inertial frames. The criterion states that a steady 2D basic flow subject to rotation Ω is unstable if there exists a streamline for which at each point 2(V/R+Ω)(W+2Ω)<0 where W is the vorticity of the streamline, R is the local algebraic radius of curvature of the streamline and V is the local norm of the velocity. If this condition is satisfied then the flow is unstable according to the geometrical optics method introduced by Lifschitz and Hameiri (1991), which consists in following wave packets along the flow trajectories using a Wentzel–Kramers–Brillouin formalism. When the streamlines are closed, it is further shown that a localized unstable normal mode can be constructed in the vicinity of a streamline. As an application, this new criterion is used to study the centrifugal-type instabilities in the Stuart vortices, which is a family of exact solutions describing a row of periodic co-rotating eddies. For each solution of that family and for each rotation parameter f = 2Ω, we give the unstable streamline interval, according to the criterion of instability. This criterion gives only a sufficient condition of centrifugal instability. The equations of the geometrical optics method are therefore numerically solved to obtain the true centrifugally unstable streamline intervals. It turns out that our criterion gives excellent results for highly concentrated vortices, i.e., the two approaches yield the same unstable streamline intervals. In less concentrated vortices, some streamlines undergo centrifugal instability although our criterion is not fulfilled. From these numerical results, another criterion of centrifugal instability for a flow with closed streamlines is conjectured which reduces to the change of sign of the absolute vorticity W+2Ω somewhere in the flow. © 2000 American Institute of Physics.
Show PACS
47.20.Ft Instability of shear flows (e.g., Kelvin-Helmholtz)
47.32.C- Vortex dynamics
47.20.Cq Inviscid instability

Receptivity to harmonic forcing of a finite-width mixing layer just downstream of the trailing edge of a flat plate

J. A. Rabchuk

Phys. Fluids 12, 1749 (2000); http://dx.doi.org/10.1063/1.870425 (13 pages)

Full Text: | Download PDF

Show Abstract
Analytical solutions for the receptivity of a piecewise linear mixing layer to harmonic disturbances at the trailing edge of a flat plate are presented. The solution is for an ideal fluid, with viscous effects at the trailing edge being accounted for by imposing an unsteady Kutta condition. The wake from the flat plate is included in the model. Initial mixing layer response for the discrete and continuous modes of instability is determined as a function of forcing strength, Strouhal number, and cross stream location. Forcing frequencies throughout the unstable regime are considered. The various components of the perturbation field and their contributions to the perturbation velocity and vorticity are investigated in the near field of the trailing edge. In the limit that the forcing frequency approaches zero, the amplification of the streamwise velocity and vorticity associated with the growing instability mode is shown to be inversely proportional to the forcing frequency, in agreement with previous results obtained using the vortex sheet approximation. In contrast, the initial cross stream velocity of the disturbance is quite small for low frequencies, and reaches a peak at the most amplified frequency of the mixing layer profile. The analysis also shows that at higher frequencies the continuous modes of instability become increasingly important channels whereby vorticity is carried away from the trailing edge. The differences in wave number between the various competing components at these higher frequencies leads to the appearance of nodes in the time-averaged perturbation quantities downstream. The entire perturbation is examined at the most amplified frequency. © 2000 American Institute of Physics.
Show PACS
47.20.-k Flow instabilities
47.27.wb Turbulent wakes
47.32.C- Vortex dynamics

Firespread through fuel beds: Modeling of wind-aided fires and induced hydrodynamics

B. Porterie, D. Morvan, J. C. Loraud, and M. Larini

Phys. Fluids 12, 1762 (2000); http://dx.doi.org/10.1063/1.870426 (21 pages) | Cited 33 times

Full Text: | Download PDF

Show Abstract
The propagation of wind-aided line fires through fuel beds is simulated by using a multiphase approach. In this approach, a gas phase flows through N solid phases which constitute an idealized reproduction of the heterogeneous combustible medium. A set of time-dependent equations is obtained for each phase and the coupling between the gas phase and the solid phases is rendered through exchange terms of mass, momentum, and energy. Turbulence is approached by using a RNG kε statistical model constructed from the Favre averaging method. The radiative transfer equation extended to multiphase media is solved using the discrete ordinates method (DOM). Soot formation is taken into account for the evaluation of the absorption coefficient of the soot/fuel particles/combustion products mixture using the gray gas assumption. First-order kinetics is incorporated to describe water vaporization, pyrolysis, and char combustion processes. The solution is performed numerically by a finite-volume method including a high-order upwind convective scheme and a flux limiter strategy along with a projection method for the pressure–velocity coupling. This model has been applied to describe the unsteady behavior of wind-aided fires spreading through a litter of dead pine needles and the induced hydrodynamics. The numerical results obtained from our model are presented and compared to measurements and predictions from other laboratory-based models. © 2000 American Institute of Physics.
Show PACS
47.70.Fw Chemically reactive flows
82.33.Vx Reactions in flames, combustion, and explosions
47.27.E- Turbulence simulation and modeling
05.60.-k Transport processes
82.30.Lp Decomposition reactions (pyrolysis, dissociation, and fragmentation)

Two-dimensional Navier–Stokes simulations of gaseous mixtures induced by Richtmyer–Meshkov instability

Claude Mügler and Serge Gauthier

Phys. Fluids 12, 1783 (2000); http://dx.doi.org/10.1063/1.870427 (16 pages) | Cited 18 times

Full Text: | Download PDF

Show Abstract
Two-dimensional numerical calculations of the fluid instability of shock-accelerated interfaces between a heavy fluid and a light one are carried out in order to simulate experiments performed by Poggi et al. [Phys. Fluids 10, 2698 (1998)]. In these experiments, the laser Doppler anemometry technique gives measurements of the fluctuating velocity. Experimental data show that a turbulent mixing zone is generated by the incident shock wave. This turbulent regime is reproduced by two-dimensional calculations. Before the first reshock, several quantities in the mixing zone, such as bubble and spike fronts, turbulent kinetic energy, enstrophy, adopt a quasi self-similar behavior versus time. In particular, we can see in numerical simulations the decay of the turbulent kinetic energy before the first reflected shock wave–mixing-zone interaction and its strong enhancement by reshocks. Furthermore, spectral analysis of the numerical results exhibit a k−3 energy spectrum. Experimental measurements also show that the turbulent boundary layers which develop on the shock-tube walls accelerate the fluid flow in the middle of the tube. Numerical simulations clearly reproduce both this acceleration and the lambda-shock structure observed in experiments. © 2000 American Institute of Physics.
Show PACS
47.10.-g General theory in fluid dynamics

Receptivity of three-dimensional boundary-layers to localized wall roughness and suction

F. P. Bertolotti

Phys. Fluids 12, 1799 (2000); http://dx.doi.org/10.1063/1.870428 (11 pages) | Cited 8 times

Full Text: | Download PDF

Show Abstract
The receptivity of three-dimensional, incompressible, boundary-layer flows to localized roughness or suction is studied theoretically and numerically. Three-dimensional boundary-layer flows have a strong “nonparallel” component due to the curvature of the potential-flow streamlines. Our theoretical model extends the Fourier-transform methodology to nonparallel flows using a Taylor expansion of the laminar mean-flow at the location of the roughness. Both the near-field and the far-field solutions are contained in the model. Additionally, the use of MacLaurin expansions in the complex plane leads quickly to the receptivity factors of eigenmodes. The theoretical results are validated with solutions to the linearized Navier–Stokes equations efficiently obtained by replacing the actual wall forcing with a smoother, but equivalent, forcing. We find that the far-field response is proportional to both math(αe) and dmath(αe)/dα, where math(α) is the Fourier transform of the forcing distribution, and αe is the eigenmode’s wave number. The receptivity coefficient for math(αe) decreases when nonparallelism is included in our models, while that for dmath(αe)/dα increases. © 2000 American Institute of Physics.
Show PACS
47.15.Cb Laminar boundary layers
47.20.Pc Flow receptivity

A physical-space version of the stretched-vortex subgrid-stress model for large-eddy simulation

Tobias Voelkl, D. I. Pullin, and Daniel C. Chan

Phys. Fluids 12, 1810 (2000); http://dx.doi.org/10.1063/1.870429 (16 pages) | Cited 33 times

Full Text: | Download PDF

Show Abstract
A physical-space version of the stretched-vortex subgrid-stress model is presented and applied to large-eddy simulations of incompressible flows. This version estimates the subgrid-kinetic energy required for evaluation of the subgrid-stress tensor using local second-order structure-function information of the resolved velocity field at separations of order the local cell size. A relation between the structure function and the energy spectrum is derived using the kinematic assumptions of the stretched-vortex model for locally homogeneous anisotropic turbulence. Results of large-eddy simulations using this model are compared to experimental and direct numerical simulation data. Comparisons are shown for the decay of kinetic energy and energy spectra of decaying isotropic turbulence and for mean velocities, root-mean-square velocity fluctuations and turbulence-kinetic energy balances of channel flow at three different Reynolds numbers. © 2000 American Institute of Physics.
Show PACS
47.11.-j Computational methods in fluid dynamics
47.32.C- Vortex dynamics
47.27.-i Turbulent flows
05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion
47.60.-i Flow phenomena in quasi-one-dimensional systems

The simulation and interpretation of free turbulence with a cognitive neural system

Francesc Giralt, A. Arenas, J. Ferre-Giné, R. Rallo, and G. A. Kopp

Phys. Fluids 12, 1826 (2000); http://dx.doi.org/10.1063/1.870430 (10 pages) | Cited 9 times

Full Text: | Download PDF

Show Abstract
An artificial neural network, based on fuzzy ARTMAP, that is capable of learning the basic nonlinear dynamics of a turbulent velocity field is presented. The neural system is capable of generating a detailed multipoint time record with the same structural characteristics and basic statistics as those of the original instantaneous velocity field used for training. The good performance of the proposed architecture is demonstrated by the generation of synthetic two-dimensional velocity data at eight different positions along the homogeneous (spanwise) direction in the far region (x/D = 420) of a turbulent wake flow generated behind a cylinder at Re = 1 200. The analysis of the synthetic velocity field, carried out with spectral techniques, POD and pattern recognition, reveals that the proposed neural system is capable of capturing the highly nonlinear dynamics of free turbulence and of reproducing the sequence of individual classes of relevant events present in turbulent wake flows. The trained neural system also yields patterns of the coherent structures embedded in the flow when presented with input data containing partial information of the instantaneous velocity maps of these events. In this way, the neural network is used as an expert system that helps in the structural interpretation of turbulence in a wake flow. © 2000 American Institute of Physics.
Show PACS
47.27.E- Turbulence simulation and modeling
07.05.Mh Neural networks, fuzzy logic, artificial intelligence

Statistics of pressure and of pressure-velocity correlations in isotropic turbulence

L. Biferale, P. Gualtieri, and F. Toschi

Phys. Fluids 12, 1836 (2000); http://dx.doi.org/10.1063/1.870431 (7 pages) | Cited 2 times

Full Text: | Download PDF

Show Abstract
Some pressure and pressure-velocity correlations in a direct numerical simulations of a three-dimensional turbulent flow at moderate Reynolds numbers have been analyzed. We have identified a set of pressure-velocity correlations which possess a good scaling behavior. Such a class of pressure-velocity correlations is determined by looking at the energy-balance across any sub-volume of the flow. According to our analysis, pressure scaling is determined by the dimensional assumption that pressure behaves as a “velocity squared,” unless finite-Reynolds effects are overwhelming. The SO(3) decompositions of pressure structure functions has also been applied in order to investigate anisotropic effects on the pressure scaling. © 2000 American Institute of Physics.
Show PACS
47.27.Gs Isotropic turbulence; homogeneous turbulence

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 164 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

Demonstration of a one-way flow of a rarefied gas induced through a pipe without average pressure and temperature gradients

Yoshio Sone and Ko Sato

Phys. Fluids 12, 1864 (2000); http://dx.doi.org/10.1063/1.870432 (5 pages) | Cited 9 times

Full Text: | Download PDF

Show Abstract
A simple experiment is performed to demonstrate that a one-way flow of a rarefied gas can be induced through a pipe without average pressure and temperature gradients by devising the shape of the pipe. A one-way flow is found to be induced through a pipe consisting of two circular pipes of different diameters connected and heated at their joint in the middle part of the pipe. The flow is from the thinner pipe to the thicker. In a pipe of a uniform diameter, no one-way flow is induced even when it is heated at a position away from the middle of the pipe. © 2000 American Institute of Physics.
Show PACS
47.45.-n Rarefied gas dynamics
47.60.-i Flow phenomena in quasi-one-dimensional systems

Non-self-similar behavior of the von Neumann reflection

S. Kobayashi, T. Adachi, and T. Suzuki

Phys. Fluids 12, 1869 (2000); http://dx.doi.org/10.1063/1.870433 (9 pages) | Cited 4 times

Full Text: | Download PDF

Show Abstract
The purpose of the present paper is to consider the von Neumann reflection (vNR) which takes place for a small incident shock Mach number and a small wedge angle. A series of experiments has been performed with ordinary smooth straight wedges and step-like wedges that simulate the former. The reflection configuration over the step-like wedge has suggested unsteady characteristics of the vNR. Contrary to the established notion of shock reflection phenomena over straight wedges, while the triple point trajectory was approximately a straight line through the apex of the wedge, the vNR showed unsteadiness in the relation between angles of incidence and reflection (ωi,ωr), and thus the flow-field proved to be non-self-similar near the triple point. Based on the measured values, the incident angle for the reflected wave and the Mach number of the flow ahead of the reflected wave were estimated. These values show that the reflected wave is not a Mach wave, but it moves on the (ωi,ωr)-plane almost along a trivial solution. In particular, the flow Mach number ahead of the reflected wave approaches unity, which leads to analytical formulas for angles of incidence and reflection as functions of the incident shock Mach number only. The reflected wave degenerates to a normal Mach wave as the incident shock proceeds. In the case of the vNR, the von Neumann paradox takes place only at an early stage of reflection, and the paradox is resolved later, since the flow properties such as angles of incidence and reflection are given by the trivial solution which is a particular solution for the three-shock theory. © 2000 American Institute of Physics.
Show PACS
47.40.Nm Shock wave interactions and shock effects
47.80.-v Instrumentation and measurement methods in fluid dynamics
47.40.-x Compressible flows; shock waves
back to top
RSS Feeds

On the stability of a stationary solution of the Tchen’s equation

O. A. Druzhinin

Phys. Fluids 12, 1878 (2000); http://dx.doi.org/10.1063/1.870434 (3 pages) | Cited 3 times

Full Text: | Download PDF

Show Abstract
This note examines the stability of a stationary velocity solution of the Tchen’s equation, describing the motion of a spherical particle with a Reynolds number much less than unity in a uniform unsteady fluid flow. The results show that the stationary solution is unstable if the ratio of the particle and fluid densities ρp/ρf<7/4. Therefore, applying the stationary solution is not justified in this case. © 2000 American Institute of Physics.
Show PACS
47.55.Kf Particle-laden flows
47.20.-k Flow instabilities
back to top
RSS Feeds
FREE

Erratum: “The motion of a falling liquid filament” [Phys. Fluids 12, 550 (2000)]

Diane Henderson, Harvey Segur, Linda B. Smolka, and Miki Wadati

Phys. Fluids 12, 1881 (2000); http://dx.doi.org/10.1063/1.870435 (1 page)

Full Text: | Download PDF

Abstract Unavailable
Show PACS
99.10.Cd Errata
47.55.D- Drops and bubbles
47.20.Cq Inviscid instability
47.20.Gv Viscous and viscoelastic instabilities
47.27.wg Turbulent jets
47.10.-g General theory in fluid dynamics
Close
Google Calendar
ADVERTISEMENT

close