Magnetically controlled spasmodic accretion during star formation. I. Formulation of the problem and method of solution

Konstantinos Tassis, Telemachos Ch Mouschovias

Research output: Contribution to journalArticle

Abstract

We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10 8 cm-3. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 1011 cm-3. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star are calculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.

Original languageEnglish (US)
Pages (from-to)769-782
Number of pages14
JournalAstrophysical Journal
Volume618
Issue number2 I
DOIs
StatePublished - Jan 10 2005

Fingerprint

star formation
accretion
formulations
magnetic field
sinks
stars
fluid
fluids
approximation
grids
magnetic fields
magnetic disks
electron
ions
ion
integrators
cells
molecular ions
detachment
gravitational fields

Keywords

  • Accretion, accretion disks
  • Dust, extinction
  • MHD -shock waves
  • Magnetic fields
  • Stars: formation

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Magnetically controlled spasmodic accretion during star formation. I. Formulation of the problem and method of solution. / Tassis, Konstantinos; Mouschovias, Telemachos Ch.

In: Astrophysical Journal, Vol. 618, No. 2 I, 10.01.2005, p. 769-782.

Research output: Contribution to journalArticle

@article{7233643f940b4c639d153a847fe8dcfb,
title = "Magnetically controlled spasmodic accretion during star formation. I. Formulation of the problem and method of solution",
abstract = "We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10 8 cm-3. The {"}central sink approximation{"} is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 1011 cm-3. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star are calculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.",
keywords = "Accretion, accretion disks, Dust, extinction, MHD -shock waves, Magnetic fields, Stars: formation",
author = "Konstantinos Tassis and Mouschovias, {Telemachos Ch}",
year = "2005",
month = "1",
day = "10",
doi = "10.1086/424479",
language = "English (US)",
volume = "618",
pages = "769--782",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "2 I",

}

TY - JOUR

T1 - Magnetically controlled spasmodic accretion during star formation. I. Formulation of the problem and method of solution

AU - Tassis, Konstantinos

AU - Mouschovias, Telemachos Ch

PY - 2005/1/10

Y1 - 2005/1/10

N2 - We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10 8 cm-3. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 1011 cm-3. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star are calculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.

AB - We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10 8 cm-3. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 1011 cm-3. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star are calculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.

KW - Accretion, accretion disks

KW - Dust, extinction

KW - MHD -shock waves

KW - Magnetic fields

KW - Stars: formation

UR - http://www.scopus.com/inward/record.url?scp=14544304967&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=14544304967&partnerID=8YFLogxK

U2 - 10.1086/424479

DO - 10.1086/424479

M3 - Article

AN - SCOPUS:14544304967

VL - 618

SP - 769

EP - 782

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 2 I

ER -