Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

Large-eddy simulations of fully developed turbulent channel and pipe flows with smooth and rough walls

Studies in turbulence often focus on two flow conditions, both of which occur frequently in real-world flows and are sought-after for their value in advancing turbulence theory. These are the high Reynolds number regime and the effect of wall surface roughness. In this dissertation, a Large-Eddy Simulation (LES) recreates both conditions over a wide range of Reynolds numbers Re_τ = O(10²)-O(10⁸) and accounts for roughness by locally modeling the statistical effects of near-wall anisotropic fine scales in a thin layer immediately above the rough surface. A subgrid, roughness-corrected wall model is introduced to dynamically transmit this modeled information from the wall to the outer LES, which uses a stretched-vortex subgrid-scale model operating in the bulk of the flow. Of primary interest is the Reynolds number and roughness dependence of these flows in terms of first and second order statistics. The LES is first applied to a fully turbulent uniformly-smooth/rough channel flow to capture the flow dynamics over smooth, transitionally rough and fully rough regimes. Results include a Moody-like diagram for the wall averaged friction factor, believed to be the first of its kind obtained from LES. Confirmation is found for experimentally observed logarithmic behavior in the normalized stream-wise turbulent intensities. Tight logarithmic collapse, scaled on the wall friction velocity, is found for smooth-wall flows when Re_τ ≥ O(10⁶) and in fully rough cases. Since the wall model operates locally and dynamically, the framework is used to investigate non-uniform roughness distribution cases in a channel, where the flow adjustments to sudden surface changes are investigated. Recovery of mean quantities and turbulent statistics after transitions are discussed qualitatively and quantitatively at various roughness and Reynolds number levels. The internal boundary layer, which is defined as the border between the flow affected by the new surface condition and the unaffected part, is computed, and a collapse of the profiles on a length scale containing the logarithm of friction Reynolds number is presented. Finally, we turn to the possibility of expanding the present framework to accommodate more general geometries. As a first step, the whole LES framework is modified for use in the curvilinear geometry of a fully-developed turbulent pipe flow, with implementation carried out in a spectral element solver capable of handling complex wall profiles. The friction factors have shown favorable agreement with the superpipe data, and the LES estimates of the Karman constant and additive constant of the log-law closely match values obtained from experiment.

Differential equation based length scales to improve DES and RANS simulations

Surprisingly expensive to compute wall distances are still used in a range of key turbulence and peripheral physics models. Potentially economical, accuracy improving differential equation based distance algorithms are considered. These involve elliptic Poisson and hyperbolic natured Eikonal equation approaches. Numerical issues relating to non-orthogonal curvilinear grid solution of the latter are addressed. Eikonal extension to a Hamilton-Jacobi (HJ) equation is discussed. Use of this extension to improve turbulence model accuracy and, along with the Eikonal, enhance Detached Eddy Simulation (DES) techniques is considered. Application of the distance approaches is studied for various geometries. These include a plane channel flow with a wire at the centre, a wing-flap system, a jet with co-flow and a supersonic double-delta configuration. Although less accurate than the Eikonal, Poisson method based flow solutions are extremely close to those using a search procedure. For a moving grid case the Poisson method is found especially efficient. Results show the Eikonal equation can be solved on highly stretched, non-orthogonal, curvilinear grids. A key accuracy aspect is that metrics must be upwinded in the propagating front direction. The HJ equation is found to have qualitative turbulence model improving properties. © 2003 by P. G. Tucker.

On the laminar-turbulent transition of the rotating-disk flow: The role of absolute instability

© 2014 Cambridge University Press. This paper describes a detailed experimental study using hot-wire anemometry of the laminar-turbulent transition region of a rotating-disk boundary-layer flow without any imposed excitation of the boundary layer. The measured data are separated into stationary and unsteady disturbance fields in order to elaborate on the roles that the stationary and the travelling modes have in the transition process. We show the onset of nonlinearity consistently at Reynolds numbers, R, of ∼ 510, i.e. at the onset of Lingwood's (J. Fluid Mech., vol. 299, 1995, pp. 17-33) local absolute instability, and the growth of stationary vortices saturates at a Reynolds number of ∼ 550. The nonlinear saturation and subsequent turbulent breakdown of individual stationary vortices independently of their amplitudes, which vary azimuthally, seem to be determined by well-defined Reynolds numbers. We identify unstable travelling disturbances in our power spectra, which continue to grow, saturating at around R=585, whereupon turbulent breakdown of the boundary layer ensues. The nonlinear saturation amplitude of the total disturbance field is approximately constant for all considered cases, i.e. different rotation rates and edge Reynolds numbers. We also identify a travelling secondary instability. Our results suggest that it is the travelling disturbances that are fundamentally important to the transition to turbulence for a clean disk, rather than the stationary vortices. Here, the results appear to show a primary nonlinear steep-fronted (travelling) global mode at the boundary between the local convectively and absolutely unstable regions, which develops nonlinearly interacting with the stationary vortices and which saturates and is unstable to a secondary instability. This leads to a rapid transition to turbulence outward of the primary front from approximately R=565 to 590 and to a fully turbulent boundary layer above 650.

Molecular dynamics simulation investigation on the plastic flow behaviour of silicon during nanometric cutting

Molecular dynamics (MD) simulation was carried out to acquire an in-depth understanding of the flow behaviour of single crystal silicon during nanometric cutting on three principal crystallographic planes and at different cutting temperatures. The key findings were that (i) the substrate material underneath the cutting tool was observed for the first time to experience a rotational flow akin to fluids at all the tested temperatures up to 1200 K. (ii) The degree of flow in terms of vorticity was found higher on the (1 1 1) crystal plane signifying better machinability on this orientation in accord with the current pool of knowledge (iii) an increase in the machining temperature reduces the springback effect and thereby the elastic recovery and (iv) the cutting orientation and the cutting temperature showed significant dependence on the location of the stagnation region in the cutting zone of the substrate.

Experimental analysis of air flow profiles in a double skin facade in a maritime climate

Glazed Double Skin Facades (DSF) offer the potential to improve the performance of all-glass building skins, common to commercial office buildings in which full facade glazing has almost become the standard. Single skin glazing results in increased heating and cooling costs over opaque walls, due to lower thermal resistance of glass, and the increased impact of solar gain through it. However, the performance benefit of DSF technology continues to be questioned and its operation poorly understood, particularly the nature of airflow through the cavity. This paper deals specifically with the experimental analysis of the air flow characteristics in an automated double skin façade. The benefit of the DSF as a thermal buffer, and to limit overheating is evaluated through analysis of an extensive set of parameters including air and surface temperatures at each level in the DSF, airflow readings in the cavity and at the inlet and outlet, solar and wind data, and analytically derived pressure differentials. The temperature and air-flow are monitored in the cavity of a DSF using wireless sensors and hot wire anemometers respectively. Automated louvre operation and building set-points are monitored via the BMS. Thermal stratification and air flow variation during changing weather conditions are shown to effect the performance of the DSF considerably and hence the energy performance of the building. The relative pressure effects due to buoyancy and wind are analysed and quantified. This research aims to developed and validate models of DSFs in the maritime climate, using multi-season data from experimental monitoring. This extensive experimental study provides data for training and validation of models.

Análise da indução de fluxo de ar por convecção livre em chaminé solar

Experiments were performed to study the effect of surface properties of a vertical channel heated by a source of thermal radiation to induce air flow through convection. Two channels (solar chimney prototype) were built with glass plates, forming a structure of truncated pyramidal geometry. We considered two surface finishes: transparent and opaque. Each stack was mounted on a base of thermal energy absorber with a central opening for passage of air, and subjected to heating by a radiant source comprises a bank of incandescent bulbs and were performed field tests. Thermocouples were fixed on the bases and on the walls of chimneys and then connected to a data acquisition system in computer. The air flow within the chimney, the speed and temperature were measured using a hot wire anemometer. Five experiments were performed for each stack in which convective flows were recorded with values ranging from 17 m³ / h and 22 m³ / h and air flow velocities ranging from 0.38 m / s and 0.56 m / s for the laboratory tests and air velocities between 0.6 m/s and 1.1m/s and convective airflows between 650 m³/h and 1150 m³/h for the field tests. The test data were compared to those obtained by semi-empirical equations, which are valid for air flow induced into channels and simulated data from 1st Thermodynamics equation. It was found that the chimney with transparent walls induced more intense convective flows than the chimney with matte finish. Based on the results obtained can be proposed for the implementation of prototype to exhaust fumes, mists, gases, vapors, mists and dusts in industrial environments, to help promote ventilation and air renewal in built environments and for drying materials, fruits and seeds

CH and C2 radicals characterization in natural gas turbulent diffusion flames

This paper reports the construction of an axisymmetric nonpremixed piloted jet burner, with well-defined initial and boundary conditions, known as the Delft burner, to assess turbulence-chemistry interaction in non-premixed turbulent flames. Detailed experimental information is described, involving hot-wire anemometry, thin-wire thermocouples and chemiluminescence visualization measurements. Radial profile of the axial mean velocity indicates excellent agreement between flow patterns developed within Delft installation and the one described herein. Chemiluminescence emissions from CH and C2 free-radicals were acquired with a CCD camera. Tomography reconstruction analysis was utilised to compare radical emissions and temperature spatial distributions. There was a strong dependence between temperature and CH/C 2 emissions. This is an indication that these radicals can be used in flame front studies.

Second law analysis and simulation techniques for the energy optimization of buildings

The research activity described in this thesis is focused mainly on the study of finite-element techniques applied to thermo-fluid dynamic problems of plant components and on the study of dynamic simulation techniques applied to integrated building design in order to enhance the energy performance of the building. The first part of this doctorate thesis is a broad dissertation on second law analysis of thermodynamic processes with the purpose of including the issue of the energy efficiency of buildings within a wider cultural context which is usually not considered by professionals in the energy sector. In particular, the first chapter includes, a rigorous scheme for the deduction of the expressions for molar exergy and molar flow exergy of pure chemical fuels. The study shows that molar exergy and molar flow exergy coincide when the temperature and pressure of the fuel are equal to those of the environment in which the combustion reaction takes place. A simple method to determine the Gibbs free energy for non-standard values of the temperature and pressure of the environment is then clarified. For hydrogen, carbon dioxide, and several hydrocarbons, the dependence of the molar exergy on the temperature and relative humidity of the environment is reported, together with an evaluation of molar exergy and molar flow exergy when the temperature and pressure of the fuel are different from those of the environment. As an application of second law analysis, a comparison of the thermodynamic efficiency of a condensing boiler and of a heat pump is also reported. The second chapter presents a study of borehole heat exchangers, that is, a polyethylene piping network buried in the soil which allows a ground-coupled heat pump to exchange heat with the ground. After a brief overview of low-enthalpy geothermal plants, an apparatus designed and assembled by the author to carry out thermal response tests is presented. Data obtained by means of in situ thermal response tests are reported and evaluated by means of a finite-element simulation method, implemented through the software package COMSOL Multyphysics. The simulation method allows the determination of the precise value of the effective thermal properties of the ground and of the grout, which are essential for the design of borehole heat exchangers. In addition to the study of a single plant component, namely the borehole heat exchanger, in the third chapter is presented a thorough process for the plant design of a zero carbon building complex. The plant is composed of: 1) a ground-coupled heat pump system for space heating and cooling, with electricity supplied by photovoltaic solar collectors; 2) air dehumidifiers; 3) thermal solar collectors to match 70% of domestic hot water energy use, and a wood pellet boiler for the remaining domestic hot water energy use and for exceptional winter peaks. This chapter includes the design methodology adopted: 1) dynamic simulation of the building complex with the software package TRNSYS for evaluating the energy requirements of the building complex; 2) ground-coupled heat pumps modelled by means of TRNSYS; and 3) evaluation of the total length of the borehole heat exchanger by an iterative method developed by the author. An economic feasibility and an exergy analysis of the proposed plant, compared with two other plants, are reported. The exergy analysis was performed by considering the embodied energy of the components of each plant and the exergy loss during the functioning of the plants.

On the onset of turbulence in natural convection on inclined plates

The problem of determination of the turbulence onset in natural convection on heated inclined plates in an air environment has been experimentally revisited. The transition has been detected by using hot wire velocity measurements. The onset of turbulence has been considered to take place where velocity fluctuations (measured through turbulence intensity) start to grow. Experiments have shown that the onset depends not only on the Grashof number defined in terms of the temperature difference between the heated plate and the surrounding air. A correlation between dimensionless Grashof and Reynolds numbers has been obtained, fitting quite well the experimental data.

Effect of cold water immersion on repeated cycling performance and limb blood flow

The purpose of the present study was to compare the effects of cold water immersion (CWI) and active recovery (ACT) on resting limb blood flow, rectal temperature and repeated cycling performance in the heat. Ten subjects completed two testing sessions separated by 1 week; each trial consisted of an initial all-out 35-min exercise bout, one of two 15-min recovery interventions (randomised: CWI or ACT), followed by a 40-min passive recovery period before repeating the 35-min exercise bout. Performance was measured as the change in total work completed during the exercise bouts. Resting limb blood flow, heart rate, rectal temperature and blood lactate were recorded throughout the testing sessions. There was a significant decline in performance after ACT (mean (SD) −1.81% (1.05%)) compared with CWI where performance remained unchanged (0.10% (0.71%)). Rectal temperature was reduced after CWI (36.8°C (1.0°C)) compared with ACT (38.3°C (0.4°C)), as was blood flow to the arms (CWI 3.64 (1.47) ml/100 ml/min; ACT 16.85 (3.57) ml/100 ml/min) and legs (CW 4.83 (2.49) ml/100 ml/min; ACT 4.83 (2.49) ml/100 ml/min). Leg blood flow at the end of the second exercise bout was not different between the active (15.25 (4.33) ml/100 ml/min) and cold trials (14.99 (4.96) ml/100 ml/min), whereas rectal temperature (CWI 38.1°C (0.3°C); ACT 38.8°C (0.2°C)) and arm blood flow (CWI 20.55 (3.78) ml/100 ml/min; ACT 23.83 (5.32) ml/100 ml/min) remained depressed until the end of the cold trial. These findings indicate that CWI is an effective intervention for maintaining repeat cycling performance in the heat and this performance benefit is associated with alterations in core temperature and limb blood flow.

Transition from radiatively inefficient to cooling dominated phase in two temperature accretion disks around black holes

We investigate the transition of a radiatively inefficient phase of a viscous two temperature accreting flow to a cooling dominated phase and vice versa around black holes. Based on a global sub-Keplerian accretion disk model in steady state, including explicit cooling processes self-consistently, we show that general advective accretion flow passes through various phases during its infall towards a black hole. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. Hence the flow governs a much lower electron temperature similar to 10(8) - 10(9.5) K compared to the hot protons of temperature similar to 10(10.2) - 10(11.8) K in the range of the accretion rate in Eddington units 0.01 less than or simiar to (M) over dot less than or similar to 100. Therefore, the solutions may potentially explain the hard X-rays and the gamma-rays emitted from AGNs and X-ray binaries. We finally compare the solutions for two different regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling dominated compared to its highly viscous counterpart which is radiatively inefficient. The flow is successfully able to reproduce the observed minosities of the under-fed AGNs and quasars (e.g. Sgr A*), ultra-luminous X-ray sources (e.g. SS433), as well as the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67) at different combinations of the mass accretion rate and ratio of specific heats.