京都大学 / 新制・課程博士 / 博士(理学) / 甲第24421号 / 理博第4920号 / 新制||理||1703(附属図書館) / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)教授 一本 潔, 准教授 浅井 歩, 教授 横山 央明 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

京都大学 / 新制・課程博士 / 博士(理学) / 甲第24415号 / 理博第4914号 / 新制||理||1702(附属図書館) / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)准教授 浅井 歩, 教授 一本 潔, 教授 横山 央明 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

Magnetic flux emergence is the subject of how magnetic fields from the solar interior can rise and expand into the atmosphere to produce active regions. It is the link that joins dynamics in the convection zone with dynamics in the atmosphere. In this thesis, we study many aspects of magnetic flux emergence through mathematical modelling and computer simulations. Our primary aim is to understand the key physical processes that lie behind emergence. The first chapter introduces flux emergence and the theoretical framework, magnetohydrodynamics (MHD), that describes it. In the second chapter, we discuss the numerical techniques used to solve the highly non-linear problems that arise from flux emergence. The third chapter summarizes the current literature. In the fourth chapter, we consider how changing the geometry and parameter values of the initial magnetic field can affect the dynamic evolution of the emerging magnetic field. For an initial toroidal magnetic field, it is found that its axis can emerge to the corona if the tube’s initial field strength is large enough. The fifth chapter describes how flux emergence models can produce large-scale solar eruptions. A 2.5D model of the breakout model, using only dynamic flux emergence, fails to produce any large scale eruptions. A 3D model of toroidal emergence with an overlying magnetic field does, however, produce multiple large-scale eruptions and the form of these is related to the breakout model. The sixth chapter is concerned with signatures of flux emergence and how to identify emerging twisted magnetic structures correctly. Here, a flux emergence model produces signatures found in observations. The signatures from the model, however, have different underlying physical mechanisms to the original interpretations of the observations. The thesis concludes with some final thoughts on current trends in theoretical magnetic flux emergence and possible future directions.

This thesis considers the properties of MHS equilibria formed through non-resistive MHD relaxation of analytical non-potential magnetic field models, which contain two null points connected by a generic separator. Four types of analytical magnetic fields are formulated, with different forms of current. The magnetic field model which has a uniform current directed along the separator, is used through the rest of this thesis to form MHS equilibria and to study reconnection. This magnetic field, which is not force-free, embedded in a high-beta plasma, relaxes non-resistively using a 3D MHD code. The relaxation causes the field about the separator to collapse leading to a twisted current layer forming along the separator. The MHS equilibrium current layer slowly becomes stronger, longer, wider and thinner with time. Its properties, and the properties of the plasma, are found to depend on the initial parameters of the magnetic field, which control the geometry of the magnetic configuration. Such a MHS equilibria is used in a high plasma-beta reconnection experiment. An anomalous resistivity ensures that only the central strong current in the separator current layer is dissipated. The reconnection occurs in two phases characterised by fast and slow reconnection, respectively. Waves, launched from the diffusion site, communicate the loss of force balance at the current layer and set up flows in the system. The energy transport in this system is dominated by Ohmic dissipation. Several methods are presented which allow a low plasma-beta value to be approached in the single-separator model. One method is chosen and this model is relaxed non-resistively to form a MHS equilibrium. A twisted current layer grows along the separator, containing stronger current than in the high plasma-beta experiments, and has a local enhancement in pressure inside it. The growth rate of this current layer is similar to that found in the high plasma-beta experiments, however, the current layer becomes thinner and narrower over time.

Ultra Low Frequency (ULF) waves are a ubiquitous feature of Earth's outer atmosphere, known as the magnetosphere, having been observed on the ground for almost two centuries, and in space over the last 50 years. These waves represent small oscillations in Earth's magnetic field, most often as a response to the external influence of the solar wind. They are important for the transfer of energy throughout the magnetosphere and for coupling different regions together. In this thesis, various features of these oscillations are considered. A detailed background on the history and previous study of ULF waves relevant to our work is given in the introductory chapter. In the following chapters, we predominantly use numerical methods to model ULF waves, which are carefully developed and thoroughly tested. We consider the application of these methods to reports on ground and spaced based observations, which allows a more in depth study of the data. In one case, the simulation results provide evidence for an alternative explanation of the data to the original report, which displays the power of theoretical modelling. An analytical model is also constructed, which is tested on simulation data, to identify the incidence and reflection of a class of ULF wave in the flank magnetosphere. This technique is developed with the aim of future applications to satellite data. Further to this, we develop models both in Cartesian and dipole geometries to investigate some of the theoretical aspects of the coupling between various waves modes. New light is shed on the coupling of compressional (fast) and transverse (Alfvén) magnetohydrodynamic (MHD) wave modes in a 3D dipole geometry. Overall, this thesis aims to develop useful numerical models, which can be used to aid in the interpretation of ULF wave observations, as well as probing new aspects of the existing wave theory.

In this dissertation, we are concerned with the numerical simulation for flows of electrically conducting fluids exposed to an external magnetic field (also known as magnetohydrodynamics or in short MHD). The aim of the present dissertation is twofold. First, the in-house CFD hydrodynamic solver SFELES is extended to MHD problems. Second, MHD turbulence is studied in the simple configuration of a MHD pipe flow within an external transverse magnetic field. Chapter 2 of this dissertation aims at reminding the physical equations that govern incompressible MHD problems. Two equivalent formulations are put forward in the particular case of quasi-static MHD. Chapter 3 is devoted to the detailed development of the hybrid spectral - stabilized finite element methods for quasi-static MHD problems. The extension of SFELES is made for both Cartesian and axisymmetric systems of coordinates. The short chapter 4 follows to provide the performances of SFELES executed by several processes in a parallel environment. The addition of a parallel direct solver is studied in regards with the memory and time requirements. The extension of SFELES is then validated in chapter 5 with test cases of increasing complexity. For this purpose, laminar flows with an existing analytical-asymptotic solution are considered. The subject of chapter 6 is the MHD turbulent pipe flow within an external transverse and uniform magnetic field. The results are partially compared with the corresponding hydrodynamic flow and with a few data available in the literature. / Le thème de cette thèse de doctorat est la simulation numérique d'écoulements de fluides conducteurs d'électricité qui sont exposés à un champ magnétique extérieur (également connu sous le nom de magnétohydrodynamique ou encore MHD). L'objectif de ce travail est double. Premièrement, le code CFD maison SFELES est étendu aux problèmes MHD. Deuxièmement, la turbulence MHD est étudiée dans la configuration de l'écoulement en conduite cylindrique à l'intérieur d'un champ magnétique transverse. Le chapitre 2 de cette thèse a pour but de rappeler les équations qui gouvernent les problèmes de MHD incompressible. Deux formulations équivalente sont mises en évidence dans le cas particulier de la MHD quasi-statique. Le chapitre 3 est dévoué au développement détaillé des méthodes spectrale - éléments finis pour la MHD quasi-statique. L'extension de SFELES est réalisée dans les systèmes de coordonnées cartésiennes et axisymétriques. Le court chapitre 4 suit pour fournir les performances de SFELES exécuté sur plusieurs processeurs dans un environnement parallèle. L'ajout d'un solveur parallèle direct est étudié en ce qui concerne les demandes en temps et mémoire. L'extension de SFELES est alors validée dans le chapitre 5 avec des cas d'étude de complexité croissante. Dans ce but, des écoulements laminaires avec solution théorique-asymptotique sont envisagés. Le sujet du chapitre 6 est l'écoulement MHD turbulent en conduite cylindrique à l'intérieur d'un champ magnétique transverse et uniforme. Les résultats sont partiellement comparés avec l'écoulement hydrodynamique correspondant et avec des données disponibles dans la littérature. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Dans cette dissertation, nous considérons l’écoulement des liquides conducteurs d’électricité dans un champ magnétique externe. De tels écoulements sont décrits par les équations de la magnétohydrodynamique (MHD) quasi-statique, et sont fréquemment rencontrés dans des applications pratiques. Il suit qu’il y a un intérêt fort pour des outils numérques qui peuvent simuler ces écoulements dans des géometries complexes.<p>La première partie de cette thèse (chapitres 2 et 3) est dédiée à la présentation de la machinerie numérique qui a été utilisée et implémentée afin de résoudre les équations de la MHD quasi-statistique (incompressible). Plus précisément, nous avons contribué au développement d’un solveur volumes finis non-structuré parallèle. La discussion sur ces méthodes est accompagnée d’une analyse numérique qui est aussi valable pour des mailles non-structurées. Dans le chapitre 3, nous vérifions notre implémentation par la simulation d’un certain nombre de cas tests avec un accent sur des écoulements dans un champ magnétique intense.<p>Dans la deuxième partie de cette thèse (chapitres 4-6), nous avons utilsé ce solveur pour étudier des écoulements MHD de proche paroi .La première géometrie considérée (chapitre 4) est celle d’une conduite circulaire infini d’axe à haut nombre de Hartmann. Nous avons investitgué la sensitivité des résultats numériques au schéma de discrétisation et à la topologie de la maille. Nos résultats permettent de caractériser in extenso l’écoulement MHD dans une conduite avec des bords bien conducteurs par moyen des lois d’échelle.<p>Le sujet du cinquième chapitre est l’écoulement dans une conduite toroïdale à section carée. Une étude du régime laminaire confirme une analyse asymptotique pour ce qui concerne les couches de cisaillement. Nous avons aussi effectué des simulations des écoulements turbulents afin d’évaluer l’effet d’un champ magnétique externe sur l’état des couches limites limites.<p>Finalement, dans le chapitre 6, nous investiguons l’écoulement MHD et dans un U-bend et dans un coude arrière. Nous expliquons comment générer une maille qui permet de toutes les couches de cisaillement à un coût computationelle acceptable. Nous comparons nos résultats aux solutions asymptotiques. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

MSc (Mathematics) / Department of Mathematcs and Applied Mathematics / Magnetohydrodynamics ows have gained signi cant attention due to their importance in engineering applications. In this study, we numerically analysed the Dufour and Soret e ects on an unsteady MHD mixed convection ow past an in nite vertical plate with viscous dissipation. The governing non-linear partial di erential equations (PDEs) are transformed into a system of ordinary di erential equations (ODEs) by the suitable similarity transformations. The resulting equations consist of the momentum, energy and mass di usion equations. These resulting equations are solved using the Spectral Local Linearization Method (SLLM). Results obtained by the SLLM are in good agreement with the bvp4c technique. The e ects of di erent physical parameters entering into the problem are displayed graphically. The values of the Skin-friction (f0(0)), Nusselt number (􀀀 0(0)) and Sherwood number (􀀀 0(0)) are shown in tabular form for di erent values of the parameters. From the results, it is noted that the Soret number (Sr) and the Dufour number (Du) have negligible e ects on temperature pro le, whereas the decrease in the Soret number (Sr) leads to a decrease in both velocity and concentration of the uid, and the increase in Dufour number (Du) reduces the velocity and also has negligilbe e ect on the concentration pro le. / NRF

MSc (Mathematics) / Department of Mathematics and Applied Mathematics / The unsteady hydromagnetic chemically reacting mixed convection MHD ow over a permeable stretching sheet embedded in a porous medium with thermal radiation and heat source/sink is investigated numerically. The original partial di erential equations are converted into ordinary di erential equations by using similarity transformation. The governing non-linear partial di erential equations of Momentum, Energy, and Concentration are considered in this study. The e ects of various physical parameters on the velocity, temperature, and species concentration have been discussed. The parameters include the Prandtl number (Pr), Magnetic parameter (M), the Schmidt number (Sc), Unsteady parameter (A), buoyancy forces ratio parameter (N), Chemical reaction (K), Radiation parameter (Nr), Eckert number (Ec), local heat source/sink parameter (Q) and buoyancy parameter due to temperature ( ). The coe cient of Skin friction and Heat transfer are investigated. The coupled non-linear partial di erential equations governing the ow eld have been solved numerically using the Spectral Relaxation Method (SRM). The results that are obtained in this study are then presented in tabular forms and on graphs and the observations are discussed. / NRF

No description available.