HEAD 2003 Meeting
Session 19. Neutrinos, Gravitational Waves and Strong-field Gravity III
Poster, Sunday-Wednesday, March 23, 2003, Duration of Meeting

[Previous] | [Session 19] | [Next]

[19.04] Finite Element Analysis for General Relativistic Flow Environments

G. A. Richardson (NRC, MSFC), T. J. Chung (University of Alabama in Huntsville)

We present our finite element methodology for numerically solving the governing equations for general relativistic environments. Our motivation for the development of such a method is to study the environment around a rotating black hole, specifically the dynamics of the accretion disk and the associated formation of relativistic jets. The numerical technique we have developed is versatile enough to perform simulations ranging from relativistic shock capturing in magnetohydrodynamics to Poisson solvers for mapping fields. This numerical technique is unique in its ability to detect small variations in the physical parameters related to the flow-field or the magnetic field. Since this information is obtained during the computations, it can be used for adaptive mesh refinement, necessary variations in the PDE solver, the detection of numerical instabilities and potentially for the detection of turbulence. This direct link between the variations in the physical parameters and the numerical parameters helps to increase the accuracy and the stability of the system. While we have developed this method for use in finite element analysis, it is also possible to adapt it for use with finite difference methods. We will present the methodology behind our numerical technique and some applications.

The National Research Council research associateship program supported this work.

[Previous] | [Session 19] | [Next]

Bulletin of the American Astronomical Society, 35#2
© 2003. The American Astronomical Soceity.