Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.Comment: 21 pages, 7 figures, to appear in Proceedings of the SPPEXA symposium, Lecture Notes in Computational Science and Engineering (LNCSE), Springe

Ab initio methods for finite temperature two-dimensional Bose gases

The stochastic Gross-Pitaevskii equation and modified Popov theory are shown to provide an ab initio description of finite temperature, weakly-interacting two-dimensional Bose gas experiments. Using modified Popov theory, a systematic approach is developed in which the momentum cut-off inherent to classical field methods is removed as a free parameter. This is shown to yield excellent agreement with the recent experiment of Hung et al. [Nature, 470, 236 (2011)], verifying that the stochastic Gross-Pitaevskii equation captures the observed universality and scale-invariance.Comment: 5 pages, 4 figure

Phase coherence in quasicondensate experiments: an ab initio analysis via the stochastic Gross-Pitaevskii equation

We perform an ab initio analysis of the temperature dependence of the phase coherence length of finite temperature, quasi-one-dimensional Bose gases measured in the experiments of Richard et al. (Phys. Rev. Lett. 91, 010405 (2003)) and Hugbart et al. (Eur. Phys. J. D 35, 155-163 (2005)), finding very good agreement across the entire observed temperature range ($0.8<T/T_{\phi}<28$). Our analysis is based on the one-dimensional stochastic Gross-Pitaevskii equation, modified to self-consistently account for transverse, quasi-one-dimensional effects, thus making it a valid model in the regime $\mu ~ few \hbar \omega_\perp$. We also numerically implement an alternative identification of $T_{\phi}$, based on direct analysis of the distribution of phases in a stochastic treatment.Comment: Amended manuscript with improved agreement to experiment, following some additional clarifications by Mathilde Hugbart and Fabrice Gerbier and useful comments by the reviewer; accepted for publication in Physical Review

Historical roots of Agile methods: where did “Agile thinking” come from?

The appearance of Agile methods has been the most noticeable change to software process thinking in the last fifteen years [16], but in fact many of the “Agile ideas” have been around since 70’s or even before. Many studies and reviews have been conducted about Agile methods which ascribe their emergence as a reaction against traditional methods. In this paper, we argue that although Agile methods are new as a whole, they have strong roots in the history of software engineering. In addition to the iterative and incremental approaches that have been in use since 1957 [21], people who criticised the traditional methods suggested alternative approaches which were actually Agile ideas such as the response to change, customer involvement, and working software over documentation. The authors of this paper believe that education about the history of Agile thinking will help to develop better understanding as well as promoting the use of Agile methods. We therefore present and discuss the reasons behind the development and introduction of Agile methods, as a reaction to traditional methods, as a result of people's experience, and in particular focusing on reusing ideas from histor

Matter-wave dark solitons: stochastic vs. analytical results

The dynamics of dark matter-wave solitons in elongated atomic condensates are discussed at finite temperatures. Simulations with the stochastic Gross-Pitaevskii equation reveal a noticeable, experimentally observable spread in individual soliton trajectories, attributed to inherent fluctuations in both phase and density of the underlying medium. Averaging over a number of such trajectories (as done in experiments) washes out such background fluctuations, revealing a well-defined temperature-dependent temporal growth in the oscillation amplitude. The average soliton dynamics is well captured by the simpler dissipative Gross-Pitaevskii equation, both numerically and via an analytically-derived equation for the soliton center based on perturbation theory for dark solitons.Comment: 4 pages, 3 figures. Added several reference

Quantitative study of quasi-one-dimensional Bose gas experiments via the stochastic Gross-Pitaevskii equation

The stochastic Gross-Pitaevskii equation is shown to be an excellent model for quasi-one-dimensional Bose gas experiments, accurately reproducing the in situ density profiles recently obtained in the experiments of Trebbia et al. [Phys. Rev. Lett. 97, 250403 (2006)] and van Amerongen et al. [Phys. Rev. Lett. 100, 090402 (2008)], and the density fluctuation data reported by Armijo et al. [Phys. Rev. Lett. 105, 230402 (2010)]. To facilitate such agreement, we propose and implement a quasi-one-dimensional stochastic equation for the low-energy, axial modes, while atoms in excited transverse modes are treated as independent ideal Bose gases.Comment: 10 pages, 5 figures; updated figures with experimental dat

Assessing the usefulness of acute physiological responses following resistance exercise: sensitivity, magnitude of change and time course of measures

A variety of strategies exist to modulate acute physiological responses following resistance exercise aimed at enhancing recovery and/or adaptation processes. To assess the true impact of these strategies, it is important to know the ability of measures to detect meaningful change. We investigated the sensitivity of measures used to quantify acute physiological responses to resistance exercise and constructed a physiological profile to characterise the magnitude of change and time course of this response. Eight males, accustomed to regular resistance exercise, performed experimental sessions during a ‘control week’, void of an exercise stimulus. Participants repeated this sequence of experimental sessions the following week, termed the ‘exercise week’, except they performed a bout of lower-limb resistance exercise following baseline assessments. Assessments were conducted at baseline, 2, 6, 24, 48, 72 and 96 h post-intervention. Based on the signal-to-noise ratio, the most sensitive measures were maximal voluntary isometric contraction, 20m sprint, countermovement jump peak force, rate of force development (100-200ms), muscle soreness, daily analysis of life demands for athletes Part B, limb girth, matrix metalloproteinase-9, interleukin-6, creatine kinase and high sensitivity C-reactive protein with ratios of >1.5. There were clear changes in these measures following resistance exercise, determined via magnitude-based inferences. These findings highlight measures that can detect real changes in acute physiological responses following resistance exercise in trained individuals. Researchers investigating strategies to manipulate acute physiological responses for recovery and/or adaptation can use these measures, as well as recommended sampling points, to be confident that their interventions are making a worthwhile impact

Fluctuating and dissipative dynamics of dark solitons in quasi-condensates

The fluctuating and dissipative dynamics of matter-wave dark solitons within harmonically trapped, partially condensed Bose gases is studied both numerically and analytically. A study of the stochastic Gross-Pitaevskii equation, which correctly accounts for density and phase fluctuations at finite temperatures, reveals dark soliton decay times to be lognormally distributed at each temperature, thereby characterizing the previously predicted long lived soliton trajectories within each ensemble of numerical realizations (S.P. Cockburn {\it et al.}, Phys. Rev. Lett. 104, 174101 (2010)). Expectation values for the average soliton lifetimes extracted from these distributions are found to agree well with both numerical and analytic predictions based upon the dissipative Gross-Pitaevskii model (with the same {\it ab initio} damping). Probing the regime for which $0.8 k_{B}T < \mu < 1.6 k_{B}T$, we find average soliton lifetimes to scale with temperature as $\tau\sim T^{-4}$, in agreement with predictions previously made for the low-temperature regime $k_{B}T\ll\mu$. The model is also shown to capture the experimentally-relevant decrease in the visibility of an oscillating soliton due to the presence of background fluctuations.Comment: 17 pages, 14 figure

Discontinuous Galerkin method for multifluid Euler equations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106454/1/AIAA2013-2595.pd