Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

A multi-scale simulation method for high Reynolds number wall-bounded turbulent flows

We present an assessment and enhancement of the hybrid two-level large-eddy simulation method (A.G. Gungor and S. Menon, A new two-scale model for large eddy simulation of wall-bounded flows, Prog. Aerosp. Sci. 46 (2010), pp. 28–45), a multi-scale formulation for simulation of high Reynolds number wall-bounded turbulent flows. The assessment of the method is performed by examining role of static and dynamic blending functions used to perform hybridisation of two-level simulation (K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows, J. Comput. Phys. 220 (2006), pp. 290–311; K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows. Part 2: Non-homogeneous flows, J. Comput. Phys. 222 (2007), pp. 673–701) and large-eddy simulation methods. The sensitivity of first- and second-order turbulence statistics to the type of blending functions is investigated by simulating a fully developed turbulent flow in a channel at a friction Reynolds number Re_τ = 395 and comparing the results with those obtained using a direct numerical simulation. The first-order statistics do not show any significant differences for different blending functions, but the second-order statistics show some minor differences. The dynamic evaluation of the hybrid region and the blending function is necessary for non-equilibrium and complex flows where use of a static blending function can lead to inaccurate results. We propose two criteria for the dynamic evaluation; first evaluates extent of the hybrid region based on the subgrid turbulent kinetic energy and the second estimates the blending function based on a characteristic length scale. The computational efficiency of the method is enhanced by incorporating a hybrid programming paradigm where a standard domain decomposition by the message-passing-interface library is combined with the open multi-processing based parallelisation. A further enhancement of the method is achieved by incorporating a closure model for the unclosed hybrid terms in the governing equations, which appear due to hybridisation of two-level- and large-eddy-simulation methods. The model is based on an order of magnitude approximation and a preliminary assessment of the model shows improvement of turbulence statistics when used to simulate turbulent flow in a periodic channel. The assessment and improvements to the multi-scale method make it more suitable for simulation of practical wall-bounded turbulent flows at higher Reynolds number than a conventional large-eddy simulation. This is demonstrated by simulating two representative cases; turbulent flow at high Reynolds number in a periodic channel and flow over a bump placed on the lower surface of a channel, where a relatively coarser computational grid is found to be sufficient for reasonably accurate results.     

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary(IB) method has been developed to simulate laminar flows. However,its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work,we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case,we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.     

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.     

Vorticity,backscatter and counter-gradient transport predictions using two-level simulation of turbulent flows

The two-level simulation (TLS) method evolves both the large-and the small-scale fields in a two-scale approach and has shown good predictive capabilities in both isotropic and wall-bounded high Reynolds number (Re) turbulent flows in the past. Sensitivity and ability of this modelling approach to predict fundamental features (such as backscatter, counter-gradient turbulent transport, small-scale vorticity, etc.) seen in high Re turbulent flows is assessed here by using two direct numerical simulation (DNS) datasets corresponding to a forced isotropic turbulence at Taylor’s microscale-based Reynolds number Re_λ ≈ 433 and a fully developed turbulent flow in a periodic channel at friction Reynolds number Re_τ ≈ 1000. It is shown that TLS captures the dynamics of local co-/counter-gradient transport and backscatter at the requisite scales of interest. These observations are further confirmed through a posteriori investigation of the flow in a periodic channel at Re_τ = 2000. The results reveal that the TLS method can capture both the large- and the small-scale flow physics in a consistent manner, and at a reduced overall cost when compared to the estimated DNS or wall-resolved LES cost.     

The self-similarity of wall-bounded temporally accelerating turbulent flows

From the study of viscous flow it is known that certain time-dependent laminar problems, such as the impulsively started flat plate and the diffusion of a vortex sheet, possess self-similar solutions. Previous studies of turbulent channel and pipe flows accelerating between two steady states have shown that the flow field evolves in three distinct stages. Furthermore, recent direct numerical simulations have shown that the perturbation velocity, i.e. the surplus velocity from the initial value, in an impulsively accelerating turbulent channel and pipe flow also possesses a self-similar distribution during the initial stage. In here, these results are developed analytically and it is shown that accelerating flows in which the centreline velocity develops as U^∧_c(t) = U₀(t/t₀)^m will possess a self-similar velocity distribution during the initial stage. The displacement thickness of the perturbation velocity is shown to be dependent only on the type of acceleration, and not on the initial Reynolds number, the acceleration rate or the change in Reynolds number. The derived formulas are verified with good agreement against measurements performed in a linearly accelerating turbulent pipe flow and with data from channel flow simulations.     

Kolmogorov scaling of turbulent flow in the vicinity of the wall

How to scale even the simplest of turbulent flows continues to be a cause for considerable controversy. In the present research, a data base compiling results from channel flow direct numerical simulations and turbulent boundary layer experiments is employed to investigate the properties of shear and normal Reynolds stresses very close to the wall. Two types of scaling based on Kolmogorov length and velocity scales are analyzed. It is shown that it is highly likely that large length scales of the order of the channel half-width or the boundary layer thickness play an important role even in the innermost regions of wall-bounded turbulent flows, which hints at the persistence of Reynolds number effects in even high Reynolds number flows.     

A one-equation turbulence model for recirculating flows

A one-equation turbulence model which relies on the turbulent kinetic energy transport equation has been developed to predict the flow properties of the recirculating flows. The turbulent eddy-viscosity coefficient is computed from a recalibrated Bradshaw’s assumption that the constant a₁ = 0.31 is recalibrated to a function based on a set of direct numerical simulation (DNS) data. The values of dissipation of turbulent kinetic energy consist of the near-wall part and isotropic part, and the isotropic part involves the von Karman length scale as the turbulent length scale. The performance of the new model is evaluated by the results from DNS for fully developed turbulence channel flow with a wide range of Reynolds numbers. However, the computed result of the recirculating flow at the separated bubble of NACA4412 demonstrates that an increase is needed on the turbulent dissipation, and this leads to an advanced tuning on the self-adjusted function. The improved model predicts better results in both the non-equilibrium and equilibrium flows, e.g. channel flows, backward-facing step flow and hump in a channel.     

Finite element-based large eddy simulation using a combination of the variational multiscale method and the dynamic Smagorinsky model

A finite element-based large eddy simulation (LES) is proposed using a combination of the residual-based variational multiscale (RBVMS) approach and the dynamic Smagorinsky eddy-viscosity model. In this combined model, the cross-stress terms are modelled using the RBVMS approach while the eddy-viscosity model is used to represent the Reynolds stresses. The eddy-viscosity is computed dynamically in a local fashion for which a localized version of the variational Germano identity is developed. To improve the robustness of the local dynamic procedure, two types of averaging schemes are considered. The first type employs spatial averaging over homogeneous direction(s) which is only applicable to turbulent flows with statistical homogeneity in at least one direction. The second type is based on Lagrangian averaging over fluid pathtubes, which is applicable to inhomogeneous turbulent flows. The predictions from the combined model are compared to the direct numerical simulation or experimental data and also to the predictions from the RBVMS model. This is done for two cases: turbulent flow in a channel (Re_τ = 590) and flow over a cylinder (Re_D = 3, 900). For the turbulent channel flow, predictions are similar between the RBVMS model and the combined model. For flow over a cylinder, the combined model provides better predictions, specifically for fluctuations in the streamwise velocity and lift.     

Effects of spanwise system rotation on turbulent stripes in a plane Poiseuille flow

In the transitional channel flow, the large-scale intermittent structure of localised turbulence, which is called the turbulent stripe pattern, can be found in the form of stripe arrangement. The structure of the turbulent stripe pattern is an oblique laminar–turbulent banded pattern and is inclined with respect to the streamwise direction. We performed direct numerical simulation at a transitional Reynolds number and very low-rotation numbers, and focused on the turbulent stripe pattern in the plane Poiseuille flow subjected to spanwise system rotation. We captured the turbulent stripe pattern in a rotating channel flow and found the augmentation and diminution of the turbulent stripe pattern were affected by the spanwise rotation. The contents of the discussion are the spatial size of the turbulent stripe pattern on the basis of the instantaneous flow fields, the energy spectra, and various statistics relating to the spanwise velocity component that characterise the turbulent stripe pattern. The turbulent stripe pattern was found to contain kinetic energy that was larger in very weakly rotating flows than in the static system. It was also found that the magnitude of the spanwise secondary flow increases, while the quasi-laminar region is wider at a very lowrotation number.     

Micro-machined super-miniature hot-film multi-array probe for field experiments with sub-Kolmogorov resolution

This is a report on the first successful experiments with a unique micro-hot-film multi-array probe. The main aim of this part of the project was to design, fabricate, implement and test in laboratory conditions a micro-hot-film multi-array probe for measurements in turbulent flows in field experiments with sub-Kolmogorov resolution. The key feature is that the probe is enabling to access, along with all three components of turbulent velocity fluctuations, also all nine components of the spatial velocity gradients tensor, including the option of obtaining the stream-wise velocity derivatives without employing the Taylor hypothesis, as well as velocity temporal derivatives. Hence the multi-array probe (typically five arrays) is built with each array consisting of four hot-film sensors, similarly to that of larger multi-array probe, based on arrays with four hot wires used in the atmospheric surface layer at Taylor micro-scale Reynolds number up to Re_λ ≈ 10⁴. This part of work relates to Reynolds numbers Re_λ < 500 and employs several laboratory flows including flow in a channel of rectangular cross section, in its bulk and in the boundary layers in the proximity of its smooth and rough walls, and in the confined jet of circular cross section.     

A mixed-time-scale low-Reynolds-number one-equation turbulence model

In this paper, a new low-Reynolds-number (LRN) one-equation turbulence model for eddy viscosity is proposed. A mixed time scale, representing a combination of three time scales: two time scales made of strain-rate parameter S and vorticity parameter Ω and the turbulent time scale k/?, is introduced into this model. The proposed model is derived from an LRN k?? two-equation model where the mixed time scale has been proved to be very effective for predicting local flows over complex terrains. In the transport equation of the model, the mixed time scale is included in the production and the dissipation terms. The new model is evaluated in channel flows at various Reynolds numbers, boundary layer flows with or without pressure gradient and backward-facing step flows with different expansion ratios and Reynolds numbers. Then the grid convergence of the model is investigated. Finally, the model performance for different values of the weighting constant C_s in the mixed time scale is assessed. The results show that the proposed model reproduces the correct wall-limiting behaviour of turbulent quantities and performs well in the near-wall region of turbulent flows. The model could be expected to be adopted in hybrid Reynolds averaged Navier–Stokes/large eddy simulation methodology for complex wall-bounded flows at high Reynolds numbers.     

On the application of the two-level large-eddy simulation method to turbulent free-shear and wake flows

We present application of the hybrid two-level large-eddy simulation (TLS-LES) method, a multi-scale simulation model, to turbulent free-shear and wake flows at moderately high Reynolds number. The TLS-LES method combines the scale-separation-based two-level simulation (TLS) model with the spatial-filtering-based conventional large-eddy simulation (LES) model in an additive manner using a normalised blending function. The additive blending can be performed in a static or a dynamic manner. We demonstrate that the method, which has been originally developed for wall-bounded flows, can be used to simulate flows in complex configurations without requiring any further adjustments to the model. In this study, three canonical flows are simulated, which are representative of free-shear and wake flows. These cases include a temporally evolving mixing layer, flow past a circular cylinder in a uniform flow and flow past a finite-span airfoil placed in a uniform flow at three different angle of attacks. We analyse the role of static and dynamic blending functions, large-scale grid resolution and the effect of small scales on the instantaneous flow features and turbulence statistics. The results obtained from these cases demonstrate robustness, accuracy and consistency of the multi-scale TLS-LES method and show that the method is suitable for investigation of turbulent flows that encompass features such as massive separation, reattachment, transition to turbulence and unsteady wake, which are challenging to model numerically.     

Fractal modelling of turbulent mixing

The aim of this work is to propose a new model for turbulent flows, called the fractal model (FM), applicable both in a Reynolds averaged Navier–Stokes (RANS) and a large-eddy simulation (LES) formulation, with the ultimate goal of applying it to simulate turbulent combustion irrelevant of its mode (premixed or non-premixed). The model is able to turn itself off in the laminar zones of the flow, and in particular near walls. It is based on the fractal theory. It describes the physics of the smaller spatial scales and therefore represents a small-scales model.

FM describes the physics of the small scales of turbulence based on the phenomenological concept of vortex cascade and on the self-similar behaviour of turbulence in the inertial range. Such a model is used in each cell of a numerical calculation. A characteristic length Δ is associated to each cell, and the local energy u ³ _Δ/Δ is distributed over a certain number of eddies, which depends on the local Reynolds number Re _Δ. Each vortex of the cascade generates N _c vortices; the recursive process of vortex generation terminates at the dissipative scale level, i.e. when the eddy Reynolds number is equal to one. FM is also able to estimate the volume fraction occupied by the dissipative fine structures of turbulence; this quantity is critical in reactive turbulent flows.

The physics of small scales is summarized by a turbulent ‘viscosity’ μ_t, to be added to the molecular one. μ_t is zero where the flow is laminar and, in particular, goes to zero at solid walls. Assuming μ_t to be isotropic, FM is applicable in a RANS formulation (IFM, isotropic fractal model). The model can be extended to the anisotropic case (AFM, anisotropic fractal model) and therefore used to close the transport equations in an LES approach. In the present paper, the model (IFM) is used in a RANS approach and is validated through a test case studied experimentally by Johnson and Bennett, and numerically (with LES) by Akselvoll and Moin. The results obtained are in good agreement both with the experimental and the numerical ones. Other tests are being performed.     

On mean flow universality of turbulent wall flows. II. Asymptotic flow analysis

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence statistics, their universality, the convergence of statistics towards their asymptotic profiles, and the overall physical flow organization. In extension of recent studies focusing on the mean flow at moderate and relatively high Re, the latter questions are addressed with respect to three canonical wall-bounded flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Main results reported here are the asymptotic logarithmic law for the mean velocity and corresponding scale-separation laws for bulk flow properties, the Reynolds shear stress, the turbulence production and turbulent viscosity. A scaling analysis indicates that the establishment of a self-similar turbulence state is the condition for the development of a strict logarithmic velocity profile. The resulting overall physical flow structure at extreme Re is discussed.     

Zonal PANS: evaluation of different treatments of the RANS–LES interface

The partially Reynolds-averaged Navier–Stokes (PANS) model can be used to simulate turbulent flows either as RANS, large eddy simulation (LES) or DNS. Its main parameter is f_k whose physical meaning is the ratio of the modelled to the total turbulent kinetic energy. In RANS f_k = 1, in DNS f_k = 0 and in LES f_k takes values between 0 and 1. Three different ways of prescribing f_k are evaluated for decaying grid turbulence and fully developed channel flow: f_k = 0.4, f_k = k^3/2 _tot/? and, from its definition, f_k = k/k_tot where k_tot is the sum of the modelled, k, and resolved, k_res, turbulent kinetic energy. It is found that the f_k = 0.4 gives the best results. In Girimaji and Wallin, a method was proposed to include the effect of the gradient of f_k. This approach is used at RANS– LES interface in the present study. Four different interface models are evaluated in fully developed channel flow and embedded LES of channel flow: in both cases, PANS is used as a zonal model with f_k = 1 in the unsteady RANS (URANS) region and f_k = 0.4 in the LES region. In fully developed channel flow, the RANS– LES interface is parallel to the wall (horizontal) and in embedded LES, it is parallel to the inlet (vertical). The importance of the location of the horizontal interface in fully developed channel flow is also investigated. It is found that the location – and the choice of the treatment at the interface – may be critical at low Reynolds number or if the interface is placed too close to the wall. The reason is that the modelled turbulent shear stress at the interface is large and hence the relative strength of the resolved turbulence is small. In RANS, the turbulent viscosity – and consequently also the modelled Reynolds shear stress – is only weakly dependent on Reynolds number. It is found in the present work that it also applies in the URANS region.     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

Direct numerical simulation of a fully developed compressible wall turbulence over a wavy wall

Compressible turbulent channel flow over a wavy surface is investigated by direct numerical simulations using high-resolution finite difference schemes. The Reynolds number considered in the present paper is 3380 based on the bulk velocity, the channel half-width and the kinetic viscosity at the wall. Four test cases are simulated and analysed at Ma_m = 0.33, 0.8, 1.2, 1.5 based on the bulk velocity and the speed of sound at the wall. We mainly focus on the curvature and the Mach number effects on the compressible turbulent flows. Numerical results show that although the wavy wall has effects on the mean and fluctuation quantities, log law still exists in the distribution of the wave-averaged streamwise velocity if the roughness effects are taken into consideration in the scaling of it. Near-wall streaks are broken by the wavy surface and near-wall quasi-streamwise vortices mostly begin at the upslope of the wave and pass over the crest of it. The wavy wall makes the turbulence more active and the flow easier to be blended. From the viewpoint of turbulent kinetic budgets, curvature effects strengthen both the diffusion terms and the dissipation terms. At the same time, they change the properties of the compressibility-related terms and promote more inner energy transferring into turbulent kinetic energy. As the Mach number increases, the reattachment of the mean flow is delayed, which indicates the mean separation bubble becomes larger. Concerning the near-wall coherent structures, the vortices are more sparsely distributed with the increasing of the Mach number. For the supersonic cases, shock waves appear. Though they have little effects on the mean turbulent quantities, they change the structures of the flow fields and induce local separations at the upper wall of the channel.     

Effects of inclined impermeable bed on the turbulent characteristics of the flow using particle image velocimetry

In this study, the turbulent characteristics of the flow in an open channel with horizontal and inclined impermeable beds were studied experimentally using two-dimensional particle image velocimetry (PIV). The experiments were conducted in a channel of 6.5 m in length, 7.5 cm in width and 25 cm in height. The slope of the channel was S = 0 for the horizontal impermeable bed and for S = ?0.002, S = ?0.004 and S = ?0.006 for the inclined impermeable bed. Hydraulic characteristics such as distributions of velocities, turbulent intensities and Reynolds stress are investigated at a fine resolution using the PIV. Velocity is measured above the horizontal and inclined impermeable bed for the same different heights (h = 5, 7, 9, 11 and 13 cm) and for the same different discharges (Q = 0.735, 0.845 and 0.970 lt/s). Results show that the channel slope influences significantly near the impermeable bed but not near the free surface the variation of turbulent characteristics of the flow and also the alteration of the channel slope from ?0.002 to ?0.006 doesn't influence the variation of turbulent characteristics of the flow, which are the longitudinal turbulent intensity u′U_*, the vertical turbulent intensity v′/U_* and the turbulent kinetic energy. The channel slope doesn't influence the Reynolds stress.