ALBERT

Description: Current research carried out with the help of the ASEE-NASA Summer Faculty Program, at NASA-Ames, is concentrated on the dynamics of nuclear regions of galaxies. From a dynamical point of view a galaxy is a collection of around 10(sup 11) stars like our Sun, each of which moves in the summed gravitational field of all the remaining stars. Thus galaxy dynamics becomes a self-consistent n-body problem with forces given by Newtonian gravitation. Strong nonlinearity in the gravitational force and the inherent nonlinearity of self-consistent problems both argue for a numerical approach. The technique of numerical experiments consis of constructing an environment in the computer that is as close as possible to the physical conditions in a real galaxy and then carrying out experiments much like laboratory experiments in physics or engineering, in this environment. Computationally, an experiment is an initial value problem, and a good deal of thought and effort goes into the design of the starting conditions that serve as initial values. Experiments are run at Ames because all the 'equipment' is in place-the programs, the necessary computational power, and good facilities for post-run analysis. Our goal for this research program is to study the nuclear regions in detail and this means replacing most of the galaxy by a suitable boundary condition to allow the full capability of numerical experiments to be brought to bear on a small region perhaps 1/1000 of the linear dimensions of an entire galaxy. This is an extremely delicate numerical problem, one in which some small feature overlook, can easily lead to a collapse or blow-up of the entire system. All particles attract each other in gravitational problems, and the 1/r(sup 2) force is: (1) nonlinear; (2) strong at short range; (3) long-range, and (4) unscreened at any distance.

Description: Trajectory design of the orbit phase of the NEAR mission involves a new process that departs significantly from those procedures used in previous missions. In most cases, a precise spacecraft ephemeris is designed well in advance of arrival at the target body. For NEAR, the uncertainty in the dynamic environment around Eros does not allow the luxury of a precise spacecraft trajectory to be defined in advance. The principal cause of this uncertainty is the limited knowledge oi' the gravity field a,-id rotational state of Eros. As a result, the concept for the NEAR trajectory design is to define a number of rules for satisfying spacecraft, mission, and science constraints, and then apply these rules to various assumptions for the model of Eros. Nominal, high, and low Eros mass models are used for testing the trajectory design strategy and to bracket the ranges of parameter variations that are expected upon arrival at the asteroid. The final design is completed after arrival at Eros and determination of the actual gravity field and rotational state. As a result of the unplanned termination of the deep space rendezvous maneuver on December 20, 1998, the NEAR spacecraft passed within 3830 km of Eros on December 23, 1998. This flyby provided a brief glimpse of Eros, and allowed for a more accurate model of the rotational parameters and gravity field uncertainty. Furthermore, after the termination of the deep space rendezvous burn, contact with the spacecraft was lost and the NEAR spacecraft lost attitude control. During the subsequent gyrations of the spacecraft, hydrazine thruster firings were used to regain attitude control. This unplanned thruster activity used Much of the fuel margin allocated for the orbit phase. Consequently, minimizing fuel consumption is now even more important.

Description: Preliminary results of a continuing search microvariations in a selected sample.

Description: Recent studies suggest that galaxies can oscillate in normal modes with essentially no damping over a Hubble time. These modes may play an important role in the structure and evolution of disk/halo systems. Motivated by the possibility that normal mode oscillations exist in real galaxies, we are investigating the response of galactic disks to halo oscillations. The goal of these investigations is to ascertain whether or not observational signatures exist for such oscillations. Our approach is to perform numerical experiments on the response of a self-gravitating disk to a time-varying halo potential. We assume that a significant fraction of the mass in a galaxy is in a dark halo. The halo oscillates and the luminous disk material responds to these oscillations. Preliminary results are reported for disks embedded in a radially oscillating gravitational potential. The equilibrium initial disk is represented by an exponential density profile. Considerable care was taken to build an initial disk model that was "stable" over long time scales. A control experiment was run with the disk in a static halo potential. The disk responds to the time-varying potential by developing a ring structure, which forms and disappears during each halo oscillation cycle. The density of stars becomes depressed in an annular region at the radius where the disk epicycle frequency is equal to the halo oscillation period. This pattern of response persists over time periods approaching a Hubble time. In the oscillating potential, a bar develops in the inner disk. This bar is absent when the halo remains static. Specific targets of this study include the implications for large-scale disk structure, the gas dynamical response of the interstellar medium in such systems, and the inflow of material into the central regions of the galaxy.

Description: Orbit calculations suggest that the sparse globular cluster, Pal 5, will pass within 7 kpc of the Galactic center the next time it crosses the plane, where it might be destroyed by tidal stresses. We study this problem, treating Pal 5 as a self-consistent dynamical system orbiting through an external potential that represents the Galaxy. The first part of the problem is to find suitable analytic approximations to the Galactic potential. They must be valid in all regions the cluster is likely to explore. Observed velocity and positional data for Pal 5 are used as initial conditions to determine the orbit. Methods we used for a different problem some 12 years ago have been adapted to this problem. Three experiments have been run, with M/L= 1, 3, and 10, for the cluster model. The cluster blew up shortly after passing through the Galactic plane (about 130 Myrs after the beginning of the run) with M/L=1. At M/L = 3 and 10 the cluster survived, although it got quite a kick in the fundamental mode on passing through the plane. But the fundamental mode oscillation died out in a couple of oscillation cycles at M/L=10. Pal 5 will probably be destroyed on its next crossing of the Galactic plane if M/L=1, but it can survive (albeit with fairly heavy damage) if NI/L=3. We haven't tried to trap the mass limits more closely than that. Pal 5 comes through pretty well unscathed at M/L=10. An interesting follow-up experiment would be to back the cluster up along its orbit to look at its previous passage through the Galactic plane, to see what kind of object it might have been at earlier times.

Description: The key to studying central regions by means of nobody numerical experiments is to concentrate on the central few parsecs of a galaxy, replacing the remainder of the galaxy by a suitable boundary condition, rather after the manner in which stellar interiors can be studied without a detailed stellar atmosphere by replacing the atmosphere with a boundary condition. Replacements must be carefully designed because the long range gravitational force means that the core region is sensitive to mass outside that region and because particles can exchange between the outer galaxy and the core region. We use periodic boundary conditions, coupled with an iterative procedure to generate initial particle loads in isothermal equilibrium. Angular momentum conservation is ensured for problems including systematic rotation by a circular reflecting boundary and by integrating in a frame that rotates with the mean flow. Mass beyond the boundary contributes to the gravitational potential, but does not participate in the dynamics. A symplectic integration scheme has been developed for rotating coordinate systems. This combination works well, leading to robust configurations. Some preliminary results with this combination show that: (1) Rotating systems are extremely sensitive to non-axisymmetric external potentials, and (2) that a second core, orbiting near the main core (like the M31 second core system), shows extremely rapid orbital decay. The experimental setups will be discussed, along with preliminary results.

Description: Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Description: The investigation of stochastic particle acceleration through resonant interactions with plasma waves that populate the magnetosphere surrounding an accreting black hole is presented. Stochastic acceleration has been successfully applied to the problem of ion and electron energization in solar flares, and is capable of accounting for a wide range of both neutral and charged particle emissions. It is also a component in diffusive shock acceleration, since pitch-angle scattering (which is necessary for multiple shock crossings) is accompanied by diffusion in momentum space, which in turn yields a net systematic energy gain; however, stochastic energization will dominate the first-order shock process only in certain parameter regimes. Although stochastic acceleration has been applied to particle energization in the lobes of radio galaxies, its application to the central regions of AGNs (active galactic nuclei) has only recently been considered, but not in detail. We proposed to systematically investigate the plasma processes responsible for stochastic particle acceleration in black hole magnetospheres along with the energy-loss processes which impede particle energization. To this end, we calculated acceleration rates and escape time scales for protons and electrons resonating with Alfven waves, and for electrons resonating with whistlers. We also considered the "hot" topic of gamma-ray line emission from the Orion complex. We proposed that the observed gamma-ray lines are produced by energetic ions that are stochastically accelerated by cascading Alfven waves in the accretion plasma near a black hole. Related research papers that were published in journals are listed.

Description: MASTER (Mainbelt Asteroid Exploration/Rendezvous), a Discovery-class orbiter, will carry out a global geological and geochemical survey of Vesta. Primary goals include determining Vesta's mean density and interior structure through detailed mapping of the gravity field, imaging surface morphology at 3-m resolution, mapping mineralogy between 0.4 and 2.5 microns, and determining abundances of key elements through X-ray and gamma-ray spectroscopy. Spectroscopic evidence indicates the presence of basaltic lava flows on Vesta and suggests that impact basins may have exposed mantle materials. These possibilities, combined with the likelihood that Vesta is the ultimate source of HED meteorites, makes this asteroid an important target for a comprehensive orbiter mission. MASTER's global survey, combined with ongoing studies of HED meteorites, will be a major step in understanding the chemical, thermal, and geological evolution of Vesta. A particularly attractive opportunity involves a launch in June 2003 with arrival at Vesta in November 2009. The orbital phase of the mission is scheduled to last one year. A flyby of another mainbelt asteroid enroute to Vesta may be possible.