RANS Modelling of Wake Induced Transition With the Dynamic Intermittency Concept

Author(s):  
Koen Lodefier ◽  
Erik Dick

A transition model for describing wake-induced transition is presented based on the SST turbulence model by Menter and two dynamic equations for intermittency: one for near-wall intermittency and one for free-stream intermittency. In the Navier-Stokes equations, the total intermittency factor, which is the sum of the two, multiplies the turbulent viscosity computed by the turbulence model. The quality of the transition model is illustrated on the T106A test cascade for different levels of inlet free-stream turbulence intensity. The unsteady results are presented in space-time diagrams of shape factor, wall shear stress, momentum thickness and intermittency on the suction side. Results show the capability of the model to capture the physics of unsteady transition. Inevitable shortcomings are also revealed.

Author(s):  
Koen Lodefier ◽  
Erik Dick

A transition model for describing wake-induced transition is presented. It is based on the SST turbulence model by Menter, with the k–ω part in low-Reynolds form according to Wilcox, and two dynamic equations for intermittency: one for near-wall-intermittency and one for free-stream-intermittency. The total intermittency factor, which is the sum of the two, multiplies the turbulent viscosity computed by the turbulence model. The quality of the transition model is illustrated on the T106a test cascade using experimental results for flow with low free-stream turbulence intensity and transition in separated state and for flow with high free-stream turbulence intensity and transition in attached state. The unsteady results are presented in S–T diagrams of the shape factor and wall shear stress on the suction side. Results show the capability of the model to capture the basics of unsteady transition.


2019 ◽  
Vol 22 (3) ◽  
pp. 25-34
Author(s):  
V. A. Vershkov ◽  
B. S. Kritsky ◽  
R. M. Mirgazov

The article considers the problem of the flow around the helicopter main rotor taking into account blades flapping in the plane of rotation and in the plane of thrust as well as the elastic blades deformation. The rotor rotation is modeled by the method of converting Navier-Stokes equations from a fixed coordinate system associated with the incoming flow into a rotating system associated with the rotor hub. For axial flow problems, this makes it possible to formulate the problem as stationary at a constant rotational speed of rotor. For a mode of skewed flow around the rotor in the terms of incident flow in this system it is necessary to solve the non-stationary problem. To solve the problem, the method of deformable grids is used, in which the equations are copied taking into account the grid nodes motion determined in accordance with the spatial blades motion, and SST turbulence model is used for closure. The results of the test calculations of the main rotor aerodynamic characteristics with and without blade flapping are presented in this paper. The coefficients of the main rotor thrust cT and the blades hinge moments mh are compared. The calculations were carried out in the CFD software ANSYS CFX (TsAGI License No. 501024). The flow around a four-bladed main rotor of a radius of 2.5 meters is modeled in the regime of skewed flow. The speed of the incoming flow came to 85 m/s under normal atmospheric conditions. The rotor was at an angle of attack of −10˚. To calculate the rotor motion without taking into account the flapping movements, we used the nonstationary system of Navier-Stokes equations with the closure with SST turbulence model. The calculation was being carried out until the change in the maximum value of the rotor thrust during one revolution became less than 1%. For modeling flapping blade movements, the control laws and equations describing the angle of blade flapping as a function from its azimuth angle obtained from the experiment were used. The procedure for reconstructing the grid according to a given law was conducted using standard grid deformation methods presented in the ANSYS CFX software. When solving the nonstationary Navier-Stokes equations, a dual time step was used. The obtained results show that accounting of the effect of flapping movements and cyclic control of the blades has an impact on the character of changing the main rotor thrust coefficient during one revolution and significantly changes the shape of the graph of the hinge moment coefficient of each blade.


1975 ◽  
Vol 69 (4) ◽  
pp. 631-656 ◽  
Author(s):  
W. Roger Briley ◽  
Henry Mcdonald

A method is presented for performing detailed computations of thin incompressible separation bubbles on smooth surfaces. The analysis consists of finite-difference solutions to the time-dependent boundary-layer or Navier-Stokes equations for the flow in the immediate vicinity of the bubble. The method employs the McDonald-Fish turbulence model, to predict the development of the time-mean flow field, as influenced by the free-stream turbulence level. It also employs a viscous-inviscid interaction model, which accounts for the elliptic interaction between the shear layer and inviscid free stream. The numerical method is based on an alternating-direction implicit scheme for the vorticity equation. It employs transformations, to allow the free-stream boundary to change in time with the shape of the computed shear layer, and to ensure an adequate resolution of the sublayer region. Numerical solutions are presented for transitional bubbles on an NACA 663-018 airfoil at zero angle of incidence with chordal Reynolds numbers of 2·0 × 106 and 1·7 × 106. These have a qualitative behaviour similar to that observed in numerous experiments; they are also in reasonable quantitative agreement with available experimental data. Little difference is found between steady solutions of the boundary-layer and Navier-Stokes equations for these flow conditions. Numerical studies based on mesh refinement suggest that the well-known singularity at separation, which is present in conventional solutions of the steady boundary-layer equations when the free-stream velocity is specified, is effectively removed when viscous-inviscid interaction is allowed to influence the imposed velocity distribution.


1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.


2014 ◽  
Vol 752 ◽  
pp. 602-625 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall

AbstractOur concern in this paper is with high-Reynolds-number nonlinear equilibrium solutions of the Navier–Stokes equations for boundary-layer flows. Here we consider the asymptotic suction boundary layer (ASBL) which we take as a prototype parallel boundary layer. Solutions of the equations of motion are obtained using a homotopy continuation from two known types of solutions for plane Couette flow. At high Reynolds numbers, it is shown that the first type of solution takes the form of a vortex–wave interaction (VWI) state, see Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666), and is located in the main part of the boundary layer. On the other hand, here the second type is found to support an equilibrium solution of the unit-Reynolds-number Navier–Stokes equations in a layer located a distance of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}O(\ln \mathit{Re})$ from the wall. Here $\mathit{Re}$ is the Reynolds number based on the free-stream speed and the unperturbed boundary-layer thickness. The streaky field produced by the interaction grows exponentially below the layer and takes its maximum size within the unperturbed boundary layer. The results suggest the possibility of two distinct types of streaky coherent structures existing, possibly simultaneously, in disturbed boundary layers.


Author(s):  
П.А. Поливанов

In this paper a numerical and experimental study of the effect of blowing/suction through a perforated surface on a turbulent boundary layer at a Mach number M = 1.4 is carried out. Most of the calculations were performed by Reynolds-averaged Navier-Stokes equations with the k-w SST turbulence model. The calculated geometry completely repeated the experimental one including the perforated surface. The numerical data were compared with experimental measurements obtained by the PIV method. Analysis of the data made it possible to find the limits of applicability of the numerical method for this flow.


2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.


1989 ◽  
Vol 111 (3) ◽  
pp. 333-340 ◽  
Author(s):  
J. F. Louis ◽  
A. Salhi

The turbulent flow between two rotating co-axial disks is driven by frictional forces. The prediction of the velocity field can be expected to be very sensitive to the turbulence model used to describe the viscosity close to the walls. Numerical solutions of the Navier–Stokes equations, using a k–ε turbulence model derived from Lam and Bremhorst, are presented and compared with experimental results obtained in two different configurations: a rotating cavity and the outflow between a rotating and stationary disk. The comparison shows good overall agreement with the experimental data and substantial improvements over the results of other analyses using the k–ε models. Based on this validation, the model is applied to the flow between counterrotating disks and it gives the dependence of the radial variation of the tangential wall shear stress on Rossby number.


Author(s):  
J. Steelant ◽  
E. Dick

Turbomachinery flows are characterized by a very high intensity turbulent mean part. As a consequence, laminar flow in boundary layer regions undergoes transition through direct excitation of turbulence. This is the so-called bypass transition. Regions form that are intermittently laminar and turbulent. In particular in accelerating flows, as on the suction side of a turbine blade, this intermittent flow can extend over a very large part of the boundary layer. Classical turbulence modelling based on global time averaging is not valid in intermittent flows. To take correctly account of the intermittency, conditioned averages are necessary. These are averages taken during the fraction of time the flow is turbulent or laminar respectively. Starting from the Navier-Stokes equations, conditioned continuity, momentum and energy equations are derived for the laminar and turbulent parts of an intermittent flow. The turbulence is described by the classical k-ε model. The supplementary parameter introduced by the conditioned averaging is the intermittency factor. In the calculations, this factor is prescribed in an algebraic way.


Author(s):  
Fahua Gu ◽  
Mark R. Anderson

The design of turbomachinery has been focusing on the improvement of the machine efficiency and the reduction of the design cost. This paper presents an integrated design system to create the machine geometry and to predict the machine performance at different levels of approximation, including one-dimensional design and analysis, quasi-three-dimensional-(blade-to-blade, throughflow) and full-three-dimensional-steady-state CFD analysis. One of the most important components, the Reynolds-averaged Navier-Stokes solver, is described in detail. It originated from the Dawes solver with numerous enhancements. They include the use of the low speed pre-conditioned full Navier-Stokes equations, the addition of the Spalart-Allmaras turbulence model and an improvement of wall functions related with the turbulence model. The latest upwind scheme, AUSM, has been implemented too. The Dawes code has been rewritten into a multi-block solver for O, C, and H grids. This paper provides some examples to evaluate the effect of grid topology on the machine performance prediction.


Sign in / Sign up

Export Citation Format

Share Document