Home News About Us Contact Contributors Disclaimer Privacy Policy Help FAQ

Quick Search
My eDoc
Session History
Support Wiki
Direct access to
document ID:

          Institute: MPI für Astronomie     Collection: Publikationen_mpia     Display Documents

ID: 660790.0, MPI für Astronomie / Publikationen_mpia
A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code
Authors:Mignone, A.; Flock, M.; Stute, M.; Kolb, S. M.; Muscianisi, G.
Date of Publication (YYYY-MM-DD):2012
Journal Abbrev.:Astronomy and Astrophysics
Start Page:id.A152
Audience:Not Specified
Abstract / Description:Context. Explicit numerical computations of hypersonic or super-fast differentially rotating disks are subject to the time-step constraint imposed by the Courant condition, according to which waves cannot travel more than a fraction of a cell during a single time-step update. When the bulk orbital velocity largely exceeds any other wave speed (e.g., sound or Alfvén), as computed in the rest frame, the time step is considerably reduced and an unusually large number of steps may be necessary to complete the computation. <BR /> Aims: We present a robust numerical scheme to overcome the Courant limitation by improving and extending the algorithm previously known as FARGO (fast advection in rotating gaseous objects) to the equations of magnetohydrodynamics (MHD) using a more general formalism. The proposed scheme conserves total angular momentum and energy to machine precision and works in cartesian, cylindrical, or spherical coordinates. The algorithm has been implemented in the next release of the PLUTO code for astrophysical gasdynamics and is suitable for local or global simulations of accretion or proto-planetary disk models. <BR /> Methods: By decomposing the total velocity into an average azimuthal contribution and a residual term, the algorithm approaches the solution of the MHD equations through two separate steps corresponding to a linear transport operator in the direction of orbital motion and a standard nonlinear solver applied to the MHD equations written in terms of the residual velocity. Since the former step is not subject to any stability restriction, the Courant condition is computed only in terms of the residual velocity, leading to substantially larger time steps. The magnetic field is advanced in time using the constrained transport method in order to fulfill the divergence-free condition. Furthermore, conservation of total energy and angular momentum is enforced at the discrete level by properly expressing the source terms in terms of upwind Godunov fluxes available during the standard solver. <BR /> Results: Our results show that applications of the proposed orbital-advection scheme to problems of astrophysical relevance provides, at reduced numerical cost, equally accurate and less dissipative results than standard time-marching schemes. Freely distributed at <A href="http://plutocode.ph.unito.it">http://plutocode.ph.unito.it</A>.
Free Keywords:methods: numerical; accretion; accretion disks; protoplanetary disks; magnetohydrodynamics (MHD); turbulence
External Publication Status:published
Document Type:Article
Communicated by:N. N.
Affiliations:MPI für Astronomie
The scope and number of records on eDoc is subject to the collection policies defined by each institute - see "info" button in the collection browse view.