GLOBAL REGISTRATION OF MULTIPLE POINT CLOUDS USING SEMIDEFINITE PROGRAMMING

Consider N points in R-d and M local coordinate systems that are related through unknown rigid transforms. For each point, we are given (possibly noisy) measurements of its local coordinates in some of the coordinate systems. Alternatively, for each coordinate system, we observe the coordinates of a subset of the points. The problem of estimating the global coordinates of the N points (up to a rigid transform) from such measurements comes up in distributed approaches to molecular conformation and sensor network localization, and also in computer vision and graphics. The least-squares formulation of this problem, although nonconvex, has a well-known closed-form solution when M = 2 (based on the singular value decomposition (SVD)). However, no closed-form solution is known for M >= 3. In this paper, we demonstrate how the least-squares formulation can be relaxed into a convex program, namely, a semidefinite program (SDP). By setting up connections between the uniqueness of this SDP and results from rigidity theory, we prove conditions for exact and stable recovery for the SDP relaxation. In particular, we prove that the SDP relaxation can guarantee recovery under more adversarial conditions compared to earlier proposed spectral relaxations, and we derive error bounds for the registration error incurred by the SDP relaxation. We also present results of numerical experiments on simulated data to confirm the theoretical findings. We empirically demonstrate that (a) unlike the spectral relaxation, the relaxation gap is mostly zero for the SDP (i.e., we are able to solve the original nonconvex least-squares problem) up to a certain noise threshold, and (b) the SDP performs significantly better than spectral and manifold-optimization methods, particularly at large noise levels.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

Characterization of friction at three contact pairs by photoelastic Isotropic point (IP)

Friction coefficient between a circular-disk periphery and V-block surface was determined by introducing the concept of isotropic point (IP) in isochromatic field of the disk under three-point symmetric loading. IP position on the symmetry axis depends on active coefficient of friction during experiment. We extend this work to asymmetric loading of circular disk in which case two frictional contact pairs out of three loading contacts, independently control the unconstrained IP location. Photoelastic experiment is conducted on particular case of asymmetric three-point loading of circular disk. Basics of digital image processing are used to extract few essential parameters from experimental image, particularly IP location. Analytical solution by Flamant for half plane with a concentrated load, is utilized to derive stress components for required loading configurations of the disk. IP is observed, in analytical simulations of three-point asymmetric normal loading, to move from vertical axis to the boundary along an ellipse-like curve. When friction is included in the analysis, IP approaches the center with increase in loading friction and it goes away with increase in support friction. With all these insights, using experimental IP information, friction angles at three contact pairs of circular disk under asymmetric loading, are determined.

A comparative study of failure features in aerospace grade unidirectional and bidirectional woven CFRP composite laminates under four-point bend fatigue loads

Damage mechanisms in unidirectional (UD) and bi-directional (BD) woven carbon fiber reinforced polymer (CFRP) laminates subjected to four point flexure, both in static and fatigue loadings, were studied. The damage progression in composites was monitored by observing the slopes of the load vs. deflection data that represent the stiffness of the given specimen geometry over a number of cycles. It was observed that the unidirectional composites exhibit gradual loss in stiffness whereas the bidirectional woven composites show a relatively quicker loss during stage II of fatigue damage progression. Both, the static and the fatigue failures in unidirectional carbon fiber reinforced polymer composites originates due to generation of cracks on compression face while in bidirectional woven composites the damage ensues from both the compression and the tensile faces. These observations are supported by a detailed fractographic analysis.

Spectral Clustering with Jensen-type kernels and their multi-point extensions

Motivated by multi-distribution divergences, which originate in information theory, we propose a notion of `multipoint' kernels, and study their applications. We study a class of kernels based on Jensen type divergences and show that these can be extended to measure similarity among multiple points. We study tensor flattening methods and develop a multi-point (kernel) spectral clustering (MSC) method. We further emphasize on a special case of the proposed kernels, which is a multi-point extension of the linear (dot-product) kernel and show the existence of cubic time tensor flattening algorithm in this case. Finally, we illustrate the usefulness of our contributions using standard data sets and image segmentation tasks.

Asymptotic Bounds on the Power-Delay Tradeoff for Fading Point-to-Point Links From Geometric Bounds on the Stationary Distribution of the Queue Length

The optimal power-delay tradeoff is studied for a time-slotted independently and identically distributed fading point-to-point link, with perfect channel state information at both transmitter and receiver, and with random packet arrivals to the transmitter queue. It is assumed that the transmitter can control the number of packets served by controlling the transmit power in the slot. The optimal tradeoff between average power and average delay is analyzed for stationary and monotone transmitter policies. For such policies, an asymptotic lower bound on the minimum average delay of the packets is obtained, when average transmitter power approaches the minimum average power required for transmitter queue stability. The asymptotic lower bound on the minimum average delay is obtained from geometric upper bounds on the stationary distribution of the queue length. This approach, which uses geometric upper bounds, also leads to an intuitive explanation of the asymptotic behavior of average delay. The asymptotic lower bounds, along with previously known asymptotic upper bounds, are used to identify three new cases where the order of the asymptotic behavior differs from that obtained from a previously considered approximate model, in which the transmit power is a strictly convex function of real valued service batch size for every fade state.

Harmonic analysis of DC capacitor current in sinusoidal and space-vector modulated neutral-point-clamped inverters

The voltage ripple and power loss in the DC-capacitor of a voltage source inverter depend on the harmonic currents flowing through the capacitor. This paper presents a double Fourier series based analysis of the harmonic contents of the DC capacitor current in a three-level neutral-point clamped (NPC) inverter, modulated with sine-triangle pulse-width modulation (SPWM) or conventional space vector pulse-width modulation (CSVPWM) schemes. The analytical results are validated experimentally on a 3-kVA three-level inverter prototype. The capacitor current in an NPC inverter has a periodicity of 120(a similar to) at the fundamental or modulation frequency. Hence, this current contains third-harmonic and triplen-frequency components, apart from switching frequency components. The harmonic components vary with modulation index and power factor for both PWM schemes. The third harmonic current decreases with increase in modulation index and also decreases with increase in power factor in case of both PWM methods. In general, the third harmonic content is higher with SPWM than with CSVPWM at a given operating condition. Also, power loss and voltage ripple in the DC capacitor are estimated for both the schemes using the current harmonic spectrum and equivalent series resistance (ESR) of the capacitor.

Part I: Physical Insights Into the Two-Stage Breakdown Characteristics of STI-Type Drain-Extended pMOS Device

In this paper, we study breakdown characteristics in shallow-trench isolation (STI)-type drain-extended MOSFETs (DeMOS) fabricated using a low-power 65-nm triple-well CMOS process with a thin gate oxide. Experimental data of p-type STI-DeMOS device showed distinct two-stage behavior in breakdown characteristics in both OFF-and ON-states, unlike the n-type device, causing a reduction in the breakdown voltage and safe operating area. The first-stage breakdown occurs due to punchthrough in the vertical structure formed by p-well, deep n-well, and p-substrate, whereas the second-stage breakdown occurs due to avalanche breakdown of lateral n-well/p-well junction. The breakdown characteristics are also compared with the STI-DeNMOS device structure. Using the experimental results and advanced TCAD simulations, a complete understanding of breakdown mechanisms is provided in this paper for STI-DeMOS devices in advanced CMOS processes.

Part II: Design of Well Doping Profile for Improved Breakdown and Mixed-Signal Performance of STI-Type DePMOS Device

Shallow-trench isolation drain extended pMOS (STI-DePMOS) devices show a distinct two-stage breakdown. The impact of p-well and deep-n-well doping profile on breakdown characteristics is investigated based on TCAD simulations. Design guidelines for p-well and deep-n-well doping profile are developed to shift the onset of the first-stage breakdown to a higher drain voltage and to avoid vertical punch-through leading to early breakdown. An optimal ratio between the OFF-state breakdown voltage and the ON-state resistance could be obtained. Furthermore, the impact of p-well/deep-n-well doping profile on the figure of merits of analog and digital performance is studied. This paper aids in the design of STI drain extended MOSFET devices for widest safe operating area and optimal mixed-signal performance in advanced system-on-chip input-output process technologies.

Determinantal point processes in the plane from products of random matrices

We show that the density of eigenvalues for three classes of random matrix ensembles is determinantal. First we derive the density of eigenvalues of product of k independent n x n matrices with i.i.d. complex Gaussian entries with a few of matrices being inverted. In second example we calculate the same for (compatible) product of rectangular matrices with i.i.d. Gaussian entries and in last example we calculate for product of independent truncated unitary random matrices. We derive exact expressions for limiting expected empirical spectral distributions of above mentioned ensembles.

Transitions and blowoff of unconfined non-premixed swirling flame

The present experimental work reports the first observations of primary and secondary transitions in the time-averaged flame topology in a non-premixed swirling flame as the geometric swirl number S-G (a non dimensional number used to quantify the intensity of imparted swirl) is varied from a magnitude of zero till flame blowout. First observations of two transition types viz. primary and secondary transitions are reported. The primary transition represents a transformation from yellow straight jet flame (at S-G = 0) to lifted flame with blue base and finally to swirling seated (burner attached) yellow flame. Time-averaged streamline plot obtained from 2D PIV in mid-longitudinal plane shows a recirculation zone (RZ) at the immediate vicinity of burner exit. The lifted flame is stabilized along the vortex core of this RZ. Further, when S-G similar to 1.4-3, the first occurrence of vortex breakdown (VB) induced internal recirculation zone (IRZ) is witnessed. The flame now stabilizes at the upstream stagnation point of the VB-IRZ, which is attached to the burner lip. The secondary transition represents a transformation from a swirling seated flame to swirling flame with a conical tailpiece and finally to a highly-swirled near blowout oxidizer-rich flame. This transition is understood to be the result of transition in vortex breakdown modes of the swirling flow field from dual-ring VB bubble to central toroidal recirculation zone (CTRZ). The physics of transition is described on the basis of modified Rossby number (Ro(m)). Finally, when the swirl intensity is very high i.e. SG similar to 10, the flame blows out due to excessive straining and due to entrainment of large amount of oxidizer due to partial premixing. The present investigation involving changes in flame topology is immensely important because any change in global flame structure causes oscillatory heat release that can couple with dynamic pressure and velocity fluctuations leading to unsteady combustion. In this light, understanding mechanisms of flame stabilization is essential to tackle the problem of thermo-acoustic instability. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Position-dependent velocity of an effective temperature point for the estimation of the thermal diffusivity of solids

A new approach is proposed to estimate the thermal diffusivity of optically transparent solids at ambient temperature based on the velocity of an effective temperature point (ETP), and by using a two-beam interferometer the proposed concept is corroborated. 1D unsteady heat flow via step-temperature excitation is interpreted as a `micro-scale rectilinear translatory motion' of an ETP. The velocity dependent function is extracted by revisiting the Fourier heat diffusion equation. The relationship between the velocity of the ETP with thermal diffusivity is modeled using a standard solution. Under optimized thermal excitation, the product of the `velocity of the ETP' and the distance is a new constitutive equation for the thermal diffusivity of the solid. The experimental approach involves the establishment of a 1D unsteady heat flow inside the sample through step-temperature excitation. In the moving isothermal surfaces, the ETP is identified using a two-beam interferometer. The arrival-time of the ETP to reach a fixed distance away from heat source is measured, and its velocity is calculated. The velocity of the ETP and a given distance is sufficient to estimate the thermal diffusivity of a solid. The proposed method is experimentally verified for BK7 glass samples and the measured results are found to match closely with the reported value.