Astrophysical Jet Simulations: Comparing Different Numerical Methods

High resolution Godunov-type methods have gained increasing popularity in the last decade due to their ability in modelling highly supersonic flows. Most of these shock-capturing schemes are based on the solution of the Riemann problem between discontinuous data at each cell's interface. The purpose of this work is to investigate the role played by different Riemann solvers in the evolution of a supersonic jet. Morphology, dynamics, conservation properties and cocoon characteristics are examined for the four different cases presented.     

Annual trace element cycles in calcite–aragonite speleothems: evidence of drought in the western Mediterranean 1200–1100 yr BP

Each of two calcitic stalagmites from Grotte de Clamouse, Herault, southern France, displays a discrete aragonite layer dated at around 1100 yr BP. The layer of fanning aragonite ray crystals is immediately preceded by calcite with Mg and Sr compositions that are uniquely high for the past 3 kyr. Trace element compositions close to the boundary between original aragonite and calcite are consistent with quasi‐equilibrium partitioning of trace elements between the phases. Study of modern dripwaters demonstrates that pronounced covariation of Mg/Ca and Sr/Ca ratios in dripwater occurs owing to large amounts of calcite precipitation upflow of the drips that fed the stalagmites. Trace element to Ca ratios are enhanced during seasonally dry periods. Ion microprobe data demonstrate a pronounced covariation of trace elements, including Mg and Sr in calcite, and Sr, U and Ba in aragonite. The mean peak spacing is close to the long‐term mean of annual growth rates determined by differences in U‐series ages and so the trace element peaks are interpreted as annual. The trace element chemistry of the stalagmites on annual to inter‐annual scales thus directly reflects the amounts of prior calcite precipitation, interpreted as an index of aridity. The longer‐term context is a multi‐decadal period of aridity (1200–1100 yr BP) possibly correlated with an analogous episode in Central America. The arid period culminated in the nucleation of aragonite, but within a decade was followed by a return to precursor conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

Orbital and mission planning constraints for the deflection of NEOs impacting on Earth

This paper is the third in a series. Paper 1 presented the results of numerical modeling of deflections of NEOs in route of collision with the Earth. The model was applied to a variety of dynamical cases including both asteroidal and cometary NEOs. Paper 2 introduced the concept of “distributed deflection,” i.e., the possibility to provide the ΔV necessary to deflect an object with a succession of maneuvers each of which would have been insufficient per se to obtain the desired result. In both papers no assumptions were made on the physical composition and structure of the NEO, nor on the details of the possible deflection maneuvers from the point of view of mission analysis. Moreover, ΔV-plots were computed assuming only along-track impulses (both in the positive and negative directions), because it is easy to demonstrate that in general this is energetically the most favorable configuration. Also in the present paper no assumptions were made on the physical composition and structure of the NEO, even if order of magnitude considerations are made on the physical feasibility of a deflection, in terms of the internal strength of the NEO. We present here the results of an investigation on the mission requirements necessary to deflect an object (or contribute to a succession of deflecting maneuvers) in terms of accessibility of the spacecraft terminal orbit from Earth with the current launchers.     

Verification of tropical cyclone deposits with oxygen isotope analyses of coeval ostracod valves

Reconstruction of prehistoric tropical cyclone (TC) activity often relies on the identification of distinctive overwash deposits (tempestites) in coastal lagoon sediments. Similar sediment deposits, however, can result from high-energy events other than TCs. In this study we assessed the utility of using the geochemistry of ostracod valves, specifically their stable oxygen isotope composition (δ¹⁸O), as a potential validation variable that could reduce the chances of misidentifying an overwash deposit as having been generated by a TC, when in fact it formed from another high-energy depositional process (type 1 error). We applied this technique to a sediment core recovered from Laguna Alejandro, Dominican Republic, which had already been analyzed for other sedimentary TC proxies. Negative δ¹⁸O anomalies identified in the ostracod valve stable isotope record are associated with TC deposits and are most easily explained by large influxes of ¹⁸O-depleted meteoric waters typical of intense tropical storms. There is potential for this technique to be used to identify TC landfalls that are not represented by overwash deposits. We, however, propose a more conservative approach and suggest this technique be used to validate the origin of a storm deposit and reduce the odds of a type 1 error.     

The Italian guidelines for seismic risk classification of constructions: technical principles and validation

The Italian “Guidelines for the seismic risk classification of constructions” approved in February 2017 define the technical principles for exploiting tax deductions with respect to seismic strengthening interventions on existing buildings (Sismabonus). Tax deductions represent a unique opportunity to improve the seismic safety of the existing Italian building stock. The guidelines are very simple and allow practitioners to deal with the sophisticated concepts behind modern seismic design, such as expected annual losses (EAL) and repair costs (expressed as a fraction of the Reconstruction Cost: %RC). The seismic risk classes of buildings and the class upgrade due to strengthening interventions can be assessed using the principles included in the guidelines. The seismic risk class is the minimum between the class defined by the building safety index at the ultimate limit state and the one related to the EAL. The latter class depends on the area under the curve of the expected losses, which is easily obtained by computing the safety index converted in the return period (annual frequency) at different limit states and the relevant %RC. This paper illustrates the technical principles at the base of the guidelines and the procedure used to calibrate the repair costs associated with the different limit states using the actual repair costs monitored in the reconstruction process following recent Italian earthquakes. Finally, simple tools to estimate the cost of the strengthening interventions to improve the seismic capacity at the life-safety limit states are provided.     

New Definition of Discovery for Solar System Objects

We propose a new formal definition of discovery for a Solar System object. It is based on an objective and mathematically rigorous algorithm to assess when a set of observations is enough to constitute a discovery. When this definition is satisfied, in almost all cases the orbit is defined well enough to establish the nature of the object discovered (Main Belt vs. Near Earth Asteroid, Trans-Neptunian vs. long period comet). The frequency of occurrence of exceptions is estimated by a set of numerical experiments. The availability of a non-subjective definition of discovery allows some rules to be adopted for the assignment of discovery credit with a minimum risk of dispute. Such rules should be fair, encourage good practice by the observers and acknowledge the contribution of the orbit computers providing the identifications and the orbits, as well as the one of all the contributing observers.