Pseudo-time particle filtering for diffuse optical tomography

We recast the reconstruction problem of diffuse optical tomography (DOT) in a pseudo-dynamical framework and develop a method to recover the optical parameters using particle filters, i.e., stochastic filters based on Monte Carlo simulations. In particular, we have implemented two such filters, viz., the bootstrap (BS) filter and the Gaussian-sum (GS) filter and employed them to recover optical absorption coefficient distribution from both numerically simulated and experimentally generated photon fluence data. Using either indicator functions or compactly supported continuous kernels to represent the unknown property distribution within the inhomogeneous inclusions, we have drastically reduced the number of parameters to be recovered and thus brought the overall computation time to within reasonable limits. Even though the GS filter outperformed the BS filter in terms of accuracy of reconstruction, both gave fairly accurate recovery of the height, radius, and location of the inclusions. Since the present filtering algorithms do not use derivatives, we could demonstrate accurate contrast recovery even in the middle of the object where the usual deterministic algorithms perform poorly owing to the poor sensitivity of measurement of the parameters. Consistent with the fact that the DOT recovery, being ill posed, admits multiple solutions, both the filters gave solutions that were verified to be admissible by the closeness of the data computed through them to the data used in the filtering step (either numerically simulated or experimentally generated). (C) 2011 Optical Society of America

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

A pseudo-dynamical systems approach to a class of inverse problems in engineering

A pseudo-dynamical approach for a class of inverse problems involving static measurements is proposed and explored. Following linearization of the minimizing functional associated with the underlying optimization problem, the new strategy results in a system of linearized ordinary differential equations (ODEs) whose steady-state solutions yield the desired reconstruction. We consider some explicit and implicit schemes for integrating the ODEs and thus establish a deterministic reconstruction strategy without an explicit use of regularization. A stochastic reconstruction strategy is then developed making use of an ensemble Kalman filter wherein these ODEs serve as the measurement model. Finally, we assess the numerical efficacy of the developed tools against a few linear and nonlinear inverse problems of engineering interest.

On the moulding of population viscosity by natural selection

In this paper, we explore the conjoint evolution of dispersal and social behaviour. The model investigated is of a population distributed over a number of sites each with a carrying capacity of two adults and an episode of dispersal in the juvenile stage. The fertilities are governed by whether an individual and its neighbour are selfish or co-operative. It is shown that the best dispersal strategy for the co-operative genotype always involves lower levels of dispersal; and further that ecological conditions favouring low levels of dispersal increase the selective advantage of a co-operative genotype. Given this positive feedback, we suggest that in any taxon viscosity and co-operativity will tend to be correlated and bimodally distributed. Hence we predict the existence of two kinds of animal societies; viscous and co-operative (e.g. quasi-social wasps such as Mischocyttarus), and non-viscous and selfish (e.g. communal sphecid wasps such as Cerceris), and relatively few social groups with intermediate levels of co-operativity and viscosity. We also suggest that when one of the two sexes disperses, it will be the sex with lower potential for co-operative behaviour.

The Origin of Regioselectivity in a-Cleavage Reactions of Cyclopropenethiones: Potential Role of Pseudo-Jahn-Teller Effect in Substituted Cyclopropenyl Systems

Arylalkylcyclopropenethiones undergo highly regioselective photochemical a-cleavage via thioketene carbene intermediates, giving rise to products derived from the less stabilized carbene. UHF MIND0/3 calculations provide an insight into this unexpected regioselectivity. The nx* triplet of cyclopropenethione is calculated to have a highly unsymmetrical geometry with an elongated C-C bond, a delocalized thiaaUyl fragment, and a pyramidal radicaloid carbon (which eventually becomes the carbene center). From this molecular electronic structure, aryl group stabilization is expected to be more effective at the thiaallyl group rather than at the pyramidal radical center. Thus, the stability of the substituted triplet thione rather than that of the thioketene carbene determines the preferred regiochemistry of cleavage. The unusual structure of the cyclopropenethione triplet is suggested to be related to one of the Jahn-Teller distorted forms of the cyclopropenyl radical. An alternative symmetrical structure is adopted by the corresponding triplet of cyclopropenone, partly accounting for its differing photobehavior. A similar structural dichotomy is demonstrated for the corresponding radical anions as well.

Seismic reliability assessment of external stability of reinforced soil walls using pseudo-dynamic method

The paper presents a method for the evaluation of external stability of reinforced soil walls subjected to earthquakes in the framework of the pseudo-dynamic method. The seismic reliability of the wall is evaluated by considering the different possible failure modes such as sliding along the base, overturning about the toe point of the wall, bearing capacity and the eccentricity of the resultant force. The analysis is performed considering properties of the reinforced backfill, foundation soil below the base of the wall, length of the geosynthetic reinforcement and characteristics of earthquake ground motions such as shear wave and primary wave velocity as random variables. The optimum length of reinforcement needed to maintain stability against four modes of failure by targeting various component reliability indices is obtained. Differences between pseudo-static and pseudo-dynamic methods are clearly highlighted in the paper. A complete analysis of pseudo-static and pseudo-dynamic methodologies shows that the pseudodynamic method results in realistic design values for the length of geosynthetic reinforcement under earthquake conditions.

Stress-induced phase transformation and pseudo-elastic/pseudo-plastic recovery in intermetallic Ni-Al nanowires

Extensive molecular dynamics (MD) simulations have been performed in a B2-NiAl nanowire using an embedded atom method (EAM) potential. We show a stress induced B2 -> body-centered-tetragonal (BCT) phase transformation and a novel temperature and cross-section dependent pseudo-elastic/pseudo-plastic recovery from such an unstable BCT phase with a recoverable strain of similar to 30% as compared to 5-8% in polycrystalline materials. Such a temperature and cross-section dependent pseudo-elastic/pseudo-plastic strain recovery can be useful in various interesting applications of shape memory and strain sensing in nanoscale devices. Effects of size, temperature, and strain rate on the structural and mechanical properties have also been analyzed in detail. For a given size of the nanowire the yield stress of both the B2 and the BCT phases is found to decrease with increasing temperature, whereas for a given temperature and strain rate the yield stress of both the B2 and the BCT phase is found to increase with increase in the cross-sectional dimensions of the nanowire. A constant elastic modulus of similar to 80 GPa of the B2 phase is observed in the temperature range of 200-500 K for nanowires of cross-sectional dimensions in the range of 17.22-28.712 angstrom, whereas the elastic modulus of the BCT phase shows a decreasing trend with an increase in the temperature.

Development of pseudo-two-dimensional turbulent wakes

Pseudotwo-dimensional wakes are generated by introducing spanwise cellular structures into an otherwise plane turbulent wake by means of the castellated blunt trailing edges of different configurations. The transverse growths of these coflowing cellular wakes are found to be independent of each other without any noticeable spanwise interaction. This wake growth is examined in the light of the plane equilibrium wake analysis. Though these wakes are not found to be exactly self-similar, their growth shows a nonmonotonous approach toward the asymptotic state appropriate to that of a plane wake. The dye emission in the wakes illustrated a coherent vortical structure in the transverse plane, similar to that of the usual two-dimensional wake, in spite of the initial spanwise irregularities.

Aspects of tautomerism 9-solvolysis of normal and pseudo phthaloyl chlorides

Solvolysis of normal and pseudo phthaloyl chlorides in aqueous acetone and in aqueous dioxane has revealed that the former solvolyses about hundred-fold faster than the latter. Contrary to the accepted belief, there is no evidence for equilibrium between the normal and the pseudo forms of phthaloyl chloride, under a variety of conditions.

Aspects of tautomerism. 6. Base-catalyzed hydrolysis of pseudo esters of .gamma.-keto acids

The rates of base-catalyzed hydrolysis of pseudo esters derived from y-keto acids show strikingly poor sensitivity to the nature of the leaving group.la The rates vary in the narrow range of about 12-fold as contrasted to a 103-fold spread of the corresponding benzoate esters. The results presented are consistent with a rate-determining formation of a tetrahedral intermediate (11) and i t s rapid collapse, by the cleavage of the lactone ring in a fast step.

Aspects of tautomerism. 8. Solvolysis of some pseudo and normal acid chlorides

Pseudo acid chlorides derived from levulinic acid ando-benzoyl-benzoic acid, solvolyse in aqueous acetone, aqueous dioxane and aqueous dimethylformamide by aS Nl process. Their reaction pattern is distinct from that of typical normal acid chlorides, viz.,p-benzoylbenzoyl chloride and fluorene-9-one-1-carboxylic acid chloride, which solvolyse by aS N2 pathway. No evidence for tautomerism could be obtained either between the normal and pseudo forms of the acid chlorides or the derived ion pairs.

Pseudo-time marching schemes for inverse problems in structural health assessment and medical imaging

We propose a self-regularized pseudo-time marching strategy for ill-posed, nonlinear inverse problems involving recovery of system parameters given partial and noisy measurements of system response. While various regularized Newton methods are popularly employed to solve these problems, resulting solutions are known to sensitively depend upon the noise intensity in the data and on regularization parameters, an optimal choice for which remains a tricky issue. Through limited numerical experiments on a couple of parameter re-construction problems, one involving the identification of a truss bridge and the other related to imaging soft-tissue organs for early detection of cancer, we demonstrate the superior features of the pseudo-time marching schemes.

Prevalence and analysis of t(14;18) and t(11;14) chromosomal translocations in healthy Indian population

Hematopoietic malignancies like leukemia and lymphoma are characteristically associated with various chromosomal translocations. Follicular lymphoma (FL) and mantle cell lymphoma (MCL) are two subtypes of non-Hodgkin's lymphoma which possess t(14;18) and t(11;14) translocations, respectively. The incidence of FL and MCL is higher in the western countries as compared to India. Interestingly, the associated translocations are also found in healthy individuals in western population, which is 50-80% for t(14;18), whereas t(11;14) occurs at a very low frequency. However, there are no studies to explore thes translocations in healthy Indian population, which could explain the lower incidence of FL and MCL. We employed Southern hybridization following nested PCR to detect above translocations in healthy individuals from India. Our results suggest that this assay can detect one t(14;18) translocation event in up to 10(7) normal cells where as one t(11;14) in 10(8) normal cells. According to our results, 87 out of 253 individuals carry t(14;18) indicating 34% prevalence in the population. The presence of this translocation was also detectable at the transcript level. Although, no gender-based difference was observed, an age-dependent increase in the prevalence of translocation was found in adults. However, even after studying 210 people, we could not detect any t(11;14) translocation, indicating that it is uncommon in Indian population. These results suggest that lower incidence of FL and MCL in India could be attributed to lower prevalence of these translocations in healthy individuals.

Convergence analysis of the Newton algorithm and a pseudo-time marching scheme for diffuse correlation tomography

We propose a self-regularized pseudo-time marching scheme to solve the ill-posed, nonlinear inverse problem associated with diffuse propagation of coherent light in a tissuelike object. In particular, in the context of diffuse correlation tomography (DCT), we consider the recovery of mechanical property distributions from partial and noisy boundary measurements of light intensity autocorrelation. We prove the existence of a minimizer for the Newton algorithm after establishing the existence of weak solutions for the forward equation of light amplitude autocorrelation and its Frechet derivative and adjoint. The asymptotic stability of the solution of the ordinary differential equation obtained through the introduction of the pseudo-time is also analyzed. We show that the asymptotic solution obtained through the pseudo-time marching converges to that optimal solution provided the Hessian of the forward equation is positive definite in the neighborhood of optimal solution. The superior noise tolerance and regularization-insensitive nature of pseudo-dynamic strategy are proved through numerical simulations in the context of both DCT and diffuse optical tomography. (C) 2010 Optical Society of America.

A new framework for solution of multidimensional population balance equations

A new framework is proposed in this work to solve multidimensional population balance equations (PBEs) using the method of discretization. A continuous PBE is considered as a statement of evolution of one evolving property of particles and conservation of their n internal attributes. Discretization must therefore preserve n + I properties of particles. Continuously distributed population is represented on discrete fixed pivots as in the fixed pivot technique of Kumar and Ramkrishna [1996a. On the solution of population balance equation by discretization-I A fixed pivot technique. Chemical Engineering Science 51(8), 1311-1332] for 1-d PBEs, but instead of the earlier extensions of this technique proposed in the literature which preserve 2(n) properties of non-pivot particles, the new framework requires n + I properties to be preserved. This opens up the use of triangular and tetrahedral elements to solve 2-d and 3-d PBEs, instead of the rectangles and cuboids that are suggested in the literature. Capabilities of computational fluid dynamics and other packages available for generating complex meshes can also be harnessed. The numerical results obtained indeed show the effectiveness of the new framework. It also brings out the hitherto unknown role of directionality of the grid in controlling the accuracy of the numerical solution of multidimensional PBEs. The numerical results obtained show that the quality of the numerical solution can be improved significantly just by altering the directionality of the grid, which does not require any increase in the number of points, or any refinement of the grid, or even redistribution of pivots in space. Directionality of a grid can be altered simply by regrouping of pivots.