High-performance direct gravitational N-body simulations on graphics processing units

We present the results of gravitational direct N-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the N-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor–corrector Hermite integrator with block time steps. We find that for a large number of particles (N  10⁴) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in 10⁶ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around 10⁻¹¹, whereas for the GPU machine precision saturates around 10⁻⁶ For N  10⁶ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.     

The primordial excitation and clearing of the asteroid belt—Revisited

We have performed new simulations of two different scenarios for the excitation and depletion of the primordial asteroid belt, assuming Jupiter and Saturn on initially circular orbits as predicted by the Nice Model of the evolution of the outer Solar System [Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A., 2005. Nature 435, 466-469; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R., 2005. Nature 435, 462-465]. First, we study the effects of sweeping secular resonances driven by the depletion of the solar nebula. We find that these sweeping secular resonances are incapable of giving sufficient dynamical excitation to the asteroids for nebula depletion timescales consistent with estimates for solar-type stars, and in addition cannot cause significant mass depletion in the asteroid belt or produce the observed radial mixing of different asteroid taxonomic types. Second, we study the effects of planetary embryos embedded in the primordial asteroid belt. These embedded planetary embryos, combined with the action of jovian and saturnian resonances, can lead to dynamical excitation and radial mixing comparable to the current asteroid belt. The mass depletion driven by embedded planetary embryos alone, even in the case of an eccentric Jupiter and Saturn, is roughly 10-20× less than necessary to explain the current mass of the main belt, and thus a secondary depletion event, such as that which occurs naturally in the Nice Model, is required. We discuss the implications of our new simulations for the dynamical and collisional evolution of the main belt.     

Thermodynamics of rare events and impulsive relaxation events in the magnetospheric substorm dynamics

The magnetosphere dynamics shows fast relaxation events following power-law distribution for many observable quantities during magnetic substorms. The emergence of such power-law distributions has been widely discussed in the framework of self-organized criticality and/or turbulence. Here, a different approach to the statistical features of these impulsive dynamical events is proposed in the framework of the thermodynamics of rare events [Lavenda, B.H., Florio, A., 1992. Thermodynamics of rare events, Int. J. Theor. Phys. 31, 1455–1475; Lavenda, B.H., 1995. Thermodynamics of Extremes. Albion]. In detail, an application of such a novel approach to the magnetospheric substorm avalanching dynamics as monitored by the auroral electroject index is discussed.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.     

Analysis of loads,motions and cavity dynamics during freefall wedges vertically entering the water surface

In this paper, theoretical models are developed and numerical methods are used to analyze the loads, motions and cavity dynamics for freefall wedges with different deadrise angles vertically entering the water surface at Froude numbers: 1 ≤ Fn < 9. The time evolutions of the penetration depth, the velocity and the acceleration are analyzed and expressed explicitly. The maximum and average accelerations are predicted. The theoretical results are compared with numerical data obtained through a single-fluid BEM model with globally satisfactory agreement. The evolution of the pressures on the impact side is investigated. Before flow separation, gravity and the acceleration of the wedge have negligible influence on the pressure on the impact side for large Froude numbers or small deadrise angles; with increasing the deadrise angle or decreasing Froude number, the effects of gravity and the acceleration of the wedge tend to become more important. Global loads, with the main emphasis on the drag coefficient, are also studied. It is found that for the light wedge, the transient drag coefficient has slow variation in the first half of the collapse stage and rapid variation in the last half of the collapse stage. For the heavy wedge, the transient drag coefficients vary slowly during the whole collapse stage and can be treated as constant. The characteristics of the transient cavity during its formation are investigated. The non-dimensional pinch-off time, pinch-off depth and submergence depth at pinch-off scale roughly linearly as the Froude number.     

Ecohydrographic constraints on biodiversity and distribution of phytoplankton and zooplankton in coral reefs of the Red Sea,Saudi Arabia

An integral concept of ecological research is the constraint of biodiversity along latitudinal and environmental gradients. The Red Sea features a natural example of a latitudinal gradient of salinity, temperature and nutrient richness. Coral reefs along the Red Sea coasts are supported with allochthonous resources such as oceanic and neritic phytoplankton and zooplankton; however, relatively little is known about how the ecohydrography correlates with plankton biodiversity and abundance. In this article we present the biodiversity of phytoplankton and zooplankton in Red Sea coral reefs. Oceanographic data (temperature, salinity), water samples for nutrient analysis, particulate organic matter, phytoplankton and zooplankton, the latter with special reference to Copepoda (Crustacea), were collected at nine coral reefs over ~1500 km distance along the Red Sea coast of Saudi Arabia. The trophic state of ambient waters [as indicated by chlorophyll a (Chl a)] changed from strong oligotrophy in the north to mesotrophy in the south and was associated with increasing biomasses of Bacillariophyceae, picoeukaryotes and Synechococcus as indicated by pigment fingerprinting (CHEMTAX) and flow cytometry. Net‐phytoplankton microscopy revealed a Trichodesmium erythraeum (Cyanobacteria) bloom north of the Farasan Islands. Several potentially harmful algae, including Dinophysis miles and Gonyaulax spinifera (Dinophyceae), were encountered in larger numbers in the vicinity of the aquaculture facilities at Al Lith. Changes in zooplankton abundance were mainly correlated to the phytoplankton biomass following the latitudinal gradient. The largest zooplankton abundance was observed at the Farasan Archipelago, despite high abundances of copepodites, veligers (Gastropoda larvae) and Chaetognatha at Al Lith. Although the community composition changed over latitude, biodiversity indices of phytoplankton and zooplankton did not exhibit a systematic pattern. As this study constitutes the first current account of the plankton biodiversity in Red Sea coral reefs at a large spatial scale, the results will be informative for ecosystem‐based management along the coastline of Saudi Arabia.     

Modelling of sediment movement in the surf and swash zones

Under the action of marine currents, non-cohesive sediments evolve by bed-load, by saltation or suspension depending on their granulometry. Several authors have considered that the movement of sediment...     

A global map to aid the identification and screening of critical habitat for marine industries

Marine industries face a number of risks that necessitate careful analysis prior to making decisions on the siting of operations and facilities. An important emerging regulatory framework on environmental sustainability for business operations is the International Finance Corporation’s Performance Standard 6 (IFC PS6). Within PS6, identification of biodiversity significance is articulated through the concept of “Critical Habitat”, a definition developed by the IFC and detailed through criteria aligned with those that support internationally accepted biodiversity designations. No publicly available tools have been developed in either the marine or terrestrial realm to assess the likelihood of sites or operations being located within PS6-defined Critical Habitat. This paper presents a starting point towards filling this gap in the form of a preliminary global map that classifies more than 13 million km² of marine and coastal areas of importance for biodiversity (protected areas, Key Biodiversity Areas [KBA], sea turtle nesting sites, cold- and warm-water corals, seamounts, seagrass beds, mangroves, saltmarshes, hydrothermal vents and cold seeps) based on their overlap with Critical Habitat criteria, as defined by IFC. In total, 5798×10³ km² (1.6%) of the analysis area (global ocean plus coastal land strip) were classed as Likely Critical Habitat, and 7526×10³ km² (2.1%) as Potential Critical Habitat; the remainder (96.3%) were Unclassified. The latter was primarily due to the paucity of biodiversity data in marine areas beyond national jurisdiction and/or in deep waters, and the comparatively fewer protected areas and KBAs in these regions. Globally, protected areas constituted 65.9% of the combined Likely and Potential Critical Habitat extent, and KBAs 29.3%, not accounting for the overlap between these two features. Relative Critical Habitat extent in Exclusive Economic Zones varied dramatically between countries. This work is likely to be of particular use for industries operating in the marine and coastal realms as an early screening aid prior to in situ Critical Habitat assessment; to financial institutions making investment decisions; and to those wishing to implement good practice policies relevant to biodiversity management. Supplementary material (available online) includes other global datasets considered, documentation and justification of biodiversity feature classification, detail of IFC PS6 criteria/scenarios, and coverage calculations.