Coherent structures in localized and global pipe turbulence

The recent discovery of unstable traveling waves (TWs) in pipe flow has been hailed as a significant breakthrough with the hope that they populate the turbulent attractor. We confirm the existence of coherent states with internal fast and slow streaks commensurate in both structure and energy with known TWs using numerical simulations in a long pipe. These only occur, however, within less energetic regions of (localized) "puff" turbulence at low Reynolds numbers (Re=2000-2400), and not at all in (homogeneous) "slug" turbulence at Re=2800. This strongly suggests that all currently known TWs sit in an intermediate region of phase space between the laminar and turbulent states rather than being embedded within the turbulent attractor itself. New coherent fast streak states with strongly decelerated cores appear to populate the turbulent attractor instead.     

Lower branch coherent states in shear flows: transition and control

Lower branch coherent states in plane Couette flow have an asymptotic structure that consists of O(1) streaks, O(R(-1)) streamwise rolls and a weak sinusoidal wave that develops a critical layer, for large Reynolds number R. Higher harmonics become negligible. These unstable lower branch states appear to have a single unstable eigenvalue at all Reynolds numbers. These results suggest that lower branch coherent states control transition to turbulence and that they may be promising targets for new turbulence prevention strategies.     

Power fluctuations in a closed turbulent shear flow

We study experimentally the power consumption P of a confined turbulent flow at constant Reynolds number Re. We analyze in details its temporal dynamics and statistical properties, in a setup that covers two decades in Reynolds numbers. We show that nontrivial power fluctuations occur over a wide range of amplitudes and that they involve coherent fluid motions over the entire system size. As a result, the power fluctuations do not result from averaging of independent subsystems and its probability density function Pi(P) is strongly non-Gaussian. The shape of Pi(P) is Reynolds number independant and we show that the relative intensity of fluctuations decreases very slowly as Re increases. These results are discussed in terms of an analogy with critical phenomena.     

Harbingers and latecomers – the order of appearance of exact coherent structures in plane Poiseuille flow

The transition to turbulence in plane Poiseuille flow (PPF) is connected with the presence of exact coherent structures. We here discuss a variety of different structures that are relevant for the transition, compare the critical Reynolds numbers and optimal wavelengths for their appearance, and explore the differences between flows operating at constant mass flux or at constant pressure drop. The Reynolds numbers quoted here are based on the mean flow velocity and refer to constant mass flux. Reynolds numbers based on constant pressure drop are always higher. The Tollmien–Schlichting (TS) waves bifurcate subcritically from the laminar profile at Re = 5772 at wavelength 6.16 and reach down to Re = 2610 at a different optimal wave length of 4.65. Their streamwise localised counter part bifurcates at the even lower value Re = 2334. Three-dimensional exact solutions appear at much lower Reynolds numbers. We describe one exact solutions that has a critical Reynolds number of 316. Streamwise localised versions of this state require higher Reynolds numbers, with the lowest bifurcation occurring near Re = 1018. The analysis shows that the various branches of TS-waves cannot be connected with transition observed near Re ≈ 1000 and that the exact coherent structures related to downstream vortices come in at lower Reynolds numbers and prepare for the transition.     

Flow separation behind a rib in a channel with laminar flow

Visualization data and results of combined measurements of flow quantities in flow with separation past a rib at nominally laminar regime of channel flow are reported. In the separation region, the flow is found to be essentially three-dimensional and unsteady, exhibiting a distinct cellular structure and flow zones with transverse motion. It is shown that the rib-induced flow separation gives rise to low-frequency fluctuations of flow velocity and initiates the turbulence transition in the channel flow. The critical Reynolds number at which flow instability starts developing in the channel is estimated. It is shown that at Reynolds numbers higher than the critical Reynolds number the linear integral scale of flow velocity fluctuations in the channel is defined by the duct size.     

转捩射流中涡结构与颗粒扩散的直接模拟

为考察转捩射流中拟序结构的空间演化过程及其对不同Stokes颗粒扩散的影响,采用有限容积方法和分步投影算法对三维气固两相射流进行了直接模拟。其中流体控制方程组的时间积分采用低存储三阶精度的Runge-Kutta格式; 颗粒的跟踪在拉格朗日框架下进行。模拟结果发现,在流场拟序结构由大尺度转化为小尺度的过程中,中、小Stokes数的颗粒能自发地调整它们的扩散方式,分别由非均匀状态向均匀状态以及由均匀状态向非均匀状态转变。     

Numerical simulation and flow visualization using soap film of the self-organized vortex structure in the wake of an array of cylinders

Abstract  

With the aid of computational fluid dynamics (CFD) and simple flow visualization technique using flowing soap-film, we present here the wake structures behind an array of cylinders for Reynolds numbers corresponding to both laminar and turbulent flow regimes. The image results illustrate interesting vortex interactions past these equally spaced cylinders; for low Reynolds number flow, well-organized wake pattern persists and manifests unsteadily to different symmetry states. An increase of free stream flow velocity causes the wake transition, resulting in the formation of asymmetric flow wake with chaotic mixing at the far wake. Observations from both the numerical simulations and soap-film are in good agreement at least qualitatively.     

Decoherence of quantum excitation of even/odd coherent states in thermal environment

In this paper, we study the decoherence of quantum excitation (photon-added) even /odd coherent states, \(((\hat a{^{\dagger }})^{m} \left | {\alpha _ \pm } \right \rangle )\), in a thermal environment by investigating the variation of negative part of the Wigner quasidistribution function vs. the rescaled time. For this purpose, at first we obtain the time-dependent Wigner function corresponding to the mentioned states in the framework of standard master equation. Then, the time evolution of the Wigner function associated with photon-added even /odd coherent states, as well as the number of added photons m are analysed. It is shown that, in both states, the negative part of the Wigner function decreases with time. By deriving the threshold value of the rescaled time for single photon-added even /odd coherent states, it is also found that, if the rescaled time exceeds the threshold value, the associated Wigner function becomes positive, i.e., the decoherence occurs completely.     

Splitting dynamics of coherent structures in a transitional round-pipe flow

Direct numerical modeling of scenarios of the initiation, development, and mutual transitions of coherent turbulent structures in a round pipe at Reynolds numbers 1800 ≤ Re ≤ 4000 is performed. High-accuracy computations were performed for very prolonged time intervals, which made is possible to reveal the fundamental long-lived and statistically stationary flow modes in the transient region between the laminar and established turbulent modes. Reclassification of coherent structures describing the splitting dynamics of the subcritical laminar-turbulent transition is proposed.     

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale(SGS) model based on coherent structures and temporal approximate deconvolution(MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation(LES) of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence(FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation(DNS) results. Compared with the LES results using the temporal approximate deconvolution model(TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.     

Reduction of quantum phase fluctuations in intermediate states

Recently we have shown that the reduction of the Carruthers-Nieto symmetric quantum phase fluctuation parameter (U) with respect to its coherent state value corresponds to an antibunched state, but the converse is not true. Consequently reduction of U is a stronger criterion of nonclassicality than the lowest order antibunching. Here we have studied the possibilities of reduction of U in intermediate states by using the Barnett-Pegg formalism. We have shown that the reduction of phase fluctuation parameter U can be seen in different intermediate states, such as binomial state, generalized binomial state, hypergeometric state, negative binomial state, and photon added coherent state. It is also shown that the depth of nonclassicality can be controlled by various parameters related to intermediate states. Further, we have provided specific examples of antibunched states, for which U is greater than its Poissonian state value.     

Origin of effectiveness degradation in active drag reduction control of turbulent channel flow at Reτ = 1000

Direct numerical simulations of turbulent channel flows are performed with opposition control at Re_τ = 180 and 1000. The drag reduction rate at the higher Reynolds number is reduced by 25% compared with that at the lower Reynolds number. In order to investigate the reason for the degradation of the control effectiveness, we examine the response of Reynolds stresses and coherent structures in both the outer and inner regions to the control and the role that large-scale motions play therein. In the outer region, the Reynolds stresses at different length scales are reduced at the same rate as the drag reduction rate, and conditionally averaged large-scale motions with spanwise scale larger than half channel width are still large-scale low-speed streaks flanked by a pair of large-scale counter-rotating streamwise vortices but with reduced velocity amplitudes. In the inner region, the effectiveness of the control in suppressing the turbulence deteriorates at the higher Reynolds number. In response to the superimposition effect of large-scale motions, the contribution to near-wall wall-parallel velocity fluctuations from large-scale motions becomes larger at the higher Reynolds number, while the suppression of large-scale motions by the control is weaker than that of near-wall coherent structures. In both controlled and uncontrolled cases, large-scale motions can modulate the amplitudes of near-wall coherent structures, and the attenuation of streamwise vortices by the control under large-scale high-speed streaks is significantly less effective than that under large-scale low-speed streaks. As a result, the effectiveness of control in suppressing near-wall coherent structures and Reynolds shear stresses becomes weaker at the higher Reynolds number. The quantitative analysis of the contributions to the drag reduction rate from outer and inner regions shows that the effectiveness of the control is mainly determined by the suppression degree of near-wall motions. Furthermore, budgets of streamwise enstrophy are analysed to reveal the interaction of large-scale motions with near-wall streamwise vorticity. The titling effect induced by large-scale motions is positive under large-scale high-speed streaks, but negative under large-scale low-speed streaks, which could be a possible way of large-scale motion to modulate streamwise vorticity. In the controlled cases, the positive titling effect induced by large-scale motions under large-scale high-speed streaks is even enhanced, while other terms in the budgets are reduced, which could explain the degradation of control effectiveness in suppressing near-wall streamwise vortices under large-scale high-speed streaks. Therefore, the loss in the drag reduction rate at the higher Reynolds number is due to the weakened control effectiveness on near-wall coherent structures, which are exposed to the modulation effect of large-scale motions.     

周期性波纹通道内换热与流动的数值研究

采用曲线坐标系下压力与速度耦合的SIMPLER算法,数值研究了一种紧凑换热器中波纹通道内周期性充分发展的层流流动与换热情况,流动Re数的范围为100~1100,Pr数为0.7.计算考察了不同波纹高度、波纹间距对流动与换热的影响,并对模型参数进行了性能评价.计算结果表明,整体Nu数及fRe数随着流动Re数的增加而增加.随着波纹高度的增加或波纹间距的减小,换热增强,特别是在高Re数下波纹高度的增加更加强化换热.最佳波纹高度和间距分别为1.15 mm和13 mm.     

交叉三角形波纹板流道在过渡流状态下的传热与阻力特性

交叉三角形波纹板流道中流动常处于过渡态。本文研究了交叉波纹板中的周期性完全发展流动及热传递。利用周期性降低几何流道的复杂性以及简化模拟对象。为了模拟这个拓扑结构中的过渡流,利用了已经验证的低雷诺数k-ε,湍流模型来说明流动中的湍流流动。得到了三维复杂计算区的温度、速度以及湍流场。计算了在恒壁温和恒热流密度两种边界条件下的摩擦系数和平均Nusselt数及其与雷诺数的关系。湍流中心从上层壁面的波纹处移向下层壁面的波纹处并逐步增强。     

Flow visualization and scanning PIV measurement of three-dimensional structure in near field of strongly buoyant jet

Abstract  

The near-field structure of strongly buoyant jet issuing from a square nozzle at low Froude and Reynolds numbers is studied by using LIF flow visualization and time-resolved scanning PIV. These experimental techniques allow the visualization of unsteady three-dimensional flow phenomenon occurring in the near-field of strongly buoyant jet. It is found that the buoyant jet is unstable to the positive buoyancy forces, which promote the inflow motion near the nozzle exit. The surrounding low temperature fluid moves into the nozzle inside along the nozzle corner and mixes with the high temperature fluid deep into the nozzle. Then, the flow pattern inside the nozzle becomes highly complex to promote the laminar to turbulent transition of the jet. The statistical flow characteristics of the strongly buoyant jet are evaluated from the scanning PIV measurement, and the result indicates the presence of axisymmetric distributions of mean flow and velocity fluctuations in the circle of diameter equal to the square side of the nozzle.     

A maximum-entropy principle for two-dimensional perfect fluid dynamics

We use Kullback entropy for Young measures to define statistical equilibrium states for a two-dimensional incompressible flow of a perfect fluid. This approach is justified, as it gives a concentration property about the equilibrium state in the phase space. It might give a statistical understanding of the appearance of coherent structures in two-dimensional turbulence.     

The distribution of headways in traffic noise simulation

In recent years researchers in jet turbulence and jet noise have become increasingly interested in what is termed “large scale coherent jet structures”. There is now considerable evidence that azimuthally coherent structures can be generated by acoustically forcing a jet from upstream. However, the evidence for such structures in unforced jets, except close to the nozzle at low Reynolds numbers, is, at best, circumstantial. The role of such structures in subsonic jet noise production is also completely unproven. In an attempt to establish a link between azimuthally coherent structures and the jet noise field a number of experimenters have recently made azimuthal cross-correlation measurements of either the near field pressure or far field noise, and used the observed coherence to infer the existence of an azimuthally coherent source field. The term azimuthally coherent is used here to infer that the source region is dominated by low order azimuthal components, with relatively little contribution coming from the higher azimuthal components. The purpose of this paper is to question the interpretation of that data. Specifically the sound field generated by a simple ring source with various types of azimuthal coherence is considered theoretically. It is shown that the azimuthal coherence of both the near and far field pressures is principally a function of the Helmholtz number and in many cases of practical interest is relatively insensitive to any coherent structure of the source.     

An Experimental Investigation of Forced Convection Heat Transfer from a Coiled Heat Exchanger Embedded in a Packed Bed

Forced convection heat transfer from a helically coiled heat exchanger embedded in a packed bed of spherical glass particles was investigated experimentally. With dry air at ambient pressure and temperature as a flowing fluid, the effect of particle size, helically coiled heat exchanger diameter, and position was studied for a wide range of Reynolds numbers. It was found that the particle diameter, the helically coiled heat exchanger diameter and position, and the air velocity are of great influence on the convective heat transfer between the helically coiled heat exchanger and air. Results indicated that the heat transfer coefficient increased with increasing the air velocity, increasing helically coiled heat exchanger diameter, and decreasing the particle size. The highest heat transfer coefficients were obtained with the packed-bed particle size of 16 mm and heat exchanger coil diameter of 9.525 mm (1/4 inch) at a Reynolds number range of 1,536 to 4,134 for all used coil positions in the conducted tests. A dimensionless correlation was proposed for Nusselt number as a function of Reynolds number, particle size, coil size, and coil position.     

Viscous near-wall flow in a wake of circular cylinder at moderate Reynolds numbers

Here we present the results of experimental investigation of a cross flow around a circular cylinder mounted near the wall of a channel with rectangular cross section. The experiments were carried out in the range of Reynolds numbers corresponding to the transition to turbulence in a wake of the cylinder. Flow visualization and SIV-measurements of instantaneous velocity fields were carried out. Evolution of the flow pattern behind the cylinder and formation of the regular vortex structures were analyzed. It is shown that in the case of flow around the cylinder, there is no spiral motion of fluid from the side walls of the channel towards its symmetry plane, typical of the flow around a spanwise rib located on the channel wall. The laminar-turbulent transition in the wake of the cylinder is caused by the shear layer instability.     

Decoherence of odd compass states in the phase-sensitive amplifying/dissipating environment

We study the evolution of odd compass states (specific superpositions of four coherent states), governed by the standard master equation with phase-sensitive amplifying/attenuating terms, in the presence of a Hamiltonian describing a parametric degenerate linear amplifier. Explicit expressions for the time-dependent Wigner function are obtained. The time of disappearance of the so called “sub-Planck structures” is calculated using the negative value of the Wigner function at the origin of phase space. It is shown that this value rapidly decreases during a short “conventional interference degradation time” (CIDT), which is inversely proportional to the size of quantum superposition, provided the anti-Hermitian terms in the master equation are of the same order (or stronger) as the Hermitian ones (governing the parametric amplification). The CIDT is compared with the final positivization time (FPT), when the Wigner function becomes positive. It appears that the FPT does not depend on the size of superpositions, moreover, it can be much bigger in the amplifying media than in the attenuating ones. Paradoxically, strengthening the Hamiltonian part results in decreasing the CIDT, so that the CIDT almost does not depend on the size of superpositions in the asymptotical case of very weak reservoir coupling. We also analyze the evolution of the Mandel factor, showing that for some sets of parameters this factor remains significantly negative, even when the Wigner function becomes positive.