Relatively few studies have been carried out on the turbulent wake structure of a finite circular cylinder and a stack partially immersed in a flat-plate turbulent boundary layer. There is a need to develop a better understanding of the wakes of these structures, since they have many important engineering applications. This thesis investigates the influence of the aspect ratio on the wake of a finite circular cylinder and the effects of the ratio of jet flow velocity to crossflow velocity (velocity ratio, R) on the wake of a stack in a cross-flow. <p>The wake characteristics of flows over a finite circular cylinder at four different aspect ratios (AR = 3, 5, 7 and 9) were investigated experimentally at a Reynolds number of ReD = 6104 using two-component thermal anemometry. Each cylinder was mounted normal to a ground plane and was either completely or partially immersed in a flat-plate turbulent boundary layer. The ratio of boundary layer thickness to the cylinder diameter was 3. <p>A similar turbulent wake structure (time-averaged velocity, turbulence intensity, and Reynolds shear stress distributions) was found for the cylinders with AR = 5, 7, and 9, while a distinctly different turbulent wake structure was found for the cylinder with AR = 3. This was consistent with the results of a previous study that focused on the time-averaged streamwise vortex structures in the wake. In addition, irrespective of the value of AR, high values were observed for the skewness and flatness factors around the free end of the cylinders, which may be attributed to the interaction of the tip vortex structures and downwash flow that dominates this region of the cylinder.<p>The wake characteristics of a stack of aspect ratio AR = 9 were investigated using both the seven-hole pressure probe and thermal anemometry. The seven-hole probe was used to measure the three components of the time-averaged velocity field, while the thermal anemometry was used to measure two components of the turbulent velocity field at various downstream locations from the stack. The stack was mounted normal to the ground plane and was partially immersed in a flat-plate turbulent boundary layer, for which the ratio of boundary layer thickness to the stack diameter was 4.5. In addition, measurements of the vortex shedding frequency were made with a single-component hot-wire probe. The cross-flow Reynolds number was ReD = 2.3 x 104, the jet Reynolds number ranged from Red = 7.6 x 103 to 4.7 x 104, and R was varied from 0 to 3. <p>In the stack study, three flow regimes were identified depending on the value of R: the downwash (R < 0.7), cross-wind-dominated (0.7 < R < 1.5), and jet-dominated (R ≥ 1.5) flow regimes. Each flow regime had a distinct structure for the time-averaged velocity and streamwise vorticity fields, and turbulence characteristics, as well as the variation of the Strouhal number and the power spectrum of the streamwise velocity fluctuations along the stack height. The turbulence structure is complex and changes in the streamwise and wall-normal directions within the near and intermediate stack and jet wakes. In the downwash and crosswind-dominated flow regimes, two pairs of counter-rotating streamwise vortex structures were identified within the stack wake. The tip-vortex pair and base-vortex pair were similar to those found in the wake of a finite circular cylinder, located close to the free end and the base of the stack (ground plane), respectively. In the jet-dominated flow regime, a third pair of streamwise vortex structures was observed, referred to as the jet-wake vortex pair, which occurred within the jet-wake region above the free end of the stack. The jet-wake vortex pair has the same orientation as the base vortex pair and is associated with the jet rise.

We explore the capability of lattice Boltzmann equation (LBE) method for complex fluid flows involving turbulence, mixing, and reaction. In the first study, LBE schemes for binary scalar mixing and multi-component reacting flow with reactions are developed. Simulations of initially non-premixed mixtures yield scalar probability distribution functions that are in good agreement with numerical data obtained from Navier-Stokes (NS) equation based computation. One-dimensional chemically-reacting flow simulation of a premixed mixture yields a flame speed that is consistent with experimentally determined value. The second study involves direct numerical simulation (DNS) and large-eddy simulation (LES) of decaying homogenous isotropic turbulence (HIT) with and without frame rotation. Three categories of simulations are performed: (i) LBE-DNS in both inertial and rotating frames; (ii) LBE-LES in inertial frame; (iii) Comparison of the LBE-LES vs. NS-LES. The LBE-DNS results of the decay exponents for kinetic energy k and dissipation rate &#949;, and the low wave-number scaling of the energy spectrum agree well with established classical results. The LBE-DNS also captures rotating turbulence physics. The LBE-LES accurately captures low-wave number scaling, energy decay and large scale structures. The comparisons indicate that the LBE-LES simulations preserve flow structures somewhat more accurately than the NS-LES counterpart. In the third study, we numerically investigate the near-field mixing features in low aspect-ratio (AR) rectangular turbulent jets (RTJ) using the LBE method. We use D3Q19 multiple-relaxation-time (MRT) LBE incorporating a subgrid Smagorinsky model for LES. Simulations of four jets which characterized by AR, exit velocity, and Reynolds number are performed. The investigated near-field behaviors include: (1) Decay of mean streamwise velocity (MSV) and inverse MSV; (2) Spanwise and lateral profiles of MSV; (3) Half-velocity width development and MSV contours; and (4) Streamwise turbulence intensity distribution and spanwise profiles of streamwise turbulence intensity. The computations are compared against experimental data and the agreement is good. We capture both unique features of RTJ: the saddle-back spanwise profile of MSV and axis-switching of long axis from spanwise to lateral direction. Overall, this work serves to establish the feasibility of the LBE method as a viable tool for computing mixing, combustion, and turbulence.

abstract: Assessments for the threats posed by volcanic eruptions rely in large part on the accurate prediction of volcanic plume motion over time. That predictive capacity is currently hindered by a limited understanding of volcanic plume dynamics. While eruption rate is considered a dominant control on volcanic plume dynamics, the effects of variable eruption rates on plume rise and evolution are not well understood. To address this aspect of plume dynamics, I conducted an experimental investigation wherein I quantified the relationship between laboratory jet development and highly-variable discharge rates under conditions analogous to those which may prevail in unsteady, short-lived explosive eruptions. I created turbulent jets in the laboratory by releasing pressurized water into a tank of still water. I then measured the resultant jet growth over time using simple video images and particle image velocimetry (PIV). I investigated jet behavior over a range of jet Reynolds numbers which overlaps with estimates of Reynolds numbers for short-duration volcanic plumes. By analysis of the jet boundary and velocity field evolution, I discovered a direct relationship between changes in vent conditions and jet evolution. Jet behavior evolved through a sequence of three stages - jet-like, transitional, and puff-like - that correlate with three main injection phases - acceleration, deceleration and off. While the source was off, jets were characterized by relatively constant internal velocity distributions and flow propagation followed that of a classical puff. However, while the source was on, the flow properties - both in the flows themselves and in the induced ambient flow - changed abruptly with changes at the source. On the basis of my findings for unsteady laboratory jets, I conclude that variable eruption rates with characteristic time scales close to eruption duration have first-order control over volcanic plume evolution. Prior to my study, the significance of this variation was largely uncharacterized as the volcanology community predominately uses steady eruption models for interpretation and prediction of activity. My results suggest that unsteady models are necessary to accurately interpret behavior and assess threats from unsteady, short-lived eruptions. / Dissertation/Thesis / Ph.D. Geological Sciences 2012

<div>Pressure-gain combustion in wave rotors offer the opportunity for substantial improvement in gas turbine efficiency and power, while controlling emissions with fuel flexibility, if provided rapid and reliable ignition of lean mixtures. In addition, tightening emission regulations and increasing availability of gas fuels for internal-combustion engines require more reliable ignition for ultra-lean operation to avoid high peak combustion temperature. Turbulent jet ignition (TJI) is able to address the ignition challenges of lean premixed combustion. Especially, the turbulent hot jet results in faster ignition penetration for wave rotor pressure-gain combustors that have high-frequency operation and fast-burn requirements. Controllability of TJI needs better understanding of the chemistry and fluid mechanics in the jet mixing region, particularly the estimation of ignition delay time and identifying the location of the ignition onset. </div><div>In the present work, numerical and analytical methods are employed to develop models capable of estimating the ignition characteristics that the turbulent hot jet exhibits as it is issued to a cold stoichiometric CH4-H2-Air mixture with varied fuel reactivity blends. Numerical models of the starting turbulent jet are developed by Reynolds-averaged and large-eddy simulation of Navier-Stokes and scalar transport equations in a high-resolution computational domain, with major focus on ignition of high-reactivity fuel blends in the jet near-field due to computational resource limitations. The chemical reactions are modeled using detailed chemistry by well-stirred and partially stirred reactor approaches. Numerical models describe the temporal evolution of jet mixture fraction, scalar dissipation rate, flow strain rate, and thermochemical quantities of the flow.</div><div>For faster estimation of ignition characteristics, analytical methods are developed to explicitly solve governing equations for the transient evolution of the near field and the leading vortex of the starting hot jet. First, the transient radial evolution of the turbulent shear-layer of a round transient jet is analytically investigated in the near-field of the nozzle, where the momentum potential core exists. The methods approximate the mixing and chemical processes in the jet shear and mixing layer. The momentum equation is integrated analytically, with a mixing-length turbulence model to represent the variation of effective viscosity due to the velocity gradients. The analytic predictions of the velocity field and mass entrainment rate of the jet are compared with numerical predictions and experimental findings. In addition, the transport equation of conserved scalars in the jet near-field is solved analytically for the history of the jet mixture fraction. This analytic solution for temperature and species is used, together with available models for instantaneous chemical induction time, to create an analytic ignition model that provides the time and radial location of the ignition onset.</div><div>Lastly, the ignition mechanism within the vortex ring, which leads the starting turbulent jet, is modeled using prior understanding about the mixing characteristics of the vortex. This mechanism is more relevant to low-reactivity fuel blends. Due to the presence of strong mixing at the large-scale, the vortex ring is treated as a homogeneous batch-reactor, which contains certain levels of the jet mixture fraction. This assumption provides the initial composition and temperature of the reactor in which ignition ensues. </div><div>This article-dissertation is developed as a collection of 4 articles published in peer-reviewed journals, one submitted article, and additional unpublished work. The study is laid out in 6 chapters with the following contributions:</div><div>Chapter 1: This chapter numerically investigates the three-dimensional behavior of a transient hot jet as modeled using the Reynolds-averaged turbulence flow. The study aims at providing an insight towards the role of mixing in the ignition progress and how the operating conditions such as fuel mixture and pre-chamber pressure ratio can influence the ignition success. An ignition prediction criterion is developed in this chapter, which helps to predict the ignition success under a broad range of operating conditions.</div><div>Chapter 2: In this chapter, the large-eddy simulation (LES) of hot jet ignition is reported in conjunction with detailed kinetics mechanism and adaptive-mesh refinement. The correlation between local values of mixture fraction gradient and ignition is discussed. Furthermore, the role of methane-hydrogen ratio on the heat release pattern is studied for two specific mixtures.</div><div>Chapter 3: The LES of CH4-H2-Air ignition is extended in this chapter to account for multivariable evaluation of ignition. Joint probability assessment of ignition explains the role of important scalars on the formation and growth of ignition. Also, the effect of CH4-H2 ratio on the spatial distribution of ignition is assessed and discussed.</div><div>Chapter 4: In this chapter, the rate of mass entrainment into the jet in the near-field region is studied. Characterization of the mass entrainment illuminates the understanding of mixing behavior of the starting turbulent jets. Through an exact solution of the momentum equation, this chapter includes a model of the diffusive transport in a round transient jet at high Reynolds numbers.</div><div>Chapter 5: This chapter proposes a method to evaluate the mass/heat exchange between a transient-turbulent jet and a quiescent environment. To analyze the transport phenomena in the jet near-field, the transient diffusion equation in cylindrical coordinates is explicitly solved and its solution is compared with the empirical findings. The transport solution then enables an ignition model to describe the spatiotemporal characteristics of ignition in the near-field.</div><div>Chapter 6: The development of ignition within the vortex ring of the transient jet is investigated in this chapter. The initiation, growth, and departure of the vortex ring are studied using the available empirical correlations and the LES. Using a perfectly-stirred, zero-dimensional representation of the vortex, chemical kinetic calculations provide estimates of ignition delay for various fuel mixtures.</div><div><br></div>

The road transport sector, dominated by internal combustion engines, accounts for as high as 23% of annual carbon emissions and is considered the major area where urgent carbon reduction strategies are required. Natural gas is considered one of the intermediate fuels to reduce carbon emissions before net carbon neutral solutions can be achieved. Methane (CH4), a major constituent of natural gas, has the highest hydrogen-to-carbon ratio among the naturally occurring hydrocarbons, and the CO2 emission from natural gas combustion is around 25% less than diesel combustion. Lean combustion shows promises for improved engine efficiency, thereby reducing carbon emissions for a given required power output. However, igniting lean natural gas mixtures requires high ignition energy, beyond the capability of spark ig nition. The pre-chamber combustion (PCC) concept can provide the required ignition energy with relatively simple components. While most pre-chamber designs found in the literature are bulky and require extensive cylinder head modifications or complete engine redesign, the narrow-throat pre-chamber design can readily fit the diesel injector pockets of most heavy-duty engines without the need for substantial hardware modifications. The unique pre-chamber design is significantly different from the contemporary pre-chamber geometries, and its engine combustion phenomena and operating characteristics are largely unknown. This thesis work investigates the effect of important pre-chamber dimensions, such as the volume, nozzle hole diameter, and throat diameter, on the engine operating characteristics and emission trends. The experiments focus on the lean operation with excess air ratios (λ) exceeding 1.6, which can be achieved by auxiliary fuel injection into the pre-chamber. The air-fuel mixture formation process inside the pre-chamber is also investigated by employing 1-D and 3-D CFD simulations, where the engine experiments provided the boundary conditions. From the simulation results, a correlation between the injected and the trapped fuel in the pre-chamber is proposed by theoretical scavenging models to estimate the air-fuel ratio in the pre-chamber with high accuracy. Although the studies largely rely on thermodynamic engine experiments, the 1-D engine simulation implements the engine studies in estimating the mixture composition and heat transfer losses from the engine.

No description available.

No description available.

Estimates of bank erosion typically require field measurements to determine the soil erodibility since soil characteristics are highly variable between sites, especially for cohesive soils. The submerged jet test device is an in situ method of determining the critical shear stress and soil erodibility of cohesive soils. A constant velocity jet, applied perpendicular to the soil surface, creates a scour hole which is measured at discrete time intervals. While the results of these tests are able to provide values of critical shear stress and soil erodibility, the results are often highly variable and do not consider certain aspects of scour phenomena found in cohesive soils. Jet test measurements taken on the lower Roanoke River showed that the results varied for samples from similar sites and bulk failures of large areas of soil were common on the clay banks. Computational Fluid Dynamics (CFD) can be used to determine the effect of scour hole shape changes on the applied shear stress. Previous calculation methods assumed that the depth of the scour hole was the only parameter that affected the applied shear stress. The analysis of the CFD models showed that depth did heavily influence the maximum shear stress applied to the soil boundary. However, the scour hole shape had an impact on the flow conditions near the jet centerline and within the scour hole. Wide, shallow holes yielded results that were similar to the flat plate, therefore it is recommended that field studies only use jet test results from wide, shallow holes to determine the coefficient of erodibility and the critical shear stress of cohesive soils. / Master of Science

The objective of this thesis is two-fold: to analyze the directional calibration of a 3-hole probe in steady flow and to develop a method for the interpretation of measurements recorded with a novel, fast-response Pitot-type probe in unsteady, turbulent flow. Calibration data for the 3-hole probe’s two side ports was taken in the steady, non-turbulent region of a free jet and was evaluated for symmetry. In addition, data that was recorded using one side port in two independent calibration runs was compared to study repeatability. Misalignment was found between the nominally symmetric data sets, which may be the result of geometric probe tip defects or a misalignment of the side ports within -2 to -10 degrees. This misalignment suggested that the two probe ports must both be calibrated. The two data sets compared for repeatability were almost indistinguishable, suggesting that probe alignment was very repeatable over multiple calibration runs. This result implied that only one calibration run may be necessary for a single probe as well as for multiple probes having nearly identical tip geometry. These methods and findings from the 3-hole probe calibration provide useful processes and considerations for the calibration of directionally sensing pressure probes. Regarding the fast-response Pitot-type probe, measurements were conducted using the same free jet as was used with the 3-hole probe. The fast-response probe, which contains a Kulite sensor fitted in the sensing orifice of a Pitot tube, was positioned at incremental centerline locations in the unsteady, turbulent region of the jet flow, and mean and dynamic pressure data were recorded. Measurements were also taken at incremental centerline locations with a standard Pitot tube and a constant temperature hot wire anemometer. The Kulite mean pressure data and standard Pitot tube data were compared directly and agreed well. The hot wire data and a relevant turbulence model was used to generate mean pressure predictions, which correlated reasonably with a slight offset from the Kulite sensor and Pitot probe mean data. Next, the dynamic pressure data from the Kulite sensor was compared with predictions generated by the hotwire data, literature static pressure fluctuation data, and a second relevant turbulence model. In the centerline region where turbulence quantities begin to stabilize, the Kulite sensor data and predictions agreed reasonably well, within 7%. Thus, while not delivering ideal results, the turbulence models used provide a plausible method for the interpretation of the fast-response Pitot-type probe pressure measurements.

La présente étude porte sur la simulation numérique de la croissance des cristaux de glace dans le sillage proche d'une aile rectangulaire munie de deux injecteurs qui modélisent les deux moteurs. Dans cette configuration, les phénomènes microphysiques interviennent lors de l'interaction du jet, issu du moteur, et du tourbillon marginal qui se développe à chaque bout d'aile. Cet écoulement, très turbulent, perturbe fortement l'air environnant. Les jets diffusent dans l'atmosphère et s'enroulent autour des deux tourbillons de bout d'aile. Ces jets contiennent de la vapeur d'eau, des suies, des gaz mais également des aérosols et particules chargées. Le modèle microphysique utilisé dans cette étude repose sur l'hypothèse que la vapeur d'eau condense uniquement sur les particules de suie. Les simulations numériques sont effectuées à l’aide du code CEDRE développé à l’ONERA. Les méthodes numériques sont basées sur une approche volume finie pour des maillages non structurés généralisés. La résolution des équations de Navier stokes compressibles pour des fluides multi-espèces se fait selon une approche de type RANS et seul le champ stationnaire, jusqu'à huit envergures en aval de la maquette, est calculé. La turbulence de l'écoulement est modélisée au moyen du système de fermeture à deux équations k-l . Cette approche permet d'obtenir une description spatiale plus réaliste de l'interaction entre le jet et le tourbillon marginal. Le champ aérodynamique du sillage est ainsi comparé aux données expérimentales existantes. Le jet est correctement enroulé autour du tourbillon à huit envergures, et la dilution du panache est bien décrite par les simulations. Le modèle microphysique est ensuite couplé au modèle aérodynamique. Une première simulation porte sur les phénomènes microphysiques intervenant dans le sillage de la maquette dans des conditions particulières, représentatives d'un avion commercial en vol de croisière dans une atmosphère saturée par rapport à la glace. L'influence de la taille initiale des particules de suies émises par les moteurs ainsi que l'humidité relative de l'atmosphère, sur les propriétés de la traînée de condensation, sont ensuite étudiées et discutées. Ce travail, de part la stratégie de calcul mise en place et notamment l'utilisation de maillages non structurés généralisés, permettra d'appréhender le rôle de certains paramètres clés liés à l'avion comme la géométrie des ailes ou bien encore la position des nacelles sur les propriétés microphysiques de la traînée de condensation. / Numerical simulations of ice particles growth, in the near wake of a rectangular wing with two injectors, are presented in this study. In this configuration, microphysical processes occur during the interaction between the engines jets and the marginal vortices developping at each wing tip. This strong turbulent flow disturbs highly the environmental flow. The jets diffuse in the atmosphere and are wrapped around the two wing tip vortices. They contain water vapour, soots, gas, aerosols and charged particles as well. One of the hypothesis of the microphysical modeling, used in this study, is that water vapour condenses on soot particles only. Numerical simulations are performed with the code CEDRE developed at ONERA. The numerical methods are based on a cell-centered finite volume approach for general unstructured grids. A Navier-Stokes solver for turbulent, compressible and multi-species flows with a RANS approach, based on the k-l turbulence model, is used. Only stationary states of the flow, until eight spans downstream the setup, are computed. This approach enables to get a better spatial description of the interaction between the jet and the marginal vortex. The numerical flow field is then compared to existing experimental data. The jet is correctly wrapped around the wing tip vortex at eight spans and the dilution of the effluents is well described by our simulations. The microphysical model is then coupled to the aerodynamics. The microphysical processes occurring in the wake of the setup in specific conditions, representative of a cruising civil aircraft in an ice-saturated atmosphere, are firstly simulated. The role, on the contrail’s properties, of soot particles initial size and of the atmospheric humidity is studied and discussed. This work, through its computational strategy, with the use of unstructured grids, will enable to understand the potential role of some key parameters such as the wings geometry or the engines position on the contrail properties.