High reynolds numbers turbulent boundary layers

Esta tesis estudia las similitudes y diferencias entre los flujos turbulentos de pared de tipo externo e interno, en régimen incompresible, y a números de Reynolds moderada¬mente altos. Para ello consideramos tanto simulaciones numéricas como experimentos de capas límites con gradiente de presiones nulo y de flujos de canal, ambos a números de Reynolds en el rango δ+ ~ 500 - 2000. Estos flujos de cortadura son objeto de numerosas investigaciones debido a la gran importancia que tienen tanto a nivel tecnológico como a nivel de física fundamental. No obstante, todavía existen muchos interrogantes sobre aspectos básicos tales como la universalidad de los perfiles medios y de fluctuación de las velocidades o de la presión, tanto en la zona cercana a la pared como en la zona logarítmica, el escalado y el efecto del número de Reynolds, o las diferencias entre los flujos internos y externos en la zona exterior. En éste estudio hemos utilizado simulaciones numéricas ya existentes de canales y capas límites a números de Reynolds δ+ ~ 2000 y δ+ ~ 700, respectivamente. Para poder comparar ambos flujos a igual número de Reynolds hemos realizado una nueva simulación directa de capa límite en el rango δ+ ~ 1000-2000. Los resultados de la misma son presentados y analizados en detalle. Los datos sin postprocesar y las estadísticas ya postprocesadas están públicamente disponibles en nuestro sitio web.162 El análisis de las estadísticas usando un único punto confirma la existencia de perfiles logarítmicos para las fluctuaciones de la velocidad transversal w'2+ y de la presión p'2+ en ambos tipos de flujos, pero no para la velocidad normal v'2+ o la velocidad longitudinal u'2+. Para aceptar o rechazar la existencia de un rango logarítmico en u'2+ se requieren números de Reynolds más altos que los considerados en éste trabajo. Una de las conse¬cuencias más importantes de poseer tales perfiles es que el valor máximo de la intensidad, que se alcanza cerca de la pared, depende explícitamente del número de Reynolds. Esto ha sido confirmado tras analizar un gran número de datos experimentales y numéricos, cor¬roborando que el máximo de u'2+, p/2+, y w'2+ aumenta proporcionalmente con el log(δ+). Por otro lado, éste máximo es más intenso en los flujos externos que en los internos. La máxima diferencia ocurre en torno a y/δ ~ 0.3-0.5, siendo esta altura prácticamente independiente del número de Reynolds considerado. Estas diferencias se originan como consecuencia del carácter intermitente de las capas límites, que es inexistente en los flujos internos. La estructura de las fluctuaciones de velocidad y de presión, junto con la de los esfuer¬zos de Reynolds, se han investigado por medio de correlaciones espaciales tridimensionales considerando dos puntos de medida. Hemos obtenido que el tamaño de las mismas es gen¬eralmente mayor en canales que en capas límites, especialmente en el caso de la correlación longitudinal Cuu en la dirección del flujo. Para esta correlación se demuestra que las es¬tructuras débilmente correladas presentan longitudes de hasta 0(75), en el caso de capas límites, y de hasta 0(185) en el caso de canales. Estas longitudes se obtienen respecti-vamente en la zona logarítmica y en la zona exterior. Las longitudes correspondientes en la dirección transversal son significativamente menores en ambos flujos, 0(5 — 25). La organización espacial de las correlaciones es compatible con la de una pareja de rollos casi paralelos con dimensiones que escalan en unidades exteriores. Esta organización se mantiene al menos hasta y ~ 0.65, altura a la cual las capas límites comienzan a organi¬zarse en rollos transversales. Este comportamiento es sin embargo más débil en canales, pudiéndose observar parcialmente a partir de y ~ 0.85. Para estudiar si estas estructuras están onduladas a lo largo de la dirección transver¬sal, hemos calculado las correlaciones condicionadas a eventos intensos de la velocidad transversal w'. Estas correlaciones revelan que la ondulación de la velocidad longitudinal aumenta conforme nos alejamos de la pared, sugiriendo que las estructuras están más alineadas en la zona cercana a la pared que en la zona lejana a ella. El por qué de esta ondulación se encuentra posiblemente en la configuración a lo largo de diagonales que presenta w'. Estas estructuras no sólo están onduladas, sino que también están inclinadas respecto a la pared con ángulos que dependen de la variable considerada, de la altura, y de el contorno de correlación seleccionado. Por encima de la zona tampón e independien¬temente del número de Reynolds y tipo de flujo, Cuu presenta una inclinación máxima de unos 10°, las correlaciones Cvv y Cm son esencialmente verticales, y Cww está inclinada a unos 35°. Summary This thesis studies the similitudes and differences between external and internal in¬compressible wall-bounded turbulent flows at moderately-high Reynolds numbers. We consider numerical and experimental zero-pressure-gradient boundary layers and chan¬nels in the range of δ+ ~ 500 — 2000. These shear flows are subjects of intensive research because of their technological importance and fundamental physical interest. However, there are still open questions regarding basic aspects such as the universality of the mean and fluctuating velocity and pressure profiles at the near-wall and logarithmic regions, their scaling and the effect of the Reynolds numbers, or the differences between internal and external flows at the outer layer, to name but a few. For this study, we made use of available direct numerical simulations of channel and boundary layers reaching δ+ ~ 2000 and δ+ ~ 700, respectively. To fill the gap in the Reynolds number, a new boundary layer simulation in the range δ+ ~ 1000-2000 is presented and discussed. The original raw data and the post-processed statistics are publicly available on our website.162 The analysis of the one-point statistic confirms the existence of logarithmic profiles for the spanwise w'2+ and pressure p'2+ fluctuations for both type of flows, but not for the wall-normal v'2+ or the streamwise u'2+ velocities. To accept or reject the existence of a logarithmic range in u'2+ requires higher Reynolds numbers than the ones considered in this work. An important consequence of having such profiles is that the maximum value of the intensities, reached near the wall, depends on the Reynolds number. This was confirmed after surveying a wide number of experimental and numerical datasets, corrob¬orating that the maximum of ul2+, p'2+, and w'2+ increases proportionally to log(δ+). On the other hand, that maximum is more intense in external flows than in internal ones, differing the most around y/δ ~ 0.3-0.5, and essentially independent of the Reynolds number. We discuss that those differences are originated as a consequence of the inter¬mittent character of boundary layers that is absent in internal flows. The structure of the velocity and pressure fluctuations, together with those of the Reynolds shear stress, were investigated using three-dimensional two-point spatial correlations. We find that the correlations extend over longer distances in channels than in boundary layers, especially in the case of the streamwise correlation Cuu in the flow direc-tion. For weakly correlated structures, the maximum streamwise length of Cuu is O(78) for boundary layers and O(188) for channels, attained at the logarithmic and outer regions respectively. The corresponding lengths for the transverse velocities and for the pressure are shorter, 0(8 — 28), and of the same order for both flows. The spatial organization of the velocity correlations is shown to be consistent with a pair of quasi-streamwise rollers that scales in outer units. That organization is observed until y ~ 0.68, from which boundary layers start to organize into spanwise rollers. This effect is weaker in channels, and it appears at y ~ 0.88. We present correlations conditioned to intense events of the transversal velocity, w', to study if these structures meander along the spanwise direction. The results indicate that the streamwise velocity streaks increase their meandering proportionally to the distance to the wall, suggesting that the structures are more aligned close to the wall than far from it. The reason behind this meandering is probably due to the characteristic organization along diagonals of w'. These structures not only meander along the spanwise direction, but they are also inclined to the wall at angles that depend on the distance from the wall, on the variable being considered, and on the correlation level used to define them. Above the buffer layer and independent of the Reynolds numbers and type of flow, the maximum inclination of Cuu is about 10°, Cvv and Cpp are roughly vertical, and Cww is inclined by 35°.

Numerical analysis of deflagration initation by a hot jet

Numerical simulations of axisymmetric reactive jets with one-step Arrhenius kinetics are used to investigate the problem of deflagration initiation in a premixed fuel–air mixture by the sudden discharge of a hot jet of its adiabatic reaction products. For the moderately large values of the jet Reynolds number considered in the computations, chemical reaction is seen to occur initially in the thin mixing layer that separates the hot products from the cold reactants. This mixing layer is wrapped around by the starting vortex, thereby enhancing mixing at the jet head, which is followed by an annular mixing layer that trails behind, connecting the leading vortex with the orifice rim. A successful deflagration is seen to develop for values of the orifice radius larger than a critical value a c in the order of the flame thickness of the planar deflagration δL. Introduction of appropriate scales provides the dimensionless formulation of the problem, with flame initiation characterised in terms of a critical Damköhler number Δc=(a d/δL)2, whose parametric dependence is investigated. The numerical computations reveal that, while the jet Reynolds number exerts a limited influence on the criticality conditions, the effect of the reactant diffusivity on ignition is much more pronounced, with the value of Δc increasing significantly with increasing Lewis numbers. The reactant diffusivity affects also the way ignition takes place, so that for reactants with the flame develops as a result of ignition in the annular mixing layer surrounding the developing jet stem, whereas for highly diffusive reactants with Lewis numbers sufficiently smaller than unity combustion is initiated in the mixed core formed around the starting vortex. The analysis provides increased understanding of deflagration initiation processes, including the effects of differential diffusion, and points to the need for further investigations corporating detailed chemistry models for specific fuel–air mixtures.

Transverse galloping at low Reynolds numbers

The possibility of transverse galloping of a square cylinder at low Reynolds numbers (Re≤200Re≤200, so that the flow is presumably laminar) is analysed. Transverse galloping is here considered as a one-degree-of-freedom oscillator subjected to fluid forces, which are described by using the quasi-steady hypothesis (time-averaged data are extracted from previous numerical simulations). Approximate solutions are obtained by means of the method of Krylov-Bogoliubov, with two major conclusions: (i) a square cylinder cannot gallop below a Reynolds number of 159 and (ii) in the range 159≤Re≤200159≤Re≤200 the response exhibits no hysteresis.

Estudo numérico do efeito de blocagem hidrodinâmica em um túnel de cavitação sobre o escoamento ao redor de um submarino.

O projeto de submarinos é um tema que vem se desenvolvendo há muitos anos. No Brasil, vivemos um momento de pleno avanço na área, sendo encabeçado pela Marinha do Brasil com o desenvolvimento do primeiro submarino de propulsão nuclear brasileiro. No decorrer desse projeto, diversos estudos devem ser desenvolvidos, dentre eles, ensaios com modelos em escala, que serão realizados em um túnel de cavitação. Durante estes experimentos, espera-se que surjam diversas fontes de incerteza, que podem comprometer os resultados obtidos, destacando-se o efeito de blocagem hidrodinâmica. Nesse contexto, esta dissertação apresenta uma investigação do fenômeno de blocagem no modelo de submarino DARPA SUBOFF na seção de testes de um túnel de cavitação. Para isso, simulações tridimensionais de média de Reynolds (RANS) foram realizadas para o modelo de submarino posicionado no túnel (sob blocagem) e sob condição de águas profundas (referência), para diferentes ângulos de ataque e número de Reynolds (Re) constante. Coeficientes de força e momento hidrodinâmicos foram obtidos e os métodos de correção de blocagem clássico e WICS foram empregados para obter fatores que pudessem corrigir esses valores. A eficiência de cada método foi avaliada comparando os valores corrigidos com os de referência. Os resultados mostraram que ambos métodos de correção são eficientes na recuperação de coeficientes de força e momento hidrodinâmicos para razões de blocagem inferiores a 10%, mas ressalvas são feitas quanto à utilização de cada um desses métodos, sendo o método clássico indicado para o estudo de geometrias simples e o método WICS para modelos mais complexos. Paralelamente, a topologia do escoamento ao redor do submarino foi estudada através do método dos pontos críticos, para avaliar se a blocagem hidrodinâmica pode causar alterações nas estruturas presentes no escoamento. A análise topológica permitiu avaliar a evolução nas estruturas do escoamento para os ângulos de ataque e Re simulados. Foram capturadas as linhas de separação no corpo e avaliada sua relação com a tensão de cisalhamento ao longo do casco. A blocagem hidrodinâmica não causou alterações drásticas na topologia do escoamento.

M=3 turbulent boundary layer measurements at very high Reynolds numbers /

"Interim report for period January 1975 to September 1976."

Constant shear along a boundary of a flume

'I'he accurate rreasurement of bed shear stress has been extremely difficult due to its changing values until white propunded a theory which would give constant shear along the bed of a flume. In this investigation a flume has been designed according to White's theory and by two separate methods proven to give constant shearing force along the bed. The first method applied the Hydrogen Bubble Technique to obtain accurate values of velocity thus allowing the velocity profile to be plotted and the momentum at the various test sections to be calculated. The use of a 16 mm Beaulieu movie camera allowed the exact velocity profiles created by the hydrogen bubbles to be recorded whilst an analysing projector gave the means of calculating the exact velocities at the various test sections. Simultaneously Preston's technique of measuring skin friction using Pitot tubes was applied. Twc banks of open ended water manometer were used for recording the static and velocity head pressure drop along the flume. This tvpe of manometer eliminated air locks in the tubes and was found to be sufficiently accurate. Readings of pressure and velocity were taken for various types and diameters of bed material both natural sands and glass spheres and the results tabulated. Graphs of particle Reynolds Number against bed shear stress were plotted and gave a linear relationship which dropped off at high values of Reynolds number. It was found that bed movement occurred instantaneously along the bed of the flume once critical velocity had been reached. On completion of this test a roof curve inappropriate to the bed material was used and then the test repeated. The bed shearing stress was now no longer constant and yet bed movement started instantaneously along the bed of the flume, showing that there are more parameters than critical shear stress to bed movement. It is concluded from the two separate methods applied that the bed shear stress is constant along the bed of the flume.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.

Evidences of turbulent mixing in multi-pumping flow systems

Multi-pumping flow systems exploit pulsed flows delivered by Solenoid pumps. Their improved performance rely on the enhanced radial mass transport inherent to the pulsed flow, which is a consequence of the establishment of vortices thus a tendency towards turbulent mixing. This paper presents several evidences of turbulent mixing in relation to pulsed flows. such as recorded peak shape, establishment of fluidized beds, exploitation of flow reversal, implementation of relatively slow chemical reactions and/or heating of the reaction medium. In addition, Reynolds number associated with the GO period of a pulsed flow is estimated and photographic images of dispersing samples flowing under laminar regime and pulsed flow conditions are presented. (C) 2009 Elsevier B.V. All rights reserved.

Estimation of Mass Transfer Velocity Based on Measured Turbulence Parameters

The aim of this study is to quantify the mass transfer velocity using turbulence parameters from simultaneous measurements of oxygen concentration fields and velocity fields. The surface divergence model was considered in more detail, using data obtained for the lower range of beta (surface divergence). It is shown that the existing models that use the divergence concept furnish good predictions for the transfer velocity also for low values of beta, in the range of this study. Additionally, traditional conceptual models, such as the film model, the penetration-renewal model, and the large eddy model, were tested using the simultaneous information of concentration and velocity fields. It is shown that the film and the surface divergence models predicted the mass transfer velocity for all the range of the equipment Reynolds number used here. The velocity measurements showed viscosity effects close to the surface, which indicates that the surface was contaminated with some surfactant. Considering the results, this contamination can be considered slight for the mass transfer predictions. (C) 2009 American Institute of Chemical Engineers AIChE J, 56: 2005-2017; 2010

Macro-Fiber Composite Actuators for Flow Control of a Variable Camber Airfoil

This research employs solid-state actuators for delay of flow separation seen in airfoils at low Reynolds numbers. The flow control technique investigated here is aimed for a variable camber airfoil that employs two active surfaces and a single four-bar (box) mechanism as the internal structure. To reduce separation, periodic excitation to the flow around the leading edge of the airfoil is induced by a total of nine piezocomposite actuated clamped-free unimorph benders distributed in the spanwise direction. An electromechanical model is employed to design an actuator capable of high deformations at the desired frequency for lift improvement at post-stall angles. The optimum spanwise distribution of excitation for increasing lift coefficient is identified experimentally in the wind tunnel. A 3D (non-uniform) excitation distribution achieved higher lift enhancement in the post-stall region with lower power consumption when compared to the 2D (uniform) excitation distribution. A lift coefficient increase of 18.4% is achieved with the identified non-uniform excitation mode at the bender resonance frequency of 125 Hz, the flow velocity of 5 m/s and at the reduced frequency of 3.78. The maximum lift (Clmax) is increased 5.2% from the baseline. The total power consumption of the flow control technique is 639 mW(RMS).

Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m(3)) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe. (C) 2010 Elsevier Inc. All rights reserved.

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex-structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Effect of the domain spanwise periodic length on the flow around a circular cylinder

We assess the effect of the choice of spanwise periodic length on simulations of the flow around a fixed circular cylinder. The Reynolds number is set to 400 because, at this value, both lift coefficient and shedding frequency show significant drop due to three-dimensional flow structures. From the analysis of the three-dimensionalities of the wake and of the integral quantities such as Strouhal number, RMS of lift coefficient and energy contained in the three-dimensional portion of the flow we obtain an estimate of the minimum spanwise length to satisfactorily represent the flow. Furthermore, we observe a distinct wake behavior when the spanwise length is approximately the mode B instability wavelength. (C) 2011 Elsevier Ltd. All rights reserved.