Engineering topics
Publications and source records attributed to Melosh, H. J..
Fulgurites result from transient high temperature processes, and some have extremely reduced phases. We performed both modeling and a microprobe analysis of natural fulgurites. The modeling suggests vapor phase C causes reduction of silicate liquid. Additional information is contained in the original extended abstract.
We have incorporated a newly upgraded version of the equation of state ANEOS that allows for the formation of molecular vapor into a smoothed-particle hydrodynamics code. Results of simulations of Moon-forming impacts with this new EOS are presented. Additional information is contained in the original extended abstract.
Europa's subsurface water ocean may be warm: that is, at the temperature of water's maximum density. This provides a natural explanation of chaos melt-through events and leads to a correct estimate of the age of its surface. Additional information is contained in the original extended abstract.
To determine the survivability of ejection from a planet, two uv-resistant spores have been tested for resistance to high acceleration and high change in acceleration, or jerk.
Natural impacts in which the projectile strikes the target vertically are virtually nonexistent. Nevertheless, our inherent drive to simplify nature often causes us to suppose most impacts are nearly vertical. Recent theoretical, observational, and experimental work is improving this situation, but even with the current wealth of studies on impact cratering, the effect of impact angle on the final crater is not well understood. Although craters' rims may appear circular down to low impact angles, the distribution of ejecta around the crater is more sensitive to the angle of impact and currently serves as the best guide to obliquity of impacts. Experimental studies established that crater dimensions depend only on the vertical component of the impact velocity. The shock wave generated by the impact weakens with decreasing impact angle. As a result, melting and vaporization depend on impact angle; however, these processes do not seem to depend on the vertical component of the velocity alone. Finally, obliquity influences the fate of the projectile: in particular, the amount and velocity of ricochet are a strong function of impact angle.
Prominent crater chains on Ganymede and Callisto are most likely the impact scars of comets tidally disrupted by Jupiter and are not secondary crater chains. We have examined the morphology of these chains in detail in order to place constraints on the properties of the comets that formed them and the disruption process. In these chains, intercrater spacing varies by no more than a factor of 2 and the craters within a given chain show almost no deviation from linearity (although the chains themselves are on gently curved small circles). All of these crater chains occur on or very near the Jupiter-facing hemisphere. For a given chain, the estimated masses of the fragments that formed each crater vary by no more than an order of magnitude. The mean fragment masses for all the chains vary by over four orders of magnitude (W. B. McKinnon and P. M. Schenk 1995, Geophys. Res. Lett. 13, 1829-1832), however. The mass of the parent comet for each crater chain is not correlated with the number of fragments produced during disruption but is correlated with the mean mass of the fragments produced in a given disruption event. Also, the larger fragments are located near the center of each chain. All of these characteristics are consistent with those predicted by disruption simulations based on the rubble pile cometary nucleus model (in which nuclei are composed on numerous small fragments weakly bound by self-gravity), and with those observed in Comet D/Shoemaker-Levy 9. Similar crater chains have not been found on the other icy satellites, but the impact record of disrupted comets on Callisto and Ganymede indicates that disruption events occur within the Jupiter system roughly once every 200 to 400 years.
A stress analysis of tectonic patterns near Venusian coronae is reported. Combined local corona stresses and uniform regional stresses are used to predict patterns of surface tectonic features. The patterns are compared to those of coronae on Magellan images to determine the regional stress and elastic lithospheric thickness about the coronae. Regional stresses of 0.1-0.6 kbar and elastic lithospheric thicknesses of 10 +/- 5 km are estimated for three specific coronae.
The thermal effects of giant impacts are studied by estimating the melt volume generated by the initial shock wave and corresponding magma ocean depths. Additionally, the effects of the planet's initial temperature on the generated melt volume are examined. The shock pressure required to completely melt the material is determined using the Hugoniot curve plotted in pressure-entropy space. Once the melting pressure is known, an impact melting model is used to estimate the radial distance melting occurred from the impact site. The melt region's geometry then determines the associated melt volume. The model is also used to estimate the partial melt volume. Magma ocean depths resulting from both excavated and retained melt are calculated, and the melt fraction not excavated during the formation of the crater is estimated. The fraction of a planet melted by the initial shock wave is also estimated using the model.
The hypervelocity impact that excavated the Stickney crater on Phobos is numerically modeled in order to understand the effects this collision had on the interior and surface of Phobos. The model assumes homogeneity prior to impact based on Fujiwara's (1991) results. Constitutive relations for water ice and basalt are applied in tandem models. Fracture energetics is used to show that the impact itself is sufficient to cause the cracks observed on Phobos. It is further shown that the Stickney impact itself could not have created large voids; thus, the low density of Phobos must either be compositional in nature or else the result of porosity at a scale sufficiently small to avoid scattering the impact energy substantially.
Vesta is a large asteroid whose crust is mostly basaltic. Spectral heterogeneity suggests a sizable olivine feature which may be explained as impact excavation. The spectral data probably show a localized feature approximately 200 km in diameter or a diffuse feature approximately 400 km in diameter. Lightcurve irregularities suggest heterogeneity on a similar scale. This heterogeneity may represent the crater bowl, the extent of its ejecta deposit, or indeed something unrelated to cratering. In any case, drawing direct inferences about the state of Vesta's surface on the basis of these observations involves substantial speculation. We are still far from understanding the boundary between cratering and catastrophic disruption, particularly on targets for which strength and self-gravity both matter. But we are now able to model the underlying physical process -- dynamic fragmentation -- accurately with fragmentation hydrocodes such as SALE 2D and SPH3D. We present several impact scenarios for Vesta; our study is similar to a previous impact model for the formation of Stickney crater on Phobos. We illustrate the effects of impactors of different sizes and velocities, and the effect of gravity and structural inhomogeneity.
Among the many effects of high-speed, giant impacts is widescale melting that can potentially trigger catastrophic core formation. If the projectile is sufficiently large, the melt pools to form an intact melt region. The dense phase then segregates from the melt, forming a density anomoly at the melt region's base. If the anomoly produces a differential stress larger than a certain minimum, it overcomes the mantle's long-term elastic strength and rapidly forms a core. It was previously shown that giant impacts effectively trigger core formation in silicate bodies by the time they grow to the mass of Mercury and in icy bodies by the time they grow larger than Triton. In order for this process to be viable, an intact melt region must be formed. Conditions under which this occurs is examined in more detail than previously published.
The present model for the timing and mechanisms of planetary core formation argues that once a planet reaches a certain minimum mass, the large impacts that are typical of late accretion can trigger core formation. This model circumvents the difficulties posed by the large-scale segregation of molten iron into diapirs, and the displacement of the cold, elastic interior of the planet by the iron. The analytical melting model used is based on the Hugoniot equations, the empirical relationship for the decline of particle velocity with distance, and the linear shock-particle velocity relationship.
The dynamic fragmentation in impacts into solids is examined using a physical model for the formation and growth of cracks in rocks. The physical model is then inserted into a numerical model (hydrocode) of stress wave propagation and interaction, from which the outcome of a given impact event can be computed. The hydrocode model predicts fragment sizes due to impact in terms of shock waves propagating in a homogeneous elastic medium containing a distribution of crack nucleation centers known as Weibull flaws.
A quantitative model in which the parabolic features are produced by the interaction of the zonal winds with material ejected ballistically from the impact crater is proposed. As the ejecta particles fall through the atmosphere, the winds transport them downwind from their entry point, smaller particles being transported a greater distance. Since the ejecta distribution is initially axially symmetric and smaller particles are thrown farther from the crater, the winds blow the particles on the upwind side back upon one another, leading to a pile-up of material. On the downwind side, the winds disperse the ejecta particles and no pile-up occurs. The resulting thickness distribution on the Venusian surface matches the observed parabolic features closely. The dual parabolic features associated with the crater Carson is also explained by this model.
The effect of giant impacts on the initial chemical and thermal states of the terrestrial planets is just now being explored. A large high speed impact creates an approximately hemispherical melt region with a radius that depends on the projectile's radius and impact speed. It is shown that giant impacts on large planets can create large, intact melt regions containing melt volumes up to a few times the volume of the projectile. These large melt regions are not created on asteroid sized bodies. If extruded to the surface, these regions contain enough melt to create a magma ocean of considerable depth, depending on the impact speed, projectile radius, and gravity of the target planet.
An evaluation is conducted of the possibility that the tidal disruption of a population of contact binary asteroids can account for terrestrial-impact 'doublet' craters. Detailed orbital integrations indicate that while such asteroids are often disrupted by tidal forces outside the Roche limit, the magnitude of the resulting separations is too small to account for the observed doublet craters. It is hypothesized that an initial population of km-scale earth-crossing objects encompassing 10-20 percent binaries must be responsible for doublet impacts, as may be verified by future observations of earth-approaching asteroids.
The SALE, or Simplified Arbitrary Lagrangian Eulerian code is presently used to calculate the antipodal pressures generated by large impacts on icy satellites and planets of various compositions and core sizes. While in a high density core/lower density mantle planet there is a refraction of the incident pressure wave away from the antipode by the core, resulting in a low antipodal pressure, a planet with a low velocity core with a higher-wave speed mantle will focus the waveform to maximize antipodal pressure. Specific impact basins on the moon, Mercury, and several icy satellites are modeled to ascertain whether the impacts were capable of producing the observed antipodal features.
Traditionally, laboratory impact experiments, designed to simulate asteroid collisions, attempted to establish a predictive capability for collisional outcomes given a particular set of initial conditions. Unfortunately, laboratory experiments are restricted to using targets considerably smaller than the modelled objects. It is therefore necessary to develop some methodology for extrapolating the extensive experimental results to the size regime of interest. Results are reported obtained through the use of two dimensional hydrocode based on 2-D SALE and modified to include strength effects and the fragmentation equations. The hydrocode was tested by comparing its predictions for post-impact fragment size distributions to those observed in laboratory impact experiments.