Finite element modelling of perturbed stellar systems
arXiv:1001.3086 · doi:10.1111/j.1365-2966.2010.16365.x
Abstract
I formulate a general finite element method (FEM) for self-gravitating stellar systems. I split the configuration space to finite elements, and express the potential and density functions over each element in terms of their nodal values and suitable interpolating functions. General expressions are then introduced for the Hamiltonian and phase space distribution functions of the stars that visit a given element. Using the weighted residual form of Poisson's equation, I derive the Galerkin projection of the perturbed collisionless Boltzmann equation, and assemble the global evolutionary equations of nodal distribution functions. The FEM is highly adaptable to all kinds of potential and density profiles, and it can deal with density clumps and initially non-axisymmetric systems. I use ring elements of non-uniform widths, choose linear and quadratic interpolation functions in the radial direction, and apply the FEM to the stability analysis of the cutout Mestel disc. I also integrate the forced evolutionary equations and investigate the disturbances of a stable stellar disc due to the gravitational field of a distant satellite galaxy. The performance of the FEM and its prospects are discussed.
11 pages, three figures, accepted for publication by MNRAS
References in corpus (8)
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
- SAPPORO: A way to turn your graphics cards into a GRAPE-6
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Units
- Unstable Disk Galaxies. I. Modal Properties
- A parallel gravitational N-body kernel
- New biorthogonal potential--density basis functions
- Finite element modelling of perturbed stellar systems
Cited by in corpus (7)
- Relaxation in self-gravitating systems
- Density waves in debris discs and galactic nuclei
- Spiral eigenmodes triggered by grooves in the phase space of disc galaxies
- Finite element modelling of perturbed stellar systems
- Global drag-induced instabilities in protoplanetary disks
- Bisymmetric normal modes in soft-centred and realistic galactic discs
- Measuring Limb Darkening of Stars in high magnification Microlensing Events by the Finite Element Method