Pattern Formation and Dynamics in Rayleigh-B\'{e}nard Convection: Numerical Simulations of Experimentally Realistic Geometries

Rayleigh-B\'{e}nard convection is studied and quantitative comparisons are made, where possible, between theory and experiment by performing numerical simulations of the Boussinesq equations for a variety of experimentally realistic situations. Rectangular and cylindrical geometries of varying aspect ratios for experimental boundary conditions, including fins and spatial ramps in plate separation, are examined with particular attention paid to the role of the mean flow. A small cylindrical convection layer bounded laterally either by a rigid wall, fin, or a ramp is investigated and our results suggest that the mean flow plays an important role in the observed wavenumber. Analytical results are developed quantifying the mean flow sources, generated by amplitude gradients, and its effect on the pattern wavenumber for a large-aspect-ratio cylinder with a ramped boundary. Numerical results are found to agree well with these analytical predictions. We gain further insight into the role of mean flow in pattern dynamics by employing a novel method of quenching the mean flow numerically. Simulations of a spiral defect chaos state where the mean flow is suddenly quenched is found to remove the time dependence, increase the wavenumber and make the pattern more angular in nature.Comment: 9 pages, 10 figure

Tree-based Coarsening and Partitioning of Complex Networks

Many applications produce massive complex networks whose analysis would benefit from parallel processing. Parallel algorithms, in turn, often require a suitable network partition. For solving optimization tasks such as graph partitioning on large networks, multilevel methods are preferred in practice. Yet, complex networks pose challenges to established multilevel algorithms, in particular to their coarsening phase. One way to specify a (recursive) coarsening of a graph is to rate its edges and then contract the edges as prioritized by the rating. In this paper we (i) define weights for the edges of a network that express the edges' importance for connectivity, (ii) compute a minimum weight spanning tree $T^m$ with respect to these weights, and (iii) rate the network edges based on the conductance values of $T^m$'s fundamental cuts. To this end, we also (iv) develop the first optimal linear-time algorithm to compute the conductance values of \emph{all} fundamental cuts of a given spanning tree. We integrate the new edge rating into a leading multilevel graph partitioner and equip the latter with a new greedy postprocessing for optimizing the maximum communication volume (MCV). Experiments on bipartitioning frequently used benchmark networks show that the postprocessing already reduces MCV by 11.3%. Our new edge rating further reduces MCV by 10.3% compared to the previously best rating with the postprocessing in place for both ratings. In total, with a modest increase in running time, our new approach reduces the MCV of complex network partitions by 20.4%

Mean flow and spiral defect chaos in Rayleigh-Benard convection

We describe a numerical procedure to construct a modified velocity field that does not have any mean flow. Using this procedure, we present two results. Firstly, we show that, in the absence of mean flow, spiral defect chaos collapses to a stationary pattern comprising textures of stripes with angular bends. The quenched patterns are characterized by mean wavenumbers that approach those uniquely selected by focus-type singularities, which, in the absence of mean flow, lie at the zig-zag instability boundary. The quenched patterns also have larger correlation lengths and are comprised of rolls with less curvature. Secondly, we describe how mean flow can contribute to the commonly observed phenomenon of rolls terminating perpendicularly into lateral walls. We show that, in the absence of mean flow, rolls begin to terminate into lateral walls at an oblique angle. This obliqueness increases with Rayleigh number.Comment: 14 pages, 19 figure

DAC-h3: A Proactive Robot Cognitive Architecture to Acquire and Express Knowledge About the World and the Self

This paper introduces a cognitive architecture for a humanoid robot to engage in a proactive, mixed-initiative exploration and manipulation of its environment, where the initiative can originate from both the human and the robot. The framework, based on a biologically-grounded theory of the brain and mind, integrates a reactive interaction engine, a number of state-of-the art perceptual and motor learning algorithms, as well as planning abilities and an autobiographical memory. The architecture as a whole drives the robot behavior to solve the symbol grounding problem, acquire language capabilities, execute goal-oriented behavior, and express a verbal narrative of its own experience in the world. We validate our approach in human-robot interaction experiments with the iCub humanoid robot, showing that the proposed cognitive architecture can be applied in real time within a realistic scenario and that it can be used with naive users

The Q2Q^2-dependence of the generalised Gerasimov-Drell-Hearn integral for the deuteron, proton and neutron

The Gerasimov-Drell-Hearn (GDH) sum rule connects the anomalous contribution to the magnetic moment of the target nucleus with an energy-weighted integral of the difference of the helicity-dependent photoabsorption cross sections. The data collected by HERMES with a deuterium target are presented together with a re-analysis of previous measurements on the proton. This provides a measurement of the generalised GDH integral covering simultaneously the nucleon-resonance and the deep inelastic scattering regions. The contribution of the nucleon-resonance region is seen to decrease rapidly with increasing $Q^2$. The DIS contribution is sizeable over the full measured range, even down to the lowest measured $Q^2$. As expected, at higher $Q^2$ the data are found to be in agreement with previous measurements of the first moment of $g_1$. From data on the deuteron and proton, the GDH integral for the neutron has been derived and the proton--neutron difference evaluated. This difference is found to satisfy the fundamental Bjorken sum rule at $Q^2 = 5$ GeV$^2$.Comment: 12 pages, 10 figure

Measurement of the Proton Spin Structure Function g1p with a Pure Hydrogen Target

A measurement of the proton spin structure function g1p(x,Q^2) in deep-inelastic scattering is presented. The data were taken with the 27.6 GeV longitudinally polarised positron beam at HERA incident on a longitudinally polarised pure hydrogen gas target internal to the storage ring. The kinematic range is 0.021<x<0.85 and 0.8 GeV^2<Q^2<20 GeV^2. The integral Int_{0.021}^{0.85} g1p(x)dx evaluated at Q0^2 of 2.5 GeV^2 is 0.122+/-0.003(stat.)+/-0.010(syst.).Comment: 7 pages, 3 figures, 1 table, RevTeX late

Parallel-in-Space-and-Time Simulation of the Three-Dimensional, Unsteady Navier-Stokes Equations for Incompressible Flow

In this paper we combine the Parareal parallel-in-time method together with spatial parallelization and investigate this space-time parallel scheme by means of solving the three-dimensional incompressible Navier-Stokes equations. Parallelization of time stepping provides a new direction of parallelization and allows to employ additional cores to further speed up simulations after spatial parallelization has saturated. We report on numerical experiments performed on a Cray XE6, simulating a driven cavity flow with and without obstacles. Distributed memory parallelization is used in both space and time, featuring up to 2,048 cores in total. It is confirmed that the space-time-parallel method can provide speedup beyond the saturation of the spatial parallelization

Determination of the Deep Inelastic Contribution to the Generalised Gerasimov-Drell-Hearn Integral for the Proton and Neutron

The virtual photon absorption cross section differences [sigma_1/2-sigma_3/2] for the proton and neutron have been determined from measurements of polarised cross section asymmetries in deep inelastic scattering of 27.5 GeV longitudinally polarised positrons from polarised 1H and 3He internal gas targets. The data were collected in the region above the nucleon resonances in the kinematic range nu < 23.5 GeV and 0.8 GeV**2 < Q**2 < 12 GeV**2. For the proton the contribution to the generalised Gerasimov-Drell-Hearn integral was found to be substantial and must be included for an accurate determination of the full integral. Furthermore the data are consistent with a QCD next-to-leading order fit based on previous deep inelastic scattering data. Therefore higher twist effects do not appear significant.Comment: 6 pages, 3 figures, 1 table, revte

Observation of a Coherence Length Effect in Exclusive Rho^0 Electroproduction

Exclusive incoherent electroproduction of the rho^0(770) meson from 1H, 2H, 3He, and 14N targets has been studied by the HERMES experiment at squared four-momentum transfer Q**2>0.4 GeV**2 and positron energy loss nu from 9 to 20 GeV. The ratio of the 14N to 1H cross sections per nucleon, known as the nuclear transparency, was found to decrease with increasing coherence length of quark-antiquark fluctuations of the virtual photon. The data provide clear evidence of the interaction of the quark- antiquark fluctuations with the nuclear medium.Comment: RevTeX, 5 pages, 3 figure

Double-Spin Asymmetry in the Cross Section for Exclusive rho^0 Production in Lepton-Proton Scattering

Evidence for a positive longitudinal double-spin asymmetry = 0.24 +-0.11 (stat) +-0.02 (syst) in the cross section for exclusive diffractive rho^0(770) vector meson production in polarised lepton-proton scattering was observed by the HERMES experiment. The longitudinally polarised 27.56 GeV HERA positron beam was scattered off a longitudinally polarised pure hydrogen gas target. The average invariant mass of the photon-proton system has a value of = 4.9 GeV, while the average negative squared four-momentum of the virtual photon is = 1.7 GeV^2. The ratio of the present result to the corresponding spin asymmetry in inclusive deep-inelastic scattering is in agreement with an early theoretical prediction based on the generalised vector meson dominance model.Comment: 10 pages, 4 embedded figures, LaTe