Titre de l'écran-titre (visionné le 9 avril 2024) / Les hydroliennes à pale oscillante (HPO) complètement passives ont été optimisées avec les années pour atteindre des efficacités allant jusqu'à 51%. Ces hautes efficacités ont été démontrées sous l'hypothèse que l'écoulement est bidimensionnel et que la taille du canal n'affecte pas les performances de la turbine. Cependant, il en ressort un besoin de démontrer la robustesse de la turbine dans différents environnements d'écoulement. En effet, l'hydrolienne réelle est d'envergure finie, donc les extrémités de la pale engendrent des tourbillons de bout qui réduisent la portance en comparaison avec une turbine 2D. Aussi, l'aire de section de n'importe quel canal a une taille finie, donc la profondeur et la largeur de l'écoulement auront elles aussi un impact sur l'hydrolienne. Ce mémoire vise à démontrer la robustesse de la turbine à des changements géométriques de la forme du canal. Cette démonstration doit être perçue comme un commencement, et non une finalité, puisque de nombreux autres effets ne sont pas abordés ici. Cette robustesse est d'abord démontrée en évaluant qualitativement et quantitativement l'impact de la profondeur d'un canal 2D sur une HPO complètement passive qui fut optimisée par d'anciens étudiants au Laboratoire de Mécanique des Fluides Numérique. Il est observé que sans modification des paramètres structuraux qui influencent les mouvements non-contraints de la pale, la turbine devient chaotique et inopérable après une faible diminution de la profondeur du canal. Pour retrouver des niveaux de performance intéressants pour des canaux encore moins profonds, une méthode simple est proposée, basée sur un ajustement adéquat du générateur électrique. Grâce à cet ajustement, l'efficacité de la turbine augmente jusqu'à 70% pour des canaux hautement restreints. Enfin, la robustesse de la turbine est démontrée pour un canal 3D où la largeur du canal est variée. Encore ici, l'efficacité de la turbine augmente avec la diminution de l'aire de section du canal, à la différence que la force du générateur ne nécessite pas de modifications. Cette démonstration fut réalisée avec le souci de reproduire la turbine et les conditions d'écoulement correspondant à une étude expérimentale réalisée à la University of Victoria en parallèle à ce mémoire afin d'également montrer que la présente méthodologie numérique se compare bien aux résultats expérimentaux. L'importance de cette démonstration n'est pas à sous-estimer puisque dans le contexte des HPO complètement passives, aucune étude n'a encore validé les simulations sous des conditions d'opérations comparables, c'est-à-dire en incluant la tridimensionnalité de la turbine. / With the passing years, the fully-passive oscillating-foil turbines (OFT) have been optimised up to 51% of efficiency. These incredible efficiencies were made possible under the constrains that the flow is bidimensional and that the size of the channel does not affect the performances of the turbine. There is however a need to demonstrate the resilience of the turbine in different flow conditions. Indeed, the turbine is made up of a three-dimensional blade, so the presence of wing-tip vortices reduce the lift produced in comparison to 2D foils. Also, the cross-sectional area of any channel is finite, so the width and depth of the flow will also impact the turbine. The goal of this master thesis is to demonstrate the resilience of the turbine to geometric changes of the channel. This resilience demonstration must be seen as a start and not a finality since numerous other effects are not discussed here. This resilience is first demonstrated by qualitatively and quantitatively evaluating the channel's depth impact on a 2D fully-passive OFT that has been optimized by previous graduate students from the Laboratoire de Mécanique des Fluides Numérique. It has been observed that without modifications of the structural parameters influencing the unconstrained blade, the turbine becomes chaotic and uninteresting after a small decrease of the channel's depth. To retrieve good performances for shallow channels, a simple method is proposed, based on the adjustment of the electric generator. Thanks to this adjustment, the efficiency of the turbine increases up to 70% for highly confined channels. Finally, the resilience of the turbine is demonstrated again in a 3D channel where where the width is varied. Again, the efficiency of the turbine increases with a reduction of the cross-sectional area of the channel. Again, the efficiency of the turbine increased with a decrease of the cross-sectional area of the channel. The difference with the previous 2D study is that the force applied by the generator did not need any modifications. This demonstration has been planned to reproduce at the same time the turbine and the flow conditions used in a experimental study. The goal is to show that the current numerical methodology used compares itself well experimental results. The importance of this demonstration must not be underestimated since in the context of fully-passive OFTs, no study has adapted to compare both approaches in comparable operating conditions, meaning to include the 3D aspect of the turbine.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Un nouveau concept d’hydrolienne basée sur des ailes oscillantes est étudié. La présente étude a pour but d’étudier l’hydrodynamique instationnaire d’une aile oscillante, d’optimiser son mouvement afin de maximiser l’extraction de puissance et de démontrer le potentiel d’une turbine à ailes oscillantes par une campagne expérimentale sur un prototype. L’analyse et l’optimisation de la turbine à ailes oscillantes ont été effectuées par simulations numériques à bas nombre de Reynolds (laminaire) ainsi qu’à haut nombre de Reynolds (Unsteady Reynolds-Averaged Navier-Stokes; URANS). Une stratégie numérique 2D et 3D impliquant l’utilisation d’interfaces de glissement a été développée spécifiquement pour cette application de corps oscillants avec de grandes amplitudes de mouvement. À l’aide de cette stratégie numérique, une étude paramétrique fut effectuée et permit l’identification des paramètres dominants en ce qui a trait à la performance hydrodynamique de la turbine à ailes oscillantes. Basé sur un grand nombre de simulations, les zones optimales de production de puissance ont été identifiées dans les espaces paramétriques pertinents. De plus, des configurations spatiales optimales ont été identifiées pour le cas de turbines à ailes oscillantes en tandem. Le potentiel de l’hydrolienne à ailes oscillantes a été formellement établi dans ce travail grâce à une campagne expérimentale sur un prototype à ailes en tandem de 2 kW. La performance de ce dernier s’avéra compétitive avec celle des hydroliennes de type rotors à axe horizontal que l’on retrouve dans la majorité des designs d’hydrolienne proposés. Les données de la campagne expérimentale ont également permis de valider les résultats des simulations numériques par leur accord avec les simulations 3D. / A new concept of hydrokinetic turbine based on oscillating hydrofoils is investigated. The objective of this study is to analyze the unsteady hydrodynamics of oscillating foils, to optimize their motions for maximum power extraction and to demonstrate in practice the potential of such a concept of turbine through experiments on a prototype. The analysis and optimization have been conducted via low Reynolds number, laminar numerical simulations as well as high Reynolds number, Unsteady Reynolds- Averaged Navier-Stokes (URANS) computations. A 2D and 3D numerical methodology relying on the use of sliding interfaces and suitable to the case of foils undergoing oscillations of large amplitudes is presented. Using that numerical strategy, a parametric study is conducted and leads to the identification of the dominant parameters impacting the hydrodynamic performance of the oscillating-foil turbine. Based on a large number of simulations, the performance of the oscillating-foil turbine has been mapped in relevant parametric spaces. In addition, optimal spatial configurations of turbines with tandem foils is also provided. The potential of the oscillating-foils hydrokinetic turbine has also been formally established in this work through field tests on a 2kW tandem-foils prototype. Its performance has been found to be competitive with the best competing technologies based on horizontal-axis rotor-blades. The experimental data have also been used here to validate the numerical models and have been found to strongly support the 3D numerical simulations.

Due à la complexité des mécanismes typiquement requis pour contraindre l’aile d’une turbine à aile oscillante à suivre des mouvements spécifiques, cette thèse étudie la possibilité de bénéficier de mouvements non contraints, dits passifs. En pratique, cela implique que l’aile est attachée à la structure de la turbine à l’aide de supports élastiques indépendants en pilonnement et en tangage, formés de ressorts et d’amortisseurs. Par conséquent, seul un contrôle indirect des mouvements est possible en ajustant adéquatement les paramètres structuraux affectant la dynamique de l’aile, tels que les paramètres d’inertie, d’amortissement et de raideur de l’aile et de ses supports élastiques. En premier lieu, un prototype ayant des mouvements passifs autant en pilonnement qu’en tangage, et donc étant complètement passif, a été conçu et testé dans un canal à surface libre. Cette première phase du présent travail de recherche a confirmé la faisabilité et le potentiel de ce concept en permettant d’extraire une quantité significative d’énergie de l’écoulement d’eau. Cependant, l’efficacité maximale atteinte est demeurée inférieure à ce qui peut être obtenu en contraignant l’aile à suivre des mouvements précis. Suite à ces expériences, un algorithme résolvant la dynamique du solide a été implémenté et couplé au logiciel résolvant la dynamique du fluide gouverné par les équations de Navier-Stokes. Des simulations numériques ont été réalisées afin d’analyser plus en détail la dynamique de chacun des deux degrés de liberté de l’aile. Plutôt que de poursuivre notre étude du concept complètement passif immédiatement, un concept de turbine semi-passive caractérisée par un mouvement de tangage passif et un mouvement de pilonnement contraint a été considéré. Des efficacités de l’ordre de 45% ont été atteintes, se comparant ainsi aux meilleures performances rapportées dans la littérature concernant les turbines à ailes oscillantes complètement contraintes. En plus de révéler le fort potentiel de ce concept de turbine semi-passive, cette étude nous a permis de nous concentrer sur certains aspects spécifiques concernant la dynamique d’une aile attachée par des ressorts en tangage. Cette analyse plus détaillée de la physique en jeu a été facilitée par le nombre réduit de paramètres structuraux en jeu par rapport à une turbine pour laquelle le mouvement de pilonnement est lui aussi passif. L’une des découvertes importantes est que le centre de masse doit être situé en aval du point de pivot afin de générer un transfert d’énergie du mouvement de pilonnement vers le mouvement de tangage par l’entremise du couplage inertiel entre les deux degrés de liberté. Ce transfert d’énergie est crucial puisque les mouvements de tangage optimaux nécessitent de l’énergie en moyenne pour être soutenus. De plus, un paramètre combinant les effets liés au moment d’inertie de l’aile par rapport à son point de pivot et à la raideur en tangage a été proposé. Ce paramètre permet de bien caractériser la dynamique du mouvement de tangage passif de la turbine semi-passive. Il permet aussi de déterminer la raideur requise pour différentes valeurs du moment d’inertie afin de maintenir une performance optimale de la turbine. Utilisant les connaissances acquises concernant la dynamique des mouvements de tangage passifs, le concept de turbine à aile oscillante complètement passive a été revisité. Les meilleures efficacités obtenues avec la turbine semi-passive ont été égalées et ont même été surpassées puisque qu’une efficacité de 53.8% a été atteinte. Les résultats ont aussi démontré qu’une performance optimale pouvait être maintenue sur de larges plages de valeurs en ce qui concerne la masse en pilonnement ainsi que le moment d’inertie par rapport au point de pivot, pourvu que les raideurs en pilonnement et en tangage soient ajustées correctement. / Due to the complexity of the mechanisms typically required when designing a flapping-foil turbine to prescribe specific heave and pitch motions, this thesis investigates the possibility of benefiting from unconstrained motions. In practice, this means that the foil is attached to the turbine structure with independent elastic supports in heave and in pitch, which consist in springs and dampers. Consequently, only an indirect control over the foil motions is possible through an adequate adjustment of the structural parameters affecting the foil dynamics, namely the inertial, damping and stiffness characteristics of the elastically-supported foil. Such motions are referred to as passive motions. As a first step, a turbine prototype with passive heave and pitch motions, thus being fully-passive, has been designed and tested in a water channel. This first phase of the present research work has confirmed the feasibility and the potential of this concept to extract a significant amount of energy from a fluid flow. However, the maximum efficiency that has been obtained is smaller than what can be achieved when prescribing specific foil motions. Following these experiments, a solid solver has been implemented and coupled with a Navier-Stokes fluid solver. Numerical simulations have been carried out to analyze the dynamics of both degrees of freedom in more details. Instead of immediately pursuing our study of the fully-passive flappingfoil turbine, a semi-passive concept, with a passive pitch motion and a prescribed heave motion, has been considered. Efficiencies of the order of 45% have been achieved, hence competing with the best performance reported in the literature for flapping-foil turbines with prescribed motions. In addition to revealing the great potential of this semi-passive turbine concept, this study has allowed us to focus on some specific aspects of the dynamics of passive pitch motions. This more detailed analysis of the physics at play has been facilitated by the reduced number of structural parameters affecting the foil dynamics compared to a turbine for which the foil is also elastically-supported in heave. One of the main findings is that the center of mass must be positioned downstream of the pitch axis in order to generate a net transfer of energy from the heave motion to the pitch motion via the inertial coupling between the two degrees of freedom. This energy transfer is crucial because optimal pitch motions require energy on average to be sustained. Moreover, a parameter combining the effects of the moment of inertia of the foil about the pitch axis and the pitch stiffness has been proposed. This parameter effectively characterizes the pitch dynamics of the semi-passive turbine. It also allows properly scaling the pitch stiffness when different moments of inertia are considered with the objective of maintaining an optimal turbine performance. Having improved our knowledge about the dynamics of passive pitch motions, the fully-passive flapping-foil turbine concept has been revisited. The best efficiencies obtained with the semi-passive concept have been matched, and even exceeded since an efficiency of 53.8% has been reached. The results have also demonstrated that an optimal performance can be maintained over large ranges of values regarding the heaving mass and the moment of inertia when the heave and pitch stiffness coefficients are adjusted adequately.

Nous étudions expérimentalement le comportement d'un film liquide, dit de Landau-Levich, lorsqu'il s'évapore dans une atmosphère constituée uniquement de sa vapeur.La dynamique de ce type de film est un paramètre-clef qui gouverne le fonctionnement des caloducs oscillant (en anglais PHP - Pulsating Heat Pipes). Les PHP sont des liens thermiques de forte conductance. Les recherches récentes leur attribuent un pouvoir de refroidissement très élevé, ce qui les rend particulièrement convoités par l'industrie. Leur géométrie est simple : il s'agit d'un tube capillaire enroulé en plusieurs branches entre une partie froide (condenseur) et une partie chaude (évaporateur). Le tube est rempli d'un fluide pur diphasique, c'est-à-dire présent sous la forme d'une succession de bulles de vapeur et de bouchons de liquide. Lorsque la différence de température entre l’évaporateur et le condenseur dépasse un certain seuil, les bulles et bouchons commencent à osciller dans le tube, entre les deux parties, ce qui permet au PHP de transférer la chaleur.Notre installation expérimentale représente un PHP dans sa configuration la plus simple, à branche unique. Une interface liquide-vapeur oscille dans un tube de section rectangulaire, et dépose un film liquide à chaque passage. Nous nous intéressons au mécanisme qui permet l'entretien de l'oscillation de l'interface, et fixe sa fréquence. L'équation de mouvement obtenue prend en compte la dissipation visqueuse engendrée par un écoulement oscillant. Dans les modèles actuels de PHP, l'hypothèse d'un écoulement de type Poiseuille est formulée. Or, notre approche montre que l'hypothèse d'un écoulement faiblement inertiel est mieux adaptée, conduit à une dissipation deux fois supérieure.Le dispositif expérimental permet l'observation du film. Une combinaison originale de méthodes optiques permet également de mesurer sa longueur et son épaisseur, et de reconstruire son profil 3D à chaque instant. Nous pouvons suivre l'évolution du film tout au long de sa durée de vie, et ainsi analyser son comportement dynamique. Le film est presque plat (pente inférieure à 0,1°). Sur toute sa longueur, qui est de quelques centimètres, cela correspond à une variation de son épaisseur de moitié, la valeur moyenne étant de 50 microns. Sous l'effet du chauffage, le film se rétracte progressivement. Dès le début de son évaporation, un bourrelet de démouillage est formé sur le pourtour du film, près de la ligne triple. L'apparition de ce bourrelet est caractéristique d'un démouillage visqueux sous conditions de non-mouillage. Ce comportement est surprenant, dans la mesure où nous avons choisi un fluide mouillant parfaitement la paroi en conditions isothermes. A l'échelle nanométrique, au plus près de la ligne triple, l'angle de contact entre le liquide et la paroi est donc très faible. Nous mesurons cependant un grand angle apparent (c'est-à-dire visible à l'échelle millimétrique), qui augmente avec la surchauffe de la paroi. Dès l'augmentation de cet angle, le bourrelet de démouillage se forme, et le film se rétracte. Ce phénomène est expliqué par l'évaporation à l'échelle microscopique. Les résultats expérimentaux sont en accord quantitatif avec la théorie développée par d'autres chercheurs. / We experimentally study the behavior of liquid films - so called Landau-Levich films - when they evaporate in their pure vapor atmosphere.The dynamics of this film is a key parameter that rules out the functioning of Pulsating Heat Pipes (PHPs). PHPs are high conductive thermal links. Their heat transfert capability is known to be extremely high. For this reason they are promising for numerous industrial applications. Their geometry is simple. It is a capillary tube bent in several branches that meander between a hot part (called evaporator) and a cold part (called condenser), and filled up with a pure two-phase fluid. When the temperature difference between evaporator and condenser exceeds a certain threshold, gas bubbles and liquid plugs begin to oscillate spontaneously back and forth inside the tube and PHP starts transferring the heat.Our experimental setup features the simplest, single branch PHP. A liquid/vapor interface oscillates in a tube. It deposits a liquid film at each passage. We focus first on the mecanism which makes possible self-sustained interface oscillations and defines its frequency. The obtained motion equation accounts for the viscous dissipation caused by oscillatory flow. In existing PHP modelling, a laminar flow is supposed. Yet, our approach shows that the assumption of weakly inertial flow is preferable and leads to a dissipation rate twice larger that the Poiseuille flow.The experimental setup allows the film visualization. An original combination of optical measurement techniques lets us measure the film length, thickness and 3D-profile at all times. The film evolution has been measured during its whole lifetime. The film is nearly flat (its slope is smaller than 0,1°). The film length is of several centimeters, and the average thickness is 50 microns. Thus, along the total length, its thickness decreases by half. Under heating conditions, the film gradually recedes. A dewetting ridge is formed, near the triple contact line. Such a behavior is typical under non-wetting conditions. At the nanometric scale the contact angle between the liquid and the solid wall is very low. However, we measure a large apparent contact angle (visible at the millimetric scale) which increases with the wall superheating. Once this angle increases, the dewetting ridge is formed and the film recedes. The large apparent contact angle is explained by evaporation in the microscopic vicinity of the contact line. The measured apparent contact angle value agrees quantitatively with theoretical results obtained by other researchers.

Les vésicules sont des bicouches fermées de molécules tensioactives, remplies de liquide, à l'intérieur d'un autre liquide. Leur taille peut être comprise entre dix et 100 microns : elles sont alors dites géantes. Nous nous intéressons à la dynamique de deux de ces objets dans un cisaillement simple, c'est à dire l'écoulement d'un liquide situé entre deux plaques planes se translatant l'une par rapport à l'autre à vitesse et distance constante. Nous commençons par une étude asymptotique, pour des vésicules quasi-sphériques en interaction lointaine. Nous utilisons ensuite un code de calcul basé sur la méthode des éléments de frontière pour étudier le cas de vésicules moins sphériques et plus proches, et comparons les résultats obtenus avec des expériences. Nous présentons enfin comment cette étude peut être utilisée pour prédire certaines propriétés de diffusion d'une suspension de vésicules, dans le régime semi-dilué, où seul le détail des interactions à deux corps est considéré.

Les vésicules sont des bicouches fermées de molécules tensioactives, remplies de liquide, à l'intérieur d'un autre liquide. Leur taille peut être comprise entre dix et 100 microns : elles sont alors dites géantes. Nous nous intéressons à la dynamique de deux de ces objets dans un cisaillement simple, c'est à dire l'écoulement d'un liquide situé entre deux plaques planes se translatant l'une par rapport à l'autre à vitesse et distance constante. Nous commençons par une étude asymptotique, pour des vésicules quasi-sphériques en interaction lointaine. Nous utilisons ensuite un code de calcul basé sur la méthode des éléments de frontière pour étudier le cas de vésicules moins sphériques et plus proches, et comparons les résultats obtenus avec des expériences. Nous présentons enfin comment cette étude peut être utilisée pour prédire certaines propriétés de diffusion d'une suspension de vésicules, dans le régime semi-dilué, où seul le détail des interactions à deux corps est considéré. / Vesicles are closed bilayers of tensioactive molecules, filled with liquid, inside another liquid. Their size can be between 10 and 100 microns : in this case, they are called giant vesicles. We study the dynamic of two of these objects in a simple shear flow, which is the one of a liquid sheared between two walls translating with respect to each other at a constant speed and distance. We begin by an asymptotic study, for quasispherical vesicles in the far field interacting regime. We then use a numerical code based on the boundary element method to study the case of less spherical and closer vesicles, and compare our results with experiments. We finish by presenting how this study can be used to predict some diffusing properties of a sheared suspension of vesicles, in the semidilute regime, where only the details of two body interactions are considered.

La texturation de surface est une thématique récente qui suscite un certain engouement pour les contacts dynamiques. Pendant de nombreuses années, les tribologues ont privilégié les surfaces lisses aux faibles rugosités pour limiter le frottement. Inspiré des rugosités de surface organisées observées dans la nature, les topologies de surfaces sont désormais axées sur la texturation et la structuration des rugosités. Fort de ce potentiel, de nombreuses études traitent, par une approche numérique, la modélisation de ces surfaces et les études expérimentales sont rares, avec une instrumentation souvent insuffisante pour appréhender tous les phénomènes physiques. Par une approche expérimentale, nous analysons le comportement des butées hydrodynamiques à faces parallèles partiellement texturées. Les 80 capteurs équipant le dispositif d'essais permettent d'apprécier avec rigueur la phénoménologie à l'interface du patin et du film lubrifiant. L'analyse met l'accent sur la capabilité de ce composant à être intégré dans un environnement industriel. Afin d'objectiver les résultats, les campagnes d'essais sont menées sur dix butées hydrodynamiques dont quatre sont munies de texturation. Une comparaison de ces butées facilite leur classement en termes de capacité de charge, de réduction de frottement ou encore de risque d'usure dans les phases de démarrage. Pour les configurations étudiées, les butées texturées ne peuvent concurrencer les butées à poches ou à plans inclinés du point de vue de la capacité de charge. En se référant à une butée à faces parallèles, les butées texturées permettent une réduction du frottement de 30% à faibles charges tandis que pour de fortes charges, les / Surface texturing is a recent topic which has raised a great interest in contact dynamics. For many years, engineers have favored smooth surfaces with low roughness in order to minimize friction losses. Inspired by textured surfaces which can be commonly found in nature, the research in surface topography is now focused on texturing and roughness characterization. Considering the great potential of surface texturing, many research studies analyze this subject, most commonly theoretically, while experimental works are often performed with inadequate equipment which does not allow a proper evaluation of the involved physical phenomena. This study uses an experimental approach in order to analyze the behavior of hydrodynamic thrust bearings with parallel textured pads. The experimental device is equipped with 80 sensors which allow a proper assessment of the phenomenology at the film/pad interface. This analysis focuses on the capability of this textured component to be integrated in an industrial environment. To objectify the results, the tests are conducted on ten hydrodynamic thrust bearings, among which four are textured. The comparison between the performances of these bearings facilitates their classification in terms of load-carrying capacity, friction loss and wears resistance during the start-up period. Results show that for the studied configurations, the textured thrust bearings cannot compete with pocketed or tapered land thrust bearings in terms of load-carrying capacity. In the case of parallel thrust bearings, surface texturing can help to reduce friction up to 30% at low loads while for heavy loads, their performance is equivalent or even lower than that of

Reproduction de : Thèse de doctorat : Mécanique : Lille 1 : 2004. / N° d'ordre (Lille 1) : 3573. Titre provenant de la page de titre du document numérisé. Bibliogr. p. 170-180.

L'instabilité de Vishniac est invoquée pour expliquer la fragmentation et la filamentation de la coquille fine de matière choquée présente dans les restes de supernova radiatifs. Toutefois l'implication et les conséquences de ce processus spécifique sur la morphologie complexe des restes de supernova ne sont pas entièrement démontrées. Nous avons réalisé des simulations numériques 2D d'une onde de souffle perturbée qui se propage dans un milieu ambiant homogène afin de mieux comprendre le mécanisme de l'instabilité de Vishniac. Le code hydrodynamique HYDRO-MUSCL 2D a été utilisé afin de déclencher la perturbation de la coquille fine dans les géométries plan-parallèle et sphérique. Dans ce travail, je montrerai que nous avons obtenu la surstabilité de Vishniac (" Vishniac overstability ") comme prédit par l'analyse théorique. Ce processus repose sur une oscillation d'amplitude et de période temporelle croissante des grandeurs fluides et spatiales. Nous avons testé l'effet de paramètres spécifiques sur la variation de masse d'une région donnée dans le cadre d'une étude paramétrique que nous avons réalisée sur le supercalculateur Titane. Nous avons trouvé que la perturbation est atténuée après quelques oscillations de la variation de masse pour tous les jeux de paramètres. Ainsi nous concluons que dans notre modèle, l'instabilité de Vishniac ne permet pas la fragmentation de la coquille fine à cause d'effets qui ne sont pas pris en compte par l'analyse théorique.

Une étude préliminaire de la chimie de l'eau de la nappe de l'Isère, établie lors d'une campagne de mesures de la résistivité et du niveau de l'eau pendant l'été 1980 pour EDF, avait permis de se faire une première idée sur les différences existant entre les conditions hydrochimiques et hydrodynamiques qui régissent l'aquifère. La présence d'une seule nappe dans l'aquifère alluvial de la vallée de l'Isère n'implique pas l'uniformité de la qualité de l'eau sur les plans verticaux et horizontaux. Le milieu alluvial semble donc étre homogène et isotrope pour la transmission des pressions, mais pas pour la transmission de matière. Il doit donc exister des conditions de mélange des eaux et des interactions eau/milieu liées à la structure de l'aquifère. Ces premiers résultats ont servi de motivation et de point de départ à notre étude. La nappe de la vallée de l'Isère, site retenu comme secteur de travail, est limitée à l'amont par la confluence Isère-Drac en aval de Grenoble, à l'aval par une remontée du mur de la nappe au seuil de Rovon (35 km à l'aval de Grenoble, en direction de Valence). Elle est alimentée par les apports de versant, par la rivière, et par la pluie; mais cette derniere est négligeable par rapport aux précédentes. Dans l'aquifère de la vallée de l'Isère, nous retrouvons deux secteurs bien définis du point de vue morphologique (largeur de la nappe), hydrodynamique (type et vitesse d'écoulement), hydrochimique (gradients de minéralisation différents): la Cluse de l'Isère et la plaine de Moirans.