[pt] Uma análise numérica de escoamento anular vertical ascendente é realizada utilizando o Modelo de Dois Fluidos unidimensional. Para escoamentos verticais, na ausência de mecanismos estabilizadores, o sistema de equações resultante é incondicionalmente mal posto. Dessa forma, visando tornar o sistema de equações bem posto, adicionou-se às equações de quantidade de movimento um termo de pressão dinâmica. Dois modelos disponíveis na literatura são investigados. O primeiro só considera a pressão dinâmica na equação de quantidade de movimento do líquido, a qual é baseada em uma simples expressão para a velocidade da onda na interface. O segundo modelo inclui uma contribuição da pressão dinâmica em ambas equações de quantidade de movimento. No presente trabalho, um terceiro modelo é proposto, o qual é baseado no primeiro modelo, com uma avaliação mais realista da velocidade da onda na interface. O efeito estabilizante da pressão dinâmica é demonstrado através de um rigoroso teste de convergência de malha. As equações de conservação são discretizadas com o método de volumes finitos, com uma integração temporal de primeira ordem e uma discretização espacial TVD de segunda ordem. Tanto o segundo quanto o terceiro modelo considerados apresentaram solução independente da malha. Parâmetros do escoamento como gradiente de pressão, espessura do filme, e variáveis características da onda obtidos numericamente são comparados com dados experimentais disponíveis na literatura, apresentando boa concordância. / [en] A numerical analysis of vertical ascending annular flow with the 1D Two-Fluid model is performed. It is well known that, in vertical flows, the resulting system of equations is unconditionally ill-posed in the absence of stabilizing mechanisms. Therefore, in the present work, to render the system of equations well-posed, modeling of dynamic pressure is included in the momentum equations. Two models available in the literature are examined. The first one only considers the dynamic pressure in the liquid momentum equation, which is based on a simple expression for the interface wave velocity. The second model includes a dynamic pressure contribution to both momentum equations. A third model is proposed based on the first, with a more realistic estimation of the interface wave velocity. A systematic grid convergence test is performed to demonstrate the stabilizing effect of the dynamic pressure. The conservation equations are discretized with the finite volume method, with a first order time integration, and a second order TVD spatial discretization. A grid-independent solution can be found when applying the second and third models considered. Flow parameters such as pressure drop, film thickness and wave characteristics numerically obtained are compared against available experimental data, presenting good agreement.

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.

Les instabilités fluides peuvent dégrader le confinement du plasma au sein des tokamaks. Étant données les échelles spatio-temporelles, on choisit les modèles fluides obtenus à partir de la dérivation des modèles cinétiques. On dérive plusieurs modèles hiérarchiques de la MagnétoHydroDynamique (MHD) et en particulier les modèles de la MHD réduite du Current Hole et de l'équilibre de Grad-Shafranov. Une des difficulté de l'ensemble de ces modèles est de respecter l'équation modélisant l'absence de monopôles magnétiques. Pour assurer cette condition en tout point du domaine, le champ magnétique est réécrit avec un potentiel vecteur. L'utilisation de potentiels fait apparaître des équations faisant intervenir des dérivées d'ordre supérieurs. La stratégie numérique développée est l'utilisation de la méthode des éléments finis avec des éléments C1 de Bell. Sur un maillage non structuré, ces éléments ont l'intérêt de présenter une base réduite définir exclusivement avec des variables aux noeuds du maillage. Les modèles de MHD réduite du Current Hole et de Grad-Shafranov ont été résolus avec ces éléments. La résolution du cas test de Grad-Shafranov avec les conditions de bords exactes a permis d'obtenir l'ordre optimale de 5. La résolution du système du Current Hole avec ces éléments, validée par l'obtention du paramètre η1/3, a permis l'observation de développement d'instabilités en dents de scies. / Fluid instabilities can degrade plasma confinement in tokamaks. Given the spatial and temporal scales, we choose the fluid models obtained from the derivation of kinetic models. We derived several hierarchical models of MagnetoHydroDynamic (MHD) and in particular models of reduced MHD like the Current Hole and the Grad-Shafranov equilibrium. One of the difficulty of all these models is to respect the absence of magnetic monopoles equation. To ensure this condition at any point, the magnetic field is rewritten with a vector potential. The use of vector portential implies that higher order derivatives appear in the equation. The numerical strategy is developed using the finite element method with C1 Bell's elements. On a unstructured mesh, these have the advantage to present a reduced basis with degrees of freedom defined exclusively on the nodes of the mesh. The reduced MHD models of the Current Hole and Grad-Shafranov have thus been resolved with these elements. The resolution of a Grad-Shafranov test case with exact boundary conditions yields the optimal order of 5. The resolution of the Current Hole system with thesse elements has been validated by obtaining physical parameter η1/3 and allowed the observation of the development of sawtooth instabilities.

Cette thèse contribue à la vérification et à la validation du modèle bifluide de Baer-Nunziato, pour modéliser les phénomènes de transitoires hydrauliques dans les réseaux de tuyauteries industrielles. Il s’agit d’abord de modéliser les écoulements de transitoires hydrauliques avec le modèle bifluide en représentation eulérienne, puis d’étendre ce modèle en formalisme ALE (Arbitrary Lagrangian Eulerian) pour prendre en compte l’interaction fluide-structure (IFS). Pour modéliser les écoulements, des lois de fermetures du modèle bifluide concernant les termes interfaciaux, les termes sources et les lois thermodynamiques (EOS) ont d’abord été étudiées. Ensuite, le système complet a été simulé avec une méthode à pas fractionnaires qui admet deux étapes, l’une pour la résolution de la partie convective, l’autre pour les termes sources. L’ensemble de schémas a été vérifié et étendu aux EOS ‘Stiffened Gas généralisées’ afin de représenter le changement de phase eau-vapeur. Après avoir retrouvé certains phénomènes typiques associés aux transitoires hydrauliques, le modèle bifluide a été validé avec l’expérience de Simpson, l’expérience Canon, et comparé avec deux modèles homogènes sur ces deux expériences. Enfin, une version ALE du modèle bifluide a été mise en œuvre et vérifiée sur un cas de propagation d’ondes de pression dans une conduite flexible. La variation de la célérité des ondes dans le fluide liée au couplage fluide/structure a été bien retrouvée. La validation a été effectuée sur un cas expérimental d'explosion dans une tuyauterie en eau. Les simulations sont en bon accord avec les données expérimentales. / This thesis contributes to the verification and the validation of the Baer-Nunziato (BN) model, to modelize water hammer phenomena in industrial piping systems. It consists of two parts, the first is to modelize water hammer flows with the BN model in Eulerian representation and the second is to extend this model to the ALE (Arbitrary Lagrangian Eulerian) formalism so as to take into account fluid-structure interaction (FSI). To modelize water hammer flows, closure laws of the BN model concerning the interfacial/source terms and the equations of state (EOS) were first studied. Then the whole system was simulated with a fractional step method including two steps, one for the resolution of the convective part, the other for the source terms. All schemes have been extended to ‘generalized Stiffened Gas’ EOS in order to represent phase-change. After regaining some typical phenomena associated with water hammer flows, the BN model was validated with the Simpson experiment, a classical water hammer test case, and the Canon experience, a rapid decompression of fluid in a high pressure duct. Moreover, the model was compared with two homogeneous models on both experiments. Finally, an ALE version of the BN model was implemented, and verified on a case of wave propagation in a ‘single’ phase flow and a two-phase flow in a flexible pipe. The variation of wave propagation speed due to the coupling between the fluid and the structure has been well retrieved. The validation was performed on an experiment which examines the response of a pipe filled with water, subjected to a violent pressure peak (140 bar). The simulations are in good agreement with the experimental data.

Reatores de leito fluidizado circulante são intensamente utilizados em aplicações de larga escala como craqueamento catalítico de petróleo e combustão de carvão. Projeto e desenvolvimento nestas áreas são fortemente baseados em plantas de demonstração, a custos elevadíssimos. Nesse contexto, tratamentos utilizando mecânica dos fluidos computacional assumem considerável relevância. Os complexos padrões de escoamento gás-sólido que se desenvolvem nos reatores de leito fluidizado circulante determinam taxas de reação e exigem, portanto, descrições hidrodinâmicas rigorosas. Os modelos Eulerianos do contínuo ou de dois fluidos são correntemente considerados a escolha mais prática na busca destas descrições. Estas formulações são baseadas na aplicação da hipótese de meio contínuo tanto para fases macroscopicamente contínuas quanto dispersas. Neste trabalho, desenvolve-se simulações numéricas para descrever processos hidrodinâmicos e reativos em reatores de leito fluidizado circulante aplicando modelagem de dois fluidos. Considera-se situações típicas de combustão em leito fluidizado circulante de carvão mineral. Como processo reativo considera-se a absorção de dióxido de enxofre por calcário. Realiza-se análises hidrodinâmicas, e de efeitos hidrodinâmicos sobre a reação de interesse. Realiza-se simulações de regime permanente e de regime transiente. Mostra-se que as simulações de regime permanente permitem análises qualitativas do processo, e provêem condições iniciais para simulações transientes diretamente no regime de escoamento estatisticamente permanente. Em relação à hidrodinâmica, conclui-se que os modelos de dois fluidos correntes são ainda bastante crus. São claramente necessárias melhores descrições reológicas e relações constitutivas sub-grade mais acuradas. Em relação à reação química, conclui-se que predições reativas acuradas somente poderão ser obtidas se descrições hidrodinâmicas rigorosas forem combinadas com descrições reativas igualmente rigorosas. / Circulating fluidized bed reactors are widely used in large scale applications such as catalytic cracking of petrol and coal combustion. Development and design in those areas are strongly based on demonstration plants, at extremely high costs. In this context, treatments applying computational fluid mechanics assume considerable relevance. The complex gas-solid flow patterns which develop inside the circulating fluidized bed reactors determine reaction rates, so that rigorous hydrodynamic descriptions are required. The continuum Eulerian or two-fluid models are currently considered the more practical choice for providing such descriptions. Those formulations are based on the application of the continuum hypothesis for both macroscopically continuous and dispersed phases. In this work numerical simulation is performed to describe both hydrodynamics and reactive processes in circulating fluidized beds applying two-fluid modeling. Typical situations of circulating fluidized bed coal combustion are considered. The reactive process considered is the absorption of sulfur dioxide by limestone. Analyses are performed of hydrodynamics, and regarding hydrodynamic effects over the concerning reaction. Both steady state and transient simulations are performed. It is shown that steady state simulations allow qualitative analyses, and do provide initial conditions for transient runs straightly inside the statistical steady state flow regime. Concerning hydrodynamics, it is concluded that the current two-fluid models are still very crude. Clearly, better rheological descriptions are required alongside with more accurate sub-grid constitutive relations. Regarding chemical reaction, it is concluded that accurate reactive predictions shall only be found if rigorous hydrodynamic descriptions are combined with equally rigorous reaction descriptions.

Reatores de leito fluidizado circulante são intensamente utilizados em aplicações de larga escala como craqueamento catalítico de petróleo e combustão de carvão. Projeto e desenvolvimento nestas áreas são fortemente baseados em plantas de demonstração, a custos elevadíssimos. Nesse contexto, tratamentos utilizando mecânica dos fluidos computacional assumem considerável relevância. Os complexos padrões de escoamento gás-sólido que se desenvolvem nos reatores de leito fluidizado circulante determinam taxas de reação e exigem, portanto, descrições hidrodinâmicas rigorosas. Os modelos Eulerianos do contínuo ou de dois fluidos são correntemente considerados a escolha mais prática na busca destas descrições. Estas formulações são baseadas na aplicação da hipótese de meio contínuo tanto para fases macroscopicamente contínuas quanto dispersas. Neste trabalho, desenvolve-se simulações numéricas para descrever processos hidrodinâmicos e reativos em reatores de leito fluidizado circulante aplicando modelagem de dois fluidos. Considera-se situações típicas de combustão em leito fluidizado circulante de carvão mineral. Como processo reativo considera-se a absorção de dióxido de enxofre por calcário. Realiza-se análises hidrodinâmicas, e de efeitos hidrodinâmicos sobre a reação de interesse. Realiza-se simulações de regime permanente e de regime transiente. Mostra-se que as simulações de regime permanente permitem análises qualitativas do processo, e provêem condições iniciais para simulações transientes diretamente no regime de escoamento estatisticamente permanente. Em relação à hidrodinâmica, conclui-se que os modelos de dois fluidos correntes são ainda bastante crus. São claramente necessárias melhores descrições reológicas e relações constitutivas sub-grade mais acuradas. Em relação à reação química, conclui-se que predições reativas acuradas somente poderão ser obtidas se descrições hidrodinâmicas rigorosas forem combinadas com descrições reativas igualmente rigorosas. / Circulating fluidized bed reactors are widely used in large scale applications such as catalytic cracking of petrol and coal combustion. Development and design in those areas are strongly based on demonstration plants, at extremely high costs. In this context, treatments applying computational fluid mechanics assume considerable relevance. The complex gas-solid flow patterns which develop inside the circulating fluidized bed reactors determine reaction rates, so that rigorous hydrodynamic descriptions are required. The continuum Eulerian or two-fluid models are currently considered the more practical choice for providing such descriptions. Those formulations are based on the application of the continuum hypothesis for both macroscopically continuous and dispersed phases. In this work numerical simulation is performed to describe both hydrodynamics and reactive processes in circulating fluidized beds applying two-fluid modeling. Typical situations of circulating fluidized bed coal combustion are considered. The reactive process considered is the absorption of sulfur dioxide by limestone. Analyses are performed of hydrodynamics, and regarding hydrodynamic effects over the concerning reaction. Both steady state and transient simulations are performed. It is shown that steady state simulations allow qualitative analyses, and do provide initial conditions for transient runs straightly inside the statistical steady state flow regime. Concerning hydrodynamics, it is concluded that the current two-fluid models are still very crude. Clearly, better rheological descriptions are required alongside with more accurate sub-grid constitutive relations. Regarding chemical reaction, it is concluded that accurate reactive predictions shall only be found if rigorous hydrodynamic descriptions are combined with equally rigorous reaction descriptions.

Ce travail a été réalisé dans le cadre du projet EcOBEx, qui consiste à réduire le bruit du groupe motopropulseur rayonné à l'extérieur par l'ajout d'écrans acoustiques dans le compartiment moteur du véhicule. Les écrans acoustiques sont fabriqués par thermocompression de matériaux poreux uniformes. Les propriétés et l'épaisseur du matériau évoluent en fonction du degré de compression subit par le matériau. L'objectif de ce travail est de proposer des lois pour prédire l'évolution des propriétés des matériaux à partir du taux de compression et de leurs valeurs initiales avant compression. Dans un premier temps, on s'intéresse aux paramètres du modèle de fluide équivalent de Johnson-Champoux-Allard-Lafarge (JCAL) : porosité, résistivité au passage d'air, tortuosité, longueurs caractéristiques visqueuse et thermique, perméabilité thermique statique. Des expressions analytiques sont proposées pour prédire la variation de ces paramètres en fonction de la compression. Elles sont développées à partir d'un modèle de matériaux fibreux à fibres cylindriques où les variations d'orientation des fibres induites par la thermocompression peuvent être prises en compte. Les résultats sont en bon accord avec les mesures effectuées sur deux types de matériaux (mousse à cellules ouvertes et fibreux). Un modèle empirique généralisé est finalement proposé pour la résistivité au passage d'air. Dans un deuxième temps, on s'attache aux paramètres élastiques dont la connaissance est essentielle pour prendre en compte la vibration du squelette. La méthode expérimentale quasistatique est d'abord appliquée pour étudier l'évolution du module de Young par rapport au taux de compression pour les fibres et les mousses. Une loi de puissance est alors proposée pour prédire ces variations. Enfin, une méthode inverse pour estimer les propriétés élastiques d'un matériau poroélastique orthotrope à partir d'une mesure vibratoire d'un écran tricouche thermo comprimé est proposée. Cette méthode permet de caractériser les propriétés élastiques du matériau poreux dans une situation proche de son application réelle / This work was carried out in the framework of the project EcOBEx, whose main objective was to reduce the passby noise by mean of acoustic shields in the engine compartment of the vehicle. The acoustic shields are manufactured by thermocompression of uniform porous materials. The material’s properties and thickness evolve according to the degree of compression experienced by the material. The objective of this work is to propose some laws to predict the evolution of the materials properties from their initial non compressed values and the compression rate. Firstly, we focus on the parameters of the Johnson-Champoux-Allard-Lafarge (JCAL) equivalent fluid model : porosity, air-flow resistivity, tortuosity, viscous characteristic lengths, thermal characteristic length, static thermal permeability. Some analytical expressions are proposed to predict the variation of these parameters as a function of compression. They are derived from a physical model of cylindrical fibres where the fibre orientation variations induced by the thermocompression can be taken into account. The results are in good agreement with the measurements made two types of materials (open cell foam and fibrous). A generalized empirical model is finally proposed for the air-flow resistivity.In a second part, we focus on the elastic parameters, which are necessary to take into account the vibration of the skeleton. The quasi-static experimental method is first applied to study the evolution of the Young’s modulus along the compression rate for fibrous and open cell foams. A power law is then proposed to predict these variations. Finally, an inverse method for estimating the elastic properties of an orthotropic poro-elastic material from a vibratory measurement of a thermocompressed three layer sandwich structure is proposed. This method allows us to characterize the elastic properties of a porous material in a situation close to its actual application.

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. / 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.

The aim of this thesis is to investigate the properties of a Radio Frequency capacitive discharge at atmospheric pressure in argon. In these conditions where the pressure x distance product is around 150 Torr.cm, the discharge usually consists of several locally hot filaments. By pulsing the RF generator with an appropriate width and period, it was found possible to control the filament to glow transition in order to obtain a diffused and stable plasma.¶ The 2 mm gap between the electrodes is open to the ambient air and fed with argon via one hundred submillimetric holes regularly spread on the surface of the top electrode. This configuration allows “on-line” surface treatment of polymer films without having to turn the discharge off between successive samples. An important and lasting improvement of the polymer wettability is quickly obtained without risk of damage. The plasma diagnostic methods are emission spectroscopy and electric measurements. The Stark broadening of the Balmer β transition line of atomic hydrogen is measured to determine a plasma density of 10^15/cm3 in the filamentary mode. The glow mode density estimation was based on power balance yielding a density of 5×10^11/cm3. Emission line ratios between neutrals and Ar+ ions are used in the Saha equation to calculate the electron temperature. It results in an approximation of 1.3 eV for the glow mode and 1.7 eV for the filaments.¶ A unidimensional self-consistent fluid model is developed to gain insight into the homogeneous discharge behaviour. Poisson’s equation for the electric field is coupled to the first moments of the Boltzmann equation (continuity equation, drift-diffusion equation and energy equation). Transport and reaction coefficients are obtained from the mean energy of the electrons.¶ The model is applied to a reduced argon kinetic with the main ionization and excitation processes. Simulation results are in agreement with experimental measurements. The atmospheric pressure RF discharge is similar to a lower pressure RF discharge for which the ionization occurs mainly inside the oscillating sheaths where electrons are the most energetic

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.