Field solver technologies for variation-aware interconnect parasitic extraction

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. 207-213).Advances in integrated circuit manufacturing technologies have enabled high density onchip integration by constantly scaling down the device and interconnect feature size. As a consequence of the ongoing technology scaling (from 45nm to 32nm, 22nm and beyond), geometrical variabilities induced by the uncertainties in the manufacturing processes are becoming more significant. Indeed, the dimensions and shapes of the manufactured devices and interconnect structures may vary by up to 40% from their design intent. The effect of such variabilities on the electrical characteristics of both devices and interconnects must be accurately evaluated and accounted for during the design phase. In the last few years, there have been several attempts to develop variation-aware extraction algorithms, i.e. algorithms that evaluate the effect of geometrical variabilities on the electrical characteristics of devices and interconnects. However, most algorithms remain computationally very expensive. In this thesis the focus is on variation-aware interconnect parasitic extraction. In the first part of the thesis several discretization-based variation-aware solver techniques are developed. The first technique is a stochastic model reduction algorithm (SMOR) The SMOR guarantees that the statistical moments computed from the reduced model are the same as those of the full model. The SMOR works best for problems in which the desired electrical property is contained in an easily defined subspace.(cont.) The second technique is the combined Neumann Hermite expansion (CNHE). The CNHE combines the advantages of both the standard Neumann expansion and the standard stochastic Galerkin method to produce a very efficient extraction algorithm. The CNHE works best in problems for which the desired electrical property (e.g. impedance) is accurately expanded in terms of a low order multivariate Hermite expansion. The third technique is the stochastic dominant singular vectors method (SDSV). The SDSV uses stochastic optimization in order to sequentially determine an optimal reduced subspace, in which the solution can be accurately represented. The SDSV works best for large dimensional problems, since its complexity is almost independent of the size of the parameter space. In the second part of the thesis, several novel discretization-free variation aware extraction techniques for both resistance and capacitance extraction are developed. First we present a variation-aware floating random walk (FRW) to extract the capacitance/resistance in the presence of non-topological (edge-defined) variations. The complexity of such algorithm is almost independent of the number of varying parameters. Then we introduce the Hierarchical FRW to extract the capacitance/resistance of a very large number of topologically different structures, which are all constructed from the same set of building blocks. The complexity of such algorithm is almost independent of the total number of structures. All the proposed techniques are applied to a variety of examples, showing orders of magnitude reduction in the computational time compared to the standard approaches. In addition, we solve very large dimensional examples that are intractable when using standard approaches.by Tarek Ali El-Moselhy.Ph.D

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

A high-performance boundary element method and its applications in engineering

As a semi-numerical and semi-analytical method, owing to the inherent advantage, of boundary-only discretisation, the boundary element method (BEM) has been widely applied to problems with complicated geometries, stress concentration problems, infinite domain problems, and many others. However, domain integrals and non-symmetrical and dense matrix systems are two obstacles for BEM which have hindered the its further development and application. This thesis is aimed at proposing a high-performance BEM to tackle the above two drawbacks and broaden the application scope of BEM. In this thesis, a detailed introduction to the traditional BEM is given and several popular algorithms are introduced or proposed to enhance the performance of BEM. Numerical examples in heat conduction analysis, thermoelastic analysis and thermoelastic fracture problems are performed to assess the efficiency and correction of the algorithms. In addition, necessary theoretical derivations are embraced for establishing novel boundary integral equations (BIEs) for specific engineering problems. The following three parts are the main content of this thesis. (1) The first part (Part II consisting of two chapters) is aimed at heat conduction analysis by BEM. The coefficient matrix of equations formed by BEM in solving problems is fully-populated which occupy large computer memory. To deal with that, the fast multipole method (FMM) is introduced to energize the line integration boundary element method (LIBEM) to performs better in efficiency. In addition, to compute domain integrals with known or unknown integrand functions which are caused by heat sources or heterogeneity, a novel BEM, the adaptive orthogonal interpolation moving least squares (AOIMLS) method enhanced LIBEM, which also inherits the advantage of boundary-only discretisation, is proposed. Unlike LIBEM, which is an accurate and stable method for computing domain integrals, but only works when the mathematical expression of integral function in domain integrals is known, the AOIMLS enhanced LIBEM can compute domain integrals with known or unknown integral functions, which ensures all the nonlinear and nonhomogeneous problems can be solved without domain discretisation. In addition, the AOIMLS can adaptively avoid singular or ill-conditioned moment matrices, thus ensuring the stability of the calculation results. (2) In the second part (Part III consisting of four chapters), the thermoelastic problems and fracture problems are the main objectives. Due to considering thermal loads, domain integrals appear in the BIEs of the thermoelastic problems, and the expression of integrand functions is known or not depending on the temperature distribution given or not, the AOIMLS enhanced LIBEM is introduced to conduct thermoelasticity analysis thereby. Besides, a series of novel unified boundary integral equations based on BEM and DDM are derived for solving fracture problems and thermoelastic fracture problems in finite and infinite domains. Two sets of unified BIEs are derived for fracture problems in finite and infinite domains based on the direct BEM and DDM respectively, which can provide accurate and stable results. Another two sets of BIEs are addressed by employing indirect BEM and DDM, which cannot ensure a stable result, thereby a modified indirect BEM is proposed which performs much more stable. Moreover, a set of novel BIEs based on the direct BEM and DDM for cracked domains under thermal stress is proposed. (3) In the third part (Part IV consisting of one chapter), a high-efficiency combined BEM and discrete element method (DEM) is proposed to compute the inner stress distribution and particle breakage of particle assemblies based on the solution mapping scheme. For the stress field computation of particles with similar geometry, a template particle is used as the representative particle, so that only the related coefficient matrices of one template particle in the local coordinate system are needed to be calculated, while the coefficient matrices of the other particles, can be obtained by mapping between the local and global coordinate systems. Thus, the combined BEM and DEM is much more effective when modelling a large-scale particle system with a small number of distinct possible particle shapes. Furthermore, with the help of the Hoek-Brown criterion, the possible cracks or breakage paths of a particle can be obtained

Modeling EMI Resulting from a Signal Via Transition Through Power/Ground Layers

Signal transitioning through layers on vias are very common in multi-layer printed circuit board (PCB) design. For a signal via transitioning through the internal power and ground planes, the return current must switch from one reference plane to another reference plane. The discontinuity of the return current at the via excites the power and ground planes, and results in noise on the power bus that can lead to signal integrity, as well as EMI problems. Numerical methods, such as the finite-difference time-domain (FDTD), Moment of Methods (MoM), and partial element equivalent circuit (PEEC) method, were employed herein to study this problem. The modeled results are supported by measurements. In addition, a common EMI mitigation approach of adding a decoupling capacitor was investigated with the FDTD method

Proceedings of the UCLA International Conference on Radiation and Remote Probing of the Atmosphere

Various articles are presented on multiple scattering problems and radiative transfer in the atmosphere. Particle size distribution and molecular absorption are also discussed

Modelling of radio wave propagation using Finite Element Analysis.

Fourth generation (4G) wireless communication systems are intended to support high data rates which requires careful and accurate modelling of the radio environment. In this thesis, for the first time finite clement based accurate and computationally efficient models of wave propagation in different outdoor and indoor environments has been developed. Three different environments were considered: the troposphere, vegetation and tunnels and wave propagation in these environments were modelled using finite element analysis. Use of finite elements in wave propagation modelling is a novel idea although many propagation models and approaches were used in past. Coverage diagrams, path loss contours and power levels were calculated using developed models in the troposphere, vegetation and tunnels. Results obtained were compared with commercially available software Advanced Refractive Effects Prediction Software (AREPS) to validate the accuracy of the developed approach and it is shown that results were accurate with an accuracy of 3dB. The developed models were very flexible in handling complex geometries and similar analysis can be easily extended to other environments. A fully vectored finite element base propagation model was developed for straight and curved tunnels. An optimum range of values of different electrical parameters for tunnels of different shapes has been derived. The thesis delivered a novel approach to modelling radio channels that provided a fast and accurate solution of radio wave propagation in realistic environments. The results of this thesis will have a great impact in modelling and characterisation of future wireless communication systems

Reconstruction algorithms for multispectral diffraction imaging

Thesis (Ph.D.)--Boston UniversityIn conventional Computed Tomography (CT) systems, a single X-ray source spectrum is used to radiate an object and the total transmitted intensity is measured to construct the spatial linear attenuation coefficient (LAC) distribution. Such scalar information is adequate for visualization of interior physical structures, but additional dimensions would be useful to characterize the nature of the structures. By imaging using broadband radiation and collecting energy-sensitive measurement information, one can generate images of additional energy-dependent properties that can be used to characterize the nature of specific areas in the object of interest. In this thesis, we explore novel imaging modalities that use broadband sources and energy-sensitive detection to generate images of energy-dependent properties of a region, with the objective of providing high quality information for material component identification. We explore two classes of imaging problems: 1) excitation using broad spectrum sub-millimeter radiation in the Terahertz regime and measure- ment of the diffracted Terahertz (THz) field to construct the spatial distribution of complex refractive index at multiple frequencies; 2) excitation using broad spectrum X-ray sources and measurement of coherent scatter radiation to image the spatial distribution of coherent-scatter form factors. For these modalities, we extend approaches developed for multimodal imaging and propose new reconstruction algorithms that impose regularization structure such as common object boundaries across reconstructed regions at different frequencies. We also explore reconstruction techniques that incorporate prior knowledge in the form of spectral parametrization, sparse representations over redundant dictionaries and explore the advantage and disadvantages of these techniques in terms of image quality and potential for accurate material characterization. We use the proposed reconstruction techniques to explore alternative architectures with reduced scanning time and increased signal-to-noise ratio, including THz diffraction tomography, limited angle X-ray diffraction tomography and the use of coded aperture masks. Numerical experiments and Monte Carlo simulations were conducted to compare performances of the developed methods, and validate the studied architectures as viable options for imaging of energy-dependent properties