We have developed a parallel, simple, and fast hydrodynamics code for multi-dimensional, self-gravitating, adiabatic flows. Our primary motivation is the study of the non-linear evolution of white dwarf oscillations excited via tidal resonances, typically over hundreds of stellar dynamical times. Consequently, we require long term stability, low diffusivity, and high algorithmic efficiency. An explicit, Eulerian, finite-difference scheme on a regular Cartesian grid fulfills these requirements. It provides uniform resolution throughout the flow, as well as simplifying the computation of the self-gravitational potential, which is done via spectral methods. In this paper, we describe the numerical scheme and present the results of some diagnostic problems. We also demonstrate the stability of a cold white dwarf in three dimensions over hundreds of dynamical times. Finally, we compare the results of the numerical scheme to the linear theory of adiabatic oscillations, finding numerical quality factors on the order of 6000, and excellent agreement with the oscillation frequency obtained by the linear analysis.

We consider the evolution of white dwarfs with compact object companions (specifically black holes with masses up to 10^6 solar masses, neutron stars, and other white dwarfs). We suppose that the orbits are initially quite elliptical and then shrink and circularise under the action of gravitational radiation. During this evolution, the white dwarfs will pass through resonances when harmonics of the orbital frequency match the stellar oscillation eigenfrequencies. As a star passes through these resonances, the associated modes will be excited and can be driven to amplitudes that are so large that there is a back reaction on the orbit which, in turn, limits the growth of the modes. A formalism is presented for describing this dynamical interaction for a non-rotating star in the linear approximation when the orbit can be treated as non-relativistic. A semi-analytical expression is found for computing the resonant energy transfer as a function of stellar and orbital parameters for the regime where back reaction may be neglected. This is used to calculate the results of passage through a sequence of resonances for several hypothetical systems. It is found that the amplitude of the l=m=2 f-mode can be driven into the non-linear regime for appropriate initial conditions. We also discuss where the no back reaction approximation is expected to fail, and the qualitative effects of back reaction.

We present a variational formalism for describing the dynamical evolution of an oscillating star with a point-mass companion in the linear, non-relativistic regime. This includes both the excitation of normal modes and the back-reaction of the modes on the orbit. The general formalism for arbitrary fluid configurations is presented, and then specialized to a homentropic potential flow. Our formalism explicitly identifies and conserves both energy and angular momentum. We also consider corrections to the orbit up to 7/2 post-Newtonian order.