Kelvin-Helmholtz-like instability of turbulent flows over riblets

The turbulent drag reduction due to riblets is a function of their size and, for different configurations, collapses well with a length scale l+g=(A+g)1/2, based in the groove cross-section Ag. The initially linear drag reduction breaks down for l+g≈11, which agrees in our DNS with the previously reported appearance of quasi-two-dimensional spanwise rollers immediately above the riblets. They are similar to those found over porous surfaces and plant canopies, and can be traced to a Kelvin-Helmholtz-like instability associated with the relaxation of the impermeability condition for the wall-normal velocity. The extra Reynolds stress associated with them accounts quantitatively for the drag degradation. An inviscid model for the instability confirms its nature, agreeing well with the observed perturbation wavelengths and shapes. The onset of the instability is determined by a length scale L+w that, for conventional riblet geometries, is proportional to l+g. The instability onset, L+w≥4, corresponds to the empirical breakdown point l+g≈11.

Study of turbulent heat transfer of aviation kerosene flows in a curved pipe at supercritical pressure

The now and heat transfer characteristics of China No. 3 aviation kerosene in a heated curved tube under supercritical pressure are numerically investigated by a finite volume method. A two-layer turbulence model, consisting of the RNG k-epsilon two-equation model and the Wolfstein one-equation model, is used for the simulation of turbulence. A 10-species kerosene surrogate model and the NIST Supertrapp software are applied to obtain the thermophysical and transport properties of the kerosene at various temperature under a supercritical pressure of 4 MPa. The large variation of thermophysical properties of the kerosene at the supercritical pressure make the flow and heat transfer more complicated, especially under the effects of buoyancy and centrifugal force. The centrifugal force enhances the heat transfer, but also increases the friction factors. The rise of the velocity caused by the variation of the density does not enhance the effects of the centrifugal force when the curvature ratios are less than 0.05. On the contrary, the variation of the density increases the effects of the buoyancy. (C) 2010 Elsevier Ltd. All rights reserved.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

The three-dimensional vortical nature of atmospheric and oceanic turbulent flows

Using a novel numerical method at unprecedented resolution, we demonstrate that structures of small to intermediate scale in rotating, stratified flows are intrinsically three-dimensional. Such flows are characterized by vortices (spinning volumes of fluid), regions of large vorticity gradients, and filamentary structures at all scales. It is found that such structures have predominantly three-dimensional dynamics below a horizontal scale LLR, where LR is the so-called Rossby radius of deformation, equal to the characteristic vertical scale of the fluid H divided by the ratio of the rotational and buoyancy frequencies f/N. The breakdown of two-dimensional dynamics at these scales is attributed to the so-called "tall-column instability" [D. G. Dritschel and M. de la Torre Juárez, J. Fluid. Mech. 328, 129 (1996)], which is active on columnar vortices that are tall after scaling by f/N, or, equivalently, that are narrow compared with LR. Moreover, this instability eventually leads to a simple relationship between typical vertical and horizontal scales: for each vertical wave number (apart from the vertically averaged, barotropic component of the flow) the average horizontal wave number is equal to f/N times the vertical wave number. The practical implication is that three-dimensional modeling is essential to capture the behavior of rotating, stratified fluids. Two-dimensional models are not valid for scales below LR. ©1999 American Institute of Physics.

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic - Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.

Pressure distribution on the frontal disk for turbulent flows in a radial diffuser

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Magnetic flux transport by turbulent reconnection in astrophysical flows

The role of magnetohydrodynamics (MHD) turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. Otherwise, the newborn stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disc that grows around the protostar. The efficiency of the mechanism that is often invoked to allow transport of magnetic fields at different stages of star formation, namely ambipolar diffusion, has recently been put in check. We discuss here an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results from three-dimensional MHD numerical simulations that indicate that this mechanism is very efficient in decoupling and transporting magnetic flux from the inner denser regions to the outskirts of collapsing clouds at different stages of star formation. We discuss this mechanism also in the context of dynamo processes and speculate that it can play a role both in solar dynamo and in accretion disc dynamo processes.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Length of the annular regime for condensing flows inside a horizontal channel - the experimental determination of its values and its trends

The experiments observe and measure the length of the annular regime in fully condensing quasi-steady (steady-in-the-mean) flows of pure FC-72 vapor in a horizontal condenser (rectangular cross-section of 2 mm height, 15 mm width, and 1 m length). The sides and top of the duct are made of clear plastic that allows flow visualization. The experimental system in which this condenser is used is able to control and achieve different quasi-steady mass flow rates, inlet pressures, and wall cooling conditions (by adjustment of the temperature and flow rate of the cooling water flowing underneath the condensing-plate). The reported correlations and measurements for the annular length are also vital information for determining the length of the annular regime and proposing extended correlation (covering many vapors and a larger parameter set than the experimentally reported version here) by ongoing independent modeling and computational simulation approach.

Laminar and turbulent nozzle-jet flows and their acoustic near-field

We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re D = 18 100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and the sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.

Numerical characterization of particle dispersion in the turbulent recirculation zones of sudden expansion pipe flows

The dispersion of solid particles in the turbulent recirculation zones of sudden expansion pipes can be characterized by different Stokes numbers and mean drift parameter and its study is important because this kind of flows appears in many technological applications.

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°.