Engineering topics
Publications and source records attributed to Duncan, M. J..
We present a progress report of a continuing study of the coupling between outer planetary system architecture and inner planetary system formation. Additional information is contained in the original extended abstract.
We have conducted a set of approximately 200 numerical experiments to test the hypothesis that a giant impact leading to the formation of Earth's Moon could have occurred tens of millions of years after most of the small debris in the inner Solar System had been incorporated into terrestrial planets or been removed from the region. More than half of these simulations ended with a giant impact between two of the five terrestrial planets that were initially present. Neglecting any rotational angular momentum prior to the collision, the merged planet typically has a rotation period of less than five hours. The mean planetary obliquity is 91.7 degrees, and the median is 87.9 degrees; thus, there is no statistically significant difference between the number of bodies with prograde rotation and the number with retrograde rotation. There is a paucity of planets with obliquity close to 90 degrees, but the total number of impacts was too small for this result to be of much significance. Several encounters leading to collisions are dominated by three-body effects, with the velocity at impact being slightly less than the free-space escape velocity of the two bodies; the obliquity distribution produced by these impacts appears to be random.
The Giant Impact Hypothesis was introduced in the mid-1970's after consideration of results from the Apollo Moon missions. This hypothesis best explains the similarity in elemental proportions in lunar and terrestrial rocks, the depletion of lunar volatiles, the lack of lunar iron. and the large angular momentum in the Earth-Moon system. Comparison between the radiometric ages of inclusions in the most primitive meteorites and those of inclusions in the oldest lunar rocks and the differentiation age of Earth suggests that the Earth-Moon system formed about 100 Myr after the oldest meteorites. In addition, the age of the famous Martian meteorite ALH84001 and an early solidification time estimated from the Martian crust, suggest that the inner Solar System was fairly clear of large bodies about 10 Myr after the oldest meteorites formed. Thus, the 'standard model' suggests that for a period of several tens of millions of years the terrestrial planet region had few. if any, lunar-sized bodies and there were five terrestrial planets, Mercury, Venus, the two progenitors of the Earth-Moon system, and Mars. To simulate the dynamics of the Solar System before the hypothesized Moon-forming impact, we are integrating the Solar System with the Earth-Moon system replaced by two bodies in heliocentric orbits between Venus and Mars. The total (orbital) angular momentum of the Earth-Moon progenitors is that of the present Earth-Moon system, and their total mass is that of the Earth-Moon system. We are looking at ranges in mass ratio and initial values for eccentricity, inclination. and semi-major axis. We are using the SYMBA integrator to integrate these systems until a collision occurs or a time of 200 Myr elapses. Results are presented in a companion paper.
The giant impact hypothesis was introduced in the mid-1970s after consideration of results from the Apollo missions. This hypothesis best explains the similarity in elemental proportions in lunar and terrestrial rocks, the depletion of lunar volatiles, the lack of lunar Fe, and the large angular momentum in the Earth-Moon system. Comparison between the radiometric ages of inclusions in the most primitive meteorites and in the oldest lunar rocks and the differentiation age of Earth suggests that the Earth-Moon system formed about100 m.y. after the oldest meteorites. In addition, the age of the famous martian meteorite ALH 84001 and an early Martian solidification time obtained by Lee and Halliday suggest that the inner solar system was fairly clear of large bodies about 10 m.y. after the oldest meteorites formed. Thus, the "standard model" suggests that for several tens of millions of years, the terrestrial planet region had few, if any, lunar-sized bodies, and there were five terrestrial planets: Mercury, Venus, the two progenitors of the Earth-Moon system, and Mars. To simulate the dynamics of the solar system before the hypothesized Moon-forming impact, we are integrating the solar system with the Earth-Moon system replaced by two bodies in heliocentric orbits between Venus and Mars. The total (orbital) angular momentum of the Earth-Moon progenitors is that of the present Earth-Moon system, and their total mass is that of the Earth-Moon System. We are looking at ranges in mass ratio and initial values for eccentricity, inclination, and semimajor axis. We are using the SYMBA integrator to integrate these systems until a collision occurs or a time of 200 m.y. elapses. Results are presented in a companion abstract, (also presented at this meeting).
Substantial evidence indicates that the Earth-Moon system formed about 100 m.y. after the oldest meteorites and that the inner solar system had five terrestrial planets for several tens of millions of years before the hypothesized Moon-forming impact. We present and discuss some results from a series of N-body integrations in which the mass ratio of the Earth-Moon progenitors is 8:1 or 1:1. We want to know if it is plausible to have the Earth-Moon progenitors collide between 8 m.y. and 200 m.y. after the other planets had formed and to have the resulting system look "similar" to the solar system. If a collision occurs, the integrations tell us which two bodies collide and the time of the collision. We also determine the angular momentum deficit (AMD) of the resulting terrestrial planets. Additionally, we calculate several parameters of the collision. We use the AMD of the terrestrial planets to compare the resulting system to our own. The AMD or a planet is the difference between its orbital angular momentum and its orbital angular momentum if it were in a circular orbit with zero inclination.
The process of planetary growth is extremely complicated, involving a myriad of physical and chemical processes, many of which are poorly understood. The ultimate configuration that a planetary system attains depends upon the properties of the disk out of which it grew, of the star at the center of the disk and, at least in some cases, of the interstellar environment. However, this dependence is poorly understood. Thus, in an effort to numerically survey the possible diversity of planetary systems, we have constructed synthetic systems of giant planets and integrated their orbits to determine the dynamical lifetimes and thus the viability of these systems. Our construction algorithm begins with 110 -- 180 planetesimals located between 4 and 40 AU from a one solar mass star; most initial planetesimals have masses several tenths that of Earth. We integrate the orbits of these bodies subject to mutual gravitational perturbations and gas drag for $10^6 - 10^7$ years, merging any pair of planetesimals which passed within one-tenth of a Hill Sphere of one another and adding "gas" to embryos larger than 10 Earth masses. Use of such large planetesimal radii provided sufficient damping to prevent the system from excessive dynamical heating. Subsequently, systems were evolved without gas drag, either with the inflated radii or with more realistic radii. Systems took from a few million years to greater than ten billion years to become stable ($10^9$ years without mergers of ejections). Some of the systems produced with the inflated radii closely resemble our Solar System. Encounters in simulations using realistic radii resulted in ejections, typically leaving only a few planets per system, most of which were in highly eccentric orbits. The structure and dynamics of the resulting "stable" systems is discussed in detail in the abstract by Levison et al.
The process of planetary growth is extremely complicated, involving a myriad of physical and chemical processes, many of which are poorly understood. The ultimate configuration that a planetary system attains depends upon the properties of the disk out of which it grew, of the star at the center of the disk and, at least in some cases, of the interstellar environment. In an effort to numerically survey the possible diversity of planetary systems, we have constructed synthetic systems of giant planets and integrated their orbits to determine the dynamical lifetimes and thus the viability of these systems. Our construction algorithm begins with 110 -- 180 planetesimals located between 4 and 40 AU from a one solar mass star; most initial planetesimals have masses several tenths that of Earth. We integrate the orbits of these bodies subject to mutual gravitational perturbations and -as drag for 10(exp 6) - 10(exp 7) years, merging any pair of planetesimals which pass within one-tenth of a Hill Sphere of one another and adding "gas" to embryos larger than 10 Earth masses. Use of such large planetesimal radii provided sufficient damping to prevent the system from excessive dynamical heating. Subsequently, systems were evolved without gas drag, either with the enlarged radii or with more realistic radii. Systems took from a few million years to greater than ten billion years to become stable (10(exp 9) years without mergers of ejections). Some of the systems produced with the enlarged radii closely resemble our outer Solar System. Many systems contained only Uranus-mass objects. Encounters in simulations using realistic radii resulted in ejections, typically leaving only a few planets per system, most of which were on very eccentric orbits. Some of the systems that we constructed were stable for at least a billion years despite undergoing macroscopic orbital changes on much shorter timescales.
Recent numerical and analytic studies of planetary orbits have demonstrated the importance of resonances and chaos in destabilizing planetary systems. Newton's "clockwork" description of regular, predictable planetary orbits has been replaced by a view in which many systems can have long but finite lifetimes. This new knowledge has altered our perceptions of the later stages of planetary growth and of the stability of planetary systems. Stability criteria are inexact and time dependent. Most previous studies have focused on the effects in initial planetary orbits on the stability of the system. We are conducting an investigation which focuses on the dependence of stability criteria on planetary mass. Synthetic systems are created by increasing the masses of the planets in our Solar System or of the moons of a particular planet; these systems are then integrated until orbit crossing occurs. We have found that over some ranges, the time until orbit crossing varies to a good approximation as a power clothe factor by which the masses of the secondaries arc increased; some scatter occurs as a consequence of vie chaotic nature of orbital evolution. The slope of this power law varies substantially from system to system, and for moons it is mildly dependent on the inclusion of the planet's quadrupole moment in the gravitational potential.
A collisionless distribution of stars with an anisotropic velocity distribution is considered as a model for the cores of spherical galaxies. The first moment of the corresponding Boltzmann equation is integrated to give the resulting density and velocity dispersion profiles. These distributions plus a central point source of light give a good fit to the data from recent observations of M87 without requiring a supermassive central object. This model requires one less fitting parameter than published models incorporating a central black hole. The same model with no point source of light also provides a good fit to the normal E1 galaxy NGC 3379. The success of this model in describing the cores of both M87 and NGC 3379 suggests that a supermassive black hole is not necessary in either.