Biquadratic pulse transmission networks

A second order transfer function with two poles and two zeros exhibits a step response characterized by a sudden rise to the steady state value, followed by oscillations around this steady state. With proper choice of the coefficients, it is possible to obtain transfer functions suitable for pulse transmission purposes.

Uncertainty visualization using HDR volume rendering

In this paper, we explore a novel idea of using high dynamic range (HDR) technology for uncertainty visualization. We focus on scalar volumetric data sets where every data point is associated with scalar uncertainty. We design a transfer function that maps each data point to a color in HDR space. The luminance component of the color is exploited to capture uncertainty. We modify existing tone mapping techniques and suitably integrate them with volume ray casting to obtain a low dynamic range (LDR) image. The resulting image is displayed on a conventional 8-bits-per-channel display device. The usage of HDR mapping reveals fine details in uncertainty distribution and enables the users to interactively study the data in the context of corresponding uncertainty information. We demonstrate the utility of our method and evaluate the results using data sets from ocean modeling.

Transition and acoustic response of recirculation structures in an unconfined co-axial isothermal swirling flow

This paper reports the first observations of transition from a pre-vortex breakdown (Pre-VB) flowreversal to a fully developed central toroidal recirculation zone in a non-reacting, double-concentric swirling jet configuration and its response to longitudinal acoustic excitation. This transition proceeds with the formation of two intermediate, critical flow regimes. First, a partially penetrated vortex breakdown bubble (VBB) is formed that indicates the first occurrence of an enclosed structure as the centre jet penetration is suppressed by the growing outer roll-up eddy; resulting in an opposed flow stagnation region. Second, a metastable transition structure is formed that marks the collapse of inner mixing vortices. In this study, the time-averaged topological changes in the coherent recirculation structures are discussed based on the non-dimensional modified Rossby number (Ro(m)) which appears to describe the spreading of the zone of swirl influence in different flow regimes. Further, the time-mean global acoustic response of pre-VB and VBB is measured as a function of pulsing frequency using the relative aerodynamic blockage factor (i.e., maximum radial width of the inner recirculation zone). It is observed that all flow modes except VBB are structurally unstable as they exhibit severe transverse radial shrinkage (similar to 20%) at the burner Helmholtz resonant modes (100-110 Hz). In contrast, all flow regimes show positional instability as seen by the large-scale, asymmetric spatial shifting of the vortex core centres. Finally, the mixing transfer function M (f) and magnitude squared coherence lambda(2)(f) analysis is presented to determine the natural couplingmodes of the system dynamic parameters (u', p'), i.e., local acoustic response. It is seen that the pre-VB flow mode exhibits a narrow-band, low pass filter behavior with a linear response window of 100-105 Hz. However, in the VBB structure, presence of critical regions such as the opposed flow stagnation region alters the linearity range with the structure showing a response even at higher pulsing frequencies (100-300 Hz). (C) 2013 AIP Publishing LLC.

Zero miss distance guidance using feedforward and periodic control

There have been attempts at obtaining robust guidance laws to ensure zero miss distance (ZMD) for interceptors with parametric uncertainties. All these laws require the plant to be of minimum phase type to enable the overall guidance loop transfer function to satisfy strict positive realness (SPR). The SPR property implies absolute stability of the closed loop system, and has been shown in the literature to lead to ZMD because it avoids saturation of lateral acceleration. In these works higher order interceptors are reduced to lower order equivalent models for which control laws are designed to ensure ZMD. However, it has also been shown that when the original system with right half plane (RHP) zeros is considered, the resulting miss distances, using such strategies, can be quite high. In this paper, an alternative approach using the circle criterion establishes the conditions for absolute stability of the guidance loop and relaxes the conservative nature of some earlier results arising from assumption of in�nite engagement time. Further, a feedforward scheme in conjunction with a lead-lag compensator is used as one control strategy while a generalized sampled hold function is used as a second strategy, to shift the RHP transmission zeros, thereby achieving ZMD. It is observed that merely shifting the RHP zero(s) to the left half plane reduces miss distances signi�cantly even when no additional controllers are used to ensure SPR conditions.

A radial differential mobility analyzer for the size-classification of gas-phase synthesized nanoparticles at low pressures

Differential mobility analyzers (DMAs) are commonly used to generate monodisperse nanoparticle aerosols. Commercial DMAs operate at quasi-atmospheric pressures and are therefore not designed to be vacuum-tight. In certain particle synthesis methods, the use of a vacuum-compatible DMA is a requirement as a process step for producing high-purity metallic particles. A vacuum-tight radial DMA (RDMA) has been developed and tested at low pressures. Its performance has been evaluated by using a commercial NANO-DMA as the reference. The performance of this low-pressure RDMA (LP-RDMA) in terms of the width of its transfer function is found to be comparable with that of other NANO-DMAs at atmospheric pressure and is almost independent of the pressure down to 30 mbar. It is shown that LP-RDMA can be used for the classification of nanometer-sized particles (5-20 nm) under low pressure condition (30 mbar) and has been successfully applied to nanoparticles produced by ablating FeNi at low pressures.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

Localized Excitations and the Morphology of Cooperatively Rearranging Regions in a Colloidal Glass-Forming Liquid

We develop a scheme based on a real space microscopic analysis of particle dynamics to ascertain the relevance of dynamical facilitation as a mechanism of structural relaxation in glass-forming liquids. By analyzing the spatial organization of localized excitations within clusters of mobile particles in a colloidal glass former and examining their partitioning into shell-like and corelike regions, we establish the existence of a crossover from a facilitation-dominated regime at low area fractions to a collective activated hopping-dominated one close to the glass transition. This crossover occurs in the vicinity of the area fraction at which the peak of the mobility transfer function exhibits a maximum and the morphology of cooperatively rearranging regions changes from stringlike to a compact form. Collectively, our findings suggest that dynamical facilitation is dominated by collective hopping close to the glass transition, thereby constituting a crucial step towards identifying the correct theoretical scenario for glass formation.

The statistics of the second order response of an FPSO in spreading seas

An expression for the probability density function of the second order response of a general FPSO in spreading seas is derived by using the Kac-Siegert approach. Various approximations of the second order force transfer functions are investigated for a ship-shaped FPSO. It is found that, when expressed in non-dimensional form, the probability density function of the response is not particularly sensitive to wave spreading, although the mean squared response and the resulting dimensional extreme values can be sensitive. The analysis is then applied to a Sevan FPSO, which is a large cylindrical buoy-like structure. The second order force transfer functions are derived by using an efficient semi-analytical hydrodynamic approach, and these are then employed to yield the extreme response. However, a significant effect of wave spreading on the statistics for a Sevan FPSO is found even in non-dimensional form. It implies that the exact statistics of a general ship-shaped FPSO may be sensitive to the wave direction, which needs to be verified in future work. It is also pointed out that the Newman's approximation regarding the frequency dependency of force transfer function is acceptable even for the spreading seas. An improvement on the results may be attained when considering the angular dependency exactly. Copyright © 2009 by ASME.

Computer generated hologram with geometric occlusion using GPU-accelerated depth buffer rasterization for three-dimensional display

We present a method of rapidly producing computer-generated holograms that exhibit geometric occlusion in the reconstructed image. Conceptually, a bundle of rays is shot from every hologram sample into the object volume.We use z buffering to find the nearest intersecting object point for every ray and add its complex field contribution to the corresponding hologram sample. Each hologram sample belongs to an independent operation, allowing us to exploit the parallel computing capability of modern programmable graphics processing units (GPUs). Unlike algorithms that use points or planar segments as the basis for constructing the hologram, our algorithm's complexity is dependent on fixed system parameters, such as the number of ray-casting operations, and can therefore handle complicated models more efficiently. The finite number of hologram pixels is, in effect, a windowing function, and from analyzing the Wigner distribution function of windowed free-space transfer function we find an upper limit on the cone angle of the ray bundle. Experimentally, we found that an angular sampling distance of 0:01' for a 2:66' cone angle produces acceptable reconstruction quality. © 2009 Optical Society of America.

Direct parameter extraction in feedthrough-embedded capacitive MEMS resonators

This paper presents a method for fast and accurate determination of parameters relevant to the characterization of capacitive MEMS resonators like quality factor (Q), resonant frequency (fn), and equivalent circuit parameters such as the motional capacitance (Cm). In the presence of a parasitic feedthrough capacitor (CF) appearing across the input and output ports, the transmission characteristic is marked by two resonances: series (S) and parallel (P). Close approximations of these circuit parameters are obtained without having to first de-embed the resonator motional current typically buried in feedthrough by using the series and parallel resonances. While previous methods with the same objective are well known, we show that these are limited to the condition where CF ≪ CmQ. In contrast, this work focuses on moderate capacitive feedthrough levels where CF > CmQ, which are more common in MEMS resonators. The method is applied to data obtained from the measured electrical transmission of fabricated SOI MEMS resonators. Parameter values deduced via direct extraction are then compared against those obtained by a full extraction procedure where de-embedding is first performed and followed by a Lorentzian fit to the data based on the classical transfer function associated with a generic LRC series resonant circuit. © 2011 Elsevier B.V. All rights reserved.

Diffraction-limited imaging on the 200-inch telescope

We have used the technique of non-redundant masking at the Palomar 200-inch telescope and radio VLBI imaging software to make optical aperture synthesis maps of two binary stars, β Corona Borealis and σ Herculis. The dynamic range of the map of β CrB, a binary star with a separation of 230 milliarcseconds is 50:1. For σ Her, we find a separation of 70 milliarcseconds and the dynamic range of our image is 30:1. These demonstrate the potential of the non-redundant masking technique for diffraction-limited imaging of astronomical objects with high dynamic range.

We find that the optimal integration time for measuring the closure phase is longer than that for measuring the fringe amplitude. There is not a close relationship between amplitude errors and phase errors, as is found in radio interferometry. Amplitude self calibration is less effective at optical wavelengths than at radio wavelengths. Primary beam sensitivity correction made in radio aperture synthesis is not necessary in optical aperture synthesis.

The effects of atmospheric disturbances on optical aperture synthesis have been studied by Monte Carlo simulations based on the Kolmogorov theory of refractive-index fluctuations. For the non-redundant masking with τ_c-sized apertures, the simulated fringe amplitude gives an upper bound of the observed fringe amplitude. A smooth transition is seen from the non-redundant masking regime to the speckle regime with increasing aperture size. The fractional reduction of the fringe amplitude according to the bandwidth is nearly independent of the aperture size. The limiting magnitude of optical aperture synthesis with τ_c-sized apertures and that with apertures larger than τ_c are derived.

Monte Carlo simulations are also made to study the sensitivity and resolution of the bispectral analysis of speckle interferometry. We present the bispectral modulation transfer function and its signal-to-noise ratio at high light levels. The results confirm the validity of the heuristic interferometric view of image-forming process in the mid-spatial-frequency range. The signal-to- noise ratio of the bispectrum at arbitrary light levels is derived in the mid-spatial-frequency range.

The non-redundant masking technique is suitable for imaging bright objects with high resolution and high dynamic range, while the faintest limit will be better pursued by speckle imaging.

Phase-space analysis of wavefront coding imaging systems

We explore the use of the Radon-Wigner transform, which is associated with the fractional Fourier transform of the pupil function, for determining the point spread function (PSF) of an incoherant defocused optical system. Then we introduce these phase-space tools to analyse the wavefront coding imaging system. It is shown that the shape of the PSF for such a system is highly invarient to the defocous-related aberrations except for a lateral shift. The optical transfer function of this system is also investigated briefly from a new understanding of ambiguity function.

Distributed Optimal Control of Cyber-Physical Systems: Controller Synthesis, Architecture Design and System Identification

The centralized paradigm of a single controller and a single plant upon which modern control theory is built is no longer applicable to modern cyber-physical systems of interest, such as the power-grid, software defined networks or automated highways systems, as these are all large-scale and spatially distributed. Both the scale and the distributed nature of these systems has motivated the decentralization of control schemes into local sub-controllers that measure, exchange and act on locally available subsets of the globally available system information. This decentralization of control logic leads to different decision makers acting on asymmetric information sets, introduces the need for coordination between them, and perhaps not surprisingly makes the resulting optimal control problem much harder to solve. In fact, shortly after such questions were posed, it was realized that seemingly simple decentralized optimal control problems are computationally intractable to solve, with the Wistenhausen counterexample being a famous instance of this phenomenon. Spurred on by this perhaps discouraging result, a concerted 40 year effort to identify tractable classes of distributed optimal control problems culminated in the notion of quadratic invariance, which loosely states that if sub-controllers can exchange information with each other at least as quickly as the effect of their control actions propagates through the plant, then the resulting distributed optimal control problem admits a convex formulation.

The identification of quadratic invariance as an appropriate means of "convexifying" distributed optimal control problems led to a renewed enthusiasm in the controller synthesis community, resulting in a rich set of results over the past decade. The contributions of this thesis can be seen as being a part of this broader family of results, with a particular focus on closing the gap between theory and practice by relaxing or removing assumptions made in the traditional distributed optimal control framework. Our contributions are to the foundational theory of distributed optimal control, and fall under three broad categories, namely controller synthesis, architecture design and system identification.

We begin by providing two novel controller synthesis algorithms. The first is a solution to the distributed H-infinity optimal control problem subject to delay constraints, and provides the only known exact characterization of delay-constrained distributed controllers satisfying an H-infinity norm bound. The second is an explicit dynamic programming solution to a two player LQR state-feedback problem with varying delays. Accommodating varying delays represents an important first step in combining distributed optimal control theory with the area of Networked Control Systems that considers lossy channels in the feedback loop. Our next set of results are concerned with controller architecture design. When designing controllers for large-scale systems, the architectural aspects of the controller such as the placement of actuators, sensors, and the communication links between them can no longer be taken as given -- indeed the task of designing this architecture is now as important as the design of the control laws themselves. To address this task, we formulate the Regularization for Design (RFD) framework, which is a unifying computationally tractable approach, based on the model matching framework and atomic norm regularization, for the simultaneous co-design of a structured optimal controller and the architecture needed to implement it. Our final result is a contribution to distributed system identification. Traditional system identification techniques such as subspace identification are not computationally scalable, and destroy rather than leverage any a priori information about the system's interconnection structure. We argue that in the context of system identification, an essential building block of any scalable algorithm is the ability to estimate local dynamics within a large interconnected system. To that end we propose a promising heuristic for identifying the dynamics of a subsystem that is still connected to a large system. We exploit the fact that the transfer function of the local dynamics is low-order, but full-rank, while the transfer function of the global dynamics is high-order, but low-rank, to formulate this separation task as a nuclear norm minimization problem. Finally, we conclude with a brief discussion of future research directions, with a particular emphasis on how to incorporate the results of this thesis, and those of optimal control theory in general, into a broader theory of dynamics, control and optimization in layered architectures.

Stability analysis of RF amplifiers based on MIMO pole-zero identification

Spurious oscillations are one of the principal issues faced by microwave and RF circuit designers. The rigorous detection of instabilities or the characterization of measured spurious oscillations is still an ongoing challenge. This project aims to create a new stability analysis CAD program that tackles this chal- lenge. Multiple Input Multiple Output (MIMO) pole-zero identification analysis is introduced on the program as a way to create new methods to automate the stability analysis process and to help designers comprehend the obtained results and prevent incorrect interpretations. The MIMO nature of the analysis contributes to eliminate possible controllability and observability losses and helps differentiate mathematical and physical quasi-cancellations, products of overmodeling. The created program reads Single Input Single Output (SISO) or MIMO frequency response data, and determines the corresponding continuous transfer functions with Vector Fitting. Once the transfer function is calculated, the corresponding pole/zero diagram is mapped enabling the designers to analyze the stability of an amplifier. Three data processing methods are introduced, two of which consist of pole/zero elimina- tions and the latter one on determining the critical nodes of an amplifier. The first pole/zero elimination method is based on eliminating non resonant poles, whilst the second method eliminates the poles with small residue by assuming that their effect on the dynamics of a system is small or non-existent. The critical node detection is also based on the residues; the node at which the effect of a pole on the dynamics is highest is defined as the critical node. In order to evaluate and check the efficiency of the created program, it is compared via examples with another existing commercial stability analysis tool (STAN tool). In this report, the newly created tool is proved to be as rigorous as STAN for detecting instabilities. Additionally, it is determined that the MIMO analysis is a very profitable addition to stability analysis, since it helps to eliminate possible problems of loss of controllability, observability and overmodeling.

The origin of the super-resolution via a nonlinear thin film

In laser applications, resolutions beyond the diffraction limit can be obtained with a thin film of strong optical nonlinear effect. The optical index of the silicon thin film is modified with the incident laser beam as a function of the local field intensity n(r) similar to E-2(r). For ultrathin films of thickness d << lambda the transmitted light through the film forms a profile of annular rings. Therefore, the device can be related to the realization of super-resolution with annular pupils. Theoretical analysis shows that the focused light spot appears significantly reduced in comparison with the diffraction limit that is determined by the laser wavelength and the numerical aperture of the converging lens. Analysis on the additional optical transfer function due to the thin film confirms that the resolving power is improved in the high spatial frequency region. (C) 2007 Published by Elsevier B.V.