Fast iterative solvers for geomechanics in a commercial FE code

There is a pressing need to improve the feasibility of three-dimensional finite element (FE) methods applied to many problems in civil engineering. This is particularly the case for static analyses in geotechnical engineering: ideally, models would be 3D, follow the actual geometry, use non-linear material formulations and allow simulation of construction sequences, and all of this with a reasonable degree of accuracy. One major obstacle to improvements in this regard is the difficulty in solving of the set of (linearised) algebraic equations which arises from a typical discretisation approach. Very large systems become cumbersome for direct techniques to solve economically. This paper describes the incorporation of iterative (rather than direct) solution techniques, developed through University research, into commercial FE software for geotechnics

Peskun–Tierney ordering for Markovian Monte Carlo: Beyond the reversible scenario

Historically time-reversibility of the transitions or processes underpinning Markov chain Monte Carlo methods (MCMC) has played a key role in their development, while the self-adjointness of associated operators together with the use of classical functional analysis techniques on Hilbert spaces have led to powerful and practically successful tools to characterise and compare their performance. Similar results for algorithms relying on nonreversible Markov processes are scarce. We show that for a type of nonreversible Monte Carlo Markov chains and processes, of current or renewed interest in the physics and statistical literatures, it is possible to develop comparison results which closely mirror those available in the reversible scenario. We show that these results shed light on earlier literature, proving some conjectures and strengthening some earlier results

Post-outburst X-ray flux and timing evolution of Swift J1822.3-1606

Swift J1822.3-1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following the detection of several bursts. The source was found to have a period of 8.4377 s and was identified as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 0.1185154343(8) s$^{-1}$ and a frequency derivative of $-4.3\pm0.3\times10^{-15}$ at MJD 55761.0, in a timing analysis that include significant non-zero second and third frequency derivatives that we attribute to timing noise. This corresponds to an estimated spin-down inferred dipole magnetic field of $B\sim5\times10^{13}$ G, consistent with previous estimates though still possibly affected by unmodelled noise. We find that the post-outburst 1--10 keV flux evolution can be characterized by a double-exponential decay with decay timescales of $15.5\pm0.5$ and $177\pm14$ days. We also fit the light curve with a crustal cooling model which suggests that the cooling results from heat injection into the outer crust. We find that the hardness-flux correlation observed in magnetar outbursts also characterizes the outburst of Swift J1822.3-1606. We compare the properties of Swift J1822.3-1606 with those of other magnetars and their outbursts.Comment: 12 pages, 8 figures, accepted for publication in Ap

Submarine glacial-landform distribution along an Antarctic Peninsula palaeo-ice stream: A shelf-slope transect through the Marguerite Trough system (66-70° S)

The Antarctic Peninsula comprises a thin spine of mountains and islands presently covered by an ice sheet up to 500 m thick that drains eastwards and westwards via outlet glaciers (Davies et al. 2012). Recently, the Peninsula has undergone rapid warming, resulting in the collapse of fringing ice shelves and the retreat, thinning and acceleration of marine-terminating outlet glaciers (e.g. Pritchard & Vaughan 2007). At the Last Glacial Maximum (LGM), the ice sheet expanded to the continental-shelf break around the Peninsula, and was organised into a series of ice streams that drained along cross-shelf bathymetric troughs (Ó Cofaigh et al. 2014). Marguerite Bay is located on the west side of the Antarctic Peninsula, at about 66° to 70° S (Fig. 1). A 12–80 km wide and 370 km long trough extends across the bay from the northern terminus of George VI Ice Shelf to the continental shelf edge. Extensive marine-geophysical surveys of the trough reveal a suite of glacial landforms which record past flow of an ice stream, which extended to the shelf edge at, or shortly after, the LGM. Subsequent retreat of the ice stream was underway by ~14 kyr ago and proceeded rapidly to the mid-shelf, where it slowed before accelerating once again to the inner shelf at ~9 kyr (Kilfeather et al. 2011).This is the author accepted manuscript. The final version is available from the Geological Society of London via https://doi.org/10.1144/M46.18

Red Noise in Anomalous X-ray Pulsar Timing Residuals

Anomalous X-ray Pulsars (AXPs), thought to be magnetars, exhibit poorly understood deviations from a simple spin-down called "timing noise". AXP timing noise has strong low-frequency components which pose significant challenges for quantification. We describe a procedure for extracting two quantities of interest, the intensity and power spectral index of timing noise. We apply this procedure to timing data from three sources: a monitoring campaign of five AXPs, observations of five young pulsars, and the stable rotator PSR B1937+21.Comment: submitted to the proceedings of the "40 Years of Pulsars" conferenc

A Survey of 56 Mid-latitude EGRET Error Boxes for Radio Pulsars

We have conducted a radio pulsar survey of 56 unidentified gamma-ray sources from the 3rd EGRET catalog which are at intermediate Galactic latitudes (5 deg. < |b| < 73 deg.). For each source, four interleaved 35-minute pointings were made with the 13-beam, 1400-MHz multibeam receiver on the Parkes 64-m radio telescope. This covered the 95% error box of each source at a limiting sensitivity of about 0.2 mJy to pulsed radio emission for periods P > 10 ms and dispersion measures < 50 pc cm-3. Roughly half of the unidentified gamma-ray sources at |b| > 5 deg. with no proposed active galactic nucleus counterpart were covered in this survey. We detected nine isolated pulsars and four recycled binary pulsars, with three from each class being new. Timing observations suggest that only one of the pulsars has a spin-down luminosity which is even marginally consistent with the inferred luminosity of its coincident EGRET source. Our results suggest that population models, which include the Gould belt as a component, overestimate the number of isolated pulsars among the mid-latitude Galactic gamma-ray sources and that it is unlikely that Gould belt pulsars make up the majority of these sources. However, the possibility of steep pulsar radio spectra and the confusion of terrestrial radio interference with long-period pulsars (P > 200 ms) having very low dispersion measures (< 10 pc cm-3, expected for sources at a distance of less than about 1 kpc) prevent us from strongly ruling out this hypothesis. Our results also do not support the hypothesis that millisecond pulsars make up the majority of these sources. Non-pulsar source classes should therefore be further investigated as possible counterparts to the unidentified EGRET sources at intermediate Galactic latitudes.Comment: 24 pages, including 4 figures and 3 tables. Accepted for publication in Ap

Are Japanese and European gastric cancer the same biological entity? An immunohistochemical study.

To examine the suggested biological difference between Japanese and British gastric cancers, immunohistochemistry was used to demonstrate eight markers of biological activity in a matched series of 40 Japanese and 33 British cases. There were no differences in the proportions of Japanese and British tumours positive to epidermal growth factor, epidermal growth factor receptor, transforming growth factor alpha, cripto or p53. A significantly greater proportion of British tumours were positive to c-erbB-2 whilst a significantly greater proportion of Japanese tumours were positive to nm23. British tumours had a significantly greater mean proliferating cell nuclear antigen proliferation index than Japanese tumours. These differences could be clinically significant

No detectable radio emission from the magnetar-like pulsar in Kes 75

The rotation-powered pulsar PSR J1846-0258 in the supernova remnant Kes 75 was recently shown to have exhibited magnetar-like X-ray bursts in mid-2006. Radio emission has not yet been observed from this source, but other magnetar-like sources have exhibited transient radio emission following X-ray bursts. We report on a deep 1.9 GHz radio observation of PSR J1846-0258 with the 100-m Green Bank Telescope in late 2007 designed to search for radio pulsations or bursts from this target. We have also analyzed three shorter serendipitous 1.4 GHz radio observations of the source taken with the 64-m Parkes telescope during the 2006 bursting period. We detected no radio emission from PSR J1846-0258 in either the Green Bank or Parkes datasets. We place an upper limit of 4.9 \mu Jy on coherent pulsed emission from PSR J1846-0258 based on the 2007 November 2 observation, and an upper limit of 27 \mu Jy around the time of the X-ray bursts. Serendipitously, we observed radio pulses from the nearby RRAT J1846-02, and place a 3\sigma confidence level upper limit on its period derivative of 1.7 * 10^{-13}, implying its surface dipole magnetic field is less than 2.6 * 10^{13} G.Comment: 15 pages, 2 figures, submitted to Ap

An automated method for mapping geomorphological expressions of former subglacial meltwater pathways (hummock corridors) from high resolution digital elevation data

Elongated tracts of hummocks or ‘hummock corridors’, exposed on palaeo-ice sheet beds, are believed to represent former subglacial meltwater pathways. Here, we present a method, coded in MATLAB, for automatically detecting and mapping hummock corridors from high-resolution digital elevation models (DEMs). Initially the DEM is filtered to remove bed roughness outside the size range of hummocks. A Fast Fourier Transform is then performed to determine the dominant orientation of hummock corridors and remove misaligned features. Finally, image segmentation is used to isolate and extract the hummock corridors as a binary mask. We tested this automated approach visually and statistically against detailed manual mapping in three areas of Canada and northern Scandinavia. Results show that while the automated method does not perfectly reproduce the manual mapping, it successfully captures the general configuration, morphometry (length, width) and location of hummock corridors, despite variation in expression across and between sites. This technique is ideally suited to take advantage of newly available high-resolution digital elevation data (e.g. the ArcticDEM), whose enormous volume makes large-scale manual mapping prohibitively time consuming. Its application will enable efficient and comprehensive mapping of the spatial distribution of hummock corridors across palaeo-beds that is necessary for deriving insights into their formation and the organisation of subglacial meltwater flow beneath ice sheets