The complex multiscale physics of nanoparticle laden functional droplets in a reacting environment is of fundamental and applied significance for a wide variety of applications ranging from thermal sprays to pharmaceutics to modern day combustors using new brands of bio-fuels. Understanding the combustion characteristics of these novel fuels (laden with energetic nanoparticle NP) is pivotal for lowering ignition delay, reducing pollutant emissions and increasing the combustion efficiency in next generation combustors. On the way to understanding the complex dynamics of sprays is to first study the behaviour of an isolated droplet. A single droplet represents a sub-grid unit of spray. In vaporizing functional droplets under high heat flux conditions, the bubble formation inside the droplet represents an unstable system. This may be either through homogenous nucleation at the superheat limit or by dispersed nanoparticle acting as heterogeneous nucleation sites. First it is shown that such self-induced boiling in burning functional pendant droplets can induce severe volumetric shape oscillations in the droplet. Internal pressure build-up due to ebullition activity force ejects bubbles from the droplet domain causing undulations on the droplet surface and oscillations in bulk thereby leading to secondary break-up of the primary droplet. Through experiments, it is established that the degree of droplet deformation depends on the frequency and intensity of these bubble expulsion events. However, in a distinct regime of single isolated bubble growing inside the droplet, pre-ejection transient time is identified by Darrieus-Landau (DL) instability at the evaporative bubble-droplet interface. In this regime the bubble-droplet system behaves as a synchronized driver-driven system with bulk bubble-shape oscillations being imposed on the droplet. However, the agglomeration of suspended anaphase additives modulates the flow structures within the droplet and also influences the bubble inception and growth leading to distinct atomization characteristics. Secondly, the secondary atomization characteristics of burning bi-component (ethanol-water) droplets containing titania nanoparticle (NPs) at both dilute (0.5% and 1% by weight) and dense particle loading rates (PLR: 5% and 7.5 wt. %) are studied experimentally at atmospheric pressure under normal gravity. It is observed that both types of nanofuel droplets undergo distinct modes of secondary break-up that are primarily responsible for transporting particles from the droplet domain to the flame zone. For dilute nanosuspensions, disruptive response is characterized by low intensity atomization modes that cause small-scale localized flame distortion. In contrast, the disruption behavior at dense concentrations is governed by high intensity bubble ejections which result in severe disruption of the flame envelope. The atomization events occur locally at the droplet surface while their cumulative effect is observed globally at the droplet scale. Apart from this, a feedback coupling between two key interacting mechanisms, namely, atomization frequency and particle agglomeration also influence the droplet deformation characteristics by regulating the effective mass fraction of NPs within the droplet. Thus, third part of the study elucidates how the initial NP concentration modulates the relative dominance of these two mechanisms thereby leading to a master-slave configuration. Secondary atomization of novel nanofuels is a crucial process since it enables an effective transport of dispersed NPs to the flame (a pre-requisite condition for NPs to burn). Contrarily, NP agglomeration at the droplet surface leads to shell formation thereby retaining NPs inside the droplet. In particular, it is shown that at dense concentrations shell formation (master process) dominates over secondary atomization (slave) while at dilute particle loading it is the high frequency bubble ejections (master) that disrupt shell formation (slave) through its rupture and continuous out flux of NPs. These results in distinct combustion residues at dilute and dense concentrations, thus, providing a method of manufacturing flame synthesized microstructures with distinct morphologies. Next, it is shown that by using external stimuli (preferential acoustic excitation) the secondary atomization of the droplet can be suppressed i.e. the external flame-acoustic interaction with bubbles inside the droplet results in controlled droplet deformation. Particularly, by exciting the droplet flame in a critical, responsive frequency range i.e. 80 Hz ≤ fP ≤ 120 Hz, the droplet deformation cycle is altered through suppression of self-excited instabilities and intensity/frequency of bubble ejection events. The acoustic tuning also enables the control of bubble dynamics, bulk droplet-shape distortion and final precipitate morphology even in burning nanoparticle laden droplets. Droplets in a non-reacting environment (heated radioactively) are also subject to instabilities. One such instability observed in drying colloidal droplets is the buckling of thin viscoelastic shell formed through consolidation of NPs. In the final part of the thesis, buckling instability driven morphology transition (sphere to ring structure) in an acoustically levitated heated nanosilica dispersion droplet is elucidated using dynamic energy balance. Droplet deformation featuring formation of symmetric cavities is initiated by the capillary pressure that is two to three orders of magnitude greater than acoustic radiation pressure, thus indicating that the standing pressure field has no influence on the buckling front kinetics. With increase in heat flux, the growth rate of surface cavities and their post-buckled volume increases while the buckling time period reduces, thereby altering the buckling pathway and resulting in distinct precipitate structures. Thus, the cavity growth is primarily driven by evaporation. However, irrespective of the heating rate, volumetric droplet deformation exhibits linear time dependence and droplet vaporization is observed to deviate from the classical D2-law. Understanding such transients of buckling phenomenon in drying colloidal suspensions is pivotal for producing new functional microstructures with tenable morphology and is particularly critical for spray drying applications that produce powders through vaporization of colloidal droplets.

Ce travail porte sur la stabilisation de flammes prémélangées et swirlées de mélanges combustibles méthane/hydrogène/air avec différents taux de dilution d’azote et de dioxyde de carbone. Une tige centrale permet de stabiliser des flammes pour de faibles nombres de swirl. Le sommet de la flamme interagît éventuellement avec les parois de la chambre de combustion. L’objectif ces travaux est d’améliorer la connaissance des mécanismes qui gouvernent la stabilisation et la topologie de ces flammes. Ces travaux démontrent que le nombre de swirl, la composition du mélange combustible, la géométrie de la chambre de combustion ainsi que les conditions aux limites thermiques ont une grande influence sur la forme prise par la flamme. Le dispositif expérimental permet de modifier la forme et la taille de la chambre de combustion, le diamètre du tube d’injection et le nombre de swirl. Des conditions opératoires propices aux transitions de forme de flamme sont ensuite étudiées pour différentes configurations de brûleur. Une caractérisation expérimentale fouillée d’un point de fonctionnement est réalisée grâce à la Fluorescence Induite par Laser sur le radical Hydroxyle (OH-PLIF), la Vélocimétrie par Images de Particules (PIV) et la Phosphorescence Induite par Laser de phosphores sensibles à la température (LIP). Une base de donnée de l’écoulement et des conditions aux limites associées est obtenue sans et avec combustion. Les mécanismes qui contrôlent les transitions de formes de flamme sont ensuite analysés lorsque la flamme interagit avec les parois de la chambre de combustion. L’influence de la composition du mélange combustible, de la vitesse débitante et du nombre de swirl est caractérisée et il est démontré que la transition d’une flamme en V vers une flamme en M est déclenchée par un retour de flamme dans la couche limite le long d’une des parois latérales de la chambre de combustion. Les nombres sans dimension contrôlant ces transitions sont identifiés et un modèle de prévision de la forme de ces flammes est développé. La physique déterminant les transitions de forme de flammes est différente lorsque celles-ci n’interagissent pas avec les parois de la chambre de combustion. En utilisant le signal de chimiluminescence OH* et la OH-PLIF, il est montré que la teneur en hydrogène dans le combustible a une grande influence sur la forme de flamme. L’utilisation de la LIP et de thermocouples a également permis de montrer que les conditions aux limites thermiques jouent un rôle prépondérant sur la forme de flamme. Les effets combinés de l’étirement et des pertes thermiques sont examinés par l’utilisation conjointe de la PIV et de la OH-PLIF. Il est montré que les limites d’extinction de flammes pauvres prémélangées sont réduites par les pertes thermiques et que la transition d’une flamme en M vers une flamme en V est consécutive à l’extinction du front de flamme situé dans la couche de cisaillement externe du jet soumis à un étirement trop important. Ces expériences sont complétées par une analyse de la dynamique de ces flammes. Des modulations de la vitesse débitante à basse fréquence et à haute amplitude modifient la forme de flamme. La stabilisation de flammes CH4/H2/air diluées par du N2 ou du CO2 est finalement examinée. La zone de recirculation produite par la tige centrale permet d’alimenter la base de la flamme avec des gaz brûlés chauds et de stabiliser des flammes fortement diluées. Augmenter la fraction molaire de diluant dans le combustible réduit l’intensité de lumière émise par le radical OH*. Il est également montré que la composition du diluant a un impact sur le champ de température des gaz brûlés et des surfaces de la chambre de combustion. La dilution par du CO2 augmente les pertes thermiques par rayonnement des gaz brûlés. Cela réduit l’efficacité de la chambre de combustion équipée de quatre parois transparentes. [...] / This work deals with the stabilization of premixed turbulent swirling flames of methane/hydrogen/air combustible mixtures with different dilution rates of nitrogen and carbon dioxide. A central bluff body helps stabilizing the flames at low swirl numbers. The flame tip eventually impinges the combustor peripheral wall. The general objective is to gain understanding of the mechanisms governing the stabilization and the topology of these flames. It is found that the swirl number, the combustible mixture composition, the geometry of the combustor, and the thermal boundary conditions have a strong impact on the shape taken by these flames. The experimental setup used to characterize flames topologies is first described. Flames prone to topology bifurcations are selected and are studied for different arrangement of the combustor when the combustion chamber shape and size, the injection tube diameter, and swirl number are varied. One operating condition is fully characterized under non-reactive and reactive conditions using Planar Hydroxyl Laser Induced Fluorescence (OH-PLIF), Particle Imaging Velocimetry (PIV), and Laser Induced Phosphorescence of thermographic phosphors (LIP) to generate a detailed database of the flow and the corresponding boundary conditions. An analysis is then conducted to understand the mechanisms controlling shape bifurcations when the flame interacts with the combustor peripheral wall. Effects of the combustible mixture composition, the bulk flow velocity, and the swirl number are analyzed. It is shown that the transition from a V to an M flame is triggered by a flashback of the V flame tip in the boundary layer of the combustor peripheral wall. Dimensionless numbers controlling these transitions are identified and a simplified model is developed to help the prediction of the flame shapes. The physics of these shape bifurcations differs when the flame does not interact with the combustor wall. The large influence of the hydrogen enrichment in the fuel on the flame shape is analyzed using flame chemiluminescence and OH-PLIF. LIP and thermocouple measurements demonstrate that the thermal boundary conditions still have a strong impact on the flame topology. The combined effects of strain and heat losses are investigated using joint OH-PLIF and PIV experiments. It is shown that flammability limits of premixed flames are reduced due to heat losses and the transitions from M to V shaped flames is consecutive to localized extinctions of flame front elements located in the outer shear layer of the jet flow that are submitted to large strain rates. These experiments are completed by an analysis of the dynamics of methane/hydrogen/air flames. It is shown that low frequency and high amplitude velocity modulations generated by a loudspeaker alter the shape taken by these flames. The stabilization of methane/hydrogen/air flames diluted by nitrogen and carbon dioxide is finally examined. It was possible to stabilize swirled flames featuring important dilution rates due to the presence of the bluff body, installed on the axis of the injection tube. The recirculation zone behind this element supplies hot burnt gases to the flame anchoring point. Using OH* chemiluminescence imaging, it is shown than increasing the molar fraction of diluent in the fuel reduces the light emission from excited OH* radicals. The influence of dilution on the flame chemistry is emphasized with experiments conducted at a fixed thermal power and fixed adiabatic flame temperature. It is also demonstrated that the composition of the diluent has a strong influence on the temperature field of the burnt gases and of the combustor wall surfaces. Dilution with carbon dioxide increases radiative heat losses from the burnt gases in comparison to dilution with nitrogen. This penalizes the combustor efficiency equipped with four transparent quartz walls. [...]

Ce travail porte sur l'étude expérimentale et numérique de jets coaxiaux de rapports de vitesses débitantes intérieure/extérieure) inférieurs à l’unité, présentant une rotation dans le jet annulaire. Dans un premier temps, des visualisations par tomographie laser ont été réalisées dans les plans méridien et transversaux permettant une description tridimensionnelle de l'écoulement. Pour différents nombres de Reynolds, de rapport de vitesses et de nombre de Swirl, un inventaire des modes dominants a pu être établi. Les champs instantanés de vitesses ont ensuite été mesurés par Vélocimétrie par Imagerie de Particules (PIV). Les résultats de mesures de vitesses longitudinales et azimutales moyennées en temps sont présentés. Une comparaison avec les structures observées par tomographie est proposée. Une décomposition de Fourier a été effectuée permettant d'identifier les modes dominants ainsi que leur position dans la direction radiale. L'approche expérimentale a été suivie par une analyse de stabilité linéaire. Une attention particulière est portée sur l'état de base stationnaire reconstruit à partir des profils de vitesses mesurées par PIV à la sortie du jet. Étant donné que l'écoulement est non parallèle, une approche globale de la stabilité est utilisée. L'étude de la stabilité est basée sur la résolution numérique des équations de Navier-Stokes par des méthodes pseudo-spectrales. L'objectif de cette analyse est de retrouver la carte en mode de Fourier azimutaux observée expérimentalement. Nous nous sommes donc intéressés au taux de croissance le plus élevé de la perturbation pour chaque mode azimutal ainsi qu'à la nature convective/ absolue des modes. Pour finir, une com / This work concerns the experimental and numerical study of coaxial jets with outer to inner velocity ratio lower than unity, presenting a rotation in the annular jet. At first, flow visualizations by tomography laser were used in the meridian and transverse plans in order to provide a spatial description of the flow. For various values of the nondimensional parameters : numbers of Reynolds, outer to inner velocity ratio and Swirl number, an inventory of the dominant modes was be established. Instantaneous velocity fields were then measured by Particle lmaging Velocimetry (PIV). The results of longitudinal and azimuthal time-averaged Velocity fields are presented. A comparison with the structures observed by tomography is proposed. A Fourier decomposition was made allowing to identify the dominant modes as well as their position in the radial direction. Experimental investigation was followed by a linear stability analysis. Special attention is paid to the steady base-flow solution reconstructed from the velocity profiles measured by PIV at the end of the nozzle. Given that the is not parallel, a global approach was used. Study of the stability is based on the numerical solution of the incompressible Navier-Stokes equations with pseudo-spectral methods. The objective of this analysis is to the map of azimuthal Fourier modes observed experimentally. We were thus interested in the most amplified growth rate of the disturbance for every azimuthal mode as well as in the absolute/convective nature of the modes. To conclude, a comparison of the results obtained in both numerical and experimental approaches is proposed.

Global renewable energy requirements rapidly increase with the transition to a fossil-free society. As a result, intermittent energy resources, such as wind- and solar power, have become increasingly popular. However, their energy production varies over time, both in the short- and long term. Hydropower plants are therefore utilized as a regulating resource more frequently to maintain a balance between production and consumption on the electrical grid. This means that they must be operated away from the design point, also known as the best-efficiency-point (BEP), and often are operated at part-load (PL) with a lower power output. Moreover, some plants are expected to provide a spinning reserve, also referred to as speed-no-load (SNL), to respond rapidly to power shortages. During this operating condition, the turbine rotates without producing any power. During the above mentioned off-design operating conditions, the flow rate is restricted by the closure of the guide vanes. This changes the absolute velocity of the flow and increases the swirl, which is unfavorable. The flow field can be described as chaotic, with separated regions and recirculating fluid. Shear layer formation between stagnant- and rotating flow regions can be an origin for rotating flow structures. Examples are the rotating-vortex-rope (RVR) found during PL operation and the vortical flow structures in the vaneless space during SNL operation, which can cause the flow between the runner blades to stall, also referred to as rotating stall. The flow structures are associated with pressure pulsations throughout the turbine, which puts high stress on the runner and other critical parts and shortens the turbine's lifetime. Numerical models of hydraulic turbines are highly coveted to investigate the detrimental flow inside the hydraulic turbines' different sections at off-design operating conditions. They enable the detailed study of the flow and the origin of the instabilities. This knowledge eases the design and assessment of mitigation techniques that expand the turbines' operating range, ultimately enabling a wider implementation of intermittent energy resources on the electrical grid and a smoother transition to a fossil-free society. This thesis presents the numerical study of the Porjus U9 model, a scaled-down version of the 10 MW prototype Kaplan turbine located along the Luleå river in northern Sweden. The distributor contains 20 guide vanes, 18 stay vanes and the runner is 6-bladed. The numerical model is a geometrical representation of the model turbine located at Vattenfall Research and Development in Älvkarleby, Sweden. The commercial software ANSYS CFX 2020 R2 is used to perform the numerical simulations. Firstly, the draft tube cone section of the U9 model is numerically studied to investigate the sensitivity of a swirling flow to the GEKO (generalized kω) turbulence model. The GEKO model aims to consolidate different eddy viscosity turbulence models. Six free coefficients are changeable to tune the model to flow conditions and obtain results closer to an experimental reference without affecting the calibration of the turbulence model to basic flow test cases. Especially, the coefficients affecting wall-bounded flows are of interest. This study aims to analyze if the GEKO model can be used to obtain results closer to experimental measurements and better predict the swirling flow at PL operation compared to other eddy viscosity turbulence models. Results show that the near-wall- and separation coefficients predict a higher swirl and give results closer to experimentally obtained ones. Secondly, a simplified version of the U9 model is investigated at BEP and PL operating conditions and includes one distributor passage with periodic boundary conditions, the runner and the draft tube. The flow is assumed axisymmetric upstream of the runner, hence the single distributor passage. Previous studies of hydraulic turbines operating at PL show difficulties predicting the flow's tangential velocity component as it is often under predicted. Therefore, a parametric analysis is performed to investigate which parameters affect the prediction of the tangential velocity in the runner domain. Results show that the model predicts the flow relatively well at BEP but has problems at PL; the axial velocity is overpredicted while the tangential is underpredicted. Moreover, the torque is overpredicted. The root cause for the deviation is an underestimation of the head losses. Another contributing reason is that the runner extracts too much swirl from the flow, hence the low tangential velocity and the high torque. Sensitive parameters are the blade clearance, blade angle and mass flow. Finally, the full version of the U9 model is analyzed at SNL operation, including the spiral casing, full distributor, runner and draft tube. During this operating condition, the flow is not axisymmetric; vortical flow structures extend from the vaneless space to the draft tube and the flow stalls between the runner blades. A mitigation technique with independent control of each guide vane is presented and implemented in the model. The idea is to open some of the guidevanes to BEP angle while keeping the remaining ones closed. The aim is to reduce the swirl and prevent the vortical flow structures from developing. Results show that the flow structures are broken down upstream the runner and the rotating stall between the runner blades is reduced, which decreases the pressure- and velocity fluctuations. The flow down stream the runner remains mainly unchanged.