Properties of stars undergoing pulsation such as the well known root-mean-density scaling relation can be useful when trying to match the observed properties of a particular star. It is often assumed that this relation is valid for p mode frequencies in rotating stars. To examine the change in frequency with rotation and mass, we have studied oscillation frequencies of two-dimensional uniformly rotating zero-age main sequence stellar models in the delta Scuti mass range. We identified 370 p and g axisymetric modes for non-rotating models and then traced the changes in their frequencies as the rotational velocity was increased. For each mass we considered a rotation sequence of ten models, with the largest rotation rate being about 200 km s$^{-1}$. We constrained the models to have the same surface shape, which can be characterized for uniform rotation by the ratio between the polar and the equatorial radii. We find that scaling relationships exist among the oscillation frequencies of the same mode for different masses when the models have the same shape. For p modes, this scaling closely follows the period root-mean-density relation found in spherical stars. The g modes also scale between models of the same shape, with the scaling reflecting the change in properties outside the convective core as the stellar mass increases. These scaling relationships can be particularly useful in finding specific stellar models to match the oscillation frequencies of individual stars.

Three dimensional hydrodynamic simulations of full amplitude RR Lyrae stars have been computed for several models across the instability strip. The three dimensional nature of the calculations allows convection to be treated without reference to a phenomenological approach such as the local mixing length theory. Specifically, the time dependent interaction of the large scale eddies and the radial pulsation is controlled by the conservation laws, while the effects of smaller convective eddies are simulated by an eddy viscosity model. The light amplitudes for these calculations are quite similar to those of our previous two dimensional calculations in the middle of the instability strip, but somewhat lower near the red edge, the fundamental blue edge, and for the one first overtone model we computed. The time dependent interaction between the radial pulsation and the convective energy transport is essentially the same in three dimensions as it is in two dimensions. There are some differences between the light curves between the two and three dimensional simulations, particularly during decreasing light. Reasons for the differences, both numerical and physical are explored.

To fully understand the Be star phenomenon, one must have a reasonable degree of knowledge about the star beneath the disk, which is often found to be rapidly rotating. Rapid rotation complicates modelling because fundamental properties like the stellar luminosity and effective temperature require knowledge of the angle of inclina- tion at which the star is observed. Furthermore our knowledge of the structure of rapidly rotating stars is on a less sure foundation than for non-rotating stars. The uncertainties in the inclination and the surface properties of a few rapidly rotating stars have been substantially reduced by interferometric observations over the last decade, and these stars can be used as tests of rotating stellar models, even if those stars themselves may not be Be stars. Vega, as an MK standard, is historically a very important star because it is used for calibration purposes. However, several studies have suggested that Vega is a rapidly rotating star viewed at a very low inclination angle, raising questions as to how well we really know its properties. Appropriate modelling has been challenging and there is still room for debate over the actual properties of Vega, as opposed to its observed properties. We have previously shown that under certain conditions both the stellar surface properties and the deduced surface properties scale from one model to another with the same surface shape. We used this scaling algorithm with realistic 2D models to compute high-resolution spectral energy distributions and interferometric visibilities to determine the best rotating model fit to Vega. Detailed comparisons between the computed and observed data will be presented.

The spectral energy distribution as a function of inclination is computed using 2D rotating stellar models and NLTE plane parallel stellar atmospheres. These models cover the range from $1.875M_{\odot}$ to $3.0M_{\odot}$. The deduced effective temperature is determined by B-V computed from the spectral energy distribution, and the deduced luminosity is computed as the integral of the spectral energy distribution over all frequencies, assuming the distance and reddening are known. These deduced quantities are obtained from the observed spectral energy distribution assuming the objects are spherically symmetric, and thus the results are dependent on the inclination. Previous work has shown that the surface properties between two rotating stellar models with the same surface shape scale, and this is also true for the deduced effective temperature and luminosity over this limited mass range.

We have developed a multidimensional radiation hydrodynamics code to simulate the interaction of radial stellar pulsation and convection for full amplitude pulsating models. Convection is computed using large eddy simulations. Here we perform three-dimensional simulations of RR Lyrae stars for comparison with previously reported two-dimensional simulations. We find that the time dependent behavior of the peak convective flux on pulsation phase is very similar in both the two-dimensional and three-dimensional calculations. The growth rates of the pulsation in the two-dimensional calculations are about $0.1\%$ higher than in the three-dimensional calculations.The amplitude of the light curve for a 6500 K RR Lyrae model is essentially the same for our 2D and 3D calculations, as is the rising light curve. There are differences in slope at various times during falling light.

We have developed a three-dimensional radiation hydrodynamics code to simulate the interaction of convection and radial pulsation in classical variable stars. One key goal is the ability to carry these simulations to full amplitude in order to compare them with observed light curves. Previous multi-dimensional calculations were prevented from reaching full amplitude because of drift in the radial coordinate system, due to the algorithm defining radial movement of the coordinate system during the pulsation cycle. We have removed this difficulty by defining our radial coordinate flow algorithm to require that the mass in a spherical shell remains constant for every time-step throughout the pulsation cycle. We have used our new code to perform 2D simulations of the interaction of radial pulsation and convection. We have made comparisons between light curves from our 2D convective simulations with observed light curves and find that our 2D simulated light curves are better able to match the observed light curve shape near the red edge of the RR Lyrae instability strip than light curves from previous 1D time dependent convective models.

Spectral energy distributions are computed using 2D rotating stellar models and NLTE plane parallel model atmospheres. A rotating, 2D stellar model has been found which matches the observed ultraviolet and visible spectrum of ? Oph. The SED match occurs for the interferometrically deduced surface shape and inclination, and is different from the SED produced by spherical models. The p mode oscillation frequencies in which the latitudinal variation is modelled by a linear combination of eight Legendre polynomials were computed for this model. The five highest and seven of the nine highest amplitude modes show agreement between computed axisymmetric, equatorially symmetric mode frequencies and the mode frequencies observed by MOST to within the observational error. Including nonaxisymmetric modes up through |m| = 2 and allowing the possibility that the eight lowest amplitude modes could be produced by modes which are not equatorially symmetric produces matches for 24 out of the 35 MOST modes to within the observational error and another eight modes to within twice the observational error. The remaining three observed modes can be fit within 4.2 times the observational error, but even these may be fit to within the observational error if the criteria for computed modes are expanded.

A rotating, two dimensional stellar model is evolved to match the approximate conditions of Alpha Oph. Both axisymmetric and nonaxisymmetric oscillation frequencies are computed for 2D rotating models which approximate the properties of Alpha Oph. These computed frequencies are compared to the observed frequencies. Oscillation calculations are made assuming the eigenfunction can be fit with six Legendre polynomials, but comparison calculations with eight Legendre polynomials show the frequencies agree to within about 0.26% on average. The surface horizontal shape of the eigenfunctions for the two sets of assumed number of Legendre polynomials agrees less well, but all calculations show significant departures from that of a single Legendre polynomial. It is still possible to determine the large separation, although the small separation is more complicated to estimate. With the addition of the nonaxisymmetric modes with |m| less than or equal 4, the frequency space becomes sufficiently dense that it is difficult to comment on the adequacy of the fit of the computed to the observed frequencies. While the nonaxisymmetric frequency mode splitting is no longer uniform, the frequency difference between the frequencies for positive and negative values of the same m remains 2m times the rotation rate.

We discuss how the p- and g-mode frequencies calculated for rotating stars are affected by the way of including the centrifugal deformations. We find that p-mode frequencies are sensitive to the accuracy in treating the centrifugal deformation, while the effect becomes smaller for smaller frequencies in the g-mode range.

Zero age main sequence models of uniformly rotating stars have been computed for ten masses between 1.625 and 8 M_\odot and 21 rotation rates from zero to nearly critical rotation. The surface shape is used to distinguish rotation rather than the surface equatorial velocity or the rotation rate. Using the surface shape is close to, but not quite equivalent to using the ratio of the rotation rate to the critical rotation rate. Using constant shape as the rotation variable means that it and the mass are separable, something that is not true for either the rotation rate or surface equatorial velocity. Thus a number of properties, including the ratio of the effective temperature anywhere on the surface to the equatorial temperature, are nearly independent of the mass of the model, as long as the rotation rate changes in such a way to keep the surface shape constant.

We are developing a 3D radiation hydrodynamics code to simulate the interaction of convection and pulsation in classical variable stars. One key goal is the ability to carry these simulations to full amplitude in order to compare them with observed light and velocity curves. Previous 2D calculations were prevented from doing this because of drift in the radial coordinate system, due to the algorithm defining radial movement of the coordinate system during the pulsation cycle. We remove this difficulty by defining our coordinate system flow algorithm to require that the mass in a spherical shell remain constant throughout the pulsation cycle. We perform adiabatic test calculations to show that large amplitude solutions repeat over more than 150 pulsation periods. We also verify that the computational method conserves the peak kinetic energy per period, as must be true for adiabatic pulsation models.

Mode splittings produced by uniform rotation and a particular form of differential rotation are computed for two-dimensional rotating 10 Mo ZAMS stellar models. The change in the character of the mode splitting is traced as a function of uniform rotation rate, and it is found that only relatively slow rotation rates are required before the mode splitting becomes asymmetric about the azimuthally symmetric (m=0) mode. Increased rotation produces a progressively altered pattern of the individual modes with respect to each other. Large mode splittings begin to overlap with the mode splittings produced by different radial and latitudinal modes at relatively low rotation rates. The mode splitting pattern for the differentially rotating stars we model is different than that for uniformly rotating stars, making the mode splitting a possible discriminant of the internal angular momentum distribution if one assumes the formidable challenge of mode identification can be overcome.

We have investigated the effects of uniform rotation and a specific model for differential rotation on the pulsation frequencies of 10 \Msun\ stellar models. Uniform rotation decreases the frequencies for all modes. Differential rotation does not appear to have a significant effect on the frequencies, except for the most extreme differentially rotating models. In all cases, the large and small separations show the effects of rotation at lower velocities than do the individual frequencies. Unfortunately, to a certain extent, differential rotation mimics the effects o f more rapid rotation, and only the presence of some specific observed frequencies with well identified modes will be able to uniquely constrain the internal rotation of pulsating stars.

Spectral energy distributions for models of arbitrarily rotating stars are computed using two dimensional rotating stellar models, NLTE plane parallel model atmospheres, and a code to integrate the appropriately weighted intensities over the visible surface of the stellar disk. The spectral energy distributions depend on the inclination angle between the observer and the rotation axis of the model. We use these curves to deduce what one would infer the model's luminosity and effective temperature to be assuming the object was nonrotating.

Radial and nonradial oscillations offer the opportunity to investigate the interior properties of stars. We use 2D stellar models and a 2D finite difference integration of the linearized pulsation equations to calculate non-radial oscillations. This approach allows us to directly calculate the pulsation modes for a distorted rotating star without treating the rotation as a perturbation. We are also able to express the finite difference solution in the horizontal direction as a sum of multiple spherical harmonics for any given mode. Using these methods, we have investigated the effects of increasing rotation and the number of spherical harmonics on the calculated eigenfrequencies and eigenfunctions and compared the results to perturbation theory. In slowly rotating stars, current methods work well, and we show that the eigenfunction can be accurately modelled using 2nd order perturbation theory and a single spherical harmonic. We use 10 Msun models with velocities ranging from 0 to 420 km/s (0.89 Omega_c) and examine low order p modes. We find that one spherical harmonic remains reasonable up to a rotation rate around 300km s^{-1} (0.69 Omega_c) for the radial fundamental mode, but can fail at rotation rates as low as 90 km/s (0.23 Omega_c) for the 2H mode or l = 2 p_2 mode, based on the eigenfrequencies alone. Depending on the mode in question, a single spherical harmonic may fail at lower rotation rates if the shape of the eigenfunction is taken into consideration. Perturbation theory, in contrast, remains valid up to relatively high rotation rates for most modes. We find the lowest failure surface equatorial velocity is 120 km/s (0.30 Omega_c) for the l = 2 p_2 mode, but failure velocities between 240 and 300 km/s (0.58-0.69 Omega_c)are more typical.

The spectral energy distribution (SED) of a non-spherical star could differ significantly from the SED of a spherical star with the same average temperature and luminosity. Calculation of the SED of a deformed star is often approximated as a composite of several spectra, each produced by a plane parallel model of given effective temperature and gravity. The weighting of these spectra over the stellar surface, and hence the inferred effective temperature and luminosity, will be dependent on the inclination of the rotation axis of the star with respect to the observer, as well as the temperature and gravity distribution on the stellar surface. Here we calculate the surface conditions of rapidly rotating stars with a 2D stellar structure and evolution code and compare the effective temperature distribution to that predicted by von Zeipel's law. We calculate the composite spectrum for a deformed star by interpolating within a grid of intensity spectra of plane parallel model atmospheres and integrating over the surface of the star. Using this method, we find that the deduced variation of effective temperature with inclination can be as much as 3000 K for an early B star, depending on the details of the underlying model.

The structure and evolution of close binary stars has been studied using the two-dimensional (2D) stellar structure algorithm developed by Deupree (1995). We have calculated a series of solar composition stellar evolution sequences of binary models, where the mass of the 2D model is 8Msun with a point-mass 5Msun companion. We have also studied the structure of the companion in 2D, by considering the zero-age main-sequence (ZAMS) structure of a 5Msun model with an 8Msun point-mass companion. In all cases the binary orbit was assumed to be circular and co-rotating with the rotation rate of the stars. We considered binary models with three different initial separations, a = 10, 14 and 20Rsun. These models were evolved through central hydrogen burning or until the more massive star expanded to fill its critical potential surface or Roche lobe. The calculations show that evolution of the deep interior quantities is only slightly modified from those of single star evolution. Describing the model surface as a Roche equipotential is also satisfactory until very close to the time of Roche lobe overflow, when the self gravity of the model about to lose mass develops a noticeable aspherical component and the surface time scale becomes sufficiently short that it is conceivable that the actual surface is not an equipotential.