A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

R. H.D. Townsend, S. P. Owocki, Asif Ud-Doula

Research output: Contribution to journalArticle

53 Citations (Scopus)

Abstract

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ∼3000 km s-1) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

Original languageEnglish (US)
Pages (from-to)139-157
Number of pages19
JournalMonthly Notices of the Royal Astronomical Society
Volume382
Issue number1
DOIs
StatePublished - Nov 1 2007

Fingerprint

massive stars
magnetospheres
magnetosphere
hydrodynamics
shock
plasma
modeling
stars
potential energy
simulation
one dimensional flow
shock heating
collision
divergence
centrifugal force
magnetic field
heating
cooling
rotating disks
loci

All Science Journal Classification (ASJC) codes

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

@article{b769aca13e0846b8a11de0ff6d276386,
title = "A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars",
abstract = "We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ∼3000 km s-1) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.",
author = "Townsend, {R. H.D.} and Owocki, {S. P.} and Asif Ud-Doula",
year = "2007",
month = "11",
day = "1",
doi = "10.1111/j.1365-2966.2007.12427.x",
language = "English (US)",
volume = "382",
pages = "139--157",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "1",

}

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars. / Townsend, R. H.D.; Owocki, S. P.; Ud-Doula, Asif.

In: Monthly Notices of the Royal Astronomical Society, Vol. 382, No. 1, 01.11.2007, p. 139-157.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

AU - Townsend, R. H.D.

AU - Owocki, S. P.

AU - Ud-Doula, Asif

PY - 2007/11/1

Y1 - 2007/11/1

N2 - We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ∼3000 km s-1) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

AB - We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ∼3000 km s-1) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

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

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

U2 - 10.1111/j.1365-2966.2007.12427.x

DO - 10.1111/j.1365-2966.2007.12427.x

M3 - Article

VL - 382

SP - 139

EP - 157

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 1

ER -