<div>The growing need to improve the power density of compact thermal systems necessitates developing new techniques to modulate the convective heat transfer efficiently. In the present research, acoustic streaming is evaluated as a potential technology to achieve this objective. Numerical simulations using the linearized and fully non-linear Navier-Stokes equations are employed to characterize the physics underlying this process. The linearized Navier-Stokes equations accurately replicate the low-frequency flow unsteadiness, which is used to find the optimal control parameters. Local and global stability analysis tools were developed to identify the modes with a global and positive heat transfer effect.</div><div><br></div><div>High-fidelity numerical simulations are performed to evaluate the effect of the excitation at selected frequencies, directed by the linear stability analysis, on the heat and momentum transport in the flow. Results indicate that, under favorable conditions, superimposing an acoustic wave, traveling along with the flow, can <i>resonate</i> within the domain and lead to a significant heat transfer enhancement with minimal skin friction losses. Two main flow configurations are considered; at the fixed Reynolds number Re_b=3000, in the supersonic case, 10.1% heat transfer enhancement is achieved by an 8.4% skin friction increase; however, in the subsonic case, 10% enhancement in heat transfer only caused a 5.3% increase to the skin friction. The deviation between these two quantities suggests a violation of the Reynolds analogy. This study is extended to include a larger Reynolds number, namely Re_b=6000 at M_b=0.75 and a similar response is observed. The effect of excitation amplitude and frequency on the resonance, limit-cycle oscillations, heat transfer, and skin friction are also investigated here.</div><div><br></div><div>Applying acoustic waves normal to the flow in the spanwise direction disrupts the near-wall turbulent structures that are primarily responsible for heat and momentum transport near the solid boundary. Direct numerical simulations were employed to investigate this technique in a supersonic channel flow at M_b=1.5 and Re_b=3000. The external excitation is applied through a periodic body force in the spanwise direction, mimicking loudspeakers placed on both walls that are operating with a 180^o phase shift. By keeping the product of forcing amplitude A_f and pulsation period (<i>T</i>) constant, spanwise velocity perturbations are generated with a similar amplitude at different frequencies. Under this condition, spanwise pulsations at <i>T</i>=20 and <i>T</i>=10 show up to 8% reduction in Nusselt number as well as the skin friction coefficient. Excitation at higher or lower frequencies fails to achieve such high level of modulations in heat and momentum transport processes near the walls.<br> <br>In configurations involving a spatially-developing boundary layer, a computational setup that includes laminar, transitional, and turbulent regions inside the domain is considered and the impact of acoustic excitation on this flow configuration has been characterized. Large-eddy simulations with dynamic Smagorinsky sub-grid scale modeling has been implemented, due to the excessive computational cost of DNS calculations at high-Reynolds numbers. The optimal excitation frequency that resembles the mode chosen for the fully-developed case has been identified via global stability analysis. Fully non-linear simulations of the spatially-developing boundary layer subjected to the excitation at this frequency reveal an interaction between the <i>pulsations</i> and the perturbations originated from the tripping which creates a re-laminarization zone traveling downstream. Such technique can locally enhance or reduce the heat transfer along the walls.<br></div>

L’urbanisation croissante fait émerger des enjeux sociétaux et environnementaux relatifs à la pollution atmosphérique et au microclimat urbain. La compréhension des phénomènes physiques de transport de quantité de mouvement, de chaleur et de masse entre la canopée urbaine et la couche limite atmosphérique est primordiale pour évaluer et anticiper les impacts négatifs de l’urbanisation. Les processus turbulents spécifiques à la couche limite urbaine sont étudiés par une approche de simulation des grandes échelles, dans une configuration urbaine représentée par un arrangement de cubes en quinconce. Le modèle de sous-maille de type Smagorinsky dynamique est implémenté pour mieux prendre en compte l’hétérogénéité de l’écoulement et les retours d’énergie des petites vers les grandes structures. Le nombre de Reynolds basé sur la hauteur du domaine et la vitesse de l’écoulement libre est de 50000. L’écoulement est résolu dans les sous-couches visqueuses et le maillage est raffiné dans la canopée. Le domaine est composé de 28 millions de cellules. Les résultats sont comparés à la littérature et aux données récentes obtenues dans la soufflerie du LHEEA. Chaque contribution au bilan d’énergie cinétique turbulente est calculée directement en tout point. Cette information, rare dans la littérature, permet d’étudier les processus dans la sous couche rugueuse. Grâce à ces résultats 3D, l’organisation complexe de l’écoulement moyen (recirculations, vorticité, points singuliers) est analysée en relation avec la production de turbulence. Enfin, une simulation où les obstacles sont remplacés par une force de traînée équivalente est réalisée à des fins d’évaluation de cette approche. / The rapid development of urbanization raises social and environmental challenges related to air pollution and urban climate. Understanding the physical processes of momentum, heat, and mass exchanges between the urban canopy and the atmospheric boundary-layer is a key to assess,predict and prevent negative impacts of urbanization. The turbulent processes occurring in the urban boundary-layer are investigated using computational fluid dynamics (CFD). The unsteady flow over an urban-like canopy modelled by a staggered arrangement of cubes is simulated using large eddy simulation (LES). Considering the highspatial and temporal in homogeneity of the flow, a dynamic Smagorinsky subgrid-scale model is implemented in the code to allow energyback scatter from small to large scales. The Reynolds number based on the domain height and free-stream velocity is 50000. The near-wall viscous sub-layers are resolved and the grid is refined in the canopy resulting in about 28 million grid cells. LES results are assessed by comparison with literature and data recently acquired in the wind tunnel of the LHEEA. The turbulent kinetic energy budget in which all contributions are independently computed is investigated. These rarely available data are used to analyse the turbulent processes in the urban canopy. By taking advantage of the three-dimensionality of the simulated flow, the complex 3D time-averaged organization of the flow (recirculation, vorticesor singular points) is analyzed in relation with production of turbulence. Finally a drag approach where obstacles are replaced by an equivalent drag force is implemented in the same domain and results are compared to obstacle-resolved data.

La nature conflictuelle des contraintes de performances, d'opérabilité et de respect des normes environnementales conduit les motoristes à optimiser finement la géométrie du brûleur afin d'identifier le meilleur design.La Simulation aux Grande Echelles (SGE) est aujourd'hui un outil performant et est déployé de manière courante dans les Bureaux d'Etudes pour la prédiction des propriétés macroscopiques de l'écoulement.Toutefois, de nombreux phénomènes influencés par les effets de chimie complexe, tels que la stabilisation, l'extinction de flamme et la formation des polluants, reste un problème crucial.En effet, la description des effets de chimie complexe nécessite l'utilisation de modèles cinétiques détaillés imposant des coûts de calculs prohibitifs, des problèmes de raideurs numérique et des difficultés de couplage avec les échelles non résolues turbulentes.Afin d'inclure une description des processus chimiques, dans les simulations numériques de chambres de combustion réelles, des modèles réduits doivent être proposés.Dans cette thèse, une méthode originale, appelée chimie virtuelle optimisée, est développée.Cette stratégie a pour objectif la description de la structure chimique de la flamme et la formation des polluants dans des configurations de flamme représentatives.Les schémas cinétiques virtuels optimisés, composés de réactions virtuelles et d'espèces virtuelles, sont construits par optimisation des paramètres réactionnels et des propriétés thermochimiques des espèces virtuelles afin de capturer les propriétés de flamme d'intérêt. / The conflicting nature of performance, operability and environmental constraints leads engine manufacturers to perform a fine optimization of the burner geometry to find the best design compromise.Large Eddy Simulation (LES) is an attractive tool to achieve this challenging task, and is routinely used in design office to capture macroscopic flow features.However, the prediction of phenomena influenced by complex kinetic effects, such as flame stabilization, extinction and pollutant formation, is still a crucial issue.Indeed, the comprehensive description of combustion chemistry effects requires the use of detailed models imposing prohibitive computational costs, numerical stiffness and difficulties related to model the coupling with unresolved turbulent scales.Reduced-cost chemistry description strategies must then be proposed to account for kinetic effects in LES of real combustion chambers.In this thesis an original modeling approach, called virtual optimized chemistry, is developed.This strategy aims at describing the chemical flame structure and pollutant formation in relevant flame configurations, at a low computational cost.Virtual optimized kinetic schemes, composed by virtual reactions and virtual species, are built through optimization of both kinetic rate parameters and virtual species thermo-chemical properties so as to capture reference target flame quantity.

Previous simulations have shown that wind farms have an impact on the near-surface atmospheric boundary layer (ABL) as turbulent wakes generated by the turbines enhance vertical mixing of momentum, heat and moisture. These changes alter downstream atmospheric properties. With the exception of a few observational data sets that focus on the impact to near-surface temperature within wind farms, little to no observational evidence exists with respect to vertical mixing. These few experimental studies also lack high spatial resolution due to their use of a limited number of meteorological sensors or remote sensing techniques. This study utilizes an instrumented small unmanned aerial system (sUAS) to gather high resolution in-situ field measurements from two state-of-the-art Midwest wind farms in order to differentially map downstream changes to relative humidity. These measurements are complemented by numerical experiments conducted using large eddy simulation (LES). Observations and numerical predictions are in good general agreement around a single wind turbine and show that downstream relative humidity is altered in the vertical, lateral, and downstream directions. A suite of LES is then performed to determine the effect of a turbine array on the relative humidity distribution in compounding wakes. In stable and neutral conditions, and in the presence of a positive relative humidity lapse rate, it is found that the humidity decreases below the turbine hub height and increases above the hub height. As the array is transitioned, the magnitude of change increases, differentially grows on the left-hand and right-hand side of the wake, and move slightly upward with downstream distance. In unstable conditions, the magnitude of near-surface decrease in relative humidity is a full order of magnitude smaller than that observed in a stable atmospheric regime.

Couette, pipe, channel, and zero-pressure gradient (ZPG) turbulent boundary layer (TBL) flows have classically been considered as canonical wall-bounded turbulent flows since their near-wall behavior is generally considered to be universal, i.e. invariant of the flow case and the Reynolds number. Nevertheless, the idea that large-scale motions, being dominant in regions further away from the wall, might interact with and influence small-scale fluctuations close to the wall has not been disregarded. This view was mainly motivated due to the observed failure of collapse of the Reynolds normal stresses in viscous scaling. While this top-down influence has been studied extensively over the last decade, the idea of a bottom-up influence (backward energy transfer) is less examined. One exception was the recent experimental work on a Couette flow by Kawata, T. &amp; Alfredsson, P. H. (Phys. Rev. Lett. 120, 244501, 2018). In the present work, a spectral representation of the Reynolds Stress transport equation is used to perform a scale-by-scale analysis of the terms in the equation. Two flow cases were studied: first, a Direct Numerical Simulation (DNS) of a Couette flow at a similar Reynolds number as Kawata and Alfredsson. The Reynolds number was ReT = 120, viscosity v. Second, a Large Eddy Simulation (LES) of a ZPG TBL at ReT = 730, 1270, and 2400. For both cases the classic interscale transport or turbulent kinetic energy was observed. However, also an inverse interscale transport of Reynolds shear stress was observed for both cases.

<p>This thesis is concerned with particle subgrid scale (SGS) modeling in large-eddy simulation (LES) of particle-laden turbulence. Although most particle-laden LES studies have neglected the effect of the subgrid scales on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM), and the stochastic models of Fukagata et al. (2004), Shotorban and Mashayek (2006) and Berrouk et al. (2007) are analyzed. The particle SGS models are assessed by conducting both a priori and a posteriori tests of a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re=74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient energy, but show a large range of results dependent on Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers. Due to the random component, the stochastic models are unable to predict preferential concentration.</p> / Master of Applied Science (MASc)

The unsteady aerodynamic and aeroelastic analysis of flapping flight under various kinematics and flow parameters is presented in this dissertation. The main motivation for this study arises from the challenges facing the development of micro air vehicles. Micro air vehicles by requirement are compact with dimensions less than 15-20 cm and flight speeds of around 10-15 m/s. These vehicles operate in low Reynolds number range of 10,000 to 100,000. At these low Reynolds numbers, the aerodynamic efficiency of conventional fixed airfoils significantly deteriorates. On the other hand, flapping flight employed by birds and insects whose flight regime coincides with that of micro air vehicles offers a viable alternate solution. For the analysis of flapping flight, a boundary fitted moving grid algorithm is implemented in a flow solver, GenIDLEST. The dynamic movement of the grid is achieved using a combination of spring analogy and trans-finite interpolation on displacements. The additional conservation equation of space required for moving grid is satisfied. The solver is validated with well known flow problems such as forced oscillation of a cylinder, a heaving airfoil, a moving indentation channel, and a hovering fruitfly. The performance of flapping flight is analyzed using Large Eddy Simulation (LES) for a wide range of Reynolds numbers and under various kinematic parameters. A spiral Leading Edge Vortex (LEV) forms during the downstroke due to the high angle of attack, which results in high force production. A strong spanwise flow of the order of the flapping velocity is observed along the core of the LEV. In addition, the formation of a negative spanwise flow is observed due to the tip vortex, which slows down the removal of vorticity from the LEV. This leads to the instability of the LEV at around mid-downstroke. Analysis with different rotation kinematics shows that a continuous rotation results in better propulsive efficiency as it generates thrust during the entire flapping cycle. Analysis with different angles of attack shows that a moderate angle of attack which results in complete shedding of the LEV offers high propulsive efficiency. The analysis of flapping flight at Reynolds numbers ranging from 100 to 100,000 shows that higher lift and thrust values are obtained for Re?100. The critical reasons are that at higher Reynolds numbers, the LEV is closer to the surface and as it sheds and convects it covers most of the upper surface. However, the Reynolds number has no or little effect on the lift and thrust as identical values are obtained for Re=10,000 and 100,000. The analysis with different tip shapes shows that tip shapes do not have a significant effect on the performance. Introduction of stroke deviation to kinematics leads to drop in average lift as wing interacts with the LEV shed during the downstroke. A linear elastic membrane model with applied aerodynamic load is developed for aeroelastic analysis. Analysis with different wing stiffnesses shows that the membrane wing outperforms the rigid wing in terms of lift, thrust and propulsive efficiency. The main reason for the increase in force production is attributed to the gliding of the LEV along the camber, which results in a high pressure difference across the surface. In addition, a high stiffness along the spanwise direction and low stiffness along the chordwise direction results in a uniform camber and high lift and thrust production. / Ph. D.

Reduced-order models are frequently used in the simulation of complex flows to overcome the high computational cost of direct numerical simulations, especially for three-dimensional nonlinear problems. Proper orthogonal decomposition, as one of the most commonly used tools to generate reduced-order models, has been utilized in many engineering and scientific applications. Its original promise of computationally efficient, yet accurate approximation of coherent structures in high Reynolds number turbulent flows, however, still remains to be fulfilled. To balance the low computational cost required by reduced-order modeling and the complexity of the targeted flows, appropriate closure modeling strategies need to be employed. In this dissertation, we put forth two new closure models for the proper orthogonal decomposition reduced-order modeling of structurally dominated turbulent flows: the dynamic subgrid-scale model and the variational multiscale model. These models, which are considered state-of-the-art in large eddy simulation, are carefully derived and numerically investigated. Since modern closure models for turbulent flows generally have non-polynomial nonlinearities, their efficient numerical discretization within a proper orthogonal decomposition framework is challenging. This dissertation proposes a two-level method for an efficient and accurate numerical discretization of general nonlinear proper orthogonal decomposition closure models. This method computes the nonlinear terms of the reduced-order model on a coarse mesh. Compared with a brute force computational approach in which the nonlinear terms are evaluated on the fine mesh at each time step, the two-level method attains the same level of accuracy while dramatically reducing the computational cost. We numerically illustrate these improvements in the two-level method by using it in three settings: the one-dimensional Burgers equation with a small diffusion parameter, a two-dimensional flow past a cylinder at Reynolds number Re = 200, and a three-dimensional flow past a cylinder at Reynolds number Re = 1000. With the help of the two-level algorithm, the new nonlinear proper orthogonal decomposition closure models (i.e., the dynamic subgrid-scale model and the variational multiscale model), together with the mixing length and the Smagorinsky closure models, are tested in the numerical simulation of a three-dimensional turbulent flow past a cylinder at Re = 1000. Five criteria are used to judge the performance of the proper orthogonal decomposition reduced-order models: the kinetic energy spectrum, the mean velocity, the Reynolds stresses, the root mean square values of the velocity fluctuations, and the time evolution of the proper orthogonal decomposition basis coefficients. All the numerical results are benchmarked against a direct numerical simulation. Based on these numerical results, we conclude that the dynamic subgrid-scale and the variational multiscale models are the most accurate. We present a rigorous numerical analysis for the discretization of the new models. As a first step, we derive an error estimate for the time discretization of the Smagorinsky proper orthogonal decomposition reduced-order model for the Burgers equation with a small diffusion parameter. The theoretical analysis is numerically verified by two tests on problems displaying shock-like phenomena. We then present a thorough numerical analysis for the finite element discretization of the variational multiscale proper orthogonal decomposition reduced-order model for convection-dominated convection-diffusion-reaction equations. Numerical tests show the increased numerical accuracy over the standard reduced-order model and illustrate the theoretical convergence rates. We also discuss the use of the new reduced-order models in realistic applications such as airflow simulation in energy efficient building design and control problems as well as numerical simulation of large-scale ocean motions in climate modeling. Several research directions that we plan to pursue in the future are outlined. / Ph. D.

<p dir="ltr">Grid convergence and scaling analyses have not been done rigorously for practical large-eddy simulations (LES). The challenge arises from the fact that there are two grid-related length scales: grid size and LES filter width. It causes the numerical and model errors in LES to be inherently coupled, making the convergence of either error difficult to analyze. This study works to overcome the challenge by developing scaling laws that can be used to guide the convergence analysis of errors in LES. Three different convergence cases are considered, and their respective scaling laws are developed by varying the ratio between grid size and filter width. A pool fire is adopted as a test case for the convergence analysis of LES. Qualitative and quantitative assessments of the LES results are made first to ensure reliable numerical solutions. In the subsequent scaling analysis, it is found that the results are consistent with their respective scaling laws. The results provide strong support to the developed scaling laws. The work is significant as it proposes a rigorous way to guide convergence analysis of LES errors. In a world where LES already has a wide range of applicability and is still becoming more prominent, it is imperative to have a thorough understanding of how it works including its convergence and scaling laws with respect to the change of grid size and filter width.</p>

In this thesis, the ribbed ducts of the internal cooling passage in turbine blading are investigated to demonstrate the effects of high speed rotation. Rotation coupled with high temperature operating conditions alters the mean flow, turbulence, and heat transfer augmentation due to Coriolis and centrifugal buoyancy forces that arises from density stratification in the domain. Gas turbine engines operate in particle laden environments (sand, volcanic ash), and particulate matter ingested by the engine can make their way into the blade internal cooling passages over thousands of operating hours. These particulates can deposit on the walls of these cooling passages and degrade performance of the turbine blade. Large-Eddy Simulations (LES) with temperature dependent properties is used for turbulent flow and heat transfer in the ribbed cooling passages and Lagrangian tracking is used to calculate the particle trajectories together with a wall deposition model. The conditions used are Re=100,000, Rotation number, Ro = 0.0 and 0.2, and centrifugal Buoyancy parameters of Bo=0, 0.5, and 1.0. First, the independent effects of Coriolis and centrifugal buoyancy forces are investigated, with a focus on the additional augmentation obtained in heat transfer with the addition of centrifugal buoyancy. Coriolis forces are known to augment heat transfer at the trailing wall and attenuate the same at the leading wall. Phenomenological arguments stated that centrifugal buoyancy augments the effects of Coriolis forces in outward flow in the first pass while opposing the effect of Coriolis forces during inward flow in the second pass. In this study, it was found that in the first pass, centrifugal buoyancy had a greater effect in augmenting heat transfer at the trailing wall than in attenuating heat transfer at the leading wall. On the contrary, it aided heat transfer in the second half of the first pass at the leading wall by energizing the flow near the wall. Also, contrary to phenomenological arguments, inclusion of centrifugal buoyancy augmented heat transfer over Coriolis forces alone on both the leading and trailing walls of the second pass. Sand ingestion is then investigated, by injecting 200,000 particles in the size range of 0.5-175μm with 65% of the particles below 10 μm. Three duct wall temperatures are investigated, 950, 1000 and 1050 °C with an inlet temperature of flow and particles at 527 °C . The impingement, deposition levels, and impact characteristics are recorded as the particles move through the domain. It was found that the Coriolis force greatly increases deposition. This was made prevalent in the first pass, as 84% of the deposits in the domain occurred in the first pass for the rotating case, whereas only 27% of deposits occurred in the first pass for the stationary case with the majority of deposits occurring in the bend region. This was due to an increased interaction with the trailing wall in the rotating case whereas particles in the stationary case were allowed to remain in the mean flow and gain momentum, making rebounding from a wall during collision more likely than deposition. In contrast, the variation of wall temperatures caused little to no change in deposition levels. This was concluded to be a result of the high Reynolds number used in the flow. At high Reynolds numbers, the particles have a short residence times in the internal cooling circuit not allowing the flow and particles to heat up to the wall temperature. Overall, 87% of the injected particles deposited in the rotating duct whereas 58% deposited in the stationary duct. / Master of Science