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

The main objective of this research is to model and analyse complex electromagnetic problems by means of a new hybridised computational technique combining the frequency domain Method of Moments (MoM), Finite-Difference Time-Domain (FDTD) method and a subgridded Finite-Difference Time-Domain (SGFDTD) method. This facilitates a significant advance in the ability to predict electromagnetic absorption in inhomogeneous, anisotropic and lossy dielectric materials irradiated by geometrically intricate sources. The Method of Moments modelling employed a two-dimensional electric surface patch integral formulation solved by independent linear basis function methods in the circumferential and axial directions of the antenna wires. A similar orthogonal basis function is used on the end surface and appropriate attachments with the wire surface are employed to satisfy the requirements of current continuity. The surface current distributions on structures which may include closely spaced parallel wires, such as dipoles, loops and helical antennas are computed. The results are found to be stable and showed good agreement with less comprehensive earlier work by others. The work also investigated the interaction between overhead high voltage transmission lines and underground utility pipelines using the FDTD technique for the whole structure, combined with a subgridding method at points of interest, particularly the pipeline. The induced fields above the pipeline are investigated and analysed. FDTD is based on the solution of Maxwell¿s equations in differential form. It is very useful for modelling complex, inhomogeneous structures. Problems arise when open-region geometries are modelled. However, the Perfectly Matched Layer (PML) concept has been employed to circumvent this difficulty. The establishment of edge elements has greatly improved the performance of this method and the computational burden due to huge numbers of time steps, in the order of tens of millions, has been eased to tens of thousands by employing quasi-static methods. This thesis also illustrates the principle of the equivalent surface boundary employed close to the antenna for MoM-FDTD-SGFDTD hybridisation. It depicts the advantage of using hybrid techniques due to their ability to analyse a system of multiple discrete regions by employing the principle of equivalent sources to excite the coupling surfaces. The method has been applied for modelling human body interaction with a short range RFID antenna to investigate and analyse the near field and far field radiation pattern for which the cumulative distribution function of antenna radiation efficiency is presented. The field distributions of the simulated structures show reasonable and stable results at 900 MHz. This method facilitates deeper investigation of the phenomena in the interaction between electromagnetic fields and human tissues. / Ministry of Higher Education Malaysia and Universiti Tun Hussein Onn Malaysia (UTHM)

No description available.

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.