Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Turbulent flows in channels with injection. results of large eddy simulation and the two-equation turbulence model

Turbulent flows in channels with intense distributed injection are modeled using the large eddy method and the two-equation k-? turbulence model. The calculations are carried out for different velocities of injection from the channel walls. For a channel with one-sided injection the results of large eddy simulation are in good agreement with the measured data, whereas the calculations in accordance with the k-? model give a less convex cross-sectional velocity profile and an appreciable error in determining the surface friction coefficient on the impermeable wall and also have certain other shortcomings. In the case of two-sided injection, the results of the calculations by the large eddy method and the k-? model are in good agreement with one another and the data of physical experiments.     

Brinkman Flow of a Viscous Fluid Through a Spherical Porous Medium Embedded in Another Porous Medium

An analytical investigation for a two-dimensional steady, viscous, and incompressible flow past a permeable sphere embedded in another porous medium is presented using the Brinkman model, assuming a uniform shear flow far away from the sphere. Semi-analytical solutions of the problem are derived and relevant quantities such as velocities and shearing stresses on the surface of the sphere are obtained. The streamlines inside and outside the sphere and the radial velocity are shown in several graphs for different values of the porous parameters \({\sigma _1 =(\mu /\tilde {\mu }) (a/\sqrt{K_1 })}\) and \({\sigma _2 =(\mu /\tilde {\mu }) (a/\sqrt{K_2 })}\) , where a is the radius of the sphere, μ is the dynamic viscosity of the fluid, \({\tilde {\mu }}\) is an effective or Brinkman viscosity, while K ₁ and K ₂ are the permeabilities of the two porous media. It is shown that the dimensionless shearing stress on the sphere is periodic in nature and its absolute value increases with an increase of both porous parameters σ ₁ and σ ₂.     

Effects of the Reynolds number on a scale-similarity model of Lagrangian velocity correlations in isotropic turbulent flows

A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η≡η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η∝ R_λ~(1.39) obtained from the numerical data.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Numerical Study of Turbulent Channel Flow over Surfaces with Variable Spanwise Heterogeneities: Topographically-driven Secondary Flows Affect Outer-layer Similarity of Turbulent Length Scales

Wall-bounded turbulent flows over surfaces with spanwise heterogeneous surface roughness – that is, spanwise-adjacent patches of relatively high and low roughness – exhibit mean flow phenomena entirely different to what would otherwise exist in the absence of spanwise heterogeneity. In the outer layer, mean counter-rotating rolls occupy the depth of the flow, and are positioned such that “upwelling” and “downwelling” occurs above the low and high roughness, respectively. It has been comprehensively shown that these secondary flows are Prandtl’s secondary flow of the second kind (Anderson et al., J. Fluid Mech. 768, 316–347 2015). This behaviour indicates that spanwise spacing, s _y , between adjacent patches of high and low roughness is, itself, a problem parameter; in this study, we have systematically assessed how s _y affects turbulence structure in high Reynolds number channel flows via two-point correlations. “High roughness” is imposed with streamwise-aligned pyramid elements with height, h, selected to be ≈ 5% of the channel half height, H. For \(s_{y}/H \gtrsim 1\), we find that the aforementioned domain-scale mean circulations exist and the surface may be regarded as a topography. For s _y /H ? 0.2, turbulence statistics show characteristics very similar to a homogeneous roughness and thus the surface may be regarded as a roughness. For 0.2 ? s _y /H ? 2, the spatial extent of the counter-rotating rolls is controlled by proximity to adjacent rows, and we define such surfaces as being intermediate. We refer to such surfaces as intermediate state.     

Numerical simulation of the interaction of a transverse jet with a supersonic flow using different turbulence models

This paper presents a numerical simulation of the flow resulting from transverse jet injection into a supersonic flow through a slot nozzle at different pressures in the injected jet and the crossflow. Calculations on grids with different resolutions use the Spalart–Allmaras turbulence model, the k–ε model, the k–ω model, and the SST model. Based on a comparison of the calculated and experimental data on the wall pressure distribution, the length of the recirculation area, and the depth of jet penetration into the supersonic flow, conclusions are made on the accuracy of the calculation results for the different turbulence models and the applicability of these models to similar problems.     

Simulation of Receptivity and Induced Transition From Discrete Roughness Elements

Simulations have been carried out to predict the receptivity and growth of crossflow vortices created by Discrete Roughness Elements (DREs) The final transition to turbulence has also been examined, including the effect of DRE spacing and freestream turbulence. Measurements by Hunt and Saric (2011) of perturbation mode shape at various locations were used to validate the code in particular for the receptivity region. The WALE sub-grid stress (SGS) model was adopted for application to transitional flows, since it allows the SGS viscosity to vanish in laminar regions and in the innermost region of the boundary layer when transition begins. Simulations were carried out for two spanwise wavelengths: λ= 12mm (critical) and λ= 6mm (control) and for roughness heights (k) from 12 μm to 42 μm. The base flow considered was an ASU (67)-0315 aerofoil with 45 ⁰ sweep at -2.9 ⁰ incidence and with onset flow at a chord-based Reynolds number Re _c= 2.4x10 ⁶. For λ= 12mm results showed, in accord with the experimental data, that the disturbance amplitude growth rate was linear for k = 12 μm and 24 μm, but the growth rate was decreased for k = 36 μm Receptivity to λ= 6mm roughness showed equally good agreement with experiments, indicating that this mode disappeared after a short distance to be replaced by a critical wavelength mode. Analysis of the development of modal disturbance amplitudes with downstream distance showed regions of linear, non-linear, saturation, and secondary instability behaviour. Examination of breakdown to turbulence revealed two possible routes: the first was 2D-like transition (probably Tollmien-Schlichting waves even in the presence of crossflow vortices) when transition occurred beyond the pressure minimum; the second was a classical crossflow vortex secondary instability, leading to the formation of a turbulent wedge.     

Combining the SSG/LRR-<Emphasis Type="Italic">ω</Emphasis> differential reynolds stress model with the detached eddy and laminar-turbulent transition models

The procedure of incorporating the detached eddy method and a model of laminar-turbulent transition into the SSG/LRR-ω turbulence model is presented. The approach proposed can be regarded as the generalization of the existing models intended to perform calculations with the SST turbulence model to the case of their use with the SSG/LRR-ω model. The advantage of the approach developed over the RANS turbulence models based on the Boussinesq hypothesis is demonstrated with respect to the problems of flow past an airfoil and cold jet outflow.     

Study of the Spectral Difference Numerical Dissipation for Turbulent Flows Using Unstructured Grids

In this paper, the numerical dissipation properties of the Spectral Difference (SD) method are studied in the context of vortex dominated flows and wall-bounded turbulence, using uniform and distorted grids. First, the validity of using the SD numerical dissipation as the only source of subgrid dissipation (the so-called Implicit-LES approach) is assessed on regular grids using various polynomial degrees (namely, p = 3, p = 4, p = 5) for the Taylor-Green vortex flow configuration at R e = 5 000. It is shown that the levels of numerical dissipation greatly depend on the order of accuracy chosen and, in turn, lead to an incorrect estimation of the viscous dissipation levels. The influence of grid distortion on the numerical dissipation is then assessed in the context of finite Reynolds number freely-decaying and wall-bounded turbulence. Tests involving different amplitudes of distortion show that highly skewed grids lead to the presence of small-scale, noisy structures, emphasizing the need of explicit subgrid modeling or regularization procedures when considering coarse, high-order SD computations on unstructured grids. Under-resolved, high-order computations of the turbulent channel flow at R e _τ = 1000 using highly-skewed grids are considered as well and present a qualitatively similar agreement to results obtained on a regular grid.     

A bootstrap mechanism for non-colloidal suspension viscosity

The role of friction in non-colloidal suspensions is examined with a model which splits the viscosity into a frictionless component (τ*) plus a frictional component which depends on the ratio of the particle pressure (P) to the shear stress (τ). The model needs the input by computation of τ* and P and a suitable choice of particle friction coefficient (μ). It can be extended to elongational flows and cases where sphere roughness is important; volume fractions up to 0.5 are considered. It is shown that friction acts in a feedback or “bootstrap” manner to increase the suspension viscosity. The analysis is also useful for deducing the friction coefficient in suspensions from experimental data. It was applied to several sets of experimental data and reasonable correlations of the viscosities were demonstrated. An example of the correlation for spheres in a silicone oil is shown for volume fractions 0.1–0.5.

Open image in new window

Graphical abstract

     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 1: Linear Stability Analysis

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Reynolds-Averaged,Scale-Adaptive and Large-Eddy Simulations of Premixed Bluff-Body Combustion Using the Eddy Dissipation Concept

A lean premixed propane/air bluff-body stabilized flame (Volvo test rig) is calculated using the Scale-Adaptive Simulation turbulence model (SAS) and Large-Eddy simulations (LES) as well as the conventional Reynolds-averaged approach (RAS). RAS and SAS are closed by the standard k-?? and the k-ω Shear Stress Transport (SST) turbulence models, respectively. The conventional Smagorinsky and the k-equation sub-grid scales models are used for the LES closure. Effects of the sub-grid scalar flux modeling using the classical gradient hypothesis and Clark’s tensor diffusivity closures both for the inert and reactive LES flows are discussed. The Eddy Dissipation Concept (EDC) is used for the turbulence-chemistry interaction. It assumes that molecular mixing and the subsequent combustion occur in the ’fine structures’ (smaller dissipative eddies, which are close to the Kolmogorov scales). Assuming the full turbulence energy cascade, the characteristic length and velocity scales of the ’fine structures’ are evaluated using different turbulence models (RAS, SAS and LES). The finite-rate chemical kinetics is taken into account by treating the ’fine structures’ as constant pressure and adiabatic homogeneous reactors, calculated as a system of ordinary-differential equations (ODEs) described by a Perfectly Stirred Reactor (PSR) concept. Several further enhancements to model the PSRs are proposed, including a new Livermore Solver (LSODA) for integrating stiff ODEs and a new correction to calculate the PSR time scales. All models have been implemented as a stand-alone application \(\text {edcPisoFoam}\) based on the OpenFOAM technology. Additionally, several RAS calculations were performed using the Turbulence Flame Speed Closure model in Ansys Fluent to assess effects of the heat losses by modeling the conjugate heat transfer between the bluff-body and the reactive flow. Effects of the turbulence Schmidt number on RAS results are discussed as well. Numerical results are compared with available experimental data. Reasonable consistency between experimental data and numerical results provided by RAS, SAS and LES is observed. In general, there is satisfactory agreement between present LES-EDC simulations, numerical results by other authors and measurements without any major modification to the EDC closure constants, which gives a quite reasonable indication on the adequacy and accuracy of the method and its further application for turbulent premixed combustion simulations.     

Does Density Ratio Significantly Affect Turbulent Flame Speed?

In order to experimentally study whether or not the density ratio σ substantially affects flame displacement speed at low and moderate turbulent intensities, two stoichiometric methane/oxygen/nitrogen mixtures characterized by the same laminar flame speed S _L = 0.36 m/s, but substantially different σ were designed using (i) preheating from T _u = 298 to 423 K in order to increase S _L , but to decrease σ, and (ii) dilution with nitrogen in order to further decrease σ and to reduce S _L back to the initial value. As a result, the density ratio was reduced from 7.52 to 4.95. In both reference and preheated/diluted cases, direct images of statistically spherical laminar and turbulent flames that expanded after spark ignition in the center of a large 3D cruciform burner were recorded and processed in order to evaluate the mean flame radius \(\bar {R}_{f}\left (t \right )\) and flame displacement speed \(S_{t}=\sigma ^{-1}{d\bar {R}_{f}} \left / \right . {dt}\) with respect to unburned gas. The use of two counter-rotating fans and perforated plates for near-isotropic turbulence generation allowed us to vary the rms turbulent velocity \(u^{\prime }\) by changing the fan frequency. In this study, \(u^{\prime }\) was varied from 0.14 to 1.39 m/s. For each set of initial conditions (two different mixture compositions, two different temperatures T _u , and six different \(u^{\prime })\), five (respectively, three) statistically equivalent runs were performed in turbulent (respectively, laminar) environment. The obtained experimental data do not show any significant effect of the density ratio on S _t . Moreover, the flame displacement speeds measured at u′/S _L = 0.4 are close to the laminar flame speeds in all investigated cases. These results imply, in particular, a minor effect of the density ratio on flame displacement speed in spark ignition engines and support simulations of the engine combustion using models that (i) do not allow for effects of the density ratio on S _t and (ii) have been validated against experimental data obtained under the room conditions, i.e. at higher σ.