scholarly journals Prandtl’s extended mixing length model applied to the two-dimensional turbulent classical far wake

Author(s):  
Ashleigh J. Hutchinson ◽  
Nicholas Hale ◽  
Kendall Born ◽  
David P. Mason

Despite its limitations, Prandtl’s mixing length model is widely applied in modelling turbulent free shear flows. Prandtl’s extended model addresses many of the shortfalls of the original model, but is not so widely used, in part due to additional mathematical complexities that arise in its derivation and implementation. Furthermore, in both models, Prandtl neglects the kinematic viscosity on the basis that it is much smaller in magnitude than the turbulent viscosity. Recent work has shown that including the kinematic viscosity in the original model has both mathematical and physical advantages. In the present work, a novel derivation of the extended model is provided, and it is demonstrated that similar advantages are again obtained when the kinematic viscosity is included. Additionally, through the use of scaling techniques, similarity mean velocity profiles of the extended model are derived, resulting in a single nonlinear ordinary differential equation that is solved numerically with a Hermite spectral method. The computed profiles for the normalized similarity mean velocity and shear stress are compared with experimental observations and shown to be in excellent agreement.

1980 ◽  
Vol 101 (1) ◽  
pp. 159-178 ◽  
Author(s):  
S. Yavuzkurt ◽  
R. J. Moffat ◽  
W. M. Kays

Hydrodynamic data are reported in the companion paper (Yavuzkurt, Moffat & Kays 1980) for a full-coverage film-cooling situation, both for the blown and the recovery regions. Values of the mean velocity, the turbulent shear stress, and the turbulence kinetic energy were measured at various locations, both within the blown region and in the recovery region. The present paper is concerned with an analysis of the recovery region only. Examination of the data suggested that the recovery-region hydrodynamics could be modelled by considering that a new boundary layer began to grow immediately after the cessation of blowing. Distributions of the Prandtl mixing length were calculated from the data using the measured values of mean velocity and turbulent shear stresses. The mixing-length distributions were consistent with the notion of a dual boundary-layer structure in the recovery region. The measured distributions of mixing length were described by using a piecewise continuous but heuristic fit, consistent with the concept of two quasi-independent layers suggested by the general appearance of the data. This distribution of mixing length, together with a set of otherwise normal constants for a two-dimensional boundary layer, successfully predicted all of the observed features of the flow. The program used in these predictions contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length. The program is a two-dimensional, finite-difference program capable of predicting the mean velocity and turbulence kinetic energy profiles based upon initial values, boundary conditions, and a closure condition.


1998 ◽  
Vol 25 (2) ◽  
pp. 232-240 ◽  
Author(s):  
Jean-Loup Robert ◽  
Mohamed Khelifi ◽  
Ahmed Ghanmi

Since the viscous analogy of turbulence was introduced by Reynolds, many formulations for turbulent viscosity have been proposed. One of them, based on the mixing length concept, is investigated here in a broader point of view. The mixing length concept was used to correctly model turbulent velocity profiles for irregular two-dimensional and three-dimensional domains. Two cases of study were investigated for this purpose: a simple two-dimensional aerodynamic problem and a more complicated three-dimensional hydraulic problem. Results showed that the use of a constant viscosity fails to correctly reproduce experimental observations. On the other hand, the use of the mixing length concept leads to a good agreement between the measured and predicted values.Key words: fluid flow, finite element method, mixing length flow theory, turbulent flow, velocity profiles.


2012 ◽  
Vol 134 (2) ◽  
Author(s):  
Rolf-Erik Keck ◽  
Dick Veldkamp ◽  
Helge Aagaard Madsen ◽  
Gunner Larsen

The work presented in this paper focuses on improving the description of wake evolution due to turbulent mixing in the dynamic wake meandering (DWM) model. From wake investigations performed with high-fidelity actuator line simulations carried out in ELLIPSYS3D, it is seen that the current DWM description, where the eddy viscosity is assumed to be constant in each cross-section of the wake, is insufficient. Instead, a two-dimensional eddy viscosity formulation is proposed to model the shear layer generated turbulence in the wake, based on the classical mixing length model. The performance of the modified DWM model is verified by comparing the mean wake velocity distribution with a set of ELLIPSYS3D actuator line calculations. The standard error (defined as the standard deviation of the difference between the mean velocity field of the DWM and the actuator line model), in the wake region extending from 3 to 12 diameters behind the rotor, is reduced by 27% by using the new eddy viscosity formulation.


Author(s):  
S. Yavuzkurt

A one-equation model (mixing-length and turbulent kinetic energy) of turbulence is used for the calculation of the full-coverage and recovery region hydrodynamics over a full-coverage film-cooled surface. The model requires a detailed description of the form and dynamics of the complex mixing-length profile encountered in this type of flow structure. This is achieved through extensive use of experimental input and physical interpretation of the data by combining equations for simple flow structures such as two-dimensional turbulent flat plate boundary layers and jets-in-cross flow. The one-equation model is used in a two-dimensional finite difference boundary layer code giving successful predictions of the spanwise averaged mean velocity and turbulent kinetic energy profiles between injection rows in the full coverage region and also in the recovery region for two blowing ratios (Ujet/U∞ = 0.4, 0.9).


1968 ◽  
Vol 90 (2) ◽  
pp. 301-312 ◽  
Author(s):  
B. S. Massey ◽  
B. R. Clayton

Equations have been developed [1] which describe the flow in a steady, two-dimensional, incompressible, laminar boundary layer on a curved surface. The method of “similar solutions” yields a fourth-order, nonlinear, ordinary differential equation which may be solved on a digital computer. Account is taken of the effects of both surface curvature and displacement thickness, and in this paper attention is given to the influence of these effects on other boundary-layer properties up to and including the position of flow separation.


2003 ◽  
Vol 47 (01) ◽  
pp. 13-23
Author(s):  
R. Muscari ◽  
A. Di Mascio

A numerical model for the simulation of two-dimensional spilling breaking waves is described. The model is derived from Cointe and Tulin's theory of steady breakers (Cointe & Tulin 1994), although some important changes have been introduced in order to obtain a stable algorithm when coupled with steady-state Reynolds averaged Navier-Stokes equations (RANSE) solvers. In particular, the shape of the breaker and its relation with the following wave height differ from the original model, and moreover, additional conditions for the tangential stress and the turbulent viscosity are proposed. The model has been implemented in a RANSE code, developed for the study of ship flows, through a modification in the free-surface boundary conditions below the breaker. This yields a simple but effective way to reproduce the breaker influence on the underlying flow. The algorithm was used for the simulation of the flow past a submerged hydrofoil. The numerical results are compared with the experimental data by Duncan (1983).


2020 ◽  
Vol 640 ◽  
pp. A53
Author(s):  
L. Löhnert ◽  
S. Krätschmer ◽  
A. G. Peeters

Here, we address the turbulent dynamics of the gravitational instability in accretion disks, retaining both radiative cooling and irradiation. Due to radiative cooling, the disk is unstable for all values of the Toomre parameter, and an accurate estimate of the maximum growth rate is derived analytically. A detailed study of the turbulent spectra shows a rapid decay with an azimuthal wave number stronger than ky−3, whereas the spectrum is more broad in the radial direction and shows a scaling in the range kx−3 to kx−2. The radial component of the radial velocity profile consists of a superposition of shocks of different heights, and is similar to that found in Burgers’ turbulence. Assuming saturation occurs through nonlinear wave steepening leading to shock formation, we developed a mixing-length model in which the typical length scale is related to the average radial distance between shocks. Furthermore, since the numerical simulations show that linear drive is necessary in order to sustain turbulence, we used the growth rate of the most unstable mode to estimate the typical timescale. The mixing-length model that was obtained agrees well with numerical simulations. The model gives an analytic expression for the turbulent viscosity as a function of the Toomre parameter and cooling time. It predicts that relevant values of α = 10−3 can be obtained in disks that have a Toomre parameter as high as Q ≈ 10.


1979 ◽  
Vol 44 (3) ◽  
pp. 700-710 ◽  
Author(s):  
Ivan Fořt ◽  
Hans-Otto Möckel ◽  
Jan Drbohlav ◽  
Miroslav Hrach

Profiles of the mean velocity have been analyzed in the stream streaking from the region of rotating standard six-blade disc turbine impeller. The profiles were obtained experimentally using a hot film thermoanemometer probe. The results of the analysis is the determination of the effect of relative size of the impeller and vessel and the kinematic viscosity of the charge on three parameters of the axial profile of the mean velocity in the examined stream. No significant change of the parameter of width of the examined stream and the momentum flux in the stream has been found in the range of parameters d/D ##m <0.25; 0.50> and the Reynolds number for mixing ReM ##m <2.90 . 101; 1 . 105>. However, a significant influence has been found of ReM (at negligible effect of d/D) on the size of the hypothetical source of motion - the radius of the tangential cylindrical jet - a. The proposed phenomenological model of the turbulent stream in region of turbine impeller has been found adequate for values of ReM exceeding 1.0 . 103.


Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1439
Author(s):  
Chaudry Masood Khalique ◽  
Karabo Plaatjie

In this article, we investigate a two-dimensional generalized shallow water wave equation. Lie symmetries of the equation are computed first and then used to perform symmetry reductions. By utilizing the three translation symmetries of the equation, a fourth-order ordinary differential equation is obtained and solved in terms of an incomplete elliptic integral. Moreover, with the aid of Kudryashov’s approach, more closed-form solutions are constructed. In addition, energy and linear momentum conservation laws for the underlying equation are computed by engaging the multiplier approach as well as Noether’s theorem.


1969 ◽  
Vol 35 (2) ◽  
pp. 225-241 ◽  
Author(s):  
M. A. Badri Narayanan ◽  
V. Ramjee

Experiments on reverse transition were conducted in two-dimensional accelerated incompressible turbulent boundary layers. Mean velocity profiles, longitudinal velocity fluctuations $\tilde{u}^{\prime}(=(\overline{u^{\prime 2}})^{\frac{1}{2}})$ and the wall-shearing stress (TW) were measured. The mean velocity profiles show that the wall region adjusts itself to laminar conditions earlier than the outer region. During the reverse transition process, increases in the shape parameter (H) are accompanied by a decrease in the skin friction coefficient (Cf). Profiles of turbulent intensity (u’2) exhibit near similarity in the turbulence decay region. The breakdown of the law of the wall is characterized by the parameter \[ \Delta_p (=\nu[dP/dx]/\rho U^{*3}) = - 0.02, \] where U* is the friction velocity. Downstream of this region the decay of $\tilde{u}^{\prime}$ fluctuations occurred when the momentum thickness Reynolds number (R) decreased roughly below 400.


Sign in / Sign up

Export Citation Format

Share Document