Les ondes internes sont produites par suite de l’équilibre dynamique entre les forces de flottabilité et la gravité quand une particule de fluide est déplacée verticalement dans un milieu stratifié stable. Les systèmes géophysiques tels que océan et l’atmosphère sont naturellement stratifiés et donc favorables à la propagation des ondes internes. En outre, ces deux environnements stockent une grande quantité de particules tant dans leur intérieur que sur les bords. Par conséquent, les ondes internes et les particules vont inévitablement interagir dans ces systèmes. Au cours de ce travail, des expériences exploratoires sont réalisées pour étudier le transport par érosion des particules, généré par les ondes internes. Afin de déterminer un seuil de transport, les propriétés particulières des réflexions d’ondes internes («réflexion critique ») sont utilisées pour augmenter l’intensité du champ d’ondes à la surface de réflexion. Une méthode a été développée en collaboration avec une équipe de traitement du signal pour améliorer la détermination des composantes de l’onde impliquées dans une réflexion quasi critique. Cela nous a permis de comparer nos résultats expérimentaux avec une théorie de la réflexion critique, montrant un bon accord et permettant d’extrapoler ces résultats à des expériences au-delà de la nôtre et à des conditions océaniques. Nous avons aussi étudié l’interaction des ondes internes avec une colonne de particules en sédimentation. Deux effets principaux ont été observés : la colonne oscille autour d’une position d’équilibre, et elle est déplacée dans son ensemble. La direction du déplacement de la colonne est expliquée par le calcul de l’effet de la dérive Lagrangienne produite pour des ondes. Cet effet pourrait également expliquer la dépendance en fréquence du déplacement. / Internal waves are produced as a consequence of the dynamic balance between buoyancy and gravity forces when a particle of fluid is vertically displaced in a stably stratified environment. Geophysical systems such as ocean and atmosphere are naturally stratified and therefore suitable for internal waves propagation. Furthermore, these two environments stock a vast amount of particles at their boundaries and in their bulk. Therefore, internal waves and particles will inexorably interact in these systems. In this work, exploratory experiments are performed to study wave generated erosive transport of particles. In order to determine a transport threshold, the peculiar properties of internal waves (“critical reflection”) are employed to increase the intensity of the wave field at the boundaries. A method was developed in collaboration with a signal processing team to improve the determination of the wave components involved in near-critical reflection. This method enabled us to compare our experimental results with a theory of critical reflection, showing good agreement and allowing to extrapolate these results to experiments beyond ours and to oceanic conditions. In addition, we study the interaction of internal waves with a column of particles in sedimentation. Two main effects are observed: the column oscillates around an equilibrium position, and it is displaced as a whole. The direction of the displacement of the column is explained by computing the effect of the Lagrangian drift of the waves. This effect could also explain the frequency dependence of the displacement.

La cascade d’énergie qui a lieu dans les océans, depuis les grandes vers les petites échelles, est capitale pour comprendre leur dynamique et le mélange irréversible associé. Les attracteurs d’ondes internes font partie des mécanismes conduisant potentiellement à une telle cascade. Dans ce manuscrit, nous étudions expérimentalement les attracteurs d’ondes internes, dans une cuve trapézoïdale remplie d’un fluide stratifié linéairement en densité. Dans cette géométrie, les ondes peuvent être focalisées vers un cycle limite : l’attracteur. Nous montrons que la formation de l’attracteur est purement linéaire : des petites échelles sont donc créées grâce à la focalisation des ondes. Les principales caractéristiques de l’attracteur dépendent uniquement de la géométrie trapézoïdale de la cuve. A l’échelle de l’océan, nous montrons que les attracteurs d’ondes internes sont très probablement instables. En effet, ceux-ci sont sujets à une instabilité de résonance triadique, qui transfère de l’énergie depuis l’attracteur vers un couple d’ondes secondaires. Cette instabilité et ses principales caractéristiques sont décrites en fonction de la géométrie du bassin. Pour des expériences de longue durée, l’instabilité produit plusieurs paires d’ondes secondaires, créant une cascade d’instabilités triadiques et transférant l’énergie injectée à grandes échelles vers des échelles plus petites. Nous montrons, pour la première fois de façon expérimentale, de très fortes signatures de turbulence d’ondes internes. Au delà de cet état, la cascade atteint un régime de mélange partiel du fluide stratifié. Cet ultime régime apparait indépendant de la géométrie trapézoïdale du bassin, et donc, universel. Cette thèse est complétée par une étude sur la masse ajoutée et l’amortissement par émission d’ondes d’objets oscillant horizontalement dans un fluide stratifié en densité. Cela a des applications concernant la conversion de l’énergie des marées en ondes internes. / A question of paramount importance in the dynamics of oceans is related to the energy cascade from large to small scales and its contribution to mixing. Internal wave attractors may be one of the possible mechanisms responsible for such a cascade. In this manuscript, we study experimentally internal wave attractors in a trapezoidal test tank filled with linearly stratified fluid. In such a geometry, the waves can form closed loops called attractors. We show that the attractor formation is purely linear: small scales are thus created by wave focusing. The attractor characteristics are found to only depend on the trapezoidal geometry of the tank. At the ocean scale, we show that attractors are very likely to be unstable. Indeed, internal wave attractors are prone to a triadic resonance instability, which transfers energy from the attractor to a pair of secondary waves. This instability and its main characteristics are described as a function of the geometry of the basin. For long-term experiments, the instability produces several pairs of secondary waves, creating a cascade of triadic interactions and transferring energy from large-scale monochromatic input to multi-scale internal-wave motion. We reveal, for the first time, experimental convincing signatures of internal wave turbulence. Beyond this cascade, we have a mixing regime, which appears to be independent of the trapezoidal geometry and, thus, universal. This manuscript is completed by a study on added mass and wave damping coefficient of bodies oscillating horizontally in a stratified fluid, with applications to tidal conversion.

Au cours de cette thèse, nous nous sommes intéressés à la modélisation des interactions fluide-structure entre un fluide visqueux, incompressible et une structure pouvant être déformable. Après avoir présenté les différentes approches possibles de modélisation, nous introduisons une nouvelle méthode de type frontière immergée : la méthode IPC ("Image Point Correction"). Combinant approches Ghost-Cell et Pénalisation, cette méthode mixte du second degré globalement et localement en vitesse, est validée sur différents cas tests (comparaisons des coefficients aérodynamiques pour des cylindres fixes ou mobiles, sédimentation 2D d'un cylindre). Nous avons ensuite appliqué la méthode IPC à la simulation de la nage. Dans un premier temps, le solveur 2D a été couplé avec un algorithme d'optimisation mathématique afin de déterminer la loi de nage optimale pour une géométrie de poisson donnée. Puis, dans un second temps, nous avons simulé la nage 3D après reconstruction approchée de la géométrie, basée sur des images du nageur. Enfin, grâce à l'outil du squelette, une reconstruction réaliste du poisson est proposée. / The aim of this thesis is to investigate the modeling of fluid-structure interactions. The fluid is viscous and incompressible, and the structure is subject to an imposed deformation. After a survey of the different existing approachs to model fluid-structure interactions, we introduce a new immersed boundary method: the IPC method (”Image Point Correction”). This ap-proach merges Ghost-Cell and Penalty concepts. It is globally and locally second order in velocity, and it is validated through several canonical simulations. Then, we apply the IPC method to fish-like swimming. First, the 2D solver is paired up with a mathematical optimization algorithm to determine the optimal swimming law for a given fish geometry. Secondly, we simulate a 3D swimmer after performing an approximated reconstruction of the geometry based on actual fish pictures. Finally, thanks to the skeleton approach, a realistic reconstruction of the fish is exposed.

Le récent développement des menaces terroristes renforce l'effort de recherche du CEA et d'EDF pour la protection des citoyens et des installations. De nombreux scénarios doivent être envisagés comme, par exemple, la chute d'un avion de ligne sur une structure de génie civil. La dispersion du carburant dans la structure, son embrasement sous forme de boule de feu et les effets thermiques associés sont des éléments essentiels du problème. L'utilisation de modèles numériques est indispensable car des expériences seraient difficiles à mettre en œuvre, coûteuses et dangereuses. Le problème type que l'on cherche à modéliser est donc l'impact d'un réservoir rempli de fluide, sa déchirure et la dispersion de son contenu. C'est un problème complexe qui fait intervenir une structure mince avec un comportement fortement non-linéaire allant jusqu'à rupture, un fluide dont la surface libre peut varier drastiquement et des interactions fluide-structure non permanentes. L'utilisation des méthodes numériques traditionnelles pour résoudre ce problème semble difficile, essentiellement parce qu'elles reposent sur un maillage. Cela complique la gestion des grandes déformations, la modélisation des interfaces variables et l'introduction de discontinuités telles que les fissures. Afin de s'affranchir de ces problèmes, la méthode sans maillage SPH (\og Smoothed Particle Hydrodynamics \fg) a été utilisée pour modéliser le fluide et la structure. Ce travail, inscrit dans la continuité de recherches précédentes, a permis d'étendre un modèle de coque SPH à la modélisation des ruptures. Un algorithme de gestion des interactions fluide-structure a également été adapté à la topologie particulière des coques. Afin de réduire les coûts de calcul importants liés à ce modèle, un couplage avec la méthode des éléments finis a également été élaboré. Il permet de n'utiliser les SPH que dans les zones d'intérêt où la rupture est attendue. Finalement, des essais réalisés par l'ONERA sont étudiés pour valider la méthode. Ces travaux ont permis de doter le logiciel de dynamique rapide Europlexus d'un outil original et efficace pour la simulation des impacts de structures minces en interaction avec un fluide. Un calcul démonstratif montre enfin la pertinence de l'approche et sa mise en œuvre dans un cadre industriel. / The recent development of terrorist threats increases the research effort of the french Atomic Energy Commission (CEA) and the French Electricity company (EDF) for the protection of citizens and facilities. Many scenarios should be considered as, for example, the fall of an airliner on a civil engineering structure. The dispersion of fuel in the structure, the formation of a fireball and associated thermal effects are essential elements of the problem. The use of numerical models is essential because experiences would be difficult to organize, costly and dangerous. The typical problem that we want to model is the impact of a tank filled with fluid, its rupture and the dispersion of its contents. It is a complex problem which involves a thin structure with a highly non-linear behavior up to rupture, a fluid with a free surface that can vary drastically and non permanent fluid-structure interactions. The use of traditional numerical methods to solve this kind of problems is difficult, mainly because they rely on a mesh. This complicates the management of large deformations, the modeling of moving interfaces and the introduction of discontinuities such as cracks. To overcome these problems, the meshfree method SPH (Smoothed Particle Hydrodynamics) was used to model both the fluid and the structure. This work, which is a continuation of previous research, has extended a model of SPH shell to the modeling of ruptures. An algorithm for managing fluid-structure interactions has also been adapted to the particular topology of shells. To reduce the important computational costs associated with this model, a coupling with the finite element method was also developed. It allows the use of SPH in areas of interest where the rupture is expected. Finally, tests performed by the french Aerospace Lab (ONERA) are studied to validate the method. This work helped to provide fast dynamic software Europlexus an original and effective tool for the simulation of the impact of thin structures interacting with fluid. A demonstrative calculation finally shows the relevance of the approach and its use within an industrial framework.

Afin d’augmenter le rendement et la puissance spécifique des turbines à gaz, la température en entrée de la turbine haute pression a été continuellement augmentée, au point d’atteindre une valeur de l’ordre de la température de fusion des matériaux. La capacité à prédire l’écoulement (température en particulier) dans la turbine est donc un élément essentiel pour la conception des turbines à gaz, particulièrement pour celle qui ne sont pas refroidies. Toutefois, cette tâche est rendue difficile par l’extrême complexité de l’environnement dans lequel la turbine évolue (interaction avec la chambre de combustion, effets technologiques, couplage thermique fluide / solide, etc.).L’objectif de ce travail de thèse est de proposer, d’une part une stratégie de simulation numérique afin de prédire au mieux les champs aérothermiques dans une turbine haute-pression, d’autre part une méthode d’analyse permettant de quantifier l’impact des différents éléments de l’environnement sur les performances. Pour répondre à ces objectifs, des simulations instationnaires d’un étage de turbine ont été réalisées avec le code elsA, prenant en compte les effets technologiques (baignoire, refroidissement externe grâce aux évents du distributeur et aux cavités moyeu et carter du rotor qui sont alimentés par le système d’air secondaire, congés de raccordement) et les transferts thermiques conjugués. Une analyse de la production locale d’entropie a aussi été menée afin de comparer les performances aérodynamiques pour plusieurs niveaux d’approximations dans la définition de la turbine (prise en compte ou non de tel ou tel effet technologique) et de localiser l’origine de ces différences. L’analyse des résultats a montré que la hauteur de jeu et le système de refroidissement externe ont l’impact le plus significatif sur les performances aérodynamiques de la turbine. La température de paroi de la pale de rotor est de son côté fortement influencée par l’écoulement issu de la chambre de combustion, le refroidissement externe et le couplage thermique fluide / solide. Ce travail est un premier pas vers la réalisation de simulations totalement intégrées qui doivent permettre d’améliorer la précision des conceptions. / In order to increase the thermodynamic efficiency of gas turbine engines, the high-pressure turbine inlet temperature has been continually increased up to reach levels of the order of magnitude of the vanes and blades melting temperatures. The ability of predicting the flow through the turbine(especially the temperature) is a key point for the design of gas turbines, especially for the uncooled ones. However, this is challenging because of the complex environment that interacts with the turbine(hot-streak migration, technological details, fluid/solid thermal coupling . . .).The aim of this work is to develop a strategy based on CFD in order to predict aerothermal fields in a high-pressure turbine as well as an analysis allowing to quantify the impact of the environment on the turbine performances. To achieve these goals, the elsA code has been used to perform unsteady simulation of a turbine stage, taking into account technological details (squealer tip, external cooling thanks to vane trailing edge cooling and rotor hub and shroud cooling cavities, fillets) and conjugate heat transfers. An analysis of the local entropy production rate was also used to compare the performances resulting from different modeling of the turbine (taking into account or not one particular technological detail) and to localize the origin of these discrepancies. The results show that the tip height and the external cooling have the greatest impact on the turbine performances. The rotor blade temperature is mainly affected by the flow coming from the combustion chamber, the external cooling and the fluid / solid thermal coupling. This work is the first step towards the realization of integrated simulations allowing to improve the accuracy of design.

Autour de nous, les systèmes granulaires sont omniprésents et rarement dissociés d'une phase liquide ou gazeuse. L'interaction entre les deux phases (solide et fluide) implique une dynamique d'écoulement complexe de ces systèmes couplés. Cette thèse expérimentale a pour but de caractériser la dynamique d'effondrement gravitaire instationnaire d'une colonne fluide-grains dans différentes configurations de mélange diphasique. Considérant uniquement l'effet des grains au sein du fluide, la dynamique d'écoulement d'une suspension isodense dans l'air est d'abord abordée avec une description macroscopique. En particulier, La rhéologie du fluide apparent est extraite à partir de l'évolution temporelle du front en utilisant des solutions auto-similaires, modèles de la dynamique de propagation aux temps longs. Sur la gamme du nombre de Reynolds étudiée, la dynamique d'écoulement de la suspension est décrite comme celle d'un fluide apparent Newtonien ou non-Newtonien (rhéofluidifiant, rhéoépaississant, viscoplastique) dépendant des paramètres considérés (fraction volumique, viscosité du fluide interstitiel, diamètre des grains, protocole de préparation). Afin de décrire la configuration opposée d'un massif sédimentaire pour lequel les interactions entre particules solides deviennent dominantes, la configuration d'une colonne granulaire dense saturée en fluide, i.e. où les grains sont plus lourds que le fluide porteur, est ensuite étudiée. Un travail préliminaire est consacré à la caractérisation de l'effondrement d'une colonne granulaire totalement immergée dans un fluide. Le rôle dissipatif de ce dernier sur le milieu granulaire est mis en évidence par une analyse de la dynamique d'effondrement et des caractéristiques du dépôt final. Cette caractérisation permet de classifier les régimes d'effondrement en fonction de la viscosité et de la masse volumique du fluide environnant, i.e. du nombre de Stokes et du rapport de masse volumique entre le fluide et les grains. Dans le cas triphasique, i.e. lorsque la colonne fluide-grains s'effondre dans l'air, la dynamique peut être fortement affectée par les effets capillaires à travers le nombre de Bond qui contrôle alors la mise en mouvement de la colonne initiale. Quand ces effets deviennent négligeables à l'échelle de la colonne et à celle du grain, le fluide interstitiel peut jouer un rôle moteur ou dissipatif vis-à-vis dumilieu granulaire conduisant à une longueur d'étalement plus ou moins importante en comparaison au cas sec. Le rôle du fluide interstitiel dépend essentiellement de sa viscosité modifiant, par la même occasion, la dynamique d'effondrement. Enfin, une étude préliminaire sur la dynamique d'écoulement d'une suspension non-isodense, initialement homogène, est réalisée. Cette configuration, à la transition des situations décrites précédemment, permet d'aborder le couplage de la dynamique de sédimentation des grains et celle du courant. En particulier, nous observons que la vitesse de sédimentation décroît avec l'augmentation de la fraction volumiqueinitiale en particules. / On the Earth's surface, granular medias are ubiquitous and they are rarely dissociated from a liquid or a gas. The fluid-solid interaction leads to a complex flow dynamics of these coupled systems. This experimental work aims at characterizing the dynamics of the unsteady gravitydrivencollapse of a granular-fluid column within different configurations of the diphasic mixture. First, the flow dynamics of a neutrally buoyant suspension in air are characterized based on a macroscopic description of the flow. In particular, the rheological parameters of the apparent fluid have been extracted using the temporal evolution of the propagating front and self-similar solutions, models of the propagating dynamics at long times. In the considered range of the Reynolds number, the flow dynamics are described as an apparent Newtonian or non-Newtonian (shear-thinning/-thickening, viscoplastic) fluid depending on the various parameters (volume fraction, viscosity of the interstitial fluid, particle diameter, mixing protocol). In order to describe the opposite case of a sedimentary environment where particle-particle interaction becomes dominant, a second part of this work investigates the case of a fluid-saturated granular collapse, i.e. for which particles are heavier than the carrier fluid, in a dense packing configuration. For this purpose, a first part of the study is dedicated to characterize the collapse of an immersed granular column. The dissipative role of the fluid on the granular media is highlighted by an analysis of the collapse dynamics and the characteristics of the final deposit. This characterization allows to classify the regimes of the collapse depending on the viscosity and the density of the surrounding fluid, i.e. the Stokes number and the fluid-grain density ratio. In the triphasic case, i.e. when the granular-fluid column collapses in air, the dynamics may be strongly affected by capillary effects through the Bond number which controls the initial dynamics of the column. When these effects can be neglected (large Bond number) at the column and grain scales, the interstitial fluid can have a driven or a dissipative role on the granular media leading to a runout length more or less extended in comparison to the dry case. The role of the interstitial fluid depends mainly on its viscosity which also modifies the collapse dynamics. Finally, a preliminary study is realized on the flow dynamics of an initially homogeneous negative buoyant suspension column. This case, which makes the transition between the above mentioned configurations, allows to study the coupling between the settling dynamics of particles and this of the current. In particular, we observe that the settling velocity decreases with the increase of the initial volume fraction of particles.

Le but de cette étude est d’étudier l’utilisation d’interactions fluide-structure (FSI) pour améliorer le transfert de chaleur et les performances de mélange dans des échangeurs-réacteurs multifonctionnels, et d’évaluer des configurations pour lesquelles l’objectif est de produire et de maintenir un régime dynamique auto-entretenu d’oscillations des générateurs de tourbillons flexibles. Dans un premier temps, deux études numériques ont été réalisées pour des écoulements laminaires bidimensionnels. Les résultats montrent qu’un minimum de trois générateurs de tourbillons alternés est nécessaire pour produire une instabilité qui engendre les oscillations de larges amplitudes. L’ajout de deux promoteurs coplanaires en amont déstabilise l’écoulement en créant des forces périodiques agissant sur les générateurs de tourbillons en aval. Il en résulte une augmentation de la vitesse réduite qui impose un blocage en fréquence des oscillations des générateurs de tourbillons en aval. Dans cette configuration, des oscillations de larges amplitudes sont obtenues pour uniquement deux générateurs de tourbillons en aval. Les oscillations des générateurs de tourbillons produisent une vorticité intense qui a une incidence positive que le transfert de chaleur et sur le mélange. Dans un second temps, une configuration tridimensionnelle HEV incluant des générateurs de tourbillons trapézoïdaux flexibles orientés a 45◦ vers l’amont est étudiée par simulations numériques. Une analyse FFT réalisée sur les coefficients issus d’une analyse POD montre un pic fréquentiel correspondant aux formations et lâchers tourbillonnaires périodiques. Cette fréquence dominante correspond bien au mode propre d’oscillation des générateurs de tourbillons et engendre ainsi de larges amplitudes d’oscillations. / The aim of this study is to investigate the use of fluid-structure interaction (FSI) to improve heat transfer and mixing performances in multi-functional heat exchangers/reactors, and to evaluate configuration designs where the main target is to produce and maintain self-sustained oscillations of flexible vortex generators. At first, two dimensional laminar flow studies are numerically investigated. The results show that a minimum of three alternating flaps is needed to produce an instability that leads to large displacement oscillations. However, the introduction of two co-planar flaps upstream destabilizes the flow by creating periodic forces that act on the alternating downstream flaps. Hence, this results in artificially increasing the reduced velocity that will induce the alternating flaps to be in a lock-in state. Thus in this case, large displacement amplitudes are created with two alternating flaps only. The free flaps oscillations produce vortices of higher strength which have a positive impact on heat transfer and mixing. Secondly, a three dimensional HEV configuration with flexible trapezoidal vortex generators inclined with an angle of 45◦ with respect to the wall and reversed opposite to the flow direction is numerically investigated. Fast Fourier Transformation is applied on the temporal variation of the Proper Orthogonal Decomposition (POD) coefficientswhich displays a dominant peak in the flow and corresponds to the vortices periodic formation and detachment. This dominant frequency synchronizes well with the structural oscillation frequency and the fundamental frequency of the tabs reaching a lock-in state and leading to large oscillation amplitudes.

Dans le domaine pétrolier, l'exploitation de pétroles bruts de plus en plus lourds nécessite de développer des procédés de conversion de ces coupes lourdes en bases carburants valorisables. Parmi ces procédés, l'hydrocraquage permet d'obtenir à partir d'une coupe distillat sous vide (DSV) des gazoles de grande qualité. Afin d'améliorer la compréhension de ce procédé catalytique, une caractérisation détaillée des charges et effluents de ce procédé est nécessaire. Les techniques existantes n'étant pas assez performantes, ce travail de thèse s'est focalisé sur l'utilisation de la chromatographie en phase gazeuse bidimensionnelle haute température (GC×GC-HT) et de la spectrométrie de masse haute résolution à transformée de Fourier (FT-ICR/MS). Il a ainsi été montré que la GC×GC-HT pouvait être utilisée pour éluer des composés ayant des points d'ébullition équivalents au nC68 (641 °C). Ce travail a également permis de proposer une méthode de caractérisation des composés soufrés par familles en GC×GC-HT-SCD dont les résultats sont cohérents avec les données obtenues par FT-MS. Pour les composés azotés, les limites de la GC×GC en terme de séparation ont été atteintes malgré la mise en place d'une pré-séparation en ligne par SFC. Une méthodologie de quantification par famille à partir de la FT-MS a alors été proposée pour étudier sélectivement les composés azotés basiques (mode ESI+) ou neutres (mode ESI-). Ce travail a permis de proposer de nouveaux outils analytiques de caractérisation des DSV et de confirmer la nature chimique des familles les plus réfractaires à l'hydrotraitement, à savoir les composés fortement alkylés de type dibenzothiophène et carbazole. / Converting heavy petroleum cuts into valuable fuels becomes a strong necessity for the refining industry as crude oil quality globally decreases. Among all the conversion processes, vacuum gas oil (VGO) hydrocracking is certainly the most suitable process to convert VGOs into high quality diesel fuel.In order to improve the understanding of this catalytic process, a detailed characterization of VGO feeds and products is needed. As existing analytical methods are not currently powerful enough, this work focused on the use of high-temperature comprehensive bidimensional gas chromatography (HT-GC×GC) and Fourier-transform high resolution mass spectrometry (FT-ICR/MS). First, it was shown that HT-GC×GC can elute heavy boiling point compounds up to nC68 (641 °C). Moreover, a HT-GC×GC-SCD method was developed to characterize sulfurcontaining compounds according to their chemical family. The obtained results were in accordance with parallel FT-ICR/MS measurements for sulfur compounds. For nitrogen-containing compounds, separation limits of HTGC×GC-NCD were reached despite the use of an additional online pre-separation by supercritical fluid chromatography (SFC). Therefore, a FT-MS based methodology was proposed to quantify nitrogen-containing chemical families and study basic and neutral compounds with ESI(+) and ESI(-) respectively.Hence, this work proposed innovative analytical tools for the characterization of vacuum gas oils. It also confirmed that highly-alkylated dibenzothiophenes and carbazoles were the most refractory sulfur and nitrogen-containing compounds towards hydrotreatment.

La méthode IIM (Immersed Interface Method) permet d'étendre certaines méthodes numériques à des problèmes présentant des discontinuités. Elle est utilisée ici pour étudier un fluide incompressible régi par les équations de Navier-Stokes, dans lequel est immergée une membrane exerçant une force singulière. Nous utilisons une méthode de projection dans une grille de différences finies de type MAC. Une dérivation très complète des conditions de saut dans le cas où la viscosité est continue est présentée en annexe. Deux exemples numériques sont présentés : l'un sans membrane, et l'un où la membrane est immobile. Le cas général d'une membrane mobile est aussi étudié en profondeur. / The Immersed Interface Method allows us to extend the scope of some numerical methods to discontinuous problems. Here we use it in the case of an incompressible fluid governed by the Navier-Stokes equations, in which a membrane is immersed, inducing a singular force. We use a projection method and staggered (MAC-type) finite difference approximations. A very complete derivation for the jump conditions is presented in the Appendix, for the case where the viscosity is continuous. Two numerical examples are shown : one without a membrane, and the other where the membrane is motionless. The general case of a moving membrane is also thoroughly studied.

La résolution analytique du problème de la stabilité d’une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s’est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l’écoulement dans le cadre d’une théorie linéaire. Ces équations sont valables dans l’approximation du faible espace annulaire et ont pour but la prévision de l’instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d’onde au voisinage de l’état critique du développement de la première instabilité. Tout d’abord, on a cherché à mettre au point les calculs dans le cas de l’écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l’on cherche à déterminer le critère d’apparition des instabilités en présence de particules en suspension et l’on détermine simultanément les paramètres de couplage entre forces d’interaction liquide-solide. L’ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l’expérience. L’étude expérimentale a permis d’analyser les effets de la concentration des particules en suspension (disques) et du rapport d’aspect radial ’ sur l’apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d’un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L’apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l’intervalle étudié, la nature des structures apparaissant dans le système d’écoulement dépend ’, alors que pour une valeur donnée de ’ dans l’intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de ’. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n’affectent pas le type de structures qui apparaissent dans le système d’écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d’apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles. / The analytical solution of the stability problem of a solid suspension (rigid spherical particles) in the system of cylindrical Taylor-Couette was conducted. We are based on the work of Ali and Lueptow (2002) to formulate the equations governing the stability of the flow in a linear theory. These equations are valid in the approximation of small gap configuration and aim to predict the primary instability. For this purpose, we used a variational method such as the Galerkin method to solve the eigenvalue problem leading to establish the stability diagram related to the wave number in the vicinity of the critical state of development of the first instability. First, we develop the calculations in the case of Taylor-Couette flow with reference to classic work of Ali and Lueptow (2002). Then, we carried out a systematic solution of the equations of motion and we search to determine the criterion of onset of instabilities in the presence of suspended particles and coupling parameters are simultaneously determined from liquid-solid interaction force. All work carried out and will remove the fundamental contradiction between theory and experiment. The experimental study has analyzed the effect of the concentration of suspended particles (disks) and radial aspect ratio ' on the occurrence of instabilities in the Taylor-Dean flow system. The experimental device used consists of a rotating inner cylinder and the outer cylinder is stationary. The rheological behavior of the fluid is viscoplastic obeying to Herschel Bulkley model. The onset of instability is examined using a visualization technique. For a given concentration in the range studied, the nature of the structure appearing in the flow system depends on ', while for a given value of ' in the range studied, the value of the critical Taylor number Tac depends on the particle concentration. We obtain that the critical Taylor number Tac corresponding to the onset of the first instability evolves nonlinearly versus '. In addition, we examined experimentally the effect of axial limitation (endwall effects) on the onset of instabilities in the Taylor-Dean flow system. The results show that the rotating endwalls do not affect the type of structures that appear in the Taylor-Dean flow system. However, they influence the threshold of appearance of instabilities which is characterized by high values of the critical Taylor number for rotating endwalls indicating a stabilizing effect of the rotating endwalls.