Material parameter estimation and hypothesis testing on a 1D viscoelastic stenosis model: Methodology

This is the post-print version of the final published paper that is available from the link below. Copyright @ 2013 Walter de Gruyter GmbH.Non-invasive detection, localization and characterization of an arterial stenosis (a blockage or partial blockage in the artery) continues to be an important problem in medicine. Partial blockage stenoses are known to generate disturbances in blood flow which generate shear waves in the chest cavity. We examine a one-dimensional viscoelastic model that incorporates Kelvin–Voigt damping and internal variables, and develop a proof-of-concept methodology using simulated data. We first develop an estimation procedure for the material parameters. We use this procedure to determine confidence intervals for the estimated parameters, which indicates the efficacy of finding parameter estimates in practice. Confidence intervals are computed using asymptotic error theory as well as bootstrapping. We then develop a model comparison test to be used in determining if a particular data set came from a low input amplitude or a high input amplitude; this we anticipate will aid in determining when stenosis is present. These two thrusts together will serve as the methodological basis for our continuing analysis using experimental data currently being collected.National Institute of Allergy and Infectious Diseases, Air Force Office of Scientific Research, Department of Education, and Engineering and Physical Sciences Research Council

Magnetic hard-direction ordering in anisotropic Kondo systems

We present a generic mechanism that explains why many Kondo materials show magnetic ordering along directions that are not favoured by the crystal-field anisotropy. Using a renormalization-group (RG) analysis of single impurity Kondo models with single-ion anisotropy, we demonstrate that strong fluctuations above the Kondo temperature drive a moment re-orientation over a wide range of parameters, e.g. for different spin values $S$ and number of Kondo channels $N$. In tetragonal systems this can happen for both easy-plane or easy axis anisotropy. The characteristic crossing of magnetic susceptibilities is not an artefact of the weak-coupling RG treatment but can be reproduced in brute-force perturbation theory. Employing numerical renormalization group (NRG), we show that for an under-screened moment ($S=1$, $N=1$) with easy-plane anisotropy, a crossing of magnetic susceptibilities can also occur in the strong-coupling regime (below the Kondo temperature). This suggests that collective magnetic ordering of such under-screened moments would develop along the magnetic hard axis

High-order space-time finite element schemes for acoustic and viscodynamic wave equations with temporal decoupling

Copyright @ 2014 The Authors. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.We revisit a method originally introduced by Werder et al. (in Comput. Methods Appl. Mech. Engrg., 190:6685–6708, 2001) for temporally discontinuous Galerkin FEMs applied to a parabolic partial differential equation. In that approach, block systems arise because of the coupling of the spatial systems through inner products of the temporal basis functions. If the spatial finite element space is of dimension D and polynomials of degree r are used in time, the block system has dimension (r + 1)D and is usually regarded as being too large when r > 1. Werder et al. found that the space-time coupling matrices are diagonalizable over inline image for r ⩽100, and this means that the time-coupled computations within a time step can actually be decoupled. By using either continuous Galerkin or spectral element methods in space, we apply this DG-in-time methodology, for the first time, to second-order wave equations including elastodynamics with and without Kelvin–Voigt and Maxwell–Zener viscoelasticity. An example set of numerical results is given to demonstrate the favourable effect on error and computational work of the moderately high-order (up to degree 7) temporal and spatio-temporal approximations, and we also touch on an application of this method to an ambitious problem related to the diagnosis of coronary artery disease

Linear magnetoresistance in commercial n-type silicon due to inhomogeneous doping

Free electron theory tells us that resistivity is independent of magnetic field. In fact, most observations match the semiclassical prediction of a magnetoresistance that is quadratic at low fields before saturating. However, a non-saturating linear magnetoresistance has been observed in exotic semiconductors such as silver chalcogenides, lightly-doped InSb, N-doped InAs, MnAs-GaAs composites, PrFeAsO, and epitaxial graphene. Here we report the observation of a large linear magnetoresistance in the ohmic regime in commonplace commercial n-type silicon wafer. It is well-described by a classical model of spatially fluctuating donor densities, and may be amplified by altering the aspect ratio of the sample to enhance current-jetting: increasing the width tenfold increased the magnetoresistance at 8 T from 445 % to 4707 % at 35 K. This physical picture may well offer insights into the large magnetoresistances recently observed in n-type and p-type Si in the non-ohmic regime.Comment: submitted to Nature Material

Session-ocaml: a session-based library with polarities and lenses

We propose session-ocaml, a novel library for session-typed concurrent/distributed programming in OCaml. Our technique solely relies on parametric polymorphism, which can encode core session type structures with strong static guarantees. Our key ideas are: ( ) polarised session types, which give an alternative formulation of duality enabling OCaml to automatically infer an appropriate session type in a session with a reasonable notational overhead; and ( ) a parameterised monad with a data structure called ‘slots’ manipulated with lenses, which can statically enforce session linearity and delegations. We show applications of session-ocaml including a travel agency usecase and an SMTP protocol

The Deep Poincare Map: A Novel Approach for Left Ventricle Segmentation

Precise segmentation of the left ventricle (LV) within cardiac MRI images is a prerequisite for the quantitative measurement of heart function. However, this task is challenging due to the limited availability of labeled data and motion artifacts from cardiac imaging. In this work, we present an iterative segmentation algorithm for LV delineation. By coupling deep learning with a novel dynamic-based labeling scheme, we present a new methodology where a policy model is learned to guide an agent to travel over the image, tracing out a boundary of the ROI – using the magnitude difference of the Poincaré map as a stopping criterion. Our method is evaluated on two datasets, namely the Sunnybrook Cardiac Dataset (SCD) and data from the STACOM 2011 LV segmentation challenge. Our method outperforms the previous research over many metrics. In order to demonstrate the transferability of our method we present encouraging results over the STACOM 2011 data, when using a model trained on the SCD dataset

The CMB and the measure of the multiverse

In the context of eternal inflation, cosmological predictions depend on the choice of measure to regulate the diverging spacetime volume. The spectrum of inflationary perturbations is no exception, as we demonstrate by comparing the predictions of the fat geodesic and causal patch measures. To highlight the effect of the measure---as opposed to any effects related to a possible landscape of vacua---we take the cosmological model, including the model of inflation, to be fixed. We also condition on the average CMB temperature accompanying the measurement. Both measures predict a 1-point expectation value for the gauge-invariant Newtonian potential, which takes the form of a (scale-dependent) monopole, in addition to a related contribution to the 3-point correlation function, with the detailed form of these quantities differing between the measures. However, for both measures both effects are well within cosmic variance. Our results make clear the theoretical relevance of the measure, and at the same time validate the standard inflationary predictions in the context of eternal inflation.Comment: 28 pages; v2: reference added, some clarification

Effect of particle size on the thermal conductivity of nanofluids containing metallic nanoparticles

A one-parameter model is presented for the thermal conductivity of nanofluids containing dispersed metallic nanoparticles. The model takes into account the decrease in thermal conductivity of metal nanoparticles with decreasing size. Although literature data could be correlated well using the model, the effect of the size of the particles on the effective thermal conductivity of the nanofluid could not be elucidated from these data. Therefore, new thermal conductivity measurements are reported for six nanofluids containing silver nanoparticles of different sizes and volume fractions. The results provide strong evidence that the decrease in the thermal conductivity of the solid with particle size must be considered when developing models for the thermal conductivity of nanofluids

Simulational study of anomalous tracer diffusion in hydrogels

In this article, we analyze different factors that affect the diffusion behavior of small tracer particles (as they are used e.g.in fluorescence correlation spectroscopy (FCS)) in the polymer network of a hydrogel and perform simulations of various simplified models. We observe, that under certain circumstances the attraction of a tracer particle to the polymer network strands might cause subdiffusive behavior on intermediate time scales. In theory, this behavior could be employed to examine the network structure and swelling behavior of weakly crosslinked hydrogels with the help of FCS.Comment: 11 pages, 11 figure