Engineering topics

Steinberg, Elad

Publications and source records attributed to Steinberg, Elad.

5 records

DOE OSTI2015

SPINS OF LARGE ASTEROIDS: A HINT OF A PRIMORDIAL DISTRIBUTION IN THEIR SPIN RATES

The Asteroid Belt and the Kuiper Belt are relics from the formation of our solar system. Understanding the size and spin distribution of the two belts is crucial for a deeper understanding of the formation of our solar system and the dynamical processes that govern it. In this paper, we investigate the effect of collisions on the evolution of the spin distribution of asteroids and KBOs. We find that the power law nature of the impactors’ size distribution leads to a Lévy distribution of the spin rates. This results in a power law tail in the spin distribution, in stark contrast to the usually quoted Maxwellian distribution. We show that for bodies larger than 10 km, collisions alone lead to spin rates peaking at 0.15–0.5 revolutions per day. Comparing that to the observed spin rates of large asteroids (R > 50 km), we find that the spins of large asteroids, peaking at ∼1–2 revolutions per day, are dominated by a primordial component that reflects the formation mechanism of the asteroids. Similarly, the Kuiper Belt has undergone virtually no collisional spin evolution, assuming current densities. Collisions contribute a spin rate of ∼0.01 revolutions per day, thus the observed fast spin rates of KBOs are also primordial in nature.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH

We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov’s method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangian simulation is better than a Eulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

THE GALACTIC CENTER CLOUD G2 AND ITS GAS STREAMER

We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late 2013 August, 2014 April, and 2014 July. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. The cloud reached minimum distance to the MBH of 1950 Schwarzschild radii. As expected for an observation near the pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to that of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (1) G1 was originally on the G2 orbit preceding G2 by 13 yr, and (2) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

BALANCING THE LOAD: A VORONOI BASED SCHEME FOR PARALLEL COMPUTATIONS

One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs. We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a two-dimensional moving mesh hydrodynamical code, but it can be extended trivially to other codes in two or three dimensions. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

BINARY YORP EFFECT AND EVOLUTION OF BINARY ASTEROIDS

The rotation states of kilometer-sized near-Earth asteroids are known to be affected by the Yarkevsky O'Keefe-Radzievskii-Paddack (YORP) effect. In a related effect, binary YORP (BYORP), the orbital properties of a binary asteroid evolve under a radiation effect mostly acting on a tidally locked secondary. The BYORP effect can alter the orbital elements over {approx}10{sup 4}-10{sup 5} years for a D{sub p} = 2 km primary with a D{sub s} = 0.4 km secondary at 1 AU. It can either separate the binary components or cause them to collide. In this paper, we devise a simple approach to calculate the YORP effect on asteroids and the BYORP effect on binaries including J{sub 2} effects due to primary oblateness and the Sun. We apply this to asteroids with known shapes as well as a set of randomly generated bodies with various degrees of smoothness. We find a strong correlation between the strengths of an asteroid's YORP and BYORP effects. Therefore, statistical knowledge of one could be used to estimate the effect of the other. We show that the action of BYORP preferentially shrinks rather than expands the binary orbit and that YORP preferentially slows down asteroids. This conclusion holds for the two extremes of thermal conductivities studied in this work and the assumption that the asteroid reaches a stable point, but may break down for moderate thermal conductivity. The YORP and BYORP effects are shown to be smaller than could be naively expected due to near cancellation of the effects at small scales. Taking this near cancellation into account, a simple order-of-magnitude estimate of the YORP and BYORP effects as a function of the sizes and smoothness of the bodies is calculated. Finally, we provide a simple proof showing that there is no secular effect due to absorption of radiation in BYORP.

79 ASTRONOMY AND ASTROPHYSICS↗