京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第22257号 / 理博第4571号 / 新制||理||1656(附属図書館) / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)教授 田口 聡, 教授 秋友 和典, 准教授 藤 浩明 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第19508号 / 理博第4168号 / 新制||理||1599(附属図書館) / 32544 / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)教授 田口 聡, 教授 家森 俊彦, 教授 余田 成男 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

This thesis describes work on building numerical models of the Earth's magnetosphere using magnetohydrodynamics (MHD) and other related modeling methods. For many years, models that solve the MHD equations have been the main tool for improving our theoretical understanding of the large-scale dynamics of the Earth's magnetosphere. While the MHD models have been very successful in capturing many large-scale features, they fail to adequately represent the important drift physics in the inner magnetosphere. Consequently, the ring current, which contains most of the particle energy in the inner magnetosphere, is not realistically represented in MHD models. In this thesis, Chapter 2 and 3 will describe in detail our effort to couple the OpenGGCM (Open Geospace General Circulation Model), one of the major MHD models, to the Rice Convection Model (RCM), an inner magnetosphere ring current model, with the goal of including energy dependent drift physics into the MHD model. In Chapter 4, we will describe an initial attempt to use a direct-integration method to calculate Birkeland currents in the MHD code. Another focus of the thesis work, presented in Chapter 5, addresses a longstanding problem on how a geomagnetic substorm can occur within the closed field line region of the tail. We find a scenario of a bubble-blob pair formation in an OpenGGCM simulation just before the expansion phase of the substorm begins and the subsequent separation of the bubble and the blob decreases the normal component of the magnetic field until finally an X-line occurs. Thus the formation of the bubble-blob pair may play an important role in changing the magnetospheric configuration from a stretched field to the X-line formation that is believed to be the major signature of a substorm.

viii Title: Dynamics of the bow shock and magnetopause Author: Karel Jelínek Department: Department of Surface and Plasma Science Supervisor: Prof. RNDr. Zdeněk Němeček, DrSc. Department of Surface and Plasma Science e-mail address: zdenek.nemecek@mff.cuni.cz Abstract: The interplanetary space is a unique laboratory which allows us to dis- cover (i) a behavior of the plasma under different conditions, (ii) origin of its insta- bilities, and (iii) its interaction with obstacles such as the Earth's magnetosphere. The present thesis analyzes the outer Earth's magnetosphere. The results are based on the in situ sensing by a variety of the spacecraft (e.g., IMP-8, INTERBALL-1, MAGION-4, Geotail, Cluster-II and Themis). The solar wind curently monitored by the WIND and ACE spacecraft near the La- grange point L1 affects by its dynamic pressure the Earth's magnetic field which acts as a counter-pressure and the boundary where these pressures are balanced is the magnetopause. Due to supersonic solar wind speed, the bow shock forms in front of the magnetopause and a region in between, where plasma flows around an obstacle is named the magnetosheath. The thesis contributes to a deaper understanding of the dependence of magnetopause and bow shock shapes and positions, especially, (1) on the orientation of the inter-...

Data visualization is a field dedicated to effectively showing large amounts of collected data. A field where data visualization has shown promising results in its ways to effectively answer questions is the fundamental research of the universe. This thesis describes how to visualize the Earth’s magnetosphere as it interacts with the solar wind, using numerical methods, semi-transparent surfaces, and contours in OpenSpace. A magnetosphere module was implemented into OpenSpace, and the OpenSpace GUI was extended to give the user a possibility to interact with the visualization. The implemented algorithm in the magnetosphere module was measured in terms of speed, robustness, and user understanding. The implementation made it possible to visualize a simple model of the Earth’s magnetosphere, both when it interacts and when it does not interact with the solar wind. The measured speed showed a trend of a linear increase when more magnetic field lines were added to the visualization, where the run time was low for all tests. The algorithm was shown to be robust in its creation of the magnetosphere. When asked about what the users thought of the implemented visualization, a majority were positive and thought it to be a good complement to learning about the Earth’s magnetosphere. For a realistic model, there is a possibility that the speed and robustness may get worse, but the results are good for a simple model. To make the visualization itself more informative, more highlights are needed for important parts of the structure, as well as different colors that differ on which surfaces belong to which geographical pole. / Creative Exploration of the Atmosphere

Ideal magnetohydrodynamics (MHD) simulations are known to have problems in satisfying the solenoidal constraint (i.e. the divergence of magnetic field should be equal to zero, $ ablacdotvec{B} = 0$). The simulations become unstable unless specific measures have been taken. In this thesis, a solenoidal constraint satisfying technique that allows discrete satisfaction of the solenoidal constraint up to the machine accuracy is presented and validated with a variety of test cases. Due to its inspiration from Chorin's artificial compressibility method developed for incompressible CFD applications, the technique was named as extit{artificial compressibility analogy (ACA)} approach. It is demonstrated that ACA is a purely hyperbolic, stable and consistent technique, which is moreover easy to implement. Unlike some other techniques, it does not pose any problems of the sort that $ ablacdotvec{B}$ errors accumulate in the vicinity of the stagnant regions of flow. With these crucial properties, ACA is thought to be a remedy to the drawbacks of the most commonly used solenoidal constraint satisfying techniques in the literature namely: Incorrect shock capturing and poor performance of the convective stabilization mechanism in regions of stagnant flow for Powell's source term method; exceedingly complex implementation for constrained transport technique due to the staggered grid representation; computationally expensive nature due to the necessity of a Poisson solver combined with hyperbolic/elliptic numerical methods for classical projection schemes. In the first chapter of the thesis, general background knowledge is given about plasmas, MHD and its history, a certain class of upwind finite volume methods, namely Riemann solvers, and their applications in MHD, the definition, constituents, formation mechanisms and effects of space weather and some of the space missions that are and will be performed in its prediction. Secondly, detailed analysis of the compressible ideal MHD equations is given in the form of the derivation of the equations, their dimensionless numbers which will be of use to specify the flows in the following chapters, and finally, the presentation of the MHD waves and discontinuities, which indicates the complexity of the system of ideal MHD equations and therefore their further numerical analysis. The next discussion is about the main subject of the thesis, namely the solenoidal constraint satisfying techniques. First of all, the definition and significance of the solenoidal constraint is given. Afterwards, the most common solenoidal constraint satisfying techniques in the literature are reviewed along with their abovementioned drawbacks. Moreover, particular emphasis is given to the Powell's source term approach which was also implemented in the upwind finite volume MHD solver developed. In addition, the hyperbolic divergence cleaning technique is presented in detail together with the resemblance and differences between it and ACA. Some other solenoidal constraint satisfying techniques are briefly mentioned at this stage. After these, ACA is presented in the following way: The point of inspiration, which is the analogy made with Chorin's artificial compressibility method developed for incompressible CFD, the introduction of the modified system of ideal MHD equations due to ACA, the derivation of the wave equation governing the propagation of $ ablacdotvec{B}$ errors and the analytical consistency proof. Having finished the core discussion of the thesis, the solver developed and its constituents are given in the fourth chapter. Furthermore, a brief overview of the platform into which this solver was implemented, namely COOLFluiD, is also given at this point. Afterwards, a thorough numerical verification of the ACA approach has been made on an increasingly complex suite of test cases. The results obtained with ACA and Powell's source term implementations are given in order to numerically analyse and verify ACA and compare the two methods and validate them with the results from literature. The sixth chapter is devoted to further validation of ACA performed with a variety of more advanced space weather-related simulations. In this chapter, also the $vec{B}_{ extrm{0}} + vec{B}_{ extrm{1}}$ splitting technique used to treat planetary magnetosphere is presented along with its application to ACA and Powell's source term approaches. This technique is utilized in obtaining the solar wind/Earth's magnetosphere interaction results and is based on suppressing the direct inclusion of the Earth's magnetic field, which is a dipole field, in the solution variables. In this way, problems are avoided with the energy equation that could arise from the drastic change of the ratio of the dipole field and the variable field computed by the solver (i.e. $frac{lvertvec{B}_{ extrm{0}}lvert}{lvertvec{B}_{ extrm{1}}lvert}$) in the computational domain. Finally, conclusions and future perspectives related to the material presented in the thesis are put forward.

Ideal magnetohydrodynamics (MHD) simulations are known to have problems in satisfying the solenoidal constraint (i.e. the divergence of magnetic field should be equal to zero, $<p>ablacdotvec{B} = 0$). The simulations become unstable unless specific measures have been taken.<p><p>In this thesis, a solenoidal constraint satisfying technique that allows discrete satisfaction of the solenoidal constraint up to the machine accuracy is presented and validated with a variety of test cases. Due to its inspiration from Chorin's artificial compressibility method developed for incompressible CFD applications, the technique was named as \ / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished