[pt] O escoamento bifásico no regime de golfadas é um padrão de escoamento que requer esforço em sua caracterização e modelagem, devido às características marcantes da distribuição espacial das fases, que gera intermitência ao escoamento. Este escoamento ocorre na presença de uma grande faixa de vazões de gás e líquido, em tubulações de diâmetro médio e pequeno, com variação periódica da densidade, fração de vazio e pressões na seção transversal da tubulação. No presente trabalho estuda-se numericamente o regime de golfadas ao longo de tubulações horizontais utilizando-se o Modelo de Dois Fluidos. Uma análise detalhada das estatísticas do escoamento é realizada com a caracterização dos principais parâmetros de uma golfada, como comprimento, velocidade da frente e cauda e freqüência de passagem das golfadas ao longo da tubulação. Os dados obtidos das simulações são analisados através de histogramas de distribuição assim como de funções densidade de probabilidade (PDF) das variáveis hidrodinâmicas, destacando-se o caráter estocástico do escoamento no padrão de golfadas. Adicionalmente, comparações com dados de trabalhos experimentais da literatura foram realizadas com concordância muito boa. / [en] Slug flow is a two-phase flow pattern that requires large effort in its characterization and modeling, due to special characteristics of the phase’s spatial distribution, which causes flow intermittency. This flow occurs in a wide range of gas and liquid flow rates in pipes of medium and small diameters, with periodic variation of density, void fraction and pressures in pipe-cross-section. This work presents a numerical study of the slug regime through horizontal pipes using the Two-Fluid Model. A detailed statistical analysis of the flow was carried out with characterization of main slug parameters, such as slug length, front and tail velocities and slug frequency along the pipeline. The numerical results were analyzed through distribution histograms as well as probability density function (PDF) of the hydrodynamic variables, showing the stochastic characteristic of slug flow pattern. Further, comparisons with experimental data from the literature were performed, showing very good agreement.

Due to the high-speed flow in a chute spillway, cavitation damages often occur. This undesired phenomenon threatens the safety of the structure. For the purpose of eliminating the damages, an aerator is often installed in the spillway. To understand its characteristics, physical model tests are a popular method. To complement the model tests, computation fluid dynamics (CFD) simulations are used to study aerator flows. To represent the two-phase flows, multiphase models should be employed. This thesis examines two of them, namely, the Volume-Of-Fluid model (VOF) and Two-Fluid model. Based on the background of the Bergeforsen dam, the aerator flow is modelled by means of the VOF model. The simulated spillway discharge capacity is in accordance with the experimental data. Compared with the results, empirical formulas fail to evaluate the air supply capacity of aerator as it is wider than the conventional width. A hypothetical vent modification is proposed. For the original and proposed layouts, the study illustrates the difference in the air-flow conditions. The results show that a larger vent area is, for a large-width aerator, preferable in the middle of the chute. To study the flip bucket-shaped aerators in the Gallejaur dam, physical model tests and prototype observations are conducted. The results lead to contradicting conclusions in terms of jet breakup and air entrainment. A CFD model is, as an option, employed to explain the reason of the discrepancy. The numerical results coincide with the prototype observations. The jet breakup and air entrainment are evaluated from air cavity profiles; the air-pressure drops are small in the cavity. The discrepancy is due to overestimation of the surface-tension effect in the physical model tests. Based on the experimental data of an aerator rig at the Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, the Two-Fluid model is used to predict air concentration distributions in the aerated flow. The model includes relevant forces governing the motion of bubbles and considers the effects of air bubble size. The numerical results are conformable to the experiments in the air cavity zone. Downstream of the cavity, the air concentration near the chute bottom is higher, which is presumably caused by the fact that the interfacial forces in the Two-Fluid model are underestimated. / <p>QC 20170224</p>

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

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.

Filtered two-fluid models used to perform large scale simulations of gas-solid fluidized flows of industrial risers require closures for filtered parameters such as filtered and residual stresses, and interphase interaction forces mainly effective drag. Closure models for those filtered parameters may be derived by averaging over results of highly resolved simulations with microscopic two-fluid modeling. This work is a contribution in that context. Recent models for filtered parameters have been written as functions of filter size, filtered solid volume fraction, and filtered slip velocity. A recent study showed that macro-scale variables like domain average solid volume fraction and gas Reynolds number also significantly affect the filtered parameters. In the current work, in addition to these filtered and macro-scale variables, the effects of two new variables over the filtered parameters are investigated: filtered solid kinetic energy and sub-grid gas turbulence. It is shown that the filtered solid kinetic energy should be accounted for in the concerning correlations, thereby improving accuracy. Regarding gas turbulence, literature shows it has no significant effects on the motion of high Stokes particles. Extending on literature, this work investigates the sub-grid gas turbulence effects on meso-scale structures formed of high Stokes particles. Results showed that sub-grid gas turbulence has no significant effects on the meso-scale structures and corresponding filtered parameters. The open source code MFIX was used for all simulations. Ranges of dilute concentration of solid and gas Reynolds number typical of riser flow regimes were considered. A modified two-fluid model with microscopic formulation was used. The sub-grid gas turbulence was generated by means of a forcing function procedure which was implemented in the physical space, over the gravitational term of the momentum conservation equation of the gas phase. First, numerical simulations of the gas phase alone were performed, accounting for literature available data, in order to set a turbulent gas field and calibrate the turbulence intensity. Then, the forcing function was introduced in to the two-fluid model and various gas-solid flows were simulated. While the current results show the necessity of accounting for additional variables in the filtered parameter correlation, they also make it clear the necessity of further developments that are required in the search for better accuracy. / Modelos filtrados de dois-fluidos usados em simulações de grandes escalas de escoamentos fluidizados de gás-sólido de risers industriais exigem fechamentos para parâmetros filtrados tais como as tensões filtradas e residuais, e forças interativas interfases, principalmente arrasto efetivo. Modelos de fechamento para estes parâmetros filtrados podem ser gerados a partir de procedimentos de media aplicados sobre resultados de simulações altamente resolvidas com modelagem microscópica de dois-fluidos. Este trabalho é uma contribuição neste contexto. Modelos de fechamento recentes para parâmetros filtrados tem sido formulados em função de tamanho de filtro, fração volumétrica de sólido filtrada, e velocidade de deslizamento filtrada. Estudo recente mostrou que variáveis de macro-escalas como fração volumétrica de sólido e número de Reynolds de gás médios no domínio também afetam significativamente os parâmetros filtrados. No presente trabalho, além dessas variáveis filtradas e de macro-escala, os efeitos de duas novas variáveis sobre os parâmetros filtrados são investigados: energia cinética filtrada do sólido e turbulência submalha do gás. Em relação à energia cinética filtrada do sólido, mostra-se que a sua consideração refina as correlações em questão, contribuindo assim para melhor acuracidade. Com relação à turbulência do gás, a literatura mostra que não tem efeitos significativos no movimento de particulados de elevados números de Stokes. Acrescentando à literatura, este trabalho investiga os efeitos da turbulência sub-malha do gás sobre estruturas de meso-escala formados de particulados de elevados números de Stokes. Os resultados mostraram que a turbulência sub-malha do gás não tem efeitos significativos sobre estruturas de meso-escalas e parâmetros filtrados correspondentes. O código aberto MFIX foi usado para todas as simulações. Faixas de concentração diluída de sólido e número de Reynolds típicos de escoamentos em risers foram considerados. Um modelo modificado de dois fluidos com formulação microscópica foi utilizado. A turbulência sub-malha do gás foi gerada por meio de um procedimento de \'forcing function\' que foi implementado no espaço físico, sobre o termo fonte gravitacional da equação de momentum da fase gás. Primeiramente, simulações numéricas da fase gás foram realizadas separadamente, levando-se em conta dados disponíveis na literatura, a fim de gerar um campo de gás turbulento e calibrar a intensidade de turbulência. Posteriormente, a \'forcing function\' foi introduzida no modelo de dois-fluidos e vários escoamentos de gás-sólido foram simulados. Enquanto os resultados obtidos mostram a necessidade de consideração de variáveis adicionais para correlação de parâmetros filtrados, também deixam claro a necessidade de desenvolvimentos mais aprofundados na busca de melhor acuracidade.

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.

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.

Dans cette thèse, on propose un nouveau modèle de couplage auto-cohérent entre des structures magnétiques cohérentes sur les échelles ioniques et des ondes dites de sifflement (whistlers, en anglais) à plus hautes fréquences, afin d'interpréter les données expérimentales recueillies par les satellites Cluster pendant un sous-orage magnétique dans la région nocturne de la magnétosphère terrestre. Le couplage fournit un mécanisme pour confiner et transporter les ondes whistlers par l'intermédiaire d' une onde nonlinéaire qui se propage obliquement par rapport au champ magnétique. Cette étude s'appuie sur une analyse des données expérimentales, sur une modélisation théorique ainsi que sur des simulations numériques. Pendant les sous-orages magnétiques, la magnétosphère est soumise à de fortes perturbations magnétiques et électriques dans une vaste gamme de fréquences, qui vont des basses fréquences, inférieures ou de l'ordre de l'échelle temporelle typique ionique, aux hautes fréquences, supérieures ou de l'ordre de l' échelle temporelle typique électronique. Afin de connaître les processus physiques qui déterminent la dynamique de la magnétosphère pendant les sous-orages, il est fondamental de comprendre si, et avec quel méchanisme, des couplages peuvent se produire entre des ondes qui se propagent sur des temps caractéristiques différents. Des structures magnétiques à basse fréquence ont déjá été obsérvées dans des régions comme la magnétogaine et le vent solaire, éventuellement associées à des ondes whistlers à plus haute fréquence. Dans cette thèse, on montre que des structures similaires sont obsérvées dans la couche de plasma à l'intérieur de la magnétosphère. On s'intérroge ensuite sur la façon dont l'inhomogénéité de telles structures peut influencer la propagation des ondes à plus haute fréquence. Grâce à ses quatre satellites en configuration tetraédrique et à ses mésures à haute résolution temporelle, la mission Cluster nous offre une occasion unique de pouvoir analyser la structure spatiale des perturbations stationnaires (ou se propageant) et d'étudier la dynamique du plasma sur des échelles temporelles plus courtes, telles que celles des ondes whistlers. Ainsi, je décrirai les émissions d'ondes whistlers détectées par les satellites Cluster à l'intérieur de structures magnétiques cohérentes situées dans un écoulement de plasma rapide pendant le sous-orage du 17 Août 2003. Au cours de cette période, les satellites Cluster sont situés dans la couche de plasma, séparés d'une distance de l'ordre des échelles spatiales typiques ioniques (le rayon de giration ou la longueur d'inertie des ions). Les ondes whistlers sont corrélées avec des structures magnétiques characterisées par un minimum du module du champ magnétique et un maximum de densité du plasma. Ces dernières ont été modélisées comme des ondes planes nonlinéaires de type lent qui piègent et transportent les ondes whistlers. A partir d'une étude théorique et numérique en utilisant une approche bi-fluide, on peut alors reproduire les données observationnelles. Le rôle possible de telles structures couplées dans la physique des sous-orages est aussi discuté. Ce nouveau mécanisme de piégeage, étudié ici en utilisant comme guide pour les whistlers une onde oblique de type magnétosonique, est d'intérêt plus général par rapport au contexte spécifique des observations présentées dans cette thèse. En effet, d'autres ondes nonlinéaires, comme par exemple les ondes d' Alfvén obliques ou d' Alfvén cinétiques dans les plasmas à beta fort (où beta est le rapport de la pression thermique du plasma sur la pression magnétique), pourraient aussi transporter les whistlers. Ce modèle de piégeage constitue aussi une explication alternative aux modèles existants qui considèrent une inhomogénéité stationnaire sous la forme d'un canal de densité. Enfin, l'étude décrite dans cette thèse concerne des problématiques fondamentales en physique des plasmas, comme la propagation d'ondes dans les milieux inhomogènes et l'interaction entre modes sur des échelles temporelles différentes.

An arterial traffic control system influences the travel time characteristics of a corridor, including the average corridor travel time and the travel time reliability. However, reliability measures have typically been outside of the focus of arterial control system performance evaluation studies. To assess the effectiveness of arterial traffic control performance evaluation studies are normally limited to average measures of travel time, speed, or delay. As an advanced traffic management system, adaptive traffic control has been developed to address real time demand variability. Thus, an evaluation of the adaptive traffic control system based on reliability may be as important as evaluation based on average travel time or delay. In addition, arterial control systems may also affect the performance of side street traffic as well as arterial corridor traffic. The performance of side street traffic is another measure that should be used in the assessment of the effectiveness of any arterial traffic control system. Finally, an arterial's operational performance often changes throughout a day and over the arterial length. Thus, a system-wide measure that reflects the range of observed operations is needed to thoroughly assess the performance. Given these issues the goal of this research is the development of procedures to evaluate adaptive traffic control's effect on arterial characteristics such as travel time distribution, reliability, side street performance, and system-wide performance. The developed procedures were applied to the evaluation of an adaptive traffic control system, SCATS (Sydney Coordinated Adaptive Traffic System) in Cobb County, Georgia that replaced a semi-actuated coordinated control system. After the procedures were applied, it was found that SCATS produced a less extreme shape of travel time distribution, possibly due to the adaptive feature, but that it did not make statistically significant changes in the selected overall analysis measures. Also, it was found that the results of the performance evaluation can vary depending on the measures selected or the study period and location.