Search NASA⌕ Search

SEARCH · Search NASA

Results for “Smooth Particle Hydrodynamics (SPH)”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗
NASA NTRS1993

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗
NASA NTRS1993

Dissipative merging of galaxies

The galaxy merging is investigated with hydrodynamical processes taken into account. For this purpose, the 3D calculations are performed by the use of a smoothed particle hydrodynamics (SPH) scheme combined with an N-body scheme. In these calculations, we find a new merging criterion and the dependence of the central phase space density of merger remnants upon the gas fraction in progenitors. It is concluded that ellipticals can be formed just by merging of fairly gas-rich primordial galaxies, not ordinary spiral galaxies.

Umemura, M.↗
NASA NTRS1993

Smoothed particle hydrodynamics with GRAPE-1A

We describe the implementation of a smoothed particle hydrodynamics (SPH) scheme using GRAPE-1A, a special-purpose processor used for gravitational N-body simulations. The GRAPE-1A calculates the gravitational force exerted on a particle from all other particles in a system, while simultaneously making a list of the nearest neighbors of the particle. It is found that GRAPE-1A accelerates SPH calculations by direct summation by about two orders of magnitudes for a ten thousand-particle simulation. The effective speed is 80 Mflops, which is about 30 percent of the peak speed of GRAPE-1A. Also, in order to investigate the accuracy of GRAPE-SPH, some test simulations were executed. We found that the force and position errors are smaller than those due to representing a fluid by a finite number of particles. The total energy and momentum were conserved within 0.2-0.4 percent and 2-5 x 10 exp -5, respectively, in simulations with several thousand particles. We conclude that GRAPE-SPH is quite effective and sufficiently accurate for self-gravitating hydrodynamics.

Umemura, Masayuki↗
NASA NTRS1993

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗
NASA NTRS1994

Shoemaker-Levy 9 and the tidal disruption of comets

The break-up of Periodic Comet Shoemaker-Levy 9 into multiple pieces following its grazing encounter with Jupiter in July 1992 can be used to study tidally-induced fracture in comets. This spectacular event allows us not only to set limits on the size, strength and density of Shoemaker-Levy 9 itself, but provides invaluable guidance to numerical modeling of such encounters. In an extensive treatment of tidal breakup which assumed self-gravitating, homogeneous, perfectly elastic bodies, Dobrovolskis derived simple analytical expressions for the tidally-induced surface and central stresses. Both can be cast in such a way that Poisson's ratio is the only material dependent constant entering these expressions. Whether both surface and central failure must be initiated as a criterion for breakup, or either one of them is sufficient, remains a subject of disagreement. To resolve this debate, we model the details of cometary breakup using a three-dimensional Smooth Particle Hydrodynamics (SPH) code modified to simulate fracture in small solid objects. At lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution.

Benz, W.↗
NASA NTRS1994

Merging neutron stars. 1. Initial results for coalescence of noncorotating systems

We present three-dimensional Newtonian simulations of the coalescence of two neutron stars, using a smoothed particle hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes (gamma = 2.4) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates result in greater mass loss and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of 2.4 solar mass, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of approximately 10 MeV, while the disk is heated by shocks to a temperature of 2 to 4 MeV.

Davies, M. B.↗
NASA NTRS1994

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗
NASA NTRS1994

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Artymowicz, Pawel↗
NASA NTRS1995

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗
NASA NTRS1994

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
NASA NTRS1994

Postexplosion hydrodynamics of supernovae in red supergiants

Shock propagation, mixing, and clumping are studied in the explosion of red supergiants as Type II supernovae using a two-dimensional smooth particle hydrodynamic (SPH) code. We show that extensive Rayleigh-Talor instabilities develop in the ejecta in the wake of the reverse shock wave. In all cases, the shell structure of the progenitor is obliterated to leave a clumpy, well-mixed supernova remnant. However, the occurrence of mass loss during the lifetime of the progenitor can significantly reduce the amount of mixing. These results are independent of the Type II supernova explosion mechanism.

Herant, Marc↗
NASA NTRS1994

Hydrodynamics of binary coalescence. 1: Polytropes with stiff equations of state

We performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations on the verge of dynamical instability. We consider synchronized equilbria only and concentrate on stiff equations of state, with adiabatic Gamma greater than 5/3. We assume that the polytropic constants (K identically equal to P/(rho(exp Gamma)) are the same for both components. These conditions apply well to models of neutron star binaries. Accordingly, we discuss our results in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, and we calculate the emission of gravitational radiation in the quadruple approximation. The fully nonlinear development of the instability is followed using SPH until a new equilibrium configuration is reached by the system. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An axisymmetric merged configuration is always produced when the adiabatic exponent Gamma approximately less than 2.3. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, triaxial merged configurations are obtained when Gamma approximately greater than 2.3, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios q not equal to 1 typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. For a mass ratio q approximately less than 0.5, however, dynamical mass transfer can temporarily retard the coalescence by causing a rapid reexpansion of the binary into a new, slightly eccentric but dynamically stable orbit. The maximum amplitude h(sub max) and peak luminosity L(sub max) of the gravitational waves emitted during the final coalescence are nearly independent of Gamma, but depend sensitively on the mass ratio q. The approximate scalings we find are h(sub max) varies as q(exp 2) and L(sub max) varies as q(exp 6) for q close to unity.

Rasio, Frederic A.↗
NASA NTRS1994

Impact on comet Shoemaker-Levy 9 on Jupiter

Three-dimensional numerical simulations of the impact of Comet Shoemaker - Levy 9 on Jupiter and the resulting vapor plume expansion were conducted using the Smoothed Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 km to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disk. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km, but no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Ahrens, Thomas J.↗
NASA NTRS1994

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗
NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗
NASA NTRS1994

Close encounters of the third-body kind

We simulated encounters involving binaries of two eccentricities: e = 0 (i.e., circular binaries) and e = 0.5. In both cases the binary contained a point mass of 1.4 solar masses (i.e., a neutron star) and a 0.8 solar masses main-sequence star modeled as a polytrope. The semimajor axes of both binaries were set to 60 solar radii (0.28 AU). We considered intruders of three masses: 1.4 solar masses (a neutron star), 0.8 solar masses (a main-sequence star or a higher mass white dwarf), and 0.64 solar masses (a more typical mass white dwarf). Our strategy was to perform a large number (40,000) of encounters using a three-body code, then to rerun a small number of cases with a three-dimensional smoothed particle hydrodynamics (SPH) code to determine the importance of hydrodynamical effects. Using the results of the three-body runs, we computed the exchange across sections, sigma(sub ex). From the results of the SPH runs, we computed the cross sections for clean exchange, denoted by sigma(sub cx); the formation of a triple system, denoted by sigma(sub trp); and the formation of a merged binary with an object formed from the merger of two of the stars left in orbit around the third star, denoted by sigma(sub mb). For encounters between either binary and a 1.4 solar masses neutron star, sigma(sub cx) approx. 0.7 sigma(sub ex) and sigma(sub mb) + sigma(sub trp) approx. 0.3 sigma(sub ex). For encounters between either binary and the 0.8 solar masses main-sequence star, sigma(sub cx) approx. 0.50 sigma(sub ex) and sigma(sub mb) + sigma(sub trp) approx. 1.0 sigma(sub ex). If the main sequence star is replaced by a main-sequence star of the same mass, we have sigma(sub cx) approx. 0.5 sigma(sub ex) and sigma(sub mb) + sigma(sub trp) approx. 1.6 sigma(sub ex). Although the exchange cross section is a sensitive function of intruder mass, we see that the cross section to produce merged binaries is roughly independent of intruder mass. The merged binaries produced have semi-major axes much larger than either those of the original binaries or those of binaries produced in clean exchanges. Coupled with their lower kick velocities, received from the encounters, their larger size will enhance their cross section, shortening the waiting time to a subsequent encounter with another single star.

Davies, M. B.↗
NASA NTRS1994

Impact of comet Shoemaker-Levy 9 on Jupiter

We have employed three-dimensional numerical simulations of the impact of Comet Shoemaker-Levy 9 (SL9) on Jupiter and the resulting vapor plume expansion using the smoothed particle hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disc. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km; however, no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Takata, Toshiko↗