首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Rotational Properties of Cometary Nuclei   总被引：1，自引：0，他引：1  
We review several techniques used to retrieve rotational parameters from observations. The spin period of a dozen of comets retrieved with these techniques are summarized. We describe how the spin period of comet Hale-Bopp (C/1995 O1) has been calculated with a high accuracy (11.30–11.34 h). Although several authors converged to a spin axis orientation at (α,δ) = (275 ± 15°, -55 ± 5°), detailed studies indicate that the dust jets morphology in 1996–1997 may be incompatible with this orientation. Comet 19P/Borrelly has been recently observed by the Deep Space 1 spacecraft. At the same time, its spin axis orientation and period have been determined by several authors to be respectively (α,δ) = (225 ± 15°, -10 ± 10°)and 26h. These two comets are likely to be in (or close to) a principal axis spin state. We discuss new modeling of the spin state of comet 46P/Wirtanen, the target of the Rosetta mission. The model involves a three-dimensional shape and thermal model, from which the torque of the non gravitational force is calculated at each time step. The moments of inertia are computed for each irregular shape. The results from numerical integrations show that this comet can remain in a principal axis spin state during more than 10 orbits if the spin period does not get above～6 h. If the spin period increases, its nucleus gets rapidly into excited spin states. It shows that even small and very active short-period comets are not necessarily in non principal axis spin states. In the last section, the consequences of recent observations and modeling of the rotational parameters of comet nuclei are discussed, and unsolved problems are presented.  相似文献   

2.
The role of non-gravitational forces in the evolution of orbitalmotion of C/1995 O1 (Hale–Bopp) has been investigated. Inorbital calculations the observational material covering theperiod from April 1993 up to August 2001 was used. To model thenon-gravitational acceleration, observed and theoretical profilesof the H2O production rates were employed. A set of forcedprecession models of a rotating cometary nucleus consistent withthe observed spin axis orientation was fitted to positionalobservations. The non-gravitational models allowed us to constrainthe mass and radius of the comet. The orbitalevolution of Comet Hale–Bopp was investigated over ±400 k yusing two sets of randomly varied orbital elements wellrepresenting all positional observations in the pure gravitationalcase, as well as in the non-gravitational case. The calculationsshowed that the comet's motion is predictable only over an interval ofa few orbital periods. The statistical conclusions changesignificantly when non-gravitational effects are included in the analysis.  相似文献   

3.
By considering model comet nuclei with a wide range of sizes, prolate ellipsoidal shapes, spin axis orientations, and surface activity patterns, constraints have been placed on the nucleus properties of the primary Rosetta target, Comet 67P/Churyumov-Gerasimenko. This is done by requiring that the model bodies simultaneously reproduce the empirical nucleus rotational lightcurve, the water production rate as function of time, and non-gravitational changes (per apparition) of the orbital period (ΔP), longitude of perihelion (Δ?), and longitude of the ascending node (ΔΩ). Two different thermophysical models are used in order to calculate the water production rate and non-gravitational force vector due to nucleus outgassing of the model objects. By requiring that the nominal water production rate measurements are reproduced as well as possible, we find that the semi-major axis of the nucleus is close to 2.5 km, the nucleus axis ratio is approximately 1.4, while the spin axis argument is either 60°±15° or 240°±15°. The spin axis obliquity can only be preliminarily constrained, indicating retrograde rotation for the first argument value, and prograde rotation for the second suggested spin axis argument. A nucleus bulk density in the range 100-370 kg m−3 is found for the nominal ΔP, while an upper limit of 500 kg m−3 can be placed if the uncertainty in ΔP is considered. Both considered thermophysical models yield the same spin axis, size, shape, and density estimates. Alternatively, if calculated water production rates within an envelope around the measured data are considered, it is no longer possible to constrain the size, shape, and spin axis orientation of the nucleus, but an upper limit on the nucleus bulk density of 600 kg m−3 is suggested.  相似文献   

4.
The nucleus of Comet 81P/Wild 2 is modeled by assuming various smooth triaxial ellipsoidal or irregular body shapes, having different rotational periods, spin axis orientations, and thermophysical properties. For these model nuclei, a large number of surface activity patterns (e.g., maps of active and inactive areas) are studied, and in each case the resulting water production rate and non-gravitational force vector versus time are calculated. By requiring that the model nuclei simultaneously reproduce certain properties of the empirical water production curve and non-gravitational changes of the orbit (focusing on changes of the orbital period and in the longitude of perihelion), constraints are placed on several properties of the nucleus. The simulations suggest that the mass of Comet 81P/Wild 2 is , resulting in a rather low bulk density, (depending on the assumed nucleus volume), and that the nucleus rotation is prograde rather than retrograde. The active area fraction is difficult to constrain, but at most 60% of the nucleus is likely to have near-surface ice.  相似文献   

5.
The nucleus of Comet 2P/Encke was detected with the Arecibo radar during the close approach of November, 2003, making this the first comet to yield radar detections at two different apparitions. Although the measured radar cross section of 1.0 km2 was close to that obtained in 1980, the Doppler bandwidth was nearly four times larger. Most of this bandwidth difference can simply be attributed to a different observing aspect relative to the spin axis proposed by Sekanina [1988, Astron. J. 95, 911] and Festou and Barale [2000, Astron. J. 119, 3119]. Comparison of the 2003 Doppler bandwidth with infrared-based size estimates supports an 11-h dominant rotation period and excludes slower 15- and 22-h periods that have also been suggested. If one assumes a short-axis-mode rotation with an 11-h period, then the Doppler bandwidth indicates that the nucleus is an oblong object with a long-axis dimension of 9 km. The estimated radar albedo of 0.05 is similar to that measured for C/IRAS-Araki-Alcock, providing further evidence that comet nuclei have relatively low surface densities of ∼0.5-1.0 g cm−3. No broadband echo component was detected from large coma grains despite predictions, based on optical/infrared models, that such a component might be detectable.  相似文献   

6.
S.J. Peale 《Icarus》2006,181(2):338-347
In determining Mercury's core structure from its rotational properties, the value of the normalized moment of inertia, C/MR2, from the location of Cassini 1 is crucial. If Mercury's spin axis occupies Cassini state 1, its position defines the location of the state, where the axis is fixed in the frame precessing with the orbit. Although tidal and core-mantle dissipation drive the spin to the Cassini state with a time scale O(105) years, the spin might still be displaced from the Cassini state if the variations in the orbital elements induced by planetary perturbations, which change the position of the Cassini state, cause the spin to lag behind as it attempts to follow the state. After being brought to the state by dissipative processes, the spin axis is expected to follow the Cassini state for orbit variations with time scales long compared to the 1000 year precession period of the spin about the Cassini state because the solid angle swept out by the spin axis as it precesses is an adiabatic invariant. Short period variations in the orbital elements of small amplitude should cause displacements that are commensurate with the amplitudes of the short period terms. The exception would be if there are forcing terms in the perturbations that are nearly resonant with the 1000 year precession period. The precision of the radar and eventual spacecraft measurements of the position of Mercury's spin axis warrants a check on the likely proximity of the spin axis to the Cassini state. How confident should we be that the spin axis position defines the Cassini state sufficiently well for a precise determination of C/MR2? By following simultaneously the spin position and the Cassini state position during long time scale orbital variations over past 3 million years [Quinn, T.R., Tremaine, S., Duncan, M., 1991. Astron. J. 101, 2287-2305] and short time scale variations for 20,000 years [JPL Ephemeris DE 408; Standish, E.M., private communication, 2005], we show that the spin axis will remain within one arcsec of the Cassini state after it is brought there by dissipative torques. In this process the spin is located in the orbit frame of reference, which in turn is referenced to the inertial ecliptic plane of J2000. There are no perturbations with periods resonant with the precession period that could cause large separations. We thus expect Mercury's spin to occupy Cassini state 1 well within the uncertainties for both radar and spacecraft measurements, with correspondingly tight constraints on C/MR2 and the extent of Mercury's molten core. Two unlikely caveats for this conclusion are: (1) an excitation of a free spin precession by an unknown mechanism or (2) a displacement by a dissipative core mantle interaction that exceeds the measurement uncertainties.  相似文献   

7.
Two classes of high-energy sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars are believed to contain slowly spinning “magnetars,” i.e. neutron stars the emission of which derives from the release of energy from their extremely strong magnetic fields (>1015 G). The enormous energy liberated in the 2004 December 27 giant flare from SGR 1806-20 (～5×1046 erg), together with the likely recurrence time of such events, points to an internal magnetic field strength of ≥1016 G. Such strong fields are expected to be generated by a coherent α?Ω dynamo in the early seconds after the Neutron Star (NS) formation, if its spin period is of a few milliseconds at most. A substantial deformation of the NS is caused by such fields and, provided the deformation axis is offset from the spin axis, a newborn millisecond-spinning magnetar would thus radiate for a few days a strong gravitational wave signal the frequency of which (～0.5–2 kHz range) decreases in time. This signal could be detected with Advanced LIGO-class detectors up to the distance of the Virgo cluster, where ≥1 yr?1 magnetars are expected to form. Recent X-ray observations revealed that SNRs around magnetar candidates do not appear to have received a larger energy input than in standard SNRs (see Vink and Kuiper, Mon. Not. Roy. Astron. Soc. 319, L14 (2006)). This is at variance with what would be expected if the spin energy of the young, millisecond NS were radiated away as electromagnetic radiation and/or relativistic particle winds. In fact, such energy would be transferred quickly and efficiently to the expanding gas shell. This may thus suggest that magnetars did not form with the expected very fast initial spin. We show here that these findings can be reconciled with the idea of magnetars being formed with fast spins, if most of their initial spin energy is radiated through GWs. In particular, we find that this occurs for essentially the same parameter range that would make such objects detectable by Advanced LIGO-class detectors up to the Virgo Cluster. If our argument holds for at least a fraction of newly formed magnetars, then these objects constitute a promising new class of gravitational wave emitters.  相似文献   

8.
Saturn’s moon, Hyperion, is subject to strongly-varying solid body torques from its primary and lacks a stable spin state resonant with its orbital frequency. In fact, its rotation is chaotic, with a Lyapunov timescale on the order of 100 days. In 2005, Cassini made three close passes of Hyperion at intervals of 40 and 67 days, when the moon was imaged extensively and the spin state could be measured. Curiously, the spin axis was observed at the same location within the body, within errors, during all three fly-bys—~ 30° from the long axis of the moon and rotating between 4.2 and 4.5 times faster than the synchronous rate. Our dynamical modeling predicts that the rotation axis should be precessing within the body, with a period of ~ 16 days. If the spin axis retains its orientation during all three fly-bys, then this puts a strong constraint on the in-body precessional period, and thus the moments of inertia. However, the location of the principal axes in our model are derived from the shape model of Hyperion, assuming a uniform composition. This may not be a valid assumption, as Hyperion has significant void space, as shown by its density of 544± 50  kg m−3 (Thomas et al. in Nature 448:50, 2007). This paper will examine both a rotation model with principal axes fixed by the shape model, and one with offsets from the shape model. We favor the latter interpretation, which produces a best-fit with principal axes offset of ~ 30° from the shape model, placing the A axis at the spin axis in 2005, but returns a lower reduced χ 2 than the best-fit fixed-axes model.  相似文献   

9.
Possible dynamical evolution of the rotation of Venus since formation   总被引：1，自引：0，他引：1  
The past evolution of the rotation of Venus has been studied by a numerical integration method using the hypothesis that only solar tidal torques and core-mantle coupling have been active since formation. It is found quite conceivable that Venus had originally a rotation similar to the other planets and has evolved in 4.5×109 years from a rapid and direct rotation (12-hour spin period and nearly zero obliquity) to the present slow retrograde one.While the solid tidal torque may be quite efficient in despinning the planet, a thermally driven atmospheric tidal torque has the capability to drive the obliquity from 0° towards 180° and to stabilize the spin axis in the latter position. The effect of a liquid core is discussed and it is shown that core-mantle friction hastens the latter part of the evolution and makes even stronger the state of equilibrium at 180°. The model assumes a nearly stable balance between solid and atmospheric tides at the current rotation rate interpreting the present 243 day spin period as being very close to the limiting value.A large family of solutions allowing for the evolution, in a few billions years, of a rapid prograde rotation to the present state have been found. Noticeably different histories of evolution are observed when the initial conditions and the values of the physical parameters are slightly modified, but generally the principal trend is maintained.The proposed evolutionary explanation of the current rotation of Venus has led us to place constraints on the solid bodyQ and on the magnitude of the atmospheric tidal torque. While the constraints seem rather severe in the absence of core-mantle friction (aQ15 at the annual frequency is required, and a dominant diurnal thermal response in the atmosphere is needed), for a large range of values of the core's viscosity, the liquid core effect allows us to relax somewhat these constraints: a solid bodyQ of the order 40 can then be allowed. ThisQ value implies that a semi-diurnal ground pressure oscillation of 2 mb is needed in the atmosphere in order for a stable balance to occur between the solid and atmospheric tides at the current rotation rate. No model of atmospheric tides on Venus has been attempted in this study, however the value of 2 mb agrees well with that predicted by the model given in Dobrovolskis (1978).  相似文献   

10.
In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.  相似文献   

11.
This paper reviews our current knowledge of Saturn's rings’ physical properties as derived from thermal infrared observations. Ring particle composition, surface structure and spin as well as the vertical structure of the main rings can be determined. These properties are the key to understand the origin and evolution of Saturn's rings. Ring composition is mainly constrained by observations in the near-infrared but the signature of some probable contaminants present in water ice may also be found at mid-infrared wavelengths. The absence of the silicate signature limits nowadays their mass fraction to 10−7±1. Recent measurements on the thermal inertia of the ring particle surface show it is very low, of the order of 5±2 Jm−2 K−1 s−1/2. New models and observations of the complete crossing of the planetary shadow are needed to attribute this low value either to compact regoliths covered by cracks due to collisions and thermal stresses or to large fluffy and irregular surfaces. Studies of the energy balance of ring particles show a preference for slowly spinning particles in the main rings. Supplementary observations at different phase angles, showing the temperature contrast between night and day sides of particles, and new models including finite spin and thermal inertia, are needed to constrain the actual spin distribution of ring particles. These results can then be compared to numerical simulations of ring dynamics. Many thermal models have been proposed to reproduce observations of the main rings, including alternative mono- or many-particles-thick layers or vertical heterogeneity, with no definitive answer. Observations on the lit and dark faces of rings as a function of longitude, at many incidence and emission angles, would provide prime information on the vertical thermal gradient due to interparticle shadowing from which constraints on the local vertical structure and dynamics can be produced. Future missions such as Cassini will provide new information to further constrain the ring thermal models.  相似文献   

12.
Extensive observations of comet 260P/McNaught were carried out between August 2012 and January 2013. The images obtained were used to analyze the comet’s inner coma morphology at resolutions ranging from 250 to about 1000 km/pixel. A deep investigation of the dust features in the inner coma allowed us to identify only a single main active source on the comet’s nucleus, at an estimated latitude of ?50°±15°. A thorough analysis of the appearance and of the motion of the morphological structures, supported by graphic simulations of the geometrical conditions of the observations, allowed us to determine a pole orientation located within a circular spot of a 15°-radius centered at RA=60°, Dec=0°. The rotation of the nucleus seems to occur on a single axis and is not chaotic, furthermore no precession effects could be estimated from our measurements. The comet’s spin axis never reached the plane of the sky from October 2012 to January 2013; during this period it did not change its direction significantly (less than 30°), thus giving us the opportunity to observe mainly structures such as bow-shaped jets departing from the single active source located on the comet’s nucleus. Only during the months of August 2012 and January 2013 the polar axis was directed towards the Earth at an angle of about 45° from the plane of the sky; this made it possible to observe the development of faint structures like fragments of shells or spirals. A possible rotation period of 0.340±0.01 days was estimated by means of differential photometric analysis.  相似文献   

13.
There are many aspects of observational evidence that cometary nuclei have irregular or nonspherical shape. The triaxial figure of the Halley's Comet nucleus is a well known fact. Therefore, the nucleus shape plays a significant role in consideration of the formation and evolution of comets and several attempts have been made to explain their nonsphericity. These studies were mainly based on the random-walk schemes for the aggregation processes. Although some results indeed lead to irregularities and deviation from sphericity, the spherical or irregular shape seem to be prevailing results. On the other hand the triaxial figure can be formed by the tidal and rotational forces. Thus, the assumption that the shape of the cometary nucleus due to some of these effects is in principle acceptable. In here assumed scenario already evolved cometary nucleus is situated as a satellite in the gravitation field of a planetary-like body. Since the rigidity of the nucleus is low, it may be easily transferred in the state of a synchronous satellite and in its shape could be imprinted the dynamical effects from this epoch. Here presented results indicate, that such a possibility should be seriously considered. The theory of this process is applied to the nucleus of comet Halley. It is shown, that the nucleus might be synchronously orbiting around a planetary-like hypothetical body with a period of 0.7 days. The minimal bulk tensile strength of the cometary material of about 102 N m–2 is estimated.  相似文献   

14.
Coupling of thermal evolution and despinning of early Iapetus   总被引：2，自引：0，他引：2  
The Cassini mission revealed two spectacular characteristics of Iapetus: (1) a geologically old and high equatorial ridge, which is unique in the Solar System and (2) a large flattening of 35 km consistent with the equilibrium figure for a hydrostatic body rotating with a period of 16 h, whereas the current spin period is 79.33 days. This study describes three-dimensional simulations of solid-state convection within an undifferentiated Iapetus. It investigates the implications for the evolution of the interior thermal structure and its spin rate and global shape using radially layered viscoelastic models. The role of the concentration in the short-lived radiogenic element [26Al], just after accretion is completed, is specifically addressed. The first result is to show that whatever the [26Al] value, convection occurs. As suggested by Castillo-Rogez et al. [Castillo-Rogez, J., Matson, D., Sotin, C., Johnson, T., Lunine, J., Thomas, P. [2007] Icarus, 190, 179-202], convection reduces the warming of the interior compared to the conductive evolution and therefore limits the conditions for despinning. In our calculations, two conceptual linear viscoelastic models are used. When considering a Maxwell rheology, the interior temperature (viscosity) never reaches a value high (low) enough to induce despinning. In order to promote dissipation at low temperature, a Burgers rheology, which includes an additional dissipation peak, is introduced. For favorable parameter values, this latter rheology leads to despinning. However, only models associated with large amounts of short-lived radiogenic elements lead to the observed flattening. This suggests that the accretion process needs to be completed shortly after the formation of CAIs (Calcium-Aluminum-rich Inclusions) (?4 Myr). For [26Al] varying between 72 and 46 ppb, the observed flattening is obtained only for a limited range of initial spin period, between 9.5 and 10.2 h. For [26Al] ranging between 30 and 15 ppb, initial spin rates smaller than 8.5 h are required. For smaller values of [26Al], the body is too cold and viscous to acquire a significant flattening even if a rotation period close to the body disruption limit is considered. Even with a thin lithosphere during the early stage, our simulations show that Iapetus never reaches the equilibrium figure for a hydrostatic body due to the non-zero rigidity of the lithosphere. The 35 km value of the flattening is the result of the partial relaxation of an ancient larger flattening ranging between 45 and 80 km, depending on the evolution of the lithosphere thickness mainly controlled by the radiogenic content. A thin lithosphere is consistent with an early building of the equatorial ridge. The lithosphere thickening due to interior cooling can explain the preservation of the ridge throughout the remaining evolution of Iapetus.  相似文献   

15.
By computer simulations, the dynamical evolution of plane triple systems of gaseous protogalaxies and galaxies with zero initial velocities has been studied. Inside the regionD of initial configurations some subregions have been revealed corresponding to a coalescence of protogalaxies on the first double approach. The average spin momenta of mergers are approximately equal to those typical of disk galaxies. In triple galaxies, a coalescence on the first double approach does not occur. The presence of significant hidden mass makes the approaches wider and prevents the coalescence of bodies in the systems without a central object. A central pair in a group of galaxies aids to coalescence. Also the change during time of the virial coefficient has been investigated.  相似文献   

16.
The new target of the Rosetta mission is comet 67P/Churyumov-Gerasimenko (hereafter 67P/C-G). In order to support the planning of the mission, in particular the strategy during the mapping and landing phases, we have performed numerical simulations of the rotational evolution of a comet in the orbit of 67P/C-G. In these simulations, the currently known observational constraints have been taken into account and a large set of initial conditions were considered. For most of the simulations, we observe that the sublimation-induced torques produce significant changes in the rotational parameters of a 67P/C-G-like comet. Typical rates of change for the spin period from the rendezvous up to the end of the nominal mission range from 0.001 to depending on different circumstances as described in the text. At perihelion, rates of change of the orientation of the angular momentum vector amount to about 0.01-. These simulations suggest that a specific strategy should be defined in order to monitor likely variations of the rotational parameters. As an example we show a possible optimized schedule for observations with the OSIRIS instrument to determine the rotational parameters of comet 67P/C-G and their possible evolution.  相似文献   

17.
It was found that certain features in the observed spin vector distribution of main belt asteroids can be explained by the differences in the dynamical spin vector evolution between objects with high and low orbital inclinations. In particular, the deficiency of high-inclination objects whose spin vectors are close to the ecliptic plane can be accounted for.The present spin vector distribution of main belt asteroids is due to several factors connected with their collisional and dynamical evolution. In this paper, the influence of the orbital evolution on the spin axis of asteroids is examined in the case of 25 objects with typical main belt orbital evolution and 125 synthetic objects, during an integration over a time period of 1 Myr. This investigation produced the following general results:• The difference between maximum and minimum obliquity increases in an approximately linear fashion with increasing orbital inclination of the studied objects.• The inclination is the major factor influencing the magnitude of the obliquity variation. This variation is generally larger for asteroids with their initial spin vectors located close to the orbital plane.• In general, the regular obliquity differences are relatively insensitive to differences in the shape, composition, and spin rate of the asteroids.The result is compared with the properties of the observed spin vectors for 73 main belt asteroids and good agreement is found between the above results and the existing spin vector distribution.  相似文献   

18.
Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s6 spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.  相似文献   

19.
Numerical simulations have been used to study high velocity two-body impacts. In this paper a two-dimensional Lagrangian finite difference hydrocode and a three-dimensional smooth particle hydrocode (SPH) are described and initial results reported.

The 2D hydrocode has successfully reproduced both the fragment size distribution and the mean fragment velocities from laboratory impact experiments using basalt and cement mortar. Further, the hydrocode calculations have determined that the energy needed to fracture a body has a much stronger dependence on target size than predicted from most scaling theories. In addition, velocity distributions obtained (using homogeneous targets at impact velocities around 2 km s−1) indicate that mean ejecta speeds resulting from large-body collisions do not generally exceed escape velocities.

The SPH model provides a fully three-dimensional framework for studying impacts, so that phenomena such as oblique collisions or impacts into non-spherical targets may be studied. The gridless code allows for arbitrary levels of distortion, and is hence appropriate for modeling the large-scale deformations which accompany most impact events. Because fragments are modeled explicitly, greater numerical accuracy is achieved in the regions of large fragments than with the purely statistical approach of the 2D model. Of course, this accuracy comes at the expense of significantly greater computational requirements.

These codes can be, and have been, used to make specific predictions about particular objects in our solar system. But more significantly, they allow us to explore a broad range of collisional events. Certain parameters (size, time) can be studied only over a very restricted range within the laboratory; other parameters (initial spin, low gravity, exotic structure or composition) are difficult to study at all experimentally. The outcomes of numerical simulations lead to a more general and accurate understanding of impacts in their many forms.  相似文献   


20.
We present the results of photometric observations of the TOPEX/Poseidon satellite performed during 2008–2016. The satellite become space debris after a failure in January, 2006, in a low Earth orbit. In the Laboratory of Space Research of Uzhhorod National University 73 light curves of the spacecraft were obtained. Standardization of photometric light curves is briefly explained. We have calculated the color indices of reflecting surfaces and the spin rate change. The general tendency of the latter is described by an exponential decay function. The satellite spin periods based on 126 light curves (including 53 light curves from the MMT-9 project operating since 2014) were taken into account. In 2016 the period of its own rotation reached its minimum of 10.6 s.A method to derive the direction of the spin axis of an artificial satellite and the angles of the light scattered by its surface has been developed in the Laboratory of Space Research of Uzhhorod National University. We briefly describe the “Orientation” program used for these purposes. The orientation of the TOPEX/Poseidon satellite in mid-2016 is given. The angle of precession β = 45°–50° and period of precession P pr = 141.5 s have been defined. The reasons for the identified nature of the satellite’s own rotation have been found. They amount to the perturbation caused by a deviation of the Earth gravity field from a central-symmetric shape and the presence of moving parts on the satellite.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号