The linear stability of high speed-viscous plane Couette and Couette-Poiseuille flows are investigated numerically. The conservation equations along with Sutherland&amp / #65533 / s viscosity law are studied using a second-order finite difference scheme. The basic velocity and temperature distributions are perturbed by a small-amplitude normalmode disturbance. The small-amplitude disturbance equations are solved numerically using a global method using QZ algorithm to find all the eigenvalues at finite Reynolds numbers, and the incompressible limit of these equations is investigated for Couette-Poiseuille flow. It is found that the instabilities occur, although the corresponding growth rates are often small. Two families of wave modes, Mode I (odd modes) and Mode II (even modes), were found to be unstable at finite Reynolds numbers, where Mode II is the dominant instability among the unstable modes for plane Couette flow. The most unstable mode for plane Couette &amp / #65533 / Poiseuille flow is Mode 0, which is not a member of the even modes. Both even and odd modes are acoustic modes created by acoustic reflections between a will and a relative sonic line. The necessary condition for the existence of such acoustic wave modes is that there is a region of locally supersonic mean flow relative to the phase speed of the instability wave. The effects of viscosity and compressibility are also investigated and shown to have a stabilizing role in all cases studied. Couette-Poiseuille flow stability is investigated in case of a choked channel flow, where the maximum velocity in the channel corresponds to sonic velocity. Neutral stability contours were obtained for this flow as a function if the wave number,Reynolds number and the upper wall Mach number. The critical Reynolds number is found as 5718.338 for an upper wall Mach number of 0.0001, corresponding to the fully Poiseuille case.

In some species, including humans, red blood cells (RBCs) under low shear stress tend to clump together and form into regular stacks called rouleaux. These stacks are not static, and constantly move and break apart. This phenomenon is referred to as red blood cell aggregation and disaggregation. When modelled as a single liquid, blood behaves as a non-Newtonian fluid. Its viscosity varies, mainly due to the aggregation of RBCs. The aim of this research is to develop a mesoscale computational model for the simulation of RBCs in plasma. This model considers RBC interaction and aggregation to predict blood-flow characteristics such as viscosity, rouleaux size and velocity distribution. In this work, the population-balance modelling (PBM) approach is utilized to model the RBC aggregation process. The PBM approach is a known method that is used for modelling agglomeration and breakage in two-phase flow fluid mechanics to find aggregate size. The PBM model is coupled to the incompressible Navier-Stokes equations for the plasma. Both models are numerically solved simultaneously. The population-balance equation has been used previously in a more restricted form, the Smoluchowski equation, to model blood viscosity, but it has never been fully coupled with the Navier-Stokes equation directly for the numerical modelling of blood flow. This approach results in a comprehensive model which aims to predict RBC aggregate size and their velocities for different flow configurations, as well as their effects on the apparent macro-scale viscosity. The PBM approach does not treat the microscopic physics of aggregation directly but rather uses experimental correlations for aggregation and disaggregation rates to account for the effects of aggregation on the bulk. To find the aggregation rate, a series of experiments on RBC sedimentation due to gravity is designed. In these tests, aggregated RBCs (rouleaux) tend to settle faster than single RBCs and, due to low shear stresses, disaggregation is very low and can be neglected. A high-speed camera is used to acquire video-microscopic pictures of the process. The size of the aggregates and their velocities are extracted using image processing techniques. For image processing, a general Matlab program is developed which can analyze all the images and report the velocity and size distribution of rouleaux. An experimental correlation for disaggregation rate is found using results from a previous steady-state Couette flow experiment. Aggregation and disaggregation rates from these experiments are used to complete the PBM model. Pressure-driven channel flow experiments are then used for the final validation of the model. Comparisons of the apparent viscosity of whole blood in previous experiments show reasonable agreement with the developed model. This model fills a gap between micro-scale and macro-scale treatments and should be more accurate than traditional macro-scale models while being cheaper than direct treatment of RBCs at the micro-scale.

In this thesis I investigate a number of problems in the nonlinear stability of density stratified plane Couette flow. I begin by describing the history of transient growth phenomena, and in particular the recent application of adjoint based optimisation to find nonlinear optimal perturbations and associated minimal seeds for turbulence, the smallest amplitude perturbations that are able to trigger transition to turbulence. I extend the work of Rabin et al. (2012) in unstratified plane Couette flow to find minimal seeds in both vertically and horizontally sheared stratified plane Couette flow. I find that the coherent states visited by such minimal seed trajectories are significantly altered by the stratification, and so proceed to investigate these states both with generalised Koopman mode analysis and by stratifying the self-sustaining process described by Waleffe (1997). I conclude with an introductory problem I considered that investigates the linear Taylor instability of layered stratified plane Couette flow, and show that the nonlinear evolution of the primary Taylor instability is not coupled to the form of the linearly unstable mode, in contrast to the Kelvin-Helmholtz instability, for example. I also include an appendix in which I describe joint work conducted with Professor Neil Balmforth of UBC during the 2015 WHOI Geophysical Fluid Dynamics summer programme, investigating stochastic homoclinic bifurcations.

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. & 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 dir="ltr">The ability to predict and manipulate how a particulate material will flow in a process is challenging for industry and researchers alike. This dissertation presents the results of a model-directed, experimental approach using a concentric cylinder rheometer titled along an axis to enable continuous chute flow of granular media. Experiments were performed using draining flows for constant and oscillatory applied shear rates. Multiple flow and stress sensors were used to investigate the interaction of mass holdup, shear rate, specific torque, particle velocity, discharge mass flow rate, and wall pressure. Depending on the flow configuration, linear ranges were observed wherein the specific torque remained steady during draining. This finding enabled systematic testing of flow behavior as a function of dimensionless shear rates. Results suggest changes in the specific torque, wall slip, and outflow variance occur with the transition from the quasi-static to dense-inertial flow regimes. A pump-curve analogy was also identified for the relationship between the outlet mass flow rate and the specific power relationship for the constant shear rate experiments. Oscillatory shear rate experiments show a significant influence of the phase shift between the applied shear rate and the specific torque. Adding an asperity to the rotor revealed rate-dependent patterns in bulk flow and force chain dynamics. Overall, the study offers valuable insights into the effects of shear rate and boundary conditions on dense granular flows. The effects of particle characteristics (e.g., size and shape distributions, friction, cohesivity) and material properties (e.g., density, modulus) remain topics for future work. </p>

Implanter une machine électrique dans un véhicule hybride pose avant tout des problèmes d’encombrement. Sous-dimensionner la machine semble légitime compte tenu de l’usage qui en est fait sur véhicule. Par contre, cela suppose que les aspects thermiques soient pris en compte non seulement lors de l’utilisation, mais aussi lors de la conception de la machine. Le phénomène majeur limitant la densité de puissance massique des machines électriques est l’échauffement interne des bobines et des aimants. La modélisation thermique de la machine est complexe compte tenu de la diversité des sources de chaleur et de la coexistence de différents modes de transferts thermiques : conduction dans la matière, convection avec l’eau de refroidissement, conduction, convection et rayonnement dans l’entrefer. En termes de géométrie, si une première approche peut être réalisée en ne considérant que des flux de chaleur radiaux, la composante axiale doit nécessairement être prise en compte dès lors qu’on veut tenir compte aussi des extrémités de machine, et notamment de la chaleur produite par les roulements et les têtes de bobines. Ainsi pour pouvoir analyser pertinemment les transferts thermiques dans la machine, des méthodes numériques de type CFD ont été utilisées pour caractériser le transfert thermique par convection. La caractérisation des échanges thermiques par convection naturelle et forcée a été réalisée à l’intérieur d’une machine synchrone à aimants permanents internes (MAPI). Des relations empiriques ont été proposées afin de prendre en compte le couplage entre la thermique et l’aérodynamique dans les cavités des machines électriques totalement fermées. Afin de valider la pertinence des modèles numériques dans le cadre de ce travail, des mesures thermiques à l’aide des moyens d’essais de l’UTC ont été réalisées. Les résultats de cette étude sont utilisés pour construire des circuits thermiques équivalents qui prennent en compte les phénomènes thermiques complexes dans les machines électriques fermées utilisées dans les véhicules hybrides et électriques. Ces recommandations de conception permettront l’optimisation de l’effort investi pour le refroidissement de la machine électrique dans ses différentes phases de fonctionnement. / The integration of an electrical machine into modern hybrid vehicles is associated with new technical constraints such as the integrability into small volume without losing certainly in performance. Therefore, the development of compacter electrical machines is a well-founded argument for car manufacturers as well as electrical machine designers. On the other hand, this finding assumes that the thermal aspects are undertaken not only during the operation of the electrical machine, but also during the design process. The internal heat generated in different areas impacts strongly the power density and the magnet health which deeply reduce the electrical machine reliability. Heat transfer modeling inside electrical machines is a tricky task because of the strong coupling between the different physics governing their operations. Indeed, the generated losses spread inside the electrical machine through three heat transfer modes which are: conduction (heat diffusion), convection(heat transport) and radiation (heat scattering). In terms of geometry, if a first approach can be carried out by considering only radial heat fluxes, the axially-transferred heat must be undertaken when it is also necessary to consider end caps effects, and particularly the heat released by the bearings. In order to carry out relevantly the thermal analysis of a permanent magnets synchronous machine, CFD based methods are used to characterize the convective heat transfer inside this machine over a large operating range. Both natural and forced convection are analyzed and the corresponding heat transfer coefficients are numerically-estimated. Empirical equations are proposed in order to take into account the coupling between thermal and fluid dynamics inside the cavities of the studied totally-enclosed machine. These correlations are integrated then into a detailed and reduced thermal network. Experimental tests are carried out using a test bench in order to measure temperature distribution in different areas of the electrical machine. Afterward, a comparison between estimated and measured temperatures shows that the results of the numerically-enhanced thermal network are in a good agreement with measurements. Thus, the proposed recommendations based on CFD modeling allow the convective heat transfer to be characterize quickly and precisely. These correlations are useful for upcoming studies dealing with convection inside automotive electrical machines as well as high speed electrical machines.

Cette thèse est une étude expérimentale sur mélange par des écoulements anisotropes turbulents en fluide stratifié. Il s'inscrit dans un projet sur l'étude globale de la structure de l'océan, principalement de sa couche supérieure. Le mélange est étudié en détail, à partir d'expériences réalisées en bicouche de densité dans un dispositif Taylor - Couette. En régime laminaire, le transport dans les zones initialement homogènes est moins intense que le flux de matière arrachée au niveau de l'interface. De multiples interfaces de densité apparaissent et leur action combinée contrôle les flux verticaux de densité. En régime turbulent, on trouve une relation entre flux et gradient de densité qui prédit un flux maximum pour des nombres de Richardson RIO d'environ 10. Un regain d'efficacité est observé pour les très grandes valeurs de RIO. Ces résultats vont dans le sens de la théorie de Balmforth (1998). Pour ces fortes valeurs de RIO, l'observation de l'écoulement par L.I.F., indique une activité intense des ondes et la présence de structures turbulentes de très petites échelles.

Ce travail porte sur l'étude théorique des interfaces entre deux fluides visqueux, soumis à un écoulement de Couette plan. Dans cette situation hors d'équilibre, les fluctuations thermiques de l'interface sont modifiées en raison du couplage par le cisaillement entre les effets visqueux et les effets de tension. Comme c'est le cas pour d'autres systèmes de matière molle (par exemple, les phases lamellaires), le cisaillement peut alors amplifier ou amortir les déformations interfaciales. On s'intéresse tout d'abord à la dynamique des fluctuations interfaciales. On montre que ces dernières vérifient une équation stochastique non-linéaire, dont la solution est contrôlée par un paramètre sans dimension qui contient toute l'information sur le système. La résolution à faible taux de cisaillement révèle que le déplacement quadratique moyen des fluctuations thermiques diminue avec l'écoulement, conformément aux observations expérimentales et numériques. Ensuite, on étudie l'influence des effets inertiels sur la stabilité de l'écoulement, dans le régime des fortes viscosités et des faibles tensions. Ce régime des grands nombres capillaires n'a été que très peu étudié, mais trouve sa pertinence par exemple dans les mélanges biphasiques de colloïdes et de polymères. Des critères de stabilité simples sont mis en évidence. Finalement, on réalise une étude numérique des propriétés des fluctuations interfaciales à grand cisaillement. Bien que les effets visqueux soient dominants, il en ressort une phénoménologie similaire à certains modèles de turbulence.

Outbreaks of food-born illness associated with the consumption of unpasteurized juice and apple cider have resulted in a rule published by the U.S. Food and Drug Administration (FDA) in order to improve the safety of juice products. The rule (21 CFR120) requires manufacturers of juice products to develop a Hazard Analysis and Critical Control Point (HACCP) plan and to achieve a 5-log reduction in the number of the most resistant pathogens. Ultraviolet (UV) disinfection is one of the promising methods to reach this 5-log reduction of pathogens. The absorption coefficients of juices typically vary from 10 to 40 1/cm and can be even higher depending on brand and processing conditions. Thin film reactors consisting of two concentric cylinders are suitable for inactivating pathogens in juices. When the two concentric cylinders are fixed, the flow pattern in the gap can be laminar Poiseuille flow or turbulent flow depending on flow rates. If the inner cylinder is rotating, and the rotating speed of the inner cylinder exceeds a certain value, the flow pattern can be either laminar or turbulent Taylor-Couette flow. UV disinfection between concentric cylinders in laminar Poiseuille flow, turbulent flow and both laminar and turbulent Taylor-Couette flow was investigated experimentally and numerically. This is the first systematic study done on UV disinfection between concentric cylinders in all three flow patterns. The present work provides new experimental data for pathogen inactivation in each of the three flow patterns. In addition, the present study constitutes the first systematic numerical CFD predictions of expected inactivation levels. Proper operating parameters and optimum gap widths for different flow patterns are suggested. It is concluded that laminar Poiseuille flow provides inferior (small) inactivation levels while laminar Taylor-Couette flow provides superior (large) inactivation levels. The relative inactivation levels are: laminar Poiseuille flow < turbulent flow < laminar Taylor-Couette flow.

Cette thèse porte sur l'étude expérimentale de la réduction de traînée par injection de bulles. Le contexte de l'étude est lié à la propulsion navale. Les bulles peuvent, lorsqu'elles sont injectées dans la couche limite en développement le long des carènes de navires, contribuer à réduire significativement la résistance de frottement, en allégeant d'une part le fluide le long de la coque et d'autre part en interagissant avec les structures turbulentes de proche paroi. La configuration expérimentale retenue pour cette étude est l'écoulement de Taylor-Couette, avec cylindre extérieur fixe et injection de bulles calibrées. Nos investigations portent sur les régimes de transition, de turbulence naissante et turbulent avec persistance des cellules de Taylor (Re=20000). La taille des bulles varie avec le capillaire et le mélange utilisés, entre 0,05 et 0,12 fois la largeur de l'entrefer. Une méthode originale de suivi de la phase dispersée dans un plan méridien, couplée à des mesures de couple global appliqué au cylindre intérieur ont permis de mettre en évidence deux régimes de réduction de traînée et plusieurs types d'arrangements des bulles, en fonction de leur taille et du nombre de Reynolds. Les bulles peuvent avoir une trajectoire glissante, oscillante, être capturées par les cellules de Taylor ou en proche paroi du cylindre intérieur dans les zones de jets sortants. La caractérisation par PIV des vitesses de la phase liquide en monophasique et diphasique a permis d'étudier les modifications induites par les bulles sur la phase liquide et de discuter des mécanismes impliqués dans la modification du couple global par la présence des bulles. Il ressort de cette étude que pour les nombres de Reynolds en deçà de la capture, les bulles contribuent à stabiliser l'écoulement en accord avec une réduction du couple visqueux pouvant atteindre -30% pour des taux de vide très faible (< 1%). Pour des nombres de Reynolds plus élevés, la capture dans les cellules conduit à une réduction de la longueur d'onde axiale et une augmentation de la vorticité des cellules, associée à une augmentation des vitesses rms. Cette configuration est favorable à une augmentation du couple visqueux. A l'inverse, la capture des bulles dans le jet sortant conduit à une augmentation de la longueur d'onde axiale, associée à une diminution de la vorticité. Cette configuration est favorable à une réduction du couple visqueux, moins marquée qu'en absence de capture.