Thesis (PhD (Electrical and Electronic Engineering))--University of Stellenbosch, 2011. / ENGLISH ABSTRACT: This work considers the acceleration of matrix-based computational electromagnetic (CEM) techniques using graphics processing units (GPUs). These massively parallel processors have gained much support since late 2006, with software tools such as CUDA and OpenCL greatly simplifying the process of harnessing the computational power of these devices. As with any advances in computation, the use of these devices enables the modelling of more complex problems, which in turn should give rise to better solutions to a number of global challenges faced at present. For the purpose of this dissertation, CUDA is used in an investigation of the acceleration of two methods in CEM that are used to tackle a variety of problems. The first of these is the Method of Moments (MOM) which is typically used to model radiation and scattering problems, with the latter begin considered here. For the CUDA acceleration of the MOM presented here, the assembly and subsequent solution of the matrix equation associated with the method are considered. This is done for both single and double precision oating point matrices. For the solution of the matrix equation, general dense linear algebra techniques are used, which allow for the use of a vast expanse of existing knowledge on the subject. This also means that implementations developed here along with the results presented are immediately applicable to the same wide array of applications where these methods are employed. Both the assembly and solution of the matrix equation implementations presented result in signi cant speedups over multi-core CPU implementations, with speedups of up to 300x and 10x, respectively, being measured. The implementations presented also overcome one of the major limitations in the use of GPUs as accelerators (that of limited memory capacity) with problems up to 16 times larger than would normally be possible being solved. The second matrix-based technique considered is the Finite Element Method (FEM), which allows for the accurate modelling of complex geometric structures including non-uniform dielectric and magnetic properties of materials, and is particularly well suited to handling bounded structures such as waveguide. In this work the CUDA acceleration of the cutoff and dispersion analysis of three waveguide configurations is presented. The modelling of these problems using an open-source software package, FEniCS, is also discussed. Once again, the problem can be approached from a linear algebra perspective, with the formulation in this case resulting in a generalised eigenvalue (GEV) problem. For the problems considered, a total solution speedup of up to 7x is measured for the solution of the generalised eigenvalue problem, with up to 22x being attained for the solution of the standard eigenvalue problem that forms part of the GEV problem. / AFRIKAANSE OPSOMMING: In hierdie werkstuk word die versnelling van matriksmetodes in numeriese elektromagnetika (NEM) deur die gebruik van grafiese verwerkingseenhede (GVEe) oorweeg. Die gebruik van hierdie verwerkingseenhede is aansienlik vergemaklik in 2006 deur sagteware pakette soos CUDA en OpenCL. Hierdie toestelle, soos ander verbeterings in verwerkings vermoe, maak dit moontlik om meer komplekse probleme op te los. Hierdie stel wetenskaplikes weer in staat om globale uitdagings beter aan te pak. In hierdie proefskrif word CUDA gebruik om ondersoek in te stel na die versnelling van twee metodes in NEM, naamlik die Moment Metode (MOM) en die Eindige Element Metode (EEM). Die MOM word tipies gebruik om stralings- en weerkaatsingsprobleme op te los. Hier word slegs na die weerkaatsingsprobleme gekyk. CUDA word gebruik om die opstel van die MOM matriks en ook die daaropvolgende oplossing van die matriksvergelyking wat met die metode gepaard gaan te bespoedig. Algemene digte lineere algebra tegnieke word benut om die matriksvergelykings op te los. Dit stel die magdom bestaande kennis in die vagebied beskikbaar vir die oplossing, en gee ook aanleiding daartoe dat enige implementasies wat ontwikkel word en resultate wat verkry word ook betrekking het tot 'n wye verskeidenheid probleme wat die lineere algebra metodes gebruik. Daar is gevind dat beide die opstelling van die matriks en die oplossing van die matriksvergelyking aansienlik vinniger is as veelverwerker SVE implementasies. 'n Verselling van tot 300x en 10x onderkeidelik is gemeet vir die opstel en oplos fases. Die hoeveelheid geheue beskikbaar tot die GVE is een van die belangrike beperkinge vir die gebruik van GVEe vir groot probleme. Hierdie beperking word hierin oorkom en probleme wat selfs 16 keer groter is as die GVE se beskikbare geheue word geakkommodeer en suksesvol opgelos. Die Eindige Element Metode word op sy beurt gebruik om komplekse geometriee asook nieuniforme materiaaleienskappe te modelleer. Die EEM is ook baie geskik om begrensde strukture soos golfgeleiers te hanteer. Hier word CUDA gebruik of om die afsny- en dispersieanalise van drie gol eierkonfigurasies te versnel. Die implementasie van hierdie probleme word gedoen deur 'n versameling oopbronkode wat bekend staan as FEniCS, wat ook hierin bespreek word. Die probleme wat ontstaan in die EEM kan weereens vanaf 'n lineere algebra uitganspunt benader word. In hierdie geval lei die formulering tot 'n algemene eiewaardeprobleem. Vir die gol eier probleme wat ondersoek word is gevind dat die algemene eiewaardeprobleem met tot 7x versnel word. Die standaard eiewaardeprobleem wat 'n stap is in die oplossing van die algemene eiewaardeprobleem is met tot 22x versnel.

This dissertation presents techniques for high-order simulation of electromagnetic fields, particularly for problems involving ships with ferromagnetic hulls and active corrosion-protection systems. A set of numerically constrained hexahedral basis functions for volume integral equation discretization is presented in a method-of-moments context. Test simulations demonstrate the accuracy achievable with these functions as well as the improvement brought about in system conditioning when compared to other basis sets. A general method for converting between a locally-corrected Nyström discretization of an integral equation and a method-of-moments discretization is presented next. Several problems involving conducting and magnetic-conducting materials are solved to verify the accuracy of the method and to illustrate both the reduction in number of unknowns and the effect of the numerically constrained bases on the conditioning of the converted matrix. Finally, a surface integral equation derived from Laplace’s equation is discretized using the locally-corrected Nyström method in order to calculate the electric fields created by impressed-current corrosion protection systems. An iterative technique is presented for handling nonlinear boundary conditions. In addition we examine different approaches for calculating the magnetic field radiated by the corrosion protection system. Numerical tests show the accuracy achievable by higher-order discretizations, validate the iterative technique presented. Various methods for magnetic field calculation are also applied to basic test cases.

<div>The design of advanced integrated circuits (ICs) and systems calls for multiphysics and large-scale modeling and simulation methods. On the one hand, novel devices and materials are emerging in next-generation IC technology, which requires multiphysics modeling and simulation. On the other hand, the ever-increasing complexity of ICs requires more efficient numerical solvers.</div><div><br></div><div>In this work, we propose a multiphysics modeling and simulation algorithm to co-simulate Maxwell's equations, dispersion relation of materials, and Boltzmann equation to characterize emerging new devices in IC technology such as Cu-Graphene (Cu-G) hybrid nano-interconnects. We also develop an unconditionally stable time marching scheme to remove the dependence of time step on space step for an efficient simulation of the multiscaled and multiphysics system. Extensive numerical experiments and comparisons with measurements have validated the accuracy and efficiency of the proposed algorithm. Compared to simplified steady-state-models based analysis, a significant difference is observed when the frequency is high or/and the dimension of the Cu-G structure is small, which necessitates our proposed multiphysics modeling and simulation for the design of advanced Cu-G interconnects. </div><div><br></div><div>To address the large-scale simulation challenge, we develop a new split-field domain-decomposition algorithm amenable for parallelization for solving Maxwell’s equations, which minimizes the communication between subdomains, while having a fast convergence of the global solution. Meanwhile, the algorithm is unconditionally stable in time domain. In this algorithm, unlike prevailing domain decomposition methods that treat the interface unknown as a whole and let it be shared across subdomains, we partition the interface unknown into multiple components, and solve each of them from one subdomain. In this way, we transform the original coupled system to fully decoupled subsystems to solve. Only one addition (communication) of the interface unknown needs to be performed after the computation in each subdomain is finished at each time step. More importantly, the algorithm has a fast convergence and permits the use of a large time step irrespective of space step. Numerical experiments on large-scale on-chip and package layout analysis have demonstrated the capability of the new domain decomposition algorithm. </div><div><br></div><div>To tackle the challenge of efficient simulation of irregular structures, in the last part of the thesis, we develop a method for the stability analysis of unsymmetrical numerical systems in time domain. An unsymmetrical system is traditionally avoided in numerical formulation since a traditional explicit simulation is absolutely unstable, and how to control the stability is unknown. However, an unsymmetrical system is frequently encountered in modeling and simulating of unstructured meshes and nonreciprocal electromagnetic and circuit devices. In our method, we reduce stability analysis of a large system into the analysis of dissembled single element, therefore provides a feasible way to control the stability of large-scale systems regardless of whether the system is symmetrical or unsymmetrical. We then apply the proposed method to prove and control the stability of an unsymmetrical matrix-free method that solves Maxwell’s equations in general unstructured meshes while not requiring a matrix solution.<br></div><div><br></div>

<p dir="ltr">This work introduces an innovative algorithm for co-simulating time-dependent Drift Diffusion (DD) equations with Maxwell\textquotesingle s equations to characterize semiconductor devices. Traditionally, the DD equations, derived from the Boltzmann transport equations, are used alongside Poisson\textquotesingle s equation to model electronic carriers in semiconductors. While DD equations coupled with Poisson\textquotesingle s equation underpin commercial TCAD software for micron-scale device simulation, they are limited by electrostatic assumptions and fail to capture time dependent high-frequency effects. Maxwell\textquotesingle s equations are fundamental to classical electrodynamics, enabling the prediction of electrical performance across frequency range crucial to advanced device fabrication and design. However, their integration with DD equations has not been studied thoroughly. The proposed method advances current simulation techniques by introducing a new broadband patch-based method to solve time-domain 3-D Maxwell\textquotesingle s equations and integrating it with the solution of DD equations. This technique is free of the low-frequency breakdown issues prevalent in conventional full-wave simulations. Meanwhile, it enables large-scale simulations with reduced computational complexity. This work extends the simulation to encompass the complete device, including metal contacts and interconnects. Thus, it captures the entire electromagnetic behavior, which is especially critical in electrically larger systems and high-frequency scenarios. The electromagnetic interactions of the device with its contacts and interconnects are investigated, providing insights into performance at the chip level. Validation through numerical experiments and comparison with results from commercial TCAD tools confirm the effectiveness of the proposed method. </p>

Dealing with the electromagnetic issue might bring a sort of discontinuous and nondifferentiable regions. Thus, it is of great interest to implement an appropriate optimisation approach, which can preserve the computational resources and come up with a global optimum. While not being trapped in local optima, as well as the feasibility to overcome some other matters such as nonlinear and phenomena of discontinuous with a large number of variables. Problems such as lengthy computation time, constraints put forward for antenna requirements and demand for large computer memory, are very common in the analysis due to the increased interests in tackling high-scale, more complex and higher-dimensional problems. On the other side, demands for even more accurate results always expand constantly. In the context of this statement, it is very important to find out how the recently developed optimization roles can contribute to the solution of the aforementioned problems. Thereafter, the key goals of this work are to model, study and design low profile antennas for wireless and mobile communications applications using optimization process over a computational electromagnetics numerical solution. The numerical solution method could be performed over one or hybrid methods subjective to the design antenna requirements and its environment. Firstly, the thesis presents the design and modelling concept of small uni-planer Ultra- Wideband antenna. The fitness functions and the geometrical antenna elements required for such design are considered. Two antennas are designed, implemented and measured. The computed and measured outcomes are found in reasonable agreement. Secondly, the work is also addressed on how the resonance modes of microstrip patches could be performed using the method of Moments. Results have been shown on how the modes could be adjusted using MoM. Finally, the design implications of balanced structure for mobile handsets covering LTE standards 698-748 MHz and 2500-2690 MHz are explored through using firefly algorithm method. The optimised balanced antenna exhibits reasonable matching performance including near-omnidirectional radiations over the dual desirable operating bands with reduced EMF, which leads to a great immunity improvement towards the hand-held. / General Secretariat of Education and Scientific Research Libya

No description available.

[pt] Este trabalho apresenta uma abordagem inovadora para calcular autopares de operadores diferenciais (OD), utilizando programação genética (PG) e computação simbólica. Na literatura atual, o Método dos Elementos Finitos (MEF) e o Método das Diferenças Finitas (MDF) são os mais utilizados. Tais métodos usam discretização para converter o OD em uma matriz finita e, por isso, apresentam limitações como perda de acurácia e presença de soluções espúrias. Além disso, se o domínio do OD fosse alterado, os autopares precisariam ser calculados novamente, pois a representação matricial do operador depende dos parâmetros do problema. Nesse contexto, este trabalho propõe evoluir autofunções analiticamente usando PG, sem discretização do domínio. Com isso, é possível incorporar parâmetros, o que torna a solução obtida válida para uma classe de problemas. Este texto descreve o modelo para OD normais, aplicando conceitos de indivíduos multi-árvore e diferenciação simbólica. O modelo evolui auto-funções e, a partir delas, calcula os autovalores empregando a razão de Rayleigh. Experimentos baseados em aplicações da Física mostram que a técnica proposta é capaz de encontrar as autofunções analíticas com a acurácia igual ou melhor que as técnicas numéricas supracitadas. Finalmente, a técnica proposta é aplicada ao problema de propagação de ondas eletromagnéticas em poços de petróleo em ULF e UHF. As soluções analíticas são dadas em função do diâmetro e do percentual de incrustações no poço. Os resultados mostram que, para um conjunto suficientemente grande de valores distintos dos parâmetros, a técnica apresenta tempo de execução inferior às técnicas clássicas, mantendo a acurácia destas. / [en] This work presents an innovative approach to calculate the eigenpairs of linear differential operators (LDO), employing genetic programming (GP) and symbolic computation. In the current literature, the Finite Element Method (FEM) and the Finite Difference Method (FDM) are more commonly used. Such methods use discretization to convert the LDO to a finite matrix, therefore causing loss of accuracy and presence of spurious solutions. Additionally, if the domain of the LDO was changed, the eigenpairs would need to be recalculated, since the matrix representation of the LDO depends on the parameters of the problem. In this context, this work proposes to evolve eigenfunctions analytically using GP, without domain discretization. Hence, it is possible to incorporate the parameter, which makes a obtained solution valid for a class of problems. This text describes the model for normal LDO, applying concepts of multi-tree individuals and symbolic differentiation. The presented model evolves eigenfunctions and, then, calculates the eigenvalues using the Rayleigh quotient. Experiments based on Physics problems show that the proposed technique is able to find the analytical eigenfunctions with the same accuracy of the numerical techniques mentioned above. Finally, the proposed technique is applied to the problem of propagation of electromagnetic waves in oil wells in ULF and UHF. The analytical solutions are given as a function of the diameter and percentage of CaCO in the well. The results show that, for a sufficiently large set of distinct values of the parameters, the technique presents execution time inferior to the FEM, while maintaining its accuracy.

Cette étude s’inscrit dans le domaine de la simulation électromagnétique des problèmes de grande taille tels que la diffraction d’ondes planes par de larges plateformes et le rayonnement d’antennes aéroportées. Elle consiste à développer une méthode combinant décomposition en sous-domaines et compression hiérarchique des équations intégrales de frontière. Pour cela, nous rappelons dans un premier temps les points importants de la méthode des équations intégrales de frontière et de leur compression hiérarchique par l’algorithme ACA (Adaptive Cross Approximation). Ensuite, nous présentons la formulation IE-DDM (Integral Equations – Domain Decomposition Method) obtenue à partir d’une représentation intégrale des sous-domaines. Les matrices résultant de la discrétisation de cette formulation sont stockées au format H-matrice (matricehiérarchique). Un solveur spécialement adapté à la résolution de la formulation IE-DDM et à sa représentation hiérarchique a été conçu. Cette étude met en évidence l’efficacité de la décomposition en sous-domaines en tant que préconditionneur des équations intégrales. De plus, la méthode développée est rapide pour la résolution des problèmes à incidences multiples ainsi que la résolution des problèmes basses fréquences / This thesis is about the electromagnetic simulation of large scale problems as the wave scattering from aircrafts and the airborne antennas radiation. It consists in the development of a method combining domain decomposition and hierarchical compression of the surface integral equations. First, we remind the principles of the boundary element method and the hierarchical representation of the surface integral equations with the Adaptive Cross Approximation algorithm. Then, we present the IE-DDM formulation obtained from a sub-domain integral representation. The matrices resulting of the discretization of the formulation are stored in the H-matrix format. A solver especially fitted with the hierarchical representation of the IE-DDM formulation has been developed. This study highlights the efficiency of the sub-domain decomposition as a preconditioner of the integral equations. Moreover, the method is fast for the resolution of multiple incidences and the resolution of low frequencies problems

La propagation à longue distance est un problème majeur dans les télécommunications, la navigation et la surveillance. L'objectif de cette thèse est de développer une méthode rapide pour simuler la propagation des ondes dans une atmosphère en 2D et 3D. Dans ce travail, deux contributions principales vers cet objectif sont obtenues. Tout d'abord, des méthodes auto-cohérentes,c'est-à-dire basées sur une théorie discrète de l'électromagnétisme, sont développées en 2D et 3D. Ensuite, une méthode rapide 2D basée sur les ondelettes est proposée. Pour simuler la propagation d'ondes électromagnétiques dans une atmosphère 2D, la méthode split-step Fourier (SSF) est largement utilisée. Le calcul est effectué itérativement en distances en tenant compte d'une réfractivité variable, du relief et des caractéristiques du sol. À chaque pas, le signal est transformé du domaine spatial au domaine spectral. La méthode des écrans de phase est appliquée pour modéliser les effets de réfraction. D'autre part, pour modéliser un sol impédant, la transformée mixte de Fourier discrète (SSF-DMFT) est utilisée. Le concept de la théorie électromagnétique auto-cohérente implique que l'utilisation d'équations de Maxwell discrètes pour la simulation numérique évite les solutions parasites. Dans la méthode couramment utilisée SSF-DMFT, la transformée spectrale est basée sur la condition aux limites d'impédance discrète, alors que le propagateur provient de l'équation continue. Pour pallier cette incohérence, une méthode auto-cohérente est proposée, notée la DSSF-DMFT. La formulation est dérivée des équations discrètes pour obtenir l'auto-cohérence. Des tests numériques montrent que SSF-DMFT présente des oscillations parasites dans certaines conditions de simulation, tandis que DSSF-DMFT reste précis. En effet, l'auto-cohérence empêche certaines instabilités numériques. Pour simuler la propagation dans des environnements en 3D, les méthodes précédentes doivent être étendues en 3D. Tout d'abord, la 3D-SSF est présentée comme une extension naturelle de la SSF. Ensuite, la 3D-DSSF est dérivée d'équations discrètes. Pour considérer un sol impédant, la 3D-DSSF-DMFT est développée conduisant à de nouvelles expressions pour les propagateurs. Ces méthodes sont testées dans plusieurs configurations incluant un profil de réfractivité extrait de mesures. Les résultats montrent une grande précision et une capacité à prendre en compte les effets latéraux. Cependant, pour la propagation dans de grand domaines, les ressources nécessaires en temps et en mémoire deviennent la préoccupation principale. Pour alléger la charge de calcul, une méthode split-step en ondelettes (SSW) est proposée en 2D comme une méthode alternative à SSF. Elle est basée sur la transformée rapide en ondelettes dont la complexité est faible et qui permet de compresser les champs. La propagation est réalisée à partir d'une combinaison linéaire d'ondelettes propagées individuellement. La compression est appliquée pour augmenter l'efficacité. Afin de considérer la réflexion sur le sol, une nouvelle méthode de source image locale dédiée à la propagation des ondelettes est proposée. Les tests numériques montrent que la SSW a une efficacité de calcul plus élevée que la SSF tout en gardant une bonne précision. / The long-range propagation of electromagnetic waves is a major issue in telecommunication, navigation, and surveillance. The objective of this Ph.D. thesis is to develop fast and accurate modeling methods for the tropospheric propagation in 2D and 3D. In this work, two main contributions towards this objective are achieved. Firstly, self-consistent methods, i.e. based on the discrete electromagnetic theory, are developed in 2D and 3D. Secondly, a fast wavelet-based 2D method is proposed. For simulating the electromagnetic wave propagation in a 2D atmosphere, the split-step Fourier method (SSF) is widely used. The computation is performed marching on in distances taking into account a variable refractivity, an irregular relief, and the electric characteristics of the ground. At each step, the signal is transformed from the spatial to the spectral domains. The phase screens method is applied to model refraction. Besides, to model an impedance ground, the discrete mixed Fourier transform (SSF-DMFT) is used. The concept of the self-consistent electromagnetic theory implies that the use of discrete Maxwell equations for numerical simulations does not lead to spurious solutions. In the widely used SSF-DMFT, the spectral transform is based on the discrete impedance boundary condition, while the propagator is derived from the continuous equation. To overcome this inconsistency, a discrete formulation of SSF-DMFT is proposed, denoted as DSSF-DMFT. The spectral transform and propagator are both derived from the discrete equations to achieve self-consistency. Numerical tests show that SSF-DMFT has spurious oscillations in certain simulation conditions, whereas DSSF-DMFT remains accurate. Indeed, the self-consistency prevents from numerical instabilities. To simulate the propagation in 3D environments, the previous methods are extended to 3D. First, 3D-SSF is presented as a natural extension of SSF. Then, 3D-DSSF is derived from discrete equations. To consider an impedance ground, 3D-DSSF-DMFT is developed leading to new expressions for the propagators. These methods are tested for several configurations, including a refractivity profile extracted from measurements. Results show that they have a high accuracy. They notably consider lateral effects. However, for the propagation in a large computation domain, time and memory occupations become the main concern. To improve the computation burden, a split-step wavelet method (SSW) is proposed in 2D as an alternative to SSF. It is based on the fast wavelet transform, which complexity is weak and which allows for data compression. The propagation is performed by means of a linear combination of wavelets that are individually propagated. Data compression is applied to increase the efficiency. A new local image source method dedicated to wavelet propagation is proposed to consider the ground reflection. Numerical tests show that this method has a higher computational efficiency than SSF while keeping a good accuracy.

With the rapid increase in available compute power and memory, and bolstered by the advent of efficient formulations and algorithms, the role of 3D full-wave computational methods for accurate modelling of complex electromagnetic (EM) structures has gained in significance. The range of problems includes Radar Cross Section (RCS) computation, analysis and design of antennas and passive microwave circuits, bio-medical non-invasive detection and therapeutics, energy harvesting etc. Further, with the rapid advances in technology trends like System-in-Package (SiP) and System-on-Chip (SoC), the fidelity of chip-to-chip communication and package-board electrical performance parameters like signal integrity (SI), power integrity (PI), electromagnetic interference (EMI) are becoming increasingly critical. Rising pin-counts to satisfy functionality requirements and decreasing layer-counts to maintain cost-effectiveness necessitates 3D full wave electromagnetic solution for accurate system modelling. Method of Moments (MoM) is one such widely used computational technique to solve a 3D electromagnetic problem with full-wave accuracy. Due to lesser number of mesh elements or discretization on the geometry, MoM has an advantage of a smaller matrix size. However, due to Green's Function interactions, the MoM matrix is dense and its solution presents a time and memory challenge. The thesis focuses on formulation and development of novel techniques that aid in fast MoM based electromagnetic solutions. With the recent paradigm shift in computer hardware architectures transitioning from single-core microprocessors to multi-core systems, it is of prime importance to parallelize the serial electromagnetic formulations in order to leverage maximum computational benefits. Therefore, the thesis explores the possibilities to expedite an electromagnetic simulation by scalable parallelization of near-linear complexity algorithms like Fast Multipole Method (FMM) on a multi-core platform. Secondly, with the best of parallelization strategies in place and near-linear complexity algorithms in use, the solution time of a complex EM problem can still be exceedingly large due to over-meshing of the geometry to achieve a desired level of accuracy. Hence, the thesis focuses on judicious placement of mesh elements on the geometry to capture the physics of the problem without compromising on accuracy- a technique called Adaptive Mesh Refinement. This facilitates a reduction in the number of solution variables or degrees of freedom in the system and hence the solution time. For multi-scale structures as encountered in chip-package-board systems, the MoM formulation breaks down for parts of the geometry having dimensions much smaller as compared to the operating wavelength. This phenomenon is popularly known as low-frequency breakdown or low-frequency instability. It results in an ill-conditioned MoM system matrix, and hence higher iteration count to converge when solved using an iterative solver framework. This consequently increases the solution time of simulation. The thesis thus proposes novel formulations to improve the spectral properties of the system matrix for real-world complex conductor and dielectric structures and hence form well-conditioned systems. This reduces the iteration count considerably for convergence and thus results in faster solution. Finally, minor changes in the geometrical design layouts can adversely affect the time-to-market of a commodity or a product. This is because the intermediate design variants, in spite of having similarities between them are treated as separate entities and therefore have to follow the conventional model-mesh-solve workflow for their analysis. This is a missed opportunity especially for design variant problems involving near-identical characteristics when the information from the previous design variant could have been used to expedite the simulation of the present design iteration. A similar problem occurs in the broadband simulation of an electromagnetic structure. The solution at a particular frequency can be expedited manifold if the matrix information from a frequency in its neighbourhood is used, provided the electrical characteristics remain nearly similar. The thesis introduces methods to re-use the subspace or Eigen-space information of a matrix from a previous design or frequency to solve the next incremental problem faster.