Cette thèse a pour but de contribuer à l’étude des aquifères karstiques en (1) identifiant les facteurs influençants le stockage et l’écoulement de l’eau souterraine et en (2) localisant et estimant de manière quantitative les ressources en eau au sein d’une série carbonatée. Cette étude a été réalisée sur l’aquifère de Dardennes, près de Toulon, en France. L’eau de l’aquifère est utilisée pour l’alimentation en eau potable et alimente le fleuve Las, susceptible de provoquer de fortes crues. Le bassin versant de l’aquifère est estimé à 70 km². Des mesures de porosité et de géochimie isotopique sur les roches carbonatées sont mises en relation dans le but de localiser et d’expliquer les vides créés au sein de la série sédimentaire. Des coupes équilibrées et des analyses structurales ont été réalisées pour évaluer la géométrie de l’aquifère. L’acquisition de données en continu (CTD) et de données hydrochimiques permet d’étudier le dynamisme et le fonctionnement de l'aquifère. Une approche de modélisation pluie-débit avec la plateforme de modélisation KarstMod, basée sur les chroniques de précipitations et de débits sur trois années a été utilisée afin d’estimer la part de l’écoulement disponible pour l’eau potable et celle générant des crues éclair. Une déconvolution du débit des sources a été effectuée à l’échelle annuelle et à l’échelle des crues. Les résultats du modèle ont été comparés avec la gestion actuelle de l’aquifère, montrant une exploitation actuelle équilibrée, et ouvrant des perspectives pour la mise en place d’une gestion active, combinant une exploitation plus importante de l’eau, un écrêtage des crues et un respect de l’écologie du cours d’eau en aval. / The aim of this PhD is to contribute to the karst aquifers understanding by (1) identifying the factors influencing the groundwater storage and flow and by (2) locating and quantitatively estimating water resources within the carbonate sedimentary series. This study was conducted on the Dardennes aquifer near Toulon, France. The aquifer water is used for the supply of drinking water and feeds the Las River, which can cause important floods. The recharge area of the aquifer is estimated at 70 km². Measurements of porosity and isotopic geochemistry on carbonate rocks are related in order to locate and explain the voids created within the sedimentary series. Balanced cross-sections and structural analyses were performed to evaluate the geometry of the aquifer at depth. Continuous data acquisition (specific conductivity, temperature and water depth-pressure) and hydrochemical data (major ions, stables isotopes of water and carbon-13) allow to study the aquifer dynamism and functioning. A daily rainfall-discharge modeling approach was used with the KarstMod platform to estimate the flow part available for drinking water and the part that generates karst flash-floods. A hydrograph separation of springs discharge was carried out at annual and floods scales. Model results have been compared with the current groundwater management, showing a balanced current exploitation. In new perspective, it would be interesting to establish active management, combining a greater groundwater exploitation, floods mitigation and stream ecology respect.

Malgré leur fort développement en assainissement domestique urbain, les bioréacteurs à membranes (BaM) pâtissent de phénomènes de colmatage, induisant des coûts énergétiques et de maintenance importants. L’aération séquencée des membranes par des grosses bulles est l’une des stratégies pour limiter le colmatage ; son impact a fait l’objet de plusieurs études mais reste globalement mal compris, notamment du fait de la complexité de ces systèmes, multiphasiques et opaques. L’hydrodynamique des réacteurs reste mal caractérisée en présence de boues. Pour apporter des éléments de compréhension aux mécanismes de limitation du colmatage par injection d’air, un pilote de filtration membranaire semi-industriel (2 m3, 3 sous modules fibres creuses Puron®) a été conçu, dimensionné et installé sur l’unité de traitement des jus (TDJ) de la station d’épuration Seine Aval (SIAAP), afin d’être alimenté en boues biologiques dans des conditions réelles de fonctionnement. Le suivi des paramètres opératoires du pilote et de ses performances de filtration sur une période de 5 mois avait un double objectif : (i) mieux caractériser la dispersion du gaz pour différentes conditions de fonctionnement (paramètres de l’aération, concentration en boues de l’alimentation), (ii) hiérarchiser les facteurs qui limitent le colmatage des membranes. Il s’agit in fine de proposer des stratégies d’aération adaptées et efficaces pour limiter le colmatage. Afin de caractériser la dispersion du gaz dans le réacteur, une méthodologie innovante basée sur la tomographie de résistivité électrique (ERT) a été adaptée au pilote. Les conditions d’utilisation de l’ERT (nombre d’électrodes de mesure, séquence de quadripôles) ont été sélectionnées à travers une étude numérique, de même que les paramètres d’inversion nécessaires pour reconstituer la cartographie des résistivités à partir des mesures expérimentales. Cette étude numérique poussée, réalisée sous COMSOL, a permis de conclure à l’intérêt de la méthode pour représenter la distribution des phases dans la géométrie considérée. L’ERT a donc été appliquée au pilote alimenté en boues, pour différentes conditions d’aération. Le jeu de données de filtration a par ailleurs été analysé par logique floue, à l’aide du logiciel FisPro. Les arbres de décision obtenus, en analysant les résultats de manière globale et en les regroupant par conditions opératoires similaires, ont mis en évidence l’impact prépondérant des variables suivantes sur la dérive de perméabilité observée (comprise entre - 9 et 2 LMH/bar) : la différence de DCO entre le surnageant des boues et le perméat (DDCO) traduisant une phase colloïdale complexe, et la concentration en matières en suspension (MES), ayant toutes deux un impact négatif sur les performances de filtration. Une augmentation du débit d’air conduirait à une limitation de la dérive de perméabilité, sauf lorsque la variable DCO est élevée (> à 500 mg/L), cette hypothèse restant cependant à vérifier sur une base de données plus conséquente. Le modèle ainsi obtenu par logique floue permet de mieux simuler les évolutions de perméabilité que les modèles obtenus par régression linéaire multivariée (erreurs de 0,61 et de 0,70 respectivement), et ce malgré une incertitude relative importante sur la mesure de perméabilité (jusqu’à 16 %). Ces résultats sont cohérents avec la dispersion du gaz observée par ERT : son homogénéité dépend de la concentration en MES et du débit d’air injecté. A forte concentration en MES (6 – 10 g/L), des zones préférentielles de passage des bulles ont été observées, en particulier à faible débit d’air, expliquant ainsi un colmatage plus important. L’utilisation nouvelle dans ce contexte de ces techniques, ERT et logique floue, donne des résultats qui confortent l’intérêt d’adapter l’aération (débit, séquençage) aux caractéristiques des boues notamment leurs concentrations, et qui permettent d’envisager des stratégies de contrôle de ces paramètres / Membrane bioreactors (MBRs) are widely used in the wastewater treatment sector. However, membrane fouling mitigation remains challenging, and leads to important maintenance and energy costs. Several strategies have been developed industrially to enhance MBR productivity, including coarse bubble sequenced aeration for fouling mitigation. The way such aeration participates in sludge hydrodynamic patterns is an important research topic. However, the methods currently used for hydrodynamic characterisation suffer from several drawbacks, mainly due to the system’s complexity (three phases with opaque deformable solids). More research is needed to characterize hydrodynamics in MBRs filled with activated sludge. A semi-industrial reactor (2.0 m3) was designed and equipped with three hollow fiber membrane modules (KMS Puron) to contribute to fill this knowledge gap. The reactor was continuously fed with activated sludge from one of the reactors of the Seine Aval (SIAAP) wastewater treatment plant. Operating parameters, activated sludge properties as well as process performances were monitored for five months with two objectives: (i) to characterize gas dispersion for different operating conditions (aeration strategy, activated sludge concentrations), (ii) to rank activated sludge properties and operating conditions, according to their impact on fouling mitigation. The ultimate goal is to propose adapted aeration strategies allowing sustainable filtration performance. In order to characterize the gas dispersion in a complex geometry such as MBR membrane tanks, an innovative approach based on electrical resistivity tomography was adapted to the pilot conditions. A numerical approach was used to define the experimental design in terms of electrode positions, quadripole sequences but also inversion parameters, used to reconstruct resistivity maps from experimental datasets. This numerical study, performed on COMSOL, demonstrated the ability of the ERT method to observe different patterns in the membrane zone. ERT was therefore applied to characterise bubble dispersion in the semi-industrial membrane bioreactor filled with activated sludge and operated using different aeration conditions. The dataset obtained from filtration monitoring was analysed using fuzzy logic. The resulting fuzzy decision trees, constructed from the whole dataset or from subsets clustering similar operating conditions, pointed out the impact on the permeability evolution of two main factors: the COD difference between activated sludge supernatant and permeate (COD), standing for complex colloidal phase, and Mixed Liquor Suspended Solid (MLSS) concentration. Both variables emphasizing fouling at high values. Also, an air flow rate augmentation seemed to mitigate fouling, except at high COD concentrations (> 500 mg/L), where its impact was reversed. This last observation needs further investigation, on a larger dataset. The model obtained from fuzzy logic allows a better simulation of permeability evolutions compared to multivariate linear regression (errors of 0.61 and 0.70, respectively), despite a relatively high measurement uncertainty of permeability (up to 16 %). Those results are consistent with the gas dispersion observed with ERT measurements: homogeneity of gas dispersion depends on MLSS concentration and air flow rate. At high concentrations of MLSS (6 g/L and 10 g/L), a high air flow rate is required for the bubbles to be distributed on the whole membrane zone and preferential flow paths have been observed, especially at low air flow rate. This gas dispersion heterogeneity explains worsened filtration performance. In this context, the new use of these technics, ERT and fuzzy logic, provided results that reinforce the interest to link aeration parameters to sludge properties, their MLSS concentration in particular, leading to consider control strategies for these operating parameters.

Ce travail de thèse s’inscrit dans une démarche d’approfondissement des connaissances du fonctionnement des aquifères karstiques. Pour cela, le site pilote des Sources du Toulon, siège d’une alimentation multiréservoirs, a été choisi et offre une fenêtre d’observation privilégiée sur les relations hydrogéologiques entre les aquifères de la marge Nord du Bassin aquitain. Elles sont utilisées depuis 1832 pour l’alimentation en eau potable de la ville de Périgueux, et constituent actuellement son unique ressource.L’objectif de ce travail a été de tester une approche multi-traceurs et un suivi haute-résolution pour caractériser le fonctionnement et la vulnérabilité d’un aquifère karstique multicouche.L’utilisation d’une approche couplée hydrodynamique et hydrochimie a mis en évidence l’alimentation des Sources du Toulon par une ressource profonde et captive en complément d’un aquifère libre plus superficiel. Ces informations ont été intégrées dans le contexte hydrogéologique global de la zone d’étude pour proposer une nouvelle délimitation du bassin d’alimentation des Sources du Toulon.L’utilisation d’un suivi haute-résolution a permis de préciser le fonctionnement hydrogéologique du système karstique étudié. Il en ressort des fonctionnements différents suivants les crues, avec des transferts de masse pouvant s’étendre de quelques jours à quelques semaines.L’apport des isotopes des nitrates (δ15N-NO3- ; δ18O-NO3-) a permis d’identifier deux origines des nitrates mesurés dans les eaux de surface et souterraines : les fertilisants minéraux de synthèse utilisés pour l’agriculture, et les rejets d’assainissement individuels.L’évaluation temporelle de la vulnérabilité spécifique des Sources du Toulon à certains marqueurs de contamination tels que les particules, les éléments bactériologiques et les nitrates a permis de distinguer plusieurs masses d’eau superficielles et subsuperficielles participant à l’alimentation pendant les crues, et jouant un rôle différent dans le transfert des contaminants. / This work is included in an approach for a better knowledge of karst aquifers. For this, the Toulon Springs pilot site was chosen, and provides the opportunity to study relationships between multilayered karst aquifers of the northern edge of the Aquitaine sedimentary basin. Toulon Springs are major regional springs and are located in Périgueux (Dordogne County, France). They have been supplying water to the metropolitan area of Périgueux since 1832.This work aim to test an innovative multitracer approach coupled with a physicochemical high-resolution auto-monitoring to characterize functioning and vulnerability of a multilayered karst aquifer.The use of a coupled hydrodynamic and hydrochemical approach has highlighted Toulon Springs supply by a deep and captive aquifer, in addition to a more subsuperficial aquifer. This information has been incorporated into the global hydrogeological situation of the study area to propose a new delineation of Toulon Springs hydrogeological cathment.The use of a high-resolution monitoring enables to specify the hydrogeological functioning of the studied karst system. The analysis performed on several flood events, identify that mass transfer can range from a few days to a few weeks.The use of nitrate isotope (δ15N-NO3- ; δ18O-NO3-) enables to identify two main nitrate sources in both surface water and groundwater: (i) inorganic fertilizer; (ii) sewage from individual house.The last part of this work was to analyse Toulon Springs temporal vulnerability during floods to marker of contamination such as particles, dissolved organic carbon, nitrate, and bacteria. Hence, several water type from surface runoff, unsaturated zone and saturated zone were identify as responsible of contaminant transfer: (i) water from saturated zone is responsible of vulnerability to mineral particles ; (ii) water from unsaturated zone is responsible of vulnerability to nitrare ; (iii) water from surface runoff is responsible of vulnerability to bacteria contamination.

Cette thèse est centrée sur l’analyse du rôle de l’adhésion cellule-substrat et des facteurs de détectionde Quorum sur la migration amibienne de Dictyostelium . Des outils pour automatiser les enregistrementsde vidéoomicroscopie et l’analyse d’image ont été développés afin de travailler avec de trèsgrands échantillons de cellules et de quantifier la migration cellulaire. Un dispositif microfluidique dedétachement cellulaire sous flux hydrodynamique combiné une platine motorisée a permis une étude statistique de l’adhésion mais aussi de la dynamique de détachement. L’analyse de la migration de Dictyostelium en milieu non nutritif met en évidence le rôle de la densité sur la différentiation des celluleset leur capacité de migration. Nous observons la présence d’une vitesse maximale de migrationaprès 6h de carence. Nous montrons que l’adhésion sur verre est deux fois plus faible en milieu carenc´equ’en milieu nutritif. Les exp´eriences de migration en milieu nutritif ont révélé la présence d’un facteur de détection de densité sécrété par les cellules et régulant leur migration aléatoire. Le coefficient de diffusion, la persistance du mouvement et la morphologie des cellules varient en fonction de la concentrationde ce facteur. Ce facteur ne modifie pas l’adhésion cellulaire mais uniquement la dynamique de d´etachement. Enfin, le protocole d’analyse développé nous a permis de faire une étude comparative de la migration en faisant varier d’autres paramètres tel que la surface ou la composition chimique du milieu expérimental. Ce travail se conclue en exposant le possible rôle de l’adhésion sur la migrationchez Dictyostelium en milieu nutritif. / This thesis focuses on the analysis of the role of adhesion between substrate and cell and factors of Quorum sensing on the migration of Dictyostelium amoeba. Tools to automate the recordings of videomicroscopy and image analysis have been developed to work with very large samples of cells and toquantify cell migration. A microfluidic device for cell detachment in hydrodynamic flow combined witha motorized stage has allowed a statistical study of adhesion but also the dynamics of detachment. The analysis of the migration of Dictyostelium in non nutritive medium highlights the role of density on celldifferentiation and migration capacity. We observe the presence of a maximum speed of migration after6 hours of starvation. We show that the adhesion to glass is twice as low in deprivation buffer as inthe nutrient medium. The experiences of migration in growth medium revealed the presence of a factorof detection of density secreted by the cells and regulating their random migration. The diffusion coefficient, the persistence of the movement and morphology of cells vary depending on the concentrationof this factor. This factor does not affect cell adhesion but only the dynamics of detachment. Finally, the testing protocol developed allowed us to make a comparative study of migration by varying otherparameters such as surface or the chemical composition of experimental medium. This work concludesby outlining the possible role of adhesion to the migration of Dictyostelium in nutrient medium.

On s'intéresse dans cette thèse à la modélisation et à la simulation multidimensionnelle de l'hydrodynamique fluviale, notamment près des estuaires. Le modèle physique de référence est le modèle 3D, mais au vu de son important coût de calcul, il est intéressant de disposer de modèles plus simples en 1D, 2D ou 2.5D, que l'on peut utiliser dans des zones adéquates du fleuve, en fonction de sa bathymétrie. Ainsi, à partir du modèle 3D basé sur les équations instationnaires et incompressibles de Navier-Stokes, des modèles plus simples sont dérivés par projection par formulations faibles du problème 3D. On obtient ainsi un modèle en 1D, écrit sur la courbe médiane de la surface libre du fleuve, ainsi que deux modèles en 2D, le 2D-vertical écrit sur la surface longitudinale médiane du fleuve et le 2D-horizontal écrit sur la surface libre. Enfin on définit un modèle en quasi-3D, le modèle 2.5D, écrit dans la somme des espaces 2D-vertical et 2D-horizontal. Tous ces modèles prennent en compte la géométrie du fleuve et fournissent une vitesse tridimensionnelle ainsi que la pression, qui n'est pas supposée hydrostatique mais qui est une inconnue entière du problème. En outre, on définit et justifie un estimateur de modèles entre le modèle 3D et chacune de ses approximations en 1D, 2D et 2.5D. Cet estimateur calcule l'erreur entre le modèle 3D et son approximation, et donne ainsi une indication sur la qualité des résultats obtenus à partir des modèles 1D, 2D ou 2.5D, dans leurs zones respectives de calcul. Tous ces modèles hydrodynamiques sont implémentés dans des codes d'éléments finis, écrits en C++. Enfin, ils sont couplés numériquement à l'aide de l'estimateur de modèles.

Ce travail est consacré à la modélisation des stries dans une couche limite laminaire et à l'étude de leur influence sur la transition. Dans une première partie, on présente une modélisation de ces structures longitudinales en utilisant les équations de Navier-Stokes paraboliques. La résolution numérique est basée sur des schémas de projection, de type Chorin-Temam, associés à des méthodes de discétisation spectrales. Ce modèle parabolique est appliqué au problème de réceptivité de la couche limite laminaire vis-à-vis de différentes perturbations extérieures. La dernière étape du processus de transition vers la turbulence est abordée par une étude de stabilité linéaire ainsi qu'une simulation numérique directe. Enfin on propose un modèle pragmatique pour reproduire la dynamique des stries et estimer le seuil de transition 'bypass'. Sur l'ensemble, on a cherché à comparer les résultats à différents travaux expérimentaux existants.

Cette thèse est une étude théorique et expérimentale de l'écoulement d'un fluide dans un cylindre en précession. Elle s'inscrit dans le domaine des écoulements tournants et trouve de nombreuses applications en aéronautique (précession des objets volants avec carburant liquide) et géophysique (influence de la précession terrestre sur la dynamique du noyau externe). Pour des nombres de Reynolds modérés, l'écoulement d'un fluide en précession est une rotation solide à laquelle se superposent des ondes appelées modes de Kelvin. Ces ondes deviennent résonnantes lorsque leur fréquence propre est égale à la fréquence de précession. Afin de prédire l'amplitude d'un mode de Kelvin résonnant, une théorie visqueuse et faiblement non-linéaire a été développée puis confirmée expérimentalement. Lorsque le nombre de Reynolds augmente, l'écoulement de base devient instable. Des mesures PIV ont montré que cette instabilité est due à une interaction entre le mode forcé par la précession et deux modes de Kelvin libres. Une analyse linéaire de stabilité basée sur un mécanisme de résonance triadique entre modes de Kelvin permet de prédire correctement le seuil et le taux de croissance de l'instabilité. Cette instabilité est saturée par un mode géostrophique créé par l'interaction non-linéaire et visqueuse des modes de Kelvin libres avec eux même. En accord avec les expériences, une théorie faiblement non-linéaire montre que les amplitudes des modes de Kelvin sont stationnaires ou oscillantes (cycles intermittents) en fonction du nombre de Reynolds. Finalement, le point fixe des équations non-linéaires (système dynamique à quatre degrés de liberté) prédit correctement l'écoulement moyen à l'intérieur du cylindre, même pour des très grands nombres de Reynolds (i.e. Re>50000).

Ce travail, réalisé dans le cadre d'une convention de thèse Cifre avec Thales Underwater Systems, et cofinancé par DCNS, vise à améliorer la prédiction du bruit propre hydrodynamique d'une antenne sonar, dû aux fluctuations de pression dans la couche limite qui se développe le long du dôme. Les estimations de bruit reposent en général sur des modèles semi-empiriques qui ne tiennent que partiellement compte de la flexibilité du dôme. La présente étude reconsidère le problème du bruit hydrodynamique dans le cas-type simplifié d'une couche limite instable le long d'une plaque plane élastique. La première partie du travail est consacrée à la simulation numérique directe d'un écoulement de couche limite bidimensionnel, caractérisé par un nombre de Reynolds hautement supercritique, le long d'une paroi élastique encastrée. La résolution des équations de Navier-Stokes incompressibles repose sur un changement de variables évolutif au cours du temps, associé à une discrétisation spatiale mixte différences finies – collocation Chebyshev. Une méthode de pas fractionnaire permet d'assurer un couplage fort entre le modèle de paroi élastique et le système fluide. Un forçage en vitesse est injecté dans ce dernier, à des fréquences instables et des amplitudes suffisamment élevées, afin de donner lieu à des instabilités convectives et saturées non linéairement, qui interagissent avec les mouvements de la paroi. Cette dernière vibre autour d'un état déformé initial résultant du couplage avec l'écoulement non perturbé. Des parois de longueurs et de matériaux variés ont été testées en vue de caractériser les vibrations, en termes de niveaux et de structures spatiales, ainsi que leur rétroaction sur les instabilités de l'écoulement, en fonction des valeurs des fréquences propres de la paroi. Dans les divers cas étudiés, on observe que le spectre de pression pariétale est enrichi par des composantes modales, à très bas nombres d'ondes et relativement hautes fréquences, qui peuvent donner lieu à un rayonnement plus important. Un post-traitement est mis en oeuvre afin d'estimer le bruit généré par les fluctuations de vitesse dans la couche limite. La pression rayonnée dans l'écoulement uniforme est évaluée à partir des données acquises au cours des simulations, dans le cadre de l'analogie de Lighthill. Le calcul, effectué dans le domaine spectral, repose sur l'utilisation d'une fonction de Green dont l'expression tient compte de la présence de la paroi souple. Dans cette approche, les vibrations de la paroi induisent à certaines fréquences, en particulier aux fréquences propres de la paroi, des niveaux de pression rayonnée supérieurs à ceux du cas rigide, en favorisant la présence de structures spatiales plus larges. D'autres modèles, élaborés à partir de la même formulation, mais tenant compte de la compressibilité de la couche limite dans le calcul des contributions vibratoires, mettent en évidence une nette augmentation des niveaux de pression rayonnée sur paroi souple dans une large gamme de fréquences, de part et d'autres des fréquences propres, lorsque ces dernières sont distinctes des fréquences de forçage, ainsi que l'apport de la prise en compte d'un couplage fort par rapport aux approches classiques basées sur une hypothèse de couplage faible. Enfin, une étude expérimentale en tunnel hydrodynamique est mise en place en vue de mesurer les vibrations d'une plaque élastique, encastrée dans un support rigide, et soumise à un écoulement transitionnel ou turbulent, ainsi que le bruit rayonné dans la cavité sous-jacente à la plaque, instrumentée avec un hydrophone. Les résultats de la campagne de mesures, qui repose sur les techniques de vibrométrie laser et de vélocimétrie laser Doppler, sont détaillés et analysés.

Ce mémoire a pour objet le calcul des collisions de gouttes dans divers écoulements laminaires et turbulents. Dans le cadre de l'approximation des équations de Stokes, les forces d'interactions hydrodynamiques entre deux gouttes sont déterminées de façon exacte au moyen d'un système de coordonnées bisphériques. Ces forces permettent de calculer les taux de collision pour des écoulements de cisaillement simple et d'élongation pure. Les sections de collision sont obtenues par trajectographie. Les résultats montrent un accroissement de l'efficacité de collision lorsque l'inertie des gouttes croît. L'expression du taux de collision en écoulement turbulent est analogue à celle utilisée en théorie cinétique des gaz. La fonction de distribution de paires est obtenue pour une turbulence homogène et isotrope et à l'aide de notions de probabilités conditionnelles. Les taux de collision sont calculés sans et avec la gravité. L'évolution granulométrique d'un nuage de gouttes est décrite à l'aide d'une équation cinétique de type Smoluchowski. Cette équation est résolue au moyen d'une méthode numérique originale. Deux cas sont traités. Pour de petites gouttes en écoulement turbulent homogène et isotrope, l'élargissement du spectre de gouttes demande plus de temps si les interactions hydrodynamiques sont prises en compte. Pour de grosses gouttes, compte tenu de la gravité, la formation par collisions de gouttes de diamètre supérieur à 200 microm dépend fortement de la turbulence. Les résultats corroborent les relevés \textit{in situ} dans des strati.

Cette thèse présente une étude numérique des interactions hydrodynamiques entre des particules et une paroi plane dans un fluide Newtonien, dans l'hypothèse d'un petit nombre de Reynolds. Les équations de Stokes sont d'abord transformées sous la forme classique d'une équation intégrale de surface donnant les contraintes sur la surface entourant le fluide. Nous avons calculé ainsi les forces et les couples exercés sur les particules. Une étude préliminaire correspond au problème axisymétrique pour une chaîne de particules sphériques arrangées suivant une ligne perpendiculaire à la paroi et en mouvement suivant cette ligne dans un fluide au repos. Puis nous avons abordé le problème général de particules non sphériques, la position des particules étant quelconque et le fluide loin des particules pouvant être en écoulement. Pour le cas axisymétrique, les points de collocation (qui sont aussi les emplacements des singularités de Stokes ou ``stokeslets'') sont choisis sur les surfaces des particules de façon que leurs distance relatives soient proportionnelles à la distance entre surfaces proches. Pour traiter le cas général, nous avons mis au point un code de calcul utilisant la méthode des éléments de frontière (BEM). Les stokeslets sont ici répartis de façon à prendre en compte les interactions hydrodynamiques et la complexité géométrique des configurations. En effet, le maillage est adapté à la variation du gradient de contraintes de façon que les zones de la surface qui entrent en forte interaction hydrodynamique soient les plus raffinées. La technique de raffinement dynamique de maillages que nous avons mise au point nous a permis de mieux détecter les zones de lubrification entre particules et paroi ainsi que les interactions entre particules. Les contraintes ont ainsi été calculées pour de nouvelles configurations géométriques. Enfin, les relations linéaires entre les forces et couples qui s'exercent sur les particules et leurs vitesses de translation et rotation sont exprimées au moyen de la « grande matrice de résistance » (ou de son inverse la « grande matrice de mobilité ») qui est alors introduite dans la relation fondamentale de la dynamique. En intégrant, nous avons déterminé les trajectoires des particules dans diverses configurations: en sédimentation dans un fluide au repos au voisinage d'une paroi, tournant librement dans un écoulement donné, ... La méthode permet aussi d'obtenir les champs de vitesse du fluide dans ces conditions.