59,232 research outputs found
Numerical Simulations of the Dark Universe: State of the Art and the Next Decade
We present a review of the current state of the art of cosmological dark
matter simulations, with particular emphasis on the implications for dark
matter detection efforts and studies of dark energy. This review is intended
both for particle physicists, who may find the cosmological simulation
literature opaque or confusing, and for astro-physicists, who may not be
familiar with the role of simulations for observational and experimental probes
of dark matter and dark energy. Our work is complementary to the contribution
by M. Baldi in this issue, which focuses on the treatment of dark energy and
cosmic acceleration in dedicated N-body simulations. Truly massive dark
matter-only simulations are being conducted on national supercomputing centers,
employing from several billion to over half a trillion particles to simulate
the formation and evolution of cosmologically representative volumes (cosmic
scale) or to zoom in on individual halos (cluster and galactic scale). These
simulations cost millions of core-hours, require tens to hundreds of terabytes
of memory, and use up to petabytes of disk storage. The field is quite
internationally diverse, with top simulations having been run in China, France,
Germany, Korea, Spain, and the USA. Predictions from such simulations touch on
almost every aspect of dark matter and dark energy studies, and we give a
comprehensive overview of this connection. We also discuss the limitations of
the cold and collisionless DM-only approach, and describe in some detail
efforts to include different particle physics as well as baryonic physics in
cosmological galaxy formation simulations, including a discussion of recent
results highlighting how the distribution of dark matter in halos may be
altered. We end with an outlook for the next decade, presenting our view of how
the field can be expected to progress. (abridged)Comment: 54 pages, 4 figures, 3 tables; invited contribution to the special
issue "The next decade in Dark Matter and Dark Energy" of the new Open Access
journal "Physics of the Dark Universe". Replaced with accepted versio
Comparison of computational codes for direct numerical simulations of turbulent Rayleigh-B\'enard convection
Computational codes for direct numerical simulations of Rayleigh-B\'enard
(RB) convection are compared in terms of computational cost and quality of the
solution. As a benchmark case, RB convection at and in a
periodic domain, in cubic and cylindrical containers is considered. A dedicated
second-order finite-difference code (AFID/RBflow) and a specialized
fourth-order finite-volume code (Goldfish) are compared with a general purpose
finite-volume approach (OpenFOAM) and a general purpose spectral-element code
(Nek5000). Reassuringly, all codes provide predictions of the average heat
transfer that converge to the same values. The computational costs, however,
are found to differ considerably. The specialized codes AFID/RBflow and
Goldfish are found to excel in efficiency, outperforming the general purpose
flow solvers Nek5000 and OpenFOAM by an order of magnitude with an error on the
Nusselt number below . However, we find that alone is not
sufficient to assess the quality of the numerical results: in fact,
instantaneous snapshots of the temperature field from a near wall region
obtained for deliberately under-resolved simulations using Nek5000 clearly
indicate inadequate flow resolution even when is converged. Overall,
dedicated special purpose codes for RB convection are found to be more
efficient than general purpose codes.Comment: 12 pages, 5 figure
A spherical shell numerical dynamo benchmark with pseudo vacuum magnetic boundary conditions
It is frequently considered that many planetary magnetic fields originate as a result of convection within planetary cores. Buoyancy forces responsible for driving the convection generate a fluid flow that is able to induce magnetic fields; numerous sophisticated computer codes are able to simulate the dynamic behaviour of such systems. This paper reports the results of a community activity aimed at comparing numerical results of several different types of computer codes that are capable of solving the equations of momentum transfer, magnetic field generation and heat transfer in the setting of a spherical shell, namely a sphere containing an inner core. The electrically conducting fluid is incompressible and rapidly rotating and the forcing of the flow is thermal convection under the Boussinesq approximation. We follow the original specifications and results reported in Harder & Hansen to construct a specific benchmark in which the boundaries of the fluid are taken to be impenetrable, non-slip and isothermal, with the added boundary condition for the magnetic field <b>B</b> that the field must be entirely radial there; this type of boundary condition for <b>B</b> is frequently referred to as ‘pseudo-vacuum’. This latter condition should be compared with the more frequently used insulating boundary condition. This benchmark is so-defined in order that computer codes based on local methods, such as finite element, finite volume or finite differences, can handle the boundary condition with ease. The defined benchmark, governed by specific choices of the Roberts, magnetic Rossby, Rayleigh and Ekman numbers, possesses a simple solution that is steady in an azimuthally drifting frame of reference, thus allowing easy comparison among results. Results from a variety of types of code are reported, including codes that are fully spectral (based on spherical harmonic expansions in angular coordinates and polynomial expansions in radius), mixed spectral and finite difference, finite volume, finite element and also a mixed Fourier-finite element code. There is good agreement among codes
- …