Source author record

Richard S. C. Cobbold

Richard S. C. Cobbold appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA physics.comp-ph Applications Computer Vision

Catalog footprint

What is connected

3works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

A second-order, perfectly matched layer formulation to model 3D transient wave propagation in anisotropic elastic media

Numerical simulation of wave propagation in an infinite medium is made possible by surrounding a finite region by a perfectly matched layer (PML). Using this approach a generalized three-dimensional (3D) formulation is proposed for time-domain modeling of elastic wave propagation in an unbounded lossless anisotropic medium. The formulation is based on a second-order approach that has the advantages of, physical relationship to the underlying equations, and amenability to be implemented in common numerical schemes. Specifically, our formulation uses three second-order equations of the displacement field and nine auxiliary equations, along with the three time histories of the displacement field. The properties of the PML, which are controlled by a complex two-parameter stretch function, are such that it acts as near perfect absorber. Using finite element method (FEM) 3D numerical results are presented for a highly anisotropic medium. An extension of the formulation to the particular case of a Kelvin-Vogit viscoelastic medium is also presented.

preprint2014arXiv

Efficient Low Dose X-ray CT Reconstruction through Sparsity-Based MAP Modeling

Ultra low radiation dose in X-ray Computed Tomography (CT) is an important clinical objective in order to minimize the risk of carcinogenesis. Compressed Sensing (CS) enables significant reductions in radiation dose to be achieved by producing diagnostic images from a limited number of CT projections. However, the excessive computation time that conventional CS-based CT reconstruction typically requires has limited clinical implementation. In this paper, we first demonstrate that a thorough analysis of CT reconstruction through a Maximum a Posteriori objective function results in a weighted compressive sensing problem. This analysis enables us to formulate a low dose fan beam and helical cone beam CT reconstruction. Subsequently, we provide an efficient solution to the formulated CS problem based on a Fast Composite Splitting Algorithm-Latent Expected Maximization (FCSA-LEM) algorithm. In the proposed method we use pseudo polar Fourier transform as the measurement matrix in order to decrease the computational complexity; and rebinning of the projections to parallel rays in order to extend its application to fan beam and helical cone beam scans. The weight involved in the proposed weighted CS model, denoted by Error Adaptation Weight (EAW), is calculated based on the statistical characteristics of CT reconstruction and is a function of Poisson measurement noise and rebinning interpolation error. Simulation results show that low computational complexity of the proposed method made the fast recovery of the CT images possible and using EAW reduces the reconstruction error by one order of magnitude. Recovery of a high quality 512$\times$ 512 image was achieved in less than 20 sec on a desktop computer without numerical optimizations.

preprint2013arXiv

Perfectly matched layer for second-order time-domain elastic wave equation: formulation and stability

A time domain system of equations is proposed to model elastic wave propagation in an unbounded two-dimensional anisotropic solid using perfectly matched layer (PML). Starting from a system of first-order frequency domain stress-velocity equations and using complex coordinate stretching approach with a two-parameter stretch function, a second-order formulation is obtained. The final system, which consists of just two second order equations along with four auxiliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. The discrete stability of the solutions for a given mesh size is examined with the help of a plane-wave analysis of the corresponding continuous problem. It is shown that increasing the scaling parameter of the stretch function leads to significant stability improvements for certain anisotropic media that have known issues. Numerical computations for different isotropic and anisotropic media are used to illustrate the results.