共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
The rotation rate distribution of small Main Belt asteroids is dominated by YORP and collisions. These mechanism act differently depending on the size of the bodies and give rise to non-linear effects when they both operate. Using a Monte Carlo method we model the formation of a steady state population of small asteroids under the influence of both mechanisms and the rotation rate distribution is compared to the observed one as derived from Pravec et al. (Pravec, P. et al. [2008]. Icarus 197, 497-504). A better match to observations is obtained with respect to the case in which only YORP is considered. In particular, an excess of slow rotators is produced in the model with both collisions and YORP because bodies driven to slow rotation by YORP have a random walk-like evolution of the spin induced by repeated collisions with small projectiles. This is a dynamical evolution different from tumbling and it lasts until a large impact takes the body to a faster rotation rate. According to our model, the rotational fission of small asteroids is a very frequent event and might explain objects like P/2010 A2 and its associated tail of millimeter-sized dust particles. The mass loss during fission of small asteroids might significantly influence the overall collisional evolution of the belt. Fission can in fact be considered as an additional erosion mechanism, besides cratering and fragmentation, acting only at small diameters. 相似文献
3.
We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune. 相似文献
4.
The effect of a neutron-proton vortex system on the rotation dynamics of neutron stars is examined. The dynamics of the motion of a two component superfluid system in the core of a neutron star yields an equation for the evolution of the pulsar's rotation period. The spin down of the star owing to energy release at the core boundary, which is associated with a contraction of the length of the neutron vortex as it moves radially and magnetic energy of the vortical cluster is released, is taken into account. Evolutionary curves are constructed for pulsars with different magnetic fields and stellar radii. For certain values of the coefficient of friction between the superfluid and normal components in the core of the neutron star, at the end of its evolution a radio pulsar may become an anomalous x-ray pulsar or a source of soft gamma radiation with a period on the order of 10 seconds. 相似文献
5.
H. Salo 《Earth, Moon, and Planets》1987,38(2):149-181
Hämeen-Anttila's (1984) analytical treatment of self-gravitating collisional particle disks is extended to include the particle spin. The equations derived for the coupled evolution of random velocities and spins indicate that friction and surface irregularity typically reduce the local velocity dispersion. Friction, and especially irregularity, also transfer significant amounts of random kinetic energy, E
kin, to rotational energy, E
rot. The equilibrium ratio E
rot/Ekin = 2/(14 – 5) if the particles are spherical, and 2(1 + )/7 if they are irregular but frictionless, and being the coefficients of restitution and friction. These results are not only exact for identical, non-gravitating mass points, but are rather accurate even if finite size, self-gravitational forces, or size distribution are included. Applications to the dynamics of Saturn's rings suggest that the inclusion of rotation is able to reduce the geometrical thickness of the layer of centimeter-sized particles to about one half, at most. Large particles are less affected.On leave from Dept. of Astronomy, University of Oulu, 90570 Oulu, Finland. 相似文献
6.
We present a numerical algorithm designed to study the evolution of a distribution of solid particles orbiting around the Sun that could be applied to similar systems in which f c ≤ Ω/2π, where f c is the frequency of collisions, and Ω is the orbital angular speed. A number of sample particles are used to represent the whole system. Binary collisions are treated in a novel way using a Monte Carlo method that works as follows. Orbits are locally sampled to compute the velocity dispersion. Then the velocity vectors of the sample particles are modified according to random collisions with virtual particles which have velocities taken from a normal distribution computed using the previously found local velocity dispersion. The energy and momentum taken up by the virtual particles are redistributed among the neighbours of the sample particle undergoing the collision, so that conservation laws are satisfied. Simulations using this model give an estimation of the final distribution of inclinations and the associated evolutionary time-scale. 相似文献
7.
We formulate a complete system of equations of two-phase multicomponent mechanics including the relative motion of the phases, coagulation processes, phase transitions, chemical reactions, and radiation in terms of the problem of reconstructing the evolution of the protoplanetary gas-dust cloud that surrounded the proto-Sun at an early stage of its existence. These equations are intended for schematized formulations and numerical solutions of special model problems on mutually consistent modeling of the structure, dynamics, thermal regime, and chemical composition of the circumsolar disk at various stages of its evolution, in particular, the developed turbulent motions of a coagulating gas suspension that lead to the formation of a dust subdisk, its gravitational instability, and the subsequent formation and growth of planetesimals. To phenomenologically describe the turbulent flows of disk material, we perform a Favre probability-theoretical averaging of the stochastic equations of heterogeneous mechanics and derive defining relations for the turbulent flows of interphase diffusion and heat as well as for the “relative” and Reynolds stress tensors needed to close the equations of mean motion. Particular attention is given to studying the influence of the inertial effects of dust particles on the properties of turbulence in the disk, in particular, on the additional generation of turbulent energy by large particles near the equatorial plane of the proto-Sun. We develop a semiempirical method of modeling the coefficient of turbulent viscosity in a two-phase disk medium by taking into account the inverse effects of the transfer of a dispersed phase (or heat) on the growth of turbulence to model the vertically nonuniform thermohydrodynamic structure of the subdisk and its atmosphere. We analyze the possible “regime of limiting saturation” of the subdisk atmosphere by fine dust particles that is responsible for the intensification of various coagulation mechanisms in a turbulized medium. For steady motion when solid particles settle to the midplane of the disk under gravity, we analyze the parametric method of moments for solving the Smoluchowski integro-differential coagulation equation for the particle size distribution function. This method is based on the fact that the sought-for distribution function a priori belongs to a certain parametric class of distributions. 相似文献
8.
K. A. Hämeen-Anttila 《Astrophysics and Space Science》1976,43(1):145-174
A dispersion in particle sizes is incorporated into the theory of partially elastic collisions on Keplerian orbits. Very different particles can never co-exist for long periods, if their elasticity parameter depends of the kinetic energy in collisions. If a system initially consists of a mixture of sizes, two possibilities for its evolution exist. If the smaller particles initially dominate, then they fall rapidly to the central body without affecting the remaining system. If the larger particles dominate, then a ringshaped system is formed. This is one possible mechanism for the origin of Saturn's rings. A dust content is also able to create a gap like the Cassini division. 相似文献
9.
Our goal is to study the regime of disk accretion in which almost all of the angular momentum and energy is carried away by the wind outflowing from the disk in numerical experiments. For this type of accretion the kinetic energy flux in the outflowing wind can exceed considerably the bolometric luminosity of the accretion disk, what is observed in the plasma flow from galactic nuclei in a number of cases. In this paper we consider the nonrelativistic case of an outflow from a cold Keplerian disk. All of the conclusions derived previously for such a system in the self-similar approximation are shown to be correct. The numerical results agree well with the analytical predictions. The inclination angle of the magnetic field lines in the disk is less than 60°, which ensures a free wind outflow from the disk, while the energy flux per wind particle is greater than the particle rotation energy in its Keplerian orbit by several orders of magnitude, provided that the ratio r A/r ? 1, where r A is the Alfvénic radius and r is the radius of the Keplerian orbit. In this case, the particle kinetic energy reaches half the maximum possible energy in the simulation region. The magnetic field collimates the outflowing wind near the rotation axis and decollimates appreciably the wind outflowing from the outer disk periphery. 相似文献
10.
Asteroids have a wide range of rotation states. While the majority spin a few times to several times each day in principal axis rotation, a small number spin so slowly that they have somehow managed to enter into a tumbling rotation state. Here we investigate whether the Yarkovsky-Radzievskii-O'Keefe-Paddack (YORP) thermal radiation effect could have produced these unusual spin states. To do this, we developed a Lie-Poisson integrator of the orbital and rotational motion of a model asteroid. Solar torques, YORP, and internal energy dissipation were included in our model. Using this code, we found that YORP can no longer drive the spin rates of bodies toward values infinitely close to zero. Instead, bodies losing too much rotation angular momentum fall into chaotic tumbling rotation states where the spin axis wanders randomly for some interval of time. Eventually, our model asteroids reach rotation states that approach regular motion of the spin axis in the body frame. An analytical model designed to describe this behavior does a good job of predicting how and when the onset of tumbling motion should take place. The question of whether a given asteroid will fall into a tumbling rotation state depends on the efficiency of its internal energy dissipation and on the precise way YORP modifies the spin rates of small bodies. 相似文献
11.
A thin gaseous disk with a nearly Keplerian rotation profile and free boundaries in the external gravitational field of a
point gravitating object does not generate any growing perturbation eigenmodes. In spite of this, a significant transient
growth of linear perturbations measured by the evolution of their total acoustic energy is possible in such a disk. This is
shown within the framework of the simplest model of an inviscid polytropic thin disk with a finite radial extent in which
small adiabatic perturbations that are a linear combination of neutral eigenmodes with a corotation radius beyond the outer
flow boundary are considered. 相似文献
12.
Christian Fendt 《Astrophysics and Space Science》2003,287(1-4):59-64
The evolution of a stellar, initially dipole type magnetosphere interacting with an accretion disk is investigated using numerical ideal MHD simulations. The simulations follow several 1000 Keplerian periods of the inner disk (for animated movies see http://www.aip.de～cfendt).Our model prescribes a Keplerian disk around a rotating star as a fixed boundary condition. The initial magnetic field distribution remains frozen into the star and the disk. The mass flow rate into the corona is fixed for both components. The initial dipole type magnetic field develops into a spherically radial outflow pattern with two main components – a disk wind and a stellar wind – both evolving into a quasi-stationary final state. A neutral field line divides both components, along which small plasmoids are ejected in irregular time intervals. The half opening angle of the stellar wind cone varies from 30° to55° depending on the ratio of the mass flow rates of disk wind and stellar wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. An axial jet forms during the first decades of rotations. However, this feature does not survive on the very long time scale and a pressure driven low velocity flow along the axis evolves. Within a cone of 15° along the axis the formation of knots may be observed if the stellar wind is weak. With the chosen mass flow rates and field strength we see almost no indication for a flow self-collimation. This is due to the weak net poloidal electric current in the magnetosphere which is in difference to typical jet models. 相似文献
13.
A. I. Neishtadt D. J. Scheeres V. V. Sidorenko P. J. Stooke A. A. Vasiliev 《Celestial Mechanics and Dynamical Astronomy》2003,86(3):249-275
Reactive torques, due to anisotropic sublimation on a comet nucleus surface, produce slow variations of its rotation. In this paper the secular effects of this sublimation are studied. The general rotational equations of motion are averaged over unperturbed fast rotation around the mass center (Euler-Poinsot motion) and over the orbital comet motion. We discuss the parameters that define typical properties of the rotational evolution and discover different classifications of the rotational evolution. As an example we discuss some possible scenarios of rotational evolution for the nuclei of the comets Halley and Borrelly. 相似文献
14.
J. Lukkari 《Earth, Moon, and Planets》1989,44(3):207-218
Quasi-equilibrium solutions for the pre-planetary disk are studied in terms of Hämeen-Anttila's theory (1984) of collisional, self-gravitating systems. The distribution of particle sizes is assumed to follow simple power-law distributions, with a power index in the range of 1.5–5.0. The treatment includes mutual impacts with a velocity dependent coefficient of restitution, as well as gravitational encounters with dynamical friction. The mean gravitational field of the disk is also taken into account. The results indicate that the energy(equi)-partition depends mainly on the index of size distribution, but is also affected by the optical thickness of the system, as well as on the vertical thickness as compared to the particle size. The vertical component of the gravitational field is found to be important, especially when the mass of the system is concentrated on the large particles. 相似文献
15.
In this paper we propose firstly a full parametrization of an impact on a target body considered as ellipsoidal, including several geometrical parameters which are generally not included. Then we construct a more detailed and complete theoretical model of the rotational changes of the target body arising from a single impact, by taking into account the various parameters above. Secondly from these theoretical studies we carry out simulations of impacts and then we evaluate the influences of the various parameters on the rotational evolution of a specific target, in particular the angular speed of rotation and the direction of the axis of rotation. For that we consider two cases: in the first one, which we call accretion, the projectile is simply stuck to the target without a significant amount of ejected mass. In the second case, which we call craterization, the target body is eroded with formation of ejecta and a crater. The physical properties of the target are close to those of the Asteroid 21 Lutetia which Rosetta mission would fly in July 2010. We obtain quite different results according to the considered mode of impact (accretion or craterization): in the case of an impact with accretion the results are intuitively foreseeable whereas those corresponding to an impact with craterization are more difficult to interpret. Our work can be applied to obtain information on the rotational effects of an impact on a given target body with well constrained physical characteristics, in particular within the framework of the Don Quijote mission project. 相似文献
16.
This paper considers, in the context of modeling the evolution of a protoplanetary cloud, the hydrodynamic aspects of the theory of concurrent processes of mass transfer and coagulation in a two-phase medium in the presence of shear turbulence in a differentially rotating gas–dust disk and of polydisperse solid particles suspended in a carrying flow of solid particles. The defining relations are derived for diffuse fluxes of particles of different sizes in the equations of turbulent diffusion in the gravitational field, which describe the convective transfer, turbulent mixing, and sedimentation of disperse dust grains onto the central plane of the disk, as well as their coagulation growth. A semiempirical method is developed for calculating the coefficients of turbulent viscosity and turbulent diffusion for particles of different kinds. This method takes into account the inverse effects of dust transfer on the turbulence evolution in the disk and the inertial differences between disperse solid particles. To solve rigorously the problem of the mutual influence of the turbulent mixing and coagulation kinetics in forming the gas–dust subdisk, the possible mechanisms of gravitational, turbulent, and electric coagulation in a protoplanetary disk are explored and the parametric method of moments for solving the Smoluchowski integro-differential coagulation equation for the particles' size distribution function is considered. This method takes into account the fact that this distribution belongs to a definite parametric class of distributions. 相似文献
17.
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. 相似文献
18.
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. 相似文献
19.
Hannes Alfvén 《Astrophysics and Space Science》1970,6(1):161-174
Ifviscosity is taken into account, Keplerian motion of a large number of grains in a gravitational field has a tendency to lead to the formation ofjet streams.In order to treat problems of this kind it is advantageous to present celestial mechanics by a simple perturbation approach which is developed in Sections 2–7.Inelastic collisions between a number of grains will tend to make their orbits similar. This leads to the formation of jet streams. Their properties are treated in Sections 8–10.Finally, in Section 11, we discuss the possible application of the jet stream theory to meteor streams, to asteroidal jet streams, and to the cosmogonic accretion process. 相似文献
20.
Song-nian Huang 《Chinese Astronomy and Astrophysics》1989,13(4):442-453
We discuss the different dynamical effects of spirals with live and rigid haloes in collisions. Our N-body simulation has 52000 particles. The target is a self-consistent three-dimensional system consisting of a disk and a live halo, the intruder galaxy has a fixed Plummer potential and moves in a straight line through the centre of the target in the plane of the disk. For comparison, the same numerical experiment is made with a rigid halo, all other conditions remaining the same. Differences in the dynamical response and the resulting morphology between the two cases are discussed. It is suggested that the torque on non-axial symmetric spirals exerted by triaxial haloes and the long term effect of dynamical friction of the halo should be further investigated. 相似文献