To gain insight into the differences between onshore and offshore atmospheric turbulence, 2 pressure fluctuations were measured for offshore wind under different environmental conditions. 3 A durable piezo-electric sensor was used to sample turbulent pressure data at 50 kHz. Offshore 4 measurements were performed at 100 m height on Germany’s FINO3 offshore platform in the 5 German Bight together with additional meteorological data provided by Deutscher Wetterdienst 6 (DWD). The statistical evaluation revealed that the stability state in the atmospheric boundary has a 7 large impact on turbulent fluctuations. Therefore, we used higher statistical properties (described 8 by so-called shape factors) to the stability state. Data was classified to be either within the unstable, 9 neutral or stable stratification. We found that in case of stable stratification, the shape factor is 10 mostly close to zero, indicating that a thermally stable environment produces closer-to Gaussian 11 distributions. Non-Gaussian distributions were found in unstable and neutral boundary layer states 12 and an occurrence probability was estimated. Possible impact on laminar-turbulent transition on the 13 blade is discussed with application of so-called laminar aerofoils on wind turbine blades. Use of a 14 cut-off frequency to separate load and aerodynamic turbulence is proposed.

The exact similarity solutions (also called as special exact solutions) of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, no similarity solutions for two dimensional turbulent boundary layers have been reported in the literature. With the help of dimensional analysis and invariance principle, Prandtl mixing length $\ell=\kappa y$ for one dimensional turbulent boundary layer is extended to $\ell(x,y)=\kappa y (1-\frac{y}{\delta})\sqrt{\frac{\nu}{U\delta}} $ for the two dimensional turbulent boundary layers, furthermore, with a similarity transformation, we successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation $f'''+ ff''+\beta(1-f'^2)+\kappa^2[\eta^2(1-\eta)^2f''|f''|]'=0$. As an application, similarity solutions of the two dimensional turbulent boundary layer on a flat plate at zero incidence have been studied in detail. To solve the ordinary differential equation numerically, a complete Maple code is provided.

The exact similarity solutions of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, for two dimensional turbulent boundary layers, no similarity solutions (special exact solutions) have ever been found. In the light of Blasius' pioneer works, we extend Blasius similarity transformation to the two dimensional turbulent boundary layers, and successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation. By author's Maple code, we numerically solve the ordinary differential equation and produce some useful quantities.

A formula for the turbulent Prandtl number is derived in terms of local variables available from two-equations turbulence models. The derivation is a direct consequence of the expected square additivity of the molecular and turbulent parameters defining the effective viscosity and the effective conductivity. The formula does not degenerate and leads to a Kays like formulation if approximated. Its utility and its limitation are discussed.

We present the results of direct numerical simulations of power spectral densities for kinetic energy, convective entropy and heat flux for unsteady Rayleigh-Bénard magnetoconvection in the frequency space. For larger values of frequency, the power spectral densities for all the global quantities vary with frequency f as f⁻². The scaling exponent is independent of Rayleigh number, Chandrasekhar’s number and thermal Prandtl number.

Polymeric drag reducers have been developed over many years due to the great number of practical applications. In all of them, the molecular stability is an essential requirement. Usually, polymers break down under turbulent flows, which causes a decrease in their efficiency as drag reducers. Besides that, some specific applications, in agro and biomedical fields, impose a specific requirement that must be fulfilled, which is the use of non-toxic materials. A suitable stable material that is elected to accomplish this necessity is the mucilage of aloe vera, which is a bio-polymer that can be used as an alternative to the synthetic ones. Here, we investigate the role played by the aging of aloe vera on its capacity to reduce drag. The results obtained by ¹H nuclear magnetic resonance indicate that the compositions of young and mature leaves of aloe vera are different and such a difference plays an important role on their efficiency as drag reducers. Tests were performed to analyse the drag reduction in a rotating apparatus and in a pipeline system and the efficiencies of leaves of different ages were compared to their composition. The main conclusion of these experiments is that the young mucilage samples, which are richer in complex polysaccharides and exhibit lower acid contents, are more efficient drag reducers.

Numerical simulations of turbulent flows in a stirred dead-end membrane bioreactor are performed by the RNG k-ε model based on finite volume method using Fluent codes. Comparisons of numerical and experimental results confirm the reliability and feasibility of the constructed model. Flow structures such as wake flows and circulation loops in stirred flows were well simulated. An increase of stirring speed is proposed to use to minimize the low velocity region. The single vane stirrer is found to be beneficial for biological separations. Results reveal that the increase of vane number can enhance the mixing effect in flow domains. However, a circular disk stirrer goes against the formation of vertical circulations. The six-vane stirrer is found to be able to provide a uniform distribution of high shear stress.

A century-old scientific conundrum is solved in this paper. The Prandtl mixing length modelled plane boundary turbulent flow is described by: $\frac{du^+}{d{y^+} }+\kappa^2{y^+} ^2\left( \frac{du^+}{d{y^+} }\right)^2=1$, together with boundary condition $ {y^+} =0:\, u^+=0$. Only an approximate solution to this nonlinear ordinary differential equation (ODE) has been sought so far, however, the exact solution to this ODE has not been obtained. By introducing a transformation, $2\kappa y^+=\sinh \xi$, I successfully find the exact solution of the ODE as follows: $u^+=\frac{1}{\kappa}\ln(2\kappa {y^+} +\sqrt{1+4\kappa^2{y^+} ^2})-\frac{2{y^+} }{1+\sqrt{1+4\kappa^2{y^+} ^2}}$.

In Shark Bay, a large hypersaline bay in Western Australia, longitudinal density gradients force gravitational circulation that is important for Bay-ocean exchange. Observations of vertical stratification and velocity are presented, confirming a steady near-bed dense water outflow from Shark Bay’s northern Geographe Channel that persisted through all stages of the tide. Outflow velocities were 2–3 times stronger than the outflows recorded previously in Naturaliste Channel (in the west) and were more resistant to breakdown by tidal mixing. Estimates of turbulent kinetic energy production derived from ADCP data using the variance method showed a more complex structure in Geographe Channel, due to shear associated with the stratified conditions, with peak levels of turbulence occurring during reversal of tidal flows. For both channels the main source of turbulence was tidal flow along the seabed, with the bottom current speed cubed, |U_b³|, providing a reasonable proxy for tidal mixing and prediction of dense water outflows from Shark Bay the majority of the time. Orientation and deeper water of Geographe Channel along the main axis of the longitudinal density gradient provided an explanation for the predominant outflow from the Bay’s northern entrance. These density-driven currents could potentially influence recruitment of commercially fished scallops and prawns through dispersal and flushing of larvae.

The study of turbulent heat transfer in liquid metal flows has gained interest because of applications in several industrial fields. The common assumption of similarity between the dynamical and thermal turbulence, namely the Reynolds analogy, has been proven to be not valid for these fluids. Many methods have been proposed in order to overcome the difficulties encountered in a proper definition of the turbulent heat flux, such as global or local correlations for the turbulent Prandtl number or four parameter turbulence models. In this work we assess a four parameter logarithmic turbulence model for liquid metals based on RANS approach. Several simulation results considering fluids with Pr = 0.01 and Pr = 0.025 are reported in order to show the validity of this approach. The Kays turbulence model is also assessed and compared with integral heat transfer correlations for a wide range of Peclet numbers.

Widely used in hydraulics, the Colebrook equation for flow friction relates implicitly to the input parameters; the Reynolds number, Re and the relative roughness of inner pipe surface, ε/D with the output unknown parameter; the flow friction factor, λ; λ=f(λ, Re, ε/D). In this paper, a few explicit approximations to the Colebrook equation; λ≈f(Re, ε/D), are generated using the ability of artificial intelligence to make inner patterns to connect input and output parameters in explicit way not knowing their nature or the physical law that connects them, but only knowing raw numbers, {Re, ε/D}→{λ}. The fact that the used genetic programming tool does not know the structure of the Colebrook equation which is based on computationally expensive logarithmic law, is used to obtain better structure of the approximations which is less demanding for calculation but also enough accurate. All generated approximations are with low computational cost because they contain a limited number of logarithmic forms used although for normalization of input parameters or for acceleration, but they are also sufficiently accurate. The relative error regarding the friction factor λ, in best case is up to 0.13% with only two logarithmic forms used. As the second logarithm can be accurately approximated by the Padé approximation, practically the same error is obtained also using only one logarithm.

The turbulent Poiseuille flow between two parallel plates is one of the simplest possible physical situations, and it has been studied intensively. In this paper, we propose a modified Prandtl-van Driest mixing length that satisfies both boundary conditions and wall damping effects. With our new formulations, we numerically solve the problem and, moreover, propose an approximate analytical solution of mean velocity. As applications of our solution, an approximate analytical friction coefficient of turbulent Poiseuille flow is proposed.

The problem of fluid dynamics can be greatly simplified if, for every point in space, the strain-rate tensor is diagonalized. This tensor is introduced into the Navier-Stokes equations via material law and divergence of the stress tensor. This article shows that local SO(3)xU(1) gauge fields can be used to locally diagonalize the diffusion components of the strain-rate tensor. The gauge fields resulting from the connection can be interpreted as convection components of the flow, they show properties of quasiparticles and can be interpreted as elementary vortices. Thus, the proposed approach not only offers new insights for the solution and situative simplification of the Navier-Stokes equations, it also uncovers hidden symmetries within the flow convection, allowing - depending on boundary conditions - further interpretation.

The aim of work is to derive an explicit expression for a function of vertical mixing induced by wind-waves. To this end, in the Navier-Stokes equations, a current is decomposed into four constituents: the mean flow, the wave-orbital motion, the wave-induced turbulent and the background turbulent currents. This decomposition allows separating the wave-induced Reynolds stress, Rw, from the background one, Rb. To make a statistical closure for Rw, the Prandtl approach for the background turbulent fluctuations is used that results in an implicit expression for the wave-induced vertical mixing function, Bv. Expression for Bv is specified based on the author’s results for the eddy viscosity found earlier in the frame of the three-layer concept for a wavy air–sea interface, used for modelling wind-drift currents [1]. Finally, the explicit parameterization for Bv(a, u*, z) is found as a linear function in both the wave amplitude at depth z, a(z), and the friction velocity in the air, u*. The linear dependence of function Bv(a) on the wave amplitude provides the enhanced vertical mixing induced by wind–waves in comparison with function Bv(a) having the cubic dependence found in [2], as far as the wind-wave amplitude a(z) decays exponentially with depth.

new subsonic blower wind tunnel design has been studied both numerically and experimentally, it is also referred to as “blower” wind tunnel. This paper, is initially aimed to address each sequential stage of the wind tunnel design process. Rather than applying the standard method of modelling solely the flow in the test section, a large-scale CFD model of the whole wind tunnel was employed. The loss of every constituent element was calculated and then all the losses are added up to determine the power needed for the wind tunnel operation which is used as “intake fan” boundary conditions in the CFD model. Then, flow uniformity and turbulent intensity measurements in an empty test section using a Pitot-static tube and hot wire anemometer (HWA) were introduced to validate the CFD results. The results showed that flow quality was significantly affected by boundary layer controllers (honeycomb and mesh screens) in the settling chamber and wide-angle diffuser. Investigations were also conducted to evaluate the flow deficit in the wake area behind a convex hump model using HWA. This was additional experimental tests carried out to validate the suitability of the wind tunnel flow aerodynamic research.

Vegetation has a significant influence on velocity distribution and turbulent energy in a confluence channel. A laboratory measurement with ADV was used to investigate the flow through a Y-shaped confluence channel partially covered with rigid vegetation on its inner bank. In this study, the flow velocities in cases with and without vegetation were measured by the ADV in a Y-shaped confluence channel. The results clearly show that the existence of non-submerged rigid plants has changed the internal flow structure, that the velocity in the non-vegetated area is greater than in the vegetated area, and that there is a large exchange of mass and momentum between the vegetated and non-vegetated areas. The velocity on both sides is significantly reduced when vegetation is present. In the vicinity of tributaries, due to the presence of vegetation, the high-velocity area moved rapidly to the middle of the non-vegetated area, and the secondary flow phenomenon disappeared. In the mainstream, when vegetation was present, circulation disappeared, and the degree of lateral mixing decreased. The presence of vegetation caused a great change in the internal flow structure and made the flow in non-vegetated areas more intense.

The Colebrook equation ζ is implicitly given in respect to the unknown flow friction factor λ ; λ=ζ(Re,ε*,λ) which cannot be expressed explicitly in exact way without simplifications and use of approximate calculus. Common approach to solve it is through the Newton-Raphson iterative procedure or through the fixed-point iterative procedure. Both requires in some case even eight iterations. On the other hand numerous more powerful iterative methods such as three-or two-point methods, etc. are available. The purpose is to choose optimal iterative method in order to solve the implicit Colebrook equation for flow friction accurately using the least possible number of iterations. The methods are thoroughly tested and those which require the least possible number of iterations to reach the accurate solution are identified. The most powerful three-point methods require in worst case only two iterations to reach final solution. The recommended representatives are Sharma-Guha-Gupta, Sharma-Sharma, Sharma-Arora, Džunić-Petković-Petković; Bi-Ren-Wu, Chun-Neta based on Kung-Traub, Neta, and Jain method based on Steffensen scheme. The recommended iterative methods can reach the final accurate solution with the least possible number of iterations. The approach is hybrid between iterative procedure and one-step explicit approximations and can be used in engineering design for initial rough, but also for final fine calculations.

The three-dimensional single fin configuration finds application in an intake geometry where the cowl-shock wave interacts with the side-wall boundary-layer. Accurate numerical simulation of such three-dimensional shock/turbulent boundary-layer interaction flows, which are characterized by the appearance of strong crossflow separation, is a challenging task. Reynolds-averaged Navier-Stokes computations using the shock-unsteadiness modified Spalart-Allmaras model is carried out at Mach of 5 at large fin angle of 23^◦. The computed results using the modified model are compared to the standard Spalart-Allmaras model and validated against the experimental data. The focus of work is to implement the modified model and to study the flow physics in detail in the complex region of swept-shock-wave turbulent boundary-layer interaction in terms of the shock structure, expansion fan, shear layer and the surface streamlines. The flow structure is correlated to the wall pressure and skin friction in detail. It is observed that the standard model predicts an initial pressure location downstream of the experiments. The modified model reduces the eddy viscosity at the shock and predicts close to the experiments. Overall, the surface pressure using modified model is predicted accurately at all the locations. The skin friction is under predicted by both the models in the reattachment region and is attributed to the poor performance of turbulence models due to flow laminarization.

The Colebrook equation is a popular model for estimating friction loss coefficients in water and gas pipes. The model is implicit in the unknown flow friction factor . To date, the captured flow friction factor can be extracted from the logarithmic form analytically only in the term of the Lambert -function. The purpose of this study is to find an accurate and computationally efficient solution based on the shifted Lambert -function also known as the Wright -function. The Wright -function is more suitable because it overcomes the problem with the overflow error by switching the fast growing term of the Lambert -function to the series expansions that further can be easily evaluated in computers without causing overflow run-time errors. Although the Colebrook equation transformed through the Lambert -function is identical to the original expression in term of accuracy, a further evaluation of the Lambert -function can be only approximate. Very accurate explicit approximations of the Colebrook equation that contains only one or two logarithms are shown. The final result is an accurate explicit approximation of the Colebrook equation with the relative error of no more than 0.0096%. The presented approximations are in the form suitable for everyday engineering use, they are both accurate and computationally efficient.

The eighty years old empirical Colebrook function widely used as an informal standard for hydraulic resistance relates implicitly the unknown flow friction factor , with the known Reynolds number and the known relative roughness of a pipe inner surface ; . It is based on logarithmic law in the form that captures the unknown flow friction factor in a way from which it cannot be extracted analytically. As an alternative to the explicit approximations or to the iterative procedures that require at least a few evaluations of computationally expensive logarithmic function or non-integer powers, this paper offers an accurate and computationally cheap iterative algorithm based on Padé polynomials with only one -call in total for the whole procedure (expensive -calls are substituted with Padé polynomials in each iteration with the exception of the first). The proposed modification is computationally less demanding compared with the standard approaches of engineering practice, but does not influence the accuracy or the number of iterations required to reach the final balanced solution.

This paper gives a new unified formula for the Newtonian fluids valid for all pipe flow regimes from laminar to the fully rough turbulent. It includes laminar, unstable sharp jump from laminar to turbulent, and all types of the turbulent regimes: smooth turbulent regime, partial non-fully developed turbulent and fully developed rough turbulent regime. The formula follows the inflectional form of curves as suggested in Nikuradse’s experiment rather than monotonic shape proposed by Colebrook and White. The composition of the proposed unified formula consists of switching functions and of the interchangeable formulas for laminar, smooth turbulent and fully rough turbulent flow. The proposed switching functions provide a smooth and a computationally cheap transition among hydraulic regimes. Thus, the here presented formulation represents a coherent hydraulic model suitable for engineering use. The model is compared to existing literature models, and shows smooth and computationally cheap transitions among hydraulic regimes.

Today, Colebrook’s equation is mostly accepted as an informal standard for modeling of turbulent flow in hydraulically smooth and rough pipes including transient zone in between. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of inner pipe surface (ε/D). It is implicit in unknown friction factor (λ). Implicit Colebrook’s equation cannot be rearranged to derive friction factor (λ) directly and therefore it can be solved only iteratively [λ=f(λ, R, ε/D)] or using its explicit approximations [λ≈f(R, ε/D)]. Of course, approximations carry in certain error compared with the iterative solution where the highest level of accuracy can be reached after enough number of iterations. The explicit approximations give a relatively good prediction of the friction factor (λ) and can reproduce accurately Colebrook’s equation and its Moody’s plot. Usually, more complex models of approximations are more accurate and vice versa. In this paper, numerical values of parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After this improvement computational burden stays unchanged while accuracy of approximations increases in some of the cases very significantly.

An overview is presented of several diverse branches of work in the area of effectively 2D fluid equilibria which have in common that they are constrained by an infinite number of conservation laws. Broad concepts, and the enormous variety of physical phenomena that can be explored, are highlighted. These span, roughly in order of increasing complexity, Euler flow, nonlinear Rossby waves, 3D axisymmetric flow, shallow water dynamics, and 2D magnetohydrodynamics. The classical field theories describing these systems bear some resemblance to perhaps more familiar fluctuating membrane and continuous spin models, but the fluid physics drives these models into unconventional regimes exhibiting large scale jet and eddy structures. From a dynamical point of view these structures are the end result of various conserved variable forward and inverse cascades. The resulting balance between large scale structure and small scale fluctuations is controlled by the competition between energy and entropy in the system free energy, in turn highly tunable through setting the values of the conserved integrals. Although the statistical mechanical description of such systems is fully self-consistent, with remarkable mathematical structure and diversity of solutions, great care must be taken because the underlying assumptions, especially ergodicity, can be violated or at minimum lead to exceedingly long equilibration times. Generalization of the theory to include weak driving and dissipation (e.g., non-equilibrium statistical mechanics and associated linear response formalism) could provide additional insights, but has yet to be properly explored.

Empirical Colebrook equation from 1939 is still accepted as an informal standard to calculate friction factor during the turbulent flow (4000 < Re < 10⁸) through pipes from smooth with almost negligible relative roughness (ε/D→0) to the very rough (up to ε/D = 0.05) inner surface. The Colebrook equation contains flow friction factor λ in implicit logarithmic form where it is, aside of itself; λ, a function of the Reynolds number Re and the relative roughness of inner pipe surface ε/D; λ = f (λ, Re, ε/D). To evaluate the error introduced by many available explicit approximations to the Colebrook equation, λ ≈ f(Re, ε/D), it is necessary to determinate value of the friction factor λ from the Colebrook equation as accurate as possible. The most accurate way to achieve that is using some kind of iterative methods. Usually classical approach also known as simple fixed point method requires up to 8 iterations to achieve the high level of accuracy, but does not require derivatives of the Colebrook function as here presented accelerated Householder’s approach (3^rd order, 2^nd order: Halley’s and Schröder’s method and 1^st order: Newton-Raphson) which needs only 3 to 7 iteration and three-point iterative methods which needs only 1 to 4 iteration to achieve the same high level of accuracy. Strategies how to find derivatives of the Colebrook function in symbolic form, how to avoid use of the derivatives (Secant method) and how to choose optimal starting point for the iterative procedure are shown. Householder’s approach to the Colebrook’s equations expressed through the Lambert W-function is also analyzed. One approximation to the Colebrook equation based on the analysis from the paper with the error of no more than 0.0617% is shown.

Downwind odor characteristics can be very different depending upon the size of the upwind point-source, interim topography, and wind conditions. At one extreme, the downwind odor plume from a relatively large, confined animal feeding operation (CAFO), located on a flat open plain and under stable, near-straight-line wind conditions can be rather broad, sustained and predictable relative to a fixed receptor site downwind. In contrast, the plume from a small point-source (e.g., a roof vent stack) located on irregular topography and under rapidly shifting wind conditions can be intermittent and fleeting. These transient odor events can be surprisingly intense and offensive, in spite of their fleeting occurrence and perception. This work reports on efforts to optimize an environmental odor sampling strategy, which is optimized for the challenges of (1) sampling of such transient odor 'spikes' and (2) the prioritization of odors/odorants from multiple, closely co-located point-sources, under such transient event conditions. Protocol refinement has emerged by way of 2 environmental odor assessment projects which have been undertaken on behalf of the Missouri Department of Natural Resources. The challenge of transient odor events has been mitigated utilizing rapid, odor cued whole-air grab capture sampling into metalized-FEP gas sampling bags, followed by immediate adsorption transfer onto SPME fibers or sorbent tubes for stabilization during the shipment and storage interval between collection and final analysis. Initial results demonstrated approximately 11 fold increases in target odorant yields for 900 mL sorbent tube transfers from 2-3 second 'burst' odor event bag-captures, as compared to equivalent direct collections at the same downwind receptor location but during perceived (stable) odor 'lull' periods. Results-to-date targeting refinement and field trials of this integrated environmental odor assessment strategy are presented. Preliminary application targeting general odor sampling and point-source differentiation utilizing tracer gases is also presented.

The aim of the studies was to test the hydraulic pressure losses along the length, local pressure losses, absolute and relative roughness values of the pipeline walls made of modern polymer material. The article presents the results of the experimental studies of hydraulic resistance of pipelines from modern polymer materials, the method of theoretical determination of resistance coefficients in the VALTEC ISO 21003 metal polymer pipeline, the classification and method of determining the values of local hydraulic resistance coefficients of T-Junctions are given, the dependence for determining the Darcy λ hydraulic friction coefficient associated with the value of the equivalent hydraulic roughness coefficient of the inner surface of pipes and the Reynolds number. The experimental data of the hydraulic coefficient of friction Darcy λ was compared with the known theoretical relationships of A.D. Altschul, Colbrooke for smooth pipes, Colbrooke and White, an alternative to the Colbrooke-White, Offengenden equation for low-shear pipes and pipes from polymer materials with Reynolds numbers Re>10^4. Dependence of hydraulic coefficient of friction on relative equivalent roughness of pipeline walls and Reynolds number, specific energy loss per one linear meter on total flow rate in system of hydraulically short pipeline is obtained. The empirical dependence of values of coefficients of local resistance ζ=f(Q2/Q3) from a ratio of the divided expenses of Q2/Q3 of the equal-pass T-Junction α=90 the pipeline to a branch up, on pass in the direct T-Junction and division of streams at various internal pipeline diameters of VALTEC.