An Unsplit, Cell-Centered Godunov Method for Ideal MHD

We present a second-order Godunov algorithm for multidimensional, ideal MHD. Our algorithm is based on the unsplit formulation of Colella (J. Comput. Phys. vol. 87, 1990), with all of the primary dependent variables centered at the same location. To properly represent the divergence-free condition of the magnetic fields, we apply a discrete projection to the intermediate values of the field at cell faces, and apply a filter to the primary dependent variables at the end of each time step. We test the method against a suite of linear and nonlinear tests to ascertain accuracy and stability of the scheme under a variety of conditions. The test suite includes rotated planar linear waves, MHD shock tube problems, low-beta flux tubes, and a magnetized rotor problem. For all of these cases, we observe that the algorithm is second-order accurate for smooth solutions, converges to the correct weak solution for problems involving shocks, and exhibits no evidence of instability or loss of accuracy due to the possible presence of non-solenoidal fields.Comment: 37 Pages, 9 Figures, submitted to Journal of Computational Physic

Porous squeeze-film flow

The squeeze-film flow of a thin layer of Newtonian fluid filling the gap between a flat impermeable surface moving under a prescribed constant load and a flat thin porous bed coating a stationary flat impermeable surface is considered. Unlike in the classical case of an impermeable bed, in which an infinite time is required for the two surfaces to touch, for a porous bed contact occurs in a finite contact time. Using a lubrication approximation an implicit expression for the fluid layer thickness and an explicit expression for the contact time are obtained and analysed. In addition, the fluid particle paths are calculated, and the penetration depths of fluid particles into the porous bed are determined. In particular, the behaviour in the asymptotic limit of small permeability, in which the contact time is large but finite, is investigated. Finally, the results are interpreted in the context of lubrication in the human knee joint, and some conclusions are drawn about the contact time of the cartilage-coated femoral condyles and tibial plateau and the penetration of nutrients into the cartilage

Squeeze-Film Flow in the Presence of a Thin Porous Bed, with Application to the Human Knee Joint

Motivated by the desire for a better understanding of the lubrication of the human knee joint, the squeeze-film flow of a thin layer of Newtonian fluid (representing the synovial fluid) filling the gap between a flat impermeable surface (representing the femoral condyles) and a flat thin porous bed (representing the articular cartilage) coating a stationary flat impermeable surface (representing the tibial plateau) is considered. As the impermeable surface approaches the porous bed under a prescribed constant load all of the fluid is squeezed out of the gap in a finite contact time. In the context of the knee, the size of this contact time suggests that when a person stands still for a short period of time their knees may be fluid lubricated, but that when they stand still for a longer period of time contact between the cartilage-coated surfaces may occur. The fluid particle paths are calculated, and the penetration depths of fluid particles into the porous bed are determined. In the context of the knee, these penetration depths provide a measure of how far into the cartilage nutrients are carried by the synovial fluid, and suggest that when a person stands still nutrients initially in the fluid layer penetrate only a relatively small distance into the cartilage. However, the model also suggests that the cumulative effect of repeated loading and unloading of the knees during physical activity such as walking or running may be sufficient to carry nutrients deep into the cartilage

‘Generation rent’ and the ability to ‘settle down’: economic and geographical variation in young people’s housing transitions

© 2016 Informa UK Limited, trading as Taylor & Francis GroupThe term ‘Generation Rent’ denotes young people who are increasingly living in the private rented sector for longer periods of their lives because they are unable to access homeownership or social housing. Drawing on qualitative data from two studies with young people and key-actors, this paper considers the phenomenon of ‘Generation Rent’ from the perspective of youth transitions and the concept of ‘home’. These frameworks posit that young people leaving the parental home traverse housing and labour markets until they reach a point of ‘settling down’. However, our data indicate that many young people face difficulties in this ‘settling’ process as they have to contend with insecure housing, unstable employment and welfare cuts which often force them to be flexible and mobile. This leaves many feeling frustrated as they struggle to remain fixed in place in order to ‘settle down’ and benefit from the positive qualities of home. Taking a Scottish focus, this paper further highlights the geographical dimension to these challenges and argues that those living in expensive and/or rural areas may find it particularly difficult to settle down

Post-Foucauldian governmentality: what does it offer critical social policy analysis?

This article considers the theoretical perspective of post-Foucauldian governmentality, especially the insights and challenges it poses for applied researchers within the critical social policy tradition. The article firstly examines the analytical strengths of this approach to understanding power and rule in contemporary society, before moving on to consider its limitations for social policy. It concludes by arguing that these insights can be retained, and some of the weaknesses overcome, by adopting a ‘realist governmentality’ approach (Stenson 2005, 2008). This advocates combining traditional discursive analysis with more ethnographic methods in order to render visible the concrete activity of governing, and unravel the messiness, complexity and unintended consequences involved in the struggles around subjectivity

Some properties of three related viruses: Andean potato latent, dulcamara mottle, and Ononis yellow mosaic

RESP-544

The beamformer and correlator for the Large European Array for Pulsars

The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.Comment: 10 pages, 6 figures, accepted for publication in Astronomy and Computin

Squeeze-film flow between a curved impermeable bearing and a flat porous bed

Axisymmetric squeeze-film flow in the thin gap between a stationary flat thin porous bed and a curved impermeable bearing moving under a prescribed constant load is analysed. The unsteady Reynolds equation is formulated and solved for the fluid pressure. This solution is used to obtain the time for the minimum fluid layer thickness to reduce to a given value, and, in particular, the finite time for the bearing and the bed to come into contact. The effect of varying the shape of the bearing and the permeability of the layer is investigated, and, in particular, it is found that both the contact time and the fluid pressure behave qualitatively differently for beds with small and large permeabilities. In addition, the paths of fluid particles initially situated in both the fluid layer and the porous bed are calculated. In particular, it is shown that, unlike in the case of a flat bearing, for a curved bearing there are fluid particles, initially situated in the fluid layer, that flow from the fluid layer into the porous bed and then re-emerge into the fluid layer, and the region in which these fluid particles are initially situated is determined

Radiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores

We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment, and whether they form a small or large number of protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from ~0.1 pc scales to ~10 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.Comment: 23 pages, 18 figures, emulateapj format. Accepted to ApJ. This version has minor typo fixes and small additions, no significant changes. Resolution of images severely degraded to fit within size limit. Download the full paper from http://www.astro.princeton.edu/~krumholz/recent.htm