Accelerating frequency-domain diffuse optical tomographic image reconstruction using graphics processing units

Diffuse optical tomographic image reconstruction uses advanced numerical models that are computationally costly to be implemented in the real time. The graphics processing units (GPUs) offer desktop massive parallelization that can accelerate these computations. An open-source GPU-accelerated linear algebra library package is used to compute the most intensive matrix-matrix calculations and matrix decompositions that are used in solving the system of linear equations. These open-source functions were integrated into the existing frequency-domain diffuse optical image reconstruction algorithms to evaluate the acceleration capability of the GPUs (NVIDIA Tesla C 1060) with increasing reconstruction problem sizes. These studies indicate that single precision computations are sufficient for diffuse optical tomographic image reconstruction. The acceleration per iteration can be up to 40, using GPUs compared to traditional CPUs in case of three-dimensional reconstruction, where the reconstruction problem is more underdetermined, making the GPUs more attractive in the clinical settings. The current limitation of these GPUs in the available onboard memory (4 GB) that restricts the reconstruction of a large set of optical parameters, more than 13, 377. (C) 2010 Society of Photo-Optical Instrumentation Engineers. DOI: 10.1117/1.3506216]

Fracture Properties of Concrete-Concrete Interfaces Using Digital Image Correlation

The mode I and mode II fracture toughness and the critical strain energy release rate for different concrete-concrete jointed interfaces are experimentally determined using the Digital Image Correlation technique. Concrete beams having different compressive strength materials on either side of a centrally placed vertical interface are prepared and tested under three-point bending in a closed loop servo-controlled testing machine under crack mouth opening displacement control. Digital images are captured before loading (undeformed state) and at different instances of loading. These images are analyzed using correlation techniques to compute the surface displacements, strain components, crack opening and sliding displacements, load-point displacement, crack length and crack tip location. It is seen that the CMOD and vertical load-point displacement computed using DIC analysis matches well with those measured experimentally.

Grinding abrasive wear and associated particle size effect

The type of abrasion that the grinding medium experiences inside a ball mill is classified as high stress or grinding abrasion, because the stress levels at the surface of the medium exceed the yield stress of the metal when hard abrasives are crushed. During dry grinding of ores the medium undergoes not only abrasion but also erosion and impact. As all three mechanisms of wear occur simultaneously, it is difficult to follow the individual components of wear. However, it is possible to show that the overall kinetics of wear follows a simple power law of the type w = at(b), where w is the weight loss of the grinding medium for a specified grinding time t and a and b are constants. Experimental data, obtained from dry grinding of quartz for a wide range of times using AISI 52100 steel balls having various microstructures in a laboratory scale batch mill, are fitted to the proposed equation and the wear rate w is calculated from the first derivative of the equation. The mean particle sizes of the quartz charge DBAR corresponding to 50 and 80% retained size are determined by mechanical sieving of the ground product after a grinding time t and thus the relationship between wear rate and particle size of the abrasive is established. It is found that w increases rapidly with DBAR up to some critical size and then increases at a much lower rate.

Ground state geometric distortions Image distal substituent effects in determining the β-facial selectivity in 7-norbornenones

The X-ray structure of Image and MNDO optimized geometries of related 7-norbornenone derivatives show a clear tilt of the carbonyl bridge away from the C=C double bond. The preferred reduction from the more hindered face of the diester reveals the electron/electrostatic origin of π - facial selectivity in these systems. X-ray structure and MNDO calculations reveal the dominance of electronic effects in determining the π-facial selectivity in 4a.

Accelerated gradient based diffuse optical tomographic image reconstruction

Purpose: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). Methods: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. Results: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. Conclusions: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]

Controlled combustion synthesis and properties of fine-particle NASICON materials

Fine-particle NASICON materials, Na1+xZr2P3-xSixO12 (where x = 0.0, 0.5, 1.0, 1.5, 2.0 and 2.5), have been prepared by controlled combustion of an aqueous solution containing stoicthiometric amounts of sodium nitrate, zirconyl nitrate, ammonium perchlorate, diammonium hydrogen phosphate, fumed silica and carbonohydrazide. Formation of NASICON has been confirmed by powder XRD, Si-29 NMR and IR spectroscopy. These NASICON powders are fine (average agglomerate size 5-12 mum) with a surface area varying from 8 to 30 m2 g-1. NASICON powders pelletized and sintered at 1100-1200-degrees-C for 5 h achieved 90-95% theoretical density and show fine-grain microstructure. The coefficient of thermal expansion of sintered NASICON compact was measured up to 500-degrees-C and changes f rom -3.4 x 10(-6) to 4.1 x 10(-6) K-1. The conductivity of Sintered Na3Zr2PSi2O12 compact at 300-degrees-C is 0.236 OMEGA-1 cm-1.

Particle filters for structural system identification using multiple test and sensor data: A combined computational and experimental study

The problem of structural system identification when measurements originate from multiple tests and multiple sensors is considered. An offline solution to this problem using bootstrap particle filtering is proposed. The central idea of the proposed method is the introduction of a dummy independent variable that allows for simultaneous assimilation of multiple measurements in a sequential manner. The method can treat linear/nonlinear structural models and allows for measurements on strains and displacements under static/dynamic loads. Illustrative examples consider measurement data from numerical models and also from laboratory experiments. The results from the proposed method are compared with those from a Kalman filter-based approach and the superior performance of the proposed method is demonstrated. Copyright (C) 2009 John Wiley & Sons, Ltd.

Symmetry Properties of the Large-Deviation Function of the Velocity of a Self-Propelled Polar Particle

A geometrically polar granular rod confined in 2D geometry, subjected to a sinusoidal vertical oscillation, undergoes noisy self-propulsion in a direction determined by its polarity. When surrounded by a medium of crystalline spherical beads, it displays substantial negative fluctuations in its velocity. We find that the large-deviation function (LDF) for the normalized velocity is strongly non-Gaussian with a kink at zero velocity, and that the antisymmetric part of the LDF is linear, resembling the fluctuation relation known for entropy production, even when the velocity distribution is clearly non-Gaussian. We extract an analogue of the phase-space contraction rate and find that it compares well with an independent estimate based on the persistence of forward and reverse velocities.

A Non-contact measurement technique to measure micro surface stress and obtain deformation profiles of the order of 1nm in Micro-cantilever based structures by single image Optical Diffraction method

A new method based on analysis of a single diffraction pattern is proposed to measure deflections in micro-cantilever (MC) based sensor probes, achieving typical deflection resolutions of 1nm and surface stress changes of 50 mu N/m. The proposed method employs a double MC structure where the deflection of one of the micro-cantilevers relative to the other due to surface stress changes results in a linear shift of intensity maxima of the Fraunhofer diffraction pattern of the transilluminated MC. Measurement of such shifts in the intensity maxima of a particular order along the length of the structure can be done to an accuracy of 0.01mm leading to the proposed sensitivity of deflection measurement in a typical microcantilever. This method can overcome the fundamental measurement sensitivity limit set by diffraction and pointing stability of laser beam in the widely used Optical Beam Deflection method (OBDM).

PVU and wave-particle splitting schemes for Euler equations of gas dynamics

A new way of flux-splitting, termed as the wave-particle splitting is presented for developing upwind methods for solving Euler equations of gas dynamics. Based on this splitting, two new upwind methods termed as Acoustic Flux Vector Splitting (AFVS) and Acoustic Flux Difference Splitting (AFDS) methods are developed. A new Boltzmann scheme, which closely resembles the wave-particle splitting, is developed using the kinetic theory of gases. This method, termed as Peculiar Velocity based Upwind (PVU) method, uses the concept of peculiar velocity for upwinding. A special feature of all these methods that the unidirectional and multidirectional parts of the flux vector are treated separately. Extensive computations done using these schemes demonstrate the soundness of the ideas.

Synthesis and characterization of barium bis(citrato) oxozirconate(iv) tetrahydrate--A new molecular precursor for fine particle BaZrO3

Barium metazirconate (BaZrO3) fine powder has been produced by thermally decomposing a molecular precursor, barium bis(citrato)oxozirconate(IV) tetrahydrate at about 700-degrees-C. The precursor, Ba[ZrO(C6H6O7)2] . 4H2O (BZO) has been synthesized and characterized by employing a combination of spectroscopic and thermoanalytical techniques. The precursor undergoes thermal decomposition in three major stages: (i) dehydration to give an anhydrous barium zirconyl citrate, (ii) decomposition of the anhydrous citrate in a multistep process to form an ionic oxycarbonate intermediate, Ba2Zr2O5CO3, and (iii) decomposition of the oxycarbonate to produce BaZrO3 fine powder. The particle size of the resultant BaZrO3 is about 0.2 mum, and the surface area is found to be 4.0 m2 g-1.

Application of scaled particle theory to alkyl-aryl surfactants: a two zone model of micellar core

Scaled Particle Theory (SPT) has been applied to predict the total free energies of micellization of ionic as well as nonionic micellar systems containing an aryl ring. A modification of the previously developed model has been made, proposing a two-zone model of micellar core which corroborates with the structural information available for such systems. The results are in good agreement with experimental data and also confirm the dictating role of cavity forming free energies for such systems

Microwave assisted synthesis and UV-Vis spectroscopic studies of silver nanoparticles synthesized using vanillin as a reducing agent

The present investigation explores the adaptability of a microwave assisted route to obtain silver nanoparticles by the reduction of AgNO3 with vanillin, an environmentally benign material. Anionic surfactants such as AOT and SDS were used separately for encapsulating AgNPs and their role was compared. The UV-Visible absorption spectra present a broad SPR band consisting of two peaks suggesting the formation of silver nanoparticle with bimodal size distribution. The TEM image shows particles with spherical and hexagonal morphologies which confirms the results of UV-Vis studies. The anisotropy in the particle morphology can be attributed to the surface oxidation which in turn produces Ag@Ag2O core-shell nanostructures. Thus an intriguing feature of this system is that the obtained colloid is a mixture of AgNPs with and without Ag2O layers. Studies on the influence of pH on the stability of the synthesized nanoparticles revealed that the presence of excess Ag2O layers has a profound influence on it. Ag2O layers can be removed from AgNPs' surface by changing the solution pH to the acidic regime. The present study attests the enhanced ability of AOT in stabilizing the AgNPs in aqueous media. (C) 2011 Elsevier B.V. All rights reserved.

Single-molecule thermodynamics: the heat distribution function of a charged particle in a static magnetic field

Interest in the applicability of fluctuation theorems to the thermodynamics of single molecules in external potentials has recently led to calculations of the work and total entropy distributions of Brownian oscillators in static and time-dependent electromagnetic fields. These calculations, which are based on solutions to a Smoluchowski equation, are not easily extended to a consideration of the other thermodynamic quantity of interest in such systems-the heat exchanges of the particle alone-because of the nonlinear dependence of the heat on a particle's stochastic trajectory. In this paper, we show that a path integral approach provides an exact expression for the distribution of the heat fluctuations of a charged Brownian oscillator in a static magnetic field. This approach is an extension of a similar path integral approach applied earlier by our group to the calculation of the heat distribution function of a trapped Brownian particle, which was found, in the limit of long times, to be consistent with experimental data on the thermal interactions of single micron-sized colloids in a viscous solvent.

Image analysis using multifractals

The efficacy of the multifractal spectrum as a tool for characterizing images has been studied. This spectrum has been computed for digitized images of the nucleus of human cervical cancer cells and it was observed that the entire spectrum is almost fully reproduced for a normal cell while only the right half (q<0) of the spectrum is reproduced for a cancerous cell. Cells in stages in between the two extremes show a shortening of the left half of the spectrum proportional to their condition. The extent of this shortening has been found to be sufficient to permit a classification between three classes of cells at varying distances from a basal cancerous layer-the superficial cells, the intermediate cells and the parabasal cells. This technique may be used for automatic screening of the population while also indicating the stage of malignancy