Dans le cadre des écoulements à l’échelle du pore d’un milieu poreux, cette thèse de calculscientifique propose une méthode itérative de résolution de Stokes incompressible basée sur lapénalisation de la matrice solide et mettant en jeu le champ de vorticité. La méthode développéeici s’inspire des méthodes de vortex habituellement utilisées pour les écoulements externes où letransport est dominant, mais adaptée aux écoulements visqueux à petites échelles (bas Reynolds).Les propriétés de cette méthode numérique itérative, qui alterne et sépare les étapes de pénalisationet de diffusion, permettent de simuler des écoulements confinés dans l’espace poral enproposant une formulation ne faisant intervenir que des opérateurs usuels, permettant de s’affranchirde l’aspect non-séparable de l’équation de diffusion-pénalisation. On assure ainsi dessimulations à haute résolution et faible stockage mémoire, aussi bien sur des géométries simplesde référence que sur des géométries poreuses réelles, et les résultats sont confrontés à des donnéesexpérimentales ou théoriques. Le code de calcul implémenté est validé et exploité pour l’estimationde la perméabilité d’échantillons réels.Dans la continuité de ce travail sur les écoulements avec adhérence à l’interface, on généralisecette approche à la problématique de la rugosité des parois solides. On présente une méthode derésolution du problème de Stokes perturbé par ce type de rugosité à l’interface pores/matricemodélisé par une condition aux limites de glissement tangentiel de type Robin-Fourier.La thèse présente également un volet sur les écoulements à viscosité variable. On présenteainsi une simulation du problème de Stokes non-linéaire couplé avec une équation de diffusionet transport résolue par méthode particulaire. On s’intéresse à une simulation originale de l’écoulementde xanthane (modélisé via une loi de Carreau comme polymère rhéofluidifiant) dans une géométrie réelle de roche de type grès sableux Bentheimer. Les résultats présentés montrent lacapacité de la méthode à gérer des variations de viscosité importante. Également à propos de lamodélisation et de la simulation des écoulement à viscosité variable, et plus précisément de l’homogénéisationdes mélanges, on présente aussi un estimateur théorique de la viscosité effectivepour les mélanges aléatoires. Ce résultat est validé numériquement grâce à une méthodologie detype Monte-Carlo sur un grand nombre de simulations de Stokes. / This scientific computing work presents an iterative method for the solution of incompressibleStokes equation at the pore scale of porous media. This method is based on solid matrixpenalization and involves vorticity field. We adapt tools from vortex methods (generally used forexternal and transport-dominant flows) to viscous flows at low scale.The pore space flow is computed with this iterative method, alternating penalization anddiffusion steps. Only usual operators are used in the proposed formulation, such that the nonseparableproperty of the penalization-diffusion initial equation is overcome. We then performhigh resolution and low memory storage simulations on various geometries, including real porousmedia. We validate the code with these results, and by estimating permeability of samples.From this work which consists in no-slip flow at solid interface, we propose a generalizationto rough solid walls at the interface. This thesis presents the modeling of such roughness usinga tangential slip Robin condition, and the ensuing method as a computation of the resultingperturbation.This work also includes a study of variable viscosity flows. With a particle-based method,we present the numerical solution of a non-linear Stokes problem coupled with diffusion andtransport for xanthan flow in a Bentheimer sandstone. The viscosity law is a Carreau law, asxanthan solution is a shear-thinning substance : the method is then able to handle high viscosityvariations. About the same topic of variable viscosity flows, and more precisely about mixturehomogenization, we introduce a theoretical estimator of effective viscosity for random mixtures.This estimator is validated on a large number of Stokes simulation, with a methodology inspiredby Monte-Carlo methods.

Le but de cette recherche consiste à la modélisation numérique et expérimentale d’écoulement sanguin dans les artères, en utilisant les nouvelles méthodes de couplage fluide structure. Plusieurs approches numériques peuvent être utilisées pour ce type de couplage, la méthode des éléments finis pour la modélisation de la structure. Pour la modélisation du domaine fluide la méthode de type éléments fnis ou méthode particulaire, Smooth Particle Hydrodynamic method, (méthode SPH) peuvent utilisées. Afin de valider ces deux méthodes numériques, Eléments Finis et SPH, une première application consiste à la modélisation numérique du gonflement d’une membrane en caoutchouc. Pour cette application des données expérimentales sont disponibles pour la validation des résultats numériques. Pour le matériau caoutchouc, une loi de comportement de type Mooney Rivlin est utilisée. Pour la structure la méthode des éléments finis, formulation coque est utilisée, pour le fluide nous avons opté pour la méthode particulaire SPH. Concernant la sensibilité par rapport au maillage, plusieurs maillages pour la modélisation des domaines fluide et structure sont testés. Une bonne corrélation en terme de déplacement et de vitesse du centre de la membrane, entre les résultats numériques et données expérimentales a été observée. Cette application est publiée dans un journal international de rang A.La seconde formulation développée dans le manuscrit consiste à la modélisation et simulation numérique du problème en vue. La méthode de couplage de type Euler Lagrange, une formulation Eulerienne à maillage fixe pour le fluide et une formulation Lagrangienne pour la structure, est utilisée pour la modélisation de l’écoulement sanguin dans les artères. Différentes méthodes de pénalisation pour le couplage fluide structure ont été testées. Pour une meilleure consistance, plusieurs types de maillages fluide et structure ont été analysés, and des pas de temps vérifiant la condition de stabilité. Comme très peu de résultats expérimentant dans le domaine biomécanique sont disponibles. Les résultats numériques sont comparés aux résultats théoriques publiés dans la littérature. Une bonne corrélation entre les résultats numériques et les données théorique a été signalée dans le manuscrit. Ces résultats sont publiés dans un Journal International de rang A. / The aim of this thesis is to investigate blood flow in arteries using Fluid-Structure Interaction (FSI) numerical approach. There are different approaches which could be used, Finite Element (FE) method for modelling artery wall and either FE or Smoothed Particles Hydrodynamic (SPH) method for blood modeling. In order to investigate the appropriate numerical method to simulate the biomedical problem. both SPH and Finite Element methods were applied. Both methods were first validated for applications where experimental data are available. In order to validate the SPH and FEM method used to simulate the fluid domain and arteries, we investigate a specific problem concerning membrane inflation. Finite Element method with shell formulation was used for the membrane made of rubber material, and SPH Particle method for the fluid. This application has been selected since experimental data are available. Mooney Rivlin constitutive material law for hyper elastic incompressible material is used for the membrane and compressible air for the fluid. For mesh sensitivity and consistency, different meshes for the membrane and different particle numbers for SPH have been investigated, Good agreement with experimental data were obtained in term of displacement and velocity of the membrane center. The application is published in an international Journal with Index Citation. The second formulation developed in the manuscript concerns the fluid structure coupling problem we are interested in. Coupled Eulerian Lagrangian (CEL) approach for the modelisation of Pulse Wave Velocity (PWV) in large arteries. Different penalty methods for the fluid structure coupling have been investigated for good accuracy and consistency of the results using different fluid and structure meshes and time step satisfying stability condition. Since experimental data are not available for these approaches, numerical results were compared to the theoretical theory and previous work published in the literature. This work is published in an International Journal with Index Citation.

L’objectif de ce travail est la prédiction numérique d’écoulements turbulents de fluides newtoniens et incompressibles à l’aide de modèles à viscosité turbulente. La présence des parois est la plupart du temps prise en compte par l’introduction des fonctions d’amortissement ou par des termes non-linéaires dans les équations du modèle, ce qui pénalise l’universalité du modèle. L’approche de relaxation elliptique permet d’éviter ces types de problèmes. Le modèle v² -f développé par Manceau et Hanjalic et basé sur la pondération elliptique a été utilisé dans cette thèse. En s’appuyant sur des procédures a priori et a posteriori, ce modèle a été comparé avec les modèles linéaire k-[epsilon] et non-linaire de Shih. Les tests a priori ont montré l’aptitude du modèle à relaxation elliptique à prédire les contraintes du tenseur de Reynolds. La comparaison des grandeurs moyennes et statistiques de l’écoulement, prédites avec les données de la DNS de la littérature pour une large gamme de nombre de Reynolds en canal, montre une bonne performance de ce modèle / The aim of this work is to predict numerically the turbulent flows of the newtonian and incompressible fluid using eddy viscosity models. The presence of walls is usually accounted for by introducing damping functions or non linear terms in the model equations. This approach does not favor its universality. The elliptic relaxation approach avoids these problems. The v² -f model developed by Manceau and Hanjalic and based on the elliptic blending was studied here. By using an a priori and a posteriori approaches, this model was compared with the linear k-[epsilon] and the non linear Shih models. The a priori tests demonstrate the capability of the elliptic relaxation model to predict the Reynolds stress tensor. The mean flow field and the turbulent statistics are compared to the existing DNS data of a turbulent channel flow for a wide range of Reynolds numbers. The model performance is shown to be satisfactory

La réduction de l’épaisseur des pièces injectées en polymères répond à des objectifs d’allègement, de productivité et de développement durable. Elle peut avoir, parce que mécanismes l’interface polymère/outillage deviennent prépondérants, des effets négatifs sur la qualité du produit fini et les moyens de fabrication à mettre en œuvre. L’optimisation du procédé nécessite d’une part la compréhension de ces mécanismes locaux, d’autre part la proposition de solutions techniques pertinentes aptes à faciliter l’écoulement. Un dispositif expérimental de rhéométrie capillaire en ligne a été développé, permettant d’accéder aux gradients de vitesse très élevés rencontrés dans le procédé de fabrication. Son utilisation, en combinaison avec celle de moyens rhéologiques plus classiques a permis de juger de la pertinence de diverses solutions envisagées pour faciliter l’injection d’un polyéthylène basse densité en parois minces (réduction de la masse moléculaire, addition d’agents glissants, ajout de nanocharges visant à augmenter la rigidité de la pièce finie). La caractérisation rhéologique isotherme a montré l’inefficacité des agents glissants considérés et les limites potentielles de la réduction de masse moléculaire. Elle a en revanche mis en évidence l’intérêt de l’utilisation de nanocharges en termes d’amélioration des propriétés mécaniques et d’aptitude à l’écoulement en parois minces. La validation sur pièces industrielles injectées (en régime anisotherme) conduit à nuancer quelque peu ces conclusions. Dans ce cas, l’utilisation d’une masse moléculaire faible nanorenforcée s’avère plus pertinente, notamment par un mécanisme de concentration de la charge dans la veine fluide / Reducing the wall thickness of the injection moulded polymeric parts meets requirements of lightness, productivity and sustainable development. This may induce poor quality of the final product and issues during processing. Those effects are due to the mechanisms taking place at the polymer/tool which prevail during the flow. Processing optimization requires a better understanding of these local mechanisms and effective technical solutions to facilitate the polymer flow. An original experimental device for online capillary rheometry has been developed, in which very high shear rates, characteristic of the injection moulding process, can be achieved. This device, combined with that data from more traditional rheological characterization allows to evaluate different solutions for improving the injection moulding of a low density polyethylene in thin walls (molecular weight reduction, slip agents, nanoparticles intended for increasing the rigidity of the final part). The isothermal rheological investigation have shown the inefficiency of the lubricating agents considered for this study and the limited potential of the reduction of polymer molecular weight. Nevertheless, it highlighted the interest of adding nanoparticles to the injected polymer to assure an improvement of the mechanical properties and for varying the flow property in thin walls. The validation of this method on industrial injected parts (anisothermal mode) offered some finer conclusions. In particular, the use of nanocomposites based on low molecular weight polymers seems more relevant, in particular due to a concentration mechanism of the nanoparticles in the fluid vein.

Les suies, particules omniprésentes dans l’atmosphère urbaine, sont d’une extrême complexité morphologique et chimique. Adsorbés à leur surface, les hydrocarbures aromatiques polycycliques (HAP) présentent une toxicité pulmonaire avérée. Durant leur transport atmosphérique, ils peuvent, par des processus d’oxydation et de photolyse, former des produits nitrés et oxygénés dont les impacts sanitaires sont peu connus. Nous proposons de caractériser le vieillissement atmosphérique de particules de suie et d’étudier l’impact de ce vieillissement sur les poumons. Ces travaux de thèse reposent sur la génération de particules de suie modèles ; la simulation en laboratoire d’un vieillissement par des expositions à l’ozone ; ainsi que l’étude des effets inflammatoires de particules vieillies ou non sur des cellules pulmonaires. Nous avons exposé, dans un premier temps, un aérosol de particules carbonées à l’ozone dans un réacteur à écoulement. Les influences du temps de contact, de la concentration initiale en ozone, de l’humidité relative, de la présence d’un rayonnement UV et de la taille des particules ont été étudiées. Le rôle de l’oxygène atomique dans le mécanisme de réaction a également été évalué. Nous nous sommes ensuite intéressés à la réactivité hétérogène entre des HAP adsorbés en surface de particules carbonées et l’ozone. Des cinétiques de dégradation de HAP particulaires ont ainsi été réalisées. Enfin, nous avons étudié la réponse inflammatoire de HAP seuls ou particulaires sur des cellules pulmonaires par le dosage d’interleukines pro-inflammatoires. L’effet du vieillissement des particules modèles, ainsi que le rôle du surfactant pulmonaire ont été évalués. / Soot particles, ubiquitous in the urban atmosphere, exhibit extreme chemical and morphological complexity. Adsorbed on their surface, polycyclic aromatic hydrocarbons (PAHs) have a proven pulmonary toxicity. During atmospheric transport, they may, through oxidation and photolysis process, form oxygenated and nitrated products about which health impacts are not well known. We propose to characterize the atmospheric aging of soot particles and to study the impact of aging on the lungs. This work deals with the generation of soot models particles; laboratory simulation of aging by exposures to ozone; and the study of the inflammatory effects of aged or not particles on lung cells. In a first step a carbonaceous aerosol has been exposed to ozone in a flow reactor. The influences of contact time, initial ozone concentration, relative humidity, UV radiations and particle size were studied. The role of the atomic oxygen in the reaction mechanism was also evaluated. Then the heterogeneous reactivity between PAHs adsorbed on carbonaceous particles and ozone has been investigated. Degradation kinetics of particulate PAHs were thus performed. Finally, we studied the inflammatory response of PAHs alone or in particulate phase on lung cells by assay of pro-inflammatory interleukins. The effect of aging models particles and the role of pulmonary surfactant were evaluated.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

La présente thèse met de l'avant le mélange turbulent d'un scalaire passif dans les écoulements cisaillés libres. L'écoulement étudié est un sillage qui se caractérise par un champ moyen scalaire antisymétrique transporté par un champ moyen dynamique symétrique. À notre connaissance ce sillage complexe est étudié pour la première fois expérimentalement. Deux plaques (2P) génératrices de sillages sont placées côte à côte, la première est faiblement refroidie et la deuxième est faiblement chauffée, ainsi la température est considérée comme un scalaire passif. Le nombre de *Reynolds* basé sur la longueur de la corde de la plaque $\it{l}$, et la vitesse externe est $\it{Re_l}$ ≈ 3.5 × (10)$^5$. La thermométrie à fil froid et l'anémométrie à fil chaud ont été mises en œuvre pour des mesures simultanées des quantités turbulentes. Les contraintes de *Reynolds*, les flux de chaleur, les moments d'ordre trois de vitesse ainsi que les moments mixtes vitesse-température, ont été mesurés et des courbes de régression ont été générées à partir de chaque profil. La contrainte de cisaillement de *Reynolds* et le flux de chaleur transversal qui sont calculés avec l'équation du mouvement moyen et l'équation d'enthalpie valident bien nos mesures expérimentales. La technique de la sonde à fils chauds croisés décalée en rotation de 45° a été utilisée pour mesurer le terme $\overline{vw^2}$. La comparaison avec les différents termes scalaires trouvés dans la littérature a requis l'utilisation de la variable autosimilaire *ξ*$_θ$ pour réconcilier les quantités et les rendre comparables. Aussi, le calcul des quantités scalaires avec les lois autosimilaires comme le flux de chaleur et les gradients longitudinaux, donne de meilleures résultats si la variable d'évolution radiale de température *ξ*$_θ$ est utilisée. Les bilans de l'énergie cinétique turbulente *k* et la demi-variance de la fluctuation du scalaire $\overline{θ^2}$/2 ont été générés dans le sillage 2P ainsi que dans le sillage d'une seule plaque (1P) légèrement chauffée pour servir de test et de validation de la technique expérimentale. La superposition des bilans de *k* et de $\overline{θ^2}$/2 du sillage 1P avec des données de la littérature, affichent le même ordre de grandeur ; les termes de production et de diffusion turbulente présentent les mêmes comportements. Par contre, le bilan de *k* du sillage 2P montre une dynamique complètement inverse par rapport au sillage 1P dans la zone −1 ≤ *ξ* ≤ 1, alors qu'à l'extérieur de cette zone, les différents termes du bilan ont les mêmes allures que dans le cas du sillage 1P. Le bilan de $\overline{θ^2}$/2 dans le sillage 2P, affiche un ordre de grandeur deux fois supérieur à celui du sillage 1P dans la zone −1 ≤ *ξ*$_θ$ ≤ 1 à cause de la forte production qui est générée par le profil moyen antisymétrique du scalaire. En effet, pour le sillage 2P, le terme de production du scalaire est prépondérant autour de l'axe, alors qu'il est pratiquement nul pour le sillage 1P. En ce qui concerne le terme de diffusion turbulente, il agit comme un terme puits important au voisinage de l'axe pour le sillage 2P alors qu'il agit comme un terme source dans la même zone pour le sillage 1P. Les mécanismes de mélange dans cette zone du sillage 2P sont par conséquent significativement différents de ceux observés pour le sillage 1P. Mais la partie extérieure du sillage semble évoluer de la même façon que la partie extérieure du sillage 1P. Le mélange du scalaire dans le sillage 1P est dominé par les structures de la vitesse longitudinale, alors que dans le sillage 2P ce sont les structures de la vitesse transversale qui dominent le mélange dans la zone −1 ≤ *ξ*$_θ$ ≤ 1. Enfin, les coefficients de diffusion de la chaleur et de la quantité de mouvement ont été calculés et comparés avec les données de la littérature, et comme attendu, les valeurs du sillage 1P se comparent bien avec les données publiées. Par ailleurs, il a été observé, en particulier, que la diffusivité turbulente dans le sillage 2P affiche une asymptote à la position *ξ*$_θ$ = 1. Ce comportement est inexistant dans le sillage 1P et aussi dans un sillage 2P avec deux plaques chauffées. Il apparait donc que ce sillage 2P avec le profil de température moyenne utilisé mérite plus d'attention. Le rapport des échelles de temps thermique et mécanique a été calculé et une courbe oscillante autour de l'unité est observée pour le sillage scalaire 2P. Ce résultat jumelé à l'analyse du plan des invariants du tenseur de *Reynolds* pourrait faire de cet écoulement du moins en ce qui concerne la région interne, un bon candidat pour la modélisation de la turbulence homogène et isotrope. La base de données collectée durant cette thèse fournira plus de résultats, en particulier les bilans des contraintes de *Reynolds* et les bilans des flux de chaleur pour les deux sillages. / The present thesis highlights the turbulent mixing of a passive scalar in free sheared flows. The flow studied is a wake that is characterized by an antisymmetric scalar mean field carried by a symmetric dynamic mean field. To our knowledge, this complex wake is being studied experimentally for the first time. Two wake generating plates (2P) are placed side by side, the first is slightly cooled and the second is slightly heated, so the temperature is considered as a passive scalar. The number of *Reynolds* based on the length of the plate chord $\it{l}$, and the external velocity is $\it{Re_l}$ ≈ 3.5 × (10)$^5$. Cold wire thermometry and hot wire anemometry were implemented for simultaneous measurements of turbulent quantities. *Reynolds* stresses, heat fluxes, third-order velocity moments as well as mixed velocity-temperature moments were measured, and regression curves were generated from each profile. The *Reynolds* shear stress and the transverse heat flux calculated with the equation of mean motion and the enthalpy equation well validate our experimental measurements. The 45° rotated cross-wire hot wire probe technique was used to measure the term $\overline{vw^2}$. Comparison with the various scalar terms found in the literature required the use of the self-similar variable *ξ*$_θ$ to reconcile the quantities and make them comparable. Also, calculating scalar quantities with self-similar laws such as heat flux and longitudinal gradients gives better results if the self-similar variable *ξ*$_θ$ is used. The turbulent kinetic energy *k* budgets and the half-variance of the scalar fluctuation $\overline{θ^2}$/2 budget were generated in the 2P wake as well as in the wake of a single plate (1P) slightly heated to serve as a test and validation of the experimental technique. The superposition of the *k* and $\overline{θ^2}$/2 budgets of the 1P wake with data from the literature show the same order of magnitude; the terms of production and turbulent diffusion exhibit the same behaviors. On the other hand, the *k* budget of the 2P wake shows a completely opposite dynamic with respect to the 1P wake in −1 ≤ *ξ* ≤ 1 area, while outside this zone, the terms of the budget look the same as in the case of the 1P wake. The budget of $\overline{θ^2}$/2 in the 2P wake displays an order of magnitude twice that of the 1P wake in the zone −1 ≤ *ξ*$_θ$ ≤ 1 because of the high production which is generated by the antisymmetric profile of the average scalar. Indeed, for the 2P wake, the scalar production term is preponderant around the axis, while it is practically zero for the 1P wake. Regarding the turbulent diffusion term, it acts as an important sink term in the vicinity of the axis for the 2P wake while it acts as a source term in the same area for the 1P wake. The mixing mechanisms in this zone of the 2P wake are therefore significantly different from those observed for the 1P wake. But the outer part of the wake seems to evolve in the same way as the outer part of the 1P wake. The scalar mixing in the 1P wake is dominated by the longitudinal velocity scales, whereas in the 2P wake it is the transverse velocity scales that dominate the mixing in the zone −1 ≤ *ξ*$_θ$ ≤ 1. Finally, the heat and momentum diffusion coefficients were calculated and compared with literature data, and as expected, the 1P wake values compare well with published data. Furthermore, it was observed in particular that the turbulent diffusivity in the 2P wake displays an asymptote at the position *ξ*$_θ$ = 1. This behavior is non-existent in the 1P wake and also in a 2P wake with two heated plates. It therefore appears that this 2P wake with the average temperature profile used deserves more attention. The ratio of thermal and mechanical time scales was calculated and an oscillating curve around unity is observed for the 2P scalar wake. This result, combined with the analysis of the plane of the invariants of the *Reynolds* stress tensor, could make this flow, at least as far as the inner region is concerned, a good candidate for the modeling of homogeneous and isotropic turbulence. The database collected during this thesis will provide more results, in particular the *Reynolds* stress balances and the heat flux balances for the two wakes.

Ce travail de thèse porte sur l’instabilité de Faraday qui résulte de l’oscillation verticale d’une couche de liquide. L’objectif est d’apporter une vérification expérimentale précise des résultats issus d’une analyse de stabilité linéaire. À basse fréquence, la longueur d'onde caractéristique de cette instabilité est de l'ordre des dimensions de la cellule expérimentale. Théoriquement les longueurs d'onde associées aux modes disponibles sont discrétisées selon le confinement latéral. Mathématiquement, la condition qui facilite les calculs est une condition limite pariétale sans contrainte, ce qui est irréaliste en raison des couches limites visqueuses et de l’hystérésis capillaire au point de contact du fluide avec les parois. Ici, nous présentons l'étude expérimentale d'un système de basse fréquence conçu pour se rapprocher d’une condition pariétale sans contrainte, par un choix judicieux de liquides. Le système a permis d’obtenir expérimentalement des seuils d’apparition de l’instabilité très proches des prédictions théoriques basées sur une analyse de stabilité linéaire. Les différences observées sont associées au ménisque formé par le liquide sur la paroi ainsi qu’au film liquide se formant à la paroi lors de l'excitation des ondes. Les écarts avec la théorie dû à ces imperfections semblent être fonction de la structure du mode et la périodicité des ondes. De plus, les ondes générées par le ménisque conduisent à des seuils inférieurs aux prévisions. Ceci indique une résonance mixte entre l'instabilité paramétrique et les oscillations forcées du ménisque où chaque contribution peut être contrôlée par un forçage caractérisé par deux fréquences. / This work is an inquiry into the current understanding of experimental single-mode Faraday waves, from the perspective of linear stability theory. At low frequency, the characteristic wavelength is of the order of the cell dimensions, and theoretically the wavelengths associated with the modes available for excitation are discretized by the sidewall boundary condition. Mathematically, the condition which eases manipulation is a stress-free sidewall condition, which is unrealistic due to viscous boundary layers and capillary hysteresis. We report the experimental investigation of a low frequency system designed to approximate this condition, by careful choice of the fluid liquids and thereby increasing the relative importance of the interior damping. The design has permitted close connection between the experimental thresholds and the predictions, contrasting the traditional approach of phenomenologically measuring the mode dissipation and treatment as a single degree of freedom system. Non-ideal behavior of the sidewall persists, evident by the contact meniscus and the sidewall film which is formed during wave excitation. The presence of these non-idealities is translated to trends in the deviation between the observed and predicted thresholds, depending upon the mode structure and wave periodicity. The presence of harmonically excited waves from the meniscus seems to result in observation of amplitudes consistently lower than the predictions, indicating a mixed resonance between the parametric instability and the forced meniscus oscillations. Furthermore, it is seen this interaction can be adjusted by excitation with two commensurate parametric frequencies.

La thèse concerne la modélisation de deux écoulements issus de problèmes géophysiques. Un premier problème a trait à la présence de tourbillons longitudinaux apparaissant dans les cours d'eau. L'origine de ces structures reste indéterminée. On aborde ce problème par une simulation numérique d'un écoulement turbulent cisaillé dans un canal ouvert utilisant un code pseudo-spectral. On tente de voir si un cisaillement de surface même faible est capable ou non de comprendre ces observations. Le second problème est lié à la resuspension des sédiments sur les bords d'une mer ou d'un lac par des ondes internes. Ces ondes existent à cause de la stratification en densité de la colonne d'eau. Lorsque elles s'approchent de la côte, elles se déstructurent générant un cisaillement sur le fond capable de resuspendre du sédiment. On a envisagé à nouveau cette étude par le biais de la simulation numérique directe. On examine comment l'onde interne se brise sur les bords en fonction de (a) la nature de la stratification, (b) la forme de la topographie du fond et (c) l'amplitude des ondes. On calcule dans chaque cas le cisaillement sur le fond. On en déduit le flux et le transport de sédiment dans toute la colonne d'eau. / This thesis studies the modeling of two problems that take their origin from a geophysical context. A first problem is related to the presence of longitudinal vortices which have been measured in rivers. The origin of these structures remains unknown. We address this problem by numerical simulations of a sheared turbulent flows in an open channel using a pseudo-spectral code. We try to determine if the presence of an imposed shear at the surface coupled with a pressure gradient is capable or not to explain these observations. The second study focuses on the sediment resuspension on shores of seas or lakes by the action of internal waves. These waves exist because of density stratification of the water column. When waves approach the shore, their patterns evolve generating a shear on the bottom capable to resuspend sediment. By direct numerical simulations, we analyse how internal waves breaking changes according to (a) stratification, (b) bottom topography and (c) wave amplitude. We compute for each case the shear exerted on the bottom, the sediment flux and transport throughout the water column.

Un réacteur à écoulement stationnaire muni d’une tête d’injection mobile a été développé au Département Chimie et Environnement de l’École des Mines de Douai afin de mieux comprendre les processus de formation d’aérosols organiques secondaires (AOS) suite à l’oxydation en phase gazeuse de certaines espèces particulièrement réactives et/ou abondantes dans l’atmosphère. Il permet d'observer la formation des produits dans les premières étapes de la réaction à pression et température ambiantes, et en conditions stationnaires, contrairement aux chambres de simulation atmosphérique. Une première partie de ce travail a consisté à caractériser l’écoulement dans le réacteur, et à développer et optimiser les méthodes utilisées pour l’analyse de la phase gazeuse. Dans un second temps, le système a été validé du point de vue cinétique par l’étude de l’ozonolyse de pentènes méthylés et de l’α-pinène. Les coefficients de vitesse ont été déterminés en excès d’alcène et ont montré un bon accord avec les valeurs de la littérature. Pour l’un des alcènes, les deux voies réactionnelles attendues de dégradation ont été identifiées et quantifiées (en bon accord avec la seule étude existante), et un intermédiaire non oxydé, provenant du réarrangement d’un biradical Criegee, a pu être mis en évidence aux faibles temps de réaction. Enfin, la formation d’AOS par ozonolyse de l’α-pinène a été investiguée, avec caractérisation simultanée des concentrations en réactifs et produits en phase gazeuse (par analyseur d’ozone et TD/GC/FID-MS), de la distribution granulométrique des particules formées (par SMPS et HR-ToF-AMS), et de la composition chimique moyenne des particules (HR-ToF-AMS). / A steady flow reactor equipped with a mobile injection head was developed at the Chemistry and Environment Department, Ecole des Mines de Douai, to better understand the formation processes of secondary organic aerosols (SOA) following the gas phase oxidation of certain species particularly reactive and/or abundant in the atmosphere. In contrast to atmospheric simulation chambers the flow reactor allows to observe in stationary conditions the product formation in the early stages of the reaction at ambient pressure and temperature.The foremost part of this work was to characterize the reactor flow, and to develop and optimize the methods used for the gas phase analysis. In a second stage a kinetic study of the ozonolysis of methylated pentenes and α-pinene was performed to validate the system. The rate coefficients were determined in pseudo first order conditions with excess alkene and showed a good agreement with literature values. The expected oxidation pathways were identified and quantified for one alkene (in good agreement with the only study available), and a nonoxidized intermediate coming from the rearrangement of a Criegee biradical, has been found for short reaction times. Finally, the SOA formation from α-pinene ozonolysis was investigated with simultaneous characterization of the gas phase concentrations of reactants and products (by ozone analyzer and TD/GC/FID-MS), of the particle size distribution (by SMPS and HR-ToF-AMS), and of the average chemical composition of particles (HR-ToF-AMS).