Resolving high Reynolds numbers in SPH simulations of subsonic turbulence

Accounting for the Reynolds number is critical in numerical simulations of turbulence, particularly for subsonic flow. For Smoothed Particle Hydrodynamics (SPH) with constant artificial viscosity coefficient alpha, it is shown that the effective Reynolds number in the absence of explicit physical viscosity terms scales linearly with the Mach number - compared to mesh schemes, where the effective Reynolds number is largely independent of the flow velocity. As a result, SPH simulations with alpha=1 will have low Reynolds numbers in the subsonic regime compared to mesh codes, which may be insufficient to resolve turbulent flow. This explains the failure of Bauer and Springel (2011, arXiv:1109.4413v1) to find agreement between the moving-mesh code AREPO and the GADGET SPH code on simulations of driven, subsonic (v ~ 0.3 c_s) turbulence appropriate to the intergalactic/intracluster medium, where it was alleged that SPH is somehow fundamentally incapable of producing a Kolmogorov-like turbulent cascade. We show that turbulent flow with a Kolmogorov spectrum can be easily recovered by employing standard methods for reducing alpha away from shocks.Comment: 5 pages, accepted to MNRAS Letters. Movies showing SPH doing the things that SPH is fundamentally incapable of doing at http://users.monash.edu.au/~dprice/pubs/reynolds/index.html#movies v2: minor changes to reflect proof-corrected version, couple of refs adde

Design, fabrication and evaluation of tungsten encapsulated silver infiltrated porous tungsten rocket nozzle inserts Final report

Rocket nozzle throat insert encapsulated in solid tungsten cas

Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics

We present a constrained formulation of Dedner et al's hyperbolic/parabolic divergence cleaning scheme for enforcing the \nabla\dot B = 0 constraint in Smoothed Particle Magnetohydrodynamics (SPMHD) simulations. The constraint we impose is that energy removed must either be conserved or dissipated, such that the scheme is guaranteed to decrease the overall magnetic energy. This is shown to require use of conjugate numerical operators for evaluating \nabla\dot B and \nabla{\psi} in the SPMHD cleaning equations. The resulting scheme is shown to be stable at density jumps and free boundaries, in contrast to an earlier implementation by Price & Monaghan (2005). Optimal values of the damping parameter are found to be {\sigma} = 0.2-0.3 in 2D and {\sigma} = 0.8-1.2 in 3D. With these parameters, our constrained Hamiltonian formulation is found to provide an effective means of enforcing the divergence constraint in SPMHD, typically maintaining average values of h |\nabla\dot B| / |B| to 0.1-1%, up to an order of magnitude better than artificial resistivity without the associated dissipation in the physical field. Furthermore, when applied to realistic, 3D simulations we find an improvement of up to two orders of magnitude in momentum conservation with a corresponding improvement in numerical stability at essentially zero additional computational expense.Comment: 28 pages, 25 figures, accepted to J. Comput. Phys. Movies at http://www.youtube.com/playlist?list=PL215D649FD0BDA466 v2: fixed inverted figs 1,4,6, and several color bar

Open plan and academe: pre- and post-hoc conversations

There now exists a strong body of evidence that creative workplaces can, in certain circumstances, exert beneficial influences on organisational cultures and outputs. Academia tends to resist such spaces and faculty buildings. The reasons are explored but the reactions of staff are not found to be different from those reported in the literature on general creative spaces. The success or failure of team oriented workspaces is in large part a socially constructed perception influenced by the manner of implementation and management. As elsewhere new workplaces are about new conversations. The cases studied lead to a model of the tensions inherent in workplace redesign.</p

Dust and gas mixtures with multiple grain species - a one-fluid approach

GL acknowledges funding from the European Research Council for the FP7 ERC advanced grant project ECOGAL. DJP is very grateful for funding via an ARC Future Fellowship, FT130100034, and Discovery Project grant DP130102078.We derive the single-fluid evolution equations describing a mixture made of a gas phase and an arbitrary number of dust phases, generalizing the approach developed by Laibe & Price. A generalization for continuous dust distributions as well as analytic approximations for strong drag regimes is also provided. This formalism lays the foundation for numerical simulations of dust populations in a wide range of astrophysical systems while avoiding limitations associated with a multiple-fluid treatment. The usefulness of the formalism is illustrated on a series of analytical problems, namely the DUSTYBOX, DUSTYSHOCK and DUSTYWAVE problems as well as the radial drift of grains and the streaming instability in protoplanetary discs. We find physical effects specific to the presence of several dust phases and multiple drag time-scales, including non-monotonic evolution of the differential velocity between phases and increased efficiency of the linear growth of the streaming instability. Interestingly, it is found that under certain conditions, large grains can migrate outwards in protoplanetary discs. This may explain the presence of small pebbles at several hundreds of astronomical units from their central star.Publisher PDFPeer reviewe