scispace - formally typeset
Search or ask a question

Showing papers on "K-epsilon turbulence model published in 2008"


Journal ArticleDOI
TL;DR: In this article, a reformulated version of the author's k-ω model of turbulence has been presented, which has been applied to both boundary layers and free shear flows and has little sensitivity to finite freestream boundary conditions on turbulence properties.
Abstract: This paper presents a reformulated version of the author'sk-ω model of turbulence. Revisions include the addition of just one new closure coefficient and an adjustment to the dependence of eddy viscosity on turbulence properties. The result is a significantly improved model that applies to both boundary layers and free shear flows and that has very little sensitivity to finite freestream boundary conditions on turbulence properties. The improvements to the k-ω model facilitate a significant expansion of its range of applicability. The new model, like preceding versions, provides accurate solutions for mildly separated flows and simple geometries such as that of a backward-facing step. The model's improvement over earlier versions lies in its accuracy for even more complicated separated flows. This paper demonstrates the enhanced capability for supersonic flow into compression corners and a hypersonic shock-wave/ boundary-layer interaction. The excellent agreement is achieved without introducing any compressibility modifications to the turbulence model.

882 citations


Journal ArticleDOI
TL;DR: In this article, an eddy-viscosity turbulence model employing three additional transport equations is presented and applied to a number of transitional flow test cases, which is based on the k- framework and represents a substantial refinement to a transition-sensitive model that has been previously documented in the open literature.
Abstract: An eddy-viscosity turbulence model employing three additional transport equations is presented and applied to a number of transitional flow test cases. The model is based on the k- framework and represents a substantial refinement to a transition-sensitive model that has been previously documented in the open literature. The third transport equation is included to predict the magnitude of low-frequency velocity fluctuations in the pretransitional boundary layer that have been identified as the precursors to transition. The closure of model terms is based on a phenomenological (i.e., physics-based) rather than a purely empirical approach and the rationale for the forms of these terms is discussed. The model has been implemented into a commercial computational fluid dynamics code and applied to a number of relevant test cases, including flat plate boundary layers with and without applied pressure gradients, as well as a variety of airfoil test cases with different geometries, Reynolds numbers, freestream turbulence conditions, and angles of attack. The test cases demonstrate the ability of the model to successfully reproduce transitional flow behavior with a reasonable degree of accuracy, particularly in comparison with commonly used models that exhibit no capability of predicting laminar-toturbulent boundary layer development. While it is impossible to resolve all of the complex features of transitional and turbulent flows with a relatively simple Reynolds-averaged modeling approach, the results shown here demonstrate that the new model can provide a useful and practical tool for engineers addressing the simulation and prediction of transitional flow behavior in fluid systems. DOI: 10.1115/1.2979230

508 citations


Journal ArticleDOI
TL;DR: A quantitative estimate of the anisotropic power and scaling of magnetic field fluctuations in inertial range magnetohydrodynamic turbulence, using a novel wavelet technique applied to spacecraft measurements in the solar wind, shows for the first time that the spacecraft-frame spectrum has a spectral index near 2.
Abstract: We present a quantitative estimate of the anisotropic power and scaling of magnetic field fluctuations in inertial range magnetohydrodynamic turbulence, using a novel wavelet technique applied to spacecraft measurements in the solar wind. We show for the first time that, when the local magnetic field direction is parallel to the flow, the spacecraft-frame spectrum has a spectral index near 2. This can be interpreted as the signature of a population of fluctuations in field-parallel wavenumbers with a k −2 k spectrum but is also consistent with the presence of a “critical balance” style turbulent cascade. We also find, in common with previous studies, that most of the power is contained in wavevectors at large angles to the local magnetic field and that this component of the turbulence has a spectral index of 5/3. Magnetised plasmas fill most of the Universe and in many regions, turbulence plays important roles in the transport of energy and momentum and the acceleration and scattering of charged particles. Many aspects of plasma turbulence remain poorly understood, however. Here we present results on one of these, the anisotropy of the energy spectrum of magnetohydrodynamic (MHD) turbulence with respect to the magnetic field. In classical hydrodynamics, velocity fluctuations δuk with a wavenumber k decay and transfer energy to smaller scales on the shear timescale, τS ≈ 1/(kδuk). Within the steady inertial range, far from the energy input (“outer”) and dissipation scales, this leads to the dimensional result (δuk) 3 ∝ ǫ/k, where ǫ is the energy dissipation rate per unit mass. This gives the familiar Kolmogorov energy spectrum P(k) ∝ k −5/3 , widely observed in hydrodynamic turbulence. In a plasma, fluc⊥ ! .

410 citations


Journal ArticleDOI
TL;DR: In this article, a new version of the SST model (SST-CC) has been extensively tested on a wide range of both wall-bounded and free shear turbulent flows with system rotation and/or streamline curvature.
Abstract: A rotation-curvature correction suggested earlier by Spalart and Shur (1997, "On the Sensitization of Turbulence Models to Rotation and Curvature, " Aerosp Sci Technol, 1(5), pp 297-302) for the one-equation Spalart-Allmaras turbulence model is adapted to the shear stress transport model This new version of the model (SST-CC) has been extensively tested on a wide range of both wall-bounded and free shear turbulent flows with system rotation and/or streamline curvature Predictions of the SST-CC model are compared with available experimental and direct numerical simulations (DNS) data, on the one hand, and with the corresponding results of the original SST model and advanced Reynolds stress transport model (RSM), on the other hand It is found that in terms of accuracy the proposed model significantly improves the original SST model and is quite competitive with the RSM, whereas its computational cost is significantly less than that of the RSM

345 citations


Journal ArticleDOI
TL;DR: In this paper, a digital filter-based generation of turbulent inflow conditions exploiting this fact is presented as a suitable technique for large eddy simulations computation of spatially developing flows.
Abstract: Using a numerical weather forecasting code to provide the dynamic large-scale inlet boundary conditions for the computation of small-scale urban canopy flows requires a continuous specification of appropriate inlet turbulence. For such computations to be practical, a very efficient method of generating such turbulence is needed. Correlation functions of typical turbulent shear flows have forms not too dissimilar to decaying exponentials. A digital-filter-based generation of turbulent inflow conditions exploiting this fact is presented as a suitable technique for large eddy simulations computation of spatially developing flows. The artificially generated turbulent inflows satisfy the prescribed integral length scales and Reynolds-stress-tensor. The method is much more efficient than, for example, Klein’s (J Comp Phys 186:652–665, 2003) or Kempf et al.’s (Flow Turbulence Combust, 74:67–84, 2005) methods because at every time step only one set of two-dimensional (rather than three-dimensional) random data is filtered to generate a set of two-dimensional data with the appropriate spatial correlations. These data are correlated with the data from the previous time step by using an exponential function based on two weight factors. The method is validated by simulating plane channel flows with smooth walls and flows over arrays of staggered cubes (a generic urban-type flow). Mean velocities, the Reynolds-stress-tensor and spectra are all shown to be comparable with those obtained using classical inlet-outlet periodic boundary conditions. Confidence has been gained in using this method to couple weather scale flows and street scale computations.

302 citations


Journal ArticleDOI
TL;DR: In this article, an extra term was added to the approximate transport equation for the turbulence energy dissipation rate, which significantly improved agreement with experimental data to validate the model, results are presented for the Nibe wind turbine ‘B’, the Danwin 180/23, and the MOD-0A.

255 citations


Journal ArticleDOI
TL;DR: In this paper, a turbulence transport model is employed to compute distributions of turbulence throughout the heliosphere, including advection, expansion, and reflection effects as well as the tendency toward dynamic alignment, and a von Karman type dissipation function that represents decay of turbulence due to cascade to small scales.
Abstract: [1] We employ a turbulence transport model to compute distributions of turbulence throughout the heliosphere. The model determines the radial dependence of three (coupled) quantities that characterize interplanetary turbulence, the energy per unit mass, the cross helicity or Alfvenicity, and a similarity length scale. A fourth integrated quantity, the plasma temperature, is modified by heat deposition due to turbulent dissipation. The model includes advection, expansion, and reflection effects as well as the tendency toward dynamic alignment, and a von Karman type dissipation function that represents decay of turbulence due to cascade to small scales. Two types of forcing are also featured, one a simple model of stream shear, and the other a driving in the outer heliosphere associated with wave energy injection due to pickup protons of interstellar origin. Parameters for the model have been tuned using observation data from Voyager and Ulysses. We analyze the constraining observations to provide boundary conditions and parameters that vary with heliocentric latitude, with some extrapolations. The fully assembled model permits the computation of the distribution of turbulence throughout the entire heliosphere, and we present solutions for several appropriate parameter sets.

209 citations


Journal ArticleDOI
TL;DR: By analyzing trajectories of solid hydrogen tracers, it is found that the distributions of velocity in decaying quantum turbulence in superfluid 4He are strongly non-Gaussian with 1/v(3) power-law tails.
Abstract: By analyzing trajectories of solid hydrogen tracers, we find that the distributions of velocity in decaying quantum turbulence in superfluid 4 He are strongly non-Gaussian with 1=v 3 power-law tails. These features differ from the near-Gaussian statistics of homogenous and isotropic turbulence of classical fluids. We examine the dynamics of many events of reconnection between quantized vortices and show by simple scaling arguments that they produce the observed power-law tails.

192 citations


Journal Article
TL;DR: In this paper, a set of direct numerical simulations of isotropic turbulence passing through a nominally normal shock wave is presented, and the instantaneous structure of the shock/turbulence interaction is examined using averages conditioned on the instantaneous shock strength.
Abstract: A set of direct numerical simulations of isotropic turbulence passing through a nominally normal shock wave is presented. Upstream of the shock, the microscale Reynolds number is 40, the mean Mach number is 1.3–6.0, and the turbulence Mach number is 0.16–0.38. It is shown that the Kolmogorov scale decreases during the shock interaction, which implies that the grid resolution needed to resolve the viscous dissipation is finer than that used in previous studies. This leads to some qualitative differences with previous work, e.g., a rapid increase in the streamwise vorticity variance behind the shock and large anisotropy of the postshock Reynolds stresses. The instantaneous structure of the shock/turbulence interaction is examined using averages conditioned on the instantaneous shock strength. For locally strong compressions, the flow is characterized by overcompression, followed by an expansion. At points where the shock is locally weak, the profiles differ qualitatively depending on the strength of the inc...

188 citations


Journal ArticleDOI
TL;DR: In this article, the performance of low-Reynolds number turbulence models in predicting mixed convection heat transfer to fluids at supercritical pressure, especially paying attention to the features which enable them to respond to the modifications of the turbulence field due to influences of flow acceleration and buoyancy, was evaluated.

178 citations


Journal ArticleDOI
TL;DR: In this article, a simulation of flow over random urban-like obstacles is presented to gain a deeper insight into the effects of randomness in the obstacle topology, e.g. spatially-averaged mean velocity, Reynolds stresses, turbulence kinetic energy and dispersive stresses.
Abstract: Further to our previous large-eddy simulation (LES) of flow over a staggered array of uniform cubes, a simulation of flow over random urban-like obstacles is presented. To gain a deeper insight into the effects of randomness in the obstacle topology, the current results, e.g. spatially-averaged mean velocity, Reynolds stresses, turbulence kinetic energy and dispersive stresses, are compared with our previous LES data and direct numerical simulation data of flow over uniform cubes. Significantly different features in the turbulence statistics are observed within and immediately above the canopy, although there are some similarities in the spatially-averaged statistics. It is also found that the relatively high pressures on the tallest buildings generate contributions to the total surface drag that are far in excess of their proportionate frontal area within the array. Details of the turbulence characteristics (like the stress anisotropy) are compared with those in regular roughness arrays and attempts to find some generality in the turbulence statistics within the canopy region are discussed.

Journal ArticleDOI
TL;DR: In this article, a numerical investigation of turbulent forced convection in a two-dimensional channel with periodic transverse grooves on the lower channel wall is conducted, where the lower wall is subjected to a uniform heat flux condition while the upper wall is insulated.

Journal ArticleDOI
TL;DR: In this paper, the authors extend the energy analysis to turbulent potential and total energies (TPE and TTE = TKE + TPE), and conclude that TTE is a conservative parameter maintained by shear in any stratification.
Abstract: Traditionally, turbulence energetics is characterized by turbulent kinetic energy (TKE) and modelled using solely the TKE budget equation In stable stratification, TKE is generated by the velocity shear and expended through viscous dissipation and work against buoyancy forces The effect of stratification is characterized by the ratio of the buoyancy gradient to squared shear, called Richardson number, Ri It is widely believed that at Ri exceeding a critical value, Ric, local shear cannot maintain turbulence, and the flow becomes laminar We revise this concept by extending the energy analysis to turbulent potential and total energies (TPE and TTE = TKE + TPE), consider their budget equations, and conclude that TTE is a conservative parameter maintained by shear in any stratification Hence there is no "energetics Ric", in contrast to the hydrodynamic-instability threshold, Ric-instability, whose typical values vary from 025 to 1 We demonstrate that this interval, 025 >1, clarify principal difference between turbulent boundary layers and free flows, and provide basis for improving operational turbulence closure models

Journal ArticleDOI
TL;DR: In this article, a large-eddy-based methodology for the simulation of turbulent sprays is discussed, where the transport equations for the spatially filtered gas phase variables, in which source terms accounting for the droplet effects are added, are solved together with a probabilistic description of the liquid phase.
Abstract: A large-eddy-based methodology for the simulation of turbulent sprays is discussed. The transport equations for the spatially filtered gas phase variables, in which source terms accounting for the droplet effects are added, are solved together with a probabilistic description of the liquid phase. The probabilistic approach for the liquid phase is based on the transport equation for the spatially filtered joint probability density function of the variables required in order to describe the state of the liquid phase. In this equation, unclosed terms representing the filtered Lagrangian rates of change of the variables describing the spray are present. General modelling ideas for subgrid-scale (SGS) effects are proposed. The capabilities of the approach and the validity of the closure models, with particular with respect to the SGS dispersion, are investigated through application to a dilute particle-laden turbulent mixing layer. It is demonstrated that the formulation is able to reproduce very closely the measured properties of both the continuous and dispersed phases. The large-eddy simulation (LES) results are also found to be entirely consistent with the experimentally observed characteristics of droplet–gas turbulence interactions. Consistent with direct numerical simulation (DNS) studies of isotropic turbulence laden with particles where the entire turbulence spectrum is found to be modulated by the presence of particles, the present investigation, which comprises the effects of particle transport upon the large-scale vortical structures of a turbulent shear flow, highlights what appears to be a selective behaviour; few large-scale frequencies gain energy whereas the remaining modes are damped.

Journal ArticleDOI
TL;DR: Lifetime measurements of turbulence in pipe flow spanning 8 orders of magnitude in time are presented, showing that no critical point exists in this regime and that in contrast to the prevailing view the turbulent state remains transient.
Abstract: The collapse of turbulence, observable in shear flows at low Reynolds numbers, raises the question if turbulence is generically of a transient nature or becomes sustained at some critical point. Recent data have led to conflicting views with the majority of studies supporting the model of turbulence turning into an attracting state. Here we present lifetime measurements of turbulence in pipe flow spanning 8 orders of magnitude in time, drastically extending all previous investigations. We show that no critical point exists in this regime and that in contrast to the prevailing view the turbulent state remains transient. To our knowledge this is the first observation of superexponential transients in turbulence, confirming a conjecture derived from low-dimensional systems.

Journal ArticleDOI
TL;DR: In this article, it is hypothesized that the dynamics of stratified turbulence explain these data and that these dynamics may also explain the behavior of strongly stratified flows in similar dynamic ranges of other geophysical flows.
Abstract: Several existing sets of smaller-scale ocean and atmospheric data appear to display Kolmogorov-Obukov-Corrsin inertial ranges in horizontal spectra for length scales up to at least a few hundred meters. It is argued here that these data are inconsistent with the assumptions for these inertial range theories. Instead, it is hypothesized that the dynamics of stratified turbulence explain these data. If valid, these dynamics may also explain the behavior of strongly stratified flows in similar dynamic ranges of other geophysical flows.

Journal ArticleDOI
TL;DR: In this paper, the axial growth of columnar structures is measured using two-point correlations and in all cases the results are consistent with structure formation via linear inertial wave propagation.
Abstract: One of the most striking features of rotating turbulence is the inevitable appearance of large-scale columnar structures. Whilst these structures are frequently observed, the processes by which they are created are still poorly understood. In this paper we consider the emergence of these structures from freely decaying, rotating turbulence with Ro ∼ 1. Our study follows the conjecture by Davidson, Staplehurst & Dalziel ( J. Fluid Mech. , vol. 557, 2006, p. 135) that the structure formation may be due to linear inertial wave propagation, which was shown to be consistent with the growth of columnar eddies in inhomogeneous turbulence. Here we extend that work and consider the case of homogeneous turbulence. We describe laboratory experiments where homogeneous turbulence is created in a rotating tank. The turbulence is generated with Ro ∼ 1, and as the energy decays, the formation of columnar vortices is observed. The axial growth of these columnar structures is then measured using two-point correlations and in all cases the results are consistent with structure formation via linear inertial wave propagation. In particular, we obtain a self-similar collapse of the two-point correlations when the axial coordinate is normalized by Ω tb , where b is a measure of the integral scale in the horizontal plane and Ω is the rotation rate. Although our results do not exclude the possibility of significant nonlinear dynamics, they are consistent with the conjecture of Davidson et al . (2006) that linear dynamics play a strong guiding hand in structure formation.

Journal ArticleDOI
TL;DR: In this paper, two-dimensional electromagnetic particle-in-cell simulations in a magnetized, homogeneous, collisionless electron-proton plasma demonstrate the forward cascade of whistler turbulence.
Abstract: Two-dimensional electromagnetic particle-in-cell simulations in a magnetized, homogeneous, collisionless electron-proton plasma demonstrate the forward cascade of whistler turbulence. The simulations represent decaying turbulence, in which an initial, narrowband spectrum of fluctuations at wavenumbers kc∕ωe≃0.1 cascades toward increased damping at kc∕ωe≃1.0, where c∕ωe is the electron inertial length. The turbulence displays magnetic energy spectra that are relatively steep functions of wavenumber and are anisotropic with more energy in directions relatively perpendicular to the background magnetic field Bo=xBo than at the same wavenumbers parallel to Bo. In the weak turbulence regime, the primary new results of the simulations are as follows: (1) Magnetic spectra of the cascading fluctuations become more anisotropic with increasing fluctuation energy; (2) the wavevector dependence of the three magnetic energy ratios, ∣δBj∣2∕∣δB∣2 with j=x,y,z, show good agreement with linear dispersion theory for whistl...

Journal ArticleDOI
TL;DR: In this article, state-of-the-art observational and numerical modeling methods for small scale turbulence and mixing with applications to coastal oceans are presented in one context, based on the approach that modern process-oriented studies should be based on both observations and models.

Journal ArticleDOI
TL;DR: In particular, the transition depends sensitively on initial conditions and the turbulent state is not persistent but has an exponential distribution of lifetimes as discussed by the authors, and coherent structures are embedded within the turbulent dynamics, which transiently show up in the temporal evolution of turbulent flow.
Abstract: Plane Couette flow and pressure-driven pipe flow are two examples of flows where turbulence sets in while the laminar profile is still linearly stable. Experiments and numerical studies have shown that the transition has features compatible with the formation of a strange saddle rather than an attractor. In particular, the transition depends sensitively on initial conditions and the turbulent state is not persistent but has an exponential distribution of lifetimes. Embedded within the turbulent dynamics are coherent structures, which transiently show up in the temporal evolution of the turbulent flow. Here we summarize the evidence for this transition scenario in these two flows, with an emphasis on lifetime studies in the case of plane Couette flow and on the coherent structures in pipe flow.

Journal ArticleDOI
TL;DR: In this paper, a detailed analysis of multiscale structures and the identification of long-lived streamer-like wavemodes in a magnetically confined plasma provides new insight into the physics of plasma turbulence.
Abstract: Detailed analysis of multiscale structures and the identification of long-lived streamer-like wavemodes in a magnetically confined plasma provides new insight into the physics of plasma turbulence.

Journal ArticleDOI
TL;DR: A reformulation of the refined similarity hypothesis in terms of the mass-weighted velocity rho1/3v yields scaling laws that are almost insensitive to the forcing of supersonic turbulent flow, implying that the most intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence.
Abstract: The statistical properties of turbulence are considered to be universal at sufficiently small length scales, i.e., independent of boundary conditions and large-scale forces acting on the fluid. Analyzing data from numerical simulations of supersonic turbulent flow driven by external forcing, we demonstrate that this is not generally true for the two-point velocity statistics of compressible turbulence. However, a reformulation of the refined similarity hypothesis in terms of the mass-weighted velocity rho1/3v yields scaling laws that are almost insensitive to the forcing. The results imply that the most intermittent dissipative structures are shocks closely following the scaling of Burgers turbulence.

Journal ArticleDOI
TL;DR: In this article, the authors study the probability distribution function (PDF) of the mass density in simulations of supersonic turbulence with properties appropriate for molecular clouds and find that the PDF is not a good measure of magnetic field strength.
Abstract: We study the probability distribution function (PDF) of the mass density in simulations of supersonic turbulence with properties appropriate for molecular clouds. For this study we use Athena, a new higher-order Godunov code. We find there are surprisingly similar relationships between the mean of the time-averaged PDF and the turbulent Mach number for driven hydrodynamic and strong-field MHD turbulence. There is, however, a large scatter about these relations, indicating a high level of temporal and spatial variability in the PDF. Thus, the PDF of the mass density is unlikely to be a good measure of magnetic field strength. We also find that the PDF of decaying MHD turbulence deviates from the mean-Mach relation found in the driven case. This implies that the instantaneous Mach number alone is not enough to determine the statistical properties of turbulence that is out of equilibrium. The scatter about the mean-Mach relation for driven turbulence, along with the large departure of decaying turbulence PDFs from those of driven turbulence, may illuminate one factor contributing to the large observed cloud-to-cloud variation in the star formation rate per solar mass.

Journal ArticleDOI
TL;DR: In this paper, the authors derived the alpha effect and turbulent magnetic diffusivity for mean magnetic fields with profiles of different length scales from simulations of isotropic turbulence, and then related these results to nonlocal formulations in which alpha and the turbulent magnetics correspond to integral kernels.
Abstract: Aims. We determine the alpha effect and turbulent magnetic diffusivity for mean magnetic fields with profiles of different length scales from simulations of isotropic turbulence. We then relate these results to nonlocal formulations in which alpha and the turbulent magnetic diffusivity correspond to integral kernels. Methods. We solve evolution equations for magnetic fields that give the response to imposed test fields. These test fields correspond to mean fields with various wavenumbers. Both an imposed fully helical steady flow consisting of a pattern of screw-like motions (Roberts flow) and time-dependent, statistically steady isotropic turbulence are considered. In the latter case the evolution equations are solved simultaneously with the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate alpha and magnetic diffusivity tensors. Results. For both, the Roberts flow under the second-order correlation approximation and the isotropic turbulence alpha and turbulent magnetic diffusivity are greatest on large scales and these values diminish toward smaller scales. In both cases, the alpha effect and turbulent diffusion kernels are approximated by exponentials, corresponding to Lorentzian profiles in Fourier space. For isotropic turbulence, the turbulent diffusion kernel is half as wide as the alpha effect kernel. For the Roberts flow beyond the second-order correlation approximation, the turbulent diffusion kernel becomes negative on large scales.

Journal ArticleDOI
TL;DR: In this article, it was shown that isotropic helical turbulence leads to an α-effect and a turbulent diffusivity whose values are independent of the magnetic Reynolds number, provided Rm exceeds unity.
Abstract: Using numerical simulations at moderate magnetic Reynolds numbers up to 220, it is shown that in the kinematic regime, isotropic helical turbulence leads to an α-effect and a turbulent diffusivity whose values are independent of the magnetic Reynolds number, Rm, provided Rm exceeds unity. These turbulent coefficients are also consistent with expectations from the first-order smoothing approximation. For small values of Rm, α and turbulent diffusivity are proportional to Rm. Over finite time-intervals, meaningful values of α and turbulent diffusivity can be obtained even when there is small-scale dynamo action that produces strong magnetic fluctuations. This suggests that the fields generated by the small-scale dynamo do not make a correlated contribution to the mean electromotive force.

Journal Article
TL;DR: Comparison of CFD results with experimental data shows that the Reynolds Stress turbulence model yields a reasonably good prediction, and the use of the Presto interpolation scheme for pressure, the Simplec algorithm for pressure-velocity coupling and the quadratic upstream interpolation for convective kinetics scheme for momentum variables gives satisfactory results for highly swirling flows in cyclones.
Abstract: The aim of this study is to investigate the suitability of various numerical schemes and turbulence models in highly complex swirling flows which occur in tangential inlet cyclones. Three-dimensional steady governing equations for incompressible turbulent flow inside a cyclone were solved numerically using Fluent CFD (computational fluid dynamics) code. The Reynolds stress turbulence model, the Standard k-e and the RNG k-e turbulence models together with various combinations of numerical schemes are used to obtain axial and tangential velocity profiles, pressure drop and turbulent quantities. Computational results were compared with experimental and numerical values given in the literature, so as to evaluate the performance of the numerical schemes and turbulent models. Comparison of CFD results with experimental data shows that the Reynolds Stress turbulence model yields a reasonably good prediction. Results obtained from the numerical tests have demonstrated that the use of the Presto interpolation scheme for pressure, the Simplec algorithm for pressure-velocity coupling and the quadratic upstream interpolation for convective kinetics (quick) scheme for momentum variables gives satisfactory results for highly swirling flows in cyclones.

Journal ArticleDOI
TL;DR: A novel dimensionless parameter, the particle moment number Pa, was derived using dimensional analysis of the particle-laden Navier-Stokes equations, in order to understand the underlying physics of turbulence modification by particles.
Abstract: A novel dimensionless parameter, the particle moment number Pa, was derived using dimensional analysis of the particle-laden Navier-Stokes equations, in order to understand the underlying physics of turbulence modification by particles. A set of 80 previous experimental measurements where the turbulent kinetic energy was modified by particles was examined, and all results could clearly be divided into three groups in Re-Pa mappings. The turbulence attenuation region was observed between the augmentation regions with two critical particle momentum numbers.

Journal ArticleDOI
TL;DR: In this paper, the authors present results from two field studies of the nocturnal stable atmospheric boundary layer (SBL) over the Great Plains of the United States, where data from a scanning remote-sensing system, a High-Resolution Doppler Lidar (HRDL), provided measurements of mean and turbulent wind components at high spatial and temporal resolution through the lowest 500-1000 m of the atmosphere.
Abstract: This paper reviews results from two field studies of the nocturnal stable atmospheric boundary layer (SBL) over the Great Plains of the United States. Data from a scanning remote-sensing system, a High-Resolution Doppler Lidar (HRDL), provided measurements of mean and turbulent wind components at high spatial and temporal resolution through the lowest 500–1000 m of the atmosphere. This data set has allowed the characteristics of the low-level jet (LLJ) maximum (speed, height, direction) to be documented through entire nights. LLJs form after sunset and produce strong shear in the layer below the LLJ maximum or nose, which is a source of turbulence and mixing in the SBL. Simultaneous HRDL measurements of turbulence quantities related to turbulence kinetic energy (TKE) has allowed the turbulence in the subjet layer to be related to LLJ properties. Turbulence structure was found to be a function of the bulk stability of the subjet layer. For the strong-LLJ (> 15 m s−1), weakly stable cases the strength of the turbulence is proportional to the strength of the LLJ. For these cases with nearly continuous turbulence in the subjet layer, low-level jet scaling, in which lengths are scaled by the LLJ height and velocity variables are scaled by the LLJ speed, was found to be appropriate. For the weak-wind (< 5 m s−1 in the lowest 200 m), very stable boundary layer (vSBL), the boundary layer was found to be very shallow (sometimes < 10 m deep), and turbulent fluxes between the earth’s surface and the atmosphere were found to be essentially shut down. For more intermediate wind speeds and stabilities, the SBL shows varying degrees of intermittency due to various mechanisms, including shearinstability and other gravity waves, density currents, and other mesoscale disturbances.

Journal ArticleDOI
TL;DR: In this paper, the authors compared the performance of a variety of turbulence models in simulating buoyancy-aided, turbulent mixed convection in vertical pipes, and found that the response to buoyancy of commonly used controlling parameters in turbulence damping functions varies significantly and that the performance can largely be correlated with the type of controlling parameter used.

Journal ArticleDOI
TL;DR: In this paper, a planar array of synthetic jets, firing upward in a spatiotemporally random pattern to create turbulence at an air-water interface, is studied.
Abstract: We report measurements of the flow above a planar array of synthetic jets, firing upwards in a spatiotemporally random pattern to create turbulence at an air–water interface. The flow generated by this randomly actuated synthetic jet array (RASJA) is turbulent, with a large Reynolds number and a weak secondary (mean) flow. The turbulence is homogeneous over a large region and has similar isotropy characteristics to those of grid turbulence. These properties make the RASJA an ideal facility for studying the behaviour of turbulence at boundaries, which we do by measuring one-point statistics approaching the air–water interface (via particle image velocimetry). We explore the effects of different spatiotemporally random driving patterns, highlighting design conditions relevant to all randomly forced facilities. We find that the number of jets firing at a given instant, and the distribution of the duration for which each jet fires, greatly affect the resulting flow. We identify and study the driving pattern that is optimal given our tank geometry. In this optimal configuration, the flow is statistically highly repeatable and rapidly reaches steady state. With increasing distance from the jets, there is a jet merging region followed by a planar homogeneous region with a power-law decay of turbulent kinetic energy. In this homogeneous region, we find a Reynolds number of 314 based on the Taylor microscale. We measure all components of mean flow velocity to be less than 10% of the turbulent velocity fluctuation magnitude. The tank width includes roughly 10 integral length scales, and because wall effects persist for one to two integral length scales, there is sizable core region in which turbulent flow is unaffected by the walls. We determine the dissipation rate of turbulent kinetic energy via three methods, the most robust using the velocity structure function. Having a precise value of dissipation and low mean flow allows us to measure the empirical constant in an existing model of the Eulerian velocity power spectrum. This model provides a method for determining the dissipation rate from velocity time series recorded at a single point, even when Taylor's frozen turbulence hypothesis does not hold. Because the jet array offers a high degree of flow control, we can quantify the effects of the mean flow in stirred tanks by intentionally forcing a mean flow and varying its strength. We demonstrate this technique with measurements of gas transfer across the free surface, and find a threshold below which mean flow no longer contributes significantly to the gas transfer velocity.