Nous avons développé une biopuce de microfluidique qui est capable de surveiller continuellement la population de cellules dans les biofilms en conditions d'écoulement laminaire bien contrôlées. Ce dispositif microfluidique est capable de modeler la formation des biofilms linéaires en utilisant une approche de flux basé sur un modèle. Les considérations de conception et méthodologie de fabrication d'un micro-bioréacteur à deux niveaux, inclus le flux basé sur un modèle (FT-μBR) qui génère un flux de croissance du biofilm entouré par les 3 côtés par un flux de confinement et inhibiteur de croissance. Grâce à une combinaison d'expériences et de simulations, nous avons évalué et exploité exhaustivement les paramètres de contrôle pour manipuler les dimensions du modèle de flux de croissance du biofilm. Ce dispositif est ensuite utilisé pour développer des modèles linéaires du biofilm avec des dimensions contrôlables. Une étude de validation de principe utilisant le dispositif démontre son utilité dans la réalisation des mesures de taux de croissance du biofilms dans différents environnements de force de cisaillement. Cela ouvre la voie à des études quantitatives sur les effets de l'environnement des cisaillements locaux sur les propriétés des biofilms et pour la synthèse d'une nouvelle génération de biomatériaux fonctionnels ayant des propriétés contrôlables. / We have developed a microfluidic biochip that is capable of continuously monitoring cell population in biofilms under well-controlled laminar flow conditions. This microfluidic device capable of patterning linear biofilm formations using a flow-templating approach. The design considerations and fabrication methodology of a two level flow-templating micro-bioreactor (FT-μBR) generates a biofilm growth stream surrounded on 3 sides by a growth inhibiting confinement stream. Through a combination of experiments and simulations we comprehensively evaluate and exploit control parameters to manipulate the biofilm growth template stream dimensions. The FT-μBR is then used to grow biofilm patterns with controllable dimensions. A proof-of-principle study using the device demonstrates its utility in conducting biofilm growth rate measurements under different shear stress environments. This opens the way for quantitative studies into the effects of the local shear environment on biofilm properties and for the synthesis of a new generation of functional biomaterials with controllable properties.

Nous étudions les phénomènes de ségrégation par taille dans des suspensions concentrées<br />de particules en écoulement. Les particules sont plus lourdes que le fluide, comme dans la<br />plupart des situations pratiques (boues, bétons, propergols solides. . .). Le matériau modèle<br />est transparent, et nous pouvons mesurer les champs de vitesse et de concentration des<br />espèces sous cisaillement.<br /><br />Dans les deux géométries étudiées (Couette et plan-plan), on observe que l'écoulement<br />principal (azimutal) est en général localisé. De plus, nous montrons l'existence systématique<br />d'une composante non azimutale des trajectoires.<br /><br />Expérimentalement, nous établissons les relations entre figure de ségrégation, écoulements<br />primaire et secondaire. Sur cette base, nous discutons les mécanismes possibles de<br />la ségrégation.

Sous la direction de SNECMA, un nouveau moteur cryotechnique pour Ariane 5 est en phase de développement. Ce moteur, VINCI, fonctionne à l'oxygène et à l'hydrogène liquides et devra être rallumable en cours de vol. Le transitoire de remplissage d'une cavité intermédiaire par laquelle transite l'oxygène liquide avant d'arriver dans la chambre de combustion s'avère être une étape critique qu'il faut s'efforcer de maîtriser. Cette cavité, appelée dôme LOX, est directement reliée à la chambre de combustion par les injecteurs à oxygène. Des outils numériques sont actuellement en cours de conception : ils permettront à terme de simuler le remplissage de cette cavité dans l'espace. Afin de valider ces outils numériques, un programme expérimental a vu le jour, impliquant SNECMA, le CNES (Centre National d'Etudes Spatiales) et le LEGI (Laboratoire des Ecoulements Géophysiques et Industriels). Il s'agit de mener à bien des expériences " simples " et reproductibles afin de disposer d'une base de données expérimentales qui servira de cas tests pour les simulations. Un banc d'essais expérimental a été progressivement mis en place et instrumenté au LEGI. Le fluide utilisé en substitution de l'oxygène liquide est de l'eau. On a choisi de respecter une similitude du nombre de Weber entre le cas " réel " et les conditions expérimentales. Ce dernier équivaut au rapport des forces d'inertie sur les forces de tension de surface. Deux campagnes d'essais ont été réalisées, dans lesquelles on s'est centré sur l'étude des aspects hydrodynamique du transitoire de remplissage : variation du débit total et des pressions au cours d'un essai, évaluation du taux de vide dans la cavité, visualisation de l'écoulement dans la cavité et en sortie des injecteurs. L'instrumentation à disposition est constituée d'un débitmètre à effet Coriolis, de capteurs de pression, d'une sonde optique, de caméras rapides et d'un laser pour l'imagerie. La première campagne d'essais a visé le remplissage de la cavité en eau seule. Le paramètre de contrôle principal était la pression génératrice de l'écoulement liquide. Dans la deuxième campagne d'essais on injecte simultanément dans la cavité un écoulement d'eau et un écoulement d'air. Cela se rapproche plus des conditions réelles du transitoire de remplissage, au cours duquel la cavité est balayée par un écoulement d'hélium. Le paramètre de contrôle supplémentaire est le débit de gaz initial. Ces campagnes ont également souligné l'importance du profil d'ouverture de la vanne de l'écoulement liquide sur le transitoire de remplissage de la cavité. Ces campagnes d'essais constituent une première étape dans la compréhension du transitoire de remplissage du dôme LOX. Elles ont permis de visualiser la forme de l'écoulement dans la cavité et en sortie des injecteurs et d'identifier certains phénomènes intervenant dans le remplissage de la cavité. En particulier, nous avons mis en évidence l'existence d'un délai de mise en place de l'écoulement par les injecteurs, qui peut être responsable d'un pic de pression dans la cavité au cours du transitoire. L'influence de la fraction gazeuse sur l'écoulement dans les injecteurs a été soulignée mais reste à quantifier de façon précise. La prochaine étape de l'étude concerne les aspects énergétiques du transitoire de remplissage, notamment les transferts thermiques ayant lieu, dans la réalité, entre l'hélium, l'oxygène et les parois du dôme LOX.

Le développement des hydroliennes se base souvent sur des conditions d'écoulement idéalisées qui ne reflètent pas entièrement le courant présent sur un site réel. Dans cette optique, ce mémoire investigue l'effet d'un écoulement non aligné ou cisaillé sur les performances et le chargement de deux types d'hydrolienne : les hydroliennes à rotor axial (HRA) et à aile oscillante (HAO). Cette étude est réalisée à l'aide de simulations numériques et les résultats démontrent que pour les deux types d'hydrolienne, le désalignement produit une réduction de la puissance et de la traînée, tandis que le cisaillement n'a qu'un effet de second ordre sur ces derniers. Le chargement additionnel sur la structure de support causé par des conditions d'écoulement non idéal est aussi similaire pour les deux technologies. Le désalignement et le cisaillement affectent cependant plus significativement le chargement en fatigue des pales de l'HRA que celui de l'aile de l'HAO. / Le développement des hydroliennes se base souvent sur des conditions d'écoulement idéalisées qui ne reflètent pas entièrement le courant présent sur un site réel. Dans cette optique, ce mémoire investigue l'effet d'un écoulement non aligné ou cisaillé sur les performances et le chargement de deux types d'hydrolienne : les hydroliennes à rotor axial (HRA) et à aile oscillante (HAO). Cette étude est réalisée à l'aide de simulations numériques et les résultats démontrent que pour les deux types d'hydrolienne, le désalignement produit une réduction de la puissance et de la traînée, tandis que le cisaillement n'a qu'un effet de second ordre sur ces derniers. Le chargement additionnel sur la structure de support causé par des conditions d'écoulement non idéal est aussi similaire pour les deux technologies. Le désalignement et le cisaillement affectent cependant plus significativement le chargement en fatigue des pales de l'HRA que celui de l'aile de l'HAO. / The development of hydrokinetic turbines is often based on idealized flow conditions which do not fully refect river or tidal currents. In this regard, this thesis investigates the effect of non-aligned or sheared flows on the performances and loading of two turbine types: the axial rotor (ART) and oscillating foil (OFT) turbines. This study was conducted with unsteady numerical simulations and the obtained results show that for both turbine types, misalignment produces a reduction in power and thrust, while the shear has a limited effect on those quantities. The additional loading on the support structure caused by the non-idealized flow is also similar for both devices. However, misalignment and shear affect more severely the fatigue loading on ART blades than they do for OFT. / The development of hydrokinetic turbines is often based on idealized flow conditions which do not fully refect river or tidal currents. In this regard, this thesis investigates the effect of non-aligned or sheared flows on the performances and loading of two turbine types: the axial rotor (ART) and oscillating foil (OFT) turbines. This study was conducted with unsteady numerical simulations and the obtained results show that for both turbine types, misalignment produces a reduction in power and thrust, while the shear has a limited effect on those quantities. The additional loading on the support structure caused by the non-idealized flow is also similar for both devices. However, misalignment and shear affect more severely the fatigue loading on ART blades than they do for OFT.

Des données de terrain et des études de modélisation ont montré que la recharge d'eau de fonte sous les calottes glaciaires peut avoir un impact important et durable sur l'écoulement des eaux souterraines. En Estonie, au nord du Bassin Artésien Balte (BAB), ce mécanisme de recharge est invoqué pour expliquer la présence d’importants volumes d'eaux souterraines marquées par un signal isotopique et géochimique glaciaire caractéristique, étant donné que la région a connu plusieurs glaciations durant le Pléistocène et a été entièrement recouverte par la calotte Fennoscandienne au cours du Dernier Maximum Glaciaire (DMG), il y a 20000 ans environ. Cette étude vise à tester cette hypothèse à l’aide de simulations numériques. En premier lieu, une étude conceptuelle a été effectuée pour déterminer quels processus sous-glaciaires doivent être représentés dans un modèle numérique qui reproduise adéquatement les écoulements souterrains et le transport de solutés. Les processus suivants ont été étudiés: la recharge sous-glaciaire d'eau de fonte, la déformation poroélastique du milieu poreux sous le poids de la glace, l’isostasie, l’évolution du drainage en surface, le permafrost et les écoulements densitaires impliquant des eaux douces de fonte et des saumures profondes. Ces processus ont été simulés dans un modèle représentant un bassin sédimentaire conceptuel, au cours d'un épisode glaciaire suivi d'une période postglaciaire. Le transport de trois traceurs d’eau glaciaire a été simulé: δ18O, solides dissouts et âge de l’eau. Les résultats montrent que la simulation de la recharge sous-glaciaire avec une condition-limite de type 1 (Dirichlet) n'est pertinente que pour des flux de faible amplitude, ce qui pourrait être le cas sous des calottes glaciaires dont la base n’est que partiellement en fusion. La compression de la matrice rocheuse diminue les surpressions, qui apparaissent uniquement dans les couches à faible diffusivité hydraulique et épaisses. Si la recharge sous-glaciaire est faible, la compression de la matrice rocheuse peut entraîner des sous-pressions après le retrait de la calotte glaciaire. L’isostasie réduit considérablement l'infiltration d'eau de fonte et les écoulements d'eau souterraine. Sous la couche de pergélisol, l'écoulement des eaux souterraines est réduit en-dessous de la calotte glaciaire mais augmente en région périglaciaire. Tenir compte des variations de densité en lien avec la salinité diminue l'infiltration d'eau de fonte en profondeur. Cette étude montre que chaque processus sous-glaciaire est potentiellement important et devrait être pris en compte dans des modèles d’écoulement des eaux souterraines et de transport de solutés en milieu sous-glaciaire. Cependant, il est raisonnable de ne représenter que la recharge sous-glaciaire si les informations manquent pour décrire correctement les autres processus. Par conséquent, ce seul processus a été simulé pour reproduire les écoulements d'eau souterraine sous la calotte Fennoscandienne dans le BAB. Les simulations ont été réalisées dans deux modèles 2D verticaux, afin de vérifier si la recharge sous-glaciaire d’eau de fonte peut expliquer la distribution particulière de δ18O (un traceur d’eau de fonte) dans les eaux souterraines de la région. L’un recoupe l’Estonie, l’autre la Lettonie et les îles estoniennes dans le Golfe de Riga. L'écoulement des eaux souterraines est simulé durant 28000 ans, depuis le DGM jusqu’à aujourd’hui, de même que le transport de δ18O pour tracer l'eau de fonte et confronter les résultats des simulations avec les données de terrain. L'espace d’incertitude de certains paramètres a été exploré, comme l’intensité et la durée de la recharge sous-glaciaire, ainsi que la composition isotopique initiale de l'eau de fonte. Les simulations fournissent un ajustement satisfaisant entre les valeurs observées et calculées de δ18O, confirmant l’hypothèse que le BAB a subi une phase de recharge sous-glaciaire durant le DMG. Elles montrent que la recharge sous-glaciaire a créé une inversion de l'écoulement des eaux souterraines dans le bassin. L’eau de fonte a infiltré tous les aquifères, en particulier les aquifères non confinés. Après le retrait de la calotte Fennoscandienne, l'eau de fonte a été entièrement remplacée par de l'eau météorique moderne, excepté dans les aquifères confinés où de l’eau de fonte a été préservée à proximité des zones de décharge. Par ailleurs, d’importants volumes d'eau de fonte sont probablement préservés sous la mer Baltique. Les simulations indiquent enfin que des épisodes de recharge sous-glaciaire antérieurs au DGM doivent être considérés afin d'expliquer les valeurs de δ18O dans la partie plus profonde du bassin. / Field evidence and modelling studies have shown that subglacial recharge of meltwater under wet-based ice-sheets can have a significant and long-lasting impact on groundwater flow. In the northern Baltic Artesian Basin (BAB), in Estonia, this mechanism of recharge is thought to be responsible of the presence of large volumes of groundwater with a characteristic glacial isotopic and geochemical signal, because the region experienced several glaciations during the Pleistocene and was entirely covered by the Fennoscandian ice-sheet during the Last Glacial Maximum (LGM), some 20 ky BP. The present study aims at testing this hypothesis by means of numerical simulations. First, a conceptual numerical study was performed to determine which glacial and subglacial processes need to be represented in numerical models for adequately capturing subglacial groundwater flow dynamics and solute transport. The relevance of the following processes was studied: subglacial recharge of meltwater, poroelastic deformation of the porous medium under ice-sheet loading, isostasy, evolution of surface drainage, permafrost, and density-dependent flow involving fresh glacial meltwater and deep brines. Simulations of these processes were conducted in a generic sedimentary basin during a single glacial event followed by a postglacial period. The transport of three common tracers of subglacial recharge was simulated: δ18O, TDS, and groundwater age. Results show that simulating subglacial recharge with a fixed flux boundary condition is relevant only for low fluxes, which could be the case under partially wet-based ice-sheets. Glacial loading decreases overpressures, which appear only in thick and low hydraulic diffusivity layers. If subglacial recharge is low, glacial loading can lead to underpressures after the retreat of the ice-sheet. Isostasy considerably reduces the infiltration of meltwater and the groundwater flow rates. Below permafrost, groundwater flow is reduced under the ice-sheet but is enhanced beyond the ice-sheet front. Accounting for salinity-dependent density reduces the infiltration of meltwater at depth. This study shows that each glacial process is potentially relevant in models of subglacial groundwater flow and solute transport. However, representing only subglacial recharge can be a reasonable assumption if information is missing to describe the other processes properly. Therefore, this single process is simulated to reproduce groundwater flow beneath the Fennoscandian ice-sheet in the northern BAB. Simulations are performed in two cross-sectional models, in order to check whether subglacial recharge of meltwater can explain the unusual distribution of δ18O in groundwater in the region, which serves as a tracer of glacial meltwater. One model crosses Estonia, the other crosses Latvia and Estonian islands in the Gulf of Riga. Groundwater flow is simulated over 28 ky, from the Last Glacial Maximum (LGM) to present-day, along with δ18O transport for tracing meltwater and to compare the results of the simulations with field data. Parameter space exploration of subglacial recharge conditions is used to tackle the uncertainty in the intensity and duration of subglacial recharge in the northern BAB, as well as in the isotopic composition of meltwater. Simulations provide a satisfying fit between the observed and the computed values of δ18O, supporting the idea that subglacial recharge happened in the northern BAB during the LGM. Simulations show that subglacial recharge created a flow reversal in the basin. Meltwater infiltrated into all aquifers, especially the shallow ones. After the retreat of the Fennoscandian ice-sheet, meltwater was entirely replaced by modern meteoric water, excepted in confined aquifers where some meltwater has been preserved close to the discharge areas. Large volumes of meltwater are also probably preserved beneath the Baltic Sea. Simulations also indicate that episodes of subglacial recharge prior to the LGM must be considered in order to explain the values of δ18O in the deeper basin.

Une étude hydrogéologique conjointe de l'Université Laval et de la Commission Géologique du Canada (CGC) a été réalisée dans le bassin versant de la rivière Chaudière. Le but de cette étude est d’accroître les connaissances sur le système hydrogéologique du bassin versant de la rivière Chaudière qui s’étend de la frontière américaine jusqu’au fleuve Saint-Laurent. La zone étudiée a une superficie de 6675 km2. Il s’agit d’une région densément peuplée (~ 390,000 hab.) à vocation principalement agricole où 65 % de la population utilise l’eau souterraine comme principale source d’eau potable. L’aquifère régional se compose de deux unités principales : le substratum rocheux de la province géologique des Appalaches et certaines unités quaternaires granulaires. Le confinement de l’unité rocheuse est contrôlé par des sédiments glaciaires, marins ou lacustres fins d’épaisseurs et de distribution variables. L’écoulement dans le roc se fait dans sa partie supérieure où la fracturation est la plus importante. La première étape des travaux a d’abord consisté à définir les propriétés hydrauliques des formations géologiques. Les données utilisées proviennent de divers rapports, de données du système d’information hydrogéologique (SIH), du Ministère du Développement Durable, de l’Environnement, de la Faune et des Parcs du Québec (MDDEFP), de rapports de consultants ainsi que d’une campagne de travaux de terrain effectuée à l’été 2007. Parmi ces données, la conductivité hydraulique (K) constitue la principale propriété à définir. Plusieurs valeurs de K ont été extraites directement ou indirectement de rapports existants. Une partie importante des valeurs de la conductivité hydraulique a été estimée à partir de la base de données SIH. Parmi celles-ci, 1333 valeurs ont été évaluées à partir de la capacité spécifique (Bradbury et al, 1985); 89 de ces mesures ont étés qualifiées de données de qualités car elles respectent toutes les conditions de la méthode, à l’exception du temps de pompage souvent trop court (94% des cas). En se basant sur les mesures recueillies, la moyenne géométrique des données de conductivité hydraulique (K) s’élève à 5 x 10-6 m/s avec un écart-type du logarithme de K égal à 0.60. La lithologie ne semble pas avoir une grande influence sur les valeurs de la conductivité hydraulique car toutes les valeurs moyennes de chacune des lithologies se trouvent à l’intérieur d’un même ordre de grandeur. Un modèle numérique a été construit pour l’ensemble de la zone d’étude à l’aide du logiciel de modélisation par éléments finis Watflow (Molson et al., 2002) afin d’interpréter le système d’écoulement de l’eau souterraine au niveau régional. Le modèle tridimensionnel est composé de 28 couches dont 18 représentent les dépôts meubles et 10 le roc. Le maillage 3D comprend 1 / A joint hydrogeological assessment study between Université Laval and the Geological Survey of Canada (GSC) has been conducted in the Chaudière River watershed. The main objective of the study was to gain further understanding into the groundwater resources within the region. The watershed extends over ~ 6700 km2 from the United States border northwards to the Saint Lawrence River in the province of Québec (Canada). Groundwater is an important source of drinking water as well as for irrigation and industrial use, however the supply is finite and there are potential risks of conflicts arising from different users. It is thus important to understand the aquifer systems throughout the watershed in order to avoid possible conflicts and to help decision makers to better manage the resource. Within the catchment, up to 65% of the population relies on groundwater as the primary drinking water resource, whereby most of the water is used for agricultural purposes. The regional aquifer system is dominated by the underlying rock substratum of the Appalachian province and Quaternary aquifer units. The degree of confinement of the fractured aquifers is controlled by the overlying glacial, marine and fine lacustrine sediments of variable thicknesses. The groundwater circulation within the fractured aquifers predominantly takes place in the uppermost part of the fractured water-bearing units, where fracture density is highest. As a first step, the hydraulic properties of the fractured basement were defined, for which the hydraulic conductivity (K) of the formations can be considered as the most important hydraulic characteristic. A literature review of the available reports as well as the SIH and the MDDEFP data bases was performed in order to assemble the data. Additional field work was completed by the GSC during the spring of 2007 to supplement the data base. Most of the values of the hydraulic conductivity K were taken from the existing reports while others were estimated from the SIH data base. The data base contains 1333 values which were evaluated using Bradbury’s method, of which 83 satisfied the “quality criteria”, since they respected suggested conditions regarding the minimum pumping (which was not met in 94% of the cases). The geometric mean value of the measured hydraulic conductivities was 5x10-6 m/s with a standard deviation of 0.60. The maximum difference between the hydraulic conductivities of the different lithological units was about a factor of 100x. A 3D numerical finite element model was built with the Watflow model (Molson et al., 2002) for simulating the regional groundwater flow system. The model domain extends over the entire watershed and is subdivided into 28 layers, 18 for the unconsolidated deposits and 10 within the basement rock. The 3D triangular prismatic mesh contains 1,694,672 nodes and 3,344,516 elements. The topographic DEM (digital elevation model) forms the uppermost surface of the model. For the model boundary conditions, the Chaudière River and major tributaries and lakes are set as first-type boundaries with constant hydraulic heads. The model base, as well as the lateral limits, are considered as no-flow boundaries. The net recharge over the model surface was estimated using the infiltration model HELP, combined with estimates of the surface-distributed local water withdrawals. The major pumping wells are represented in the model domain as point sinks. The horizontal hydraulic conductivity together with the vertical anisotropy values for each formation were used as the model calibration parameters, using the observation data from 68 observation wells. The results of the calibration of the model showed that the mass balance error is close to zero which implies that the sum of inflows is equal to the sum of outflows. Different scenarios with recharge variations were simulated with the calibrated model. The simulations allowed identifying the principal sensitive zones which need piezometric and water quality control. It was also possible to determine the areas most sensitive to climate variations. Changes in the rate of pumping and recharge (± 20%) did not have a significant impact on the modeled system mainly in relation to groundwater level (observed variations were ≤5m).

Un modèle numérique a été utilisé pour simuler l'écoulement transitoire à saturation variable dans une croûte argileuse fracturée afin d'étudier les fluctuations des pressions d'eau dans une pente susceptible aux glissements superficiels. Le modèle hydrogéologique considère l'écoulement de surface et l'évapotranspiration. La caractérisation du site expérimental a servis à définir les propriétés hydrauliques d'une croûte argileuse. Les méthodes du profil instantané et de Köhne et al. (2002) ont été utilisées pour calibrer les courbes hydrauliques non saturées avec le modèle de double continuum de Gerke et van Genuchten (1993). Les résultats ont démontré que le réseau de fractures d'un sol argileux influence grandement la courbe K([psi]), mais peu [theta]([psi]). Les simulations montrent que les plus grandes variations du niveau d'eau sont généralement observées à la fin de la pluie torrentielle vis-à-vis la zone où se développent les glissements superficiels et correspondent à l'épaisseur de la frange capillaire.

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids.<p>We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a subtriangulation of the triangular Pk elements where we apply the distribution used for a P1 element.<p>This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes.<p><p>The first part is subdivided in four chapters:<p>The aim of the second chapter is to present the multidimensional upwind residual distributive schemes and to explain what was the status of their development at the beginning of this work.<p>The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes.<p>The fourth chapter is composed of two parts: we start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation.<p>This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion.<p>Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems.<p>In the second part, we extend the schemes of the first part to systems of equations as follows:<p>The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries.<p>Then, we look at the performance of 2nd and 3rd order schemes on viscous flow.<p>Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Les particules soumises à un écoulement peuvent manifester des orientations préférées et une variété de déformations en fonction de leur géométrie et élasticité et du champ de vitesse de l'écoulement. A l’inverse, les flux peuvent être modifiés lorsque les contraintes des particules sont non négligeables. Cette thèse présente des résultats théoriques et numériques sur cette relation bilatérale en deux parties. La Ière partie commence par une analyse de stabilité et des simulations numériques qui montrent une particule brownienne semi-flexible dans un écoulement élongationnel effectuant un retournement, un phénomène associé aux flux de cisaillement. Le Chap. 2 étend les outils analytiques dédiés aux modèles simples ou aux flux indépendante du temps pour les modèles perle-barre-bond généraux dans les flux aléatoires. En partant des résultats des chapitres précédents, Le Chap. 3 aboutit à l'étude d'un degré de liberté lagrangien inexploré dans un écoulement turbulent : la flexion. Une particule semi-flexible courbe différemment dans les flux aléatoires bidimensionnels et tridimensionnels. La IIème partie concerne la turbulence élastique, un régime chaotique créé dans un écoulement de faibles forces inertielles par l'ajout de polymères élastiques. Le Chap. 4 estime le nombre de degrés de liberté d'une solution de ce régime via la dimension Lyapunov de l'attracteur du modèle Oldroyd-B bidimensionnel, un modèle connu de reproduire la turbulence élastique dans les simulations numériques. Le Chap. 5 pose des questions sur la nécessité d'élasticité pour produire un régime chaotique et conclut qu'une solution de polymère de barres peut créer un régime comparable à la turbulence élastique. / Particles, when subjected into a flow, may display preferred orientations and a variety of deformations depending on their geometry and elasticity and the flow velocity field. Flows can conversely be modified when the particle stresses are sufficiently large. This thesis presents theoretical and numerical results on this two-way relationship between particles and flows in two parts. Part I starts with a stability analysis and numerical simulations that show a simple semiflexible Brownian particle in an extensional flow undergo tumbling, a phenomenon normally associated to shear flows. Chapter 2 extends analytical tools available only for elementary polymer models or for steady flows to general bead-rod-spring models in random flows. By building on the results from the previous chapters, Chap. 3 culminates with the study of an unexplored Lagrangian degree of freedom in a turbulent flow: bending. A semiflexible particle is shown to display different bending behaviours in two- and three-dimensional random flows. This prediction is confirmed via direct numerical simulations of the particle in a turbulent flow. Part II concerns “elastic turbulence", a chaotic regime created in a flow with low inertial forces by the addition of elastic polymers. Chapter 4 provides an estimate for the number of degrees of freedom of a solution of this chaotic system via the Lyapunov dimension of the attractor of the two-dimensional Oldroyd-B model, a model known to reproduce elastic turbulence in numerical simulations. Chapter 5 questions the necessity of elasticity in producing a chaotic regime and concludes that a rodlike polymer solution can create a regime similar to elastic turbulence.

Cette thèse porte sur l'écoulement à l'intérieur d'une turbine de type Hélice. Deux méthodes de mesures ont été utilisées afin de remplir les objectifs de l'étude. La précision des instruments et l'incertitude des mesures effectuées à l'aide d'outils de vélocimétrie laser tels que la vélocimétrie Doppler laser (LDV) et la vélocimétrie par images de particules (PIV) ont été analysées. Trois phénomènes principaux sont étudiés à l'aide du LDV et du PIV. Le premier porte sur le débalancement circonférentiel de l'écoulement dans toute la machine. Le second porte sur les interactions rotor-stator se produisant entre les directrices et les aubes de la roue. Le troisième englobe les écoulements à l'entrée de l'aspirateur, le vortex à charge partielle (torche) et à forte charge ainsi que le sillage des directrices. Une attention particulière est portée sur les fréquences caractéristiques trouvées à l'aide des mesures LDV à cet endroit. Parallèlement aux mesures expérimentales, des simulations numériques ont été effectuées et comparées avec les données expérimentales disponibles. Cette comparaison permet d'améliorer les simulations numériques en posant différentes conditions limites et modèles de turbulence. Une meilleure utilisation et une amélioration des modèles RANS, URANS et LES avec l'utilisation de la base de données expérimentales est discutée dans l'étude. En conclusion à la recherche, le débalancement circonférentiel de l'écoulement est identifié avec les méthodes expérimentales et numériques. Le débalancement est de l'ordre de 10% de la vitesse de référence crête à crête aux avants directrices et est amorti à 5% à l'interface. Le débalancement maximal est de 12% à l'intérieur du cône divergent de l'aspirateur. L'intensité turbulente dans les canaux d'admission a été estimée à 8%. Les résultats expérimentaux ont démontré que les interactions rotor-stator sont faibles à l'intérieur de ce type de machine car il y a une dissipation rapide du sillage des directrices avant la roue. La torche à charge partielle et les recirculations se produisant dans l'aspirateur sont présentées avec leurs effets sur l'efficacité hydraulique et la configuration de l'écoulement. Plusieurs structures et vortex sont identifiés tels que le vortex se produisant dans l'espace en porte-à-faux des directrices, le vortex de bout d'aubes, la torche et les vortex contrarotatifs des canaux d'admission. Une base de donnée a été crée pour des recherches futures portant sur l'écoulement d'une turbine Hélice.