A Perturbation Scheme for Passivity Verification and Enforcement of Parameterized Macromodels

This paper presents an algorithm for checking and enforcing passivity of behavioral reduced-order macromodels of LTI systems, whose frequency-domain (scattering) responses depend on external parameters. Such models, which are typically extracted from sampled input-output responses obtained from numerical solution of first-principle physical models, usually expressed as Partial Differential Equations, prove extremely useful in design flows, since they allow optimization, what-if or sensitivity analyses, and design centering. Starting from an implicit parameterization of both poles and residues of the model, as resulting from well-known model identification schemes based on the Generalized Sanathanan-Koerner iteration, we construct a parameter-dependent Skew-Hamiltonian/Hamiltonian matrix pencil. The iterative extraction of purely imaginary eigenvalues ot fhe pencil, combined with an adaptive sampling scheme in the parameter space, is able to identify all regions in the frequency-parameter plane where local passivity violations occur. Then, a singular value perturbation scheme is setup to iteratively correct the model coefficients, until all local passivity violations are eliminated. The final result is a corrected model, which is uniformly passive throughout the parameter range. Several numerical examples denomstrate the effectiveness of the proposed approach.Comment: Submitted to the IEEE Transactions on Components, Packaging and Manufacturing Technology on 13-Apr-201

Physics-based passivity-preserving parameterized model order reduction for PEEC circuit analysis

The decrease of integrated circuit feature size and the increase of operating frequencies require 3-D electromagnetic methods, such as the partial element equivalent circuit (PEEC) method, for the analysis and design of high-speed circuits. Very large systems of equations are often produced by 3-D electromagnetic methods, and model order reduction (MOR) methods have proven to be very effective in combating such high complexity. During the circuit synthesis of large-scale digital or analog applications, it is important to predict the response of the circuit under study as a function of design parameters such as geometrical and substrate features. Traditional MOR techniques perform order reduction only with respect to frequency, and therefore the computation of a new electromagnetic model and the corresponding reduced model are needed each time a design parameter is modified, reducing the CPU efficiency. Parameterized model order reduction (PMOR) methods become necessary to reduce large systems of equations with respect to frequency and other design parameters of the circuit, such as geometrical layout or substrate characteristics. We propose a novel PMOR technique applicable to PEEC analysis which is based on a parameterization process of matrices generated by the PEEC method and the projection subspace generated by a passivity-preserving MOR method. The proposed PMOR technique guarantees overall stability and passivity of parameterized reduced order models over a user-defined range of design parameter values. Pertinent numerical examples validate the proposed PMOR approach

A Two-Step Global Maximum Error Controller-Based TPWL MOR with POD Basis Vectors and Its Applications to MEMS

In our previous study, we have proposed a linearization point (LP) selection method based on a global maximum error controller for the trajectory piecewise-linear (TPWL) method. It has been demonstrated that this method has many advantages over other existing methods. In this paper, a more efficient version of this method is presented, which introduces a preliminary LP selection procedure and constructs projection matrix by the proper orthogonal decomposition (POD) method. Compared with the original method, the improved method takes much less time for extracting a reduced-order model (ROM) of similar quality and gets some other benefits (such as being easier to implement, having lower memory requirement, and enhanced flexibility). The effectiveness of the new method is fully demonstrated by a diode transmission line RLC circuit. And then, the method is applied to three more complicated microelectromechanical systems (MEMS) devices, which are a micromachined switch, an electrostatic micropump diaphragm, and a thermomechanical in-plane microactuator

Theoretical and practical aspects of linear and nonlinear model order reduction techniques

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2008.This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.Includes bibliographical references (p. 133-142).Model order reduction methods have proved to be an important technique for accelerating time-domain simulation in a variety of computer-aided design tools. In this study we present several new techniques for model reduction of the large-scale linear and nonlinear systems. First, we present a method for nonlinear system reduction based on a combination of the trajectory piecewise-linear (TPWL) method with truncated-balanced realizations (TBR). We analyze the stability characteristics of this combined method using perturbation theory. Second, we describe a linear reduction method that approximates TBR model reduction and takes advantage of sparsity of the system matrices or available accelerated solvers. This method is based on AISIAD (approximate implicit subspace iteration with alternate directions) and uses low-rank approximations of a system's gramians. This method is shown to be advantageous over the common approach of independently approximating the controllability and observability gramians, as such independent approximation methods can be inefficient when the gramians do not share a common dominant eigenspace. Third, we present a graph-based method for reduction of parameterized RC circuits. We prove that this method preserves stability and passivity of the models for nominal reduction. We present computational results for large collections of nominal and parameter-dependent circuits. Finally, we present a case study of model reduction applied to electroosmotic flow of a marker concentration pulse in a U-shaped microfluidic channel, where the marker flow in the channel is described by a three-dimensional convection-diffusion equation. First, we demonstrate the effectiveness of the modified AISIAD method in generating a low order models that correctly describe the dispersion of the marker in the linear case; that is, for the case of concentration-independent mobility and diffusion constants.(cont) Next, we describe several methods for nonlinear model reduction when the diffusion and mobility constants become concentration-dependent.by Dmitry Missiuro Vasilyev.Ph.D

Fast Simulation of Analog Circuit Blocks under Nonstationary Operating Conditions

This paper proposes a black-box behavioral modeling framework for analog circuit blocks operating under small-signal conditions around non-stationary operating points. Such variations may be induced either by changes in the loading conditions or by event-driven updates of the operating point for system performance optimization, e.g., to reduce power consumption. An extension of existing data-driven parameterized reduced-order modeling techniques is proposed that considers the time-varying bias components of the port signals as non-stationary parameters. These components are extracted at runtime by a lowpass filter and used to instantaneously update the matrices of the reduced-order state-space model realized as a SPICE netlist. Our main result is a formal proof of quadratic stability of such Linear Parameter Varying (LPV) models, enabled by imposing a specific model structure and representing the transfer function in a basis of positive functions whose elements constitute a partition of unity. The proposed quadratic stability conditions are easily enforced through a finite set of small-size Linear Matrix Inequalities (LMI), used as constraints during model construction. Numerical results on various circuit blocks including voltage regulators confirm that our approach not only ensures the model stability, but also provides speedup in runtime up to 2 orders of magnitude with respect to full transistor-level circuits

Tensor Computation: A New Framework for High-Dimensional Problems in EDA

Many critical EDA problems suffer from the curse of dimensionality, i.e. the very fast-scaling computational burden produced by large number of parameters and/or unknown variables. This phenomenon may be caused by multiple spatial or temporal factors (e.g. 3-D field solvers discretizations and multi-rate circuit simulation), nonlinearity of devices and circuits, large number of design or optimization parameters (e.g. full-chip routing/placement and circuit sizing), or extensive process variations (e.g. variability/reliability analysis and design for manufacturability). The computational challenges generated by such high dimensional problems are generally hard to handle efficiently with traditional EDA core algorithms that are based on matrix and vector computation. This paper presents "tensor computation" as an alternative general framework for the development of efficient EDA algorithms and tools. A tensor is a high-dimensional generalization of a matrix and a vector, and is a natural choice for both storing and solving efficiently high-dimensional EDA problems. This paper gives a basic tutorial on tensors, demonstrates some recent examples of EDA applications (e.g., nonlinear circuit modeling and high-dimensional uncertainty quantification), and suggests further open EDA problems where the use of tensor computation could be of advantage.Comment: 14 figures. Accepted by IEEE Trans. CAD of Integrated Circuits and System

System- and Data-Driven Methods and Algorithms

An increasing complexity of models used to predict real-world systems leads to the need for algorithms to replace complex models with far simpler ones, while preserving the accuracy of the predictions. This two-volume handbook covers methods as well as applications. This first volume focuses on real-time control theory, data assimilation, real-time visualization, high-dimensional state spaces and interaction of different reduction techniques

Regelungstheorie

The workshop “Regelungstheorie” (control theory) covered a broad variety of topics that were either concerned with fundamental mathematical aspects of control or with its strong impact in various fields of engineering

Stability-preserving model reduction for linear and nonlinear systems arising in analog circuit applications

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2010.Cataloged from PDF version of thesis.Includes bibliographical references (p. 221-229).Despite the increasing presence of RF and analog components in personal wireless electronics, such as mobile communication devices, the automated design and optimization of such systems is still an extremely challenging task. This is primarily due to the presence of both parasitic elements and highly nonlinear elements, which makes simulation computationally expensive and slow. The ability to generate parameterized reduced order models of analog systems could serve as a first step toward the automatic and accurate characterization of geometrically complex components and subcircuits, eventually enabling their synthesis and optimization. This thesis presents techniques for reduced order modeling of linear and nonlinear systems arising in analog applications. Emphasis is placed on developing techniques capable of preserving important system properties, such as stability, and parameter dependence in the reduced models. The first technique is a projection-based model reduction approach for linear systems aimed at generating stable and passive models from large linear systems described by indefinite, and possibly even mildly unstable, matrices. For such systems, existing techniques are either prohibitively computationally expensive or incapable of guaranteeing stability and passivity. By forcing the reduced model to be described by definite matrices, we are able to derive a pair of stability constraints that are linear in terms of projection matrices.(cont.) These constraints can be used to formulate a semidefinite optimization problem whose solution is an optimal stabilizing projection framework. The second technique is a projection-based model reduction approach for highly nonlinear systems that is based on the trajectory piecewise linear (TPWL) method. Enforcing stability in nonlinear reduced models is an extremely difficult task that is typically ignored in most existing techniques. Our approach utilizes a new nonlinear projection in order to ensure stability in each of the local models used to describe the nonlinear reduced model. The TPWL approach is also extended to handle parameterized models, and a sensitivity-based training system is presented that allows us to efficiently select inputs and parameter values for training. Lastly, we present a system identification approach to model reduction for both linear and nonlinear systems. This approach utilizes given time-domain data, such as input/output samples generated from transient simulation, in order to identify a compact stable model that best fits the given data. Our procedure is based on minimization of a quantity referred to as the 'robust equation error', which, provided the model is incrementally stable, serves as up upper bound for a measure of the accuracy of the identified model termed 'linearized output error'. Minimization of this bound, subject to an incremental stability constraint, can be cast as a semidefinite optimization problem.by Bradley Neil Bond.Ph.D