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.     

Reduced complexity strategies for modelling urban floodplain inundation

Significant advances in flood inundation modelling have been made in the last decade through the use of a new generation of 2D hydraulic numerical models. These offer the potential to predict the local pattern and timing of flood depth and velocity, enabling informed flood risk zoning and improved emergency planning. With the availability of high resolution DEMs derived from airborne lidar, these models can theoretically now be routinely parameterized to represent considerable topographic complexity, even in urban areas where the potential exists to represent flows at the scale of individual buildings. Currently, however, computational constraints on conventional finite element and volume codes typically require model discretization at scales well below those achievable with lidar and are thus unable to make optimal use of this emerging data stream.In this paper we review two strategies that attempt to address this mismatch between model and data resolution in an effort to improve urban flood forecasts. The first of these strives for a solution by simplifying the mathematical formulation of the numerical model by using a computationally efficient 2D raster storage cell approach coupled to a 1D channel model. This parsimonious model structure enables simulations over large model domains offering the opportunity to employ a topographic discretization strategy which explicitly represents the built environment. The second approach seeks to further reduce the computational overhead of this raster method by employing a subgrid parameterization to represent the effect of buildings and micro-relief on flow pathways and floodplain storage. This multi-scale methodology enables highly efficient model applications at coarse spatial resolutions while retaining information about the complex geometry of the built environment.These two strategies are evaluated through numerical experiments designed to reconstruct a flood in the small town of Linton in southern England, which occurred in response to a 1 in 250 year rainfall event in October 2001. Results from both approaches are encouraging, with the spatial pattern of inundation and flood wave propagation matching observations well. Both show significant advantages over a coarse resolution model without subgrid parameterisation, particularly in terms of their ability to reproduce both hydrograph and inundation depth measurements simultaneously, without need for recalibration. The subgrid parameterization is shown to achieve this without contributing significant computational complexity and reduces model run-times by an order of magnitude.     

An inexpensive system for continuous lake core photography

Representative images of split sediment cores document geological records, serve as visual aids for multimedia, and provide samples for image analysis. Existing core imaging systems, often containing integrated hardware and software, are capable of providing excellent images; unfortunately these systems are also bulky and expensive, and are therefore limited to labs that process large volumes of core and those actively involved in image analysis. For most research groups such systems are unfeasible to buy and maintain. Producing good quality core images with a hand-held digital camera is very difficult without a consistent lighting source, and changes in camera angle (in all three spatial dimensions) between photographs may prevent accurate image compositing. Tripods allow for camera stability, but typically do not accommodate downward-facing photography. Presented here is an easily constructed core-imaging system that minimizes many of the drawbacks of personal digital cameras. Software necessary for using this system is readily available and can be run from a personal computer.     

Neutron background in large-scale xenon detectors for dark matter searches

Simulations of the neutron background for future large-scale particle dark matter detectors are presented. Neutrons were generated in rock and detector elements via spontaneous fission and (α,n) reactions, and by cosmic-ray muons. The simulation techniques and results are discussed in the context of the expected sensitivity of a generic liquid xenon dark matter detector. Methods of neutron background suppression are investigated. A sensitivity of 10⁻⁹–10⁻¹⁰ pb to WIMP-nucleon interactions can be achieved by a tonne-scale detector.     

Connectivity clues from short-term variability in settlement and geochemical tags of mytilid mussels

The use of geochemical tags in calcified structures of fish and invertebrates is an exciting tool for investigating larval population connectivity. Tag evaluation over relatively short intervals (weeks) may detect environmental and ecological variability at a temporal scale highly relevant to larval transport and settlement. We collected newly settled mussels (Mytilus californianus and M. galloprovincialis) weekly during winter/spring of 2002 along the coast of San Diego, CA, USA, at sites on the exposed coast (SIO) and in a protected coastal bay (HI), to investigate temporal patterns of geochemical tags in mussel shells. Analyses of post-settlement shell via LA-ICP-MS revealed statistically significant temporal variability for all elements we examined (Mg, Mn, Cu, Sr, Cd, Ba, Pb and U). Despite this, our ability to distinguish multielemental signatures between sites was largely conserved. Throughout our 13-week study, SIO and HI mussels could be chemically distinguished from one another in 78–87% of all cases. Settlement varied between 2 and 27 settlers gram-byssus^?1 week^?1 at SIO and HI, and both sites were characterized by 2–3 weeks with “high” settlement. Geochemical tags recorded in early larval shell of newly settled mussels differed between “high” and “low” settlement weeks at both sites (MANOVA), driven by Mg and Sr at SIO (p = 0.013) and Sr, Cd, Ba and Pb at HI (p < 0.001). These data imply that shifts in larval sources or transport corridors were responsible for observed settlement variation, rather than increased larval production. In particular, increased settlement at HI was observed concurrent with the appearance of geochemical tags (e.g., elevated Cd), suggesting that those larvae were retained in upwelled water near the mouth of the bay. Such shifts may reflect short-term changes in connectivity among sites due to altered transport corridors, and influence the demography of local populations.     

Short-term survival of severe propeller strike injuries and observations on wound progression in a bottlenose dolphin

Until recently, vessel collisions with small cetaceans were presumed rare, mainly as a consequence of limited reporting. Observations on dolphin wound healing from propeller strike injuries also remain scarce. We present an extreme case of a bottlenose dolphin (Tursiops truncatus) with multiple propeller wounds, including a penetration to the bone, where survival was possible for at least 23 days post injury. We used photographic records in conjunction with field observations to describe wound progression in the absence of treatment. Considering the severity of the wounds, it was surprising that the injuries were not immediately fatal. A practical solution remains to be found for the problem of odontocete vessel collisions, in particular for small highly mobile species using neritic waters. In view of the ongoing problem in the Hauraki Gulf and, indeed, worldwide, we recommend continued evaluation of cases and reporting of incidents, in addition to investigations into mechanisms that may reduce the risk of occurrence of vessel collisions with cetaceans.     

Study and suppression of anomalous fast events in inorganic scintillators for dark matter searches

The status of dark matter searches with inorganic scintillator detectors at Boulby mine is reviewed and the results of tests with a CsI(Tl) crystal are presented. The objectives of the latter experiment were to study anomalous fast events previously observed and to identify ways to remove this background. Clear indications were found that these events were due to surface contamination of crystals by alphas, probably from radon decay. A new array of unencapsulated NaI(Tl) crystals immersed either in liquid paraffin or pure nitrogen atmosphere is under construction at Boulby. Such an approach allows complete control of the surface of the crystals and the ability to remove any surface contamination. First data from the unencapsulated NaI(Tl) do not show the presence of anomalous fast events.