Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

Ultrasound Tomography has seen a revival of interest in the past decade, especially for breast imaging, due to improvements in both ultrasound and computing hardware. In particular, three-dimensional ultrasound tomography, a fully tomographic method in which the medium to be imaged is surrounded by ultrasound transducers, has become feasible. In this paper, a comprehensive derivation and study of a robust framework for large-scale bent-ray ultrasound tomography in 3D for a hemispherical detector array is presented. Two ray-tracing approaches are derived and compared. More significantly, the problem of linking the rays between emitters and receivers, which is challenging in 3D due to the high number of degrees of freedom for the trajectory of rays, is analysed both as a minimisation and as a root-finding problem. The ray-linking problem is parameterised for a convex detection surface and three robust, accurate, and efficient ray-linking algorithms are formulated and demonstrated. To stabilise these methods, novel adaptive-smoothing approaches are proposed that control the conditioning of the update matrices to ensure accurate linking. The nonlinear UST problem of estimating the sound speed was recast as a series of linearised subproblems, each solved using the above algorithms and within a steepest descent scheme. The whole imaging algorithm was demonstrated to be robust and accurate on realistic data simulated using a full-wave acoustic model and an anatomical breast phantom, and incorporating the errors due to time-of-flight picking that would be present with measured data. This method can used to provide a low-artefact, quantitatively accurate, 3D sound speed maps. In addition to being useful in their own right, such 3D sound speed maps can be used to initialise full-wave inversion methods, or as an input to photoacoustic tomography reconstructions

Lithium-ion battery thermal-electrochemical model-based state estimation using orthogonal collocation and a modified extended Kalman filter

This paper investigates the state estimation of a high-fidelity spatially resolved thermal- electrochemical lithium-ion battery model commonly referred to as the pseudo two-dimensional model. The partial-differential algebraic equations (PDAEs) constituting the model are spatially discretised using Chebyshev orthogonal collocation enabling fast and accurate simulations up to high C-rates. This implementation of the pseudo-2D model is then used in combination with an extended Kalman filter algorithm for differential-algebraic equations to estimate the states of the model. The state estimation algorithm is able to rapidly recover the model states from current, voltage and temperature measurements. Results show that the error on the state estimate falls below 1 % in less than 200 s despite a 30 % error on battery initial state-of-charge and additive measurement noise with 10 mV and 0.5 K standard deviations.Comment: Submitted to the Journal of Power Source

Accurate Modelling of Left-Handed Metamaterials Using Finite-Difference Time-Domain Method with Spatial Averaging at the Boundaries

The accuracy of finite-difference time-domain (FDTD) modelling of left-handed metamaterials (LHMs) is dramatically improved by using an averaging technique along the boundaries of LHM slabs. The material frequency dispersion of LHMs is taken into account using auxiliary differential equation (ADE) based dispersive FDTD methods. The dispersive FDTD method with averaged permittivity along the material boundaries is implemented for a two-dimensional (2-D) transverse electric (TE) case. A mismatch between analytical and numerical material parameters (e.g. permittivity and permeability) introduced by the time discretisation in FDTD is demonstrated. The expression of numerical permittivity is formulated and it is suggested to use corrected permittivity in FDTD simulations in order to model LHM slabs with their desired parameters. The influence of switching time of source on the oscillation of field intensity is analysed. It is shown that there exists an optimum value which leads to fast convergence in simulations.Comment: 17 pages, 7 figures, submitted to Journal of Optics A Nanometa special issu

Solving the Poisson equation on small aspect ratio domains using unstructured meshes

We discuss the ill conditioning of the matrix for the discretised Poisson equation in the small aspect ratio limit, and motivate this problem in the context of nonhydrostatic ocean modelling. Efficient iterative solvers for the Poisson equation in small aspect ratio domains are crucial for the successful development of nonhydrostatic ocean models on unstructured meshes. We introduce a new multigrid preconditioner for the Poisson problem which can be used with finite element discretisations on general unstructured meshes; this preconditioner is motivated by the fact that the Poisson problem has a condition number which is independent of aspect ratio when Dirichlet boundary conditions are imposed on the top surface of the domain. This leads to the first level in an algebraic multigrid solver (which can be extended by further conventional algebraic multigrid stages), and an additive smoother. We illustrate the method with numerical tests on unstructured meshes, which show that the preconditioner makes a dramatic improvement on a more standard multigrid preconditioner approach, and also show that the additive smoother produces better results than standard SOR smoothing. This new solver method makes it feasible to run nonhydrostatic unstructured mesh ocean models in small aspect ratio domains.Comment: submitted to Ocean Modellin

Computational convergence of the path integral for real dendritic morphologies

Neurons are characterised by a morphological structure unique amongst biological cells, the core of which is the dendritic tree. The vast number of dendritic geometries, combined with heterogeneous properties of the cell membrane, continue to challenge scientists in predicting neuronal input-output relationships, even in the case of sub-threshold dendritic currents. The Green’s function obtained for a given dendritic geometry provides this functional relationship for passive or quasi-active dendrites and can be constructed by a sum-over-trips approach based on a path integral formalism. In this paper, we introduce a number of efficient algorithms for realisation of the sum-over-trips framework and investigate the convergence of these algorithms on different dendritic geometries. We demonstrate that the convergence of the trip sampling methods strongly depends on dendritic morphology as well as the biophysical properties of the cell membrane. For real morphologies, the number of trips to guarantee a small convergence error might become very large and strongly affect computational efficiency. As an alternative, we introduce a highly-efficient matrix method which can be applied to arbitrary branching structures

Numerical Efficiency of Inverse Simulation Methods Applied to a Wheeled Rover

Extending the navigational capability of planetary rovers is essential for increasing the scientific outputs from such exploratory missions. In this paper a navigation method based on Inverse Simulation is applied to a four wheel rover. The method calculates the required control inputs to achieve a desired, specified response. Here this is a desired trajectory defined as a series of waypoints. Inverse Simulation considers the complete system dynamics of the rover to calculate the control input using an iterative, numerical Newton - Raphson scheme. The paper provides an insight into the numerical parameters that affect the performance of the method. Also, the influence of varying the timestep and the convergence tolerance is examined in terms of the quality of the calculated control input and the resulting trajectory, as well as the execution time. From this analysis a set of parameters and recommendations to successfully apply Inverse Simulation to a rover is presented