Guidelines for Modeling a 2D Rough Wall Channel Flow

The effects of 2-D roughness elements on the Reynolds stress anisotropy tensor, and of the energy dissipation rate anisotropy tensor of a turbulent channel flow are investigated using data obtained from direct numerical simulations (DNSs). The roughness elements consist of transverse square rods of size , placed on one wall of the channel only. While is kept constant (, is the half-width of the channel), the spacing between the rods is varied from to . The results show that the variation in can dramatically change the structure of the wall region flow. The modeling of the near-wall region needs to reflect the structural changes caused by the variation in . On the basis of the Reynolds stress budgets, attention needs to be given to the turbulent energy and pressure diffusion terms while local isotropy may be a reasonable approximation for the energy dissipation rate, especially over a range of for which the drag is near its maximum.     

Direct Numerical Simulation of Flow in a Ribbed Channel

The issue of turbine lifetime is an important one, particularly for modern turbines operating at high temperature regimes. A cooling design such as ribs may achieve an improved lifetime and complex mechanisms of heat transfer need to be well studied. In this paper, a Direct Numerical Simulation (DNS) is presented for a 3-D channel flow with two square ribs on the lower wall. The full unsteady compressible Navier-Stokes equations are solved with an original hybrid finite difference/finite element scheme. The Reynolds number of the simulation is 7 000 based on the bulk velocity at the inlet and the channel height. The present study is mainly devoted to understand the mechanism of heat transfer at the wall through the topological analysis of the flow and the temperature flux. Results show that the large-scale structures generated by obstacles splash onto the lower surface and induce longitudinal vortices which enhance heat transfer at the wall. A comprehensive data base including 56 correlations was set up for testing and improving turbulence models for this complex, separated flow.     

Passive Scalar in a Turbulent Channel Flow with Wall Velocity Disturbances

Direct Numerical Simulations (DNS) of a passive scalar in a turbulent channel flow with a normal velocity disturbance on the lower wall are presented for high and low Reynolds numbers. The aim is to reproduce the complex physics of turbulent rough flows without dealing with the geometric complexity. In addition, isothermal walls that cannot be easily assigned in an experiment, are considered. The paper explains the increase of heat transfer through the changes of the velocity and thermal structures. As in real rough flows, the transpiration produces an isotropization of the turbulence near the wall.     

Numerical Simulation of a Turbulent Flow in a Channel with Surface Mounted Cubes

In this paper we report on a fourth-order, spectro-consistent simulation of a complex turbulent flow. A spatial discretization of a convection-diffusion equation is termed spectro-consistent if the spectral properties of the convective and diffusive operators are preserved, i.e. convection skew-symmetric; diffusion symmetric positive definite. We consider a fully developed flow in a channel, where a matrix of cubes is placed at a wall of the channel. The Reynolds number (based on the channel width and the mean bulk velocity) is equal to Re = 13,000. The three-dimensional flow around the surface mounted cubes has served at a test case at the 6th ERCOFTAC/IAHR/COST workshop on refined flow modeling (Delft, June 1997). Here, mean velocity profiles as well as Reynolds stresses at various locations in the channel have been computed without using any turbulence models. The results agree well with the available experimental data.     

Temporal Evolution of Vortical Structures in the Wall Region of Turbulent Channel Flow

Hairpin-like vortical structures that form in the wall region of turbulent channel flow are investigated. The analysis is performed by following a procedure in which the Navier-Stokes equations are first integrated by means of a computational code based on a mixed spectral-finite difference technique in the case of the flow in a plane channel. A DNS turbulent-flow database, representing the turbulent statistically steady state of the velocity field through 10 viscous time units, is computed and the vortex-detection method of the imaginary part of the complex eigenvalue pair of the velocity-gradient tensor is applied to the velocity field. As a result, hairpin-like vortical structures are educed. Flow visualizations are provided of the processes of evolution that characterize hairpin vortices in the wall region of turbulent channel flow. The relationship is investigated between vortex dynamics and 2nd- and 4th- quadrant events, showing that ejections and sweeps play a fundamental role in the way the morphological evolution of a hairpin vortex develops with time.     

Direct Numerical Simulation of an Accelerating Channel Flow

Direct Numerical Simulation of a linearly accelerating channel flow starting from an initially statistically steady turbulent flow has been performed. It is shown that the response of the accelerating flow is fundamentally the same as that of the step-change transient flow described in He and Seddighi (J Fluid Mech 715:60–102, 2013). The flow structure again behaves like a boundary layer bypass transition undergoing three distinct phases, namely, (i) initially (pre-transition), the flow is laminar-like and the pre-existing turbulent structures are modulated resulting in elongated streaks leading to a strong and continuous increase in the streamwise fluctuating velocity but little changes in the other two components; (ii) it then undergoes transition when isolated turbulent spots are generated which spread and merge with each other, and (iii) they eventually cover the entire surface of the wall when the flow is fully turbulent. The similarity between the turbulence responses in the two flows is significant noting the contrasting features of the two types of mean flow unsteadiness: in the step-change flow, a sharp boundary layer is resulted in nearly instantly on the wall which closely resembles the spatially developing boundary layer, whereas the linear flow acceleration causes a continuing change of velocity gradient adjacent to the wall which propagates into the flow field with time, resulting in a gradually-developing boundary layer. There are, however, quantitative differences in the detailed behavior of the two flows and especially the transition is much delayed in the accelerating flow. It is also shown that the late pre-transition and early transition stages in both flows are characterised by significantly increased inwards sweep events in the wall region and ejection events in the outer layer. The flatness of the wall-normal velocity increases markedly near the wall around the time of onset of transition as a consequence of the huge intermittency of the velocity fluctuations. That is, there are long periods of quiescent flow coupled with occasional turbulent bursts.     

Numerical Modeling of Two-Dimensional Flow of a Nonhomogeneous Fluid in a Confined Domain

The pattern of the two-dimensional vortex flow of a nonhomogeneous fluid in a confined domain is studied using two-dimensional numerical calculations. It is found that in the case of a nonhomogeneous initial density distribution the kinetic energy decay rates are proportional to the square root of viscosity at the active stage of flow restructuring. The correlation functions of the velocity and the density are derived for different moments of time in the inertial range. All these results indicate the choice of the two-dimensional turbulence development scenario in a nonhomogeneous fluid.     

Numerical Simulation of Turbulent Flow and Heat Transfer in a Three-Dimensional Channel Coupled with Flow Through Porous Structures

This study investigates numerically the turbulent flow and heat transfer characteristics of a T-junction mixing, where a porous media flow is vertically discharged in a 3D fully developed channel flow. The fluid equations for the porous medium are solved in a pore structure level using an Speziale, Sarkar and Gatski turbulence model and validated with open literature data. Overall, two types of porous structures, consisted of square pores, are investigated over a wide range of Reynolds numbers: an in-line and a staggered pore structure arrangement. The flow patterns, including the reattachment length in the channel, the velocity field inside the porous medium as well as the fluctuation velocity at the interface, are found to be strongly affected by the velocity ratio between the transversely interacting flow streams. In addition, the heat transfer examination of the flow domain reveals that the temperature distribution in the porous structure is more uniform for the staggered array. The local heat transfer distributions inside the porous structure are also studied, and the general heat transfer rates are correlated in terms of area-averaged Nusselt number accounting for the effects of Reynolds number, velocity ratio as well as the geometrical arrangement of the porous structures.     

Direct Numerical Simulation of a Fully Developed Turbulent Flow over a Wavy Wall

Turbulent flow over a sinusoidal solid wavy surface was investigated by a direct numerical simulation using a spectral element technique. The train of waves has an amplitude to wavelength ratio of 0.05. For the flow conditions (Re=hU _b/2ν= 3460) considered, adverse pressure gradients were large enough to cause flow separation. Numerical results compare favorably with those of Hudson's (1993) measurements. Instantaneous flow fields show a large variation of the flow pattern in the spanwise direction in the separated bubble at a given time. A surprising result is the discovery of occasional velocity bursts which originate in the separated region and extend over large distances away from the wavy wall. Turbulence in this region is very different from that near a flat wall in that it is associated with a shear layer which is formed by flow separation. Received 17 April 1996 and accepted 19 November 1997     

Numerical Modeling of a Viscous-Gas Flow in a Turbine Cascade with Air Ejection

The gas flow pattern, the local friction coefficients, the profile losses, and the outlet flow angles in a plane turbine cascade are modeled numerically using the 2-D Reynolds equations. For describing the turbulence characteristics, a two-parameter q- turbulence model is used. The flow pattern behind the cascade trailing edge is studied. The calculated values of the local gas pressure and the friction coefficient on the profile contour, the profile losses, and the outlet flow angles are compared with the experimental data for a transonic flow in a nozzle cascade at various reduced gas velocities behind the cascade and relative mass flows of the air ejected.     

Steady Vortex Ring with Isochronous Flow in the Vortex Core

Steady-state solutions to the problem of a thin vortex ring in an inviscid incompressible fluid in infinite space are investigated. The Fraenkel procedure is used to construct the steady-state solutions. In this procedure a given vorticity distribution in plane flow with circular streamlines is transformed into a steady vortex ring using an expansion in the ring thinness parameter. For example, a two-dimensional vortex of constant vorticity is transformed into a steady vortex ring with the uniform distribution in which the absolute value of vorticity is proportional to the distance from the axis of symmetry. The principal aim of our study is to construct the algorithm of finding the flow for an isochronous vortex ring in which the periods of revolution are the same for all the liquid particles in the vortex core. The problem is that the two-dimensional distribution which goes over in the isochronous ring in accordance with the Fraenkel procedure is unknown in advance. In particular, the ring with the uniform distribution is not isochronous despite the isochronism of the initial two-dimensional flow. In this connection the Fraenkel procedure is significantly modified so that the initial two-dimensional vorticity distribution is determined in each of the steps of the iteration procedure. The solution for the vortex ring with the uniform distribution obtained in the present study is significantly used to construct the isochronous solution. The necessary corrections to the former solution are calculated in each step. Obtaining of the isochronous flow is the key step for the investigation of stability of three-dimensional oscillations of the vortex ring since the oscillation spectrum of this flow is discrete.     

Thermal Instability in a Horizontal Porous Channel with Horizontal Through Flow and Symmetric Wall Heat Fluxes

The linear thermoconvective instability of the basic parallel flow in a plane and horizontal porous channel is investigated. The boundary walls are assumed to be impermeable and subject to symmetric and uniform heat fluxes. The wall heat fluxes produce either a net heating or a net cooling of the fluid saturated porous medium. A horizontal mass flow rate is externally impressed leading to a stationary basic state with a temperature gradient inclined to the vertical. A region of possibly unstable thermal stratification exists either in the lower half-channel (boundary heating), or in the upper half-channel (boundary cooling). The convective instability of the basic flow is governed by the Rayleigh number and by the Péclet number. In the case of boundary heating, the thermal instability arises when the Rayleigh number exceeds its critical value, that depends on the Péclet number. The change of the critical Rayleigh number as a function of the Péclet number is determined numerically for arbitrary normal modes oblique to the basic flow direction. The most dangerous modes are the longitudinal rolls, with a wave vector perpendicular to the basic velocity. There exists a minimum value of the Péclet number, 19.1971, below which no linear instability is detected.     

床面固体颗粒圆柱绕流的试验现象观察

本文介绍了床面固体颗粒随水流绕过圆柱体时，将在圆柱周围的床面上形成一个无粒子运动区的试验现象。水槽试验结果表明，当固体颗粒的粒径减小时，无粒子区的范围将增大；无粒子区的范围随圆柱直径的增大而增大；水流条件的变化直接影响着床面固体颗粒的运动情况，同无粒子区的形成、消失及范围大小有密切的关系。根据试验资料，结合量纲分析，建立了无粒子区的无量纲经验关系式     

Large-Scale Streamwise Vortices in Turbulent Channel Flow Induced by Active Wall Actuations

Direct numerical simulations of turbulent flow in a plane channel using spanwise alternatively distributed strips (SADS) are performed to investigate the characteristics of large-scale streamwise vortices (LSSVs) induced by small-scale active wall actuations, and their role in suppressing flow separation. SADS control is obtained by alternatively applying out-of-phase control (OPC) and in-phase control (IPC) to the wall-normal velocity component of the lower channel wall, in the spanwise direction. Besides the non-controlled channel flow simulated as a reference, four controlled cases with 1, 2, 3 and 4 pairs of OPC/IPC strips are studied at M =?0.2 and R e =?6,000, based on the bulk velocity and the channel half height. The case with 2 pairs of strips, whose width is Δz ⁺ =?264 based on the friction velocity of the non-controlled case, is the most effective in terms of generating large-scale motions. It is also found that the OPC (resp. IPC) strips suppress (resp. enhance) the coherent structures and that leads to the creation of a vertical shear layer, which is responsible for the LSSVs presence. They are in a statistically steady state and their cores are located between two neighbouring OPC and IPC strips. These motions contribute significantly to the momentum transport in the wall-normal and spanwise directions showing potential for flow separation suppression.     

Modeling the Interaction of a Vortex Pair with the Ground

The results of a numerical investigation of the interaction between an aircraft vortex wake (vortex pair) and the ground during the takeoff and landing phases are presented. The calculations were performed within the framework of the time-dependent two-dimensional Reynolds-averaged Navier-Stokes equations using the Spalart-Shur turbulence model generalizing the turbulent viscosity transport model of Spalart and Allmaras to the case of the flows with curved streamlines and rotation. Similar calculations were carried out on the basis of the original Spalart-Allmaras model and the k- model of Menter. Some new qualitative and quantitative data on the distinctive features of the phenomenon under consideration are obtained.     

Modeling the Turbulent Vortex Wake Behind a High-Lift Wing

A mathematical model and the results of a numerical simulation of the multiple vortex structure behind the high-lift wing (with deflected flaps) of a civil aircraft model are presented. The calculations are performed within the framework of the simplified (with allowance for the specific features of the flow under consideration) three-dimensional Reynolds equations using three turbulence models. The results obtained are compared with the experimental data; a considerable influence of the streamline curvature, typical of vortex flows, not only on the turbulence parameters but also on the parameters of the average flow is demonstrated. As a result, only one of the turbulence models considered, namely, that taking the streamline curvature and rotation effects on the turbulence into account, ensures acceptable accuracy of the main wake characteristics.     

Axisymmetric Vortex Fluid Flow in a Long Elastic Tube

A mathematical model for axisymmetric eddy motion of a perfect incompressible fluid in a long tube with thin elastic walls is proposed. Necessary and sufficient conditions for hyperbolicity of the system of equations of motion for flows with monotonic radial velocity profiles are formulated. The propagation velocities of the characteristics of the system under study and the characteristic shape of this system are calculated. The existence of simple waves continuously attached to a given steady shear flow is proved. The group of transformations admitted by the system is found, and submodels that determine invariant solutions are given. By integrating factorsystems, new classes of exact solutions of equations of motion are found.     

Vortex in a Supersonic Flow and its Influence on Blunt Body Flow and Heat Transfer

The dissipation of a tip vortex generated by a rectangular wing and the stagnation pressure and temperature distributions within the vortex are experimentally investigated. The interaction of the vortex with the bow shock ahead of a sphere and the heat flux distribution over the spherical surface are studied. The experiments were carried out at the Mach number M=3 and the unit Reynolds numbers Re₁=1.1·10⁷ and 3.7·10⁷ 1/m.