Experimental study of premixed flames in intense isotropic turbulence
Bedat, B.; Cheng, R.K.
1994-04-01
A methodology for investigating premixed turbulent flames propagating in intense isotropic turbulence has been developed. The burner uses a turbulence generator developed by Videto and Santavicca and the flame is stabilized by weak-swirl generated by air injectors. This set-up produces stable premixed turbulent flames under a wide range of mixture conditions and turbulence intensities. The experiments are designed to investigate systematically the changes in flame structures for conditions which can be classified as wrinkled laminar flames, corrugated flames and flames with distributed reaction zones. Laser Doppler anemometry and Rayleigh scattering techniques are used to determine the turbulence and scalar statistics. In the intense turbulence, the flames are found to produce very little changes in the mean and rams velocities. Their flame speed increase linearly with turbulence intensity as for wrinkled laminar flames. The Rayleigh scattering pdfs for flames within the distributed reaction zone regime are distinctly bimodal. The probabilities of the reacting states (i.e. contributions from within the reaction zone) is not higher than those of wrinkled laminar flame. These results show that there is no drastic changes in flame structures at Karlovitz number close to unity. This suggest that the Klimov-Williams criterion under-predicts the resilience of wrinkled flamelets to intense turbulence.
The structure of intense vorticity in homogeneous isotropic turbulence
NASA Technical Reports Server (NTRS)
Jimenez, J.; Wray, A. A.; Saffman, P. G.; Rogallo, R. S.
1992-01-01
The structure of the intense vorticity regions is studied in numerically simulated homogeneous, isotropic, equilibrium turbulent flow fields at four different Reynolds numbers in the range Re(sub lambda) = 36-171. In accordance with previous investigators, this vorticity is found to be organized in coherent, cylindrical or ribbon-like, vortices ('worms'). A statistical study suggests that they are just especially intense features of the background, O(omega'), vorticity. Their radii scale with the Kolmogorov microscale and their lengths with the integral scale of the flow. An interesting observation is that the Reynolds number based on the circulation of the intense vortices, gamma/nu, increases monotonically with Re(sub lambda), raising the question of the stability of the structures in the limit of Re(sub lambda) approaching infinity. One and two-dimensional statistics of vorticity and strain are presented; they are non-gaussian, and the behavior of their tails depends strongly on the Reynolds number. There is no evidence of convergence to a limiting distribution in our range of Re(sub lambda), even though the energy spectra and the energy dissipation rate show good asymptotic properties in the higher Reynolds number cases. Evidence is presented to show that worms are natural features of the flow and that they do not depend on the particular forcing scheme.
Decay Power Law in, High Intensity, Isotropic Turbulent Flow
NASA Astrophysics Data System (ADS)
Koster, Timothy; Puga, Alejandro; Larue, John
2014-11-01
In the study reported here, isotropy is determined using the measure proposed by George (1992), where isotropy corresponds to those downstream positions where the product of the Taylor Reynolds number and the skewness of the velocity derivative is a constant. Straight forward approach can be used which is based on the observation of Batchelor (1953), that the square of the Talor micorscale is linearly related to downstream distance relative to the virtual origin. The fact that the decay of downstream velocity variance is described by a power law is shown to imply power law behavior for various other parameters such as the dissipation, the integral length scale, the Taylor microscale, the Kolmogorov microscale and the Taylor Reynolds number and that there is an algebraic relationship between the various power law exponents. Results are presented for mean velocities of 6 and 8 m/s for the downstream decay of the parameters listed in the preceding. The corresponding values of the Taylor Reynolds number at the start of the isotropic region are 290 and 400, and the variance decay exponent and virtual origin are found to be respectively -1.707 and -1.298 and -27.95 and -5.757. The exponents in the decay law for the other parameters are found to be within +/- 3% of the expected values. University of California Irvine Research Funds.
Linearly Forced Isotropic Turbulence
NASA Technical Reports Server (NTRS)
Lundgren, T. S.
2003-01-01
Stationary isotropic turbulence is often studied numerically by adding a forcing term to the Navier-Stokes equation. This is usually done for the purpose of achieving higher Reynolds number and longer statistics than is possible for isotropic decaying turbulence. It is generally accepted that forcing the Navier-Stokes equation at low wave number does not influence the small scale statistics of the flow provided that there is wide separation between the largest and smallest scales. It will be shown, however, that the spectral width of the forcing has a noticeable effect on inertial range statistics. A case will be made here for using a broader form of forcing in order to compare computed isotropic stationary turbulence with (decaying) grid turbulence. It is shown that using a forcing function which is directly proportional to the velocity has physical meaning and gives results which are closer to both homogeneous and non-homogeneous turbulence. Section 1 presents a four part series of motivations for linear forcing. Section 2 puts linear forcing to a numerical test with a pseudospectral computation.
DNS of Shock / Isotropic Turbulence Interaction
NASA Astrophysics Data System (ADS)
Grube, Nathan; Taylor, Ellen; Martín, Pino
2010-11-01
We discuss DNS of Shock / Isotropic Turbulence Interactions (SITI). We vary the incoming turbulence Mach number up to 0.8 and the convective Mach number up to 5 in order to determine their effects on the interaction. These cases are challenging due to the presence of shocklets in the incoming turbulence as well as significant motion of the main shock. Shock-capturing must be used at all points while still maintaining low enough numerical dissipation to preserve the turbulent fluctuations. We use the linearly- and nonlinearly-optimized Weighted Essentially Non-Oscillatory (WENO) method[1,2]. Particular attention is paid to the inflow boundary condition, where we find the use of snapshots of "frozen" turbulence from decaying isotropic box simulations to be unsatisfactory. We instead use time-varying inflow data generated by a separate forced isotropic turbulence simulation with a specified convection speed. This allows us to access flow conditions where the assumptions of Taylor's Hypothesis are not met. 1.) Mart'in, M.P., Taylor, E.M., Wu, M., and Weirs, V.G., JCP 220(1) 270-89, 2006. 2.) Taylor, E.M., Wu, M., and Mart'in, M.P., JCP 223(1) 384-97, 2007.
Some Basic Laws of Isotropic Turbulent Flow
NASA Technical Reports Server (NTRS)
Loitsianskii, L. G.
1945-01-01
An Investigation is made of the diffusion of artificially produced turbulence behind screens or other turbulence producers. The method is based on the author's concept of disturbance moment as a certain theoretically well-founded measure of turbulent disturbances.
The signature of initial production mechanisms in isotropic turbulence decay
NASA Astrophysics Data System (ADS)
Meldi, M.
2016-03-01
In the present work the quantification of the time-lasting effects of production mechanisms in homogeneous isotropic turbulence decay is addressed. The analysis is developed through the use of theoretical tools as well as numerical calculations based on the eddy damped quasinormal Markovian (EDQNM) model. In both cases a modified Lin equation is used, which accounts for production mechanisms as proposed by Meldi, Lejemble, and Sagaut ["On the emergence of non-classical decay regimes in multiscale/fractal generated isotropic turbulence," J. Fluid Mech. 756, 816-843 (2014)]. The approaches used show that an exponential decay law can be observed if the intensity of the forcing is strong enough to drive the turbulence dynamics, before a power-law decay is eventually attained. The EDQNM numerical results indicate that the exponential regime can persist for long evolution times, longer than the observation time in grid turbulence experiments. A rigorous investigation of the self-similar behavior of the pressure spectrum has been performed by a comprehensive comparison of EDQNM data with direct numerical simulation (DNS)/experiments in the literature. While DNS and free decay EDQNM simulations suggest the need of a very high Reλ threshold in order to observe a clear -7/3 slope of the pressure inertial range, experimental data and forced EDQNM calculations indicate a significantly lower value. This observation suggests that the time-lasting effects of production mechanisms, which cannot be excluded in experiments, play a role in the lack of general agreement with classical numerical approaches. These results reinforce the urge to evolve the numerical simulation state of the art towards the prediction of realistic physical states.
Sudden relaminarisation and lifetimes in forced isotropic turbulence
NASA Astrophysics Data System (ADS)
Linkmann, Moritz; Morozov, Alexander
2015-11-01
We demonstrate an unexpected connection between isotropic turbulence and wall-bounded shear flows. We perform direct numerical simulations of isotropic turbulence forced at large scales at moderate Reynolds numbers and observe sudden transitions from chaotic dynamics to a spatially simple flow, analogous to the laminar state in wall bounded shear flows. We find that the survival probabilities of turbulence are exponential and the typical lifetimes increase super-exponentially with the Reynolds number, similar to results on relaminarisation of localised turbulence in pipe and plane Couette flow. Results from simulations subjecting the observed large-scale flow to random perturbations of variable amplitude demonstrate that it is a linearly stable simple exact solution that can be destabilised by a finite-amplitude perturbation, like the Hagen-Poiseuille profile in pipe flow. Our results suggest that both isotropic turbulence and wall-bounded shear flows qualitatively share the same phase-space dynamics.
The energy decay in self-preserving isotropic turbulence revisited
NASA Technical Reports Server (NTRS)
Speziale, Charles G.; Bernard, Peter S.
1991-01-01
The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.
The energy decay in self-preserving isotropic turbulence revisited
NASA Technical Reports Server (NTRS)
Speziale, Charles G.; Bernard, Peter S.
1992-01-01
The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup -alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.
Diffusion of Heat from a Line Source in Isotropic Turbulence
NASA Technical Reports Server (NTRS)
Uberoi, Mahinder S; Corrsin, Stanley
1953-01-01
An experimental and analytical study has been made of some features of the turbulent heat diffusion behind a line heated wire stretched perpendicular to a flowing isotropic turbulence. The mean temperature distributions have been measured with systematic variations in wind speed, size of turbulence-producing grid, and downstream location of heat source. The nature of the temperature fluctuation field has been studied. A comparison of Lagrangian and Eulerian analyses for diffusion in a nondecaying turbulence yields an expression for turbulent-heat-transfer coefficient in terms of turbulence velocity and a Lagrangian "scale." the ratio of Eulerian to Lagrangian microscale has been determined theoretically by generalization of a result of Heisenberg and with arbitrary constants taken from independent sources, shows rough agreement with experimental results. A convenient form has been deduced for the criterion of interchangeability of instantaneous space and time derivatives in a flowing turbulence.
Studies of Shock Wave Interactions with Homogeneous and Isotropic Turbulence
NASA Technical Reports Server (NTRS)
Briassulis, G.; Agui, J.; Watkins, C. B.; Andreopoulos, Y.
1998-01-01
A nearly homogeneous nearly isotropic compressible turbulent flow interacting with a normal shock wave has been studied experimentally in a large shock tube facility. Spatial resolution of the order of 8 Kolmogorov viscous length scales was achieved in the measurements of turbulence. A variety of turbulence generating grids provide a wide range of turbulence scales. Integral length scales were found to substantially decrease through the interaction with the shock wave in all investigated cases with flow Mach numbers ranging from 0.3 to 0.7 and shock Mach numbers from 1.2 to 1.6. The outcome of the interaction depends strongly on the state of compressibility of the incoming turbulence. The length scales in the lateral direction are amplified at small Mach numbers and attenuated at large Mach numbers. Even at large Mach numbers amplification of lateral length scales has been observed in the case of fine grids. In addition to the interaction with the shock the present work has documented substantial compressibility effects in the incoming homogeneous and isotropic turbulent flow. The decay of Mach number fluctuations was found to follow a power law similar to that describing the decay of incompressible isotropic turbulence. It was found that the decay coefficient and the decay exponent decrease with increasing Mach number while the virtual origin increases with increasing Mach number. A mechanism possibly responsible for these effects appears to be the inherently low growth rate of compressible shear layers emanating from the cylindrical rods of the grid.
Charge pariticle transport in the non-isotropic turbulences
NASA Astrophysics Data System (ADS)
Sun, P.; Jokipii, J. R.
2015-12-01
The scattering and diffusion of energetic charged particles is not only important for understanding phenomena such as diffusive shock acceleration but it also is a natural probe of the statistical characteristics of magnetohydrodynamic (MHD) turbulence. Although Parker's transport equation (Parker 1965) allows us to describe the propagation of charged particles, the transport coefficients needed in the equation must be determined. Using Quasi-Linear Theory (QLT, e.g. Jokipii (1966)), one finds that coefficients can be related to the correlation function or power spectrum of homogeneous magnetic turbulence. However, different turbulence models will generally have a different influence on particle's scattering and diffusion. Among those models developed in MHD Turbulence, such as isotropic, Slab plus 2D (Tu & Marsch 1993; Gray et al 1996; Bieber et al 1996), etc. Here, using test-particle orbit simulations to calculate the transport coefficients, we study particle transport in synthesized asymmetric turbulence using the form first proposed by Goldreich & Sridhar (1995). We developed and introduce a systematic method to synthesize scale-dependent non-isotropic magnetic turbulences. We also developed and introduce a method to synthesize the 3d turbulent magnetic field from the observed solar wind time series dataset. We present the comparison of their effects on charge particle transport with previous theories and models.
Laboratory Study of Homogeneous and Isotropic Turbulence at High Reynolds Number
NASA Astrophysics Data System (ADS)
Pecenak, Zachary; Dou, Zhongwang; Yang, Fan; Cao, Lujie; Liang, Zach; Meng, Hui
2013-11-01
To study particle dynamics modified by isotropic turbulence at high Reynolds numbers and provide experimental data for DNS validation, we have developed a soccer-ball-shaped truncated icosahedron turbulence chamber with 20 adjoining hexagon surfaces, 12 pentagon surfaces and twenty symettrically displaced fans, which form an enclosed chamber of 1m diameter. We use Particle Image Velocimetry (PIV) technique to characterize the base turbulent flow, using different PIV set ups to capture various characteristic scales of turbulence. Results show that the stationary isotropic turbulence field is a spherical domain with diameter of 40 mm with quasi-zero mean velocities. The maximum rms velocity is ~1.5 m/s, corresponding to a Taylor microscale Re of 450. We extract from the PIV velocity field the whole set of turbulent flow parameters including: turbulent kinetic energy, turbulent intensity, kinetic energy dissipation rate, large eddy length and time scales, the Kolmogorov length, time and velocity scales, Taylor microscale and Re, which are critical to the study of inter-particle statistics modified by turbulence. This research is funded by an NSF grant CBET-0967407.
Computation of the sound generated by isotropic turbulence
NASA Technical Reports Server (NTRS)
Sarkar, S.; Hussaini, M. Y.
1993-01-01
The acoustic radiation from isotropic turbulence is computed numerically. A hybrid direct numerical simulation approach which combines direct numerical simulation (DNS) of the turbulent flow with the Lighthill acoustic analogy is utilized. It is demonstrated that the hybrid DNS method is a feasible approach to the computation of sound generated by turbulent flows. The acoustic efficiency in the simulation of isotropic turbulence appears to be substantially less than that in subsonic jet experiments. The dominant frequency of the computed acoustic pressure is found to be somewhat larger than the dominant frequency of the energy-containing scales of motion. The acoustic power in the simulations is proportional to epsilon (M(sub t))(exp 5) where epsilon is the turbulent dissipation rate and M(sub t) is the turbulent Mach number. This is in agreement with the analytical result of Proudman (1952), but the constant of proportionality is smaller than the analytical result. Two different methods of computing the acoustic power from the DNS data bases yielded consistent results.
Computation of large-scale statistics in decaying isotropic turbulence
NASA Technical Reports Server (NTRS)
Chasnov, Jeffrey R.
1993-01-01
We have performed large-eddy simulations of decaying isotropic turbulence to test the prediction of self-similar decay of the energy spectrum and to compute the decay exponents of the kinetic energy. In general, good agreement between the simulation results and the assumption of self-similarity were obtained. However, the statistics of the simulations were insufficient to compute the value of gamma which corrects the decay exponent when the spectrum follows a k(exp 4) wave number behavior near k = 0. To obtain good statistics, it was found necessary to average over a large ensemble of turbulent flows.
A non-isotropic multiple-scale turbulence model
NASA Technical Reports Server (NTRS)
Chen, C. P.
1990-01-01
A newly developed non-isotropic multiple scale turbulence model (MS/ASM) is described for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts for modeling multiple scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model perform significantly better than the single scale k-epsilon model. The present model is relatively inexpensive in terms of CPU time which makes it suitable for broad engineering flow applications.
Bulk viscosity effect on freely decaying compressible homogeneous isotropic turbulence
NASA Astrophysics Data System (ADS)
Pan, Shaowu; Johnsen, Eric
2015-11-01
Despite growing interests in compressible turbulence, the effect of bulk viscosity has been long ignored. For certain gases, the bulk viscosity may be 1000 times greater than the shear viscosity and thus modify energy transfer and dissipation mechanisms. In this study, we use direct numerical simulations to investigate the role of bulk viscosity on decaying isotropic compressible turbulence. Our results show that bulk viscosity exhibits a negligible decrease on enstrophy, but moderate and significant increases on the turbulent kinetic energy and Taylor-scale Reynolds number, respectively. A Helmholtz decomposition of the velocity field indicates that the bulk viscosity has a negligible effect on the solenoidal part, but exhibits a cross-scale effect on the dilatational component.
Modification of homogeneous and isotropic turbulence by solid particles
NASA Astrophysics Data System (ADS)
Hwang, Wontae
2005-12-01
Particle-laden flows are prevalent in natural and industrial environments. Dilute loadings of small, heavy particles have been observed to attenuate the turbulence levels of the carrier-phase flow, up to 80% in some cases. We attempt to increase the physical understanding of this complex phenomenon by studying the interaction of solid particles with the most fundamental type of turbulence, which is homogeneous and isotropic with no mean flow. A flow facility was developed that could create air turbulence in a nearly-spherical chamber by means of synthetic jet actuators mounted on the corners. Loudspeakers were used as the actuators. Stationary turbulence and natural decaying turbulence were investigated using two-dimensional particle image velocimetry for the base flow qualification. Results indicated that the turbulence was fairly homogeneous throughout the measurement domain and very isotropic, with small mean flow. The particle-laden flow experiments were conducted in two different environments, the lab and in micro-gravity, to examine the effects of particle wakes and flow structure distortion caused by settling particles. The laboratory experiments showed that glass particles with diameters on the order of the turbulence Kolmogorov length scale attenuated the fluid turbulent kinetic energy (TKE) and dissipation rate with increasing particle mass loadings. The main source of fluid TKE production in the chamber was the speakers, but the loss of potential energy of the settling particles also resulted in a significant amount of production of extra TKE. The sink of TKE in the chamber was due to the ordinary fluid viscous dissipation and extra dissipation caused by particles. This extra dissipation could be divided into "unresolved" dissipation caused by local velocity disturbances in the vicinity of the small particles and dissipation caused by large-scale flow distortions from particle wakes and particle clusters. The micro-gravity experiments in NASA's KC-135
Hindered Energy Cascade in Highly Helical Isotropic Turbulence.
Stepanov, Rodion; Golbraikh, Ephim; Frick, Peter; Shestakov, Alexander
2015-12-01
The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade to small scales has never been observed. We suggest, for the first time, a generalized phenomenology for isotropic turbulence with an arbitrary spectral distribution of the helicity. We discuss various scenarios of direct turbulent cascades with new helicity effect, which can be interpreted as a hindering of the spectral energy transfer. Therefore, the energy is accumulated and redistributed so that the efficiency of nonlinear interactions will be sufficient to provide a constant energy flux. We confirm our phenomenology by high Reynolds number numerical simulations based on a shell model of helical turbulence. The energy in our model is injected at a certain large scale only, whereas the source of helicity is distributed over all scales. In particular, we found that the helical bottleneck effect can appear in the inertial interval of the energy spectrum. PMID:26684120
Hindered Energy Cascade in Highly Helical Isotropic Turbulence
NASA Astrophysics Data System (ADS)
Stepanov, Rodion; Golbraikh, Ephim; Frick, Peter; Shestakov, Alexander
2015-12-01
The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade to small scales has never been observed. We suggest, for the first time, a generalized phenomenology for isotropic turbulence with an arbitrary spectral distribution of the helicity. We discuss various scenarios of direct turbulent cascades with new helicity effect, which can be interpreted as a hindering of the spectral energy transfer. Therefore, the energy is accumulated and redistributed so that the efficiency of nonlinear interactions will be sufficient to provide a constant energy flux. We confirm our phenomenology by high Reynolds number numerical simulations based on a shell model of helical turbulence. The energy in our model is injected at a certain large scale only, whereas the source of helicity is distributed over all scales. In particular, we found that the helical bottleneck effect can appear in the inertial interval of the energy spectrum.
Interaction of a converging spherical shock wave with isotropic turbulence
NASA Astrophysics Data System (ADS)
Bhagatwala, Ankit; Lele, Sanjiva K.
2012-08-01
Simulations of converging spherical shock waves propagating through a region of compressible isotropic turbulence are carried out. Both converging and reflected phases of the shock are studied. Effect of the reflected phase of the shock is found to be quite different from the expanding shock in the Taylor blast wave-turbulence interaction problem. Vorticity and turbulent kinetic energy are amplified due to passage of the shock. Similar to the latter problem, the vorticity-dilatation term is primarily responsible for the observed behavior. This is confirmed through Eulerian and Lagrangian statistics. Transverse vorticity amplification is compared with linear planar shock-turbulence theory. The smallest eddies, represented by the Kolmogorov scale, decrease in size after passing through the converging shock and this is shown to be related to a decrease in kinematic viscosity and increase in dissipation behind the converging shock. Distortion of the shock due to turbulence is also investigated and quantified. Turbulence also affects maximum compression achieved at the point of shock reflection, when the shock radius is at a minimum. This decrease in compression is quantified by comparing with pure shock simulations.
Asymptotic behavior of curvature of surface elements in isotropic turbulence
NASA Technical Reports Server (NTRS)
Girimaji, S. S.
1991-01-01
The asymptotic behavior of the curvature of material elements in turbulence is investigated using Lagrangian velocity-gradient time series obtained from direct numerical simulations of isotropic turbulence. Several material-element ensembles of different initial curvatures and shapes are studied. It is found that, at long times, the (first five) moments of the logarithm of characteristic curvature and shape factor asymptote to values that are independent of the initial curvature or shape. This evidence strongly suggests that the asymptotic pdf's of the curvature and shape of material elements are stationary and independent of initial conditions. Irrespective of initial curvature or shape, the asymptotic shape of a material surface is cylindrical with a high probability.
Subfilter scalar-flux vector orientation in homogeneous isotropic turbulence.
Verma, Siddhartha; Blanquart, G
2014-06-01
The geometric orientation of the subfilter-scale scalar-flux vector is examined in homogeneous isotropic turbulence. Vector orientation is determined using the eigenframe of the resolved strain-rate tensor. The Schmidt number is kept sufficiently large so as to leave the velocity field, and hence the strain-rate tensor, unaltered by filtering in the viscous-convective subrange. Strong preferential alignment is observed for the case of Gaussian and box filters, whereas the sharp-spectral filter leads to close to a random orientation. The orientation angle obtained with the Gaussian and box filters is largely independent of the filter width and the Schmidt number. It is shown that the alignment direction observed numerically using these two filters is predicted very well by the tensor-diffusivity model. Moreover, preferred alignment of the scalar gradient vector in the eigenframe is shown to mitigate any probable issues of negative diffusivity in the tensor-diffusivity model. Consequentially, the model might not suffer from solution instability when used for large eddy simulations of scalar transport in homogeneous isotropic turbulence. Further a priori tests indicate poor alignment of the Smagorinsky and stretched vortex model predictions with the exact subfilter flux. Finally, strong filter dependence of subfilter scalar-flux orientation suggests that explicit filtering may be preferable to implicit filtering in large eddy simulations. PMID:25019887
Large Deviation Statistics of Vorticity Stretching in Isotropic Turbulence
NASA Astrophysics Data System (ADS)
Johnson, Perry; Meneveau, Charles
2015-11-01
A key feature of 3D fluid turbulence is the stretching/re-alignment of vorticity by the action of the strain-rate. It is shown using the cumulant-generating function that cumulative vorticity stretching along a Lagrangian path in isotropic turbulence behaves statistically like a sum of i.i.d. variables. The Cramer function for vorticity stretching is computed from the JHTDB isotropic DNS (Reλ = 430) and compared to those of the finite-time Lyapunov exponents (FTLE) for material deformation. As expected the mean cumulative vorticity stretching is slightly less than that of the most-stretched material line (largest FTLE), due to the vorticity's preferential alignment with the second-largest eigenvalue of strain-rate and the material line's preferential alignment with the largest eigenvalue. However, the vorticity stretching tends to be significantly larger than the second-largest FTLE, and the Cramer functions reveal that the statistics of vorticity stretching fluctuations are more similar to those of largest FTLE. A model Fokker-Planck equation is constructed by approximating the viscous destruction of vorticity with a deterministic non-linear relaxation law matching conditional statistics, while the fluctuations in vorticity stretching are modelled by stochastic noise matching the statistics encoded in the Cramer function. The model predicts a stretched-exponential tail for the vorticity magnitude PDF, with good agreement for the exponent but significant error (30-40%) in the pre-factor. Supported by NSF Graduate Fellowship (DGE-1232825) and NSF Grant CMMI-0941530.
Taylor length-scale size particles in Isotropic Turbulence
NASA Astrophysics Data System (ADS)
Lucci, Francesco
The present study investigates the two-way coupling effects of finite-size solid spherical particles on decaying isotropic turbulence using an immersed boundary method. The conventional point particle assumption is valid only in the case of particles with a diameter, dp, much smaller than the Kolmogorov length scale, eta. In a simulation with particles of diameter dp > eta the flow around each particle needs to be resolved. In this study, we use a method similar to that of Uhlmann(2005) [55] that adapts the Immersed Boundary(IB) Method developed by Peskin [38] to simulate the flow around suspended spherical solid particles. The main idea of the method is to distribute a number of Lagrangian points uniformly over the surface of the particle. A force is applied at each Lagrangian point to represent the momentum exchange between the particle and the surrounding fluid. An analytic three-point delta function is used to distribute the force to the Eulerian grid points saddling the spherical surface to satisfy the no-slip condition at each Lagrangian point. Decaying turbulence is simulated in a periodic box with a uniform mesh of up to (512)3 grid points and an initial microscale Reynolds number of up to Relambda0 = 110. We compare the single phase flow (SPF) with particle-laden flows with particles of different diameters. The density of the particle varies from 2.56 to 10 times that of the fluid. The effects of the particles on the temporal development of turbulence kinetic energy E(t), its dissipation rate epsilon( t), its two-way coupling rate of change Ψp( t) and frequency spectra E(o) are discussed. In this study, in contrast to particles with dp < eta [15], particles with dp > eta always increase the dissipation rate of turbulence kinetic energy, epsilon( t). In addition, Ψp(t) is always positive, whereas it can be positive or negative for particles with dp < eta. The balance between these two effects caused E(t) to be smaller than that of the single-phase flow
Streamlines in stationary homogeneous isotropic turbulence and fractal-generated turbulence
NASA Astrophysics Data System (ADS)
Boschung, J.; Peters, N.; Laizet, S.; Vassilicos, J. C.
2016-04-01
We compare streamline statistics in stationary homogeneous isotropic turbulence and in turbulence generated by a fractal square grid. We examine streamline segments characterised by the velocity difference {{Δ }}u and the distance l between extremum points. We find close agreement between the stationary homogeneous isotropic turbulence and the decay region of the fractal-generated turbulence as well as the production region of the fractal flow for small segments. The statistics of larger segments are very similar for the isotropic turbulence and the decay region, but differ for the production region. Specifically, we examine the first, second and third conditional mean < {[{{Δ }}u]}n| l> . Noticeably, non-vanishing < {[{{Δ }}u]}n| l> for n=1,3 are due to an asymmetry of positive and negative segments, i.e. those for which {{Δ }}u\\gt 0 and {{Δ }}u\\lt 0, respectively. This asymmetry is not only kinematic, but is also due to dissipative effects and therefore < {[{{Δ }}u]}n| l> contains cascade information.
Coherent Vortex Simulations of 3D isotropic turbulence
NASA Astrophysics Data System (ADS)
Goldstein, Daniel E.; Vasilyev, Oleg V.; Kevlahan, Nicholas K.-R.
2006-11-01
This is the first of three talks on the wavelet filter based dynamically adaptive eddy capturing computational methodology that unifies variable fidelity simulation approaches such as wavelet-based DNS, Coherent Vortex Simulation (CVS), and Stochastic Coherent Adaptive Large Eddy Simulation. The commonality of these approaches is their ability to identify and ``track" on an adaptive mesh energetic coherent vortical structures. In CVS the velocity field is decomposed into two orthogonal parts: a coherent, inhomogeneous, non-Gaussian component and an incoherent, homogeneous, Gaussian component. This separation of coherent and incoherent components is achieved by wavelet thresholding which can be viewed as a non-linear filter that depends on each flow realization. The essence of the CVS approach is to solve for the coherent non-Gaussian component of a turbulent flow field. It has been shown previously that second generation bi-orthogonal wavelet threshold filtering is able to decompose a turbulent velocity field such that the total resulting SGS dissipation is approximately zero. This physically allows a CVS simulation to recover low order statistics with no SGS model. In this work CVS simulations of decaying incompressible 3D isotropic turbulence are compared to DNS results. -6pt
Local topology of energy transport in isotropic turbulence
NASA Astrophysics Data System (ADS)
Boschung, Jonas; Meneveau, Charles
2012-11-01
Similar to the velocity vector field, whose tangent (stream) lines represent how fluid volume (or mass in constant density flows) is transported in the flow, it is of interest to consider the vector field corresponding to the transport of mechanical energy (Meyers & Meneveau, 2012). The transport includes advection and viscous diffusion. In order to characterize the local topology of this vector field in turbulence, we examine statistical properties of its gradient field. This energy transport field is not divergence-free, due to dissipation and unsteady changes of kinetic energy. Therefore, the first invariant (the trace) of its gradient tensor is not zero, as in compressible flow. The three invariants PE, QE and RE of the energy transport gradient tensor are analyzed using concepts developed earlier for analysis of compressible flows. Data from DNS of isotropic turbulence is used, from the JHU database (Li et al. 2008, JoT), as well as other sources. Contracting node-like topology occurs very frequently, consistent with the dissipative nature of fluid turbulence. Further topological properties are established based on conditional PDFs of the invariants, and flow visualizations are used to develop insights into the local structure of the energy transport vector field. This work is supported by project CMMI-0941530. The authors also thank Prof. J. Meyers, Prof. N. Peters and Mr. P. Schaefer for interesting discussions on this topic.
Energy transfer and constrained simulations in isotropic turbulence
NASA Technical Reports Server (NTRS)
Jimenez, Javier
1993-01-01
The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.
Dynamics of Aerosol Particles in Stationary, Isotropic Turbulence
NASA Technical Reports Server (NTRS)
Collins, Lance R.; Meng, Hui
2004-01-01
A detailed study of the dynamics of sub-Kolmogorov-size aerosol particles in stationary isotropic turbulence has been performed. The study combined direct numerical simulations (DNS; directed by Prof. Collins) and high-resolution experimental measurements (directed by Prof. Meng) under conditions of nearly perfect geometric and parametric overlap. The goal was to measure the accumulation of particles in low-vorticity regions of the flow that arises from the effect commonly referred to as preferential concentration. The grant technically was initiated on June 13, 2000; however, funding was not available until July 11, 2000. The grant was originally awarded to Penn State University (numerical simulations) and SUNY-Buffalo (experiments); however, Prof. Collins effort was moved to Cornell University on January 2002 when he joined that university. He completed the study there. A list of the specific tasks that were completed under this study is presented.
Fully developed isotropic turbulence: Symmetries and exact identities
NASA Astrophysics Data System (ADS)
Canet, Léonie; Delamotte, Bertrand; Wschebor, Nicolás
2015-05-01
We consider the regime of fully developed isotropic and homogeneous turbulence of the Navier-Stokes equation with a stochastic forcing. We present two gauge symmetries of the corresponding Navier-Stokes field theory and derive the associated general Ward identities. Furthermore, by introducing a local source bilinear in the velocity field, we show that these symmetries entail an infinite set of exact and local relations between correlation functions. They include in particular the Kármán-Howarth relation and another exact relation for a pressure-velocity correlation function recently derived in G. Falkovich, I. Fouxon, and Y. Oz [J. Fluid Mech. 644, 465 (2010)], 10.1017/S0022112009993429 that we further generalize.
Non-isotropic turbulence effects on spray combustion
NASA Technical Reports Server (NTRS)
Kim, Y. M.; Shang, H. M.; Chen, C. P.
1991-01-01
A numerical model for the prediction of local properties of statistically stationary spray-combusting flows is evaluated by comparison with experimental data. To appraise the relative performance of turbulence models, computations were carried out by the k-epsilon model and the algebraic stress model. The present numerical results show the qualitative agreement with experimental data. In terms of overall local flow properties, the algebraic stress model improves a degree of conformity to the experimental data due to its ability to introduce the nonisotropic turbulence effects. Two swirl numbers are considered to investigate the influence of swirl on the droplet evaporation and trajectories, and the effects of droplet/turbulence interactions in flow properties. It is found that the large swirl produces a higher evaporationn rate, and more intensive turbulent mixing and burning. The discrepancies observed in the results are attributed mainly to uncertainties in the initial spray size and velocity distributions, the droplet/wall impingement interaction, the combustion model with the fast chemistry and the turbulence models dealing with the strong streamline curvature and complex interactions between the dispersed droplets and the continuous gas-phase flows.
NASA Technical Reports Server (NTRS)
Eaton, John; Hwang, Wontae; Cabral, Patrick
2002-01-01
This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact
Clustering of vertically constrained passive particles in homogeneous isotropic turbulence
NASA Astrophysics Data System (ADS)
De Pietro, Massimo; van Hinsberg, Michel A. T.; Biferale, Luca; Clercx, Herman J. H.; Perlekar, Prasad; Toschi, Federico
2015-05-01
We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale η , that maximizes the clustering of the particles.
DNS of fully resolved spherical particles dispersed in isotropic turbulence
NASA Astrophysics Data System (ADS)
Lucci, Francesco; Ferrante, Antonino; Elghobashi, Said
2008-11-01
Our DNS study concerns the interactions between decaying isotropic turbulence and solid spherical particles with diameter, d, larger than the Kolmogorov length scale, η. We employ an Immersed Boundary method similar to that of Uhlmann (JCP, 2005) to resolve the flow around 6400 spherical particles with a volume fraction of φv=0.1. The monosize particles have a diameter, d = 16 ηo. Our simulations, with 256^3 mesh points and Reλ0= 75, cover a range of 38 <=τp/τKo<=149, for the ratio of the particle response time to the initial Kolmogorov time scale. A Lagrangian approach is used to compute the frequency spectrum of the turbulence kinetic energy (TKE) of the fluid phase. The effects of varying τp/τKo on the spectrum and the decay rate of TKE are discussed. The effects of the formation of the particle boundary layer on the viscous dissipation rate of TKE are also discussed.
Preferential concentration of heavy particles in compressible isotropic turbulence
NASA Astrophysics Data System (ADS)
Zhang, Qingqing; Liu, Han; Ma, Zongqiang; Xiao, Zuoli
2016-05-01
Numerical simulations of particle-laden compressible isotropic turbulence with Taylor Reynolds number Reλ ˜ 100 are conducted by using a high-order turbulence solver, which is based on high-order compact finite difference method in the whole flow domain and localized artificial diffusivities for discontinuities. For simplicity, only one-way coupling (i.e., the influence of fluid on particles) between the carrier flow and particles is considered. The focus is on the study of the preferential concentration of heavy particles in dissipative scale of turbulence and the underlying mechanisms. Firstly, the effect of Stokes number (St) on the particle distribution in flow of Mach 1.01 (referred to as high-Mach-number case in this study) is investigated as a necessary supplementation for the previous studies in incompressible and weakly compressible flows. It turns out that heavy particles with Stokes number close to unity exhibit the strongest preferential concentration, which is in agreement with the observation in incompressible flow. All types of heavy particles have a tendency to accumulate in high-density regions of the background flow. While all kinds of particles dominantly collect in low-vorticity regions, intermediate and large particles (St = 1 and St = 5) are also found to collect in high-vorticity regions behind the randomly formed shocklets. Secondly, the impact of turbulent Mach number (Mt) (or the compressibility) of the carrier flow on the spatial distribution of the particles with St = 1 is discussed using the simulated compressible flows with Mt being 0.22, 0.68, and 1.01, respectively. In low-Mach-number flow, particles tend to concentrate in regions of low vorticity due to the centrifuge effect of vortices and particle concentration decreases monotonically with the increasing vorticity magnitude. As Mach number increases, the degree of particle clustering is slightly weakened in low-vorticity regions but is enhanced in high-vorticity regions, which
The role of wall confinement on the decay rate of an initially isotropic turbulent field
NASA Astrophysics Data System (ADS)
Dowling, David R.; Movahed, Pooya; Johnsen, Eric
2014-11-01
The problem of freely decaying isotropic turbulence has been the subject of intensive research during the past few decades due to its importance for modeling purposes. While isotropy and periodic boundary conditions assumptions simplify the analysis, large-scale anisotropy (e.g., caused by rotation, shear, acceleration or walls) is in practice present in most turbulent flows and affects flow dynamics across different scales, as well as the kinetic energy decay. We investigate the role of wall confinement and viscous dissipation on the decay rate of an initially isotropic field for confining volumes of different aspect ratios. We first generate an isotropic velocity field in a cube with periodic boundary conditions. Next, using this field, we change the boundary conditions to no-slip walls on all sides. These walls restrict the initial field to a confined geometry and also provide an additional viscous dissipation mechanism. The problem is considered for confining volumes of different aspect ratios by adjusting the initial field. The change in confining volume introduces an additional length scale to the problem. Direct numerical simulation of the proposed set-up is used to verify the scaling arguments for the decay rate of kinetic energy. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation Grant Number ACI-1053575.
Symmetries and the approach to statistical equilibrium in isotropic turbulence
NASA Astrophysics Data System (ADS)
Clark, Timothy T.; Zemach, Charles
1998-11-01
The relaxation in time of an arbitrary isotropic turbulent state to a state of statistical equilibrium is identified as a transition to a state which is invariant under a symmetry group. We deduce the allowed self-similar forms and time-decay laws for equilibrium states by applying Lie-group methods (a) to a family of scaling symmetries, for the limit of high Reynolds number, as well as (b) to a unique scaling symmetry, for nonzero viscosity or nonzero hyperviscosity. This explains why a diverse collection of turbulence models, going back half a century, arrived at the same time-decay laws, either through derivations embedded in the mechanics of a particular model, or through numerical computation. Because the models treat the same dynamical variables having the same physical dimensions, they are subject to the same scaling invariances and hence to the same time-decay laws, independent of the eccentricities of their different formulations. We show in turn, by physical argument, by an explicitly solvable analytical model, and by numerical computation in more sophisticated models, that the physical mechanism which drives (this is distinct from the mathematical circumstance which allows) the relaxation to equilibrium is the cascade of turbulence energy toward higher wave numbers, with the rate of cascade approaching zero in the low wave-number limit and approaching infinity in the high wave-number limit. Only the low-wave-number properties of the initial state can influence the equilibrium state. This supplies the physical basis, beyond simple dimensional analysis, for quantitative estimates of relaxation times. These relaxation times are estimated to be as large as hundreds or more times the initial dominant-eddy cycle times, and are determined by the large-eddy cycle times. This mode of analysis, applied to a viscous turbulent system in a wind tunnel with typical initial laboratory parameters, shows that the time necessary to reach the final stage of decay is
Large-deviation statistics of vorticity stretching in isotropic turbulence
NASA Astrophysics Data System (ADS)
Johnson, Perry L.; Meneveau, Charles
2016-03-01
A key feature of three-dimensional fluid turbulence is the stretching and realignment of vorticity by the action of the strain rate. It is shown in this paper, using the cumulant-generating function, that the cumulative vorticity stretching along a Lagrangian path in isotropic turbulence obeys a large deviation principle. As a result, the relevant statistics can be described by the vorticity stretching Cramér function. This function is computed from a direct numerical simulation data set at a Taylor-scale Reynolds number of Reλ=433 and compared to those of the finite-time Lyapunov exponents (FTLE) for material deformation. As expected, the mean cumulative vorticity stretching is slightly less than that of the most-stretched material line (largest FTLE), due to the vorticity's preferential alignment with the second-largest eigenvalue of strain rate and the material line's preferential alignment with the largest eigenvalue. However, the vorticity stretching tends to be significantly larger than the second-largest FTLE, and the Cramér functions reveal that the statistics of vorticity stretching fluctuations are more similar to those of the largest FTLE. In an attempt to relate the vorticity stretching statistics to the vorticity magnitude probability density function in statistically stationary conditions, a model Kramers-Moyal equation is constructed using the statistics encoded in the Cramér function. The model predicts a stretched-exponential tail for the vorticity magnitude probability density function, with good agreement for the exponent but significant difference (35%) in the prefactor.
Carroll, Jonathan J.; Frank, Adam; Blackman, Eric G.
2010-10-10
Feedback from protostellar outflows can influence the nature of turbulence in star-forming regions even if they are not the primary source of velocity dispersion for all scales of molecular clouds. For the rate and power expected in star-forming regions, we previously (Carroll et al.) demonstrated that outflows could drive supersonic turbulence at levels consistent with the scaling relations from Matzner although with a steeper velocity power spectrum than expected for an isotropically driven supersonic turbulent cascade. Here, we perform higher resolution simulations and combine simulations of outflow driven turbulence with those of isotropically forced turbulence. We find that the presence of outflows within an ambient isotropically driven turbulent environment produces a knee in the velocity power spectrum at the outflow scale and a steeper slope at sub-outflow scales than for a purely isotropically forced case. We also find that the presence of outflows flattens the density spectrum at large scales effectively reducing the formation of large-scale turbulent density structures. These effects are qualitatively independent of resolution. We have also carried out Principal Component Analysis (PCA) for synthetic data from our simulations. We find that PCA as a tool for identifying the driving scale of turbulence has a misleading bias toward low amplitude large-scale velocity structures even when they are not necessarily the dominant energy containing scales. This bias is absent for isotropically forced turbulence but manifests strongly for collimated outflow driven turbulence.
Energy transfer and dissipation in forced isotropic turbulence
NASA Astrophysics Data System (ADS)
McComb, W. D.; Berera, A.; Yoffe, S. R.; Linkmann, M. F.
2015-04-01
A model for the Reynolds-number dependence of the dimensionless dissipation rate Cɛ was derived from the dimensionless Kármán-Howarth equation, resulting in Cɛ=Cɛ ,∞+C /RL+O (1 /RL2) , where RL is the integral scale Reynolds number. The coefficients C and Cɛ ,∞ arise from asymptotic expansions of the dimensionless second- and third-order structure functions. This theoretical work was supplemented by direct numerical simulations (DNSs) of forced isotropic turbulence for integral scale Reynolds numbers up to RL=5875 (Rλ=435 ), which were used to establish that the decay of dimensionless dissipation with increasing Reynolds number took the form of a power law RLn with exponent value n =-1.000 ±0.009 and that this decay of Cɛ was actually due to the increase in the Taylor surrogate U3/L . The model equation was fitted to data from the DNS, which resulted in the value C =18.9 ±1.3 and in an asymptotic value for Cɛ in the infinite Reynolds-number limit of Cɛ ,∞=0.468 ±0.006 .
Energy transfer and dissipation in forced isotropic turbulence.
McComb, W D; Berera, A; Yoffe, S R; Linkmann, M F
2015-04-01
A model for the Reynolds-number dependence of the dimensionless dissipation rate C(ɛ) was derived from the dimensionless Kármán-Howarth equation, resulting in C(ɛ)=C(ɛ,∞)+C/R(L)+O(1/R(L)(2)), where R(L) is the integral scale Reynolds number. The coefficients C and C(ɛ,∞) arise from asymptotic expansions of the dimensionless second- and third-order structure functions. This theoretical work was supplemented by direct numerical simulations (DNSs) of forced isotropic turbulence for integral scale Reynolds numbers up to R(L)=5875 (R(λ)=435), which were used to establish that the decay of dimensionless dissipation with increasing Reynolds number took the form of a power law R(L)(n) with exponent value n=-1.000±0.009 and that this decay of C(ɛ) was actually due to the increase in the Taylor surrogate U(3)/L. The model equation was fitted to data from the DNS, which resulted in the value C=18.9±1.3 and in an asymptotic value for C(ɛ) in the infinite Reynolds-number limit of C(ɛ,∞)=0.468±0.006. PMID:25974586
Interacting scales and energy transfer in isotropic turbulence
NASA Technical Reports Server (NTRS)
Zhou, YE
1993-01-01
The dependence of the energy transfer process on the disparity of the interacting scales is investigated in the inertial and far-dissipation ranges of isotropic turbulence. The strategy for generating the simulated flow fields and the choice of a disparity parameter to characterize the scaling of the interactions is discussed. The inertial range is found to be dominated by relatively local interactions, in agreement with the Kolmogorov assumption. The far-dissipation is found to be dominated by relatively non-local interactions, supporting the classical notion that the far-dissipation range is slaved to the Kolmogorov scales. The measured energy transfer is compared with the classical models of Heisenberg, Obukhov, and the more detailed analysis of Tennekes and Lumley. The energy transfer statistics measured in the numerically simulated flows are found to be nearly self-similar for wave numbers in the inertial range. Using the self-similar form measured within the limited scale range of the simulation, an 'ideal' energy transfer function and the corresponding energy flux rate for an inertial range of infinite extent are constructed. From this flux rate, the Kolmogorov constant is calculated to be 1.5, in excellent agreement with experiments.
Inverse energy cascade in three-dimensional isotropic turbulence.
Biferale, Luca; Musacchio, Stefano; Toschi, Federico
2012-04-20
We study the statistical properties of homogeneous and isotropic three-dimensional (3D) turbulent flows. By introducing a novel way to make numerical investigations of Navier-Stokes equations, we show that all 3D flows in nature possess a subset of nonlinear evolution leading to a reverse energy transfer: from small to large scales. Up to now, such an inverse cascade was only observed in flows under strong rotation and in quasi-two-dimensional geometries under strong confinement. We show here that energy flux is always reversed when mirror symmetry is broken, leading to a distribution of helicity in the system with a well-defined sign at all wave numbers. Our findings broaden the range of flows where the inverse energy cascade may be detected and rationalize the role played by helicity in the energy transfer process, showing that both 2D and 3D properties naturally coexist in all flows in nature. The unconventional numerical methodology here proposed, based on a Galerkin decimation of helical Fourier modes, paves the road for future studies on the influence of helicity on small-scale intermittency and the nature of the nonlinear interaction in magnetohydrodynamics. PMID:22680722
Statistics of pressure and pressure gradient in homogeneous isotropic turbulence
NASA Technical Reports Server (NTRS)
Gotoh, T.; Rogallo, R. S.
1994-01-01
The statistics of pressure and pressure gradient in stationary isotropic turbulence are measured within direct numerical simulations at low to moderate Reynolds numbers. It is found that the one-point pdf of the pressure is highly skewed and that the pdf of the pressure gradient is of stretched exponential form. The power spectrum of the pressure P(k) is found to be larger than the corresponding spectrum P(sub G)(k) computed from a Gaussian velocity field having the same energy spectrum as that of the DNS field. The ratio P(k)/P(sub G)(k), a measure of the pressure-field intermittence, grows with wavenumber and Reynolds number as -R(sub lambda)(exp 1/2)log(k/k(sub d)) for k less than k(sub d)/2 where k(sub d) is the Kolmogorov wavenumber. The Lagrangian correlations of pressure gradient and velocity are compared and the Lagrangian time scale of the pressure gradient is observed to be much shorter than that of the velocity.
Pressure and higher-order spectra for homogeneous isotropic turbulence
NASA Technical Reports Server (NTRS)
Pullin, D. I.; Rogallo, R. S.
1994-01-01
The spectra of the pressure, and other higher-order quantities including the dissipation, the enstrophy, and the square of the longitudinal velocity derivative are computed using data obtained from direct numerical simulation of homogeneous isotropic turbulence at Taylor-Reynolds numbers R(sub lambda) in the range 38 - 170. For the pressure spectra we find reasonable collapse in the dissipation range (of the velocity spectrum) when scaled in Kolmogorov variables and some evidence, which is not conclusive, for the existence of a k(exp -7/3) inertial range where k = absolute value of K, is the modulus of the wavenumber. The power spectra of the dissipation, the enstrophy, and the square of the longitudinal velocity derivative separate in the dissipation range, but appear to converge together in the short inertial range of the simulations. A least-squares curve-fit in the dissipation range for one value of R(sub lambda) = 96 gives a form for the spectrum of the dissipation as k(exp 0)exp(-Ck eta), for k(eta) greater than 0.2, where eta is the Kolmogorov length and C is approximately equal to 2.5.
Testing a similarity theory for isotropic turbulence on DNS data.
NASA Astrophysics Data System (ADS)
Melander, Mogens; Fabijonas, Bruce
2006-11-01
Using direct numerical simulations, we consider the issue of self-similarity in 3D incompressible isotropic turbulence. The starting point for our investigation is a similarity theory we have developed on the basis of high Reynolds number shell model calculations. Like Kolmogorov's 1941 theory, our theory calls for similarity across all scales in the inertial range. Unlike K41, our theory does not fail on account of intermittency, but is developed to blossom in that environment. To observe self-similarity, it is essential that the correct variables are used, otherwise one sees only intermittency. The correct variables are reasonably easy to spot for the shell model, but they are more difficult to identify for the full Navier-Stokes equations. Moreover, one has to overcome the fact that the DNS has lower Reynolds numbers than in the shell model simulations so that the inertial range is shorter. Using the technique ESS, we clear this obstacle with only a minor modification to the theory. The DNS data then collapse on the theoretical pdf at all scales.
On the decay of homogeneous nearly isotropic turbulence behind active fractal grids
NASA Astrophysics Data System (ADS)
Thormann, Adrien; Meneveau, Charles
2012-11-01
The study of decaying isotropic turbulent flow is an important point of reference for turbulence theories and numerical simulations. For the past several decades, most experimental results appear to favor power-law decay with exponents between -1.2 and -1.4, approximately. More recently, fractal-generated turbulence (Hurst & Vassilicos, PoF 2007, and subsequent papers) using multi-scale passive grids suggest possible faster decay, and non-trivial behavior especially near the grid, where the mean velocity is spatially evolving. In order to generate spatially homogeneous flow using multi-scale injection of kinetic energy at high Reynolds numbers, we use a new type of active-grid consisting of winglets with various fractal shapes. We test space-filling fractal shaped winglets as well as Sierpisky-carpet and Apollonian packing type fractal shapes. Data are acquired using X-wire thermal anemometry. Tests of homogeneity of mean flow and turbulence intensity will be presented as well as decay of kinetic energy and spectral characteristics of the flow. This research is supported by NSF-CBET-1033942. The assistance of Ms. Imbi Salasoo and Mr. Nathan Greene in designing and building the fractal winglets is much appreciated. The authors also thank Mr. Vince Rolin for his assistance with the active grid.
Numerical studies of microscopic oil droplets under intense turbulence
NASA Astrophysics Data System (ADS)
Snyder, Murray; Knio, Omar
2008-11-01
The rise of small oil droplets in water experiencing isotropic turbulence conditions is analyzed computationally under four different turbulence intensities. The computational method combines DNS of the turbulent flow with Lagrangian tracking of the slightly buoyant droplets using a dynamical equation with buoyancy, virtual mass, pressure, drag, lift and history forces. In our recent work, Snyder et al. (2008), we showed that the puzzling behavior observed by Friedman and Katz (2002), where the rise velocity of droplets smaller than 800 μm in diameter is enhanced by turbulence whereas the rise of larger droplets is retarded, could be explained by significant drop then enhancement of the droplet drag coefficient, and corresponding drop in the virtual mass coefficient. In this study we use the same technique to explain the recent experimental results of Gopalan and Katz (2008), who also showed both suppression and enhancement of droplet rise velocities. Using drag and virtual mass coefficients which vary with Reynolds number, our computations approximate the experimental behavior observed by Gopalan and Katz in isotropic turbulence with 79, 100 and 151 μm Kolmogorov length scales. Combined with close agreement with the Friedman and Katz results, with an 88 μm Kolmogorov scale, our results provide further evidence that both the quasi-steady drag and virtual mass coefficients may be heavily modified under intense turbulence.
Local energy transfer and nonlocal interactions in homogeneous, isotropic turbulence
NASA Technical Reports Server (NTRS)
Domaradzki, J. Andrzej; Rogallo, Robert S.
1990-01-01
Detailed computations were made of energy transfer among the scales of motion in incompressible turbulent fields at low Reynolds numbers, generated by direct numerical simulations. It was observed that, although the transfer resulted from triad interactions that were nonlocal in k space, the energy always transferred locally. The energy transfer calculated from the eddy-damped quasi-normal Markovian (EDQNM) theory of turbulence at low Reynolds numbers is in excellent agreement with the results of the numerical simulations. At high Reynolds numbers, the EDQNM theory predicts the same transfer mechanism in the inertial range that is observed at low Reynolds numbers.
A detailed look at turbulence intensity
Technology Transfer Automated Retrieval System (TEKTRAN)
The effect of turbulence intensity on energy capture by small wind turbines has been a point of debate in the last few years. Claims of 25% de-rating of the power curve for turbines installed at sites with high turbulence are not uncommon. Over the years, many attempts have been made to model the ef...
Sweeping and straining effects in sound generation by high Reynolds number isotropic turbulence
NASA Technical Reports Server (NTRS)
Zhou, YE; Rubinstein, Robert
1995-01-01
The sound radiated by isotropic turbulence is computed using inertial range scaling expressions for the relevant two time and two point correlations. The result depends on whether the decay of Eulerian time correlations is dominated by large scale sweeping or by local straining: the straining hypothesis leads to an expression for total acoustic power, whereas the sweeping hypothesis leads to a more recent result.
Settling of almost neutrally buoyant particles in homogeneous isotropic turbulence
NASA Astrophysics Data System (ADS)
van Hinsberg, Michel; Clercx, Herman; Toschi, Federico
2015-11-01
Settling of particles in a turbulent flow occurs in various industrial and natural phenomena, examples are clouds and waste water treatment. It is well known that turbulence can enhance the settling velocity of particles. Many studies have been done, numerically and experimentally to investigate this behavior for the case of ``heavy'' particles, with particle to fluid density ratios above 100. Here we investigate the case of almost neutrally buoyant particles, i.e. density ratios between 1 and 100. In the case of light particles the Maxey-Riley equations cannot be simplified to only the Stokes drag and gravity force as pressure gradient, added mass and Basset history force are important as well. We investigate the influence of these forces on the settling velocity of particles and show that the extra forces can both increase or decrease the settling velocity, depending on the combination of the Stokes number and gravity applied.
Energy transfer in isotropic turbulence at low Reynolds numbers
NASA Technical Reports Server (NTRS)
Domaradzki, J. A.; Rogallo, R. S.
1988-01-01
Detailed measurements were made of energy transfer among the scales of motion in incompressible turbulent fields at low Reynolds numbers generated by direct numerical simulation. It was observed that although the transfer resulted from triad interactions that were non-local in k space, the energy always transferred locally. The results are consistent with the notion of non-uniform advection of small weak eddies by larger and stronger ones, similar to transfer processes in the far dissipation range at high Reynolds numbers.
Grid-generated isotropic homogeneous turbulence at high Reynolds numbers
NASA Technical Reports Server (NTRS)
Rosen, G.
1981-01-01
Consideration is given to an empirical formula for the longitudinal correlation function for grid-generated incompressible fluid turbulence at Reynolds numbers above 12,800. The formula, which relates the longitudinal correlation function to the inverse cube of a dimensionless geometrical ratio, is shown to minimize the global correlation integrals into which the two-point velocity correlation tensor has been substituted subject to a global constraint on the Sobolev concomitent of the longitudinal correlation function. Furthermore, the energy spectrum function associated with the empirical formula is shown to satisfy a tertiary Helmholtz-type linear condition throughout the initial period of decay.
Mixing of a passive scalar in isotropic and sheared homogeneous turbulence
NASA Technical Reports Server (NTRS)
Shirani, E.; Ferziger, J. H.; Reynolds, W. C.
1981-01-01
In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.
Theoretical and computational studies of isotropic homogeneous turbulence
NASA Astrophysics Data System (ADS)
Huang, Mei-Jiau
Numerical simulations are presented for viscous incompressible homogeneous turbulent flows with periodic boundary conditions. Our numerical method is based on the spectral Fourier method. Rogallo's code is modified and extended to trace fluid particles and simulate the evolution of material line elements. The first part of the thesis is about modifying and applying the code to simulate a passive vector field convected and stretched by the so called ABC flows in the presence of viscosity. The correlation of the geometry of the physical structures of the passive vector with the external straining is investigated. It is observed that most amplifications either occur in the neighborhoods of local unstable manifolds of the stagnation points of the ABC flows, if they exist, especially those with only one positive eigenvalue, or they are confined within the chaotic regions of the ABC flows if there is no stagnation point. The second part of the thesis is an investigation of the power-law energy decay of turbulence. Two decay exponents, 1.24 and 1.54, are measured from simulations. A new similarity form for the double and triple velocity autocorrelation functions using the Taylor microscale as the scaling, consistent with the Karman-Howarth equation and a power-law energy decay, is proposed and compared with numerical results. The proposed similarity form seems applicable at small to intermediate Reynolds number. For flows with very large Reynolds number, an expansion form of energy spectrum is proposed instead. The third part of the thesis is a presentation of the Lagrangian data collected by tracking fluid particles in decaying turbulent flows. The mean growth rates of the magnitudes of material line elements, that of the vorticity due to nonlinear forces, and the mean principal rates of strain tensors are found to be proportional to the square root of the mean enstrophy. The proportional coefficients remain constant during the decay. The mean angles between material
Theoretical and Computational Studies of Isotropic Homogeneous Turbulence
NASA Astrophysics Data System (ADS)
Huang, Mei-Jiau
Numerical simulations are presented for viscous incompressible homogeneous turbulent flows with periodic boundary conditions. Our numerical method is based on the spectral Fourier method. Rogallo's code is modified and extended to trace fluid particles and simulate the evolution of material line elements. The first part of the thesis is about modifying and applying the code to simulate a passive vector field convected and stretched by the so-called ABC flows in the presence of viscosity. The correlation of the geometry of the physical structures of the passive vector with the external straining is investigated. It is observed that most amplifications either occur in the neighborhoods of local unstable manifolds of the stagnation points of the ABC flows, if they exist, especially those with only one positive eigenvalue, or they are confined within the chaotic regions of the ABC flows if there is no stagnation point. The second part of the thesis is an investigation of the power-law energy decay of turbulence. Two decay exponents, 1.24 and 1.54, are measured from simulations. A new similarity form for the double and triple velocity autocorrelation functions using the Taylor microscale as the scaling, consistent with the Karman-Howarth equation and a power-law energy decay, is proposed and compared with numerical results. The proposed similarity form seems applicable at small to intermediate Reynolds number. For flows with very large Reynolds number, an expansion form of energy spectrum is proposed instead. The third part of the thesis is a presentation of the Lagrangian data collected by tracking fluid particles in decaying turbulent flows. The mean growth rates of the magnitudes of material line elements, that of the vorticity due to nonlinear forces, and the mean principal rates of strain tensors are found to be proportional to the square root of the mean enstrophy. The proportional coefficients remain constant during the decay. The mean angles between material
The radiated noise from isotropic turbulence and heated jets
NASA Technical Reports Server (NTRS)
Lilley, G. M.
1995-01-01
Our understanding of aerodynamic noise has its foundations in the work of Sir James Lighthill (1952), which was the first major advance in acoustics since the pioneering work of Lord Rayleigh in the last century. The combination of Lighthill's theory of aerodynamic noise as applied to turbulent flows and the experimental growing database from the early 1950's was quickly exploited by various jet propulsion engine designers in reducing the noise of jet engines at takeoff and landing to levels marginally acceptable to communities living in the neighborhoods of airports. The success in this noise containment led to the rapid growth of fast economical subsonic civil transport aircraft worldwide throughout the 1960's and has continued to the present day. One important factor in this success story has been the improvements in the engine cycle that have led to both reductions in specific fuel consumption and noise. The second is the introduction of Noise Certification, which specifies the maximum noise levels at takeoff and landing that all aircraft must meet before they can be entered on the Civil Aircraft Register. The growing interest in the development of a new supersonic civil transport to replace 'Concorde' in the early years of the next century has led to a resurgence of interest in the more challenging problem of predicting the noise of hot supersonic jets and developing means of aircraft noise reduction at takeoff and landing to meet the standards now accepted for subsonic Noise Certification. The prediction of aircraft noise to the accuracy required to meet Noise Certification requirements has necessitated reliance upon experimental measurements and empirically derived laws based on the available experimental data bases. These laws have their foundation in the results from Lighthill's theory, but in the case of jet noise, where the noise is generated in the turbulent mixing region with the external ambient fluid, the complexity of the turbulent motion has
Short-time evolution of Lagrangian velocity gradient correlations in isotropic turbulence
NASA Astrophysics Data System (ADS)
Fang, L.; Bos, W. J. T.; Jin, G. D.
2015-12-01
We show by direct numerical simulation (DNS) that the Lagrangian cross correlation of velocity gradients in homogeneous isotropic turbulence increases at short times, whereas its auto-correlation decreases. Kinematic considerations allow to show that two invariants of the turbulent velocity field determine the short-time velocity gradient correlations. In order to get a more intuitive understanding of the dynamics for longer times, heuristic models are proposed involving the combined action of local shear and rotation. These models quantitatively reproduce the effects and disentangle the different physical mechanisms leading to the observations in the DNS.
How isotropic are turbulent flows generated by using periodic conditions in a cube?
NASA Astrophysics Data System (ADS)
Qin, Z. C.; Fang, L.; Fang, J.
2016-03-01
In numerical simulations, "isotropic" turbulent flows are always generated by using periodic conditions. We show that these periodic conditions mathematically lead to large-scale anisotropy which can be about 10% of the mean values, and thus prevent existing post-processing results from being accurate. A decomposition method by employing spherical harmonics is then proposed to distinguish this scale-dependent anisotropy effect from others and to minimize the related post-processing error.
New forcing scheme to sustain particle-laden homogeneous and isotropic turbulence
NASA Astrophysics Data System (ADS)
Mallouppas, G.; George, W. K.; van Wachem, B. G. M.
2013-08-01
This paper evaluates the effect of forcing to sustain turbulence on the transfer function of the fluid with particles suspended in a homogeneous and isotropic flow. As mentioned by Lucci et al. ["Modulation of isotropic turbulence by particles of Taylor length-scale size," J. Fluid Mech. 650, 5-55 (2010), 10.1017/S0022112009994022], there are three limitations of forcing particle-laden homogeneous and isotropic turbulence: (a) large fluctuations on the temporal evolution of the kinetic energy are created when forcing is active at low wavenumbers, (b) the redistribution of energy is affected when forcing is performed over all wavenumbers, and (c) the nonlinear transfer function of the fluid due to the triadic interactions is affected when forcing is active over a wavenumber range. These limitations make the interpretation of the effects of particles on the energy spectrum of the fluid difficult. A new forcing scheme in physical space has been designed which avoids these limitations in wavenumber space, so the spectral effects of particles can be evaluated. The performance of this forcing scheme is tested using Direct Numerical Simulations. It is shown that the nonlinear transfer function of the fluid with the current forcing scheme is only affected at the wavenumbers it is acting, consistent with the theory. Even so, the spatial coherence and phase spectra between the two-way coupling and the fluid computed from the simulations show that the new forcing scheme is only moderately correlated even for the forced wavenumbers, with correlation coefficient typically about 10%.
NASA Technical Reports Server (NTRS)
Ling, S. C.; Saad, A.
1977-01-01
The energetic isotropic turbulence generated by a waterfall of low head was found to be developed in part through the unstable two-phase flow of entrained air bubbles. The resulting turbulent field had a turbulent Reynolds number in excess of 20,000 and maintained a self-similar structure throughout the decay period studied. The present study may provide some insight into the structure of turbulence produced by breaking waves over the ocean.
NASA Astrophysics Data System (ADS)
Canet, Léonie; Delamotte, Bertrand; Wschebor, Nicolás
2016-06-01
We investigate the regime of fully developed homogeneous and isotropic turbulence of the Navier-Stokes (NS) equation in the presence of a stochastic forcing, using the nonperturbative (functional) renormalization group (NPRG). Within a simple approximation based on symmetries, we obtain the fixed-point solution of the NPRG flow equations that corresponds to fully developed turbulence both in d =2 and 3 dimensions. Deviations to the dimensional scalings (Kolmogorov in d =3 or Kraichnan-Batchelor in d =2 ) are found for the two-point functions. To further analyze these deviations, we derive exact flow equations in the large wave-number limit, and show that the fixed point does not entail the usual scale invariance, thereby identifying the mechanism for the emergence of intermittency within the NPRG framework. The purpose of this work is to provide a detailed basis for NPRG studies of NS turbulence; the determination of the ensuing intermittency exponents is left for future work.
Suppression of self-organized structure coarsening in homogenous isotropic turbulence
NASA Astrophysics Data System (ADS)
Takagi, Youhei
2014-11-01
Self-organized structure by spinodal decomposition is often seen in quenched binary mixture. Complex network structure is formed through coarsening process of self-organized structure when the phase separation due to spinodal decomposition proceeds. The phase separation governed by the Cahn-Hilliard equation have been well investigated for stationary fluid in previous studies, however, the turbulent effect on the formation of structures was not fully discussed. In this study, we carried out a numerical simulation for homogenous isotropic turbulence with phase separation, the relation between turbulent vortex formation and self-organized structure coarsening. The governing equations are incompressible Navier-Stokes equation considering phase separation force and Cahn-Hilliard equation with the chemical potential based on the Landau-Ginzburg free energy. From the identification and visualization of turbulent structures, it was found that the local entrainment of small eddy structure suppressed the coarsening process of self-organized structure. The energy used in phase separation was related to the initial process of vortex sheet-tube transition in turbulent flow, and the energy cascade from large turbulent structure to small eddy was different from that without phase separation.
Turbulence generation through intense kinetic energy sources
NASA Astrophysics Data System (ADS)
Maqui, Agustin F.; Donzis, Diego A.
2016-06-01
Direct numerical simulations (DNS) are used to systematically study the development and establishment of turbulence when the flow is initialized with concentrated regions of intense kinetic energy. This resembles both active and passive grids which have been extensively used to generate and study turbulence in laboratories at different Reynolds numbers and with different characteristics, such as the degree of isotropy and homogeneity. A large DNS database was generated covering a wide range of initial conditions with a focus on perturbations with some directional preference, a condition found in active jet grids and passive grids passed through a contraction as well as a new type of active grid inspired by the experimental use of lasers to photo-excite the molecules that comprise the fluid. The DNS database is used to assert under what conditions the flow becomes turbulent and if so, the time required for this to occur. We identify a natural time scale of the problem which indicates the onset of turbulence and a single Reynolds number based exclusively on initial conditions which controls the evolution of the flow. It is found that a minimum Reynolds number is needed for the flow to evolve towards fully developed turbulence. An extensive analysis of single and two point statistics, velocity as well as spectral dynamics and anisotropy measures is presented to characterize the evolution of the flow towards realistic turbulence.
Optical intensity interferometry through atmospheric turbulence
NASA Astrophysics Data System (ADS)
Tan, P. K.; Chan, A. H.; Kurtsiefer, C.
2016-04-01
Conventional ground-based astronomical observations suffer from image distortion due to atmospheric turbulence. This can be minimized by choosing suitable geographic locations or adaptive optical techniques, and avoided altogether by using orbital platforms outside the atmosphere. One of the promises of optical intensity interferometry is its independence from atmospherically induced phase fluctuations. By performing narrow-band spectral filtering on sunlight and conducting temporal intensity interferometry using actively quenched avalanche photodiodes, the Solar g(2)(τ) signature was directly measured. We observe an averaged photon bunching signal of g(2)(τ) = 1.693 ± 0.003 from the Sun, consistently throughout the day despite fluctuating weather conditions, cloud cover and elevation angle. This demonstrates the robustness of the intensity interferometry technique against atmospheric turbulence and opto-mechanical instabilities, and the feasibility to implement measurement schemes with both large baselines and long integration times.
On the dynamics of small-scale vorticity in isotropic turbulence
NASA Technical Reports Server (NTRS)
Jimenez, Javier; Wray, A. A.
1994-01-01
It was previously shown that the strong vorticity in isotropic turbulence is organized into tubular vortices ('worms') whose properties were characterized through the use of full numerical simulations at several Reynolds numbers. At the time most of the observations were kinematic, and several scaling laws were discovered for which there was no theoretical explanation. In the meantime, further analysis of the same fields yielded new information on the generation of the vortices, and it was realized that even if they had to be formed by stretching, they were at any given moment actually compressed at many points of their axes. This apparent contradiction was partially explained by postulating axial inertial waves induced by the nonuniformity of the vortex cores, which helped to 'spread' the axial strain and allowed the vortices to remain compact even if not uniformly stretched. The existence of such solutions was recently proved numerically. The present report discusses a set of new numerical simulations of isotropic turbulence, and a reanalysis of the old ones, in an effort to prove or disprove the presence of these waves in actual turbulent flows and to understand the dynamics, as opposed to the kinematics, of the vortices.
The isotropic nature of the background turbulence spectra in the solar wind
NASA Astrophysics Data System (ADS)
Wang, X.; Tu, C. Y.; He, J.; Marsch, E.; Wang, L.
2014-12-01
At the high-frequency end of the inertial range, the solar wind turbulence power spectrum was recently found to be anisotropic with respect to the direction of local magnetic field, as an evidence for the presence of a "critical balance" style turbulence cascade. However, we find that the spectral anisotropy seems to result from intermittent structures. The following two independent studies corroborate this statement by showing that the power spectra of the background turbulence, in which there are no intermittent structures, have an isotropic nature. In Study 1, we remove the wavelet coefficients of the local intermittency with large partial variance increment (PVI), and find the spectral indices of the magnetic field are 1.63±0.02, independent of the angle θRB between the direction of the local background magnetic field and the radial direction. In Study 2, we make a statistical study on the magnetic field spectral indices obtained by using Fast Fourier Transform on 40 time series, in which no intermittent structures appear. We find that for the time series with 0o<θRB <6o, the probability distribution of the observed spectral indices peaks at -1.7, while the -2 index predicted by the "critical balance" theory rarely appears. For the time series with 84 o <θRB <90 o, the probability distribution of the indices peaks at -1.5. Considering the uncertainty of the statistics, these results show that the background-turbulence spectra are nearly isotropic with respect to θRB, which may be consistent with some explanations based on hydrodynamic turbulence theory.
NASA Astrophysics Data System (ADS)
Xanthos, Savvas; Gong, Minwei; Andreopoulos, Yiannis
2010-01-01
Further analysis of the experimental data of the velocity gradient tensor first published by Xanthos et al. [J. Fluid Mech. 584, 301 (2007)] has been carried out and new results are reported here to provide additional insights on the effects of expansion waves interacting with isotropic turbulence. The flow field was generated by the reflection of an incoming shock wave at the open end of a large scale shock tube facility which interacted with the induced flow behind the incident shock wave which passed through a turbulence generating grid. In the present configuration the interaction is free from streamline curvature effects, which cause additional effects on turbulence. The strength of the applied expansive straining was 240 s-1. Rectangular pattern grids of different mesh sizes were used to generate isotropic and homogeneous turbulence with turbulent Reynolds number Reλ based on Taylor's microscale between 450 and 488. Lateral vorticity fluctuations and fluctuations of enstrophy and all stretching vector components are drastically reduced during the interaction. Residual attenuation in the postinteraction flow field was found only in the lateral vorticity fluctuations and in the longitudinal stretching term S11Ω1. Helicity and the helicity angle were computed from the data and the orientation angle of the vorticity vector in reference to the velocity vector was determined. Large fluctuations of the helicity angle were observed which extend from 0° to 180° with most probable values close to 30° and 130° and a mean value of 85°. Rotational dissipation rate was found to be high at these angles. The time-dependent signals of enstrophy, vortex stretching/tilting vector, and dissipation rate were found to exhibit a rather strong intermittent behavior which is characterized by high amplitude bursts followed by low level activities. It was found that the observed strong dissipative events are mostly associated with strong activities in the longitudinal stretching
Fang, L; Zhang, Y J; Fang, J; Zhu, Y
2016-08-01
We show by direct numerical simulations (DNSs) that in different types of isotropic turbulence, the fourth-order statistical invariants have approximately a linear relation, which can be represented by a straight line in the phase plane, passing two extreme states: the Gaussian state and the restricted Euler state. Also, each DNS case corresponds to an equilibrium region that is roughly Reynolds-dependent. In addition, both the time reversal and the compressibility effect lead to nonequilibrium transition processes in this phase plane. This observation adds a new restriction on the mean-field theory. PMID:27627399
The upper atmosphere of Uranus - A critical test of isotropic turbulence models
NASA Technical Reports Server (NTRS)
French, R. G.; Elliot, J. L.; Sicardy, B.; Nicholson, P.; Matthews, K.
1982-01-01
Observations of the August 15, 1980, Uranus occultation of KM 12, obtained from Cerro Tololo InterAmerican Observatory, European Southern Observatory, and Cerro Las Campanas Observatory, are used to compare the atmospheric structure at points separated by approximately 140 km along the planetary limb. The results reveal striking, but by no means perfect correlation of the light curves, ruling out isotropic turbulence as the cause of the light curve spikes. The atmosphere is strongly layered, and any acceptable turbulence model must accommodate the axial ratios of greater than about 60 which are observed. The mean temperature of the atmosphere is 150 plus or minus 15 K for the region near number density 10 to the 14th per cu cm. Derived temperature variations of vertical scale approximately 130 km and amplitude plus or minus 5 K are in agreement for all stations, and correlated spikes correspond to low-amplitude temperature variations with a vertical scale of several kilometers.
Modulation of isotropic turbulence by deformable droplets of Taylor lengthscale size
NASA Astrophysics Data System (ADS)
Dodd, Michael; Ferrante, Antonino
2014-11-01
We investigate the effects of finite-size deformable droplets on decaying isotropic turbulence via direct numerical simulation (DNS). DNS is performed using the two-fluid pressure-correction method by Dodd and Ferrante [J. Comput. Phys. 273, 416 (2014)] coupled with the volume of fluid method by Baraldi et al. [Comput. & Fluids 96, 322 (2014)]. We fully-resolve the flow around and inside 3130 droplets of Taylor lengthscale size, resulting in a droplet volume fraction of 0.05. The initial Taylor lengthscale Reynolds number is Reλ0 = 75 , and the computational mesh has 10243 grid points. We analyze the effects on turbulence modulation of varying the droplet- to carrier-fluid viscosity ratio (1 <=μd /μc <= 100) and the droplet Weber number based on the r.m.s velocity of turbulence (0 . 1 <=Werms <= 5). We discuss how varying these parameters affects the turbulence kinetic energy budget, and explain the physical mechanisms for such modulation. This work was supported by the National Science Foundation CAREER Award, Grant Number OCI-1054591.
Model for drop coalescence in a locally isotropic turbulent flow field.
Narsimhan, Ganesan
2004-04-01
The proposed model views drop coalescence in a turbulent flow field as a two-step process consisting of formation of a doublet due to drop collisions followed by coalescence of the individual droplets in a doublet due to the drainage of the intervening film of continuous phase under the action of colloidal (van der Waals and electrostatic) and random turbulent forces. The turbulent flow field was assumed to be locally isotropic. A first-passage-time analysis was employed for the random process of intervening continuous-phase film thickness between the two drops of a doublet in order to evaluate the first two moments of coalescence-time distribution of the doublet. The average drop coalescence time of the doublet was dependent on the barrier for coalescence due to the net repulsive force (net effect of colloidal repulsive and turbulent attractive forces). The predicted average drop coalescence time was found to be smaller for larger turbulent energy dissipation rates, smaller surface potentials, larger drop sizes, larger ionic strengths, and larger drop size ratios of unequal-sized drop pairs. The predicted average drop coalescence time was found to decrease whenever the ratio of average turbulent force to repulsive force barrier became larger. The calculated coalescence-time distribution was broader, with a higher standard deviation, at lower energy dissipation rates, higher surface potentials, smaller drop sizes, and smaller size ratios of unequal drop pairs. The model predictions of average coalescence-rate constants for tetradecane-in-water emulsions stabilized by sodium dodecyl sulfate (SDS) in a high-pressure homogenizer agreed fairly well with the inferred experimental values as reported by Narsimhan and Goel (J. Colloid Interface Sci. 238 (2001) 420-432) at different homogenizer pressures and SDS concentrations. PMID:14985038
Magnetic Field Line Random Walk in Isotropic Turbulence with Varying Mean Field
NASA Astrophysics Data System (ADS)
Sonsrettee, W.; Subedi, P.; Ruffolo, D.; Matthaeus, W. H.; Snodin, A. P.; Wongpan, P.; Chuychai, P.; Rowlands, G.; Vyas, S.
2016-08-01
In astrophysical plasmas, the magnetic field line random walk (FLRW) plays an important role in guiding particle transport. The FLRW behavior is scaled by the Kubo number R=(b/{B}0)({{\\ell }}\\parallel /{{\\ell }}\\perp ) for rms magnetic fluctuation b, large-scale mean field {{\\boldsymbol{B}}}0, and coherence scales parallel ({{\\ell }}\\parallel ) and perpendicular ({{\\ell }}\\perp ) to {{\\boldsymbol{B}}}0. Here we use a nonperturbative analytic framework based on Corrsin’s hypothesis, together with direct computer simulations, to examine the R-scaling of the FLRW for varying B 0 with finite b and isotropic fluctuations with {{\\ell }}\\parallel /{{\\ell }}\\perp =1, instead of the well-studied route of varying {{\\ell }}\\parallel /{{\\ell }}\\perp for b \\ll {B}0. The FLRW for isotropic magnetic fluctuations is also of astrophysical interest regarding transport processes in the interstellar medium. With a mean field, fluctuations may have variance anisotropy, so we consider limiting cases of isotropic variance and transverse variance (with b z = 0). We obtain analytic theories, and closed-form solutions for extreme cases. Padé approximants are provided to interpolate all versions of theory and simulations to any B 0. We demonstrate that, for isotropic turbulence, Corrsin-based theories generally work well, and with increasing R there is a transition from quasilinear to Bohm diffusion. This holds even with b z = 0, when different routes to R\\to ∞ are mathematically equivalent; in contrast with previous studies, we find that a Corrsin-based theory with random ballistic decorrelation works well even up to R = 400, where the effects of trapping are barely perceptible in simulation results.
Effect of freestream isotropic turbulence on heat transfer from a sphere
NASA Astrophysics Data System (ADS)
Bagchi, Prosenjit; Kottam, Kirit
2008-07-01
We consider direct numerical simulation (DNS) based on pseudospectral methods to study the heat transfer around a stationary sphere held at a constant temperature and subject to an ambient turbulent velocity and temperature condition. The sphere Reynolds number is in the range of 63-400, and the sphere diameter (d) varies from one to eight times the Kolmogorov scale (η). The ambient turbulent field is isotropic, and the Taylor microscale Reynolds number Rλ varies from 38 to 240. Results from two sets of DNS are presented. In the first set, the ambient velocity field is turbulent, but the ambient temperature is held constant. In the second set of simulations, both the ambient velocity and the temperature fields are turbulent. These two sets of simulations allow us to isolate the role of freestream velocity fluctuations and temperature fluctuations in modifying the mean and time-dependent heat transfer from the sphere. The mean Nusselt number is observed to be independent of Rλ. It is shown that the freestream turbulence does not have any significant effect on the mean Nusselt number, and the available correlations for a steady and uniform ambient can predict the mean Nusselt number under the turbulent ambient condition. The instantaneous Nusselt number, however, can differ significantly from the mean, and can be negative in case of large temperature fluctuation in the far field. The instantaneous Nusselt number obtained from the DNS is analyzed and compared with the analytical expression for the unsteady thermal response of a sphere. It is shown that the thermal added-mass effect is small for d /η≈1 but introduces spurious oscillation at higher d. The thermal history effect is shown to be insignificant for all d /η. Properties of the thermal wake in the presence of the turbulent velocity and temperature fields are studied. The mean thermal wake is observed to be shorter in streamwise direction and wider in crossflow direction in a turbulent ambient than that
Creating Only Isotropic Homogeneous Turbulence in Liquid Helium near Absolute Zero
NASA Astrophysics Data System (ADS)
Ihas, G. G.; Thompson, K. J.; Labbe, G.; McClintock, P. V. E.
2012-02-01
Flow through a grid is a standard method to produce isotropic, homogeneous turbulence for laboratory study. This technique has been used to generate quantum turbulence (QT) above 1 K in superfluid heliumootnotetextS. R. Stalp, L. Skrbek, and R. J. Donnelly, Phys. Rev. Lett. 82, 4831 (1999). where QT seems to mimic classical turbulence. Efforts have been made recentlyootnotetextG. G. Ihas, G. Labbe, S-c. Liu, and K. J. Thompson, J. Low Temp. Phys. 150, 384 (2008). to make similar measurements near absolute zero, where there is an almost total absence of normal fluid and hence classical viscosity. This presents the difficulty that most motive force devices produce heat which overwhelms the phenomena being investigated. The process of designing and implimenting a ``dissipation-free'' motor for pulling a grid through superfluid helium at millikelvin temperatures has resulted in the development of new techniques which have broad application in low temperature research. Some of these, such as Meissner-affect magnetic drives, capacitive and inductive position sensors, and magnetic centering devices will be described. Heating results for devices which can move in a controlled fashion from very low speed up to 10 cm/s will be presented. Acknowledgement: We thank W.F. Vinen for many useful discussions.
NASA Astrophysics Data System (ADS)
Hammond, Adam; Dou, Zhongwang; Tripathi, Anjan; Liang, Zach; Meng, Hui
2015-11-01
Study of droplet collision and cloud formation should consider the effects of both turbulence and electrostatic charge on particle dynamics. We present the first experimental observation of radial relative velocity (RV) of charged particles in homogeneous and isotropic turbulence (HIT). Charges on particles were generated through triboelectric effect between the inner wall of the chamber and the particles. To measure charge distribution, a particle-laden head-on impinging flow mimicking our HIT chamber conditions was built and holographic particle tracking was applied to quantify particle charges by measuring their displacements in an electric field. Particles were observed to have opposite charges. Next, in our HIT chamber, we measured particle RV by a novel 4-frame particle tracking velocimetry technique with and without charges on particles, wherein charges were neutralized by coating the interior of the HIT chamber with conductive carbon paint. We compared RV under the same turbulence conditions between charged particles and neutral particles and observed that when particles were oppositely charged, their mean inward RV increased at small separation distances. This result is consistent with recent theory and simulations (Lu and Shaw, Physics of Fluids, 2015). This work was supported by the National Science Foundation through a Collaborative Research Grant CBET-0967407.
NASA Astrophysics Data System (ADS)
Farge, Marie; Roussel, Olivier; Schneider, Kai
2004-11-01
We compare the extraction of coherent vortices in 3D homogeneous isotropic turbulence computed by DNS using either orthogonal or biorthogonal wavelets. The method is based on a wavelet decomposition of the vorticity field and a subsequent thresholding of the wavelet coefficients (PRL, 87(5), 2001, Phys. Fluids 15(10), 2003). The coherent vorticity is reconstructed from few strong wavelet coefficients while the incoherent vorticity is reconstructed from the remaining weak coefficients. In the orthogonal case the choice of the threshold is motivated from statistical denoising theory and has no adjustable parameters. Using 3% of the coefficients we show that both decompositions extract the coherent vortices out of the turbulent flow. They contain 99.6% of the energy and retain 74% and 68% of the enstrophy in the orthogonal and biorthogonal case, respectively. Concerning the incoherent background flow, it is structureless and decorrelated for the orthogonal decomposition, with a Gaussian velocity PDF. In contrast, the biorthogonal decomposition yields a background flow which exhibits quasi-2D sheet-like structures with an exponetial velocity PDF instead. In conclusion, modeling the incoherent background flow might be more difficult using biorthogonal wavelets for the CVS (Coherent Vortex Simulation, cf. Flow, Turbulence and Combustion 66(4), 2001).
Presumed PDF Modeling of Early Flame Propagation in Moderate to Intense Turbulence Environments
NASA Technical Reports Server (NTRS)
Carmen, Christina; Feikema, Douglas A.
2003-01-01
The present paper describes the results obtained from a one-dimensional time dependent numerical technique that simulates early flame propagation in a moderate to intense turbulent environment. Attention is focused on the development of a spark-ignited, premixed, lean methane/air mixture with the unsteady spherical flame propagating in homogeneous and isotropic turbulence. A Monte-Carlo particle tracking method, based upon the method of fractional steps, is utilized to simulate the phenomena represented by a probability density function (PDF) transport equation. Gaussian distributions of fluctuating velocity and fuel concentration are prescribed. Attention is focused on three primary parameters that influence the initial flame kernel growth: the detailed ignition system characteristics, the mixture composition, and the nature of the flow field. The computational results of moderate and intense isotropic turbulence suggests that flames within the distributed reaction zone are not as vulnerable, as traditionally believed, to the adverse effects of increased turbulence intensity. It is also shown that the magnitude of the flame front thickness significantly impacts the turbulent consumption flame speed. Flame conditions studied have fuel equivalence ratio s in the range phi = 0.6 to 0.9 at standard temperature and pressure.
NASA Astrophysics Data System (ADS)
Borgas, Michael S.; Yeung, P. K.
2004-03-01
A new model for Lagrangian particle-pair separation in turbulent flows is developed and compared with data from direct numerical simulations (DNS) of isotropic turbulence. The model formulation emphasizes (i) non-Gaussian behaviour in Eulerian and Lagrangian statistics, (ii) the occurrence of large separation velocities, (iii) the role of straining and streaming flow structure as recognized in kinematic simulations of turbulence, and (iv) the role of conditionally averaged accelerations in stochastic modelling of turbulent relative dispersion. Previous stochastic models of relative dispersion have produced unrealistic behaviour, particularly in the dissipation subrange where viscous effects are important, which have led to questions on the adequacy of stochastic modelling. However, this failure can now be recognized as inadequate detail in formulation, which is explained and rectified in this paper. The model is quasi-one-dimensional in nature, and is focused on the statistics of particle-pair separation distance and its rate of change, referred to as the separation speed. Detailed comparisons are presented at several Reynolds numbers using the DNS database reported in a companion paper (Part 1). Up to fourth-order moments for these quantities are examined, as well as the separation-distance probability density function, which is discussed in the context of recent claims of Richardson scaling in the literature. The model is able to account for the spatial representation of straining regions as well as incompressibility of the flow, and is shown to reproduce strong non-Gaussianity and intermittency in the Lagrangian separation statistics observed in DNS. Comparisons with recent physical experiments are also remarkably consistent. This work demonstrates that stochastic models when properly formulated are effective and efficient representations of the dispersion process and this general class of models therefore possess great utility for calculations of both one
Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation
NASA Technical Reports Server (NTRS)
Swanson, R. C.; Rumsey, Christopher L.; Rubinstein, Robert; Balakumar, Ponnampalam; Zang, Thomas A.
2012-01-01
Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.
Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence
NASA Astrophysics Data System (ADS)
Bragg, Andrew D.; Ireland, Peter J.; Collins, Lance R.
2016-01-01
In this paper, we investigate both theoretically and numerically the Forward-In-Time (FIT) and Backward-In-Time (BIT) dispersion of fluid and inertial particle-pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dispersion is in general much stronger than that for fluid particles. For example, the ratio of the BIT to FIT mean-square separation can be up to an order of magnitude larger for the inertial particles than for the fluid particles. We also find that for both the inertial and fluid particles, the irreversibility becomes stronger as the scale of their separation decreases. Regarding the physical mechanism for the irreversibility, we argue that whereas the irreversibility of fluid particle-pair dispersion can be understood in terms of a directional bias arising from the energy transfer process in turbulence, inertial particles experience an additional source of irreversibility arising from the non-local contribution to their velocity dynamics, a contribution that vanishes in the limit St → 0, where St is the particle Stokes number. For each given initial (final, in the BIT case) separation, r0, there is an optimum value of St for which the dispersion irreversibility is strongest, as such particles are optimally affected by both sources of irreversibility. We derive analytical expressions for the BIT, mean-square separation of inertial particles and compare the predictions with numerical data obtained from a Reλ ≈ 582 (where Reλ is the Taylor Reynolds number) Direct Numerical Simulation (DNS) of particle-laden isotropic turbulent flow. The small-time theory, which in the dissipation range is valid for times ≤max[Stτη, τη] (where τη is the Kolmogorov time scale), is in excellent agreement with the DNS. The theory for long-times is in good agreement with the DNS
Isotropic boundary adapted wavelets for coherent vorticity extraction in turbulent channel flows
NASA Astrophysics Data System (ADS)
Farge, Marie; Sakurai, Teluo; Yoshimatsu, Katsunori; Schneider, Kai; Morishita, Koji; Ishihara, Takashi
2015-11-01
We present a construction of isotropic boundary adapted wavelets, which are orthogonal and yield a multi-resolution analysis. We analyze DNS data of turbulent channel flow computed at a friction-velocity based Reynolds number of 395 and investigate the role of coherent vorticity. Thresholding of the wavelet coefficients allows to split the flow into two parts, coherent and incoherent vorticity. The statistics of the former, i.e., energy and enstrophy spectra, are close to the ones of the total flow, and moreover the nonlinear energy budgets are well preserved. The remaining incoherent part, represented by the large majority of the weak wavelet coefficients, corresponds to a structureless, i.e., noise-like, background flow and exhibits an almost equi-distribution of energy.
NASA Astrophysics Data System (ADS)
Goto, Susumu; Kida, Shigeo
1999-07-01
The Lagrangian direct-interaction approximation developed previously by the present authors [S. Kida and S. Goto, J. Fluid Mech. 345, 307 (1997)] is applied to a passive scalar field in isotropic turbulence. We examine the behavior of solutions to the resultant closure equations for the correlation function of the scalar field for arbitrary values of the Schmidt number, and show systematically that the solutions are completely consistent with the phenomenological theories on the scalar spectral function by Obukhov (1949), Corrsin (1951), Batchelor et al. (1959), and Batchelor (1959). The universal forms of the function in the statistically stationary state are obtained by solving the closure equations numerically in the whole wave number range for each case of moderate, extremely large, and small values of the Schmidt number.
He, G.; Doolen, G.D.; Chen, S.
1999-12-01
The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations. {copyright} {ital 1999 American Institute of Physics.}
Passive scalar convective-diffusive subrange for low Prandtl numbers in isotropic turbulence
NASA Astrophysics Data System (ADS)
Briard, A.; Gomez, T.
2015-01-01
In this Rapid Communication, we study the behavior of a strongly diffusive passive scalar field T submitted to a freely decaying, homogeneous and isotropic turbulence with eddy-damped quasinormal Markovian simulations. We present a new subrange located between the k-17 /3 inertial-diffusive subrange and the Kolmogorov wave number kη. This subrange is generated by small-scale convection linked to kη that balances diffusion effects. Thus, we build a typical length scale kCD -1 based on convection and diffusion and give an expression for the shape of the passive scalar spectrum in this subrange ET˜√{Pr}k-11 /3 using physical arguments. This result unifies two different theories coming from Batchelor [G. K. Batchelor, J. Fluid. Mech. 5, 113 (1959), 10.1017/S002211205900009X] and Chasnov [J. Chasnov et al., Phys. Fluids A 1, 1698 (1989), 10.1063/1.857535] and explains results previously obtained experimentally.
NASA Astrophysics Data System (ADS)
Gotoh, Toshiyuki
2012-11-01
Spectrum of passive scalar variance at very high Schmidt number up to 1000 in isotropic steady turbulence has been studied by using very high resolution DNS. Gaussian random force and scalar source which are isotropic and white in time are applied at low wavenumber band. Since the Schmidt number is very large, the system was integrated for 72 large eddy turn over time for the system to forgot the initial state. It is found that the scalar spectrum attains the asymptotic k-1 spectrum in the viscous-convective range and the constant CB is found to be 5.7 which is larger than 4.9 obtained by DNS under the uniform mean scalar gradient. Reasons for the difference are inferred as the Reynolds number effect, anisotropy, difference in the scalar injection, duration of time average, and the universality of the constant is discussed. The constant CB is also compared with the prediction by the Lagrangian statistical theory for the passive scalar. The scalar spectrum in the far diffusive range is found to be exponential, which is consistent with the Kraichnan's spectrum. However, the Kraichnan spectrum was derived under the assumption that the velocity field is white in time, therefore theoretical explanation of the agreement needs to be explored. Grant-in-Aid for Scientific Research No. 21360082, Ministry of Education, Culture, Sports, Science and Technology of Japan.
Large fluctuations of the nonlinearities in isotropic turbulence. Anisotropic filtering analysis
NASA Astrophysics Data System (ADS)
Tordella, D.; Di Savino, S.; Sitzia, L.
2014-09-01
Using a Navier-Stokes isotropic turbulent field numerically simulated in a box with a discretization of 10243 (Biferale et al., 2005), we show that the probability of having a stretching-tilting larger than a few times the local enstrophy is low. By using an anisotropic kind of filter in the Fourier space, where wavenumbers that have at least one component below a threshold or inside a range are removed, we analyze these survival statistics when the large, the small inertial or the small inertial and dissipation scales are filtered out. By considering a flow obtained by randomizing the phases of the Fourier modes, and applying our filtering techniques, we identified clearly the properties attributable to turbulence. It can be observed that, in the unfiltered isotropic Navier-Stokes field, the probability of the ratio (|ωṡ∇U|/|) being higher than a given threshold is higher than in the fields where the large scales were filtered out. At the same time, it is lower than in the fields where the small inertial and dissipation range of scales is filtered out. This is basically due to the suppression of compact structures in the ranges that have been filtered in different ways. The partial removal of the background of filaments and sheets does not have a first order effect on these statistics. These results are discussed in the light of a hypothesized relation between vortical filaments, sheets and blobs in physical space and in Fourier space. The study in fact can be viewed as a kind of test for this idea and tries to highlight its limits. We conclude that a qualitative relation in physical space and in Fourier space can be supposed to exist for blobs only. That is for the near isotropic structures which are sufficiently described by a single spatial scale and do not suffer from the disambiguation problem as filaments and sheets do. Information is also given on the filtering effect on statistics concerning the inclination of the strain rate tensor eigenvectors with
NASA Astrophysics Data System (ADS)
Dou, Zhongwang; Pecenak, Zachary; Liang, Zach; Cao, Lujie; Ireland, Peter; Collins, Lance; Meng, Hui
2015-11-01
Particle-pair radial relative velocity (RV) in turbulence plays a critical role in droplet collision and cloud formation. Both simulations and experiments are performed to better understand RV of inertial particles in homogeneous and isotropic turbulence (HIT). However, past experimental measurement of particle RV statistics exhibited large deviations from DNS results (de Jong et al., 2010). In the current study, we identified intrinsic limitations in our previous study and devised a 4-frame particle tracking velocimetry technique to measure particle RV. In a second-generation, enclosed, fan-driven HIT chamber, both tracer and inertial particles were studied at R_ λ of 366. The experimentally measured RV statistics were compared with DNS with excellent agreement. Additionally, for both kinds of particles, the mean inward RV vs. particle separation distance r also matched very well with DNS, but at near-zero r, experimental values were slightly higher. To investigate the cause of this discrepancy, we compared DNS of both mono- and tri-dispersed particles. We found that the tri-dispersed particles exhibited higher mean inward RV at small r than any mono-dispersed particles. This suggests that the increase of mean inward RV in the experiment could be due to the Stokes number (St) distribution present in the particles, while DNS employed single St values. This work was supported by the National Science Foundation through a Collaborative Research Grant CBET-0967407.
NASA Astrophysics Data System (ADS)
Collins, Lance; Bragg, Andrew; Ireland, Peter
2014-11-01
In this talk, we consider the physical mechanism for the clustering of inertial particles in the inertial range of turbulence. By comparisons with DNS data we demonstrate that the mechanism in the theory of Zaichik et al. (Phys. Fluids 19, 113308, 2007) quantitatively describes the clustering of particles in the inertial range. We then analyze the theory for isotropic turbulence in the limit Reλ --> ∞ . For arbitrary St (Stokes number), there exists a separation in the inertial range beyond which Str << 1 , where Str is the Stokes number based on the eddy turnover timescale at separation r. The inertial-range clustering in this limit can be understood to be due to the preferential sampling of the coarse-grained velocity gradient tensor at that scale. At smaller separations, there may be transitions to Str ~ 1 , where a path history symmetry breaking effect dominates the clustering mechanism, and in some cases Str >> 1 , which implies ballistic behavior and a flat RDF. The scaling for each of these regimes is derived and compared to DNS, where applicable. Finally, we compare the results with the ``sweep-stick'' mechanism by Coleman and Vassilicos (Phys. Fluids 21, 113301, 2009) and discuss the similarities and differences between the two theories.
Long-time behavior of material-surface curvature in isotropic turbulence
NASA Technical Reports Server (NTRS)
Girimaji, S. S.
1992-01-01
The behavior at large times of the curvature of material elements in turbulence is investigated using Lagrangian velocity-gradient time series obtained from direct numerical simulations of isotropic turbulence. The main objectives are: to study the asymptotic behavior of the pdf curvature as a function of initial curvature and shape; and to establish whether the curvature of an initially plane material element goes to a stationary probability distribution. The evidence available in the literature about the asymptotic curvature-pdf of initially flat surfaces is ambiguous, and the conjecture is that it is quasi-stationary. In this work several material-element ensembles of different initial curvatures and shapes are studied. It is found that, at long times the moments of the logarithm of curvature are independent of the initial pdf of curvature. This, it is argued, supports the view that the curvature attains a stationary distribution at long times. It is also shown that, irrespective of initial shape or curvature, the shape of any material element at long times is cylindrical with a high probability.
The rotation and translation of non-spherical particles in homogeneous isotropic turbulence
NASA Astrophysics Data System (ADS)
Byron, Margaret
The motion of particles suspended in environmental turbulence is relevant to many scientific fields, from sediment transport to biological interactions to underwater robotics. At very small scales and simple shapes, we are able to completely mathematically describe the motion of inertial particles; however, the motion of large aspherical particles is significantly more complex, and current computational models are inadequate for large or highly-resolved domains. Therefore, we seek to experimentally investigate the coupling between freely suspended particles and ambient turbulence. A better understanding of this coupling will inform not only engineering and physics, but the interactions between small aquatic organisms and their environments. In the following pages, we explore the roles of shape and buoyancy on the motion of passive particles in turbulence, and allow these particles to serve as models for meso-scale aquatic organisms. We fabricate cylindrical and spheroidal particles and suspend them in homogeneous, isotropic turbulence that is generated via randomly-actuated jet arrays. The particles are fabricated with agarose hydrogel, which is refractive-index-matched to the surrounding fluid (water). Both the fluid and the particle are seeded with passive tracers, allowing us to perform Particle Image Velocimetry (PIV) simultaneously on the particle and fluid phase. To investigate the effects of shape, particles are fabricated at varying aspect ratios; to investigate the effects of buoyancy, particles are fabricated at varying specific gravities. Each particle type is freely suspended at a volume fraction of F=0.1%, for which four-way coupling interactions are negligible. The suspended particles are imaged together with the surrounding fluid and analyzed using stereoscopic PIV, which yields three velocity components in a two-dimensional measurement plane. Using image thresholding, the results are separated into simultaneous fluid-phase and solid-phase velocity
NASA Astrophysics Data System (ADS)
Onishi, Ryo; Vassilicos, J. C.
2014-11-01
This study investigates the collision statistics of inertial particles in inverse-cascading 2D homogeneous isotropic turbulence by means of a direct numerical simulation (DNS). A collision kernel model for particles with small Stokes number (St) in 2D flows is proposed based on the model of Saffman & Turner (1956) (ST56 model). The DNS results agree with this 2D version of the ST56 model for St < 0.1. It is then confirmed that our DNS results satisfy the 2D version of the spherical formulation of the collision kernel. The fact that the flatness factor stays around 3 in our 2D flow confirms that the present 2D turbulent flow is nearly intermittency-free. Collision statistics for St = 0.1, 0.4 and 0.6, i.e. for St <1, are obtained from the present 2D DNS and compared with those obtained from the three-dimensional (3D) DNS of Onishi et al. (2013). We have observed that the 3D radial distribution function at contact (g(R), the so-called clustering effect) decreases for St = 0.4 and 0.6 with increasing Reynolds number, while the 2D g(R) does not show a significant dependence on Reynolds number. This observation supports the view that the Reynolds-number dependence of g(R) observed in three dimensions is due to internal intermittency of the 3D turbulence. We have further investigated the local St, which is a function of the local flow strain rates, and proposed a plausible mechanism that can explain the Reynolds-number dependence of g(R).
Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence
NASA Astrophysics Data System (ADS)
Johnson, Perry L.; Meneveau, Charles
2015-08-01
One of the hallmarks of turbulent flows is the chaotic behavior of fluid particle paths with exponentially growing separation among them while their distance does not exceed the viscous range. The maximal (positive) Lyapunov exponent represents the average strength of the exponential growth rate, while fluctuations in the rate of growth are characterized by the finite-time Lyapunov exponents (FTLEs). In the last decade or so, the notion of Lagrangian coherent structures (which are often computed using FTLEs) has gained attention as a tool for visualizing coherent trajectory patterns in a flow and distinguishing regions of the flow with different mixing properties. A quantitative statistical characterization of FTLEs can be accomplished using the statistical theory of large deviations, based on the so-called Cramér function. To obtain the Cramér function from data, we use both the method based on measuring moments and measuring histograms and introduce a finite-size correction to the histogram-based method. We generalize the existing univariate formalism to the joint distributions of the two FTLEs needed to fully specify the Lyapunov spectrum in 3D flows. The joint Cramér function of turbulence is measured from two direct numerical simulation datasets of isotropic turbulence. Results are compared with joint statistics of FTLEs computed using only the symmetric part of the velocity gradient tensor, as well as with joint statistics of instantaneous strain-rate eigenvalues. When using only the strain contribution of the velocity gradient, the maximal FTLE nearly doubles in magnitude, highlighting the role of rotation in de-correlating the fluid deformations along particle paths. We also extend the large-deviation theory to study the statistics of the ratio of FTLEs. The most likely ratio of the FTLEs λ1 : λ2 : λ3 is shown to be about 4:1:-5, compared to about 8:3:-11 when using only the strain-rate tensor for calculating fluid volume deformations. The results
Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence
Johnson, Perry L. Meneveau, Charles
2015-08-15
One of the hallmarks of turbulent flows is the chaotic behavior of fluid particle paths with exponentially growing separation among them while their distance does not exceed the viscous range. The maximal (positive) Lyapunov exponent represents the average strength of the exponential growth rate, while fluctuations in the rate of growth are characterized by the finite-time Lyapunov exponents (FTLEs). In the last decade or so, the notion of Lagrangian coherent structures (which are often computed using FTLEs) has gained attention as a tool for visualizing coherent trajectory patterns in a flow and distinguishing regions of the flow with different mixing properties. A quantitative statistical characterization of FTLEs can be accomplished using the statistical theory of large deviations, based on the so-called Cramér function. To obtain the Cramér function from data, we use both the method based on measuring moments and measuring histograms and introduce a finite-size correction to the histogram-based method. We generalize the existing univariate formalism to the joint distributions of the two FTLEs needed to fully specify the Lyapunov spectrum in 3D flows. The joint Cramér function of turbulence is measured from two direct numerical simulation datasets of isotropic turbulence. Results are compared with joint statistics of FTLEs computed using only the symmetric part of the velocity gradient tensor, as well as with joint statistics of instantaneous strain-rate eigenvalues. When using only the strain contribution of the velocity gradient, the maximal FTLE nearly doubles in magnitude, highlighting the role of rotation in de-correlating the fluid deformations along particle paths. We also extend the large-deviation theory to study the statistics of the ratio of FTLEs. The most likely ratio of the FTLEs λ{sub 1} : λ{sub 2} : λ{sub 3} is shown to be about 4:1:−5, compared to about 8:3:−11 when using only the strain-rate tensor for calculating fluid volume
The length distribution of streamline segments in homogeneous isotropic decaying turbulence
NASA Astrophysics Data System (ADS)
Schaefer, P.; Gampert, M.; Peters, N.
2012-04-01
by Schaefer et al. ["Fast and slow changes of the length of gradient trajectories in homogenous shear turbulence," in Advances in Turbulence XII, edited by B. Eckhardt (Springer-Verlag, Berlin, 2009), pp. 565-572] we will refer to the morphological part of the evolution of streamline segments as slow changes while the topological part of the evolution is referred to as fast changes. This separation yields a transport equation for the probability density function (pdf) P(l) of the arclength l of streamline segments in which the slow changes translate into a convection and a diffusion term when terms up to second order are included and the fast changes yield integral terms. The overall temporal evolution (morphological and topological) of the arclength l of streamline segments is analyzed and associated with the motion of the above isosurface. This motion is diffusion controlled for small segments, while large segments are mainly subject to strain and pressure fluctuations. The convection velocity corresponds to the first order jump moment, while the diffusion term includes the second order jump moment. It is concluded, both theoretically and from direct numerical simulations (DNS) data of homogeneous isotropic decaying turbulence at two different Reynolds numbers, that the normalized first order jump moment is quasi-universal, while the second order one is proportional to the inverse of the square root of the Taylor based Reynolds number Re_{λ }^{-1/2}. Its inclusion thus represents a small correction in the limit of large Reynolds numbers. Numerical solutions of the pdf equation yield a good agreement with the pdf obtained from the DNS data. The interplay of viscous drift acting on small segments and linear strain acting on large segments yield, as it has already been concluded for dissipation elements, that the mean length of streamline segments should scale with Taylor microscale.
Kazansky, Peter G; Shimotsuma, Yasuhiko; Sakakura, Masaaki; Beresna, Martynas; Gecevičius, Mindaugas; Svirko, Yuri; Akturk, Selcuk; Qiu, Jianrong; Miura, Kiyotaka; Hirao, Kazuyuki
2011-10-10
We present the first experimental evidence of anisotropic photosensitivity of an isotropic homogeneous medium under uniform illumination. Our experiments reveal fundamentally new type of light induced anisotropy originated from the hidden asymmetry of pulsed light beam with a finite tilt of intensity front. We anticipate that the observed phenomenon, which enables employing mutual orientation of a light polarization plane and pulse front tilt to control interaction of matter with ultrashort light pulses, will open new opportunities in material processing. PMID:21997076
Magnetic Field Line Random Walk in Isotropic Turbulence with Zero Mean Field
NASA Astrophysics Data System (ADS)
Sonsrettee, W.; Subedi, P.; Ruffolo, D.; Matthaeus, W. H.; Snodin, A. P.; Wongpan, P.; Chuychai, P.
2015-01-01
In astrophysical plasmas, magnetic field lines often guide the motions of thermal and non-thermal particles. The field line random walk (FLRW) is typically considered to depend on the Kubo number R = (b/B 0)(l∥/l) for rms magnetic fluctuation b, large-scale mean field B 0, and parallel and perpendicular coherence scales l∥ and l, respectively. Here we examine the FLRW when R → ∞ by taking B 0 → 0 for finite bz (fluctuation component along B 0), which differs from the well-studied route with bz = 0 or bz Lt B 0 as the turbulence becomes quasi-two-dimensional (quasi-2D). Fluctuations with B 0 = 0 are typically isotropic, which serves as a reasonable model of interstellar turbulence. We use a non-perturbative analytic framework based on Corrsin's hypothesis to determine closed-form solutions for the asymptotic field line diffusion coefficient for three versions of the theory, which are directly related to the k -1 or k -2 moment of the power spectrum. We test these theories by performing computer simulations of the FLRW, obtaining the ratio of diffusion coefficients for two different parameterizations of a field line. Comparing this with theoretical ratios, the random ballistic decorrelation version of the theory agrees well with the simulations. All results exhibit an analog to Bohm diffusion. In the quasi-2D limit, previous works have shown that Corrsin-based theories deviate substantially from simulation results, but here we find that as B 0 → 0, they remain in reasonable agreement. We conclude that their applicability is limited not by large R, but rather by quasi-two-dimensionality.
MAGNETIC FIELD LINE RANDOM WALK IN ISOTROPIC TURBULENCE WITH ZERO MEAN FIELD
Sonsrettee, W.; Ruffolo, D.; Snodin, A. P.; Wongpan, P.; Subedi, P.; Matthaeus, W. H.; Chuychai, P. E-mail: david.ruf@mahidol.ac.th E-mail: pat.wongpan@postgrad.otago.ac.nz E-mail: prasub@udel.edu
2015-01-01
In astrophysical plasmas, magnetic field lines often guide the motions of thermal and non-thermal particles. The field line random walk (FLRW) is typically considered to depend on the Kubo number R = (b/B {sub 0})(ℓ{sub ∥}/ℓ ) for rms magnetic fluctuation b, large-scale mean field B {sub 0}, and parallel and perpendicular coherence scales ℓ{sub ∥} and ℓ , respectively. Here we examine the FLRW when R → ∞ by taking B {sub 0} → 0 for finite b{sub z} (fluctuation component along B {sub 0}), which differs from the well-studied route with b{sub z} = 0 or b{sub z} << B {sub 0} as the turbulence becomes quasi-two-dimensional (quasi-2D). Fluctuations with B {sub 0} = 0 are typically isotropic, which serves as a reasonable model of interstellar turbulence. We use a non-perturbative analytic framework based on Corrsin's hypothesis to determine closed-form solutions for the asymptotic field line diffusion coefficient for three versions of the theory, which are directly related to the k {sup –1} or k {sup –2} moment of the power spectrum. We test these theories by performing computer simulations of the FLRW, obtaining the ratio of diffusion coefficients for two different parameterizations of a field line. Comparing this with theoretical ratios, the random ballistic decorrelation version of the theory agrees well with the simulations. All results exhibit an analog to Bohm diffusion. In the quasi-2D limit, previous works have shown that Corrsin-based theories deviate substantially from simulation results, but here we find that as B {sub 0} → 0, they remain in reasonable agreement. We conclude that their applicability is limited not by large R, but rather by quasi-two-dimensionality.
A priori study of subgrid-scale flux of a passive scalar in isotropic homogeneous turbulence
Chumakov, Sergei
2008-01-01
We perform a direct numerical simulation (DNS) of forced homogeneous isotropic turbulence with a passive scalar that is forced by mean gradient. The DNS data are used to study the properties of subgrid-scale flux of a passive scalar in the framework of large eddy simulation (LES), such as alignment trends between the flux, resolved, and subgrid-scale flow structures. It is shown that the direction of the flux is strongly coupled with the subgrid-scale stress axes rather than the resolved flow quantities such as strain, vorticity, or scalar gradient. We derive an approximate transport equation for the subgrid-scale flux of a scalar and look at the relative importance of the terms in the transport equation. A particular form of LES tensor-viscosity model for the scalar flux is investigated, which includes the subgrid-scale stress. Effect of different models for the subgrid-scale stress on the model for the subgrid-scale flux is studied.
NASA Astrophysics Data System (ADS)
Miura, H.; Araki, K.
2014-07-01
Hall effects on local structures in homogeneous, isotropic, and incompressible magnetohydrodynamic turbulence are studied numerically. The transition of vortices from sheet-like to tubular structures induced by the Hall term is found, while the kinetic energy spectrum does not distinguish the two types of structures. It is shown by the use of the sharp low-pass filter that the transition occurs not only in the scales smaller than the ion skin depth but also in a larger scale. The transition is related with the forward energy transfer in the spectral space. Analyses by the use of the sharp low-pass filter show that the nonlinear energy transfer associated with the Hall term is dominated by the forward transfer and relatively local in the wave number space. A projection of the simulation data to a Smagorinsky-type sub-grid-scale model shows that the high wave number component of the Hall term may possibly be replaced by the model effectively.
Miura, H.; Araki, K.
2014-07-15
Hall effects on local structures in homogeneous, isotropic, and incompressible magnetohydrodynamic turbulence are studied numerically. The transition of vortices from sheet-like to tubular structures induced by the Hall term is found, while the kinetic energy spectrum does not distinguish the two types of structures. It is shown by the use of the sharp low-pass filter that the transition occurs not only in the scales smaller than the ion skin depth but also in a larger scale. The transition is related with the forward energy transfer in the spectral space. Analyses by the use of the sharp low-pass filter show that the nonlinear energy transfer associated with the Hall term is dominated by the forward transfer and relatively local in the wave number space. A projection of the simulation data to a Smagorinsky-type sub-grid-scale model shows that the high wave number component of the Hall term may possibly be replaced by the model effectively.
Clustering of vertically constrained passive particles in homogeneous and isotropic turbulence
NASA Astrophysics Data System (ADS)
van Hinsberg, Michel; de Pietro, Massimo; Biferale, Luca; Clercx, Herman; Toschi, Federico
2014-11-01
We analyze the dynamics of small particles confined within a horizontal fluid slab in a three-dimensional (3D) homogenous isotropic turbulent velocity field. Particles can freely move horizontally as fluid tracers but are vertically confined around a given horizontal plane via a simple linear restoring force. The present model may be considered as the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy or other mechanisms, are capable to maintain themselves in a shallow horizontal layer somewhere below the free surface of oceans or lakes. In the model varying the strength of the restoring force can control the thickness of the fluid slab in which the particles can move. Whenever some confinement is present, particle trajectories deviate from fluid tracers and experience an effectively compressible velocity field. We report a quantification of this effective compressibility as well as a quantification of preferential concentration of tracer particles in terms of the correlation dimension. We found that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale, that maximizes the clustering of the particles. This work is part of the research programmes 11PR2841 and FP112 of the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO). The work was partially funded by ERC Grant No 339032.
A Stochastic Model for the Relative Motion of High Stokes Number Particles in Isotropic Turbulence
NASA Astrophysics Data System (ADS)
Dhariwal, Rohit; Rani, Sarma; Koch, Donald
2014-11-01
In the current study, a novel analytical closure for the diffusion current in the PDF equation is presented that is applicable to high-inertia particle pairs with Stokes numbers Str >> 1 . Here Str is a Stokes number based on the time-scale τr of eddies whose size scales with pair separation r. Using this closure, Langevin equations were solved to evolve particle-pair relative velocities and separations in stationary isotropic turbulence. The Langevin equation approach enables the simulation of the full PDF of pair relative motion, instead of only the first few moments of the PDF as is the case in a moments-based approach. Accordingly, PDFs Ω (U | r) and Ω (Ur | r) are computed for various separations r, where the former is the PDF of relative velocity U and the latter is the PDF of the radial component of relative velocity Ur, both conditioned upon the separation r. Consistent with the DNS study of Sundaram & Collins, the Langevin simulations capture the transition of Ω (U | r) from being Gaussian at integral-scale separations to an exponential PDF at Kolmogorov-scale separations. The radial distribution functions (RDFs) computed from these simulations also show reasonable quantitative agreement with those from the DNS of Fevrier et al.
Power and Nonpower Laws of Passive Scalar Moments Convected by Isotropic Turbulence
NASA Astrophysics Data System (ADS)
Gotoh, Toshiyuki; Watanabe, Takeshi
2015-09-01
The scaling behavior of the moments of two passive scalars that are excited by two different methods and simultaneously convected by the same isotropic steady turbulence at Rλ=805 and Sc=0.72 is studied by using direct numerical simulation with N =40963 grid points. The passive scalar θ is excited by a random source that is Gaussian and white in time, and the passive scalar q is excited by the mean uniform scalar gradient. In the inertial convective range, the n th-order moments of the scalar increment δ θ (r ) do not obey a simple power law, but have the local scaling exponents ξnθ+βnlog (r /r*) with βn>0 . In contrast, the local scaling exponents of q have well-developed plateaus and saturate with increasing order. The power law of passive scalar moments is not trivial. The universality of passive scalars is found not in the moments, but in the normalized moments.
Power and nonpower laws of passive scalar moments convected by isotropic turbulence.
Gotoh, Toshiyuki; Watanabe, Takeshi
2015-09-11
The scaling behavior of the moments of two passive scalars that are excited by two different methods and simultaneously convected by the same isotropic steady turbulence at R_{λ}=805 and Sc=0.72 is studied by using direct numerical simulation with N=4096^{3} grid points. The passive scalar θ is excited by a random source that is Gaussian and white in time, and the passive scalar q is excited by the mean uniform scalar gradient. In the inertial convective range, the nth-order moments of the scalar increment δθ(r) do not obey a simple power law, but have the local scaling exponents ξ_{n}^{θ}+β_{n}log(r/r_{*}) with β_{n}>0. In contrast, the local scaling exponents of q have well-developed plateaus and saturate with increasing order. The power law of passive scalar moments is not trivial. The universality of passive scalars is found not in the moments, but in the normalized moments. PMID:26406833
NASA Astrophysics Data System (ADS)
Meneveau, Charles; Yang, Yunke; Perlman, Eric; Wan, Minpin; Burns, Randal; Szalay, Alex; Chen, Shiyi; Eyink, Gregory
2008-11-01
A public database system archiving a direct numerical simulation (DNS) data set of isotropic, forced turbulence is used for studying basic turbulence dynamics. The data set consists of the DNS output on 1024-cubed spatial points and 1024 time-samples spanning about one large-scale turn-over timescale. This complete space-time history of turbulence is accessible to users remotely through an interface that is based on the Web-services model (see http://turbulence.pha.jhu.edu). Users may write and execute analysis programs on their host computers, while the programs make subroutine-like calls that request desired parts of the data over the network. The architecture of the database is briefly explained, as are some of the new functions such as Lagrangian particle tracking and spatial box-filtering. These tools are used to evaluate and compare subgrid stresses and models.
Absorption of intense microwaves and ion acoustic turbulence due to heat transport
De Groot, J.S.; Liu, J.M.; Matte, J.P.
1994-02-04
Measurements and calculations of the inverse bremsstrahlung absorption of intense microwaves are presented. The isotropic component of the electron distribution becomes flat-topped in agreement with detailed Fokker-Planck calculations. The plasma heating is reduced due to the flat-topped distributions in agreement with calculations. The calculations show that the heat flux at high microwave powers is very large, q{sub max} {approx} 0.3 n{sub e}v{sub e}T{sub e}. A new particle model to, calculate the heat transport inhibition due to ion acoustic turbulence in ICF plasmas is also presented. One-dimensional PIC calculations of ion acoustic turbulence excited due to heat transport are presented. The 2-D PIC code is presently being used to perform calculations of heat flux inhibition due to ion acoustic turbulence.
Lipkens, B; Blackstock, D T
1998-09-01
A model experiment was reported to be successful in simulating the propagation of sonic booms through a turbulent atmosphere [B. Lipkens and D. T. Blackstock, J. Acoust. Soc. Am. 103, 148-158 (1998)]. In this study the effect on N wave characteristics of turbulence intensity and propagation distance through turbulence are investigated. The main parameters of interest are the rise time and the peak pressure. The effect of turbulence intensity and propagation distance is to flatten the rise time and peak pressure distributions. Rise time and peak pressure distributions always have positive skewness after propagation through turbulence. Average rise time grows with turbulence intensity and propagation distance. The scattering of rise time data is one-sided, i.e., rise times are almost always increased by turbulence. Average peak pressure decreases slowly with turbulence intensity and propagation distance. For the reported data a threefold increase in average rise time is observed and a maximum decrease of about 20% in average peak pressure. Rise times more than ten times that of the no turbulence value are observed. At most, the maximum peak pressure doubles after propagation through turbulence, and the minimum peak pressure values are about one-half the no-turbulence values. Rounded waveforms are always more common than peaked waveforms. PMID:9745733
NASA Astrophysics Data System (ADS)
Dou, Zhongwang; Pecenak, Zachary K.; Cao, Lujie; Woodward, Scott H.; Liang, Zach; Meng, Hui
2016-03-01
Enclosed flow apparatuses with negligible mean flow are emerging as alternatives to wind tunnels for laboratory studies of homogeneous and isotropic turbulence (HIT) with or without aerosol particles, especially in experimental validation of Direct Numerical Simulation (DNS). It is desired that these flow apparatuses generate HIT at high Taylor-microscale Reynolds numbers ({{R}λ} ) and enable accurate measurement of turbulence parameters including kinetic energy dissipation rate and thereby {{R}λ} . We have designed an enclosed, fan-driven, highly symmetric truncated-icosahedron ‘soccer ball’ airflow apparatus that enables particle imaging velocimetry (PIV) and other whole-field flow measurement techniques. To minimize gravity effect on inertial particles and improve isotropy, we chose fans instead of synthetic jets as flow actuators. We developed explicit relations between {{R}λ} and physical as well as operational parameters of enclosed HIT chambers. To experimentally characterize turbulence in this near-zero-mean flow chamber, we devised a new two-scale PIV approach utilizing two independent PIV systems to obtain both high resolution and large field of view. Velocity measurement results show that turbulence in the apparatus achieved high homogeneity and isotropy in a large central region (48 mm diameter) of the chamber. From PIV-measured velocity fields, we obtained turbulence dissipation rates and thereby {{R}λ} by using the second-order velocity structure function. A maximum {{R}λ} of 384 was achieved. Furthermore, experiments confirmed that the root mean square (RMS) velocity increases linearly with fan speed, and {{R}λ} increases with the square root of fan speed. Characterizing turbulence in such apparatus paves the way for further investigation of particle dynamics in particle-laden homogeneous and isotropic turbulence.
Analysis of subgrid models using direct and large-eddy simulations of isotropic turbulence
NASA Astrophysics Data System (ADS)
Menon, S.; Yeung, P. K.
1994-12-01
Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.
Flow structure interaction between a flexible cantilever beam and isotropic turbulence
NASA Astrophysics Data System (ADS)
Vogel, Andrew; Morvan, Thomas; Goushcha, Oleg; Andreopoulos, Yiannis; Elvin, Niell
2015-11-01
In the present experimental work we consider the degree of distortion of isotropy and homogeneity of grid turbulence caused by the presence of a thin flexible cantilever beam immersed in the flow aligned in the longitudinal direction. Beams of various rigidities and lengths were used in the experiments. Piezoelectric patches were attached to the beams which provided an output voltage proportional to the strain and therefore proportional to the beam's deflection. The experiments were carried out in a large scale wind tunnel and hot-wires were used to measure turbulence intensity in the vicinity of the beams for various values of the ratio of aerodynamic loading to beam's rigidity. It was found that the flow field distortion depends on the rigidity of the beam. For very rigid beams this distortion is of the order of the boundary layer thickness developing over the beam while for very flexible beams the distorted region is of the order of the beam's tip deflection. Analysis of the time-dependent signals indicated some correlation between the frequency of beam's vibration and flow structures detected. Supported by NSF Grant: CBET #1033117.
NASA Astrophysics Data System (ADS)
Aliseda, Alberto; Bourgoin, Mickael; Eswirp Collaboration
2014-11-01
We present preliminary results from a recent grid turbulence experiment conducted at the ONERA wind tunnel in Modane, France. The ESWIRP Collaboration was conceived to probe the smallest scales of a canonical turbulent flow with very high Reynolds numbers. To achieve this, the largest scales of the turbulence need to be extremely big so that, even with the large separation of scales, the smallest scales would be well above the spatial and temporal resolution of the instruments. The ONERA wind tunnel in Modane (8 m -diameter test section) was chosen as a limit of the biggest large scales achievable in a laboratory setting. A giant inflatable grid (M = 0.8 m) was conceived to induce slowly-decaying homogeneous isotropic turbulence in a large region of the test section, with minimal structural risk. An international team or researchers collected hot wire anemometry, ultrasound anemometry, resonant cantilever anemometry, fast pitot tube anemometry, cold wire thermometry and high-speed particle tracking data of this canonical turbulent flow. While analysis of this large database, which will become publicly available over the next 2 years, has only started, the Taylor-scale Reynolds number is estimated to be between 400 and 800, with Kolmogorov scales as large as a few mm . The ESWIRP Collaboration is formed by an international team of scientists to investigate experimentally the smallest scales of turbulence. It was funded by the European Union to take advantage of the largest wind tunnel in Europe for fundamental research.
Intensity fluctuations of ultrasonic scattering in a highly turbulent flow.
Shen, C; Lemmin, U
2000-05-01
Aspects of ultrasound intensity fluctuations backscattered from additive microstructures in a turbulent flow have been investigated theoretically and experimentally for the conditions of a small insonified volume, a high sound frequency and strong turbulence. These conditions are typically found in high resolution Doppler sonar applications. An easily applicable expression for the auto-correlation of scattering intensity fluctuations is obtained by introducing open-channel turbulence theory, a semi-empirical scalar spectrum (including a Batchelor spectrum) and a Gaussian window function. Experiments carried out in a laboratory-clear water, open-channel flow for different turbulence levels verify the underlying assumptions. A good agreement is found with the predictions made with the above-derived expression. The feasibility of extracting flow information from the backscattered intensity fluctuations is discussed. PMID:10857575
Intense sediment transport: Collisional to turbulent suspension
NASA Astrophysics Data System (ADS)
Berzi, Diego; Fraccarollo, Luigi
2016-02-01
A recent simple analytical approach to the problem of steady, uniform transport of sediment by a turbulent shearing fluid dominated by interparticle collisions is extended to the case in which the mean turbulent lift may partially or totally support the weight of the sediment. We treat the granular-fluid mixture as a continuum and make use of constitutive relations of kinetic theory of granular gases to model the particle phase and a simple mixing-length approach for the fluid. We focus on pressure-driven flows over horizontal, erodible beds and divide the flow itself into layers, each dominated by different physical mechanisms. This permits a crude analytical integration of the governing equations and to obtain analytical expressions for the distribution of particle concentration and velocity. The predictions of the theory are compared with existing laboratory measurements on the flow of glass spheres and sand particles in water. We also show how to build a regime map to distinguish between collisional, turbulent-collisional, and fully turbulent suspensions.
Gomez, T; Sagaut, P; Schilling, O; Zhou, Y
2006-07-05
A spectral subggrid-scale eddy viscosity and magnetic resisitivity model based on the eddy-damped quasi-normal Markovian (EDQNM) spectral kinetic and magnetic energy transfer presented in [12] is used in large-eddy simulation (LES) of large kinetic and magnetic Reynold number magneto-hydrodynamic (MHD) turbulence. The proposed model is assessed via a posteri tests on three-dimensional, incompressible, isotropic, non-helical, freely-decaying MHD turbulence at asymptotically large Reynolds numbers. Using LES with an initial condition characterized by an Alfv{acute e}n ratio of kinetic to magnetic energy {tau}{sub A} equal to unity, it is shown that at the kinetic energy spectrum E{sub K}(k) and magnetic energy spectrum E{sub M}(k) exhibit Kolmogorov -5/3 inertial subrange scalings in the LES, consistent with the EDQNM model.
NASA Astrophysics Data System (ADS)
Yue, Ping; Zhang, Qiang; Wang, Runyuan; Li, Yaohui; Wang, Sheng
2015-09-01
A deep understanding of turbulence structure is important for investigating the characteristics of the atmospheric boundary layer, especially over heterogeneous terrain. In the present study, turbulence intensity and turbulent kinetic energy (TKE) parameters are analyzed for different conditions with respect to stability, wind direction and wind speed over a valley region of the Loess Plateau of China during December 2003 and January 2004. The purpose of the study is to examine whether the observed turbulence intensity and TKE parameters satisfy Monin-Obukhov similarity theory (MOST), and analyze the wind shear effect on, and thermal buoyancy function of, the TKE, despite the terrain heterogeneity. The results demonstrate that the normalized intensity of turbulence follows MOST for all stability in the horizontal and vertical directions, as well as the normalized TKE in the horizontal direction. The shear effect of the wind speed in the Loess Plateau region is strong in winter and could enhance turbulence for all stability conditions. During daytime, the buoyancy and shear effect together constitute the generation of TKE under unstable conditions. At night, the contribution of buoyancy to TKE is relatively small, and mechanical shearing is the main production form of turbulence.
Power Law Decay in High Intensity Turbulence
NASA Astrophysics Data System (ADS)
Koster, Timothy; Puga, Alejandro; Nguyen, Baolong; Larue, John
2015-11-01
In the study reported herein, the region where the power decay law is applicable for active grid generated turbulence is found by an iterative approach which determines the largest range where the ratio of the dissipation from the power law and the dissipation from the temporal velocity derivative are unity. The square of the Taylor microscale, as noted by Batchelor (1953), is linearly related to downstream distance relative to the virtual origin and can be used in a straightforward manner to find the virtual origin. The fact that the decay of downstream velocity variance is described by a power law is shown to imply power law behavior for various other parameters such as the dissipation, the integral length scale, the Taylor microscale, the Kolmogorov microscale and the Taylor Reynolds number and that there is an algebraic relationship between the various power law exponents. Results are presented for various mean velocities to show the decay exponent as a function of the Taylor Reynolds number.
Comparison of hot wire/laser velocimeter turbulence intensity measurements
NASA Technical Reports Server (NTRS)
Meyers, J. F.; Wilkinson, S. P.
1982-01-01
The question of whether a random measure of particle velocities yields a good statistical estimate of the stationary condition of the turbulence flow field was investigated by comparing hot-wire and laser velocimeter turbulence intensity measurements. Great care was taken to insure that the instrument precision of both the laser velocimeter and hot wire was maximized. In this attempt to reduce the measurement uncertainties in the hot wire, direct digitization of the analog output signal was performed with point-by-point conversion to velocity through a spline fit calibration curve and the turbulence intensity function was calculated statistically. Frequent calibrations of the hot wire were performed using the laser velocimeter as the velocity standard to account for the presence of the small seed particles in the air flow and signal drift in the hot wire.
NASA Astrophysics Data System (ADS)
Joly, A.; Moulin, F.; Violeau, D.; Astruc, D.
2012-10-01
The prediction of solid bodies transport (such as algae, debris, sediment grains, or corrosion deposits) is a necessary requirement in many industrial or environmental processes. The physical processes involved cover a wide range of processes, from tidal flow to turbulent eddies and particle drag. A stochastic model was therefore developed to link the different scales of the physical processes where it was assumed that the particles are dilute enough that they do not affect the flow or the motion of other particles while being large enough that each particle does not follow exactly the fluid motions (i.e., macro-particles). The stochastic model is built in such a way that it uses Reynolds-averaged fluid properties to predict trajectories of individual particles. This model was then tested using experimental measurements obtained for isotropic particles released in semi-homogeneous turbulence. The turbulent flow was generated using a pair of oscillating grids and was characterized using particle image velocimetry measurements. The trajectories of the particles were measured using a pair of high resolution cameras. The comparison between the experimental data and different numerical models gives satisfactory results.
LATERAL TURBULENCE INTENSITY AND PLUME MEANDERING DURING STABLE CONDITIONS
There is much evidence in the literature for the presence of mesoscale lateral meanders in the stable nighttime boundary layer. These meanders result in relatively high lateral turbulence intensities and diffusion rates when averaged over an hour. Anemometer data from 17 overnigh...
NASA Astrophysics Data System (ADS)
Bassenne, Maxime; Urzay, Javier; Park, George I.; Moin, Parviz
2016-03-01
This study investigates control-based forcing methods for incompressible homogeneous-isotropic turbulence forced linearly in physical space which result in constant turbulent kinetic energy, constant turbulent dissipation (also constant enstrophy), or a combination of the two based on a least-squares error minimization. The methods consist of proportional controllers embedded in the forcing coefficients. During the transient, the controllers adjust the forcing coefficients such that the controlled quantity achieves very early a minimal relative error with respect to its target stationary value. Comparisons of these forcing methods are made with the non-controlled approaches of Rosales and Meneveau ["Linear forcing in numerical simulations of isotropic turbulence: Physical space implementations and convergence properties," Phys. Fluids 17, 095106 (2005)] and Carroll and Blanquart ["A proposed modification to Lundgren's physical space velocity forcing method for isotropic turbulence," Phys. Fluids 25, 105114 (2013)], using direct numerical simulations (DNS) and large-eddy simulations (LES). The results indicate that the proposed constant-energetics forcing methods shorten the transient period from a user-defined artificial flow field to Navier-Stokes turbulence while maintaining steadier statistics. Additionally, the proposed method of constant kinetic-energy forcing behaves more robustly in coarse LES when initial conditions are employed that favor the occurrence of subgrid-scale backscatter, whereas the other approaches fail to provide physical turbulent flow fields. For illustration, the proposed forcing methods are applied to dilute particle-laden homogeneous-isotropic turbulent flows; the results serve to highlight the influences of the forcing strategies on the disperse-phase statistics.
NASA Technical Reports Server (NTRS)
Bennett, Floyd V.; Yntema, Robert T.
1959-01-01
Several approximate procedures for calculating the bending-moment response of flexible airplanes to continuous isotropic turbulence are presented and evaluated. The modal methods (the mode-displacement and force-summation methods) and a matrix method (segmented-wing method) are considered. These approximate procedures are applied to a simplified airplane for which an exact solution to the equation of motion can be obtained. The simplified airplane consists of a uniform beam with a concentrated fuselage mass at the center. Airplane motions are limited to vertical rigid-body translation and symmetrical wing bending deflections. Output power spectra of wing bending moments based on the exact transfer-function solutions are used as a basis for the evaluation of the approximate methods. It is shown that the force-summation and the matrix methods give satisfactory accuracy and that the mode-displacement method gives unsatisfactory accuracy.
Turbulence generation through intense localized sources of energy
NASA Astrophysics Data System (ADS)
Maqui, Agustin; Donzis, Diego
2015-11-01
Mechanisms to generate turbulence in controlled conditions have been studied for nearly a century. Most common methods include passive and active grids with a focus on incompressible turbulence. However, little attention has been given to compressible flows, and even less to hypersonic flows, where phenomena such as thermal non-equilibrium can be present. Using intense energy from lasers, extreme molecule velocities can be generated from photo-dissociation. This creates strong localized changes in both the hydrodynamics and thermodynamics of the flow, which may perturb the flow in a way similar to an active grid to generate turbulence in hypersonic flows. A large database of direct numerical simulations (DNS) are used to study the feasibility of such an approach. An extensive analysis of single and two point statistics, as well as spectral dynamics is used to characterize the evolution of the flow towards realistic turbulence. Local measures of enstrophy and dissipation are studied to diagnose the main mechanisms for energy exchange. As commonly done in compressible flows, dilatational and solenoidal components are separated to understand the effect of acoustics on the development of turbulence. Further results for cases that assimilate laboratory conditions will be discussed. The authors gratefully acknowledge the support of AFOSR.
NASA Astrophysics Data System (ADS)
Tóth-Katona, Tibor; Gleeson, James T.
2004-01-01
Fluctuations of the injected electric power during electroconvection (EHC) of liquid crystals are reported in both the conductive and the dielectric regime of convection. The amplitude and the frequency of the fluctuations, as well as the probability density functions have been compared in these two regimes and substantial differences have been found both in defect turbulence of EHC and at the DSM1→DSM2 transition.
NASA Astrophysics Data System (ADS)
Lee, J. H.; Kevin; Monty, J. P.; Hutchins, N.
2016-08-01
The discrepancy between measured turbulence intensity obtained from experiments in wall-bounded turbulence and the fully resolved reference results (usually from DNS datasets) are often attributed to spatial resolution issues, especially in PIV measurements due to the presence of spatial averaging within the interrogation region/volume. In many cases, in particular at high Reynolds numbers (where there is a lack of DNS data), there is no attempt to verify that this is the case. There is a risk that attributing unexpected PIV statistics to spatial resolution, without careful checks, could mask wider problems with the experimental setup or test facility. Here, we propose a robust technique to validate the under-resolved PIV obtained turbulence intensity profiles for canonical wall-bounded turbulence. This validation scheme is independent of Reynolds number and does not rely on empirical functions. It is based on arguments that (1) the viscous-scaled small-scale turbulence energy is invariant with Reynolds number and that (2) the spatially under-resolved measurement is sufficient to capture the large-scale energy. This then suggests that we can estimate the missing energy from volume-filtered DNS data at much lower Reynolds numbers. Good agreement is found between the experimental results and estimation profiles for all three velocity components, demonstrating that the estimation tool successfully computes the missing energy for given spatial resolutions over a wide range of Reynolds numbers. A database for a canonical turbulent boundary layer and associated MATLAB function are provided that enable this missing energy to be calculated across a range of interrogation volume sizes, so that users do not require access to raw DNS data. This methodology and tool will provide PIV practitioners, investigating canonical wall-bounded turbulent flow with a convenient check of the effects of spatial resolution on a given experiment.
NASA Astrophysics Data System (ADS)
Coleman, S. W.; Vassilicos, J. C.
2009-11-01
Our work focuses on the sweep-stick mechanism of particle clustering in turbulent flows introduced by Chen et al. [L. Chen, S. Goto, and J. C. Vassilicos, "Turbulent clustering of stagnation points and inertial particles," J. Fluid Mech. 553, 143 (2006)] for two-dimensional (2D) inverse cascading homogeneous, isotropic turbulence (HIT), whereby heavy particles cluster in a way that mimics the clustering of zero-acceleration points. We extend this phenomenology to three-dimensional (3D) HIT, where it was previously reported that zero-acceleration points were extremely rare. Having obtained a unified mechanism we quantify the Stokes number dependency of the probability of the heavy particles to be at zero-acceleration points and show that in the inertial range of Stokes numbers, the sweep-stick mechanism is dominant over the conventionally proposed mechanism of heavy particles being centrifuged from high vorticity regions to high strain regions. Finally, having a clustering coincidence between particles and zero-acceleration points, both in 2D and 3D HIT, motivates us to demonstrate the sweep and stick parts of the mechanism in both dimensions. The sweeping of regions of low acceleration regions by the local fluid velocity in both flows is demonstrated by introducing a velocity of the acceleration field. Finally, the stick part is demonstrated by showing that heavy particles statistically move with the same velocity as zero-acceleration points, while moving away from any nonzero-acceleration region, irrespective of their Stokes number. These results explain the clustering of inertial particles given the clustering of zero-acceleration points.
Estimating Overwater Turbulence Intensity from Routine Gust-Factor Measurements.
NASA Astrophysics Data System (ADS)
Hsu, S. A.; Blanchard, Brian W.
2004-12-01
For overwater diffusion estimates the Offshore and Coastal Dispersion (OCD) model is preferred by the U.S. Environmental Protection Agency. The U.S. Minerals Management Service has recommended that the OCD model be used for emissions located on the outer continental shelf. During southerly winds over the Gulf of Mexico, for example, the pollutants from hundreds of offshore platforms may affect the gulf coasts. In the OCD model, the overwater plume is described by the Gaussian equation, which requires the computation of σy and σz, which are, in turn, related to the turbulence intensity, overwater trajectory, and atmospheric stability. On the basis of several air sea interaction experiments [the Barbados Oceanographic and Meteorological Experiment (BOMEX), the Air-Mass Transformation Experiment (AMTEX), and, most recently, the Tropical Ocean and Global Atmosphere Coupled Ocean Atmosphere Response Experiment (TOGA COARE)] and the extensive datasets from the National Data Buoy Center (NDBC), it is shown that under neutral and stable conditions the overwater turbulence intensities are linearly proportional to the gust factor (G), which is the ratio of the wind gust and mean wind speed at height z (Uz) as reported hourly by the NDBC buoys. Under unstable conditions, it is first shown that the popular formula relating the horizontal turbulence intensity (σu,/u, where u is the friction velocity) to the ratio of the mixing height (h) and the buoyancy length (L) (i.e., h/L) suffers from a self-correlation problem and cannot be used in the marine environment. Then, alternative formulas to estimate the horizontal turbulence intensities (σu,/Uz) using G are proposed for practical applications. Furthermore, formulas to estimate u and z/L are fundamentally needed in air sea interaction studies, in addition to dispersion meteorology.
NASA Astrophysics Data System (ADS)
Donzis, Diego; Yakhot, Victor; Sreenivasan, K. R.
2015-11-01
Most approaches to understand turbulence have sought universal behavior believed to manifest at high Reynolds numbers (Rλ). However, recent theory and simulations suggest that universal characteristics, such as the non-trivial anomalous scaling exponents of moments of velocity gradients, emerge even at very low Rλ at which no inertial range exists. Furthermore, with decreasing Reynolds numbers, a transition occurs from fully intermittent turbulence to (approximately) Gaussian behavior at an apparently universal critical Rλ. A potential implication of these observations is that significant information concerning the inertial range (e.g. scaling exponents) is already manifest in the dissipation range at very low Rλ. Thus, high Rλ properties can be studied with well-resolved low-Rλ simulations instead of marginally resolved high-Reynolds flows. The focus of this talk is to explore signatures of universality at high-Reynolds numbers in the dissipation range of highly resolved DNS (kmax η ~ O (20)) for Rλ up to 90, and decaying simulations close to the critical Rλ. In addition to statistics of velocity gradients and dissipation we explore evidence of Beltramization as suggested in past theoretical work.
NASA Astrophysics Data System (ADS)
Ireland, Peter J.; Bragg, Andrew D.; Collins, Lance R.
2016-06-01
In Part I of this study, we analyzed the motion of inertial particles in isotropic turbulence in the absence of gravity using direct numerical simulation (DNS). Here, in Part II, we introduce gravity and study its effect over a wide range of flow Reynolds numbers, Froude numbers, and particle Stokes numbers. We see that gravity causes particles to sample the flow more uniformly and reduces the time particles can spend interacting with the underlying turbulence. We also find that gravity tends to increase inertial particle accelerations, and we introduce a model to explain that effect. We then analyze the particle relative velocities and radial distribution functions (RDFs), which are generally seen to be independent of Reynolds number for low and moderate Kolmogorov-scale Stokes numbers $St$. We see that gravity causes particle relative velocities to decrease, and that the relative velocities have higher scaling exponents with gravity. We observe that gravity has a non-trivial effect on clustering, acting to decrease clustering at low $St$ and to increase clustering at high $St$. By considering the effect of gravity on the clustering mechanisms described in the theory of Zaichik & Alipchenkov (New J. Phys., 11:103018, 2009), we provide an explanation for this non-trivial effect of gravity. We also show that when the effects of gravity are accounted for in the theory of Zaichik & Alipchenkov, the results compare favorably with DNS. The relative velocities and RDFs exhibit considerable anisotropy at small separations, and this anisotropy is quantified using spherical harmonic functions. We use the relative velocities and the RDFs to compute the particle collision kernels, and find that the collision kernel remains as it was for the case without gravity, namely nearly independent of Reynolds number for low and moderate $St$.
Two-color correlation between intensity fluctuations in atmospheric turbulence
NASA Astrophysics Data System (ADS)
Luo, Meilan; Zhao, Daomu
2016-06-01
The correlation between intensity fluctuations generated by two varying wavelengths through a turbulent medium is investigated, where the influences arising from source correlation and perturbation of atmosphere are mainly emphasized. It is demonstrated that the correlation between intensity fluctuations can be enhanced or reduced by modulating the difference of two incident wavelengths. For shorter wavelength, the correlation between intensity fluctuations is stronger at the far field. In addition, in the case of single wavelength, a relationship λ1z1 =λ2z2 =λnzn holding in free space could be found, from which the distance where the peak value occurs may be inferred. However, it can be destroyed by increasing the strength of atmosphere.
Rotational motion of elongated particles in isotropic turbulent flow: statistical perspective
NASA Astrophysics Data System (ADS)
Zhao, Lihao; Andersson, Helge; Variano, Evan
2014-11-01
We consider the rotational motion of non-spherical particles in turbulent flow, comparing the statistics of particles' angular velocity to the corresponding quantities computed in the fluid phase. We use numerical (DNS) and laboratory measurements for particles that are both larger and smaller than the Kolmogorov lengthscale. The particles are spheroids or rods, with aspect ratios between 1 and 10. We will discuss the subtleties of defining a meaningful Stokes number for these particles, focusing on the effect of asphericity and the fact that our interest is in rotation and not translation. Comparing the probability density function of angular velocity between fluid and particle phase indicates that the angular velocity of particles has a narrower distribution than that of the fluid phase, and that. particles do respond to extreme events in the fluid phase. The first four moments of the PDFs are analyzed, and these show that the ``filtering'' effect is very similar between DNS and lab experiments, despite differences in particle sizes and mass. We propose a nondimensional curve for predicting the magnitude of the filtering effect, and discuss the implications of this curve for the definition of Stokes number, as discussed earlier. This work has been supported by grants from the Peder Sather Center for Advanced Study at UC Berkeley and from the Research Council of Norway (Contract No. 213917/F20).
NASA Astrophysics Data System (ADS)
Jin, Guodong; He, Guo-Wei
2015-11-01
Clustering and intermittency in radial relative velocity (RRV) of heavy particles of same size settling in turbulent flows can be remarkably changed due to gravity. Clustering is monotonically reduced at Stokes number less than 1 under gravity due to the disability of the centrifugal mechanism, however it is non-monotonically enhanced at Stokes number greater than 1 due to the multiplicative amplification in the case that the proposed effective Kubo number is less than 1. Although gravity causes monotonical reduction in the rms of RRV of particles at a given Stokes number with decreasing Froude number, the variation tendency in the tails of standardized PDF of RRV versus Froude number is obviously different: the tails become narrower at a small Stokes number, while they become broader at a large Stokes number. The mechanism of this variation stems from the compromise between the following two competing factors. The mitigation of correlation of particle positions and the regions of high strain rate which are more intermittent reduces the intermittency in RRV at small Stokes numbers, while the significant reduction in the backward-in-time relative separations will make particle pairs see small-scale structures, leading to a higher intermittency in RRV at large Stokes numbers. NSAF of China (grant number U1230126); NSFC (grant numbers 11072247 and 11232011).
The structure of partially premixed methane flames in high-intensity turbulent flows
Yaldizli, Murat; Mehravaran, Kian; Mohammad, Hyderuddin; Jaberi, Farhad A.
2008-09-15
Direct numerical simulations (DNS) are conducted to study the structure of partially premixed and non-premixed methane flames in high-intensity two-dimensional isotropic turbulent flows. The results obtained via ''flame normal analysis'' show local extinction and reignition for both non-premixed and partially premixed flames. Dynamical analysis of the flame with a Lagrangian method indicates that the time integrated strain rate characterizes the finite-rate chemistry effects and the flame extinction better than the strain rate. It is observed that the flame behavior is affected by the ''pressure-dilatation'' and ''viscous-dissipation'' in addition to strain rate. Consistent with previous studies, high vorticity values are detected close to the reaction zone, where the vorticity generation by the ''baroclinic torque'' was found to be significant. The influences of (initial) Reynolds and Damkoehler numbers, and various air-fuel premixing levels on flame and turbulence variables are also studied. It is observed that the flame extinction occurs similarly in flames with different fuel-air premixing. Our simulations also indicate that the CO emission increases as the partial premixing of the fuel with air increases. Higher values of the temperature, the OH mass fraction and the CO mass fraction are observed within the flame zone at higher Reynolds numbers. (author)
NASA Astrophysics Data System (ADS)
Ireland, Peter J.; Bragg, Andrew D.; Collins, Lance R.
2016-06-01
In this study, we analyze the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number $R_\\lambda$ on the particle statistics is examined over the largest range to date (from $R_\\lambda = 88-597$). We first explore the effect of preferential sampling on the single-particle statistics, and use our understanding of preferential sampling to provide a physical explanation for many of the trends in the particle velocity gradients, kinetic energies, and accelerations at low $St$. As $St$ increases, inertial filtering effects become more important, causing the particle kinetic energies and accelerations to decrease. We then consider particle-pair statistics, and focus our attention on the relative velocities and radial distribution functions (RDFs) of the particles. The relative velocity statistics indicate that preferential-sampling effects are important for $St \\lesssim 0.1$ and that path-history/non-local effects become increasingly important for $St \\gtrsim 0.2$. The lower-order relative velocity statistics are only weakly sensitive to changes in Reynolds number at low $St$. We find that the RDFs peak near $St$ of order unity, that they exhibit power-law scaling for low and intermediate $St$, and that they are largely independent of Reynolds number for low and intermediate $St$. We also observe that at large $St$, changes in the RDF are related to changes the scaling exponents of the relative velocity variances. The particle collision kernel is found to be largely insensitive to the flow Reynolds number, suggesting that relatively low-Reynolds-number simulations may be able to capture much of the relevant physics of droplet collisions and growth in the adiabatic cores of atmospheric clouds.
Pawar, Shashikant S; Arakeri, Jaywant H
2016-08-01
Frequency spectra obtained from the measurements of light intensity and angle of arrival (AOA) of parallel laser light propagating through the axially homogeneous, axisymmetric buoyancy-driven turbulent flow at high Rayleigh numbers in a long (length-to-diameter ratio of about 10) vertical tube are reported. The flow is driven by an unstable density difference created across the tube ends using brine and fresh water. The highest Rayleigh number is about 8×10^{9}. The aim of the present work is to find whether the conventional Obukhov-Corrsin scaling or Bolgiano-Obukhov (BO) scaling is obtained for the intensity and AOA spectra in the case of light propagation in a buoyancy-driven turbulent medium. Theoretical relations for the frequency spectra of log amplitude and AOA fluctuations developed for homogeneous isotropic turbulent media are modified for the buoyancy-driven flow in the present case to obtain the asymptotic scalings for the high and low frequency ranges. For low frequencies, the spectra of intensity and vertical AOA fluctuations obtained from measurements follow BO scaling, while scaling for the spectra of horizontal AOA fluctuations shows a small departure from BO scaling. PMID:27505375
NASA Astrophysics Data System (ADS)
Sahli, Ouissem; Adjlout, Lahouari; Ladjedel, Omar; Amine Ghazi, Mohamed
2016-03-01
The present work is an experimental investigation on the effect of the turbulence intensity variation in a staggered tube bundle equipped with grooves at 90° and 270°.The experiments were carried out in a subsonic wind tunnel. Three Reynolds numbers and three turbulence levels were tested. The pressure distributions and drag forces were measured. Surface visualizations were also performed. The obtained results show that the turbulence intensity for different Reynolds number has an influence on the reduction of the drag coefficient.
NASA Astrophysics Data System (ADS)
Kestoras, M. D.; Simon, T. W.
1995-04-01
In an attempt to characterize the turbulence characteristics of the high free-stream turbulence (TI approximately 8%) flow, several experiments are conducted on a flat recovery wall downstream of sustained concave curvature in the presence of such flow. A turbulent boundary layer that grows from the leading edge of a concave wall, then passes onto a downstream flat wall is considered. The results indicate that turbulence intensities increase profoundly in the outer region of the boundary layer over the recovery wall. The recovery wall is found to be lifted off by near-wall turbulent eddies while a 'stabilized' region forms near the wall. Contrary to the low-free-stream turbulence intensity flow, turbulent eddies penetrate the outer parts of the 'stabilized' region. The behavior of the Stanton numbers as well as the velocity distribution on the core of the flow are also accounted for.
NASA Astrophysics Data System (ADS)
Kulikov, Victor A.
2016-03-01
The inner scale and the refractive-index structure parameter Cn2 are estimated by the intensity scintillations of the laser beam propagated through the turbulent water layer. The laser beam scintillations caused by underwater propagation can be described in the framework of a turbulence model with accounting of the inner scale. An intensity field has been observed at a 2 m distance when the laser beam passed through a water convective cell with a characteristic Rayleigh number of about 108. Similar intensity fields were simulated by using the split-step method. Characteristics of the experimentally obtained and numerically simulated intensity fields are compared and analyzed. A simple method of the turbulent parameter estimation is proposed.
Intensity of thunderstorm-generated turbulence revealed by large-eddy simulation
NASA Astrophysics Data System (ADS)
Lane, Todd P.; Sharman, Robert D.
2014-03-01
Thunderstorms are characterized by turbulent processes that constitute an important aviation hazard and cause vertical transport of atmospheric constituents. Turbulence occurs within cloud and in the surrounding clear air, but, despite its importance, the characteristics of thunderstorm-generated turbulence and its spatial distribution are poorly understood, especially outside of cloud. Here we use large-eddy simulation to characterize turbulence generated by a canonical thunderstorm. The simulation identifies regions of notable three-dimensional anisotropic turbulence more than 5 km above the storm, in a shallow layer above the storm's anvil, and a horizontally asymmetric pattern of weaker turbulence that extends more than 50 km horizontally away from the cloud. Our results provide the first continuous estimate of turbulence intensity in and around thunderstorms and represent a major step toward improved turbulence avoidance methods. The results have broader implications for understanding the fundamental aspects of how thunderstorms affect their environment through vertical exchange processes.
Fichtl, G.H.
1983-09-01
When designing a wind energy converison system (WECS), it may be necessary to take into account the distribution of wind across the disc of rotation. The specific engineering applications include structural strength, fatigue, and control. This wind distribution consists of two parts, namely that associated with the mean wind profile and that associated with the turbulence velocity fluctuation field. The work reported herein is aimed at the latter, namely the distribution of turbulence velocity fluctuations across the WECS disk of rotation. A theory is developed for the two-time covariance matrix for turbulence velocity vector components for wind energy conversion system (WECS) design. The theory is developed for homogeneous and iotropic turbulance with the assumption that Taylor's hypothesis is valid. The Eulerian turbulence velocity vector field is expanded about the hub of the WECS. Formulae are developed for the turbulence velocity vector component covariance matrix following the WECS blade elements. It is shown that upon specification of the turbulence energy spectrum function and the WECS rotation rate, the two-point, two-time covariance matrix of the turbulent flow relative to the WECS bladed elements is determined. This covariance matrix is represented as the sum of nonstationary and stationary contributions. Generalized power spectral methods are used to obtain two-point, double frequency power spectral density functions for the turbulent flow following the blade elements. The Dryden turbulence model is used to demonstrate the theory. A discussion of linear system response analysis is provided to show how the double frequency turbulence spectra might be used to calculate response spectra of a WECS to turbulent flow. Finally the spectrum of the component of turbulence normal to the WECS disc of rotation, following the blade elements, is compared with experimental results.
How is turbulence intensity determined by macroscopic variables in a toroidal plasma?
NASA Astrophysics Data System (ADS)
Inagaki, S.; Tokuzawa, T.; Tamura, N.; Itoh, S.-I.; Kobayashi, T.; Ida, K.; Shimozuma, T.; Kubo, S.; Tanaka, K.; Ido, T.; Shimizu, A.; Tsuchiya, H.; Kasuya, N.; Nagayama, Y.; Kawahata, K.; Sudo, S.; Yamada, H.; Fujisawa, A.; Itoh, K.; the LHD Experiment Group
2013-11-01
We report observations of the dynamic response of micro-fluctuations and turbulent flux to a low-frequency heating power modulation in the Large Helical Device. The responses of heat flux and micro-fluctuation intensity differ from that of the change in temperature gradient. This result violates the local transport model, where turbulence is determined by the local temperature gradient. A new relationship between flux, gradient and turbulence is found. In addition to the temperature gradient, the heating rate is proposed as a new, direct controlling parameter of turbulence to explain the fast response of turbulence against periodic modulation of heating power.
NASA Technical Reports Server (NTRS)
Arnold, S. M.; Binienda, W. K.; Tan, H. Q.; Xu, M. H.
1992-01-01
Analytical derivations of stress intensity factors (SIF's) of a multicracked plate can be complex and tedious. Recent advances, however, in intelligent application of symbolic computation can overcome these difficulties and provide the means to rigorously and efficiently analyze this class of problems. Here, the symbolic algorithm required to implement the methodology described in Part 1 is presented. The special problem-oriented symbolic functions to derive the fundamental kernels are described, and the associated automatically generated FORTRAN subroutines are given. As a result, a symbolic/FORTRAN package named SYMFRAC, capable of providing accurate SIF's at each crack tip, was developed and validated. Simple illustrative examples using SYMFRAC show the potential of the present approach for predicting the macrocrack propagation path due to existing microcracks in the vicinity of a macrocrack tip, when the influence of the microcrack's location, orientation, size, and interaction are taken into account.
NASA Technical Reports Server (NTRS)
Binienda, W. K.; Arnold, S. M.; Tan, H. Q.
1992-01-01
An essential part of describing the damage state and predicting the damage growth in a multicracked plate is the accurate calculation of stress intensity factors (SIF's). Here, a methodology and rigorous solution formulation for SIF's of a multicracked plate, with fully interacting cracks, subjected to a far-field arbitrary stress state is presented. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. This analytical derivation and numerical solution are obtained by using intelligent application of symbolic computations and automatic FORTRAN generation capabilities (described in the second part of this paper). As a result, a symbolic/FORTRAN package, named SYMFRAC, that is capable of providing accurate SIF's at each crack tip was developed and validated.
NASA Technical Reports Server (NTRS)
Boldman, D. R.; Brinich, P. F.
1977-01-01
To explain the increase in noise when a jet impinges on a large flat plate, mean velocity, turbulence intensity, and scale were measured at nominal nozzle-exit velocities of 61, 138, and 192 meters per second with the plate located 7.1 nozzle-exit diameters from the nozzle. The maximum turbulence intensities in free and impinging jets were about the same; however, the integral length scale near the plate surface was only about one-half the free jet scale. The measured intensities and length scales, in conjunction with a contemporary theory of aerodynamic noise, provided a good explanation for the observed increase in noise associated with the impinging jet. An increase in the volume of highly turbulent flow could be the principal reason for the increase in noise.
Yi, S.; Kwon, J. M.; Rhee, T.; Diamond, P. H.
2012-11-15
This paper reports the results of gyrokinetic simulation studies of ion temperature gradient driven turbulence which investigate the role of non-resonant modes in turbulence spreading, turbulence regulation, and self-generated plasma rotation. Non-resonant modes, which are those without a rational surface within the simulation domain, are identified as nonlinearly driven, radially extended convective cells. Even though the amplitudes of such convective cells are much smaller than that of the resonant, localized turbulence eddies, we find from bicoherence analysis that the mode-mode interactions in the presence of such convective cells increase the efficiency of turbulence spreading associated with nonlocality phenomena. Artificial suppression of the convective cells shows that turbulence spreading is reduced, and that the turbulence intensity profile is more localized. The more localized turbulence intensity profile produces stronger Reynolds stress and E Multiplication-Sign B shear flows, which in turn results in more effective turbulence self-regulation. This suggests that models without non-resonant modes may significantly underestimate turbulent fluctuation levels and transport.
NASA Astrophysics Data System (ADS)
Furukawa, Junichi; Noguchi, Yoshiki; Hirano, Toshisuke; Williams, Forman A.
2002-07-01
The density change across premixed flames propagating in turbulent flows modifies the turbulence. The nature of that modification depends on the regime of turbulent combustion, the burner design, the orientation of the turbulent flame and the position within the flame. The present study addresses statistically stationary turbulent combustion in the flame-sheet regime, in which the laminar-flame thickness is less than the Kolmogorov scale, for flames stabilized on a vertically oriented cylindrical burner having fully developed upward turbulent pipe flow upstream from the exit. Under these conditions, rapidly moving wrinkled laminar flamelets form the axisymmetric turbulent flame brush that is attached to the burner exit. Predictions have been made of changes in turbulence properties across laminar flamelets in such situations, but very few measurements have been performed to test the predictions. The present work measures individual velocity changes and changes in turbulence across flamelets at different positions in the turbulent flame brush for three different equivalence ratios, for comparison with theory.
Kestoras, M.D.; Simon, T.W.
1995-04-01
Experiments are conducted on a flat recovery wall downstream of sustained concave curvature in the presence of high free-stream turbulence (TI {approximately} 8%). This flow simulates some of the features of the flow on the latter parts of the pressure surface of a gas turbine airfoil. The combined effects of concave curvature and TI, both present in the flow over a turbine airfoil, have so far been little studied. Computation of such flows with standard turbulence closure models has not been particularly successful. This experiment attempts to characterize the turbulence characteristics of this flow. In the present study, a turbulent boundary layer grows from the leading edge of a concave wall, then passes onto a downstream flat wall. Results show that turbulence intensities increase profoundly in the outer region of the boundary layer over the recovery wall. Near-wall turbulent eddies appear to lift off the recovery wall and a stabilized region forms near the wall. In contrast to a low-free-stream turbulence intensity flow, turbulent eddies penetrate the outer parts of the stabilized region where sharp velocity and temperature gradients exist. These eddies can more readily transfer momentum and heat. As a result, skin friction coefficients and Stanton numbers on the recovery wall are 20 and 10%, respectively, above their values in the low-free-stream turbulence intensity case. Stanton numbers do not undershoot flat-wall expectations at the same Re{sub {Delta}2} values as seen in the low-TI case. Remarkably, the velocity distribution in the core of the flow over the recovery wall exhibits a negative gradient normal to the wall under high-free-stream turbulence intensity conditions. This velocity distribution appears to be the result of two effects: (1) cross transport of kinetic energy by boundary work in the upstream curved flow and (2) readjustment of static pressure profiles in response to the removal of concave curvature.
Gaussian Multiplicative Chaos for Symmetric Isotropic Matrices
NASA Astrophysics Data System (ADS)
Chevillard, Laurent; Rhodes, Rémi; Vargas, Vincent
2013-02-01
Motivated by isotropic fully developed turbulence, we define a theory of symmetric matrix valued isotropic Gaussian multiplicative chaos. Our construction extends the scalar theory developed by J.P. Kahane in 1985.
Measurements in film cooling flows: Hole L/D and turbulence intensity effects
Burd, S.W.; Kaszeta, R.W.; Simon, T.W.
1996-12-31
Hot-wire anemometry of simulated film cooling was used to study the influence of freestream turbulence intensity and film cooling hole length-to-diameter ratio on mean velocity and turbulence intensity. Measurements were made in the zone where the coolant and freestream flows mix. Flow from one row of film cooling holes with a streamwise injection of 35{degree} and no lateral injection and with a coolant- to-freestream flow velocity ratio of 1.0 was investigated under freestream turbulence levels of 0.5 and 12%. Coolant-to-freestream density ratio was unity. Two length-to-diameter ratios for the film cooling holes, 2.3 and 7.0, are tested. Results show that under low freestream turbulence conditions, pronounced differences exist in the flowfield between L/D=7.0 and 2.3; the differences are less prominent at high freestream turbulence intensities. Generally, short-L/D injection results in ``jetting`` of the coolant further into the freestream flow and enhanced mixing. Other changes in the flowfield attributable to a rise in freestream turbulence intensity to engine- representative conditions are documented. 15 figs, 2 tabs, refs.
Free-Stream Turbulence Intensity in the Langley 14- by 22-Foot Subsonic Tunnel
NASA Technical Reports Server (NTRS)
Neuhart, Dan H.; McGinley, Catherine B.
2004-01-01
An investigation was conducted using hot-wire anemometry to determine the turbulence intensity levels in the test section of the Langley 14- by 22-Foot Subsonic Tunnel in the closed or walls-down configuration. This study was one component of the three-dimensional High-Lift Flow Physics experiment designed to provide code validation data. Turbulence intensities were measured during two stages of the study. In the first stage, the free-stream turbulence levels were measured before and after a change was made to the floor suction surface of the wind tunnel s boundary layer removal system. The results indicated that the new suction surface at the entrance to the test section had little impact on the turbulence intensities. The second stage was an overall flow quality survey of the empty tunnel including measurements of the turbulence levels at several vertical and streamwise locations. Results indicated that the turbulence intensity is a function of tunnel dynamic pressure and the location in the test section. The general shape of the frequency spectrum is fairly consistent throughout the wind tunnel, changing mostly in amplitude (also slightly with frequency) with change in condition and location.
Probe systems for measuring static pressure and turbulence intensity in fluid streams
NASA Technical Reports Server (NTRS)
Rossow, Vernon J. (Inventor)
1993-01-01
A method and an apparatus for measuring time-averaged static or ambient pressure and turbulence intensity in a turbulent stream are discussed. The procedure involves placing a plurality of probes in the stream. Each probe responds in a different manner to characteristics of the fluid stream, preferably as a result of having varying cross sections. The responses from the probes are used to eliminate unwanted components in the measured quantities for accurate determination of selected characteristics.
Effect of turbulent atmosphere on the on-axis average intensity of Pearcey–Gaussian beam
NASA Astrophysics Data System (ADS)
F, Boufalah; L, Dalil-Essakali; H, Nebdi; A, Belafhal
2016-06-01
The propagation characteristics of the Pearcey–Gaussian (PG) beam in turbulent atmosphere are investigated in this paper. The Pearcey beam is a new kind of paraxial beam, based on the Pearcey function of catastrophe theory, which describes diffraction about a cusp caustic. By using the extended Huygens–Fresnel integral formula in the paraxial approximation and the Rytov theory, an analytical expression of axial intensity for the considered beam family is derived. Some numerical results for PG beam propagating in atmospheric turbulence are given by studying the influences of some factors, including incident beam parameters and turbulence strengths.
NASA Technical Reports Server (NTRS)
Kerr, R. A.
1983-01-01
In a three dimensional simulation higher order derivative correlations, including skewness and flatness factors, are calculated for velocity and passive scalar fields and are compared with structures in the flow. The equations are forced to maintain steady state turbulence and collect statistics. It is found that the scalar derivative flatness increases much faster with Reynolds number than the velocity derivative flatness, and the velocity and mixed derivative skewness do not increase with Reynolds number. Separate exponents are found for the various fourth order velocity derivative correlations, with the vorticity flatness exponent the largest. Three dimensional graphics show strong alignment between the vorticity, rate of strain, and scalar-gradient fields. The vorticity is concentrated in tubes with the scalar gradient and the largest principal rate of strain aligned perpendicular to the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single curve and a short minus 5/3 spectral regime is observed.
A new statistical tool to study the geometry of intense vorticity clusters in turbulence
NASA Astrophysics Data System (ADS)
Vela-Martin, Alberto; Ishihara, Takashi
2016-04-01
Recent large-scale direct numerical simulations (DNS) of high-Reynolds number (high-Re) turbulence, suggest that strong micro-scale tube-like vortices form clusters in localized thin regions of space. However, to this date no thorough quantitative and statistical analysis of the geometry of such vortical clusters has been conducted. This study is intended to generate new statistical tools to study the shape and dynamics of these intense vorticity and strain regions. We first propose a new method for locating and analysing the geometrical properties of thresholded vortical clusters contained inside boxes of a given size. Second, we use this new tool to investigate the natural presence of intense shear layers and their relevance as geometrical features of high-Re homogeneous turbulence. This new method is applied to the DNS of homogeneous incompressible turbulence with up to 40963 grid points, showing that the geometry of high vorticity regions varies strongly depending on the threshold and on the size of the clusters. In particular for sizes in the inertial range of scales and high thresholds, approximately layer-like structures of vortices are extracted and visualized. Agreement of results with previous observations and known features of turbulence supports the validity of the proposed method to characterize the geometry of intense vorticity and strain regions in high-Re turbulence.
Sectoral Analysis of Wind Speed and Turbulence Intensity over Forest and Open Land
NASA Astrophysics Data System (ADS)
Beyer, Elisabeth
2015-04-01
More and more wind turbines are built onshore and reduce the available areas for wind energy. Forests are an additional potential for wind energy priority areas. But the high roughness of wooden areas and the resulting turbulences make it difficult to assess sites in forests. In order to cope with this problem some measurements, using Lidar and met masts were done inside and outside wooden areas. Therefore, two posters were created, dealing with the dependence of turbulence intensity on wind speed and wind direction. This poster summarises the existing results, concentrating on the sectoral analysis of wind speed and turbulence intensity. Furthermore it investigates the differences of inflow over landcover with low and high surface roughnesses.
NASA Astrophysics Data System (ADS)
Dhariwal, Rohit; Rani, Sarma; Koch, Donald
2015-11-01
In an earlier work, Rani, Dhariwal, and Koch (JFM, Vol. 756, 2014) developed an analytical closure for the diffusion current in the PDF transport equation describing the relative motion of high-Stokes-number particle pairs in isotropic turbulence. In this study, an improved closure was developed for the diffusion coefficient, such that the motion of the particle-pair center of mass is taken into account. Using the earlier and the new analytical closures, Langevin simulations of pair relative motion were performed for four particle Stokes numbers, Stη = 10 , 20 , 40 , 80 and at two Taylor micro-scale Reynolds numbers Reλ = 76 , 131 . Detailed comparisons of the analytical model predictions with those of DNS were undertaken. It is seen that the pair relative motion statistics obtained from the improved theory show excellent agreement with the DNS statistics. The radial distribution functions (RDFs), and relative velocity PDFs obtained from the improved-closure-based Langevin simulations are found to be in very good agreement with those from DNS. It was found that the RDFs and relative velocity RMS increased with Reλ for all Stη . The collision kernel also increased strongly with Reλ , since it depended on the RDF and the radial relative velocities.
NASA Astrophysics Data System (ADS)
Nguyen van yen, Romain; Farge, Marie; Schneider, Kai
2012-02-01
Classical statistical theories of turbulence have shown their limitations, in that they cannot predict much more than the energy spectrum in an idealized setting of statistical homogeneity and stationarity. We explore the applicability of a conditional statistical modeling approach: can we sort out what part of the information should be kept, and what part should be modeled statistically, or, in other words, “dissipated”? Our mathematical framework is the initial value problem for the two-dimensional (2D) Euler equations, which we approximate numerically by solving the 2D Navier-Stokes equations in the vanishing viscosity limit. In order to obtain a good approximation of the inviscid dynamics, we use a spectral method and a resolution going up to 8192 2. We introduce a macroscopic concept of dissipation, relying on a split of the flow between coherent and incoherent contributions: the coherent flow is constructed from the large wavelet coefficients of the vorticity field, and the incoherent flow from the small ones. In previous work, a unique threshold was applied to all wavelet coefficients, while here we also consider the effect of a scale by scale thresholding algorithm, called scale-wise coherent vorticity extraction. We study the statistical properties of the coherent and incoherent vorticity fields, and the transfers of enstrophy between them, and then use these results to propose, within a maximum entropy framework, a simple model for the incoherent vorticity. In the framework of this model, we show that the flow velocity can be predicted accurately in the L2 norm for about 10 eddy turnover times.
Bumblebee Flight in Heavy Turbulence
NASA Astrophysics Data System (ADS)
Engels, T.; Kolomenskiy, D.; Schneider, K.; Lehmann, F.-O.; Sesterhenn, J.
2016-01-01
High-resolution numerical simulations of a tethered model bumblebee in forward flight are performed superimposing homogeneous isotropic turbulent fluctuations to the uniform inflow. Despite tremendous variation in turbulence intensity, between 17% and 99% with respect to the mean flow, we do not find significant changes in cycle-averaged aerodynamic forces, moments, or flight power when averaged over realizations, compared to laminar inflow conditions. The variance of aerodynamic measures, however, significantly increases with increasing turbulence intensity, which may explain flight instabilities observed in freely flying bees.
Effects of q-profile structure on turbulence spreading: A fluctuation intensity transport analysis
Yi, S.; Kwon, J. M.; Diamond, P. H.; Hahm, T. S.
2014-09-15
This paper studies effects of q-profile structure on turbulence spreading. It reports results of numerical experiments using global gyrokinetic simulations. We examine propagation of turbulence, triggered by an identical linear instability in a source region, into an adjacent, linearly stable region with variable q-profile. The numerical experiments are designed so as to separate the physics of turbulence spreading from that of linear stability. The strength of turbulence spreading is measured by the penetration depth of turbulence. Dynamics of spreading are elucidated by fluctuation intensity balance analysis, using a model intensity evolution equation which retains nonlinear diffusion and damping, and linear growth. It is found that turbulence spreading is strongly affected by magnetic shear s, but is hardly altered by the safety factor q itself. There is an optimal range of modest magnetic shear which maximizes turbulence spreading. For high to modest shear values, the spreading is enhanced by the increase of the mode correlation length with decreasing magnetic shear. However, the efficiency of spreading drops for sufficiently low magnetic shear even though the mode correlation length is comparable to that for the case of optimal magnetic shear. The reduction of spreading is attributed to the increase in time required for the requisite nonlinear mode-mode interactions. The effect of increased interaction time dominates that of increased mode correlation length. Our findings of the reduction of spreading and the increase in interaction time at weak magnetic shear are consistent with the well-known benefit of weak or reversed magnetic shear for core confinement enhancement. Weak shear is shown to promote locality, as well as stability.
Turbulence in planetary occultations. IV - Power spectra of phase and intensity fluctuations
NASA Technical Reports Server (NTRS)
Haugstad, B. S.
1979-01-01
Power spectra of phase and intensity scintillations during occultation by turbulent planetary atmospheres are significantly affected by the inhomogeneous background upon which the turbulence is superimposed. Such coupling is particularly pronounced in the intensity, where there is also a marked difference in spectral shape between a central and grazing occultation. While the former has its structural features smoothed by coupling to the inhomogeneous background, such features are enhanced in the latter. Indeed, the latter power spectrum peaks around the characteristic frequency that is determined by the size of the free-space Fresnel zone and the ray velocity in the atmosphere; at higher frequencies strong fringes develop in the power spectrum. A confrontation between the theoretical scintillation spectra computed here and those calculated from the Mariner 5 Venus mission by Woo et al. (1974) is inconclusive, mainly because of insufficient statistical resolution. Phase and/or intensity power spectra computed from occultation data may be used to deduce characteristics of the turbulence and to distinguish turbulence from other perturbations in the refractive index. Such determinations are facilitated if observations are made at two or more frequencies (radio occultation) or in two or more colors (stellar occultation).
Slant path average intensity of finite optical beam propagating in turbulent atmosphere
NASA Astrophysics Data System (ADS)
Zhang, Yixin; Wang, Gaogang
2006-10-01
The average intensity of finite laser beam propagating through turbulent atmosphere is calculated from the extended Huygens Fresnel principle. Formulas are presented for the slant path average intensity from an arbitrarily truncated Gaussian beam. The new expressions are derived from the modified von Karman spectrum for refractive-index fluctuations, quadratic approximation of the structure function, and Gaussian approximation for the product of Gaussian function and Bessel function. It is shown that the form of average intensity is not a Gaussian function but a polynomial of the power of the binomial function, Gaussian function, and the incomplete gamma function. The results also show that the mean irradiance of a finite optical beam propagating in slant path turbulent atmosphere not only depends on the effective beam radius at the transmitting aperture plane, propagation distance, and long-term lateral coherence length of spherical wave, but also on the radius of emit aperture.
NASA Astrophysics Data System (ADS)
Kamardin, A. P.; Nevzorova, I. V.; Odintsov, S. L.
2015-11-01
In the work, we consider estimates of the height of layer of intense turbulent heat exchange in stably stratified atmospheric boundary layer, made with the use of meteorological acoustic radar (sodar). Dependence of this height on temperature gradient is analyzed. Current temperature stratification of the atmosphere in the layer with height up to 1 000 m was determined with the help of MTP-5 meteorological temperature profiler.
NASA Technical Reports Server (NTRS)
Dryden, Hugh L; Schubauer, G B; Mock, W C , Jr; Skramstad, H K
1937-01-01
The investigation of wind-tunnel turbulence, conducted at the National Bureau of Standards with the cooperation of the National Advisory Committee for Aeronautics, has been extended to include a new variable, namely, the scale of the turbulence. This report presents the results of a study of this new variable together with the intensity of the turbulence, and the effect of both on the critical Reynolds number of spheres.
Sankaran, Ramanan; Hawkes, Evatt R.; Yoo, Chun Sang; Chen, Jacqueline H.
2015-06-22
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methane–air chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methane–air mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall on the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.
Sankaran, Ramanan; Hawkes, Evatt R.; Yoo, Chun Sang; Chen, Jacqueline H.
2015-06-22
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methane–air chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methane–air mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall onmore » the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.« less
NASA Astrophysics Data System (ADS)
Hari Prasad, K. B. R. R.; Srinivas, C. V.; Satyanarayana, A. N. V.; Naidu, C. V.; Baskaran, R.; Venkatraman, B.
2015-08-01
Season- and stability-dependent turbulence intensity ( σ u / u *, σ v / u *, σ w / u *) relationships are derived from experimental turbulence measurements following surface layer scaling and local stability at the tropical coastal site Kalpakkam, India for atmospheric dispersion parameterization. Turbulence wind components ( u', v', w') measured with fast response UltraSonic Anemometers during an intense observation campaign for wind field modeling called Round Robin Exercise are used to formulate the flux-profile relationships using surface layer similarity theory and Fast Fourier Transform technique. The new relationships (modified Hanna scheme) are incorporated in a Lagrangian Particle Dispersion model FLEXPART-WRF and tested by conducting simulations for a field tracer dispersion experiment at Kalpakkam. Plume dispersion analysis of a ground level hypothetical release indicated that the new turbulent intensity formulations provide slightly higher diffusivity across the plume relative to the original Hanna scheme. The new formulations for σ u , σ v , σ w are found to give better agreement with observed turbulent intensities during both stable and unstable conditions under various seasonal meteorological conditions. The simulated concentrations using the two methods are compared with those obtained from a classical Gaussian model and the observed SF6 concentration. It has been found that the new relationships provide comparatively higher diffusion across the plume relative to the model default Hanna scheme and provide downwind concentration results in better agreement with observations.
Gross separation approaching a blunt trailing edge as the turbulence intensity increases.
Scheichl, B
2014-07-28
A novel rational description of incompressible two-dimensional time-mean turbulent boundary layer (BL) flow separating from a bluff body at an arbitrarily large globally formed Reynolds number, Re, is devised. Partly in contrast to and partly complementing previous approaches, it predicts a pronounced delay of massive separation as the turbulence intensity level increases. This is bounded from above by a weakly decaying Re-dependent gauge function (hence, the BL approximation stays intact locally), and thus the finite intensity level characterizing fully developed turbulence. However, it by far exceeds the moderate level found in a preceding study which copes with the associated moderate delay of separation. Thus, the present analysis bridges this self-consistent and another forerunner theory, proposing extremely retarded separation by anticipating a fully attached external potential flow. Specifically, it is shown upon formulation of a respective distinguished limit at which rate the separation point and the attached-flow trailing edge collapse as [Formula: see text] and how on a short streamwise scale the typical small velocity deficit in the core region of the incident BL evolves to a large one. Hence, at its base, the separating velocity profile varies generically with the one-third power of the wall distance, and the classical triple-deck problem describing local viscous-inviscid interaction crucial for moderately retarded separation is superseded by a Rayleigh problem, governing separation of that core layer. Its targeted solution proves vital for understanding the separation process more close to the wall. Most importantly, the analysis does not resort to any specific turbulence closure. A first comparison with the available experimentally found positions of separation for the canonical flow past a circular cylinder is encouraging. PMID:24936020
Wake Measurements at alpha ventus - Dependency on Stability and Turbulence Intensity
NASA Astrophysics Data System (ADS)
Westerhellweg, Annette; Cañadillas, Beatriz; Kinder, Friederike; Neumann, Tom
2014-12-01
Wind and power deficit in the wake are assessed for the offshore wind farm Alpha Ventus. Operational data are evaluated for the power deficit in the wake of a single wind turbine and in a row of wind turbines. The wake of a single wind turbine is described by the maximum power deficit and expansion width of the wake. The impact of atmospheric stability in respect to vertical wind shear and turbulence intensity is assessed showing that wake effects are more pronounced under stable conditions.
NASA Astrophysics Data System (ADS)
Ghassemlooy, Zabih; Popoola, Wasiu O.; Ahmadi, Vahid; Leitgeb, Erich
In this paper, we analyse the error performance of transmitter/receiver array free-space optical (FSO) communication system employing binary phase shift keying (BPSK) subcarrier intensity modulation (SIM) in clear but turbulent atmospheric channel. Subcarrier modulation is employed to eliminate the need for adaptive threshold detector. Direct detection is employed at the receiver and each subcarrier is subsequently demodulated coherently. The effect of irradiance fading is mitigated with an array of lasers and photodetectors. The received signals are linearly combined using the optimal maximum ratio combining (MRC), the equal gain combining (EGC) and the selection combining (SelC). The bit error rate (BER) equations are derived considering additive white Gaussian noise and log normal intensity fluctuations. This work is part of the EU COST actions and EU projects.
Wu, Yuqian; Zhang, Yixin; Zhu, Yun
2016-08-01
We studied Gaussian beams with three different partially coherent models, including the Gaussian-Schell model (GSM), Laguerre-Gaussian Schell model (LGSM), and Bessel-Gaussian Schell model (BGSM), propagating through oceanic turbulence. The expressions of average intensity, beam spreading, and beam wander for GSM, LGSM, and BGSM beams in the paraxial channel are derived. We make a contrast for the three models in numerical simulations and find that the GSM beam has smaller spreading than the others, and the LGSM beam needs longer propagation distance to transform into a well-like profile of average intensity than the BGSM beam in the same conditions. The salinity fluctuation has a greater contribution to the wander of LGSM and BGSM beams than that of the temperature fluctuation. Our results can be helpful in the design of an optical wireless communication link operating in oceanic environment. PMID:27505642
NASA Technical Reports Server (NTRS)
Mccolgan, C. J.; Larson, R. S.
1977-01-01
The effect of flight on the mean flow and turbulence properties of a 0.056m circular jet were determined in a free jet wind tunnel. The nozzle exit velocity was 122 m/sec, and the wind tunnel velocity was set at 0, 12, 37, and 61 m/sec. Measurements of flow properties including mean velocity, turbulence intensity and spectra, and eddy convection velocity were carried out using two linearized hot wire anemometers. This report contains the raw data and graphical presentations. The final technical report includes a description of the test facilities, test hardware, along with significant test results and conclusions.
NASA Astrophysics Data System (ADS)
Sutton, David M.
The effect of freestream turbulence intensities ranging from Tu = 1.26% to Tu = 3.2% is studied. Skin friction measurements made on the surface of the airfoil using oil film interferometry (OFI) show that, in general, the effect of the increased Tu is to inhibit separation of the laminar boundary layer. With increased Tu, the near-wall flow experiences strong deceleration in the adverse pressure gradient, but does not reverse as it does in the baseline case where Tu = 0.05%. The Cp distribution resulting from this decelerated fluid is similar in appearance to that of a laminar separation bubble. OFI results also show that laminar separation initiates a more rapid transition process than does higher turbulence intensity: transition of the boundary layer occurs over a shorter distance with Tu = 1.26% than it does with Tu = 2.19% due to the presence of a LSB at the lower turbulence intensity.
NASA Astrophysics Data System (ADS)
Stevens, Richard J. A. M.; Gayme, Dennice; Meneveau, Charles
2015-06-01
We use the recently introduced coupled wake boundary layer (CWBL) model to predict the effect of turbulence intensity on the performance of a wind farm. The CWBL model combines a standard wake model with a “top-down” approach to get improved predictions for the power output compared to a stand-alone wake model. Here we compare the CWBL model results for different turbulence intensities with the Horns Rev field measurements by Hansen et al., Wind Energy 15, 183196 (2012). We show that the main trends as function of the turbulence intensity are captured very well by the model and discuss differences between the field measurements and model results based on comparisons with LES results from Wu and Porté-Agel, Renewable Energy 75, 945-955 (2015).
NASA Astrophysics Data System (ADS)
Goleneva, N. V.; Lavrinov, V.; Lavrinova, L. N.
2015-11-01
The wavefront sensor of Hartmann type consists of two parts: the optical and algorithmic. The parameters of the optical part of the sensor may vary. Since the time of "frozen" turbulence due to the Fried's length and to the cross wind transport turbulent distortion speed, the measurement Shack-Hartmann sensor depend on the intensity of turbulent distortions. In this paper are presented the results of the analysis of the measurements of the sensor according to the size of lens array and to the intensity of turbulent distortions. The analysis is performed on basis of a numerical model of the Shack-Hartmann wavefront sensor and on Kolmogorov's turbulence model.
NASA Astrophysics Data System (ADS)
Frisch, Uriel
1996-01-01
Written five centuries after the first studies of Leonardo da Vinci and half a century after A.N. Kolmogorov's first attempt to predict the properties of flow, this textbook presents a modern account of turbulence, one of the greatest challenges in physics. "Fully developed turbulence" is ubiquitous in both cosmic and natural environments, in engineering applications and in everyday life. Elementary presentations of dynamical systems ideas, probabilistic methods (including the theory of large deviations) and fractal geometry make this a self-contained textbook. This is the first book on turbulence to use modern ideas from chaos and symmetry breaking. The book will appeal to first-year graduate students in mathematics, physics, astrophysics, geosciences and engineering, as well as professional scientists and engineers.
Contribution to the numerical study of turbulence in high intensity discharge lamps
Kaziz, S.; Ben Ahmed, R.; Helali, H.; Gazzah, H.; Charrada, K.
2011-07-15
We present in this paper a comparison between results obtained with a laminar and turbulent models for high-pressure mercury arc. The two models are based on the resolution of bidimensional time-dependent equations by a semi-implicit finite-element code. The numerical computation of turbulent model is solved with large eddy simulation model; this approach takes into account the various scales of turbulence by a filtering method on each scale. The results show the quantitative influence of turbulence on the flow fields and also the difference between laminar and turbulent effects on the dynamic thermal behaviour and on the characteristics of the discharge.
NASA Technical Reports Server (NTRS)
Hoffmann, Jon A.
1988-01-01
The influence of near isotropic free-stream turbulence on the shape factors and skin friction coefficients of turbulent bounday layers is presented for the cases of zero and mild adverse pressure gradients. With free-stream turbulence, improved fluid mixing occurs in boundary layers with adverse pressure gradients relative to the zero pressure gradient condition, with the same free-stream turbulence intensity and length scale. Stronger boundary layers with lower shape factors occur as a result of a lower ratio of the integral scale of turbulence to the boundary layer thickness, and to vortex stretching of the turbulent eddies in the free stream, both of which act to improve the transmission of momentum from the free stream to the boundary layers.
NASA Technical Reports Server (NTRS)
Hoffmann, J. A.; Kassir, S. M.; Larwood, S. M.
1989-01-01
The influence of near isotropic free-stream turbulence on the shape factors and skin friction coefficients of turbulent boundary layers is presented for the cases of zero and mild adverse pressure gradients. With free-stream turbulence, improved fluid mixing occurs in boundary layers with adverse pressure gradients relative to the zero pressure gradient condition, with the same free-stream turbulence intensity and length scale. Stronger boundary layers with lower shape factors occur as a result of a lower ratio of the integral scale of turbulence to the boundary layer thickness, and to vortex stretching of the turbulent eddies in the free-stream, both of which act to improve the transmission of momentum from the free-stream to the boundary layers.
Effects of free-stream turbulence on turbulent boundary layers with mild adverse pressure gradients
NASA Technical Reports Server (NTRS)
Hoffmann, J. A.; Kassir, S. M.
1988-01-01
The influence of near isotropic free-stream turbulence on the shape factors and skin friction coefficients of turbulent boundary layers is presented for the cases of zero and mild adverse pressure gradients. With free-stream turbulence, improved fluid mixing occurs in boundary layers with adverse pressure gradients relative to the zero pressure gradient condition, with the same free-stream turbulence intensity and length scale. Stronger boundary layers with lower shape factors occur as a result of a lower ratio of the integral scale of turbulence to the boundary layer thickness, and to vortex stretching of the turbulent eddies in the free-stream, both of which act to improve the transmission of momentum from the free-stream to the boundary layers.
NASA Astrophysics Data System (ADS)
Ivanova, I. V.; Dmitriev, D. I.; Sirazetdinov, V. S.
2007-02-01
In this paper we analyze some results of natural and numerical experiments on probability density of intensity fluctuations on an axis for 1,06 microns and 0,53 microns laser beams in comparison with theoretical dependences (lognormal, exponential and K-distribution). Beams were propagated in aviation engine exhaust at various angles between the jet and beam axes. It has been shown that for a beam with a wavelength of 0,53 microns experimental data can be approximated as exponential and K-distribution, while for radiation with a wavelength of 1,06 microns good conformity to K-distribution has been observed. Optimum conditions for image registration with CCD-cameras of laser beams distorted by turbulence have been chosen. For this purpose transfer characteristics of several same type samples of CCD-cameras have been studied at various irradiation modes and registration tunings. It has been shown that the dynamic range of the cameras is used to maximum capacity for image recording when gamma-correction is applied.
Experiments on the interaction between hydrodynamic turbulence and surface waves
NASA Astrophysics Data System (ADS)
Jamin, Timothee; Berhanu, Michael; Falcon, Eric
2014-11-01
Different regimes of interaction between hydrodynamic turbulence and a free surface are investigated in a meter scale basin. A homogeneous and isotropic turbulence is generated by an 8×8 array of jets pointing upwards at the bottom of the tank. The 64 jets are driven individually to reach a random spatiotemporal forcing pattern and produce an intense turbulence. Using fluid velocity measurements, we characterize the turbulence obtained with this setup, then we investigate free-surface deformations induced by hydrodynamic turbulence. In a second stage an electromechanical shaker will generate gravity-capillary waves at the free surface. We aim to study reduction or amplification of surface waves and then measure energy exchange between hydrodynamic turbulence and wave turbulence. This work was supported by the DGA-CNRS Ph.D program and ANR Turbulon 12-BS04-0005.
Experimental assessment of helicopter rotor turbulence ingestion noise in hover
NASA Technical Reports Server (NTRS)
Simonich, J. C.; Schlinker, R. H.; Amiet, R. K.
1989-01-01
An experiment was conducted to assess the accuracy of a theory for non-isotropic turbulence ingestion. In order to generate non-isotropic turbulence in a controlled environment, a scale model rotor in a closed chamber was used so that the turbulence generated by the rotor was reingested by the recirculating flow. Simultaneous measurements of turbulence inflow properties and far field acoustics were acquired. Measurements confirmed that the inflow turbulence was highly non-isotropic. The measured aerodynamic properties were used as inputs for the noise prediction procedure. The general agreement between the non-isotropic noise prediction procedure and the experiment was good, although the procedure generally overpredicts the quasi-tonal low to mid range frequencies and underpredicts the higher broadband signals. The predicted sound power level as a function of polar angle was in close agreement with measurements, except near the rotor plane, which is not modeled by the present analysis. It is determined that the most sensitive parameter influencing the predicted noise was the turbulence intensity.
A region of intense plasma wave turbulence on auroral field lines
NASA Technical Reports Server (NTRS)
Gurnett, D. A.; Frank, L. A.
1976-01-01
This report presents a detailed study of the plasma wave turbulence observed by HAWKEYE-1 and IMP-6 on high latitude auroral field lines and investigates the relationship of this turbulence to magnetic field and plasma measurements obtained in the same region.
Measurement of width and intensity of particle streaks in turbulent flows
NASA Astrophysics Data System (ADS)
Håkansson, Karl M. O.; Kvick, Mathias; Lundell, Fredrik; Prahl Wittberg, Lisa; Söderberg, L. Daniel
2013-06-01
Fibre streaks are observed in experiments with fibre suspensions in a turbulent half-channel flow. The preferential concentration methods, most commonly used to quantify preferential particle concentration, are in one dimension found to be concentration dependent. Two different new streak quantification methods are evaluated, one based on Voronoi analysis and the other based on artificial particles with an assigned fixed width. The width of the particle streaks and a measure of the intensity of the streaks, i.e. streakiness, are sought. Both methods are based on the auto-correlation of a signal, generated by summing images in the direction of the streaks. Common for both methods is a severe concentration dependency, verified in experiments keeping the flow conditions constant while the (very dilute) concentration of fibres is altered. The fixed width method is shown to be the most suitable method, being more robust and less computationally expensive. By assuming the concentration dependence to be related to random noise, an expression is derived, which is shown to make the streak width and the streakiness independent of the concentration even at as low concentrations as 0.05 particles per pixel column in an image. The streakiness is obtained by applying an artificial particle width equal to 20 % of the streak width. This artificial particle width is in this study found to be large enough to smoothen the correlation without altering the streakiness nor the streak width. It is concluded that in order to make quantitative comparisons between different experiments or simulations, the evaluation has to be performed with care and be very well documented.
Active control for turbulent premixed flame simulations
Bell, John B.; Day, Marcus S.; Grcar, Joseph F.; Lijewski, Michael J.
2004-03-26
Many turbulent premixed flames of practical interest are statistically stationary. They occur in combustors that have anchoring mechanisms to prevent blow-off and flashback. The stabilization devices often introduce a level of geometric complexity that is prohibitive for detailed computational studies of turbulent flame dynamics. As a result, typical detailed simulations are performed in simplified model configurations such as decaying isotropic turbulence or inflowing turbulence. In these configurations, the turbulence seen by the flame either decays or, in the latter case, increases as the flame accelerates toward the turbulent inflow. This limits the duration of the eddy evolutions experienced by the flame at a given level of turbulent intensity, so that statistically valid observations cannot be made. In this paper, we apply a feedback control to computationally stabilize an otherwise unstable turbulent premixed flame in two dimensions. For the simulations, we specify turbulent in flow conditions and dynamically adjust the integrated fueling rate to control the mean location of the flame in the domain. We outline the numerical procedure, and illustrate the behavior of the control algorithm. We use the simulations to study the propagation and the local chemical variability of turbulent flame chemistry.
Hierarchical structures in a turbulent pipe flow
NASA Astrophysics Data System (ADS)
She, Zhen-Su; Zou, Zhengping; Zhu, Yuanjie; Zhou, Mingde
2003-11-01
Statistical structures of a series of longitudinal velocity fluctuation signals at different distances (10
Mixing and chemical reaction in sheared and nonsheared homogeneous turbulence
NASA Technical Reports Server (NTRS)
Leonard, Andy D.; Hill, James C.
1992-01-01
Direct numerical simulations were made to examine the local structure of the reaction zone for a moderately fast reaction between unmixed species in decaying, homogeneous turbulence and in a homogeneous turbulent shear flow. Pseudospectral techniques were used in domains of 64 exp 3 and higher wavenumbers. A finite-rate, single step reaction between non-premixed reactants was considered, and in one case temperature-dependent Arrhenius kinetics was assumed. Locally intense reaction rates that tend to persist throughout the simulations occur in locations where the reactant concentration gradients are large and are amplified by the local rate of strain. The reaction zones are more organized in the case of a uniform mean shear than in isotropic turbulence, and regions of intense reaction rate appear to be associated with vortex structures such as horseshoe vortices and fingers seen in mixing layers. Concentration gradients tend to align with the direction of the most compressive principal strain rate, more so in the isotropic case.
NASA Astrophysics Data System (ADS)
Pearson, Juli K.
The growing demand for increased efficiency in turbine engine designs has sparked a growing interest for research of air flow around curved surfaces. The turbine's operating conditions result in material property constraints, especially in the first stage turbine vanes and blades. These turbine vane components experience extreme loading conditions of both high temperature and high turbulence intensities exiting the combustor. The surface of the turbine blades has cylindrical leading edges that promote stabilizing flow accelerations. These convex surfaces can cause a reduced eddy diffusivity across the boundary layer. This thesis reviews measurements of velocity and turbulence intensities taken just shy of the thirty degrees offset from the stagnation line of a two-dimensional cylindrical leading edge under a wide range of turbulence and flow conditions flow conditions. Flow conditions and velocity measurements were gathered with respect to the distance to the surface. The length of the measurements extended from the surface to beyond the boundary layer's edge. The instrumentation used to collect data was a single wire driven by a constant temperature anemometer bridge. The hot wire is specially modified to measure data near the cylindrical leading edges curved surface. The traversing system allowed the acquisition of high-resolution boundary layer data. The traversing system was installed internally to the cylindrical leading edge to reduce probe blockage.
Hierarchical structures in a turbulent pipe flow
NASA Astrophysics Data System (ADS)
Zou, Zhengping; Zhu, Yuanjie; Zhou, Mingde; She, Zhen-Su
2003-12-01
A hierarchical structure (HS) analysis ( β-test and γ-test) is applied to a fully developed turbulent pipe flow. Velocity signals are measured at two cross sections in the pipe and at a series of radial locations from the pipe wall. Particular attention is paid to the variation of turbulent statistics at wall units 10< y+<3000. It is shown that at all locations the velocity fluctuations satisfy the She-Leveque hierarchical symmetry (Phys. Rev. Lett. 72 (1994) 336). The measured HS parameters, β and γ, are interpreted in terms of the variation of fluid structures. Intense anisotropic fluid structures generated near the wall appear to be more singular than the most intermittent structures in isotropic turbulence and appear to be more outstanding compared to the background fluctuations; this yields a more intermittent velocity signal with smaller γ and β. As turbulence migrates into the logarithmic region, small-scale motions are generated by an energy cascade and large-scale organized structures emerge which are also less singular than the most intermittent structures of isotropic turbulence. At the center, turbulence is nearly isotropic, and β and γ are close to the 1994 She-Leveque predictions. A transition is observed from the logarithmic region to the center in which γ drops and the large-scale organized structures break down. We speculate that it is due to the growing eddy viscosity effects of widely spread turbulent fluctuations in a similar way as in the breakdown of the Taylor vortices in a turbulent Couette-Taylor flow at high Reynolds numbers.
Turbulent intensity and Reynolds number effects on an airfoil at low Reynolds numbers
NASA Astrophysics Data System (ADS)
Wang, S.; Zhou, Y.; Alam, Md. Mahbub; Yang, H.
2014-11-01
This work investigates the aerodynamics of a NACA 0012 airfoil at the chord-based Reynolds numbers (Rec) from 5.3 × 103 to 2.0 × 104. The lift and drag coefficients, CL and CD, of the airfoil, along with the flow structure, were measured as the turbulent intensity Tu of oncoming flow varies from 0.6% to 6.0%. The analysis of the present data and those in the literature unveils a total of eight distinct flow structures around the suction side of the airfoil. Four Rec regimes, i.e., the ultra-low (<1.0 × 104), low (1.0 × 104-3.0 × 105), moderate (3.0 × 105-5.0 × 106), and high Rec (>5.0 × 106), are proposed based on their characteristics of the CL-Rec relationship and the flow structure. It has been observed that Tu has a more pronounced effect at lower Rec than at higher Rec on the shear layer separation, reattachment, transition, and formation of the separation bubble. As a result, CL, CD, CL/CD and their dependence on the airfoil angle of attack all vary with Tu. So does the critical Reynolds number Rec,cr that divides the ultra-low and low Rec regimes. It is further noted that the effect of increasing Tu bears similarity in many aspects to that of increasing Rec, albeit with differences. The concept of the effective Reynolds number Rec,eff advocated for the moderate and high Rec regimes is re-evaluated for the low and ultra-low Rec regimes. The Rec,eff treats the non-zero Tu effect as an addition of Rec and is determined based on the presently defined Rec,cr. It has been found that all the maximum lift data from both present measurements and previous reports collapse into a single curve in the low and ultra-low Rec regimes if scaled with Rec,eff.
Turbulence generation in homogeneous dilute particle- laden flows
NASA Astrophysics Data System (ADS)
Chen, Jeng-Horng
Homogeneous turbulence generated by the motion of particles in dispersed multiphase flows was studied both theoretically and experimentally, motivated by applications to sprays, particle-laden jets, bubble plumes and rainstorms, among others. The experiments involved uniform fluxes of monodisperse spherical particles falling through a slow upflow of air. Particle fluxes and phase velocities were measured by sampling and phase-discriminating laser Doppler velocimetry (LDV), respectively. Measured particle velocities included mean and fluctuating streamwise and cross-stream velocities and probability density functions (PDF's). Measured continuous-phase velocities included mean and fluctuating streamwise and cross-stream velocities, PDF's and the higher moments of velocity fluctuations such as skewness and kurtosis, energy spectra of velocity fluctuations and integral length scales based on streamwise velocity fluctuations. Continuous-phase velocity measurements included conditional averages for particle wake disturbances and the turbulent inter-wake region surrounding these disturbances as well as overall flow properties. Present and earlier results in the literature provided particle Reynolds numbers of 38-990, particle volume fractions less than 0.01% and turbulence intensities (normalized by mean particle relative velocities) of 0.1-10.0%. Theory included characterization of particle wake disturbances as laminar-like turbulent wakes observed for intermediate particle Reynolds numbers in turbulent environments, characterization of the turbulent inter-wake region by analogy to grid-generated isotropic turbulence, and estimation of overall flow properties by conditional averaging of the properties of the wake disturbances and the turbulent inter-wake region. Present measurements showed that particle wake disturbances during turbulence generation were properly characterized by the properties of laminar-like turbulent wakes. The turbulent inter-wake region was
Probe shapes that measure time-averaged streamwise momentum and cross-stream turbulence intensity
NASA Technical Reports Server (NTRS)
Rossow, Vernon J. (Inventor)
1993-01-01
A method and apparatus for directly measuring the time-averaged streamwise momentum in a turbulent stream use a probe which has total head response which varies as the cosine-squared of the angle of incidence. The probe has a nose with a slight indentation on its front face for providing the desired response. The method of making the probe incorporates unique design features. Another probe may be positioned in a side-by-side relationship to the first probe to provide a direct measurement of the total pressure. The difference between the two pressures yields the sum of the squares of the cross-stream components of the turbulence level.
Wu, Pey-Shey; Tsai, Shen-Ta; Jhuo, Yue-Hua
2014-01-01
This study is concerned with a film cooling technique applicable to the protection of the endwalls of a gas turbine vane. In the experiments, cross-injection coolant flow from two-row, paired, inclined holes with nonintersecting centerlines was utilized. The test model is a scaled two-half vane. The levels of turbulence intensity used in the experiments are T.I. = 1.8%, 7%, and 12%. Other parameters considered in the film cooling experiments include three inlet Reynolds numbers (9.20 × 10(4), 1.24 × 10(5), and 1.50 × 10(5)), three blowing ratios (0.5, 1.0, and 2.0), and three endwall conditions (smooth endwall and stepped endwall with forward-facing or backward-facing step). Thermochromic liquid crystal (TLC) technique with steady-state heat transfer experiments was used to obtain the whole-field film cooling effectiveness. Results show that, at low turbulence intensity, increasing Reynolds number decreases the effectiveness in most of the vane passage. There is no monotonic trend of influence by Reynolds number at high turbulence intensity. The effect of blowing ratio on the effectiveness has opposite trends at low and high turbulence levels. Increasing turbulent intensity decreases the effectiveness, especially near the inlet of the vane passage. With a stepped endwall, turbulence intensity has only mild effect on the film cooling effectiveness. PMID:24592153
Tsai, Shen-Ta; Jhuo, Yue-Hua
2014-01-01
This study is concerned with a film cooling technique applicable to the protection of the endwalls of a gas turbine vane. In the experiments, cross-injection coolant flow from two-row, paired, inclined holes with nonintersecting centerlines was utilized. The test model is a scaled two-half vane. The levels of turbulence intensity used in the experiments are T.I. = 1.8%, 7%, and 12%. Other parameters considered in the film cooling experiments include three inlet Reynolds numbers (9.20 × 104, 1.24 × 105, and 1.50 × 105), three blowing ratios (0.5, 1.0, and 2.0), and three endwall conditions (smooth endwall and stepped endwall with forward-facing or backward-facing step). Thermochromic liquid crystal (TLC) technique with steady-state heat transfer experiments was used to obtain the whole-field film cooling effectiveness. Results show that, at low turbulence intensity, increasing Reynolds number decreases the effectiveness in most of the vane passage. There is no monotonic trend of influence by Reynolds number at high turbulence intensity. The effect of blowing ratio on the effectiveness has opposite trends at low and high turbulence levels. Increasing turbulent intensity decreases the effectiveness, especially near the inlet of the vane passage. With a stepped endwall, turbulence intensity has only mild effect on the film cooling effectiveness. PMID:24592153
Banakh, V A; Marakasov, D A
2008-04-30
An algorithm for the wind profile recovery from spatiotemporal spectra of a laser beam reflected in a turbulent atmosphere is presented. The cases of a spherical wave incident on a diffuse reflector of finite size and a spatially limited beam reflected from an infinite random surface are considered. (laser applications and other topics in quantum electronics)
NASA Astrophysics Data System (ADS)
Melnick, M. Blake; Thurow, Brian S.
2014-02-01
Simultaneous particle image velocimetry (PIV) and flow visualization measurements were performed in a turbulent boundary layer in an effort to better quantify the relationship between the velocity field and the image intensity typically observed in a classical flow visualization experiment. The freestream flow was lightly seeded with smoke particles to facilitate PIV measurements, whereas the boundary layer was densely seeded with smoke through an upstream slit in the wall to facilitate both PIV and classical flow visualization measurements at Reynolds numbers, Re θ , ranging from 2,100 to 8,600. Measurements were taken with and without the slit covered as well as with and without smoke injection. The addition of a narrow slit in the wall produces a minor modification of the nominal turbulent boundary layer profile whose effect is reduced with downstream distance. The presence of dense smoke in the boundary layer had a minimal effect on the observed velocity field and the associated proper orthogonal decomposition (POD) modes. Analysis of instantaneous images shows that the edge of the turbulent boundary layer identified from flow visualization images generally matches the edge of the boundary layer determined from velocity and vorticity. The correlation between velocity deficit and smoke intensity was determined to be positive and relatively large (>0.7) indicating a moderate-to-strong relationship between the two. This notion was extended further through the use of a direct correlation approach and a complementary POD/linear stochastic estimation (LSE) approach to estimate the velocity field directly from flow visualization images. This exercise showed that, in many cases, velocity fields estimated from smoke intensity were similar to the actual velocity fields. The complementary POD/LSE approach proved better for these estimations, but not enough to suggest using this technique to approximate velocity measurements from a smoke intensity image. Instead, the
Effects of axisymmetric contractions on turbulence of various scales
NASA Technical Reports Server (NTRS)
Tan-Atichat, J.; Nagib, H. M.; Drubka, R. E.
1980-01-01
Digitally acquired and processed results from an experimental investigation of grid generated turbulence of various scales through and downstream of nine matched cubic contour contractions ranging in area ratio from 2 to 36, and in length to inlet diameter ratio from 0.25 to 1.50 are reported. An additional contraction with a fifth order contour was also utilized for studying the shape effect. Thirteen homogeneous and nearly isotropic test flow conditions with a range of turbulence intensities, length scales and Reynolds numbers were generated and used to examine the sensitivity of the contractions to upstream turbulence. The extent to which the turbulence is altered by the contraction depends on the incoming turbulence scales, the total strain experienced by the fluid, as well as the contraction ratio and the strain rate. Varying the turbulence integral scale influences the transverse turbulence components more than the streamwise component. In general, the larger the turbulence scale, the lesser the reduction in the turbulence intensity of the transverse components. Best agreement with rapid distortion theory was obtained for large scale turbulence, where viscous decay over the contraction length was negligible, or when a first order correction for viscous decay was applied to the results.
NASA Technical Reports Server (NTRS)
Kiock, R.
1978-01-01
Turbulence intensity (Tu) measurements were made in two-dimensional and rotating cascades of blades in a low-speed cascade wind tunnel using hot-wire probes as sensors. The local Tu at Re = 1.6x100000 was determined in the wake zone behind a two-dimensional cascade. Then the values were recomputed for a rotating cascade, giving a mean turbulence intensity of 6.5% at 1/10 chord downstream and 2.9% at one chord. Fans were used for measurements on the rotating cascade. Re was equal to 7x100000. Frequency analysis was employed to separate the actual Tu of the entry flow from the effects caused by interaction with the rotor blades, showing that the true Tu increased from a few tenths of a percent to 6.6% in the 1st rotor, and from 7.2 to 9.3% in the 2d rotor. The Tu behind the 3d rotor was equal to 8.9%.
NASA Astrophysics Data System (ADS)
Ji, Hongjun; Wang, Qiang; Li, Mingyu; Wang, Chunqing
2014-07-01
With the miniaturization of portable electronic devices, the size of solder joint interconnects is decreasing to micrometer levels. These joints possess only several or even one or two grains, resulting in anisotropy and failure issues. Direct ultrasound-assisted solidification of Cu/SAC305/Cu interconnects for grain refinement and fabrication of isotropic solder joints is presented herein. These joints consist of many β-Sn grains. The average cross-sectional area of the Sn-rich phase is significantly reduced by up to 99% when compared with conventional as-reflowed samples. The ultrasonic power density exhibits a threshold value for affecting the microstructures. Below 200 W cm-2, the β-Sn grains were refined and had circular shape. The Ag3Sn phase grew in a manner similar to branched coral to sizes reaching 30 μm, or as rods aggregated together with Cu6Sn5 tube fragments. Above 200 W cm-2, the microstructures were coarsened and Ag3Sn had plate-like shape. The thickness of Cu6Sn5 intermetallic layers at the Cu/solder interfaces was reduced by more than 26%. The relationships among the ultrasonic power, nucleation rate, local temperature drop, and pressure were identified. At the highest power density of 267 W cm-2, the nucleation rate was about 4.05 × 1014 m-3 s-1, the local temperature drop was 248 K, and the local pressure was on the order of several GPa.
Noise of a model helicopter rotor due to ingestion of turbulence
NASA Technical Reports Server (NTRS)
Paterson, R. W.; Amiet, R. K.
1979-01-01
A theoretical and experimental investigation of the noise of a model helicoper rotor due to ingestion of turbulence was conducted. Experiments were performed with a 0.76 m dia, articulated model rotor for a range of inflow turbulence and rotor operating conditions. Inflow turbulence levels varied from approximately 2 to 19 percent and tip Mach number was varied from 0.3 to 0.52. Test conditions included ingestion of a atmospheric turbulence in outdoor hover as well as ingestion of grid generated isotropic turbulence in the wind tunnel airstream. In wind tunnel testing, both forward flight and vertical ascent (climb) were simulated. Far field noise spectra and directivity were measured in addition to incident turbulence intensities, length scales, and spectra. Results indicate that ingestion of atmospheric turbulence is the dominant helicopter rotor hover noise mechanism at the moderate to high frequencies which determine perceived noise level.
A Jet-Stirred Apparatus for Turbulent Combustion Experiments
NASA Astrophysics Data System (ADS)
Davani, Abbasali; Ronney, Paul
2015-11-01
A novel jet-stirred combustion chamber is designed to study turbulent premixed flames. In the new approach, multiple impinging turbulent jets are used to stir the mixture. It is well known that pair of counterflowing turbulent jets produces nearly a constant intensity (u') along the jet axes. In this study, different numbers of impinging jets in various configurations are used to produce isotropic turbulence intensity. FLUENT simulations have been conducted to assess the viability of the proposed chamber. In order to be able to compare different configurations, three different non dimensional indices are introduces. Mean flow index; Homogeneity index, and Isotropicity index. Using these indices one can compare various chambers including conventional Fan-stirred Reactors. Results show that a concentric inlet/outlet chamber (CAIO) with 8 inlets and 8 outlets with inlet velocity of 20 m/s and initial intensity of 15% produces near zero mean flow and 2.5 m/s turbulence intensity which is much more higher than reported values for Fan-stirred chamber. This research was sponsored by National Science Foundation.
Turbulence Intensity at Inlet of 80- by 120-Foot Wind Tunnel Caused by Upwind Blockage
NASA Technical Reports Server (NTRS)
Salazar, Denise; Yuricich, Jillian
2014-01-01
In order to estimate the magnitude of turbulence in the National Full-Scale Aerodynamics Complex (NFAC) 80- by 120-Foot Wind Tunnel (80 x 120) caused by buildings located upwind from the 80 x 120 inlet, a 150th-scale study was performed that utilized a nominal two-dimensional blockage placed ahead of the inlet. The distance of the blockage ahead of the inlet was varied. This report describes velocity measurements made in the plane of the 80 x 120 model inlet for the case of zero ambient (atmospheric) wind.
Analysis of two-equation turbulence models for recirculating flows
NASA Technical Reports Server (NTRS)
Thangam, S.
1991-01-01
The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.
NASA Astrophysics Data System (ADS)
Lapenta, Giovanni; Kemel, Koen; Henri, Pierre; Califano, Francesco; Markidis, Stefano
2015-11-01
Using the full kinetic implicit PIC code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Kelvin-Helmholtz vortex chain via secondary instabilities, vortex pairing and non-linear interactions. We found that the magnetic energy spectral cascade between ion and electron inertial scales is in agreement with satellite observations and previous numerical simulations; however, in our case the spectrum ends with a peak beyond de due to the occurrence of the lower hybrid drift instability. The electric energy spectrum is influenced by secondary instabilities: anomalous resistivity, fed by the development of the lower hybrid drift instability, steepens the spectral decay and, depending on the alignment of B and the shear vorticity, peaks due to ion-Bernstein waves may dominate the spectrum around di. A key conclusion of the study is that the anomalous resistivity produced by these complex wave and instabilities can indeed very accurately be described in terms of a proportionality with the current. This research used resources of NERSC, a DOE Office of Science User Facility supported by the Office of Science of the U.S. DOE under Contract No. DE-AC02-05CH11231.
NASA Technical Reports Server (NTRS)
Mielke, Amy F.; Seasholtz, Richard G.; Elam, Kristie A.; Panda, Jayanta
2004-01-01
A molecular Rayleigh scattering based flow diagnostic is developed to measure time average velocity, density, temperature, and turbulence intensity in a 25.4-mm diameter nozzle free jet facility. The spectrum of the Rayleigh scattered light is analyzed using a Fabry-Perot interferometer operated in the static imaging mode. The resulting fringe pattern containing spectral information of the scattered light is recorded using a low noise CCD camera. Nonlinear least squares analysis of the fringe pattern using a kinetic theory model of the Rayleigh scattered light provides estimates of density, velocity, temperature, and turbulence intensity of the gas flow. Resulting flow parameter estimates are presented for an axial scan of subsonic flow at Mach 0.95 for comparison with previously acquired pitot tube data, and axial scans of supersonic flow in an underexpanded screeching jet. The issues related to obtaining accurate turbulence intensity measurements using this technique are discussed.
Extreme events in computational turbulence
Yeung, P. K.; Zhai, X. M.; Sreenivasan, Katepalli R.
2015-01-01
We have performed direct numerical simulations of homogeneous and isotropic turbulence in a periodic box with 8,1923 grid points. These are the largest simulations performed, to date, aimed at improving our understanding of turbulence small-scale structure. We present some basic statistical results and focus on “extreme” events (whose magnitudes are several tens of thousands the mean value). The structure of these extreme events is quite different from that of moderately large events (of the order of 10 times the mean value). In particular, intense vorticity occurs primarily in the form of tubes for moderately large events whereas it is much more “chunky” for extreme events (though probably overlaid on the traditional vortex tubes). We track the temporal evolution of extreme events and find that they are generally short-lived. Extreme magnitudes of energy dissipation rate and enstrophy occur simultaneously in space and remain nearly colocated during their evolution. PMID:26424452
Structure of wind-shear turbulence
NASA Technical Reports Server (NTRS)
Trevino, G.; Laituri, T. R.
1989-01-01
The statistical characteristics of wind shear turbulence are modelled. Isotropic turbulence serves as the basis of comparison for the anisotropic turbulence which exists in wind shear. The question of turbulence scales in wind shear is addressed from the perspective of power spectral density.
Structure of wind-shear turbulence
NASA Technical Reports Server (NTRS)
Trevino, G.; Laituri, T. R.
1988-01-01
The statistical characteristics of wind-shear turbulence are modelled. Isotropic turbulence serves as the basis of comparison for the anisotropic turbulence which exists in wind shear. The question of how turbulence scales in a wind shear is addressed from the perspective of power spectral density.
NASA Technical Reports Server (NTRS)
Flegel, Ashlie B.; Giel, Paul W.; Welch, Gerard E.
2014-01-01
The effects of high inlet turbulence intensity on the aerodynamic performance of a variable speed power turbine blade are examined over large incidence and Reynolds number ranges. These results are compared to previous measurements made in a low turbulence environment. Both high and low turbulence studies were conducted in the NASA Glenn Research Center Transonic Turbine Blade Cascade Facility. The purpose of the low inlet turbulence study was to examine the transitional flow effects that are anticipated at cruise Reynolds numbers. The current study extends this to LPT-relevant turbulence levels while perhaps sacrificing transitional flow effects. Assessing the effects of turbulence at these large incidence and Reynolds number variations complements the existing database. Downstream total pressure and exit angle data were acquired for 10 incidence angles ranging from +15.8deg to -51.0deg. For each incidence angle, data were obtained at five flow conditions with the exit Reynolds number ranging from 2.12×10(exp 5) to 2.12×10(exp 6) and at a design exit Mach number of 0.72. In order to achieve the lowest Reynolds number, the exit Mach number was reduced to 0.35 due to facility constraints. The inlet turbulence intensity, Tu, was measured using a single-wire hotwire located 0.415 axial-chord upstream of the blade row. The inlet turbulence levels ranged from 8 to 15 percent for the current study. Tu measurements were also made farther upstream so that turbulence decay rates could be calculated as needed for computational inlet boundary conditions. Downstream flow field measurements were obtained using a pneumatic five-hole pitch/yaw probe located in a survey plane 7 percent axial chord aft of the blade trailing edge and covering three blade passages. Blade and endwall static pressures were acquired for each flow condition as well. The blade loading data show that the suction surface separation that was evident at many of the low Tu conditions has been eliminated. At
NASA Technical Reports Server (NTRS)
Flegel, Ashlie B.; Giel, Paul W.; Welch, Gerard E.
2014-01-01
The effects of high inlet turbulence intensity on the aerodynamic performance of a variable speed power turbine blade are examined over large incidence and Reynolds number ranges. These results are compared to previous measurements made in a low turbulence environment. Both high and low turbulence studies were conducted in the NASA Glenn Research Center Transonic Turbine Blade Cascade Facility. The purpose of the low inlet turbulence study was to examine the transitional flow effects that are anticipated at cruise Reynolds numbers. The current study extends this to LPT-relevant turbulence levels while perhaps sacrificing transitional flow effects. Assessing the effects of turbulence at these large incidence and Reynolds number variations complements the existing database. Downstream total pressure and exit angle data were acquired for 10 incidence angles ranging from +15.8deg to -51.0deg. For each incidence angle, data were obtained at five flow conditions with the exit Reynolds number ranging from 2.12×10(exp 5) to 2.12×10(exp 6) and at a design exit Mach number of 0.72. In order to achieve the lowest Reynolds number, the exit Mach number was reduced to 0.35 due to facility constraints. The inlet turbulence intensity, Tu, was measured using a single-wire hotwire located 0.415 axial-chord upstream of the blade row. The inlet turbulence levels ranged from 8 to 15 percent for the current study. Tu measurements were also made farther upstream so that turbulence decay rates could be calculated as needed for computational inlet boundary conditions. Downstream flow field measurements were obtained using a pneumatic five-hole pitch/yaw probe located in a survey plane 7 percent axial chord aft of the blade trailing edge and covering three blade passages. Blade and endwall static pressures were acquired for each flow condition as well. The blade loading data show that the suction surface separation that was evident at many of the low Tu conditions has been eliminated. At
The effects of anisotropic free-stream turbulence on turbulent boundary layer behavior
NASA Technical Reports Server (NTRS)
Liang-Wei, F.; Hoffman, J. A.
1985-01-01
The effects of near-isotropic and highly anisotropic free-stream turbulence on mean flow properties of the turbulence structure of turbulent boundary layers in a near zero pressure gradient flow has been experimentally evaluated. Turbulence levels vary from 0.5% to 8.0% and the momentum thickness Reynolds number varies from 800 to 1100. The results indicate that the effects of free-stream turbulence on the classical boundary layer properties for near-isotropic turbulence which have been published by other investigators are similar to the case of highly anisotropic turbulence fields, while the effects of free-stream turbulence on the properties of the turbulent structure within the boundary layer for the case of near-isotropic turbulence are quite different compared to the highly anisotropic case.
Direct simulation of compressible turbulence in a shear flow
NASA Technical Reports Server (NTRS)
Sarkar, S.; Erlebacher, G.; Hussaini, M. Y.
1991-01-01
Compressibility effects on the turbulence in homogeneous shear flow are investigated. The growth of the turbulent kinetic energy was found to decrease with increasing Mach number: a phenomenon which is similar to the reduction of turbulent velocity intensities observed in experiments on supersonic free shear layers. An examination of the turbulent energy budget shows that both the compressible dissipation and the pressure-dilatation contribute to the decrease in the growth of kinetic energy. The pressure-dilatation is predominantly negative in homogeneous shear flow, in contrast to its predominantly positive behavior in isotropic turbulence. The different signs of the pressure-dilatation are explained by theoretical consideration of the equations for the pressure variance and density variance. Previously, the following results were obtained for isotropic turbulence: (1) the normalized compressible dissipation is of O(M(sub t)(exp 2)); and (2) there is approximate equipartition between the kinetic and potential energies associated with the fluctuating compressible mode. Both of these results were substantiated in the case of homogeneous shear. The dilatation field is significantly more skewed and intermittent than the vorticity field. Strong compressions seem to be more likely than strong expansions.
Coherence in Turbulence: New Perspective
NASA Astrophysics Data System (ADS)
Levich, Eugene
2009-07-01
It is claimed that turbulence in fluids is inherently coherent phenomenon. The coherence shows up clearly as strongly correlated helicity fluctuations of opposite sign. The helicity fluctuations have cellular structure forming clusters that are actually observed as vorticity bands and coherent structures in laboratory turbulence, direct numerical simulations and most obviously in atmospheric turbulence. The clusters are named BCC - Beltrami Cellular Clusters - because of the observed nearly total alignment of the velocity and vorticity fields in each particular cell, and hence nearly maximal possible helicity in each cell; although when averaged over all the cells the residual mean helicity in general is small and does not play active dynamical role. The Beltrami like fluctuations are short-lived and stabilize only in small and generally contiguous sub-domains that are tending to a (multi)fractal in the asymptotic limit of large Reynolds numbers, Re → ∞. For the model of homogeneous isotropic turbulence the theory predicts the leading fractal dimension of BCC to be: DF = 2.5. This particular BCC is responsible for generating the Kolmogorov -5/3 power law energy spectrum. The most obvious role that BCC play dynamically is that the nonlinear interactions in them are relatively reduced, due to strong spatial alignment between the velocity field v(r, t) and the vorticity field ω(r, t) = curlv(r, t), while the physical quantities typically best characterizing turbulence intermittency, such as entrophy, vorticity stretching and generation, and energy dissipation are maximized in and near them. The theory quantitatively relates the reduction of nonlinear inter-actions to the BCC fractal dimension DF and subsequent turbulence intermittency. It is further asserted that BCC is a fundamental feature of all turbulent flows, e.g., wall bounded turbulent flows, atmospheric and oceanic flows, and their leading fractal dimension remains invariant and universal in these flows
NASA Technical Reports Server (NTRS)
Flegel, Ashlie Brynn; Giel, Paul W.; Welch, Gerard E.
2014-01-01
The effects of inlet turbulence intensity on the aerodynamic performance of a variable speed power turbine blade are examined over large incidence and Reynolds number ranges. Both high and low turbulence studies were conducted in the NASA Glenn Research Center Transonic Turbine Blade Cascade Facility. The purpose of the low inlet turbulence study was to examine the transitional flow effects that are anticipated at cruise Reynolds numbers. The high turbulence study extends this to LPT-relevant turbulence levels while perhaps sacrificing transitional flow effects. Downstream total pressure and exit angle data were acquired for ten incidence angles ranging from +15.8 to 51.0. For each incidence angle, data were obtained at five flow conditions with the exit Reynolds number ranging from 2.12105 to 2.12106 and at a design exit Mach number of 0.72. In order to achieve the lowest Reynolds number, the exit Mach number was reduced to 0.35 due to facility constraints. The inlet turbulence intensity, Tu, was measured using a single-wire hotwire located 0.415 axial-chord upstream of the blade row. The inlet turbulence levels ranged from 0.25 - 0.4 for the low Tu tests and 8- 15 for the high Tu study. Tu measurements were also made farther upstream so that turbulence decay rates could be calculated as needed for computational inlet boundary conditions. Downstream flow field measurements were obtained using a pneumatic five-hole pitchyaw probe located in a survey plane 7 axial chord aft of the blade trailing edge and covering three blade passages. Blade and endwall static pressures were acquired for each flow condition as well. The blade loading data show that the suction surface separation that was evident at many of the low Tu conditions has been eliminated. At the extreme positive and negative incidence angles, the data show substantial differences in the exit flow field. These differences are attributable to both the higher inlet Tu directly and to the thinner inlet endwall
NASA Technical Reports Server (NTRS)
Haugstad, B. S.
1978-01-01
The nature and magnitude of turbulence-induced errors in atmospheric profiles derived from Doppler measurements made during radio occultations are investigated. It is found that turbulence in planetary atmospheres induces both fluctuating and systematic errors in derived profiles, but the errors of both types are very small. Consideration of the occultation of Mariner 10 by Venus and of the Pioneer occultations by Jupiter shows that the rms fractional errors in the atmospheric profiles derived from these observations were less than 0.01 in both temperature and pressure, while the fractional systematic errors were typically of the order of 1 millionth. The extent to which atmospheric profiles derived from radio and optical intensity measurements are affected by turbulence is also examined. The results indicate that turbulence in planetary atmospheres has only a marginal effect on derived profiles in the weak-scattering limit and that the turbulence-induced errors in this case are always much larger than the corresponding errors in profiles derived from radio Doppler measurements.
NASA Astrophysics Data System (ADS)
Madala, Srikanth; Satyanarayana, A. N. V.; Srinivas, C. V.
2015-12-01
Accurate representation of air pollutant dispersion is essential for environmental management and planning purposes. In this study, semi-empirical relationships of turbulence intensity (σu/u*, σv/u* and σw/u*) as a function of surface layer scaling and local stability are developed following boundary layer similarity concepts at Ranchi, a complex terrain in Jharkhand, Eastern India for various seasons. The impact of the new turbulence parameterization for air pollution dispersion simulation is studied by incorporating the same in the Hanna scheme of FLEXPART-WRF Lagrangian Particle dispersion model over study region. The model is used to estimate the ground level concentrations of nitrogen oxides (NOx) due to industrial and vehicular sources in study region. The meteorological parameters needed in air-quality simulation are simulated using the Advanced Research WRF (ARW) mesoscale model at high resolution (3 km). Three turbulence schemes (YSU, MYNN2 and ACM2) in ARW are alternatively tested in dispersion simulation and comparisons are made with available air quality data for eight days in different seasons (winter, pre-monsoon, monsoon and post-monsoon). Simulations with FLEXPART revealed distinct seasonal variation of dispersion patterns. It has been found that the new turbulence intensity relationships in FLEXPART improved the NOx concentration estimates by reducing the negative bias seen with default Hanna scheme. Further, the ARW simulated meteorological parameters using ACM2 and MYNN2 significantly reduced the bias in modeled pollutant concentrations. The study demonstrates the utility of high quality seasonal turbulence measurements in pollution dispersion model for better diffusion parameterization needed in air quality modeling.
Scaling of turbulent flame speed for expanding flames with Markstein diffusion considerations
NASA Astrophysics Data System (ADS)
Chaudhuri, Swetaprovo; Wu, Fujia; Law, Chung K.
2013-09-01
In this paper we clarify the role of Markstein diffusivity, which is the product of the planar laminar flame speed and the Markstein length, on the turbulent flame speed and its scaling, based on experimental measurements on constant-pressure expanding turbulent flames. Turbulent flame propagation data are presented for premixed flames of mixtures of hydrogen, methane, ethylene, n-butane, and dimethyl ether with air, in near-isotropic turbulence in a dual-chamber, fan-stirred vessel. For each individual fuel-air mixture presented in this work and the recently published iso-octane data from Leeds, normalized turbulent flame speed data of individual fuel-air mixtures approximately follow a ReT,f0.5 scaling, for which the average radius is the length scale and thermal diffusivity is the transport property of the turbulence Reynolds number. At a given ReT,f, it is experimentally observed that the normalized turbulent flame speed decreases with increasing Markstein number, which could be explained by considering Markstein diffusivity as the leading dissipation mechanism for the large wave number flame surface fluctuations. Consequently, by replacing thermal diffusivity with the Markstein diffusivity in the turbulence Reynolds number definition above, it is found that normalized turbulent flame speeds could be scaled by ReT,M0.5 irrespective of the fuel, equivalence ratio, pressure, and turbulence intensity for positive Markstein number flames.
Scaling of turbulent flame speed for expanding flames with Markstein diffusion considerations.
Chaudhuri, Swetaprovo; Wu, Fujia; Law, Chung K
2013-09-01
In this paper we clarify the role of Markstein diffusivity, which is the product of the planar laminar flame speed and the Markstein length, on the turbulent flame speed and its scaling, based on experimental measurements on constant-pressure expanding turbulent flames. Turbulent flame propagation data are presented for premixed flames of mixtures of hydrogen, methane, ethylene, n-butane, and dimethyl ether with air, in near-isotropic turbulence in a dual-chamber, fan-stirred vessel. For each individual fuel-air mixture presented in this work and the recently published iso-octane data from Leeds, normalized turbulent flame speed data of individual fuel-air mixtures approximately follow a Re_{T,f}^{0.5} scaling, for which the average radius is the length scale and thermal diffusivity is the transport property of the turbulence Reynolds number. At a given Re_{T,f}^{}, it is experimentally observed that the normalized turbulent flame speed decreases with increasing Markstein number, which could be explained by considering Markstein diffusivity as the leading dissipation mechanism for the large wave number flame surface fluctuations. Consequently, by replacing thermal diffusivity with the Markstein diffusivity in the turbulence Reynolds number definition above, it is found that normalized turbulent flame speeds could be scaled by Re_{T,M}^{0.5} irrespective of the fuel, equivalence ratio, pressure, and turbulence intensity for positive Markstein number flames. PMID:24125342
Isotropic sequence order learning.
Porr, Bernd; Wörgötter, Florentin
2003-04-01
In this article, we present an isotropic unsupervised algorithm for temporal sequence learning. No special reward signal is used such that all inputs are completely isotropic. All input signals are bandpass filtered before converging onto a linear output neuron. All synaptic weights change according to the correlation of bandpass-filtered inputs with the derivative of the output. We investigate the algorithm in an open- and a closed-loop condition, the latter being defined by embedding the learning system into a behavioral feedback loop. In the open-loop condition, we find that the linear structure of the algorithm allows analytically calculating the shape of the weight change, which is strictly heterosynaptic and follows the shape of the weight change curves found in spike-time-dependent plasticity. Furthermore, we show that synaptic weights stabilize automatically when no more temporal differences exist between the inputs without additional normalizing measures. In the second part of this study, the algorithm is is placed in an environment that leads to closed sensor-motor loop. To this end, a robot is programmed with a prewired retraction reflex reaction in response to collisions. Through isotropic sequence order (ISO) learning, the robot achieves collision avoidance by learning the correlation between his early range-finder signals and the later occurring collision signal. Synaptic weights stabilize at the end of learning as theoretically predicted. Finally, we discuss the relation of ISO learning with other drive reinforcement models and with the commonly used temporal difference learning algorithm. This study is followed up by a mathematical analysis of the closed-loop situation in the companion article in this issue, "ISO Learning Approximates a Solution to the Inverse-Controller Problem in an Unsupervised Behavioral Paradigm" (pp. 865-884). PMID:12689389
Inlet free-stream turbulence effects on diffuser performance
NASA Technical Reports Server (NTRS)
Hoffman, J. A.; Gonzales, G.
1983-01-01
The performance of a subsonic two dimensional diffuser was experimentally evaluated as a function of inlet free-stream turbulence parameters. Anisotropic inlet free-stream turbulence with the eddy axis perpendicular to the flow and parallel to the diverging walls of the diffuser appears to be more effective at transmitting energy to the diverging walls of the diffuser, thereby improving diffuser performance, as compared to isotropic turbulence or anisotropic turbulence with the eddy axis perpendicular to the diverging walls of the diffuser. The pressure recovery of the diffuser was found to be strongly dependent upon the inlet free-stream total turbulence intensity, was independent of eddy size for large eddy dimensions, and was dependent upon eddy size for small eddy dimensions. The improvement in the diffuser's static pressure recovery coefficient at a total included divergence angle of 20 deg, compared to the low inlet turbulence case, was found to be as much as 21 times larger than the pressure loss across the turbulence generators.
NASA Astrophysics Data System (ADS)
Madala, Srikanth; Satyanarayana A. N., V.; Srinivas C., V.; Boadh, Rahul; Pinaka Pani V. V. S., N.; Kumar, Manoj
2015-04-01
The complex terrain region of Patratu, Jharkhand in southern Chota Nagpur of eastern India has high air pollution problems besides complex mesoscale flow and meteorology. The FLEXPART-WRF mesoscale Lagrangian Particle dispersion model is used to simulate the dispersion of elevated effluent releases of nitrogen dioxide (NO2) and suspended particulate matter (SPM) from Patratu thermal power plant over Patratu at a high resolution of 1 km. The WRF is integrated with nested domains (27, 9, 3 km resolutions, 51 vertical levels). The relationships for turbulent intensities in the default diffusion parameterization of the Hanna scheme of FLEXPART is modified with new empirical relationships derived as a function of atmospheric stability from one year fast response turbulence measurements from a nearby observational site at Ranchi. The pollutant dispersion simulated by FLEXPART is evaluated with modified version of the model and using the WRF simulated atmospheric flow field and thermodynamical structure with three alternative PBL schemes [Yonsei University (YSU), Asymmetric Convective Model version 2 (ACM2) and Mellor- Yamada Nakanishi and Niino Level 2.5 PBL (MYNN2]. Results indicate that the new turbulence intensity relationships in FLEXPART provide better comparisons for concentrations of NO2 and SPM with available observations relative to the default relationships. Further, the meteorological parameters simulated using YSU significantly reduces the bias in modeled pollutant concentrations in terms of lesser mean absolute error (MAE), root mean square error (RMSE), normalized mean square error (NMSE), fractional bias (FB) and FAC2 (Factor of 2). These parametric tests enabled to fine tune and validate the FLEXPART-WRF dispersion model with YSU PBL physics and improved Hanna relationships to realistically simulate pollution dispersion over complex terrain of the study region. The study demonstrates the utility of high quality turbulence measurements in pollution
NASA Astrophysics Data System (ADS)
Zhang, Qiang
The effects of surface roughness, turbulence intensity, Mach number, and streamline curvature-airfoil shape on the aerodynamic performance of turbine airfoils are investigated in compressible, high speed flows. The University of Utah Transonic Wind Tunnel is employed for the experimental part of the study. Two different test sections are designed to produce Mach numbers, Reynolds numbers, passage mass flow rates, and physical dimensions, which match values along turbine blades in operating engines: (i) a nonturning test section with a symmetric airfoil, and (ii) a cascade test section with a cambered turbine vane. The nonuniform, irregular, three-dimensional surface roughness is characterized using the equivalent sand grain roughness size. Changing the airfoil surface roughness condition has a substantial effect on wake profiles of total pressure loss coefficients, normalized Mach number, normalized kinetic energy, and on the normalized and dimensional magnitudes of Integrated Aerodynamic Losses produced by the airfoils. Comparisons with results for a symmetric airfoil and a cambered vane show that roughness has more substantial effects on losses produced by the symmetric airfoil than the cambered vane. Data are also provided that illustrate the larger loss magnitudes are generally present with flow turning and cambered airfoils, than with symmetric airfoils. Wake turbulence structure of symmetric airfoils and cambered vanes are also studied experimentally. The effects of surface roughness and freestream turbulence levels on wake distributions of mean velocity, turbulence intensity, and power spectral density profiles and vortex shedding frequencies are quantified one axial chord length downstream of the test airfoils. As the level of surface roughness increases, all wake profile quantities broaden significantly and nondimensional vortex shedding frequencies decrease. Wake profiles produced by the symmetric airfoil are more sensitive to variations of surface
Turbulent Sediment Suspension and Induced Ripple Dynamics Absent Mean Shear
NASA Astrophysics Data System (ADS)
Johnson, B. A.; Cowen, E.
2014-12-01
The uprush and backwash phases in the swash zone, the region of the beach that is alternately covered and uncovered by wave run-up, are fundamentally different events. Backwash is dominated by a growing boundary layer where the turbulence is set by the bed shear stress. In this phase traditional boundary layer turbulence models and Shields-type critical stress pickup functions work well. However, the uprush phase, while often viewed in the context of traditional boundary layer turbulence models, has little in common with the backwash phase. During uprush, the entire water column is turbulent, as it rapidly advects well-stirred highly turbulent flow generated offshore from breaking waves or collapsing bores. Turbulence levels in the uprush are several times higher than turbulent boundary layer theory would predict and hence the use of a boundary layer model to predict turbulence levels during uprush grossly under predicts the turbulence and subsequent sediment suspension in the swash zone. To study the importance of this advected turbulence to sediment suspension we conduct experiments in a water tank designed to generate horizontally homogeneous isotropic turbulence absent mean shear using randomly actuated synthetic jet arrays suspended above both a solid glass plate and a narrowly graded sediment bed. Using jet arrays with different jet spacings allows the generation of high Reynolds number turbulence with variable integral length scales, which we hypothesize control the characteristic length scales in the induced ripple field. Particle image velocimetry and acoustic Doppler velocimetry measurements are used to characterize the near-bed flow and this unique turbulent boundary layer. Metrics include the mean flow and turbulence intensities and stresses, temporal and spatial spectra, dissipation of turbulent kinetic energy, and integral length scales of the turbulence. We leverage our unique dataset to compare the flows over impermeable fixed and permeable mobile
Oren, P.E.; Lee, R.L.; Tek, M.R.
1988-09-01
This paper represents the first contribution in qualifying the effects of wellbore storage along with those of skin and turbulence intensity on deliverability from gas wells. A decay factor, W/sub D/, has been introduced to account for the contribution of sandface flow rate to early-time transient production from the gas wells. Generalized charts for type-curve studies on gas wells, based on pseudopressure and pseudotime concepts, along with wellbore storage, are presented with applications and support from field data.
NASA Technical Reports Server (NTRS)
Seasholtz, R. G.; Goldman, L. J.
1982-01-01
A technique for measuring a small optical axis velocity component in a flow with a large transverse velocity component is presented. Experimental results are given for a subsonic free jet operating in a laboratory environment, and for a 0.508 meter diameter turbine stator cascade. Satisfactory operation of the instrument was demonstrated in the stator cascade facility with an ambient acoustic noise level during operation of about 105 dB. In addition, the turbulence intensity measured with the interferometer was consistent with previous measurements taken with a fringe type laser anemometer.
NASA Astrophysics Data System (ADS)
Brunk, Brett Kenneth
1997-11-01
Pollutant and particle transport in estuaries is affected by a multitude of physical, chemical and biological processes. In this research the importance of equilibrium sorption and turbulent coagulation were studied. Sorption in estuaries was modeled using phenanthrene, bacterial extracellular polymer and kaolinite clay as surrogates for a hydrophobic organic pollutant, dissolved organic matter and inorganic suspended sediment, respectively. Experiments over a range of estuarine salinities showed that ionic strength had the largest effect on the extent of sorption, while the effect of extracellular polymer coatings on the mineral surfaces was insignificant. Further calculations using typical estuarine suspended sediment concentrations indicated that equilibrium sorption could not fully account for the solid/solution phase distribution of hydrophobic organic compounds in the estuarine water column. For particles that are small compared to the length scales of turbulence, the rate of coagulation is related to the dynamics of the smallest turbulent eddies since they have the highest shear rate. Experimental and theoretical effort focused on determining the coagulation rate of spherical particles in isotropic turbulence. A pair diffusion approximation valid for rapidly fluctuating flows was used to calculate the rate of coagulation in a randomly varying isotropic linear flow field. Dynamic simulations of particle coagulation in Gaussian turbulence were computed over a range of representative values of particle-particle interactions (i.e, hydrodynamic interactions and van der Waals attraction) and total strain (i.e., the product of the strain rate and its time scale). The computed coagulation rates for isotropic turbulence differed from analytical approximations valid at large and small total strain. As expected, particle interactions were found to be significant. Experimental measurements of coagulation in grid-generated turbulence were obtained by measuring the loss
NASA Astrophysics Data System (ADS)
Amato, Alberto; Fortini, Stefania; Watteaux, Romain; Diano, Marcello; Espa, Stefania; Esposito, Serena; Ferrante, Maria I.; Peters, Francesc; Iudicone, Daniele; Ribera d'Alcalà, Maurizio
2016-03-01
In recent years, there has been a renewed interest in the impact of turbulence on aquatic organisms. In response to this interest, a novel instrument has been constructed, TURBOGEN, that generates turbulence in water volumes up to 13 l. TURBOGEN is fully computer controlled, thus, allowing for a high level of reproducibility and for variations of the intensity and characteristics of turbulence during the experiment. The calibration tests, carried out by particle image velocimetry, showed TURBOGEN to be successful in generating isotropic turbulence at the typical relatively low levels of the marine environment. TURBOGEN and its sizing have been devised with the long-term scope of analyzing in detail the molecular responses of plankton to different mixing regimes, which is of great importance in both environmental and biotechnological processes.
Amato, Alberto; Fortini, Stefania; Watteaux, Romain; Diano, Marcello; Espa, Stefania; Esposito, Serena; Ferrante, Maria I; Peters, Francesc; Iudicone, Daniele; Ribera d'Alcalà, Maurizio
2016-03-01
In recent years, there has been a renewed interest in the impact of turbulence on aquatic organisms. In response to this interest, a novel instrument has been constructed, TURBOGEN, that generates turbulence in water volumes up to 13 l. TURBOGEN is fully computer controlled, thus, allowing for a high level of reproducibility and for variations of the intensity and characteristics of turbulence during the experiment. The calibration tests, carried out by particle image velocimetry, showed TURBOGEN to be successful in generating isotropic turbulence at the typical relatively low levels of the marine environment. TURBOGEN and its sizing have been devised with the long-term scope of analyzing in detail the molecular responses of plankton to different mixing regimes, which is of great importance in both environmental and biotechnological processes. PMID:27036831
Suppression of turbulent resistivity in turbulent Couette flow
Si, Jiahe Sonnenfeld, Richard G.; Colgate, Arthur S.; Westpfahl, David J.; Romero, Van D.; Martinic, Joe; Colgate, Stirling A.; Li, Hui; Nornberg, Mark D.
2015-07-15
Turbulent transport in rapidly rotating shear flow very efficiently transports angular momentum, a critical feature of instabilities responsible both for the dynamics of accretion disks and the turbulent power dissipation in a centrifuge. Turbulent mixing can efficiently transport other quantities like heat and even magnetic flux by enhanced diffusion. This enhancement is particularly evident in homogeneous, isotropic turbulent flows of liquid metals. In the New Mexico dynamo experiment, the effective resistivity is measured using both differential rotation and pulsed magnetic field decay to demonstrate that at very high Reynolds number rotating shear flow can be described entirely by mean flow induction with very little contribution from correlated velocity fluctuations.
Suppression of turbulent resistivity in turbulent Couette flow
NASA Astrophysics Data System (ADS)
Si, Jiahe; Colgate, Stirling A.; Sonnenfeld, Richard G.; Nornberg, Mark D.; Li, Hui; Colgate, Arthur S.; Westpfahl, David J.; Romero, Van D.; Martinic, Joe
2015-07-01
Turbulent transport in rapidly rotating shear flow very efficiently transports angular momentum, a critical feature of instabilities responsible both for the dynamics of accretion disks and the turbulent power dissipation in a centrifuge. Turbulent mixing can efficiently transport other quantities like heat and even magnetic flux by enhanced diffusion. This enhancement is particularly evident in homogeneous, isotropic turbulent flows of liquid metals. In the New Mexico dynamo experiment, the effective resistivity is measured using both differential rotation and pulsed magnetic field decay to demonstrate that at very high Reynolds number rotating shear flow can be described entirely by mean flow induction with very little contribution from correlated velocity fluctuations.
Orientation statistics and settling velocity of ellipsoids in decaying turbulence
NASA Astrophysics Data System (ADS)
Siewert, C.; Kunnen, R. P. J.; Meinke, M.; Schröder, W.
2014-06-01
Motivated by applications in technology as well as in other disciplines where the motion of particles in a turbulent flow field is important, the orientation and settling velocity of ellipsoidal particles in a spatially decaying isotropic turbulent flow are numerically investigated. With respect to cloud microphysics ellipsoidal particles of various shapes are interpreted as archetypes of regular ice crystals, i.e., plates and columns approximated by oblate and prolate ellipsoids. The motion of 19 million small and heavy ellipsoidal particles is tracked by a Lagrangian point-particle model based on Stokes flow conditions. Five types of ellipsoids of revolution such as prolates, spheres, and oblates are considered. The orientation and settling velocity statistics are gathered at six turbulence intensities characterized by the turbulent kinetic energy dissipation rate ranging from 30 to 250 cm2s- 3. It is shown that the preferential orientation of ellipsoids is disturbed by the turbulent fluctuations of the fluid forces and moments. As the turbulence intensity increases the orientation probability distribution becomes more and more uniform. That is, the settling velocity of the ellipsoids is influenced by the turbulence level since the drag force is dependent on the orientation. The effect is more pronounced, the longer the prolate or the flatter the oblate is. The theoretical settling velocity based on the orientation probability of the non-spherical particles is smaller than that found in the simulation. The results show the existence of the preferential sweeping phenomenon also for non-spherical particles. These two effects of turbulence on the motion of ellipsoids change the settling velocity and as such the swept volume, that is expected to result in modified collision probabilities of ellipsoid-shaped particles.
Pulsating instability and self-acceleration of fast turbulent flames
NASA Astrophysics Data System (ADS)
Poludnenko, Alexei Y.
2015-01-01
A series of three-dimensional numerical simulations is used to study the intrinsic stability of high-speed turbulent flames. Calculations model the interaction of a fully resolved premixed flame with a highly subsonic, statistically steady, homogeneous, isotropic turbulence. The computational domain is unconfined to prevent the onset of thermoacoustic instabilities. We consider a wide range of turbulent intensities and system sizes, corresponding to the Damköhler numbers Da = 0.1 - 6.0. These calculations show that turbulent flames in the regimes considered are intrinsically unstable. In particular, we find three effects. (1) Turbulent flame speed, ST, develops pulsations with the observed peak-to-peak amplitude ST max / ST min > 10 and a characteristic time scale close to a large-scale eddy turnover time. Such variability is caused by the interplay between turbulence, which continuously creates the flame surface, and highly intermittent flame collisions, which consume the flame surface. (2) Unstable burning results in the periodic pressure build-up and the formation of pressure waves or shocks, when ST approaches or exceeds the speed of a Chapman-Jouguet deflagration. (3) Coupling of pressure gradients formed during pulsations with density gradients across the flame leads to the anisotropic amplification of turbulence inside the flame volume and flame acceleration. Such process, which is driven by the baroclinic term in the vorticity transport equation, is a reacting-flow analog of the mechanism underlying the Richtmyer-Meshkov instability. With the increase in turbulent intensity, the limit-cycle instability discussed here transitions to the regime described in our previous work, in which the growth of ST becomes unbounded and produces a detonation.
Kentfield, J.A.C.
1996-10-01
A study is presented of apparent drag anomalies based on the results of measurements made by many workers. In some cases the addition, to an airfoil, of a Gurney flap adds little or nothing to the airfoil drag coefficient; yet in other cases, for a prescribed Gurney-flap height-to-chord ratio, the addition of the Gurney flap can increase the airfoil drag coefficient,m at low lift coefficients, by a factor as large as three. The lift increment obtained is essentially independent of the drag level recorded and is, to a first approximation, proportional to the height-to-chord ratio of the Gurney flap. It is, as a consequence of the study, shown that the major adverse influence on the drag coefficient of Gurney-flap equipped airfoils is the thinness of the pressure-surface boundary layer at the airfoil trailing edge. It has been shown by other workers that provided the height of a Gurney flap is less than the boundary layer thickness, the drag penalty due to the Gurney flap is small or can be essentially zero. The influence of the chord-based Reynolds number on the pressure-surface boundary layer thickness appears to parallel that associated with a turbulent flow over a flat plate. Hence the boundary layer thickness diminishes with increasing chord-based Reynolds number. A factor that is shown to have a major impact on the boundary layer thickness is the inherent turbulence level of the approaching flow. Generally, the higher the turbulence level and the lower the Reynolds number, the thicker, for a prescribed airfoil chord, the boundary layer.
NASA Astrophysics Data System (ADS)
Chang, Keh-Chin; Chen, Ching-Shun
1993-03-01
A hybrid k-epsilon turbulence model, based on the concept that the modification of anisotropic effects should not be made in the flow regions inherent to small streamline curvatures, has been developed and examined with the swirling recirculating flows, with the swirl levels ranging from 0.6 to 1.23 in an abrupt pipe expansion. A fairly satisfactory agreement of model predictions with the experimental data shows that this hybrid k-epsilon model can perform better simulation of swirling recirculating flows as compared to the standard k-epsilon model and the modified k-epsilon model proposed by Abujelala and Lilley (1984).
NASA Technical Reports Server (NTRS)
Kao, S. K.; Lordi, N. J.
1977-01-01
Analyses of the meteorological rocket data obtained from an experiment conducted at 3-hour intervals at 8 western meridional rocket stations are presented. Large variations in the meridional wind contribute substantially to overall turbulence in the tropical stratosphere. The solar semidiurnal component of wind oscillations in the tropics was observed to be much higher than predicted by theory, often exceeding the magnitude of the diurnal amplitude throughout the stratosphere. The observed value of the solar diurnal amplitude in the stratosphere was in line with theoretical prediction. The solar terdiurnal amplitudes for temperature, meridional and zonal winds were non-negligible and must be considered in any harmonic analysis. Phase angle variation with height was rapid for all harmonics; however, there was general agreement between predicted and observed phase angles. Because of large changes in the mean winds in the mesosphere with season, harmonic determinations are difficult. There appear to be large zonal wind changes even within the same season as mentioned previously. Turbulence diffusivity in the upper stratosphere is greater near the equator than in the mid-latitudes.
Isotropic Monte Carlo Grain Growth
Energy Science and Technology Software Center (ESTSC)
2013-04-25
IMCGG performs Monte Carlo simulations of normal grain growth in metals on a hexagonal grid in two dimensions with periodic boundary conditions. This may be performed with either an isotropic or a misorientation - and incliantion-dependent grain boundary energy.
NASA Astrophysics Data System (ADS)
Smith, Craig M.; Barthelmie, R. J.; Pryor, S. C.
2013-09-01
Observations of wakes from individual wind turbines and a multi-megawatt wind energy installation in the Midwestern US indicate that directly downstream of a turbine (at a distance of 190 m, or 2.4 rotor diameters (D)), there is a clear impact on wind speed and turbulence intensity (TI) throughout the rotor swept area. However, at a downwind distance of 2.1 km (26 D downstream of the closest wind turbine) the wake of the whole wind farm is not evident. There is no significant reduction of hub-height wind speed or increase in TI especially during daytime. Thus, in high turbulence regimes even very large wind installations may have only a modest impact on downstream flow fields. No impact is observable in daytime vertical potential temperature gradients at downwind distances of >2 km, but at night the presence of the wind farm does significantly decrease the vertical gradients of potential temperature (though the profile remains stably stratified), largely by increasing the temperature at 2 m.
NASA Astrophysics Data System (ADS)
Dmitrenko, Artur V.
2016-07-01
The stochastic equations of continuum are used for determining the hydraulic drag coefficients. As a result, the formulas for the hydraulic drag coefficients dependent on the turbulence intensity and scale instead of only on the Reynolds number are proposed for the classic flows of an incompressible fluid along a smooth flat plate and a round smooth tube. It is shown that the new expressions for the classical drag coefficients, which depend only on the Reynolds number, should be obtained from these new general formulas if to use the well-known experimental data for the initial turbulence. It is found that the limitations of classical empirical and semiempirical formulas for the hydraulic drag coefficients and their deviation from the experimental data depend on different parameters of initial fluctuations in the flow for different experiments in a wide range of Reynolds numbers. On the basis of these new dependencies, it is possible to explain that the differences between the experimental results for the fixed Reynolds number are caused by the difference in the values of flow fluctuations for each experiment instead of only due to the systematic error in the processing of experiments. Accordingly, the obtained general dependencies for the smooth flat plate and the smooth round tube can serve as the basis for clarifying the results of experiments and the experimental formulas, which used for continuum flows in different devices.
NASA Technical Reports Server (NTRS)
Sewell, Jesse; Chew, Larry
1994-01-01
In recent years, the interest in developing a high-speed civil transport has increased. This has led to an increase in research activity on compressible supersonic flows, in particular the boundary layer. The structure of subsonic boundary layers has been extensively documented using conditional sampling techniques which exploit the knowledge of both u and v velocities. Researchers using these techniques have been able to explore some of the complex three-dimensional motions which are responsible for Reynolds stress production and transport in the boundary layer. As interest in turbulent structure has grown to include supersonic flows, a need for simultaneous multicomponent velocity measurements in these flows has developed. The success of conditional analysis in determining the characteristics of coherent motions and structures in the boundary layer relies on accurate, simultaneous measurement of two instantaneous velocity components.
Transversely isotropic poroelasticity arising from thin isotropic layers
Berryman, J.G.
1996-11-01
Percolation phenomena play central roles in the field of poroelasticity, where two distinct sets of percolating continua intertwine. A connected solid frame forms the basis of the elastic behavior of a poroelastic medium in the presence of confining forces, while connected pores permit a percolating fluid (if present) to influence the mechanical response of the system from within. The present paper discusses isotropic and anisotropic poroelastic media and establishes general formulas for the behavior of transversely isotropic poroelasticity arising from laminations of isotropic components. The Backus averaging method is shown to provide elementary means of constructing general formulas. The results for confined fluids are then compared with the more general Gassmann formulas that must be satisfied by any anisotropic poroelastic medium and found to be in complete agreement.
Charged Particle Diffusion in Isotropic Random Static Magnetic Fields
NASA Astrophysics Data System (ADS)
Subedi, P.; Sonsrettee, W.; Matthaeus, W. H.; Ruffolo, D. J.; Wan, M.; Montgomery, D.
2013-12-01
Study of the transport and diffusion of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider Diffusion of charged particles in fully three dimensional statistically isotropic magnetic field turbulence with no mean field which is pertinent to many astrophysical situations. We classify different regions of particle energy depending upon the ratio of Larmor radius of the charged particle to the characteristic outer length scale of turbulence. We propose three different theoretical models to calculate the diffusion coefficient each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing very good agreement.
Turbulent solutions of equations of fluid motion
NASA Technical Reports Server (NTRS)
Deissler, R. G.
1985-01-01
Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.
Turbulence generation by waves
Kaftori, D.; Nan, X.S.; Banerjee, S.
1995-12-31
The interaction between two-dimensional mechanically generated waves, and a turbulent stream was investigated experimentally in a horizontal channel, using a 3-D LDA synchronized with a surface position measuring device and a micro-bubble tracers flow visualization with high speed video. Results show that although the wave induced orbital motion reached all the way to the wall, the characteristics of the turbulence wall structures and the turbulence intensity close to the wall were not altered. Nor was the streaky nature of the wall layer. On the other hand, the mean velocity profile became more uniform and the mean friction velocity was increased. Close to the free surface, the turbulence intensity was substantially increased as well. Even in predominantly laminar flows, the introduction of 2-D waves causes three dimensional turbulence. The turbulence enhancement is found to be proportional to the wave strength.
Turbulence modelling of flow fields in thrust chambers
NASA Astrophysics Data System (ADS)
Chen, C. P.; Kim, Y. M.; Shang, H. M.
1993-02-01
Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.
Turbulence modelling of flow fields in thrust chambers
NASA Technical Reports Server (NTRS)
Chen, C. P.; Kim, Y. M.; Shang, H. M.
1993-01-01
Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.
Interaction of a free flame front with a turbulence field
NASA Technical Reports Server (NTRS)
Tucker, Maurice
1956-01-01
Small-perturbation spectral-analysis techniques are used to obtain the root-mean-square flame-generated turbulence velocities and the attenuating pressure fluctuations stemming from interaction of a constant-pressure flame front with a field of isotropic turbulence in the absence of turbulence decay processes.
Subtraction threshold for an isotropic fluorescence emission difference microscope
NASA Astrophysics Data System (ADS)
Wang, Nan; Kobayashi, Takayoshi
2015-12-01
Isotropic fluorescence emission difference microscopy proposed recently provides a simple method to enhance the spatial resolution in three-dimensions (3D) for fluorescence imaging. However, the subtraction threshold to achieve the condition for appropriately resolving the sample in 3D have not been studied. Then the subtraction factors used in this type of microscopes are still experientially chosen. Based on vector diffraction theory and a 3D numerical model developed here, the subtraction threshold is numerically investigated for the isotropic fluorescence subtraction microscopy. The subtraction factors and peak intensities at the threshold are obtained and comparied both in lateral and axial planes for achieving most appropriate subtraction and inspecting the isotropic characteristic. The effects of radius ratios of implemented 0-π annular phase plate for generating three dimensional donut spot on the subtracted resolution, peak intensity and negative sidebands are also discussed.