Analysis of pipe flow transition. Part I. Direct numerical simulation

We have developed an accurate hybrid finite-difference code for the simulation of unsteady incompressible pipe flow. The numerical scheme uses compact finite differences of at least eighth-order accuracy for the axial coordinate, and Chebyshev and Fourier polynomials for the radial and azimuthal coordinates, respectively. Boundary conditions for the incompressible flow are enforced using an influence-matrix technique, and the Poisson equation for pressure is solved using a fast direct method. The code has been used to simulate and analyze the spatial transition process in developed laminar pipe flow at a Reynolds number of Re=2350. Results of the simulation are compared to experimental data given by Han, Tumin and Wygnanski [18]. PACS 47.11.+j, 47.20.Ft, 47.27.Cn     

NUMERICAL SIMULATION OF COMPLICATED PIPE FLOW WITH k-εMODEL

In this paper, turbulence in a complicated pipe is simulated by using the k-ε model. The ladder-like mesh approximation is used to solve the problem of complicated boundary with the result of numerical simulation favorable. Two computational examples are given to validate the strong adaptability and stability of k-ε model.     

Numerical simulation of boundary-layer transition at Mach two

The linear and early nonlinear stages of boundary-layer transition at free-stream Mach numberM ==2.0 are investigated by direct numerical simulation of the compressible Navier-Stokes equations. Results from simulations with a large computational box and small-amplitude random initial conditions are compared with linear stability theory. The growth rates of oblique waves are reproduced correctly. Two-dimensional waves show a growth that is modulated in time, indicating the presence of an extra unstable mode which moves supersonically relative to the free stream. Further simulations are conducted to investigate the nonlinear development of two- and three-dimensional disturbances The transition due to oblique disturbance waves is the most likely cause of transition at this Mach number, and is found to lead to the development of strong streamwise vortices.     

Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

Very large‐scale computation of very small‐scale motions in a fully developed pipe flow

Local isotropy theory is examined using direct numerical simulation in a fully developed pipe flow at two Reynolds numbers Re_τ=1285.6 and 684.8. The approach to local isotropy is assessed with reference to the two Kolmogorov classical equations for longitudinal and transverse velocity structure functions. The results for the second‐order longitudinal structure functions in both the dissipative and inertial ranges indicate an improved agreement with the local isotropy hypothesis as the centreline is approached. However, the transverse structure functions satisfy isotropy neither in the dissipative range or in the inertial range. The distribution of the longitudinal and transverse structure functions also shows a substantial Reynolds number dependance in the logarithmic region of the flow and beyond. The results for the third‐order longitudinal structure function demonstrate an increased Reynolds number influence, and a deteriorating tendency to local isotropy for large separations. Contour images of axial velocity differences in the dissipative and inertial ranges have exhibited interesting patterns in relation to those of the instantaneous axial velocity. Finally, the results obtained in this investigation are in very good agreement with other published experimental and numerical data on channel and duct flows. Copyright © 2010 John Wiley & Sons, Ltd.     

Universal threshold of the transition to localized turbulence in shear flows

Experimental and numerical studies have shown similarities between localized turbulence in channel and pipe flows. By scaling analysis of a disturbed-flow model, this paper proposes a local Reynolds number Re_M to characterize the threshold of transition triggered by finite-amplitude disturbances. The Re_M represents the maximum contribution of the basic flow to the momentum ratio between the nonlinear convection and the viscous diffusion. The lower critical Re_M observed in experiments of plane Poiseuille flow, pipe Poiseuille flow and plane Couette flow are all close to 323, indicating the uniformity of mechanism governing the transition to localized turbulence.     

Subharmonic transition to turbulence in channel flow

A direct numerical simulation of subharmonic transition to turbulence in channel flow has been performed. The stages of primary and secondary instability have been identified in the results leading to a staggered pattern of A-shaped vortices. The associated staggered symmetry has been found to persist up to rather late stages of the breakdown process. This symmetry rapidly disappears in the final stage of transition and a developed turbulent flow is attained. A pronounced difference in the development between the two channel halves is observed which is consistent with the predictions of secondary instability theory.     

Research on the motion of particles in the turbulent pipe flow of fiber suspensions

The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence. The orientation distribution of fibers in the computational area for different Re numbers was computed. The results which were consistent with the experimental ones show that the fluctuation velocity of turbulence cause fibers to orient randomly. The orientation distributions become broader as the Re number increases. Then the fluctuation velocity and angular velocity of fibers were obtained. Both are affected by the fluctuation velocity of turbulence. The fluctuation velocity intensity of fiber is stronger at longitudinal than at lateral, while it was opposite for the fluctuation angular velocity intensity of fibers. Finally, the spatial distribution of fiber was given. It is obvious that the fiber dispersion is strenghened with the increase of Re numbers.     

Laminar-turbulent transition in the boundary layer on cones in a hypersonic flow at high Reynolds numbers per meter

The laminar-turbulent transition is experimentally studied in boundary-layer flows on cones with a rectangular axisymmetric step in the base part of the cone and without the step. The experiments are performed in an A-1 two-step piston-driven gas-dynamic facility with adiabatic compression of the working gas with Mach numbers at the nozzle exit M _∞ = 12–14 and pressures in the settling chamber P₀ = 60–600 MPa. These values of parameters allow obtaining Reynolds numbers per meter near the cone surface equal to Re _1e = (53–200) · 10⁶ m ⁻¹. The transition occurs at Reynolds numbers Re _tr = (2.3–5.7) · 10⁶. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 76–83, May–June, 2007.     

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

An experimental investigation of heat transfer to pulsating pipe air flow with different amplitudes

Heat transfer characteristics to both laminar and turbulent pulsating pipe flows under different conditions of Reynolds number, pulsation frequency, pulsator location and tube diameter were experimentally investigated. The tube wall of uniform heat flux condition was considered for both cases. Reynolds number varied from 750 to 12,320 while the frequency of pulsation ranged from 1 to 10 Hz. With locating the pulsator upstream of the inlet of the test section tube, results showed an increase in heat transfer rate due to pulsation by as much as 30% with flow Reynolds number of 1,643 and pulsation frequency of 1 Hz, depending on the upstream location of the pulsator valve. Closer the valve to the tested section inlet, the better improvement in the heat transfer coefficient is achieved. Upon comparing the heat transfer results of the upstream and the downstream pulsation, at Reynolds number of 1,366 and 1,643, low values of the relative mean Nusselt number were obtained with the upstream pulsation. Comparing the heat transfer results of the two studied test sections tubes for Reynolds number range from 8,000 to 12,000 and pulsation frequency range from 1.0 to 10 Hz showed that more improvement in heat transfer rate was observed with a larger tube diameter. For Reynolds number ranging from 8,000 to 12,000 and pulsation frequency of 10 Hz, an improvement in the relative mean Nusselt number of about 50% was obtained at Reynolds number of 8,000 for the large test section diameter of 50 mm. While, for the small test section diameter of 15 mm, at same conditions of Reynolds number and frequency, a reduction in the relative mean Nusselt number of up to 10% was obtained.     

内管自转且公转时环管幂律流的受力分析

数值模拟了环管中内管偏心自转且公转时由轴向压力所驱动的幂律流体充分发展层流,分析了内管上的流体作用力。结果表明,内管偏心自转时流体作用力具有推动内管作和自转同向公转的效果。当只有外力矩驱动内管自转时,由于流体的作用,随内管线密度的不同,内管能达到的受力平衡态也不同:线密度较小时内管仅能在同心自转时达到受力平衡;线密度较大时内管能在作具有不变角速度和偏心率公转时达到受力平衡,且内管线密度越大,对应的受力平衡的公转的偏心率也越大。     

Delay of finite-span excrescence-induced transition using plasma-based control

Direct numerical simulations are carried out to explore the use of flow control that delays transition generated by excrescence on a plate-like geometry in subsonic flow. Both forward-facing and rearward-facing steps of small roughness heights are considered in the investigation. These are representative of joints and other surface imperfections on wing sections that disrupt laminar flow, thereby increasing skin friction and configuration drag. Unlike previous studies, the steps have a finite lateral extent, such that sharp edges occur in both the spanwise and streamwise directions, and provide a more realistic characterisation of misaligned panels in aerodynamic configurations. The effect of spanwise corners upon transition is examined, and dielectric barrier discharge plasma-based flow control is applied to delay transition and increase the extent of the laminar flow region. Solutions are obtained to the Navier– Stokes equations that were augmented by source terms used to represent body forces imparted by plasma actuators on the fluid. A simple phenomenological model provided these forces resulting from the electric field generated by the plasma. The numerical method is based upon a high-fidelity scheme and an implicit time-marching approach, on an overset mesh system that is used to represent the finite-span steps. Very small-amplitude numerical forcing is employed to generate perturbations, which are amplified by the geometric disturbances and result in transition, similar to the physical situation. Both continuous and pulsed operations of actuators are considered, and the effectiveness of the control is quantified. Transition with the forward-facing step is considerably exacerbated by the presence of a spanwise edge. Plasma control is minimally effective, even with the use of multiple actuators and increased applied force. For the rearward-facing step, transition is substantially delayed by plasma control with small force application.     

Direct simulation of breakdown to turbulence following oblique instability waves in a supersonic boundary layer

The late stages of transition to turbulence in a Mach two boundary layer are investigated by direct numerical simulation of the compressible Navier-Stokes equations. The primary instability at this Mach number consists of oblique waves, which are known to form a pattern of quasi-streamwise vortices. It is found that breakdown does not follow immediately from these vortices, which decay in intensity. The generation of new vortices is observed by following the evolution of the pressure and vorticity in the simulation, and analysed by consideration of vorticity stretching. It is found that the slight inclined and skewed nature of the quasi-streamwise vortices leads to a production of oppositely signed streamwise vorticity, which serves as a strong localised forcing of the shear layer alongside the original vortices, formed by convection and stretching of spanwise vorticity. The shear layer rolls up into many new vortices, and is followed by a sharp increase in the energy of higher frequencies and in the skin friction.     

Heat transfer to air–water annular flow in a horizontal pipe

Two-phase air–water flow and heat transfer in a 25 mm internal diameter horizontal pipe were investigated experimentally. The water superficial velocity varied from 24.2 m/s to 41.5 m/s and the air superficial velocity varied from 0.02 m/s to 0.09 m/s. The aim of the study was to determine the heat transfer coefficient and its connection to flow pattern and liquid film thickness. The flow patterns were visualized using a high speed video camera, and the film thickness was measured by the conductive tomography technique. The heat transfer coefficient was calculated from the temperature measurements using the infrared thermography method. It was found that the heat transfer coefficient at the bottom of the pipe is up to three times higher than that at the top, and becomes more uniform around the pipe for higher air flow-rates. Correlations on local and average Nusselt number were obtained and compared to results reported in the literature. The behavior of local heat transfer coefficient was analyzed and the role of film thickness and flow pattern was clarified.     

Effect of liquid viscosity on the stratified-slug transition in horizontal pipe flow

The effect of liquid viscosity on the initiation of slug flow was studied in horizontal 2.52 and 9.53 cm pipelines. The results show the stabilizing effect of viscosity predicted by Lin & Hanratty, and are at variance with analyses which use a long-wavelength inviscid approximation. For very viscous liquids a stability analysis which recognizes that slugs originate from a train of small-wavelength sinusoidal waves seems consistent with the measurements.     

Numerical solutions of pulsating flow and heat transfer characteristics in a pipe

Numerical studies are made of flow and heat transfer characteristics of a pulsating flow in a pipe. Complete time-dependent laminar boundary-layer equations are solved numerically over broad ranges of the parameter spaces, i.e., the frequency parameter β and the amplitude of oscillation A. Recently developed numerical solution procedures for unsteady boundary-layer equations are utilized. The capabilities of the present numerical model are satisfactorily tested by comparing the instantaenous axial velocities with the existing data in various parameters. The time-mean axial velocity profiles are substantially unaffected by the changes in β and A. For high frequencies, the prominent effect of pulsations is felt principally in a thin layer near the solid wall. Skin friction is generally greateer than that of a steady flow. The influence of oscillation on skin friction is appreciable both in terms of magnitude and phase relation. Numerical results for temperature are analyzed to reveal significant heat transfer characteristics. In the downstream fully established region, the Nusselt number either increases or decreases over the steady-flow value, depending on the frequency parameter, although the deviations from the steady values are rather small in magnitude for the parameter ranges computed. The Nusselt number trend is amplified as A increases and when the Prandtl number is low below unity. These heat transfer characteristics are qualitatively consistent with previous theoretical predictions.     

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.