Th. univ.--Math. appl.--Nice-Sophia Antipolis, 1996. / Contient des chapitres en anglais. Bibliogr. p. 161-166. Résumés en français et en anglais. 1996 d'après la déclaration de dépôt légal.

Th. 3e cycle--Méc.-rhéol.--Toulouse--I.N.P., 1978. N°: 22.

Th. 3e cycle--Méc.--Toulouse--I.N.P., 1983. N°: 175.

Th. doct.-ing.--Bioméc.--Toulouse--I.N.P., 1983. N°: 238.

Th.--Méc. des fluides--Toulouse--I.N.P., 1979. N°: 43.

Trente-quatre essais d'écoulement provenant de quatre groupes d'auteurs ont été sélectionnés dans la littérature afin d'étudier l'instabilité interne de matériaux de faible plasticité soumis à l'écoulement de l'eau et initiée par suffosion. Les essais sélectionnés ont tous été effectués dans un perméamètre en appliquant les contraintes hydrodynamiques de manière ascendante et par palier. L'objectif était de déterminer la vitesse critique (vcr(l))p qui initie la mobilisation des particules fines dans les pores de la structure prédominante des matériaux, soit celle qui forme les constrictions. La méthode d'analyse présentée par Kovacs et Ujfaludi (1983), mettant en relation la résistance à l'écoulement ([Delta]h/Q) en fonction de la charge totale appliquée ([Delta]h), a été utilisée afin de déterminer les vitesses critiques (vcr(¹))p. Les essais ont tous été réinterprétés avec cette méthode et celle-ci a permis d'extraire des conditions critiques à des stades inférieurs durant les essais comparativement aux méthodes d'analyses utilisées par leurs auteurs respectifs. Ainsi, par cette approche, un nouveau critère de stabilité interne, basé sur celui développé par Kovacs et Ujfaludi (1983) pour l'application sur les puits, a été proposé pour l'étude des barrages en remblai. Puis, toujours suite à l'analyse des essais avec la méthode [Delta]h/Q - [Delta]h, les matériaux utilisés lors des essais ont été classés en trois catégories selon leur type de comportement sous écoulement après la vitesse critique (vcr(¹))p dans le repère [Delta]h/Q - [Delta]h. Il a ensuite été discuté que la sélection de matériaux stables pour l'étude des barrages doit tenir compte du type de comportement sous écoulement des matériaux et doit inclure un critère admissible de perte de fines, tel que présenté par l'approche russe. Enfin, via l'étude de certains paramètres qui influencent la vitesse critique, il a été démontré que la seule façon d'obtenir une tendance dans les résultats était de représenter dans un repère cartésien logarithmique les vitesses d'écoulement critiques (vcr(¹))p par rapport à la conductivité hydraulique (K) des matériaux.

Ce mémoire présente les résultats de la première étape d'une étude sur un jet circulaire turbulent faiblement chauffé et confiné. Ces travaux, réalisés d'août 2004 à mai 2006, ont mené au développement d'une nouvelle approche expérimentale pour l'étude des jets circulaires confinés. L'étude a été réalisée dans trois enceintes de surfaces différentes, permettant ainsi d'étudier le comportement d'un jet faiblement confiné (pour (An)05 = 42 D) à celui d'un jet moyennement confiné (pour (AR)0,5 = 30Det20D). Les profils de la vitesse et de la température moyennes du jet ont été mesurés et comparés avec des résultats expérimentaux issus de la littérature. Ces essais ont permis d'établir que le jet est autosimilaire à une distance d'environ x/D = 20 en aval de la sortie. La plupart des mesures ont été effectuées dans la région auto-similaire, c'est-à-dire à des distances supérieures à x/D = 20. L'objectif était de vérifier l'effet du nombre de Reynolds sur la diffusion du jet dans la région où la turbulence est pleinement développée. Le mémoire présente aussi une étude expérimentale dans la zone la plus proche de la sortie du jet dénommée "cône potentiel", ainsi que la mesure de l'ordre de grandeur des échelles de Kolmogorov dans la région auto-similaire du jet.

Des écoulements instationnaires au niveau du distributeur des turbines à réaction ont été expérimentalement étudiés dans ce projet doctoral. Le but est de développer une expertise intégrée pour mesurer des écoulements instationnaires au niveau des distributeurs des turbines à réaction en utilisant une vélocimétrie par image de particules (PIV) et d'autres techniques de mesure. Les expériences ont été effectuées en deux phases appelées configuration 1 et configuration 2. La configuration 1 est composée d'une directrice, installée dans un canal bidimensionnel d'eau, examinée aux nombres de Reynolds comparable à ceux dans les modèles standards de turbines Francis. La configuration 2 consiste en un distributeur d'un modèle réduit d'une turbine Francis. Dans la configuration 1, les champs d'écoulement ont été mesurés avec un PIV autour d'une directrice asymétrique dans un canal d'eau, à différents angles d'attaque entre 0° et 20°. Les données obtenues comprennent les champs de vitesse instantanée et moyenne. Pour les mêmes conditions, les pertes de charge dues à la directrice ont été mesurées et comparées. Dans la configuration 2, les champs d'écoulement entre deux directrices voisines ont été mesurés avec le PIV au travers d'un accès optique spécifique. Des résultats instantanés tels que les champs de vitesse et les résultats moyens tels que la fluctuation de vitesse et les iso-contours de vitesse ont été présentés. La distribution de pression du côté intrados et du côté extrados de deux directrices voisines a été mesurée avec des capteurs de pression miniatures montés sur les surfaces des directrices. 'Basé sur les mesures, les effets des points d'opération sur le comportement des écoulements sont analysés.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2005-2006 / Un modèle numérique tridimensionnel a été développé pour la simulation du système chimique quartz-eau couplé avec l’écoulement à densité et viscosité variable dans les milieux poreux discrètement fracturés. Le nouveau modèle simule aussi le transfert de chaleur dans les milieux poreux fracturés en supposant que l’expansion thermique du milieu est négligeable. Les propriétés du fluide, densité et viscosité, ainsi que les constantes chimiques (constant de taux de dissolution, constant d’équilibre, coefficient d’activité) sont calculées en fonction de la concentration des ions majeurs et de la température. Des paramètres de réaction et d’écoulement, comme la surface spécifique du minéral et la perméabilité sont mis jour à la fin de chaque pas de temps avec des taux de réaction explicitement calculés. Le modèle suppose que des changements de la porosite et des ouvertures de fractures n’ont pas d’impact sur l’emmagasinement spécifique. Des pas de temps adaptatifs sont utilisés pour accélérer et ralentir la simulation afin d’empêcher des résultats non physiques. Les nouveaux incréments de temps dépendent des changements maximum de la porosité et/ou de l’ouverture de fracture. Des taux de réaction au niveau temporel L+1 (schéma de pondération temporelle implicite) sont utilisés pour renouveler tous les paramètres du modèle afin de garantir la stabilité numérique. Le modèle a été vérifié avec des problèmes analytiques, numériques et physiques de l’écoulement à densité variable, transport réactif et transfert de chaleur dans les milieux poreux fracturés. La complexité de la formulation du modèle permet d’étudier des réactions chimiques et l’écoulement à densité variable d’une façon plus réaliste qu’auparavant possible. En premier lieu, cette étude adresse le phénomène de l’écoulement et du transport à densité variable dans les milieux poreux fracturés avec une seule fracture à inclinaison arbitraire. Une formulation mathématique générale du terme de flottabilité est dérivée qui tient compte de l’écoulement et du transport à densité variable dans des fractures de toute orientation. Des simulations de l’écoulement et du transport à densité variable dans une seule fracture implanté dans une matrice poreuse ont été effectuées. Les simulations montrent que l’écoulement à densité variable dans une fracture cause la convection dans la matrice poreuse et que la fracture à perméabilité élevée agit comme barrière pour la convection. Le nouveau modèle a été appliqué afin de simuler des exemples, comme le mouvement horizontal d’un panache de fluide chaud dans un milieu fracturé chimiquement réactif. Le transport thermohalin (double-diffusif) influence non seulement l’écoulement à densité variable mais aussi les réactions chimiques. L’écoulement à convection libre dépend du contraste de densité entre le fluide (panache chaud ou de l’eau salée froide) et le fluide de référence. Dans l’exemple, des contrastes de densité sont généralement faibles et des fractures n’agissent pas comme des chemins préférés mais contribuent à la dispersion transverse du panache. Des zones chaudes correspondent aux régions de dissolution de quartz tandis que dans les zones froides, la silice mobile précipite. La concentration de silice est inversement proportionnelle à la salinité dans les régions à salinité élevée et directement proportionnelle à la température dans les régions à salinité faible. Le système est le plus sensible aux inexactitudes de température. Ceci est parce que la température influence non seulement la cinétique de dissolution (équation d’Arrhenius), mais aussi la solubilité de quartz. / A three-dimensional numerical model is developed that couples the quartz-water chemical system with variable-density, variable-viscosity flow in fractured porous media. The new model also solves for heat transfer in fractured porous media, under the assumption of negligible thermal expansion of the rock. The fluid properties density and viscosity as well as chemistry constants (dissolution rate constant, equilibrium constant and activity coefficient) are calculated as a function of the concentrations of major ions and of temperature. Reaction and flow parameters, such as mineral surface area and permeability, are updated at the end of each time step with explicitly calculated reaction rates. The impact of porosity and aperture changes on specific storage is neglected. Adaptive time stepping is used to accelerate and slow down the simulation process in order to prevent physically unrealistic results. New time increments depend on maximum changes in matrix porosity and/or fracture aperture. Reaction rates at time level L+1 (implicit time weighting scheme) are used to renew all model parameters to ensure numerical stability. The model is verified against existing analytical, numerical and physical benchmark problems of variable-density flow, reactive solute transport and heat transfer in fractured porous media. The complexity of the model formulation allows chemical reactions and variable-density flow to be studied in a more realistic way than previously possible. The present study first addresses the phenomenon of variable-density flow and transport in fractured porous media, where a single fracture of an arbitrary incline can occur. A general mathematical formulation of the body force vector is derived, which accounts for variable-density flow and transport in fractures of any orientation. Simulations of variable-density flow and solute transport are conducted for a single fracture, embedded in a porous matrix. The simulations show that density-driven flow in the fracture causes convective flow within the porous matrix and that the highpermeability fracture acts as a barrier for convection. The new model was applied to simulate illustrative examples, such as the horizontal movement of a hot plume in a chemically reactive fractured medium. Thermohaline (double-diffusive) transport impacts both buoyancy-driven flow and chemical reactions. Free convective flow depends on the density contrast between the fluid (hot brine or cool saltwater) and the reference fluid. In the example, density contrasts are generally small and fractures do not act like preferential pathways but contribute to transverse dispersion of the plume. Hot zones correspond to areas of quartz dissolution while in cooler zones, precipitation of imported silica prevails. The silica concentration is inversely proportional to salinity in high-salinity regions and directly proportional to temperature in low-salinity regions. The system is the most sensitive to temperature inaccuracy. This is because temperature impacts both the dissolution kinetics (Arrhenius equation) and the quartz solubility.

Tableau d'honneur de la Faculté des études supérieures et postdorales, 2014-2015 / Cette thèse présente une étude de la turbulence et du scalaire passif dans la zone autosimilaire d’un jet rond libre faiblement chauffé, caractérisé par un haut nombre de Reynolds ([symbol]). D’abord, une analyse des limitations des sondes à fils chauds croisés stationnaires basée sur des mesures LDV est présentée. Celle-ci met en évidence l’importance des événements de vitesse à grands angles d’incidence dans les valeurs mesurées. Ensuite, la turbulence et la chaleur employée comme scalaire passif dans le jet sont étudiées à l’aide du LDV, de la thermométrie à fil froid et d’une combinaison de ces deux méthodes. On mesure ainsi les distributions radiales des valeurs moyennes et des moments de vitesse, de température et de leurs corrélations croisées. Ces statistiques incluent notamment [symbol]. Ces mesures sont caractérisées par un contrôle précis des conditions aux limites et une grande convergence statistique. On remarque un bon accord entres les valeurs mesurées de [symbol] , [symbol] ainsi que [symbol] et leur expressions analytiques, dérivées des équations de continuité, du mouvement moyen et d’enthalpie, respectivement. Les bilans de transports de l’énergie cinétique turbulente (k), du scalaire passif ([symbol] ) ainsi que des flux de chaleur turbulents [symbol] sont ensuite établis, dans le même écoulement de jet rond faiblement chauffé. Pour le bilan de k, le modèle proposé par Lumley (1978) est employé pour représenter les termes de diffusion de pression. Pour les deux bilans, l’intégration à zéro sur un volume de contrôle des termes de diffusion est vérifiée et la portée de ce critère de validité est analysée. Aussi, les restes des deux bilans constituent leurs termes de dissipation respectifs, ek et eq . Dans le cas de ek, le modèle prédictif mis de l’avant par Thiesset et al (2014) ainsi qu’une analyse originale des bases de données de Panchapakesan & Lumley (1993a) et de Hussein et al (1994) laisse entrevoir une possibilité de réconciliation en ce qui a trait aux résultats produits par ces deux expériences réputées. Les budgets des composantes du tenseur de Reynolds ([symbol] ) sont aussi établis et analysés. La valeur de eq sur l’axe du jet est employée comme référence dans une étude exhaustive des techniques de mesures directes de ce taux de dissipation. Cette étude expose clairement les difficultés de ces mesures, notamment par des analyses spectrales. Les données amassées permettent d’établir les distributions radiales de paramètres d’intérêt pour la modélisation, tels que la viscosité tourbillonnaire ([symbol] ), la diffusivité tourbillonnaire ([symbol] ) et le nombre de Prandtl turbulent ([symbol] ). L’assomption de [symbol] constant est examinée. Le ratio entre les échelles de temps de dissipation de température et d’énergie cinétique turbulente (R) est aussi calculé et analysé. / This thesis presents a study of the turbulence generated by a slightly heated turbulent round jet, evolving in still air. The study focuses on the self-preserving region of the high Reynolds number jet ([symbol] ). The jet is heated to only 20°C above ambient, so that temperature is treated as a passive scalar. In the first chapter of the thesis, a new and novel analysis of the stationary crossed hot wire probes based on unbiased LDV data is exposed, thus revealing the importance of the wide incidence angle events on the measured values. LDV, cold wire thermometry and a combination of both methods is used to study the slightly heated jet. Radial distributions of mean velocities and heating level, second and third order moments of fluctuating velocity and temperature ([symbol]) along with cross correlations, such as [symbol] are made available. During data acquisition, great attention was devoted to the close control of boundary conditions and statistical convergence. The measured values of V, uv as well as vq are in good agreement with their analytical expressions, provided by continuity, mean momentum and enthalpy equations, respectively. The transport budgets for the turbulent kinetic energy, ([symbol] ), passive scalar ([symbol] ), and temperature turbulent fluxes ([symbol] ) are then established using the data acquired in the very same, slightly heated, jet flow. The model introduced by Lumley (1978) for the pressure diffusion, is used as a surrogate for this term, in the case of the k budget. The integration to zero of the diffusion terms of both budgets over a control volume is verified. The significance of this validity criterion is analysed. The remainders of both budgets, are their respective dissipation terms, ek and eq . In the case of ek, the predictive model recently put forward by Thiesset et al (2014) along with a novel treatment of the databases of Panchapakesan & Lumley (1993a) and Hussein et al (1994), may help bring the outcome of the two famous experiments to a reconciliation. The Reynolds stress tensor component budgets were also established, leading to an analysis of the importance of the redistribution terms. The [symbol] value provided by the [symbol] budget is used as a reference in an exhaustive study of the means of dissipation measurement, on the centerline of the jet. This study reveals the shortcoming of the direct measurements of eq , namely thought spectral analysis. Radial distributions of turbulence modeling parameters, such as turbulent momentum viscosity ([symbol] ), turbulent thermal diffusivity ([symbol] ) and turbulent Prandtl number ([symbol] ) were established. The assumption of a constant PrT is scrutinized. The radial distribution of the thermal-to-mechanical timescale ratio (R) is calculated and analysed.