共查询到20条相似文献，搜索用时 31 毫秒
1.
Keiji Ohtsuki 《Icarus》2004,172(2):432-445
We examine the rotation of a small moonlet embedded in planetary rings caused by impacts of ring particles, using analytic calculation and numerical orbital integration for the three-body problem. Taking into account the Rayleigh distribution of particles' orbital eccentricities and inclinations, we evaluate both systematic and random components of rotation, where the former arises from an average of a large number of small impacts and the latter is contribution from large impacts. Calculations for parameter values corresponding to inner parts of Saturn's rings show that a moonlet would spin slowly in the prograde direction if most impactors are small particles whose velocity dispersion is comparable to or smaller than the moonlet's escape velocity. However, we also find that the effect of the random component can be significant, if the velocity dispersion of particles is larger and/or impacts of large particles comparable to the moonlet's size are common: in this case, both prograde and retrograde rotations can be expected. In the case of a small moonlet embedded in planetary rings of equal-sized particles, we find that the systematic component dominates the moonlet rotation when m/M?0.1 (m and M are the mass of a particle and a moonlet, respectively), while the random component is dominant when m/M?0.3. We derive the condition for the random component to dominate moonlet rotation on the basis of our results of three-body orbital integration, and confirm agreement with N-body simulation. 相似文献
2.
We investigate the spin rates of moonlets embedded in planetary rings, subject to collisions with surrounding small particles, using three-body integrations including friction and spins. All successive impacts of the particle with the moonlet are followed, including a possible sliding phase after the initial inelastic rebounds. Two methods for treating impacts, (1) as instantaneous velocity changes and (2) using an impact force model, are applied after Salo (1995, Icarus 117, 287). Conducting a series of integrations with various initial summed spin velocity of the moonlet and the particle, we determine the equilibrium spin rate for which the averaged torque vanishes. This equilibrium spin rate corresponds to the final spin rate of the moonlet if the moonlet is much larger than the surrounding particles; it also corresponds to the mean spin rate for a ring composed of identical particles. We find that the equilibrium spin rate is enhanced by sliding orbits as compared with the spin rate determined by considering only the first impacts of the particles with the moonlet. If the random velocities of incident particles are small enough, the resulting equilibrium spin rate of the moonlet can be larger than the synchronous rotation rate, for rp∼1, where rp denotes the sum of radii of the colliding pair normalized by their mutual Hill radius. In this special case aggregates without internal strength may become rotationally unstable. However, the equilibrium spin rate decreases with increasing random velocity, and aggregates are always rotationally stable in the more likely case where the relative velocities are comparable to the mutual escape velocity. We also compare our results with the mean spin rates found in previous N-body simulations, and find a good agreement for optically thin rings; however the spin rates for optically thick rings are significantly larger than those predicted by our three-body calculations. 相似文献
3.
As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 107 years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System. 相似文献
4.
This paper presents the results of N-body simulations of moonlets embedded in broad rings, focusing specifically on the saturnian A ring. This work adds to previous efforts by including particle self-gravity and particle size distributions. The discussion here focuses primarily on the features that form in the background particles as a result of the moonlet. Particle self-gravity tends to damp out features produced by embedded moonlets and this damping is enhanced if the moonlet is simply the largest member of a continuous size distribution. Observable features around an embedded moonlet appear to require that the largest ring particles be no more massive than 1/30 the mass of the moonlet. These results, compared with current and future Cassini observations, will provide insight into the nature of the particle population in the saturnian rings. Some time is also spent analyzing the way in which the background particles cluster around the moonlet. The accretion of small particles onto the moonlet can be limited by disruptive collisions with the largest ring particles in the particle size distribution. 相似文献
5.
We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 105 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation. 相似文献
6.
We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular “atmosphere” that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate. 相似文献
7.
We present a numerical method that incorporates particle sticking in simulations using the N-body code pkdgrav to study motions in a local rotating frame, such as a patch of a planetary ring. Particles stick to form non-deformable but breakable aggregates that obey the (Eulerian) equations of rigid-body motion. Applications include local simulations of planetary ring dynamics and planet formation, which typically feature hundreds of thousands or more colliding bodies. Bonding and breaking thresholds are tunable parameters that can approximately mimic, for example, van der Waals forces or interlocking of surface frost layers. The bonding and breaking model does not incorporate a rigorous treatment of internal fracture; rather the method serves as motivation for first-order investigation of how semi-rigid bonding affects the evolution of particle assemblies in high-density environments.We apply the method to Saturn’s A ring, for which laboratory experiments suggest that interpenetration of thin, frost-coated surface layers may lead to weak cohesive bonding. These experiments show that frost-coated icy bodies can bond at the low impact speeds characteristic of the rings. Our investigation is further motivated by recent simulations that suggest a very low coefficient of restitution is needed to explain the amplitude of the azimuthal brightness asymmetry in Saturn’s A ring, and the hypothesis that fine structure in Saturn’s B ring may in part be caused by large-scale cohesion.This work presents the full implementation of our model in pkdgrav, as well as results from initial tests with a limited set of parameters explored. We find a combination of parameters that yields aggregate size distribution and maximum radius values in agreement with Voyager data for ring particles in Saturn’s outer A ring. We also find that the bonding and breaking parameters define two strength regimes in which fragmentation is dominated either by collisions or other stresses, such as tides. We conclude our study with a discussion of future applications of and refinements to our model. 相似文献
8.
Douglas C. Heggie 《Celestial Mechanics and Dynamical Astronomy》1974,10(2):217-241
The work of Aarseth and Zare (1974) is extended to provide aglobal regularisation of the classical gravitational three-body problem: by transformation of the variables in a way that does not depend on the particular configuration, we obtain equations of motion which are regular with respect to collisions between any pair of particles. The only cases excepted are those in which collisions between more than one pair occur simultaneously and those in which at least one of the masses vanishes. However, by means of the same principles the restricted problem is regularised globally if collisions between the two primaries are excluded. Results of numerical tests are summarised, and the theory is generalised to provide global regularisations, first, for perturbed three-body motion and, second, for theN-body problem. A way of increasing the number of degrees of freedom of a dynamical system is central to the method, and is the subject of an Appendix. 相似文献
9.
Gravitational accretion in the rings of Saturn is studied with local N-body simulations, taking into account the dissipative impacts and gravitational forces between particles. Common estimates of accretion assume that gravitational sticking takes place beyond a certain distance (Roche distance) where the self-gravity between a pair of ring particles exceeds the disrupting tidal force of the central object, the exact value of this distance depending on the ring particles' internal density. However, the actual physical situation in the rings is more complicated, the growth and stability of the particle groups being affected also by the elasticity and friction in particle impacts, both directly via sticking probabilities and indirectly via velocity dispersion, as well as by the shape, rotational state and the internal packing density of the forming particle groups. These factors are most conveniently taken into account via N-body simulations. In our standard simulation case of identical 1 m particles with internal density of solid ice, ρ=900 kg m−3, following the Bridges et al., 1984 elasticity law, we find accretion beyond a=137,000-146,000 km, the smaller value referring to a distance where transient aggregates are first obtained, and the larger value to the distance where stable aggregates eventually form in every experiment lasting 50 orbital periods. Practically the same result is obtained for a constant coefficient of restitution εn=0.5. In terms of rp parameter, the sum of particle radii normalized by their mutual Hill radius, the above limit for perfect accretion corresponds to rp<0.84. Increased dissipation (εn=0.1), or inclusion of friction (tangential force 10% of normal force) shifts the accretion region inward by about 5000 km. Accretion is also more efficient in the case of size distribution: with a q=3 power law extending over a mass range of 1000, accretion shifts inward by almost 10,000 km. The aggregates forming in simulations via gradual accumulation of particles are synchronously rotating. 相似文献
10.
The viscosity (the angular momentum flux) in the disk of mutually gravitating particles of Saturn's rings is investigated. The hydrodynamic theory of the gravitational Jeans-type instability of small gravity perturbations (e.g., those produced by spontaneous disturbances) of the disk is developed. It is suggested that in such a system the hydrodynamic turbulence may arise as a result of the instability. The turbulence is related to stochastic motions of “fluid” elements. The objective of this paper is to show that in the Jeans-unstable Saturnian ring disk the turbulent viscosity may exceed the ordinary microscopic viscosity substantially. The main result of local N-body simulations of planetary rings by Daisaka et al. (2001. Viscosity in a dense planetary ring with self-gravitating particles. Icarus 154, 296-312) is explained: in the presence of the gravitationally unstable density waves, the effective turbulent viscosity νeff is given as νeff=CG2Σ2/Ω3, where G, Σ, and Ω are the gravitational constant, the surface mass density of a ring, and the angular velocity, respectively, and the nondimensional correction factor C≈10. We argue that both Saturn's main rings and their irregular of the order of 100 m or even less fine-scale structure (being recurrently created and destroyed on the time scale of an order of Keplerian period ) are not likely much younger than the solar system. 相似文献
11.
Keiji Ohtsuki 《Icarus》2006,183(2):384-395
We examine rotation rates of gravitating particles in low optical depth rings, on the basis of the evolution equation of particle rotational energy derived by Ohtsuki [Ohtsuki, K., 2006. Rotation rate and velocity dispersion of planetary ring particles with size distribution. I. Formulation and analytic calculation. Icarus 183, 373-383]. We obtain the rates of evolution of particle rotation rate and velocity dispersion, using three-body orbital integration that takes into account distribution of random velocities and rotation rates. The obtained stirring and friction rates are used to calculate the evolution of velocity dispersion and rotation rate for particles in one- and two-size component rings as well as those with a narrow size distribution, and agreement with N-body simulation is confirmed. Then, we perform calculations to examine equilibrium rotation rates and velocity dispersion of gravitating ring particles with a broad size distribution, from 1 cm up to 10 m. We find that small particles spin rapidly with 〈ω2〉1/2/Ω?102-103, where ω and Ω are the particle rotation rate and its orbital angular frequency, respectively, while the largest particles spin slowly, with 〈ω2〉1/2/Ω?1. The vertical scale height of rapidly rotating small particles is much larger than that of slowly rotating large particles. Thus, rotational states of ring particles have vertical heterogeneity, which should be taken into account in modeling thermal infrared emission from Saturn's rings. 相似文献
12.
In this paper we present a new semianalytical model of oligarchic growth of planets considering a distribution of planetesimal sizes, fragmentation of planetesimals in mutual collisions, sublimation of ices through the snow line, random velocities out of equilibrium and merging of planetary embryos. We show that the presence of several planetary embryos growing simultaneously at different locations in the protoplanetary disk affects the whole accretion history, specially for the innermost planets. The results presented here clearly indicate the relevance of considering a distribution of planetesimal sizes. Fragmentation occurring during planetesimal-planetesimal collisions represent only a marginal effect in shaping the surface density of solid material in the protoplanetary disc. 相似文献
13.
Zoë M. LeinhardtDerek C. Richardson 《Icarus》2002,159(2):306-313
We present results from direct N-body simulations of collisions between gravitational aggregates of varying size as part of a study to parameterize planetesimal growth in the Solar System. We find that as the ratio of projectile to target mass departs from unity, the impact angle has less effect on the outcome. At the same time, the probability of planetesimal growth increases. Conversely, for a fixed impact energy, collisions between impactors with mass ratio near unity are more dispersive than those with impactor mass ratio far from unity. We derive an expression for the accretion probability as a function of mass ratio. For an average mass ratio of 1:5, we find an accretion probability of ∼60% over all impact parameters. We also compute the critical specific dispersal energy Q*D as a function of projectile size. Extrapolating to a projectile size of 1 m with a 1-km target, we find Q*D=103−104 J kg−1, in agreement with several other collision models that use fundamentally different techniques. Our model assumes that the components of each gravitational aggregate are identical and indestructible over the range of sampled impact speeds. In future work we hope to incorporate a simple fracture model to extend the range of applicable speeds and we plan to implement our results in a large-scale planetesimal evolution code. 相似文献
14.
Alexander V. Krivov 《中国天文和天体物理学报》2010,10(5)
Debris disks are optically thin, almost gas-free dusty disks observed arounda significant fraction of main-sequence stars older than about 10 Myr. Since the circumstellar dust is short-lived, the very existence of these disks is considered as evi-dence that dust-producing planetesimals are still present in mature systems, in whichplanets have formed – or failed to form – a long time ago. It is inferred that theseplanetesimals orbit their host stars at asteroid to Kuiper-belt distances and continuallysupply ... 相似文献
15.
Z. M. Leinhardt D. C. Richardson G. Lufkin J. Haseltine 《Monthly notices of the Royal Astronomical Society》2009,396(2):718-728
In this paper, we extend our numerical method for simulating terrestrial planet formation to include dynamical friction from the unresolved debris component. In the previous work, we implemented a rubble pile planetesimal collision model into direct N -body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10 per cent or more of the total disc mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead, they grow concurrently, similar to oligarchic growth. The dynamical friction from the unresolved debris damps the eccentricities of the planetesimals, reducing the mean impact speeds and causing all collisions to result in merging with no mass loss. As a result, there is no debris production. The mass in debris slowly decreases with time. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions (for which the initial mass in unresolved debris is at least equal to the mass in resolved planetesimals). 相似文献
16.
John Chambers 《Icarus》2006,180(2):496-513
A new semi-analytic model for the oligarchic growth phase of planetary accretion is developed. The model explicitly calculates damping and excitation of planetesimal eccentricities e and inclinations i due to gas drag and perturbations from embryos. The effects of planetesimal fragmentation, enhanced embryo capture cross sections due to atmospheres, inward planetesimal drift, and embryo-embryo collisions are also incorporated. In the early stages of oligarchic growth, embryos grow rapidly as e and i fall below their equilibrium values. The formation of planetesimal collision fragments also speeds up embryo growth as fragments have low-e, low-i orbits, thereby optimizing gravitational focussing. At later times, the presence of thick atmospheres captured from the nebula aids embryo growth by increasing their capture cross sections. Planetesimal drift due to gas drag can lead to substantial inward transport of solid material. However, inward drift is greatly reduced when embryo atmospheres are present, as the drift timescale is no longer short compared to the accretion timescale. Embryo-embryo collisions increase embryo growth rates by 50% compared to the case where growth is solely due to accretion of planetesimals. Formation of 0.1-Earth-mass protoplanets at 1 AU and 10-Earth-mass cores at 5 AU requires roughly 0.1 and 1 million years respectively, in a nebula where the local solid surface density is 7 g cm−2 at each of these locations. 相似文献
17.
We obtain the viscous stirring and dynamical friction rates of planetesimals with a Rayleigh distribution of eccentricities and inclinations, using three-body orbital integration and the procedure described by Ohtsuki (1999, Icarus137, 152), who evaluated these rates for ring particles. We find that these rates based on orbital integrations agree quite well with the analytic results of Stewart and Ida (2000, Icarus 143, 28) in high-velocity cases. In low-velocity cases where Kepler shear dominates the relative velocity, however, the three-body calculations show significant deviation from the formulas of Stewart and Ida, who did not investigate the rates for low velocities in detail but just presented a simple interpolation formula between their high-velocity formula and the numerical results for circular orbits. We calculate evolution of root mean square eccentricities and inclinations using the above stirring rates based on orbital integrations, and find excellent agreement with N-body simulations for both one- and two-component systems, even in the low-velocity cases. We derive semi-analytic formulas for the stirring and dynamical friction rates based on our numerical results, and confirm that they reproduce the results of N-body simulations with sufficient accuracy. Using these formulas, we calculate equilibrium velocities of planetesimals with given size distributions. At a stage before the onset of runaway growth of large bodies, the velocity distribution calculated by our new formulas are found to agree quite well with those obtained by using the formulas of Stewart and Ida or Wetherill and Stewart (1993, Icarus106, 190). However, at later stages, we find that the inclinations of small collisional fragments calculated by our new formulas can be much smaller than those calculated by the previously obtained formulas, so that they are more easily accreted by larger bodies in our case. The results essentially support the previous results such as runaway growth of protoplanets, but they could enhance their growth rate by 10-30% after early runaway growth, where those fragments with low random velocities can significantly contribute to rapid growth of runaway bodies. 相似文献
18.
In this chapter, we will give a brief overview on our current theoretical understanding how planets form from the solid material in circumstellar disks in the core accretion-gas capture model. This chapter will not be as concise and complete as a review on this matter, yet will serve as an introductory text to generate interest in the subject. Students are referred to comprehensive text books and some important reviews.This chapter will discuss “dusty storms”, e.g. the dust transport in turbulent protoplanetary disks, followed by the latest model of planetesimal formation, e.g. gravoturbulent planetesimal formation, which deals with particle concentration in turbulence and N-body simulations thereof. We also briefly describe the core accretion-gas capture process and talk about nascent planets, e.g. the observability of planet–disk interaction concluding with the migration of young planets and the final arrangement of planetary systems. 相似文献
19.
Raine Karjalainen 《Icarus》2007,189(2):523-537
Ring particle aggregates are formed in the outer parts of Saturn's main rings. We study how collisions between aggregates can lead to destruction or coalescence of these aggregates, with local N-body simulations taking into account the dissipative impacts and gravitational forces between particles. Impacts of aggregates with different mass ratios are studied, as well as aggregates that consist of particles with different physical properties. We find that the outcome of the collision is very sensitive to the shape of the aggregate, in the sense that more elongated aggregates are more prone to be destroyed. We were interested in testing the accretion criterion Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171] used in their F ring simulations, according to which accretion requires that the masses of the colliding bodies differ at least by a factor of 100. We confirm that such a critical mass ratio exists. In particular, simulations indicate that the exact critical mass ratio depends on the internal density and elasticity of particles, and the distance from the planet. The zone of transition, defined by the distance where individual particles or small aggregates first start to stick on the larger aggregate, and by the distance where two similar sized aggregates on the average eventually coalesce is only about 5000 km wide, if fixed particle properties are used. The rotational state of the aggregates that form via aggregate collision rapidly reaches synchronous rotation, similarly to the aggregates that form via gradual growth. 相似文献
20.
When preplanetary bodies reach proportions of ∼1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. The numerical model is novel as it attempts to preserve the individual particle nature of the bodies (like N-body codes); yet it is statistical in nature since it must incorporate the very large number of planetesimals. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R0. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We associate the initial runaway growth phase with increasing GF-factors for the largest body. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high-mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two-component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. We identify the peak of GF with the transition between the runaway growth and oligarchy accretion stages. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius, especially at the outer disk. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion, which could lead to a very fast growth. 相似文献