Toroidal figures of equilibrium from a 2nd-order accurate, accelerated SCF-method with subgrid approach

Abstract

We compute the structure of a self-gravitating torus with polytropic equation-of-state (EOS) rotating in an imposed centrifugal potential. The Poisson-solver is based on isotropic multigrid with optimal covering factor (fluid section-to-grid area ratio). We work at $2$nd-order in the grid resolution for both finite difference and quadrature schemes. For soft EOS (i.e. polytropic index $n\ge 1$), the underlying $2$nd-order is naturally recovered for Boundary Values (BVs) and any other integrated quantity sensitive to the mass density (mass, angular momentum, volume, Virial Parameter, etc.), i.e. errors vary with the number $N$ of nodes per direction as $\sim 1/N^2$. This is, however, not observed for purely geometrical quantities (surface area, meridional section area, volume), unless a subgrid approach is considered (i.e. boundary detection). Equilibrium sequences are also much better described, especially close to critical rotation. Yet another technical effort is required for hard EOS ($n<1$), due to infinite mass density gradients at the fluid surface. We fix the problem by using kernel splitting. Finally, we propose an accelerated version of the SCF-algorithm based on a node-by-node pre-conditionning of the mass density at each step. The computing time is reduced by a factor $2$ typically, regardless of the polytropic index. There is a priori no obstacle to applying these results and techniques to ellipsoidal configurations and even to $3$D-configurations.

arXiv commentaccepted for publication; 18 pages; 26 figures

References in corpus (16)

• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• Mapping the radial structure of AGN tori
• An algorithm for computing the 2D structure of fast rotating stars
• Pseudo-Newtonian and general relativistic barotropic tori in Schwarzschild-de Sitter spacetimes
• A multi-domain spectral method for scalar and vectorial Poisson equations with non-compact sources
• Composite Polytrope Models of Molecular Clouds. I. Theory
• Uniformly rotating axisymmetric fluid configurations bifurcating from highly flattened Maclaurin spheroids
• Charged tori in spherical gravitational and dipolar magnetic fields
• Formation of the Ring-like Structure in the SN 1987A Nebula due to the Magnetic Pressure of the Toroidal Field
• Slicing the Torus: Obscuring Structures in Quasars
• Equilibrium sequences of irrotational binary polytropic stars : The case of double polytropic stars
• Effect of the toroidal magnetic field on the runaway instability of relativistic tori
• Uniformly Rotating Polytropic Rings in Newtonian Gravity
• Uniformly Rotating Homogeneous and Polytropic Rings in Newtonian Gravity
• Analytical approximation for the structure of differentially rotating barotropes
• Relativistic Singular Isothermal Toroids

Cited by in corpus (5)

• The equilibrium of over-pressurised polytropes
• Rigidly rotating, incompressible spheroid-ring systems: new bifurcations, critical rotations and degenerate states
• Analytic construction of baroclinic tori by non-linear summation
• Nested spheroidal figures of equilibrium -- IV. On heterogeneous configurations
• Impact of rotation on synthetic mass-radius relationships of two-layer rocky planets and water worlds