Simulation of Comet-Solar Wind Interaction in the AMROC Framework with the IGLM-MHD Scheme
摘要
In this work we present 2-dimensional, 1-fluid magnetohydrodynamic (MHD) simulations of the interaction of a comet with the solar wind. Even though our model can in principle run simulations with time-dependent physical parameters representing the comet travelling along its elliptic orbit, here we focus on “snapshots” at four heliocentric distances, namely: 1.3, 2.0, 2.7 and 3.25 AU. This strategy allows us to analyse some key features of the comet at relatively low computational cost. Our aim here is to present our implementation of such a scenario with the AMROC computational code and study the main features of the comet, which has the physical characteristics of 67P/Churyumov-Gerasimenko (67P/CG). The dependence of the physical parameters on the heliocentric distance is given by means of Parker’s model and of functions obtained from curve fitting of data collected in the literature. In addition, we model comets with three forms: circular, elliptic and bilobed (or “peanut-shaped”.) AMROC is in development and the model shown here is an valuable contribution to these efforts, as it represents a new kind of problem implemented in this code. Furthermore, our model provides the possibility of creation of new and interesting scenarios with regard to comets and other sources of plasma.