Dissertations / Theses on the topic 'Modèle eulerien à deux fluides'

Ce travail vise à l'amélioration des modèles eulériens à deux fluides pour l'hydrodynamique et le transfert de masse en écoulement turbulent à bulles. Nous étudions des zones de mélange oxygène-eau pour des taux de vide 2<br>In this work we develop a Eulerian two-fluids model for hydrodynamics and mass transfer in turbulent bubbly flows. We study oxygen-water mixing layers for void fractions 2

Un modele a deux fluides applique a la simulation des lits fluidises denses (des milieux diphasiques gaz-particules qui font intervenir une large gamme de fraction volumique de solide (0-70 %)) est valide d'une part microscopiquement en etudiant le mouvement de particules dans ces milieux, d'autre part macroscopiquement en comparant les resultats du modele avec des resultats experimentaux. Le couplage entre les phases intervient sur les equations de transport de la quantite de mouvement et des energies fluctuantes de chaque phase. La modelisation s'appuie non seulement sur la theorie cinetique des milieux granulaires dans laquelle l'action du fluide a ete incluse de facon a obtenir les expressions correctes des coefficients de transport ; mais egalement sur une modelisation des interactions hydrodynamiques entre particules a grand nombre de stokes. La validation microscopique de la modelisation a ete effectuee a l'aide de simulations numeriques lagrangiennes considerees comme des experiences numeriques. Le premier cas de validation utilise un ecoulement gaz-particules en cisaillement simple, il permet de valider la competition entre les mecanismes dissipatifs issus et des chocs entre particules et de la trainee. Le second cas de validation fait intervenir un milieu dense de particules en sedimentation dans un ecoulement a grand nombre de stokes, et permet de valider la modelisation des interactions hydrodynamiques. Ce dernier cas nous a conduit a developper une methode originale de simulation lagrangienne avec couplage reciproque entre le fluide et les particules. La validation macroscopique utilise les resultats de la simulation eulerienne des lits fluidises denses. La hauteur du lit fluidise dense, la circulation des particules, et les caracteristiques des bulles qui apparaissent spontanement ont ete etudiees et comparees a des mesures experimentales en regardant l'influence des parametres physiques et numeriques.

Nous avons mis au point un système expérimental modèle d'adhésion spécifique entre deux surfaces fluides sur lesquelles diffusent des protéines de type ligand-récepteur. Nous avons notamment voulu relier la mesure d'une énergie d'adhésion macroscopique à son moteur microscopique, la formation de liens spécifiques. Pour mesurer une énergie macroscopique d'adhésion, nous avons utilisé des gouttelettes d'émulsion. Lorsqu'elles adhèrent, les gouttelettes se déforment. Le gain de surface permet d'accéder à l'énergie d'adhésion. Comme moteur de l'adhésion spécifique, nous avons utilisé le couple modèle streptavidine/biotine. Les gouttelettes d'émulsion sont fonctionnalisées avec ces protéines ligantes. Amenées par gravité en présence de bicouches lipidiques supportées fluides dont une fraction des lipides porte les protéines réceptrices, les gouttelettes adhèrent spécifiquement aux bicouches lipidiques. Nous avons mené une étude statique de cette adhésion. Grâce à différentes techniques de fluorescence (épifluorescence, FCS (Fluorescence Correlation Spectroscopy), dosages en retour) nous avons quantifié le nombre de protéines mises en jeu dans nos patches adhésifs afin de le comparer à l'énergie de déformation de la goutte. Nous avons pu observer que l'énergie d'adhésion est proportionnelle à l'énergie d'un lien. Par ailleurs, des mesures en RICM (Reflection Iterference Contrast Microscopy) nous ont permis d'étudier la dynamique de croissance du patch adhésif. Nous avons pu mettre en évidence une croissance diffusive et quasistatique de la zone de contact. En jouant sur les quantités de protéines mises en jeu ainsi que sur leur nature (longueur des liens), nous avons exploré les différentes situations statiques et dynamiques de l'adhésion dans notre système modèle. Nous avons notamment mis en évidence les effets de compressibilité des films adhésifs qui apportent de nouvelles clés pour envisager l'adhésion cellulaire.

La simulation de réacteurs à bulles en régime industriel est un grand défi. L'objectif principal de ce travail est la prédiction de la taille des bulles à l’aide d’un modèle numérique de bilan de population, basé sur la modélisation des phénomènes de brisure et de coalescence, et pouvant être couplé aux conditions hydrodynamiques présentes dans les réacteurs. Différentes données expérimentales sont obtenues pour valider le modèle. La taille des bulles est mesurée à l'aide d'une technique innovante de corrélation croisée. Les essais, réalisés en eau du réseau (partiellement contaminée) et en eau déminéralisée avec ajout éventuel d'éthanol, montrent que les additifs réduisent la coalescence et diminuent la taille moyenne des bulles. Deux distributeurs du gaz différents sont utilisés pour découpler l'étude de la brisure et de la coalescence. Les données expérimentales sont utilisées initialement pour valider des simulations CFD 3D transitoires Eulériennes-Eulériennes. La loi de traînée est corrigée par un facteur de swarm pour intégrer l’effet d’une fraction de gaz élevée. Différents modèles de turbulence sont testés. La contribution de la turbulence induite par les sillages de bulles au mélange de scalaires est évaluée. Enfin, pour prédire la taille des bulles, un bilan de population est couplé au modèle hydrodynamique préalablement validé et est résolu par la méthode de quadrature des moments (QMOM). Un set original de kernels de brisure et coalescence est proposé, capable de prédire la taille des bulles pour différentes conditions opératoires. Le comportement du modèle lors de l’extrapolation des réacteurs est également examiné<br>The simulation of bubble column reactors under industrial operating conditions is an exciting challenge. The main objective of this work is to predict the bubble size, in turn interconnected to the reactor hydrodynamic conditions, with computational models, by modelling bubble breakage and coalescence. Experimental data is collected for model validation, including bubble size measurements with an innovative cross-correlation technique. Experiments are carried out with tap water and demineralized water, with or without the addition of ethanol, and gathered results show that additives reduce coalescence and lower the mean bubble size. Two different spargers are used, in order to decouple the investigation of breakage and coalescence. The experimental data set is used to validate out unsteady three-dimensional Eulerian-Eulerian CFD simulations. A drag law for oblate bubbles is considered, together with a swarm factor, that accounts for the swarm effect. Several turbulence models are tested. The contribution of bubble induced turbulence (BIT) to scalar mixing is assessed. To predict bubble size, a population balance model is coupled to the hydrodynamic model and is solved with the quadrature method of moments. A set of breakage and coalescence kernels is proposed, capable of predicting the bubble size for different operating conditions. Scale-up effects are also investigated

Les milieux diphasiques mettent en jeu des interactions aux interfaces qui modifient d’une manière significative la structure des champs moyens et fluctuants des écoulements. La conception de modèles à deux fluides adaptés aux écoulements industriels nécessite la prise en compte de l’effet de ces interactions au niveau des relations de fermeture adoptées. Le travail développé dans cette thèse porte sur le développement de modèles à deux fluides au premier ordre déduits par réduction des fermetures au second ordre. La démarche adoptée, basée sur le principe de décomposition du tenseur des contraintes de Reynolds en contributions turbulente et pseudo-turbulente statistiquement indépendantes, permet de préserver le contenu physique des fermetures au second ordre. L’analyse de la structure de la turbulence dans deux configurations d’écoulements de référence : écoulements homogènes à bulles uniforme derrière une grille et à cisaillement constant permet de déduire une formulation de la viscosité turbulente en milieux diphasiques faisant intervenir les échelles caractéristiques de la turbulence à bulles, ainsi qu’une description analytique de la modification induite par la présence des bulles de la structure de la turbulence homogène. Le modèle eulérien à deux fluides a été généralisé ensuite au cas des écoulements inhomogènes à faibles taux de vide. Les résultats numériques obtenus par l’application de ce modèle intégré dans le code de calcul MELODIF dans le cas de l’écoulement turbulent à bulles cisaillé libre de sillage ont montré une concordance satisfaisante avec les données expérimentales et ont permis d’analyser la modification des échelles caractéristiques de cet écoulement par les interactions interfaciales. Le modèle à deux fluides au premier ordre est généralisé enfin au cas des écoulements à forts taux de vide où les interactions hydrodynamiques entre les bulles ne sont plus négligeables.

Les milieux diphasiques mettent en jeu des interactions aux interfaces qui modifient d'une manière significative la structure des champs moyens et fluctuants des écoulements. La conception de modèles à deux fluides adaptés aux écoulements industriels nécessite la prise en compte de l'effet de ces interactions au niveau des relations de fermeture adoptées. Le travail développé dans cette thèse porte sur le développement de modèles à deux fluides au premier ordre déduits par réduction des fermetures au second ordre. La démarche adoptée, basée sur le principe de décomposition du tenseur des contraintes de Reynolds en contributions turbulente et pseudo-turbulente statistiquement indépendantes, permet de préserver le contenu physique des fermetures au second ordre. L'analyse de la structure de la turbulence dans deux configurations d'écoulements de référence : écoulements homogènes à bulles uniforme derrière une grille et à cisaillement constant permet de déduire une formulation de la viscosité turbulente en milieux diphasiques faisant intervenir les échelles caractéristiques de la turbulence à bulles, ainsi qu'une description analytique de la modification induite par la présence des bulles de la structure de la turbulence homogène. Le modèle eulérien à deux fluides a été généralisé ensuite au cas des écoulements inhomogènes à faibles taux de vide. Les résultats numériques obtenus par l'application de ce modèle intégré dans le code de calcul MELODIF dans le cas de l'écoulement turbulent à bulles cisaillé libre de sillage ont montré une concordance satisfaisante avec les données expérimentales et ont permis d'analyser la modification des échelles caractéristiques de cet écoulement par les interactions interfaciales. Le modèle à deux fluides au premier ordre est généralisé enfin au cas des écoulements à forts taux de vide où les interactions hydrodynamiques entre les bulles ne sont plus négligeables<br>The work developed in this thesis concerns the development of first order two-fluid models deduced by reduction of second order closures. The adopted reasoning allows to preserve the physical contents of the second order relations closure. Analysis of the turbulence structure in homogeneous uniformly sheared bubbly flow allows to deduce a formulation of the two-phase turbulent viscosity. The eulerian two-fluid model was then generalized with the case of the inhomogeneous flows with low void fractions. The numerical results obtained by the application of this model integrated in the computer code MELODIF in the case of free sheared turbulent bubbly flow of wake made it possible to analyze the modification of the characteristic scales of such flow. The two-fluid first order model is generalized finally with the case of high void fractions bubbly flows where the hydrodynamic interactions between the bubbles are not negligible any more

L'objet de cette étude est d'utiliser la décomposition en modes propres (POD) pour établir un modèle de la composante déterministe et aléatoire d'un écoulement turbulent présentant une instationnarité à grande échelle pseudo périodique.<br />Une étude expérimentale de l'écoulement bidimensionnel en moyenne en aval d'un cylindre semi-circulaire, par vélocimétrie par image de particules en deux temps (PIV2T) caractérise l'écoulement<br />Une analyse POD du champ de vitesse permet d'extraire les modes spatiaux et de définir un paramètre de phase décrivant l'instationnarité à grande échelle qui régit la partie déterministe. La modélisation de l'évolution temporelle des coefficients associés aux modes s'effectue par des fonctions soit harmoniques pour la partie déterministe, soit stochastiques pour la partie aléatoire.<br />La modélisation est en accord avec les mesures expérimentales des premiers moments statistiques en un point et des fonctions de corrélation spatio-temporelle du champ de vitesse.

Le but de cette thèse est de proposer une extension de la modélisation statistique edqnm (ou modele quasi-normal markovianise à amortissement tourbillonnaire) au cas d'une turbulence inhomogène, donc d'élargir le champ d'application des fermetures en deux points, encore appelées modèles spectraux, à des écoulements plus réalistes, se rapprochant plus de ceux rencontrés et etudiés par un ingénieur. Ce travail se situe dans la suite logique des travaux amorcés au laboratoire de mécanique des fluides et d'acoustique, URA CNRS 263, pour aboutir à une description spectrale de la turbulence inhomogène. Des comparaisons numériques (avec des calculs de simulation des grandes échelles et avec le modèle scit) sont effectuées dans un cas test très simple qui est la couche de mélange non cisaillée. Il s'agit de l'étude de l'interaction de deux champs turbulents situés c^cte a côte, soumis à aucun gradient de vitesse moyenne, et possédant des échelles caractéristiques et des niveaux d'énergie différents. Nous nous sommes également intéressés aux implications que notre modèle pouvaient avoir sur les modèles en un point. De nouvelles modélisations sont en particulier proposées à partir d'une analyse spectrale de nos équations, puis comparées aux formulations classiques<br>The aim of the present work is to extend a two-point closure theory (namely the Eddy Damped Quasi Normal Markovian theory) to inhomogeneous turbulence and to propose a model that may be applied to more realistic flows than the usual two-point theories (usually restricted to homogeneous turbulence). The present study follows the earlier works carried out at the "Laboratoire de Mécanique des Fluides et d'Acoustique" of the Ecole Centrale de Lyon, in order to device two-point models for inhomogeneous turbulence. Numerical comparisons are carried with Large Eddy Simulations and the SCIT model, in the simple case of the shearless turbulence mixing layer. The adopted configuration consists in the interaction of two nearly homogeneous turbulences, situated side by side and having different energies and length scales. Spectral analysis of the present model also leads to propose improved one-point formulations that are compared to classical ones

Notre étude porte sur la caractérisation d’un jet d’air turbulent chargé de particules ou de gouttelettes. Elle se place dans le cadre d’une coopération entre le CEMAGREF et des partenaires industriels. Le but ultime de l’étude est le transfert vers ces partenaires industriels du maximum de connaissances utilisables pour le développement et l’amélioration de ses produits. En ce sens, l’étude se veut la plus complète possible. Nous avons également le souci de produire des résultats applicables et d’utiliser et de développer des techniques d’étude susceptibles d’être adoptées par les partenaires à l’issue de notre collaboration. Notre étude porte sur la dynamique pure de l’écoulement. De plus, nous nous intéressons à un modèle réduit et simplifié des produits de nos partenaires qui offre des débits, et donc une portée, plus faibles, ce qui permet le montage d’une installation expérimentale en laboratoire où on contrôle parfaitement l’environnement dans lequel le jet fonctionne. Ce mémoire consiste en trois parties. La première présente les bases fondamentales des phénomènes qui entrent en jeu dans les écoulements turbulents diphasiques à phase dispersée, puis dresse un état de l’art succinct des connaissances et moyens d’étude pour ce type d’écoulements. La seconde partie détaille notre étude expérimentale, les techniques employées et les résultats obtenus. Une troisième partie présente le travail effectué sur le développement d’un outil numérique de prédiction de l’écoulement, notre base de départ, les modèles passés en revue pour les différents phénomènes et les premiers résultats obtenus<br>We present here the study we have conducted on a gaseous turbulent jet loaded with solid particles or liquid droplets. This study is performed in the framework of the collaboration between the Cemagref and industrial partners. The final aim of our work is the transfer to our industrial partners of a maximum of useful knowledge for the development and improvement of their products. Therefore we have tried to perform the most complete study. We also had the commitment to produce utilizable results and to use and develop study techniques likely to be adopted by the partners at the end of our collaboration. Our study focuses on the dynamic side of the flow. Furthermore we have built a reduced-scale model with smaller flow rates and therefore smaller flow range. The use of this scaled model enabled us to design an indoor laboratory experimental facility where we control exactly the conditions in which the flow operates. This report consists of three parts. The first one presents the fundamentals for the phenomena that govern dispersed two-phase turbulent flows. Then we review the state of the art for the knowledge and study techniques for such flows. In the second part we detail the experimental part of our study, the techniques we used and the results we have obtained. Then the third part presents the work performed for the development of a numerical prediction tool, where from we started, the models we reviewed to account for the physics involved in the flows and give an outline of the first results

Les écoulements polyphasiques sont constitués de différentes tailles d'inclusions et de structures turbulentes. Les formalismes classiques diphasiques sont généralement dédiés à une gamme restreinte de tailles d'interface. Les travaux entrepris ici proposent un formalisme à deux fluides compressible permettant de prendre en compte différentes tailles d'interface, appellé Simulation des Grandes Interfaces (SGI). Une opération de filtrage a été réalisée, s'inspirant de la LES monophasique, afin de simuler les grandes interfaces (GI) et de modéliser les petites (PI). Le terme de force capillaire de GI apparaissant dans les équations de quantité de mouvement est modélisé à l'aide de la méthode CSF. Les termes interfaciaux associés aux petites inclusions sont modélisés par les lois de fermeture classique utilisées par l'approche moyennée. Pour calculer les termes GI uniquement à la présence d'une GI, un algorithme de reconnaissance des GI est développé à partir d'un critère adimensionnel. Le modèle SGI est appliqué à l'aide d'un code diphasique moyenné. La reconnaissance des GI et la modélisation des forces capillaires de GI sont validées sur des cas analytiques et expérimentaux. Le modèle SGI est appliqué à la fragmentation d'une bulle soumise à un champ liquide turbulent; à la fragmentation d'une bulle (passage de GI et de PI) soumise à la gravité et traversant une grille. Puis impactant sur un plafond (passage de PI à GI)<br>Different ranges of size of interfaces and eddies are involved in multiphase flow phenomena. Classical formalisms focus on a specifie range of size. This study presents a Large Interface Simulation (LIS) two-fluid compressible formalism taking into account different sizes of interfaces. As in the single-phase Large Eddy Simulation, a filtering process is used to point out Large Interface (LI) simulation and Smallinterface (SI) modelisation. The LI surface tension force is modelled adaptating the well-known CSF method. The modelling of SI transfer terms is done calling for classical closure laws of the averaged approach. To simulate accurately LI transfer terms, we develop a LI recognition algorithm based on a dimensionless criterion. The LIS model is applied in a classical averaged two-fluid code. The LI transfer terms modelling and the LI recognition are validated on analytical and experimental tests. A square base basin excited by a horizontal periodic movement is studied with the LIS mode!. The capability of the model is also shown on the case of the break-up of a bubble in a turbulent liquid flow. The break-up of a large bubble at a gridimpact performed regime transition between two different scales of interface from LI to SI and from PI to LI

L'objectif de cette thèse est d'étudier différents aspects microscopiques et macroscopiques liés à la présence de superfluidité dans les étoiles à neutrons. Dans un premier temps, nous avons calculé des configurations stationnaires d'étoiles à neutrons superfluides en rotation, en relativité générale, basées sur l'utilisation d'équations d'état réalistes. A l'aide de ces configurations d'équilibre, nous avons ensuite développé un modèle simple de glitch, en relativité générale, vu comme un transfert de moment cinétique entre les neutrons superfluides et les particules chargées constituant l'étoile. Cela nous a permis d'obtenir des temps caractéristiques de montée qui pourront être comparés à de futures observations précises de glitches afin d'apporter de meilleures contraintes sur l'intérieur de ces étoiles. Enfin, nous nous sommes également intéressés à la dynamique des vortex superfluides, en présence de tubes de flux, dans le cas où les protons dans le coeur des étoiles formeraient un supraconducteur de type II<br>The aim of this thesis is to study different aspects, both microscopic and macroscopic, associated with the presence of a large amount of superfluid matter inside neutron stars. First, we computed stationary configurations of rotating superfluid neutron stars, in general relativity, using realistic equations of state. Based on these equilibrium configurations, we then developed a simple model of pulsar glitches, in general relativity, seen as angular momentum transfers between the superfluid neutrons and the charged particles composing the star. This enables us to infer spin-up time scales that could be compared with future accurate glitch observations, in order to get some constraints on the interior of neutron stars. Finally, we also focused on the dynamics of superfluid vortex lines, accounting for the presence of fluxtubes, if the protons are forming a type II superconductor in the core of neutron stars

Ce travail traite des aspects macroscopiques et microscopiques de l'écorce interne d'une étoile à neutrons, formée d'un solide de noyaux plongé dans un superfluide de neutrons. Une première partie expose une formulation quadridimensionnelle covariante de l'hydrodynamique non relativiste d'un mélange de fluides parfaits, basée sur un principe variationnel convectif. Ce formalisme est appliqué à la description de l'écorce, comme un mélange de deux fluides, un superfluide de neutrons et un plasma de noyaux et d'électrons, couplés par un entraînement non dissipatif. La seconde partie est dédiée à l'étude microscopique de cet entraînement.<br />Appliquant des méthodes de champ moyen au-delà de l'approximation de Wigner-Seitz, nous montrons que cet entraînement résulte de la diffraction de Bragg des neutrons libres sur les noyaux. Celle-ci se traduit par une masse de neutron effective "mésoscopique", qui, contrairement à la masse effective "microscopique", est très grande devant la masse "nue", dans les couches intermédiaires.

Cette travail de thèse est lié au besoin de modélisation des écoulements diphasiques en générateurs de vapeur (entrée liquide et sortie vapeur). La démarche proposée consiste à faire le choix d'une modélisation hybride de l'écoulement, en scindant la phase gaz en deux champs, modélisés de manières différentes. Ainsi, les petites bulles sphériques sont modélisées avec une approche dispersée classique avec le modèle eulérien à deux fluides, et les bulles déformées sont simulées à l'aide d'une méthode de localisation d'interface. Le travail effectué porte sur la mise en place, la vérification et la validation du modèle dédié aux larges bulles déformées, ainsi que le couplage entre les deux approches pour le gaz gaz, permettant des premiers calculs de démonstration utilisant l'approche hybride complète.

Dans les réacteurs nucléaires, divers régimes d’écoulement gaz-liquide peuvent apparaître avec des transitions entre eux. La modélisation de ces transitions dans les codes CFD 3D requière le traitement d’interfaces déformables de différentes tailles, la prise en compte d’interactions par coalescence et fragmentation ainsi que le développement de lois de fermeture indépendantes du régime. Ce travail vise la modélisation et la simulation de l’hydrodynamique des écoulements gaz-liquide adiabatiques grâce à un modèle bi-fluide à trois champs dans Neptune CFD. Dans une première étape, un modèle avec un champ liquide continu et deux champs de gaz dispersés représentant petites et grandes bulles est utilisé pour simuler des écoulements cap et churn avec un taux de vide jusqu’à 0.5 et une emphase est mise sur la prédiction de l’aire interfaciale. Dans une seconde étape, le deuxième champ dispersé est remplacé par un champ hybride continu/dispersé représentant les grandes bulles et les régions continues de gaz. Le modèle est validé sur plusieurs régimes d’écoulements en tuyaux de large diamètre et dans des canaux rectangulaires confinés<br>In nuclear reactors, several regimes of gas-liquid flows may occur with some transitions between them. The main challenges associated with simulating these transitions in 3D CFD codes are associated with deformable interfaces of different sizes, accounting for coalescence and breakup interactions between gas structures and developing flow regime independent closure relations. This work aims at modelling and simulating the hydrodynamics of adiabatic gas-liquid flows thanks to a three-field two-fluid model in Neptune CFD. In a first step, a model with one continuous liquid field and two dispersed gas fields for small and large bubbles is used to simulate cap and churn flows with a void fraction up to 0.5 and a focus is put on the interfacial area prediction. In a second step, the second dispersed field is replaced by a hybrid continuous/dispersed field representing both large bubbles and continuous gas regions. The model is validated on several flow regimes in large diameter pipes and in confined rectangular channels

Motivées par des données expérimentales sur la magnétisation de réseau de nanofils de plomb, les résolutions numériques des équations stationnaires de Ginzburg-Landau (GL) se sont focalisées sur les géométries à symétrie axiale. L'effet Meissner, les états représentant un vortex d'Abrikosov ou encore des Vortex Géants (``GiantVortex') centrés à l'origine du cylindre ont alors pu être identifiés sous l’hypothèse d’invariance sous rotation selon l’axe de symétrie du cylindre étudié (modèle à une dimension, 1D). En identifiant le type de transition par le caractère continu ou non du paramètre d'ordre autour du changement de phase, une frontière à l'échelle mésoscopique a également pu être identifiée au travers du modèle 1D. Plus spécifiquement, la limite entre les deux types de transitions décrite par le paramètre phénoménologique κ = λ /ξ ( =1/√2 à l’échelle macroscopique) devient une fonction non constante dépendant à la fois du rayon normalisé, u=R/λ, et de la vorticité L: κ =f(u,L). Les deux longueurs caractéristiques λ et ξ représentent respectivement les longueurs de pénétration et de cohérence d’un échantillon supraconducteur. Une comparaison avec les résultats obtenus par Zharkov permet de valider notre démarche numérique employée pour la résolution numérique des équations de GL à une dimension. En employant un modèle à deux dimensions (2D), la symétrie sous rotation des solutions a également été relâchée. Basée sur le principe de moindre action, la résolution propose alors un schéma numérique indépendant du type d'équations du mouvement à solutionner. Les configurations du type MultiVortex ont alors pu être identifiées, et comparées aux solutions du groupe du Professeur F. Peeters. Ces différents accords ont confirmé la démarche développée. Une modélisation de la magnétisation expérimentale d'un réseau de nanofils a également été développée. De par la taille réduite des nanofils, l'interaction magnétique entre ceux-ci a pu être négligée. La magnétisation totale du réseau est alors construite par une sommation incluant la contribution individuelle en magnétisation de chaque fil, pondérée par un poids reflétant une distribution gaussienne pour les rayons des fils constituant le réseau. La magnétisation individuelle est évidemment obtenue par résolution des équations du mouvement de GL précédemment étudiées avec les modèles 1D et 2D. En ajustant les paramètres libres associés à ce modèle décrivant la magnétisation totale du réseau, les données expérimentales ont pu être reproduites endéans 10% de marge d'erreur, l'intervalle d'incertitude caractéristique de la théorie effective de Ginzburg-Landau. Ces variables attachées au modèle de la magnétisation totale, reprennent la valeur moyenne m et l'écart-type s de la distribution gaussienne, ainsi que les longueurs caractéristiques λ(T) et ξ(T) présentes dans la théorie de GL. Un test totalement indépendant de l'analyse des magnétisations a permis de valider les valeurs déterminées pour la distribution des rayons. Les grandeurs ajustées pour les longueurs λ(T) et ξ(T) ont fait l'objet d'une analyse supplémentaire en termes de leur dépendance en température et du libre parcours moyen des électrons. Malgré l'accord entre les données expérimentales et la magnétisation théorique, il est important de mentionner qu'un paramètre libre supplémentaire, associé à l'apparition de configurations décrivant un vortex magnétique, a dû être introduit. Il modifie empiriquement la métastabilité trop longue en mode champ externe décroissant de l'état décrivant un vortex d'Abrikosov. La correction expulse donc le vortex avant sa prédiction théorique liée à la disparition de la barrière de Bean-Linvingston. Une étude plus approfondie de cette barrière de potentiel fut donc également réalisée. Cependant, elle n'est pas concluante en regard des données expérimentales analysées. Il n'en demeure pas moins que la transition apparaît dans un domaine en champ magnétique cohérent vis-à-vis de la description en énergie libre des états de vorticités voisines d'une unité de quantum de flux magnétique. La correspondance entre les longueurs caractéristiques du modèle phénoménologique de GL et les longueurs issues des théories microscopiques de Pippard et BCS a également abordée. Cette étude permet entre autre de comparer les différentes dépendances possibles en température avec les longueurs obtenues de l'analyse de magnétisation des nanofils en plomb. Au delà de l'accord avec le modèle des deux-fluides de Gorter et Casimir, une extrapolation bien en deçà de la température critique Tc est proposée pour les paramètres phénoménologiques λ(T) et ξ(T) de Ginzburg-Landau. Même si la correspondance entre les magnétisations expérimentales et théoriques semblait déjà l'indiquer, il est possible d'appliquer les équations de Ginzburg-Landau pour décrire le comportement magnétique du plomb bien en deçà de sa température critique. De plus, les paramètres associés possèdent une dépendance tout à fait conforme à une autre théorie empirique, le modèle des deux-fluides. Basée sur le modèle de Pippard, une détermination de la valeur du libre parcours moyen des normaux a également été isolée. Elle justifie alors une distinction entre les deux échantillons analysés en terme de leur degré d'impureté. Les résultats électrons obtenus étant en accord avec les procédures de fabrication des nanofils de plomb, cette nouvelle constatation, positive avec l'expérience, confirme une fois de plus la cohérence du modèle développé pour la magnétisation totale, et justifie l'emploi des équations de GL à toutes les températures en dessous de Tc. / Mesoscopic superconductors are described within the framework of the nonlinear Ginzburg-Landau theory. The two coupled nonlinear equations are solved numerically and we investigate the properties, in particular the order of the transition and the vortex configurations, of cylinders submitted to an external magnetic field. Meissner state, Abrikosov vortices, GiantVortex and MultiVortex solutions are described. The Bean-Livingston barrier in mesoscopic cylinders is also numerically studied. This theoretical work was applied to understand experimental magnetizations of lead nanowires in an array well below the superconducting transition temperature Tc. By freely adjusting the GL phenomenological lengths λ (T) and ξ (T), the experimental magnetization curves are reproduced to within a 10% error margin. The Meissner and the Abrikosov state were also experimentally observed in this apparently type-I superconductor. This fact is a consequence of the non-trivial behaviour of the critical boundary κ _c ($=1/√2 in bulk materials) between type-I and type-II phase transition at mesoscopic scales. Beyond the experimental-theoretical agreement, the question whether the GL model remains valid far below Tc is also addressed. The temperature dependence of the adjusted characteristic lengths is compared with different theoretical and empirical laws. The best agreement is achieved for the Gorter-Casimir two-fluid model. A comparison between lead nanowire arrays electrodeposited under constant and pulsed voltage conditions allows us to distinguish both samples in terms of their electronic mean free paths. The characterisation of the latter quantities concurs perfectly with the experimental expectation given the different electrodeposition techniques.

Des analyses de sensibilité et des estimations de paramètres sont étudiées sur deux études de cas de transfert de masse en milieu poreux. La première partie est consacrée à la sensibilité des écoulements souterrains dans une modélisation des échanges drain-aquifère pour mettre en évidence les différences entre les deux méthodes de discrétisation mises en œuvre. La seconde partie est dédiée à la modélisation de l’écoulement en milieu poreux variablement saturé dans une zone humide artificielle, au calage des paramètres du modèle de van Genuchten-Mualem et à l’évaluation de son efficacité à reproduire des données piézométriques collectées sur le site de l’Ostwaldergraben. La variabilité temporelle des paramètres hydrodynamiques, incluant l’effet d’hystérésis, montre que ceux de la couche active du filtre changent au cours du temps. Ces deux études sont conduites à l’aide de la différenciation automatique<br>Sensitivity analyses and parameter estimation are applied to mass transfer in porous media for two remediation facilities. The first part is devoted to the sensitivity analysis of groundwater flows in a modeling of drain-aquifer exchanges to highlight the differences between the two implemented methods of discretization. The second part is dedicated to the modeling of the flow in a variably saturated porous medium in a stormwater constructed wetland, to the calibration of van Genuchten-Mualem parameters and to the evaluation of its efficiency in the reproduction of piezometric data collected on the Ostwaldergraben site. The temporal variability of the hydrodynamic parameters, including the hysteresis effect, shows that the characteristics of the filter layer alters along time. Both studies are carried using automatic differentiation

Les simulations numériques des équations d’Euler deux-fluides réalisé sur des maillages grossiers éliminent les structures fins d’écoulement gaz-solide dans les lits fluidisés. Pour précisément estimer l’hydrodynamique globale de lit, il faut proposer une modélisation qui prend en compte les effets de structure non-résolue. Dans ce but, les maillages sont raffinés pour obtenir le résultat de simulation pleinement résolue ce que les grandeurs statistiques ne modifient plus avec un autre raffinement pour le lit fluidisé périodique dilué gaz-particules sur une géométrie 3D cartésienne et ce résultat est utilisé pour tests "a priori". Les résultats de tests "a priori" montrent que l’équation filtrée de la quantité de mouvement est effectuée mais il faut prendre en compte le flux de la fraction volumique de solide de sous-maille en raison de l’interaction locale de la vitesse du gaz et la fraction volumique de solide pour la force traniée. Nous proposons les modèles fonctionnels et structurels pour le flux de la fraction volumique de solide de sous-maille. En plus, les modèles fermetures du tenseur de sous-maille de la phase dispersée sont similaires aux modèles classiquement utilisés en écoulement turbulent monophasique. Tous les modèles sont validés par test "a priori" et "a posteriori"<br>Eulerian two fluid approach is generally used to simulate gas-solid flows in industrial circulating fluidized beds. Because of limitation of computational resources, simulations of large vessels are usually performed by using too coarse grid. Coarse grid simulations can not resolve fine flow scales which can play an important role in the dynamic behaviour of the beds. In particular, cancelling out the particle segregation effect of small scales leads to an inadequate modelling of the mean interfacial momentum transfer between phases and particulate shear stresses by secondary effect. Then, an appropriate modelling ac counting for influences of unresolved structures has to be proposed for coarse-grid simu-lations. For this purpose, computational grids are refined to get mesh-independent result where statistical quantities do not change with further mesh refinement for a 3-D peri-odic circulating fluidized bed. The 3-D periodic circulating fluidized is a simple academic configuration where gas-solid flow conducted with A-type particles is periodically driven along the opposite direction of the gravity. The particulate momentum and agitation equations are filtered by the volume averaging and the importance of additional terms due to the averaging procedure are investigated by budget analyses using the mesh independent result. Results show that the filtered momentum equation of phases can be computed on coarse grid simulations but sub-grid drift velocity due to the sub-grid correlation between the local fluid veloc- ity and the local particle volume fraction and particulate sub-grid shear stresses must be taken into account. In this study, we propose functional and structural models for sub- grid drift velocity, written in terms of the difference between the gas velocity-solid volume fraction correlation and the multiplication of the filtered gas velocity with the filtered solid volume fraction. Particulate sub-grid shear stresses are closed by models proposed for single turbulent flows. Models’ predictabilities are investigated by a priori tests and they are validated by coarse-grid simulations of 3-D periodic circulating, dense fluidized beds and experimental data of industrial scale circulating fluidized bed in manner of a posteriori tests

Nous avons étudié des méthodes de volumes finis pour la simulation numérique d'un flux impliquant deux phases incompressibles ou deux phases générales compressibles en déséquilibre mécanique. Les principales difficultés du régime où il y a une apparition de phase ou une disparition de phase est la singularité de la vitesse. Nous montrons que l'utilisation du l'entropie correction améliorer beaucoup ces problèmes. Enfin, nous simulons certains tests numériques importants pour vérifier les méthodes numériques, telles que la séparation de phase par gravité ou un canal bouillant<br>We investigated some finite volume methods for the numerical simulation of a flow involving two incompressible phases or general two compressible phases in mechanical disequilibrium. The main difficulties of the regime where there is either a phase appearance or a phase disappearance is the singularity of the velocity. We show that using the entropy fix will much improve these problems. Finally, we perfom some important numerical tests to verify the numerical methods, such as a phase separation by gravity or a boiling channel