We carry out 3-D numerical simulations of the dynamical instability in rapidly rotating stars initially modeled as polytropes with n = 1.5, 1.0, and 0.5. The calculations are done with a SPH code using Newtonian gravity, and the gravitational radiation is calculated in the quadrupole limit. All models develop the global m=2 bar mode, with mass and angular momentum being shed from the ends of the bar in two trailing spiral arms. The models then undergo successive episodes of core recontraction and spiral arm ejection, with the number of these episodes increasing as n decreases: this results in longer-lived gravitational wave signals for stiffer models. This instability may operate in a stellar core that has expended its nuclear fuel and is prevented from further collapse due to centrifugal forces. The actual values of the gravitational radiation amplitudes and frequencies depend sensitively on the radius of the star R_{eq} at which the instability develops.

We have carried out 3-D numerical simulations of the dynamical bar instability in a rotating star and the resulting gravitational radiation using both an Eulerian code written in cylindrical coordinates and a smooth particle hydrodynamics (SPH) code. The star is modeled initially as a polytrope with index $n = 3/2$ and $T_{\rm rot}/|W| \approx 0.30$, where $T_{\rm rot}$ is the rotational kinetic energy and $|W|$ is the gravitational potential energy. In both codes the gravitational field is purely Newtonian, and the gravitational radiation is calculated in the quadrupole approximation. We have run 3 simulations with the Eulerian code, varying the number of angular zones and the treatment of the boundary between the star and the vacuum. Using the SPH code we did 7 runs, varying the number of particles, the artificial viscosity, and the type of initial model. We compare the growth rate and rotation speed of the bar, the mass and angular momentum distributions, and the gravitational radiation quantities. We highlight the successes and difficulties of both methods, and make suggestions for future improvements.

We present the first calculations of the gravitational radiation produced by nonaxisymmetric dynamical instability in a rapidly rotating compact star. The star deforms into a bar shape, shedding $\sim 4\%$ of its mass and $\sim 17\%$ of its angular momentum. The gravitational radiation is calculated in the quadrupole approximation. For a mass $M \sim 1.4$ M$_{\odot}$ and radius $R \sim 10$ km, the gravitational waves have frequency $\sim 4$ kHz and amplitude $h \sim 2 \times 10^{-22}$ at the distance of the Virgo Cluster. They carry off energy $\Delta E/M \sim 0.1\%$ and radiate angular momentum $\Delta J/J \sim 0.7\%$.