Engineering topics
Publications and source records attributed to Krauss-Varban, Dietmar.
In this report, we provide a summary of our most significant research accomplishments resulting from this contract. For the sake of brevity, most of the projects are explained in a paragraph length, highlighting only pertinent results.
This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic-simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost supercomputers, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enables us to make a number of important discoveries predictions, which need to be looked for in satellite data. Examples include the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH), (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs), and (3) this finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.
This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost desktop superconductors, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enabled us to make a number of important discoveries/predictions, which need to be looked for in satellite data. Examples include: (1) the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH;) (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs;) and (3) the finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.
The work carried out under this contract attempts a better understanding of whistler wave generation and associated scattering of electrons in the solar wind. This task is accomplished through simulations using a particle-in-cell code and a Vlasov code. In addition, the work is supported by the utilization of a linear kinetic dispersion solver. Previously, we have concentrated on gaining a better understanding of the linear mode properties, and have tested the simulation codes within a known parameter regime. We are now in a new phase in which we implement, execute, and analyze production simulations. This phase is projected to last over several reporting periods, with this being the second cycle. In addition, we have started to research to what extent the evolution of the pertinent instabilities is two-dimensional. We are also continuing our work on the visualization aspects of the simulation results, and on a code version that runs on single-user Alpha-processor based workstations.
This report describes the work done by SciberNet, Inc. during the month of January. During this time, we primarily worked on further analysis of the results presented at the AGU as well as writing them up for publication. Using large scale simulations, we showed that the magnetopause during the southward IMF case is quite irregular with varying thickness, and has a complex flow pattern owing to the nonlinear effects of the convective flow superimposed on the flows generated in the reconnection layer. We used inflow-outflow boundary conditions to examine the kinetic nature of the discontinuities that are formed in the reconnection layer and concluded that nonlocal effects play a major role in the formation of such discontinuities and can alter their properties from the usual structures expected from 1-D simulations or from fluid theories. Finally, we used our 3-D simulations to examine the nonlinear interaction of the tearing mode with the Kelvin-Helmholtz instability. We showed that this interaction leads to the generation of a large core field which is observed both in the magnetotail as well as the magnetopause.
This report describes the work done by SciberNet, Inc. during the month of October. We are working on the further refinement of the model used in our large-scale hybrid simulations of the magnetopause. Specifically, we are experimenting with several ways of modeling the effects of cold magnetospheric ions into our simulations. In addition, we are preparing two presentations for the upcoming Fall AGU highlighting the results of these simulations. We have also made progress in our development of a new kinetic linear code which we are using to study the linear properties of the Kelvin-Helmholtz instability at the magnetopause. We have extended the code from the electrostatic limit to the fully electromagnetic regime and are currently in the process of debugging and testing the code. Finally, we have made several test runs with our 2-D hybrid code for the magnetopause. The inflow-outflow boundary conditions are working properly. However, there are issues related to the setup and evolution of the original equilibrium that we are still trying to resolve. Finally, we are preparing several presentations for the upcoming Fall AGU.
This report describes the work done by SciberNet, Inc. during the month of August. We have resolved the issues associated with the implementation of the dipole field in our large scale hybrid simulations of the magnetopause. We have setup several runs and will spend the next several months analyzing the data. The results will be presented at the Fall AGU. We are also continuing our analysis of the 3-D simulations of thin current sheets at the magnetopause, paying special attention to the conditions under which Kelvin-Helmholtz would lead to sizable perturbations of the magnetopause. In a related study, we are in the process of developing a new kinetic linear code that would for the first time enable us to examine the linear properties of the Kelvin-Helmholtz instability in the fully kinetic regime. Finally, we are continuing our code development to include inflow-outflow boundary conditions in our 2-D and 3-D hybrid codes. We are also comparing the different methods of code parallelization in order to extend the limits of our calculations.