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

Novel auditory motivated subband temporal envelope based fundamental frequency estimation algorithm

We address the problem of estimating the fundamental frequency of voiced speech. We present a novel solution motivated by the importance of amplitude modulation in sound processing and speech perception. The new algorithm is based on a cumulative spectrum computed from the temporal envelope of various subbands. We provide theoretical analysis to derive the new pitch estimator based on the temporal envelope of the bandpass speech signal. We report extensive experimental performance for synthetic as well as natural vowels for both realworld noisy and noise-free data. Experimental results show that the new technique performs accurate pitch estimation and is robust to noise. We also show that the technique is superior to the autocorrelation technique for pitch estimation.

3D Image Reconstruction From Truncated Helical Cone Beam Projection Data - A Linear Prediction Approach

With the introduction of 2D flat-panel X-ray detectors, 3D image reconstruction using helical cone-beam tomography is fast replacing the conventional 2D reconstruction techniques. In 3D image reconstruction, the source orbit or scanning geometry should satisfy the data sufficiency or completeness condition for exact reconstruction. The helical scan geometry satisfies this condition and hence can give exact reconstruction. The theoretically exact helical cone-beam reconstruction algorithm proposed by Katsevich is a breakthrough and has attracted interest in the 3D reconstruction using helical cone-beam Computed Tomography.In many practical situations, the available projection data is incomplete. One such case is where the detector plane does not completely cover the full extent of the object being imaged in lateral direction resulting in truncated projections. This result in artifacts that mask small features near to the periphery of the ROI when reconstructed using the convolution back projection (CBP) method assuming that the projection data is complete. A number of techniques exist which deal with completion of missing data followed by the CBP reconstruction. In 2D, linear prediction (LP)extrapolation has been shown to be efficient for data completion, involving minimal assumptions on the nature of the data, producing smooth extensions of the missing projection data.In this paper, we propose to extend the LP approach for extrapolating helical cone beam truncated data. The projection on the multi row flat panel detectors has missing columns towards either ends in the lateral direction in truncated data situation. The available data from each detector row is modeled using a linear predictor. The available data is extrapolated and this completed projection data is backprojected using the Katsevich algorithm. Simulation results show the efficacy of the proposed method.

Thermal Weight Functions and Stress Intensity Factors for Bonded Dissimilar Media Using Body Analogy

In this study, an analytical method is presented for the computation of thermal weight functions in two dimensional bi-material elastic bodies containing a crack at the interface and subjected to thermal loads using body analogy method. The thermal weight functions are derived for two problems of infinite bonded dissimilar media, one with a semi-infinite crack and the other with a finite crack along the interface. The derived thermal weight functions are shown to reduce to the already known expressions of thermal weight functions available in the literature for the respective homogeneous elastic body. Using these thermal weight functions, the stress intensity factors are computed for the above interface crack problems when subjected to an instantaneous heat source.

Provider-Customer Coalitional Games

Efficacy of commercial wireless networks can be substantially enhanced through large-scale cooperation among involved entities such as providers and customers. The success of such cooperation is contingent upon the design of judicious resource allocation strategies that ensure that the individuals' payoffs are commensurate to the resources they offer to the coalition. The resource allocation strategies depend on which entities are decision-makers and whether and how they share their aggregate payoffs. Initially, we consider the scenario where the providers are the only decision-makers and they do not share their payoffs. We formulate the resource allocation problem as a nontransferable payoff coalitional game and show that there exists a cooperation strategy that leaves no incentive for any subset of providers to split from the grand coalition, i.e., the core of the game is nonempty. To compute this cooperation strategy and the corresponding payoffs, we subsequently relate this game and its core to an exchange market setting and its equilibrium, which can be computed by several efficient algorithms. Next, we investigate cooperation when customers are also decision-makers and decide which provider to subscribe to based on whether there is cooperation. We formulate a coalitional game in this setting and show that it has a nonempty core. Finally, we extend the formulations and results to the cases where the payoffs are vectors and can be shared selectively.

The phenomenon of notch in an NEMP simulator and its reduction

The effect of variation in the switching instant of the output switch of the pulser circuit used in energizing an NEMP simulator on the voltage fed to the simulator and hence the electric field within the working volume of the simulator has been studied. Depending upon the instant at which the output switch closes, the amplitude and the wave shape of the voltage that is fed to the illuminator varies. This wave shape of the output voltage from the pulser circuit determines the shape and characteristics of the electric field within the working volume of the simulator. To study the effect of variation in the switching instant on the vertical electric field within the working volume, the vertical electric field has been computed in time and frequency domains. For certain switching instants, the electric field shows a sharp reduction in its amplitude after the peak which is called the notch. The presence of notch results in the test object not getting illuminated with all the frequencies of interest. The notch has been successfully reduced by suitably modifying the pulser circuit.

Coupling of High Voltage AC Power Line Fields to Metallic Pipelines

In the recent years, there has been a trend to run metallic pipelines carrying petroleum products and high voltage AC power lines parallel to each other in a relatively narrow strip of land. Due to this sharing of the right-of-way, verhead AC power line electric field may induce voltages on the metallic pipelines running in close vicinity leading to serious adverse effects. In this paper, the induced voltages on metallic pipelines running in close vicinity of high voltage power transmission lines have been computed. Before computing the induced voltages, an optimum configuration of the phase conductors based on the lowest conductor surface gradient and field under transmission line has been arrived at. This paper reports the conductor surface field gradients calculated for the various configurations. Also the electric fields under transmission line, for single circuit and double circuit (various phase arrangements) have been analyzed. Based on the above results, an optimum configuration giving the lowest field under the power line as well as the lowest conductor surface gradient has been arrived at and for this configuration, induced voltage on the pipeline has been computed using the Charge Simulation Method (CSM). For comparison, induced voltages on the pipeline has been computed for the various other phase configurations also.

Radiated Lightning Electromagnetic fields at different heights above ground

In this paper, the radiated electric and magnetic fields above a perfectly conducting ground at different heights from 10 m to 10 km and for lateral distances varying from 20 m to 10 km from a lightning return stroke channel are computed and the field waveforms are presented. It has been observed that the vertical electric field reverses its polarity with height and this height depends on the radial distance from the lightning channel. The magnitude of the horizontal electric field, on the other hand,increases with height up to a certain height and then reduces. The effect of variation in the rate of rise of lightning current (di/dt) and the velocity of return stroke current on the radiated electric and magnetic fields for the above heights and distances have also been studied. It is seen that the variation in maximum current derivative does not have a significant influence on the electric field when ground is assumed as a perfect conductor but it influences significantly the horizontal electric field when ground has finite conductivity. The velocity of propagation of return stroke current on the other hand has significant influence for both perfectly as well as finitely conducting ground conditions.

Thermodynamic study of Fe2O3 Fe2(SO4)3 equilibrium using an oxyanionic electrolyte (Na2SO4 I)

The emf of the cell, Pt, Ar + O2 + SO2 + SO3/Na2SO4-I/Fe2O2 + Fe2(SO4)3, Pt, has been measured in the temperature range 800 to 1000 K, using a gas mixture of known input composition as the reference electrode. The equilibrium composition of the reference gas at the measuring temperatures was computed using the thermodynamic data on the gaseous species reported in the literature. A mixture of ferric oxide and sulfate was kept in a closed system to ensure establishment of equilibrium partial pressure at the electrode. The cell was designed to avoid physical contact between Fe2(SO4)3 and Na2SO4 electrolyte. Uncertainties arising from the formation of sulfate solid solution were thus eliminated. The Gibbs’ energy of formation of ferric sulfate calculated from the emf is discussed in comparison with data reported in the literature. There is no evidence for the formation of oxysulfates in the Fe-S-0 system. Based on the results obtained in the present study for Fe2(SO4)3 and literature data for other phases, chemical potential diagrams have been constructed for the Fe-S-O system at 900 and 1100 K.

Spinel deoxidation equilibria in the Fe Mn Al O system - experiments and computation

The concentration and chemical potential of oxygen in liquid Fe--Mn alloys equilibrated with the spinel solution, (Fe, Mn)Al sub 2+2x O sub 4+3x , and alpha -Al sub 2 O sub 3 have been determined at 1873K as a function of manganese concentration. The composition of the spinel phase has been determined using electron probe microanalysis. The results are compared with data reported in the literature. The deoxidation equilibrium has been computed using data on free energy of solution of oxygen in liquid iron, free energies of formation of hercynite and galaxite, and interaction parameters reported in the literature. The activity--composition relationship in spinel solution was derived from a cation distribution model. The model is in excellent agreement with the experimental data on oxygen concentration and potential and the composition of the spinel phase. 23 ref.--AA

Output-Sensitive Construction of Reeb Graphs

The Reeb graph of a scalar function represents the evolution of the topology of its level sets. This paper describes a near-optimal output-sensitive algorithm for computing the Reeb graph of scalar functions defined over manifolds or non-manifolds in any dimension. Key to the simplicity and efficiency of the algorithm is an alternate definition of the Reeb graph that considers equivalence classes of level sets instead of individual level sets. The algorithm works in two steps. The first step locates all critical points of the function in the domain. Critical points correspond to nodes in the Reeb graph. Arcs connecting the nodes are computed in the second step by a simple search procedure that works on a small subset of the domain that corresponds to a pair of critical points. The paper also describes a scheme for controlled simplification of the Reeb graph and two different graph layout schemes that help in the effective presentation of Reeb graphs for visual analysis of scalar fields. Finally, the Reeb graph is employed in four different applications-surface segmentation, spatially-aware transfer function design, visualization of interval volumes, and interactive exploration of time-varying data.

Tension-Softening Properties for Concrete-Concrete Interfaces

The tension-softening parameters for different concrete-concrete interfaces are determined using the bimaterial cracked hinge model. Beams of different sizes having a jointed interface between two different strengths of concrete are tested under three-point bending (TPB). The load versus crack mouth opening displacement (CMOD) results are used to obtain the stress-crack opening relation through an inverse analysis. In addition, the fracture energy, tensile strength, and modulus of elasticity are also computed from the inverse analysis. The fracture properties are used in the nonlinear fracture mechanics analysis of a concrete patch-repaired beam to determine its load-carrying capacity when repaired with concrete of different strengths.

Measurement of the Gibbs energy of formation of URh3 using an oxide solid electrolyte

The Gibbs' energy offormation of the intermetallic compound URh3has been measured in the temperature range 980 to 1320 K using an oxide solid state cell incorporating yttria-doped thoria as the solid electrolyte and a mixture of manganese and manganese oxide as the reference electrode. The cell can be represented as Pt, Mn + MnO I (Y203)Th02 I Rh + URh3 + U02 + x' Rh, Pt The reversible emf of the cell was a linear function of temperature E = 15.60 +0.0237 T (±0.8) mY. Using auxiliary thermodynamic data for MnO and U02+ x the Gibbs' energy of formation of URh3 from component metals has been computed. The results can be expressed by the equation L'.G?< URh3 > = -316240 + 13.22 T (± 3000) J mol-1. The "third-law" enthalpy of formation of URh3at 298 K is -293.2 (± 4) kJ mol-1, significantly more negative than the value of -181.5 kJ mol-1 calculated using Miedema's model.

Gibbs’ energy of formation of YBa2Cu3O7-x (tetragonal)

The high temperature ceramic oxide superconductor YBa2Cu3O7-x (1–2–3 compound) is generally synthesized in an oxygen-rich environment. Hence any method for determining its thermodynamic stability should operate at a high oxygen partial pressure. A solid-state cell incorporating CaF2 as the electrolyte and functioning under pure oxygen at a pressure of 1·01 × 105 Pa has been employed for the determination of the Gibbs’ energy of formation of the 1–2–3 compound. The configuration of the galvanic cell can be represented by: Pt, O2, YBa2Cu3O7−x , Y2BaCuO5, CuO, BaF2/CaF2/BaF2, BaZrO3, ZrO2, O2, Pt. Using the values of the standard Gibbs’ energy of formation of the compounds BaZrO3 and Y2BaCuO5 from the literature, the Gibbs’ energy of formation of the 1–2–3 compound from the constituent binary oxides has been computed at different temperatures. The value ofx at each temperature is determined by the oxygen partial pressure. At 1023 K for O content of 6·5 the Gibbs’ energy of formation of the 1–2–3 compound is −261·7 kJ mol−1.

Thermodynamics of Aluminium -Strontium Alloys

The activity of strontium in liquid Al-Sr alloys (X(Sr) less-than-or-equal-to 0.17) at 1323 K has been determined using the Knudsen effusion-mass loss technique. At higher concentrations (X(Sr) greater-than-or-equal-to 0.28), the activity of strontium has been determined by the pseudoisopiestic technique. Activity of aluminium has been derived by Gibbs-Duhem integration. The concentration - concentration structure factor of Bhatia and Thornton at zero wave vector has been computed from the thermodynamic data. The behaviour of the mean square thermal fluctuation in composition and the thermodynamic mixing functions suggest association tendencies in the liquid state. The associated solution model with Al2Sr as the predominant complex can account for the properties of the liquid alloy. Thermodynamic data for the intermetallic compunds in the Al-Sr system have been derived using the phase diagram and the Gibbs' energy and enthalpy of mixing of liquid alloys. The data indicate the need for redetermination of the phase diagram near the strontium-rich corner.