共查询到20条相似文献，搜索用时 6 毫秒
1.
Numerical simulation of impact cratering on granular material 总被引：1，自引：0，他引：1
A new numerical code based on the Distinct Element Method (DEM) is developed to study the impact cratering processes on granular material. This code has a potential advantage to simulate the cratering process on granular material, since the movement of discrete particles can be treated. To show the physical plausibility of this code, we conduct 3-D numerical simulations of vertical impact into granular material targets that consist of 384,000 particles, and compare the results with those from experimental studies. It is shown that the excavation stage of cratering derived from experimental studies is represented well by our simulation: the size of the crater cavity, and the ejecta velocity and angle distributions are consistent with those obtained in laboratory experiments. The impact simulation code developed in this study is thus suggested to be useful for the analysis of the impact cratering process on granular material. 相似文献
2.
We investigate planetary fly-bys of asteroids using an approximate volume-averaged method that offers a relatively simple, but very flexible, approach to study the rotational dynamics of ellipsoids. The asteroid is considered to be a deformable, prolate ellipsoid, with its interior being modeled as a rigid-granular material. Effects due to the asteroid's rotation, its self-gravity and gravitational interaction with the planet are included. Using a simplified approach allows us to explore in detail the mechanics of asteroid's deformations and disruptions during planetary encounters. We also compare our results with those obtained by Richardson et al. [Richardson, D.C., Bottke Jr., W.F., Love, S.G., 1998. Icarus 134, 47-76] who used a large numerical code. We find that many of the features reported by them can indeed be captured by our rather simple methodology, and we discuss the reasons why some of our results differ from theirs. 相似文献
3.
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. 相似文献
4.
《New Astronomy》2016
We present a GPU accelerated CUDA-C implementation of the Barnes Hut (BH) tree code for calculating the gravitational potential on octree adaptive meshes. The tree code algorithm is implemented within the FLASH4 adaptive mesh refinement (AMR) code framework and therefore fully MPI parallel. We describe the algorithm and present test results that demonstrate its accuracy and performance in comparison to the algorithms available in the current FLASH4 version. We use a MacLaurin spheroid to test the accuracy of our new implementation and use spherical, collapsing cloud cores with effective AMR to carry out performance tests also in comparison with previous gravity solvers. Depending on the setup and the GPU/CPU ratio, we find a speedup for the gravity unit of at least a factor of 3 and up to 60 in comparison to the gravity solvers implemented in the FLASH4 code. We find an overall speedup factor for full simulations of at least factor 1.6 up to a factor of 10. 相似文献
5.
《Icarus》1998,132(1):113-124
We present results of two-dimensional gravitationalN-body simulations of the late stage of planetary formation. This stage is characterized by the direct accretion of hundreds of lunar-sized planetesimals into planetary bodies. Our simulation code is based on the Hermite Individual Timestep integration algorithm, and gravitational interactions among all bodies are included throughout the simulations. We compare our simulation with earlier works that do not include all interactions, and we find very good agreement. A previously published collisional fragmentation model is included in our simulation to study the effects of the production of fragments on the subsequent evolution of the larger planetary bodies. It is found that for realistic two-body collisions that, according to this model, both bodies will suffer fragmentation, and that the outcome of the collision will be a relatively large core containing most of the mass and a few small fragments. We present the results of simulations that include this simple fragmentation model. They indicate that the presence of small fragments have only a small effect on the growth or orbital evolution of the large planet-sized bodies. 相似文献
6.
Krzysztof Goz´dziewski Eric Bois rzej J. Maciejewski 《Monthly notices of the Royal Astronomical Society》2002,332(4):839-855
The Gliese 876 planetary system consists of two Jupiter-like planets having a nearly commensurate 2:1 orbital periods ratio. Because the semimajor axes of the planets are very small (of the order 0.1 au and 0.2 au, respectively), and the eccentricity of the inner companion is ≃0.3, the mutual perturbations are extremely large. However, many authors claim the long-term orbital stability of the system, at least over 500 Myr for initial conditions found by Rivera & Lissauer. Results of investigations of a migration of initially separated planets into the close 2:1 mean motion resonance lock from Lee & Peale also support the conclusion that the system should be stable for the lifetime of the parent star. Initial conditions of the system, found from non-linear N -body fits by Laughlin & Chambers and Rivera & Lissauer, to the radial velocity curve, formally allow for a variety of orbital configurations of the GJ 876 system, e.g. coplanar, with planetary inclinations in the range [≃30°, 90°], and with relative inclinations of orbital planes as high as 80°. Our work is devoted to the stability investigation of the systems originating from the fitted initial conditions. We study neighbourhoods of these initial states in the orbital parameter space. We found estimations of the 2:1 mean motion resonance width and dynamical limitations on the planetary masses. We also obtain a global representation of the domains of the orbital parameters space in which initial conditions leading to stable evolutions can be found. Our results can be useful in localization of the best, stable fits to the observational data. In our investigations we use the MEGNO technique (the Mean Exponential Growth factor of Nearby Orbits) invented by Cincotta & Simó. It allows us to distinguish efficiently and precisely between chaotic and regular behaviour of a planetary system. 相似文献
7.
Impacts into quartz sand: Crater formation,shock metamorphism,and ejecta distribution in laboratory experiments and numerical models 
下载免费PDF全文
下载免费PDF全文 We investigated the ejection mechanics by a complementary approach of cratering experiments, including the microscopic analysis of material sampled from these experiments, and 2‐D numerical modeling of vertical impacts. The study is based on cratering experiments in quartz sand targets performed at the NASA Ames Vertical Gun Range. In these experiments, the preimpact location in the target and the final position of ejecta was determined by using color‐coded sand and a catcher system for the ejecta. The results were compared with numerical simulations of the cratering and ejection process to validate the iSALE shock physics code. In turn the models provide further details on the ejection velocities and angles. We quantify the general assumption that ejecta thickness decreases with distance according to a power‐law and that the relative proportion of shocked material in the ejecta increase with distance. We distinguish three types of shock metamorphic particles (1) melt particles, (2) shock lithified aggregates, and (3) shock‐comminuted grains. The agreement between experiment and model was excellent, which provides confidence that the models can predict ejection angles, velocities, and the degree of shock loading of material expelled from a crater accurately if impact parameters such as impact velocity, impactor size, and gravity are varied beyond the experimental limitations. This study is relevant for a quantitative assessment of impact gardening on planetary surfaces and the evolution of regolith layers on atmosphereless bodies. 相似文献
8.
Abstract— We present a numerical simulation technique for the retrieval of the electric properties relative permittivity and conductivity of planetary, asteroid, and cometary regolith. Our simulation techniques aim at accompanying hardware development and conducting virtual experiments, e.g., to assess the response of arbitrary heterogeneous conductivity and permittivity distributions or to scrutinize possibilities for spatial reconstruction methods using inverse schemes. In a first step, we have developed a finite element simulation code on the basis of unstructured, adaptive triangular grids for arbitrary two‐dimensional axisymmetric distributions of conductivity and permittivity. The code is able to take into account the spatial geometry of the probe and allows for possible inductive effects. In previous studies, the non‐inductive approach has been used to convert potential and phase data into apparent material properties. By our simulations, we have shown that this approach is valid for the frequency range from 102 Hz to 107 Hz and electric conductivities of 10?8 S/m that are typical for the near‐surface region of asteroids and comets composed of chondritic materials and/or frozen volatiles such as H2O and CO2 ice. We prove the accuracy of our code to be better than 10%, using mixed types of boundary conditions and present a simulated vertical log through a horizontally stratified subsurface layer as a representative example of a heterogeneous distribution of the electrical properties. Resolution studies for the given electrode separation reveal that the material parameters of layers having thicknesses of less than about half the electrode spread are not reconstructible if only apparent quantities are considered. Therefore, spatial distributions of the complex sensitivity are presented having in mind a future data inversion concept that will permit the multi‐dimensional reconstruction of material parameters in heterogeneous environments. 相似文献
9.
Joshua E. Colwell Stein Sture Dan Durda Tyler Goudie Daniel J. Ashcom Thomas Keohane Michael Lupton 《Icarus》2008,195(2):908-917
Collisions between planetary ring particles and in some protoplanetary disk environments occur at speeds below 10 m/s. The particles involved in these low-velocity collisions have negligible gravity and may be made of or coated with smaller dust grains and aggregates. We undertook microgravity impact experiments to better understand the dissipation of energy and production of ejecta in these collisions. Here we report the results of impact experiments of solid projectiles into beds of granular material at impact velocities from 0.2 to 2.3 m/s performed under near-weightless conditions on the NASA KC-135 Weightless Wonder V. Impactors of various densities and radii of 1 and 2 cm were launched into targets of quartz sand, JSC-1 lunar regolith simulant, and JSC-Mars-1 martian regolith simulant. Most impacts were at normal or near-normal incidence angles, though some impacts were at oblique angles. Oblique impacts led to much higher ejection velocities and ejecta masses than normal impacts. For normal incidence impacts, characteristic ejecta velocities increase with impactor kinetic energy, KE, as approximately KE0.5. Ejecta masses could not be measured accurately due to the nature of the experiment, but qualitatively also increased with impactor kinetic energy. Some experiments were near the threshold velocity of 0.2 m/s identified in previous microgravity impact experiments as the minimum velocity needed to produce ejecta [Colwell, J.E., 2003. Icarus 164, 188-196], and the experimental scatter is large at these low speeds in the airplane experiment. A more precise exploration of the transition from low-ejecta-mass impacts to high-ejecta-mass impacts requires a longer and smoother period of reduced gravity. Coefficient of restitution measurements are not possible due to the varying acceleration of the airplane throughout the experiment. 相似文献
10.
Adrián Brunini Héctor R. Viturro 《Monthly notices of the Royal Astronomical Society》2003,346(3):924-932
We present a tree code for simulations of collisional systems dominated by a central mass. We describe the implementation of the code and the results of some test runs with which the performance of the code was tested. A comparison between the behaviour of the tree code and a direct hybrid integrator is also presented. The main result is that tree codes can be useful in numerical simulations of planetary accretion, especially during intermediate stages, where possible runaway accretion and dynamical friction lead to a population with a few large bodies in low-eccentricity and low-inclination orbits embedded in a large swarm of small planetesimals in rather excited orbits. Some strategies to improve the performance of the code are also discussed. 相似文献
11.
The nonlinear motion (collapse/explosion) of a homogeneous, rapidly rotating, gaseous ellipsoid has been studied with the effects of an external gravitational field incorporated. The non-axisymmetric motion has been followed using an efficient numerical code based on certain properties of the potential functions. The gravitational radiation associated with the phenomenon was calculated, and the wave forms studied to extract information on the dynamics of the source. Possible implications of our results for various astrophysical processes are discussed. 相似文献
12.
Yi Zhao-hua 《Chinese Astronomy and Astrophysics》2001,25(4):399-408
This is a discussion and personal views on current important topics and methods in solar system dynamics. The topics include dynamical models, orbit resonance, planetary rings, chaos and secular evolution, motion of near-earth asteroids, the Kuiper belt, gravitational theory of the solar system and other related problems. 相似文献
13.
We investigate the excitation of magnetoacoustic–gravity waves generated from localized pulses in the gas pressure as well as in the vertical component of velocity. These pulses are initially launched at the top of the solar photosphere, which is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical, and oblique to the gravitational force. We numerically model magnetoacoustic–gravity waves by implementing a realistic (VAL-C) model of the solar temperature. We solve the two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic–gravity waves, which are later reflected from the transition region due to the large-temperature gradient. The cavities in the lower solar atmosphere are found to have the best conditions to act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics. 相似文献
14.
Fathi Namouni 《Astrophysics and Space Science》2011,331(2):575-595
Dynamical friction arises from the interaction of a perturber and the gravitational wake it excites in the ambient medium.
We study the effects of the presence of a boundary on dynamical friction by studying analytically the interaction of perturber
with uniform rectilinear motion in a uniform homogeneous medium with a reflecting planar boundary. Wake reflection at a medium’s
boundary may occur at the edges of truncated disks perturbed by planetary or stellar companions as well as in numerical simulations
of planet-disk interaction with no-outflow boundary conditions. In this paper, we show that the presence of the boundary modifies
the behaviour of dynamical friction significantly. We find that perturbers are invariably pushed away from the boundary and
reach a terminal subsonic velocity near Mach 0.37 regardless of initial velocity. Dynamical friction may even be reversed
for Mach numbers less than 0.37 thereby accelerating instead of decelerating the perturber. Perturbers moving parallel to
the boundary feel additional friction orthogonal to the direction of motion that is much stronger than the standard friction
along the direction of motion. These results indicate that the common use of the standard Chandrasekhar formula as a short
hand estimate of dynamical friction may be inadequate as observed in various numerical simulations. 相似文献
15.
W. Kley 《Celestial Mechanics and Dynamical Astronomy》2003,87(1-2):85-97
We study the evolution of a system consisting of two protoplanets still embedded in a protoplanetary disk. Results of two different numerical approaches are presented. In the first kind of model the motion of the disk material is followed by fully viscous hydrodynamical simulations, and the planetary motion is determined by N-body calculations including exactly the gravitational potential from the disk material. In the second kind we only solve the N-body part and add additional analytically given forces which model the effect of the torques of the disk. This type of modeling is of course orders of magnitudes faster than the full hydro-model. Another advantage of this two-fold approach is the possibility of adjusting the otherwise unknown parameters of the simplified model.The results give very good agreement between the methods. Using two different initial setups for the planets and disk, we obtain in the first case a resonant trapping into the 3:1 resonance, and in the second case a trapping into the 2:1 resonance. Resonant capture leads to a rise in the eccentricity and to an alignment of the spatial orientation of orbits. The characteristics of the numerical results agree very favorably with those of three observed planetary systems (GJ 876, HD 82943, and 55 Cnc) known to be in mean motion resonances. 相似文献
16.
17.
Modeling of migrating grains on asteroid’s surface 总被引：1，自引：0，他引：1
18.
Inge Loes TEN KATE James R. C. GARRY Zan PEETERS Richard QUINN Bernard FOING Pascale EHRENFREUND 《Meteoritics & planetary science》2005,40(8):1185-1193
Abstract— In the framework of international planetary exploration programs, several space missions are planned to search for organics and bio‐signatures on Mars. Previous attempts have not detected any organic compounds in the Martian regolith. It is therefore critical to investigate the processes that may affect organic molecules on and below the planet's surface. Laboratory simulations can provide useful data about the reaction pathways of organic material at Mars' surface. We have studied the stability of amino acid thin films against ultraviolet (UV) irradiation and use those data to predict the survival time of these compounds on and in the Martian regolith. We show that thin films of glycine and D‐alanine are expected to have half‐lives of 22 ± 5 hr and of 3 ± 1 hr, respectively, when irradiated with Mars‐like UV flux levels. Modelling shows that the half‐lives of the amino acids are extended to the order of 107 years when embedded in regolith. These data suggest that subsurface sampling must be a key component of future missions to Mars dedicated to organic detection. 相似文献
19.
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. 相似文献
20.
Ved Mitra 《Astrophysics and Space Science》1976,39(2):387-396
The formation of the solar nebula and the distribution of mass in its planetary system is studied. The underlying idea is that the protosun, fragmented out from an interstellar cloud as a result of cluster formation, gathered the planetary material and, hence, spin angular momentum by gravitational accretion during its orbital motion around the centre of the Galaxy. The study gives the initial angular momentum of the solar nebula nearly equal to the present value of the solar system. 相似文献