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

J.-M. Huré 1 F. Hersant 1
1 ECLIPSE 2017
LAB - Laboratoire d'Astrophysique de Bordeaux [Pessac]
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.
Type de document :
Article dans une revue
Monthly Notices of the Royal Astronomical Society, Oxford University Press (OUP): Policy P - Oxford Open Option A, 2017, 464 (4), pp.4761 - 4776. <10.1093/mnras/stw2569>
Liste complète des métadonnées

https://hal.archives-ouvertes.fr/hal-01384481
Contributeur : Marie-Paule Pomies <>
Soumis le : jeudi 20 octobre 2016 - 09:24:43
Dernière modification le : mardi 6 décembre 2016 - 09:37:32

Identifiants

Collections

Citation

J.-M. Huré, F. Hersant. Toroidal figures of equilibrium from a second-order accurate, accelerated SCF method with subgrid approach. Monthly Notices of the Royal Astronomical Society, Oxford University Press (OUP): Policy P - Oxford Open Option A, 2017, 464 (4), pp.4761 - 4776. <10.1093/mnras/stw2569>. <hal-01384481>

Partager

Métriques

Consultations de la notice

33