Control of two-dimensional turbulent wall jet on a coanda surface

The Coanda effect that is the tendency of a fluid jet to stay attached to an adjacent curved surface that is very well shaped has been employed to improve the performance of various devices. The main objective of this paper is to investigate ways of keeping the flow attached to a larger length of a Coanda surface. There are considered two possibilities: one passive, which uses a slot that connects the low pressure and high pressure points on the Coanda surface and an active one, based on the principle of synthetic jet, created through an orifice located near the point of detachment of the jet. Reynolds averaged Navier-Stokes simulations (RANS) with shear stress transport k-ω (SST model) of Menter have been used to compute the two-dimensional turbulent wall jet. The numerical results are presented for the two methods considered. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A dynamic subgrid-scale model for the large eddy simulation of stratified flow

A new dynamic subgrid-scale (SGS) model, including subgrid turbulent stress and heat flux models for stratified shear flow is proposed by using Yoshizawa’s eddy viscosity model as a base model. Based on our calculated results, the dynamic subgrid-scale model developed here is effective for the large eddy simulation (LES) of stratified turbulent channel flows. The new SGS model is then applied to the large eddy simulation of stratified turbulent channel flow under gravity to investigate the coupled shear and buoyancy effects on the near-wall turbulent statistics and the turbulent heat transfer at different Richardson numbers. The critical Richardson number predicted by the present calculation is in good agreement with the value of theoretical analysis     

Asymptotic behaviour of commutation errors and the divergence of the Reynolds stress tensor near the wall in the turbulent channel flow

The derivation of the space averaged Navier–Stokes equations for the large eddy simulation (LES) of turbulent incompressible flows introduces two groups of terms which do not depend only on the space averaged flow field variables: the divergence of the Reynolds stress tensor and commutation errors. Whereas the former is studied intensively in the literature, the latter terms are usually neglected. This note studies the asymptotic behaviour of these terms for the turbulent channel flow at a wall in the case that the commutation errors arise from the application of a non‐uniform box filter. To perform analytical calculations, the unknown flow field is modelled by a wall law (Reichardt law and 1/αth power law) for the mean velocity profile and highly oscillating functions model the turbulent fluctuations. The asymptotics show that near the wall, the commutation errors are at least as important as the divergence of the Reynolds stress tensor. Copyright © 2006 John Wiley & Sons, Ltd.     

Laws of the wall for velocity and temperature in supersonic turbulent boundary layers

Scaling laws for velocity and temperature profiles in the near-wall region of sub- and supersonic turbulent boundary layers have been developed, which allow us to represent velocity and temperature profiles in compressible gas stream in terms of those in an incompressible boundary layer. They are obtained as asymptotic expansions of the solutions to the Reynolds equations in a small parameter — the Mach number based on the friction velocity and gas enthalpy on the wall. The leading term of the expansion for velocity corresponds to known Van Driest's formula. However, the obtained solution contains additional terms of order unity, which explains the contradiction between Van Driest's formula and experimental data. The law of the wall for temperature, which has been formulated for the first time, has an analogous structure. Besides the von Kármán constant and the turbulent Prandtl number in the logarithmic region, known for incompressible flow, the obtained relations contain three new universal constants, which do not depend on gas molecular properties and the specific heat ratio. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A stochastic method to account for the ambient turbulence in Lagrangian Vortex computations

This paper describes a detailed implementation of the Synthetic Eddy Method (SEM) initially presented in Jarrin et al. (2006) applied to the Lagrangian Vortex simulation. While the treatment of turbulent diffusion is already extensively covered in scientific literature, this is one of the first attempts to represent ambient turbulence in a fully Lagrangian framework. This implementation is well suited to the integration of PSE (Particle Strength Exchange) or DVM (Diffusion Velocity Method), often used to account for molecular and turbulent diffusion in Lagrangian simulations. The adaptation and implementation of the SEM into a Lagrangian method using the PSE diffusion model is presented, and the turbulent velocity fields produced by this method are then analysed. In this adaptation, SEM turbulent structures are simply advected, without stretching or diffusion of their own, over the flow domain. This implementation proves its ability to produce turbulent velocity fields in accordance with any desired turbulent flow parameters. As the SEM is a purely mathematical and stochastic model, turbulent spectra and turbulent length scales are also investigated. With the addition of variation in the turbulent structures sizes, a satisfying representation of turbulent spectra is recovered, and a linear relation is obtained between the turbulent structures sizes and the Taylor macroscale. Lastly, the model is applied to the computation of a tidal turbine wake for different ambient turbulence levels, demonstrating the ability of this new implementation to emulate experimentally observed tendencies.     

An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.     

Turbulent Poiseuille flow with near-critical wall transpiration

An incompressible, pressure-driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress with the first and the second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary-value problem for a second-order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. The case of near-critical transpiration, when the shear stress at the injection wall vanishes, is considered. It is shown that the maximum point on the mean velocity profile lies in a thin sublayer near the suction wall in this case. A formula for the position of the maximum point as a function of the transpiration factor is obtained. The mean velocity profiles near the suction wall are calculated. A friction law for Poiseuille flow with near-critical transpiration is found, which makes it possible to describe the relation between the shear stress at the wall, the Reynolds number, and the transpiration velocity by a single function of one variable. Direct numerical simulation of the flow for some transpiration factors is performed. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic study of turbulent Poiseuille flow with wall transpiration

An incompressible, pressure–driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress to the first and second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary–value problem for a second–order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. A limiting transpiration velocity is obtained, such that the shear stress at the injection wall vanishes, while the maximum point on the velocity profile approaches the suction wall. In this case, a sublayer near the suction wall appears where the mean velocity is proportional to the square root of the distance from the wall. A friction law for Poiseuille flow with transpiration is found, which makes it possible to describe the relation between the wall shear stress, the Reynolds number, and the transpiration velocity by a function of one variable. A velocity defect law, which generalizes the classical law for the core region in a channel with impermeable walls to the case of transpiration, is also established. In similarity variables, the mean velocity profiles across the whole channel width outside viscous sublayers can be described by a one–parameter family of curves. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Calculation of Acoustic Sources for the Landing Gear of Aeroplane during Take-off and Landing

In this paper we present Orthogonal Sub-grid scale stabilized method with static subscales applied on turbulent flow around landing gear of aircraft during take-off and landing. Numerical results of proposed method are compared with famous LES Smagorinsky method for turbulent flow and it is shown how obtained results affect allocation of acoustic sources and then propagation of acoustic waves in computational domain. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical modeling of turbulence-induced interfacial instability in two-phase flow with moving interface

The present paper introduces a new interfacial marker-level set method (IMLS) which is coupled with the Reynolds averaged Navier–Stokes (RANS) equations to predict the turbulence-induced interfacial instability of two-phase flow with moving interface. The governing RANS equations for time-dependent, axisymmetric and incompressible two-phase flow are described in both phases and solved separately using the control volume approach on structured cell-centered collocated grids. The transition from one phase to another is performed through a consistent balance of kinematic and dynamic conditions on the interface separating the two phases. The topological changes of the interface are predicted by applying the level set approach. By fitting a number of interfacial markers on the intersection points of the computational grids with the interface, the interfacial stresses and consequently, the interfacial driving forces are easily estimated. Moreover, the normal interface velocity, calculated at the interfacial markers positions, can be extended to the higher dimensional level set function and used for the interface advection process. The performance of linear and non-linear two-equation k–ε turbulence models is investigated in the context of the considered two-phase flow impinging problem, where a turbulent gas jet impinging on a free liquid surface. The numerical results obtained are evaluated through the comparison with the available experimental and analytical data. The nonlinear turbulence model showed superiority in predicting the interface deformation resulting from turbulent normal stresses. However, both linear and nonlinear turbulence models showed a similar behavior in predicting the interface deformation due to turbulent tangential stresses. In general, the developed IMLS numerical method showed a remarkable capability in predicting the dynamics of the considered two-phase immiscible flow problems and therefore it can be applied to quite a number of interface stability problems.     

Langevin equation of a fluid particle in wall-induced turbulence

We derive the Langevin equation describing the stochastic process of fluid particle motion in wall-induced turbulence (turbulent flow in pipes, channels, and boundary layers including the atmospheric surface layer). The analysis is based on the asymptotic behavior at a large Reynolds number. We use the Lagrangian Kolmogorov theory, recently derived asymptotic expressions for the spatial distribution of turbulent energy dissipation, and also newly derived reciprocity relations analogous to the Onsager relations supplemented with recent measurement results. The long-time limit of the derived Langevin equation yields the diffusion equation for admixture dispersion in wall-induced turbulence.     

Incipient separation over wall irregularities in transonic flow

Incipient separation over wall irregularities in a steady two dimensional flow field of a perfect fluid which has transonic speed characteristics has been investigated considering viscous-inviscid interactions at high Reynolds number. The aim of this work is to investigate dependence of the critical hump height (when a well attached flow over rigid body surface turns into a separated one) on the Karman–Guderley parameter which characterizes of the local flow field. The analysis of the flow field starts with the so-called inspection analysis of the flow properties and then the interaction problem has been constructed using the asymptotic analysis of triple-deck structure of interaction region. Finally, a method based on a semi-direct solution of governing equations of the transonic interaction problem has been used to obtain the numerical solution of the problem.     

Numerical simulation of pollution and oil spill spreading by the stochastic discrete particle method

The features of the stochastic discrete particle method are discussed as applied to the simulation of pollutant advection and diffusion in a turbulent flow and to the spread of a thin film of a viscous substance (oil) on the surface of water. The diffusion tensor in the former problem depends on the scale of the pollution cloud, and the diffusivity in the latter problem depends nonlinearly on the desired function. For pollution dispersion by a turbulent flow, a stochastic discrete particle algorithm is constructed in the case when the diffusion tensor corresponds to the Richardson 4/3 law. The numerical and analytical results are shown to agree well. The problem of oil film spreading is described by a quasilinear advection-diffusion equation. For this problem, a random walking algorithm is constructed in which the variance of the walking particle step depends on the desired function. For both instantaneous and time-continuous sources of pollutants, the solution produced by the stochastic discrete particle method agrees well with the analytical and/or numerical solutions to the test problems under consideration.     

Robin-type wall functions and their numerical implementation

The paper is devoted to numerical implementation of the wall functions of Robin-type for modeling near-wall turbulent flows. The wall functions are based on the transfer of a boundary condition from a wall to some intermediate boundary near the wall. The boundary conditions on the intermediate boundary are of Robin-type and represented in a differential form. The wall functions are formulated in an analytical easy-to-implement form, can take into account the source terms, and do not include free parameters. The relation between the wall functions of Robin type and the theory of Calderon–Ryaben'kii's potentials is demonstrated. A universal robust approach to the implementation of the Robin-type wall functions in finite-volume codes is provided. The example of an impinging jet is considered.     

Exergetic aspects of high-speed boundary layers

Exergy is an established concept for evaluating heat cycles. It quantifies the available work that can be extracted from a system. The extension of the classical definition to non-uniform flow field properties is applied here to high-speed boundary layers. The destruction of exergy is quantified by the corresponding loss thickness and is evaluated locally. Turbulent fluctuations in this approach are described by their turbulent exergy, which is a generalisation of the concept of turbulent kinetic energy for incompressible flow. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Physical Modelling of Point Source Dispersion Inside Atmospheric Boundary Layer

One of the most important problems in the turbulence theory is the process of turbulent mixing. Although many theoretical and experimental studies on the structure functions of a passive scalar as well as the probability density functions of its increment and gradient in a turbulent fluid have been published so far, the instantaneous feature of an advected scalar has not been explored very well. There are only a few experimental studies on a scalar concentration fluctuation analysis in turbulent flows. In this study the time‐series of both longitudinal velocity component and concentration analysis are to be presented. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stationary single waves in turbulent open-channel flow

Steady two-dimensional turbulent open-channel flow is considered. Stationary single-wave solutions are investigated. The fully-developed oncoming flow is slightly supercritical. The Reynolds number is very large. The analysis is kept free of turbulence modelling. As stationary solitary waves cannot exist in turbulent flow for a plane bottom with constant roughness [1], two particular perturbations of the conditions at the channel bottom are examined: 1) We revisit the case [1] where the friction coefficient locally differs slightly by a constant from the reference value upstream; 2) An unevenness of very small height in the channel bottom (bump, ramp) is admitted, with the bottom roughness taken constant. An analogy between these cases is presented. In both cases, three stationary solutions for the surface elevation are found: A stable and an unstable solitary wave, respectively, and a single wave of a second kind with smaller amplitude. For the latter, an analysis for weak dissipation yields a uniformly valid solution that is in good agreement with the numerical results for various parameters. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hybrid LES-RANS Modeling of Complex Turbulent Flows

The paper describes a state-of-the-art hybrid LES-URANS method for the simulation of complex internal and external turbulent flows. Relying on a unified LES-URANS approach with a soft interface the methodology is designed for wall-bounded non-equilibrium flows. The unsteady Reynolds-averaged Navier-Stokes (URANS) mode within the hybrid approach is taken into account by an explicit algebraic Reynolds stress model (EARSM), which guarantees an appropriate representation of the anisotropic near-wall turbulence. All non-closed terms in the transport equation of the turbulent kinetic energy are modeled by enhanced formulations using the EARSM (production and diffusion term) and the splitting of the dissipation rate into a homogeneous and an inhomogeneous contribution. The former is expressed analytically by a Taylor series expansion of the homogeneous lateral Taylor microscale in the vicinity of the wall guaranteeing the correct asymptotic behavior. The latter is incorporated into the diffusion term. The interface location between the large-eddy simulation (LES) mode and the URANS mode is determined automatically on-the-fly based on the modeled turbulent kinetic energy and the distance to the wall. For transitional (external) flows an additional dynamic transition criterion is applied which determines the laminar and the turbulent flow regimes and is combined with the existing interface criterion. An internal flow over a periodic arrangement of hills and an external flow past a SD7003 airfoil with a laminar separation bubble are taken into account for a detailed evaluation of the method. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Integration of infinite chain of transport equations for cumulants

The method of compatible differential constraints is applied to integrate an infinite chain of transport equations for cumulants. The dynamics of a stratified momentumless turbulent planar wake is considered as an example. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)