5,202 research outputs found
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 and a frequency derivative of
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 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 and
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
- …