EvoL: The new Padova T-SPH parallel code for cosmological simulations - I. Basic code: gravity and hydrodynamics

We present EvoL, the new release of the Padova N-body code for cosmological simulations of galaxy formation and evolution. In this paper, the basic Tree + SPH code is presented and analysed, together with an overview on the software architectures. EvoL is a flexible parallel Fortran95 code, specifically designed for simulations of cosmological structure formation on cluster, galactic and sub-galactic scales. EvoL is a fully Lagrangian self-adaptive code, based on the classical Oct-tree and on the Smoothed Particle Hydrodynamics algorithm. It includes special features such as adaptive softening lengths with correcting extra-terms, and modern formulations of SPH and artificial viscosity. It is designed to be run in parallel on multiple CPUs to optimize the performance and save computational time. We describe the code in detail, and present the results of a number of standard hydrodynamical tests.Comment: 33 pages, 49 figures, accepted on A&

Angular Signatures of Annihilating Dark Matter in the Cosmic Gamma-Ray Background

The extragalactic cosmic gamma-ray background (CGB) is an interesting channel to look for signatures of dark matter annihilation. In particular, besides the imprint in the energy spectrum, peculiar anisotropy patterns are expected compared to the case of a pure astrophysical origin of the CGB. We take into account the uncertainties in the dark matter clustering properties on sub-galactic scales, deriving two possible anisotropy scenarios. A clear dark matter angular signature is achieved when the annihilation signal receives only a moderate contribution from sub-galactic clumps and/or cuspy haloes. Experimentally, if galactic foregrounds systematics are efficiently kept under control, the angular differences are detectable with the forthcoming GLAST observatory, provided that the annihilation signal contributes to the CGB for a fraction >10-20%. If, instead, sub-galactic structures have a more prominent role, the astrophysical and dark matter anisotropies become degenerate, correspondingly diluting the DM signature. As complementary observables we also introduce the cross-correlation between surveys of galaxies and the CGB and the cross-correlation between different energy bands of the CGB and we find that they provide a further sensitive tool to detect the dark matter angular signatures.Comment: 13 pages, 8 figures; improved discussion; matches published versio

Devonian and Carboniferous stratigraphical correlation and interpretation in the Central North Sea, Quadrants 25 – 44

This report details the stratigraphy and palaeogeography of Devonian and Carboniferous rocks of the UK Central North Sea for the 21CXRM Palaeozoic project. The 21CXRM Palaeozoic project results are delivered as a series of reports and digital datasets for each area. This report describes the stratigraphical correlation of Devonian and Carboniferous strata south of, and over, the Mid North Sea High (Quadrants 25-44) using both well and seismic data. It builds on the work of Cameron (1993a) and others, and uses lithostratigraphy to understand how facies vary across the area. A major outcome of this work is the detailed description of the Cleveland Group, the basinal correlative of the Scremerston, Yoredale and Millstone Grit formations within the offshore extension of the Cleveland Basin

SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch

We present a novel implementation of Smoothed Particle Hydrodynamics (SPHS) that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch -- which is second order accurate -- detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (for example, the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single-valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with `classic' SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that -- treating the particle mass similarly to the entropy -- we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov-Taylor blast wave, Kelvin-Helmholtz Instability, the `blob test', and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.Comment: 21 pages; 15 Figures. Submitted to MNRAS. Comments welcom

Transport coefficients of multi-particle collision algorithms with velocity-dependent collision rules

Detailed calculations of the transport coefficients of a recently introduced particle-based model for fluid dynamics with a non-ideal equation of state are presented. Excluded volume interactions are modeled by means of biased stochastic multiparticle collisions which depend on the local velocities and densities. Momentum and energy are exactly conserved locally. A general scheme to derive transport coefficients for such biased, velocity dependent collision rules is developed. Analytic expressions for the self-diffusion coefficient and the shear viscosity are obtained, and very good agreement is found with numerical results at small and large mean free paths. The viscosity turns out to be proportional to the square root of temperature, as in a real gas. In addition, the theoretical framework is applied to a two-component version of the model, and expressions for the viscosity and the difference in diffusion of the two species are given.Comment: 31 pages, 8 figures, accepted by J. Phys. Cond. Matte

SPH simulations of accretion disks and narrow rings

We model a massless viscous disk using Smoothed Particle Hydrodynamics (SPH) and note that it evolves according to the Lynden-Bell \& Pringle theory (1974) until a non-axisymmetric instability develops at the inner edge of the disk. This instability may have the same origin as the instability of initially axisymmetric viscous disks discussed by Lyubarskij et al. (1994). To clarify the evolution we evolved single and double rings of particles. It is actually inconsistent with the SPH scheme to set up a single ring as an initial condition because SPH assumes a smoothed initial state. As would be expected from an SPH simulation, the ring rapidly breaks up into a band. We analyse the stability of the ring and show that the predictions are confirmed by the simulation

Enabling internal electronic circuitry within additively manufactured metal structures - The effect and importance of inter-laminar topography

Purpose: This paper aims to explore the potential of ultrasonic additive manufacturing (UAM) to incorporate the direct printing of electrical materials and arrangements (conductors and insulators) at the interlaminar interface of parts during manufacture to allow the integration of functional and optimal electrical circuitries inside dense metallic objects without detrimental effect on the overall mechanical integrity. This holds promise to release transformative device functionality and applications of smart metallic devices and products. Design/methodology/approach: To ensure the proper electrical insulation between the printed conductors and metal matrices, an insulation layer with sufficient thickness is required to accommodate the rough interlaminar surface which is inherent to the UAM process. This in turn increases the total thickness of printed circuitries and thereby adversely affects the integrity of the UAM part. A specific solution is proposed to optimise the rough interlaminar surface through deforming the UAM substrates via sonotrode rolling or UAM processing. Findings: The surface roughness (Sa) could be reduced from 4.5 to 4.1 µm by sonotrode rolling and from 4.5 to 0.8 µm by ultrasonic deformation. Peel testing demonstrated that sonotrode-rolled substrates could maintain their mechanical strength, while the performance of UAM-deformed substrates degraded under same welding conditions ( approximately 12 per cent reduction compared with undeformed substrates). This was attributed to the work hardening of deformation process which was identified via dual-beam focussed ion beam–scanning electron microscope investigation. Originality/value: The sonotrode rolling was identified as a viable methodology in allowing printed electrical circuitries in UAM. It enabled a decrease in the thickness of printed electrical circuitries by ca. 25 per cent

Magnetic field structure due to the global velocity field in spiral galaxies

We present a set of global, self-consistent N-body/SPH simulations of the dynamic evolution of galactic discs with gas and including magnetic fields. We have implemented a description to follow the evolution of magnetic fields with the ideal induction equation in the SPH part of the Vine code. Results from a direct implementation of the field equations are compared to a representation by Euler potentials, which pose a div(B)-free description, an constraint not fulfilled for the direct implementation. All simulations are compared to an implementation of magnetic fields in the Gadget code which includes also cleaning methods for div(B). Starting with a homogeneous seed field we find that by differential rotation and spiral structure formation of the disc the field is amplified by one order of magnitude within five rotation periods of the disc. The amplification is stronger for higher numerical resolution. Moreover, we find a tight connection of the magnetic field structure to the density pattern of the galaxy in our simulations, with the magnetic field lines being aligned with the developing spiral pattern of the gas. Our simulations clearly show the importance of non-axisymmetry for the evolution of the magnetic field.Comment: 17 pages, 18 figure

Flow-History-Dependent Behavior in Entangled Polymer Melt Flow with Multiscale Simulation

Polymer melts represent the flow-history-dependent behavior. To clearly show this behavior, we have investigated flow behavior of an entangled polymer melt around two cylinders placed in tandem along the flow direction in a two dimensional periodic system. In this system, the polymer states around a cylinder in downstream side are different from the ones around another cylinder in upstream side because the former ones have a memory of a strain experienced when passing around the cylinder in upstream side but the latter ones do not have the memory. Therefore, the shear stress distributions around two cylinders are found to be different from each other. Moreover, we have found that the averaged flow velocity decreases accordingly with increasing the distance between two cylinders while the applied external force is constant. While this behavior is consistent with that of the Newtonian fluid, the flow-history-dependent behavior enhances the reduction of the flow resistance.Comment: 6 pages, 3 figures, Proceedings of 5th International Mini-Symposium on Liquid