A pseudo-time EnKF incorporating shape based reconstruction for diffuse optical tomography

Purpose: The authors aim at developing a pseudo-time, sub-optimal stochastic filtering approach based on a derivative free variant of the ensemble Kalman filter (EnKF) for solving the inverse problem of diffuse optical tomography (DOT) while making use of a shape based reconstruction strategy that enables representing a cross section of an inhomogeneous tumor boundary by a general closed curve. Methods: The optical parameter fields to be recovered are approximated via an expansion based on the circular harmonics (CH) (Fourier basis functions) and the EnKF is used to recover the coefficients in the expansion with both simulated and experimentally obtained photon fluence data on phantoms with inhomogeneous inclusions. The process and measurement equations in the pseudo-dynamic EnKF (PD-EnKF) presently yield a parsimonious representation of the filter variables, which consist of only the Fourier coefficients and the constant scalar parameter value within the inclusion. Using fictitious, low-intensity Wiener noise processes in suitably constructed ``measurement'' equations, the filter variables are treated as pseudo-stochastic processes so that their recovery within a stochastic filtering framework is made possible. Results: In our numerical simulations, we have considered both elliptical inclusions (two inhomogeneities) and those with more complex shapes (such as an annular ring and a dumbbell) in 2-D objects which are cross-sections of a cylinder with background absorption and (reduced) scattering coefficient chosen as mu(b)(a)=0.01mm(-1) and mu('b)(s)=1.0mm(-1), respectively. We also assume mu(a) = 0.02 mm(-1) within the inhomogeneity (for the single inhomogeneity case) and mu(a) = 0.02 and 0.03 mm(-1) (for the two inhomogeneities case). The reconstruction results by the PD-EnKF are shown to be consistently superior to those through a deterministic and explicitly regularized Gauss-Newton algorithm. We have also estimated the unknown mu(a) from experimentally gathered fluence data and verified the reconstruction by matching the experimental data with the computed one. Conclusions: The PD-EnKF, which exhibits little sensitivity against variations in the fictitiously introduced noise processes, is also proven to be accurate and robust in recovering a spatial map of the absorption coefficient from DOT data. With the help of shape based representation of the inhomogeneities and an appropriate scaling of the CH expansion coefficients representing the boundary, we have been able to recover inhomogeneities representative of the shape of malignancies in medical diagnostic imaging. (C) 2012 American Association of Physicists in Medicine. [DOI: 10.1118/1.3679855]

Microstructure and texture evolution during accumulative roll bonding of aluminium alloys AA2219/AA5086 composite laminates

Accumulative roll bonding of two aluminium alloys, AA2219 and AA5086 was carried out up to 8 passes. During the course of ARB, the deformation inhomogeneity between the two alloy layers results in interfacial instability after the 4th pass, necking of the AA5086 layers after the 6th pass and fracture along the necked regions after the 7th and 8th pass. The EBSD analysis shows deformation bands along the interfaces after 8 passes of ARB. The ARB-processed materials predominantly show characteristic deformation texture components. The weak texture after the 2nd pass results from the combination of a weakly-textured starting AA2219 layer and a strongly-textured starting AA5086 layer. A strong deformation texture forms due to the high imposed strain after a higher number of ARB passes. Subgrain formation and related shear banding induces copper/S components in the case of the small elongated grains, while planar slip leads to the formation of brass component in the large elongated grains.

Evidence of nanoscale structural phase separation in large bandwidth La0.2Sr0.8MnO3

We have investigated the structural evolution of La0.2Sr0.8MnO3 using temperature dependent high resolution synchrotron x-ray diffraction technique. In a wide temperature range, La0.2Sr0.8MnO3 reveals nanoscale structural inhomogeneity consisting of cubic and tetragonal phases. The present results suggest that domains of nanometer size of the tetragonal (low temperature) phase start nucleating in the cubic (high temperature) phase even above the Neel temperature (T-N). The tetragonal phase fraction increases substantially below T-N. Detailed analysis suggests that the twinned phase is tetragonal, orbital ordered, and insulating. At temperatures below 170 K, a small amount of the cubic phase is retained. The present results reveal the significance of the connectivity between the nanoscale structural phase separation with the physical properties.

Common Ground Method of Current Injection in Electrical Impedance Tomography

Surface electrodes are essentially required to be switched for boundary data collection in electrical impedance tomography (Ell). Parallel digital data bits are required to operate the multiplexers used, generally, for electrode switching in ELT. More the electrodes in an EIT system more the digital data bits are needed. For a sixteen electrode system. 16 parallel digital data bits are required to operate the multiplexers in opposite or neighbouring current injection method. In this paper a common ground current injection is proposed for EIT and the resistivity imaging is studied. Common ground method needs only two analog multiplexers each of which need only 4 digital data bits and hence only 8 digital bits are required to switch the 16 surface electrodes. Results show that the USB based data acquisition system sequentially generate digital data required for multiplexers operating in common ground current injection method. The profile of the boundary data collected from practical phantom show that the multiplexers are operating in the required sequence in common ground current injection protocol. The voltage peaks obtained for all the inhomogeneity configurations are found at the accurate positions in the boundary data matrix which proved the sequential operation of multiplexers. Resistivity images reconstructed from the boundary data collected from the practical phantom with different configurations also show that the entire digital data generation module is functioning properly. Reconstructed images and their image parameters proved that the boundary data are successfully acquired by the DAQ system which in turn indicates a sequential and proper operation of multiplexers.

Criteria for thermally induced atomization and catastrophic breakup of acoustically levitated droplet

We report thermally induced instability leading to catastrophic breakup in acoustically levitated vaporizing fuel droplets. Change in surface tension and viscosity with increase in droplet temperature causes wide fluctuations in droplet aspect ratio. If the viscous damping of aspect ratio oscillation is not strong enough, the droplet goes through unbounded stretching. If the droplet exceeds a critical Weber number locally, a bag type and capillary wave induced atomization can occur, which leads to catastrophic breakup. A stability criterion has been established based on the inhomogeneity of Bernoulli (acoustic) pressure and surface tension of the droplet in terms of a local Weber number and Ohnesorge number. This instability is thermally induced in a droplet which does not experience instabilities without heating. (C) 2012 Elsevier Ltd. All rights reserved.

Investigation of the origin of glassiness in La0.5Sr0.5CoO3

A comprehensive magnetic study has been carried out on the two sets of La0.5Sr0.5CoO3 samples with a view to understand the origin of low temperature glassiness in the ferromagnetic state. The samples prepared by the conventional solid-state synthesis method show a low temperature shoulder in both dc magnetization as well as in the ac susceptibility measurements, which exhibit characteristics of glassiness such as the frequency dependence and memory effect. These observations suggest the existence of a distinct low temperature cluster-glass like phase within dominant ferromagnetic phase. But, once the same sample is properly homogenized by repeated grinding and annealing process, the low temperature glassy phase disappears, and it shows a pure ferromagnetic behavior. Our comparative study clearly reveals that the reentrant spin-glass like nature is not intrinsic to La0.5Sr0.5CoO3 system, in fact this is an outcome of the compositional inhomogeneity.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Correlated Conductance Fluctuations Close to the Berezinskii-Kosterlitz-Thouless Transition in Ultrathin NbN Films

We probe the presence of long-range correlations in phase fluctuations by analyzing the higher-order spectrum of resistance fluctuations in ultrathin NbN superconducting films. The non-Gaussian component of resistance fluctuations is found to be sensitive to film thickness close to the transition, which allows us to distinguish between mean field and Berezinskii-Kosterlitz-Thouless (BKT) type superconducting transitions. The extent of non-Gaussianity was found to be bounded by the BKT and mean field transition temperatures and depends strongly on the roughness and structural inhomogeneity of the superconducting films. Our experiment outlines a novel fluctuation-based kinetic probe in detecting the nature of superconductivity in disordered low-dimensional materials.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

Temperature dependent electrical behaviour of Cu2SnS3 films

The temperature dependent electrical properties of the dropcasted Cu2SnS3 films have been measured in the temperature range 140 K to 317 K. The log I versus root V plot shows two regions. The region at lower bias is due to electrode limited Schottky emission and the higher bias region is due to bulk limited Poole Frenkel emission. The ideality factor is calculated from the ln I versus V plot for different temperatures fitted with the thermionic emission model and is found to vary from 6.05 eV to 12.23 eV. This large value is attributed to the presence of defects or amorphous layer at the Ag / Cu2SnS3 interface. From the Richardson's plot the Richardson's constant and the barrier height were calculated. Owing to the inhomogeneity in the barrier heights, the Richardson's constant and the barrier height were also calculated from the modified Richardson's plot. The I-V-T curves were also fitted using the thermionic field emission model. The barrier heights were found to be higher than those calculated using thermionic emission model. From the fit of the I-V-T curves to the field emission model, field emission was seen to dominate in the low temperature range of 140 K to 177 K. The temperature dependent current graphs show two regions of different mechanisms. The log I versus 1000/T plot gives activation energies E-a1 = 0.367095 - 0.257682 eV and E-a2 = 0.038416 - 0.042452 eV. The log ( I/T-2) versus 1000/T graph gives trap depths Phi(o1) = 0.314159 - 0.204752 eV and Phi(o2) = 0.007425- 0.011163 eV. With increasing voltage the activation energy E-a1 and the trap depth Phi(o1) decrease. From the ln (IT1/ 2) versus 1/T-1/ 4 graph, the low temperature region is due to variable range hopping mechanism and the high temperature region is due to thermionic emission. (C) 2014 Author(s).

Epitaxial Co metal thin film grown by pulsed laser deposition using oxide target

We report here the growth of epitaxial Co metal thin film on c-plane sapphire by pulsed laser deposition (RD) using Co:ZnO target utilizing the composition inhomogeneity of the corresponding plasma. Two distinct plasma composition regions have been observed using heavily alloyed Co0.6Zn0.4O target. The central and intense region of the plasma grows Co:ZnO film; the extreme tail grows only Co metal with no trace of either ZnO or Co oxide In between the two extremes, mixed phases (Co +Co-oxides +Co:ZnO) were observed. The Co metal thin film grown in this way shows room temperature ferromagnetism with large in plane magnetization similar to 1288 emu cm(-3) and a coerciviLy of similar to 230 Oe with applied field parallel to the film-substrate interface. Carrier density of the film is similar to 10(22) cm(-3). The film is epiLaxial single phase Co metal which is confirmed by both X-ray diffraction and transmission electron microscopy characierizaLions. Planar Hall Effect (PHE) and Magneto Optic Kerr Effect (MOKE) measurements confirm that the film possesses similar attributes of Co metal. The result shows that the epiLaxial Co metal thin film can be grown from its oxides in the PLD. (C) 2014 Elsevier B.V. All rights reserved.

Monte Carlo simulation of light transport in turbid medium with embedded object-spherical, cylindrical, ellipsoidal, or cuboidal objects embedded within multilayered tissues

Monte Carlo modeling of light transport in multilayered tissue (MCML) is modified to incorporate objects of various shapes (sphere, ellipsoid, cylinder, or cuboid) with a refractive-index mismatched boundary. These geometries would be useful for modeling lymph nodes, tumors, blood vessels, capillaries, bones, the head, and other body parts. Mesh-based Monte Carlo (MMC) has also been used to compare the results from the MCML with embedded objects (MCML-EO). Our simulation assumes a realistic tissue model and can also handle the transmission/reflection at the object-tissue boundary due to the mismatch of the refractive index. Simulation of MCML-EO takes a few seconds, whereas MMC takes nearly an hour for the same geometry and optical properties. Contour plots of fluence distribution from MCML-EO and MMC correlate well. This study assists one to decide on the tool to use for modeling light propagation in biological tissue with objects of regular shapes embedded in it. For irregular inhomogeneity in the model (tissue), MMC has to be used. If the embedded objects (inhomogeneity) are of regular geometry (shapes), then MCML-EO is a better option, as simulations like Raman scattering, fluorescent imaging, and optical coherence tomography are currently possible only with MCML. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Studies and Evaluation of EIT Image Reconstruction in EIDORS with Simulated Boundary Data

Simulated boundary potential data for Electrical Impedance Tomography (EIT) are generated by a MATLAB based EIT data generator and the resistivity reconstruction is evaluated with Electrical Impedance Tomography and Diffuse Optical Tomography Reconstruction Software (EIDORS). Circular domains containing subdomains as inhomogeneity are defined in MATLAB-based EIT data generator and the boundary data are calculated by a constant current simulation with opposite current injection (OCI) method. The resistivity images reconstructed for different boundary data sets and images are analyzed with image parameters to evaluate the reconstruction.

Conductivity noise study of the insulator-metal transition and phase coexistence in epitaxial samarium nickelate thin films

Interaction between the lattice and the orbital degrees of freedom not only makes rare-earth nickelates unusually ``bad metal,'' but also introduces a temperature-driven insulator-metal phase transition. Here we investigate this insulator-metal phase transition in thin films of SmNiO3 using the slow time-dependent fluctuations (noise) in resistivity. The normalized magnitude of noise is found to be extremely large, being nearly eight orders of magnitude higher than thin films of common disordered metallic systems, and indicates electrical conduction via classical percolation in a spatially inhomogeneous medium. The higher-order statistics of the fluctuations indicate a strong non-Gaussian component of noise close to the transition, attributing the inhomogeneity to the coexistence of the metallic and insulating phases. Our experiment offers insight into the impact of lattice-orbital coupling on the microscopic mechanism of electron transport in the rare-earth nickelates.

Instabilities in a fluid overlying an inclined anisotropic and inhomogeneous porous layer

In this paper, linear stability analysis on a Newtonian fluid film flowing under the effect of gravity over an inclined porous medium saturated with the same fluid in isothermal condition is carried out. The focus is placed on the effect of the anisotropic and inhomogeneous variations in the permeability of the porous medium on the shear mode and surface mode instabilities. The fluid-porous system is modelled by a coupled two-dimensional Navier-Stokes/Darcy problem. The perturbation equations are solved numerically using the Chebyshev collocation method. Detailed stability characteristics as a function of the depth ratio (the ratio of the depth of the fluid layer to that of the porous layer), the anisotropic parameter (the ratio of the permeability in the direction of the basic flow to that in the direction transverse to the basic flow) and the inhomogeneity functions are presented.