Spelling suggestions: "subject:"computational electromagnetics"" "subject:"computational electromagneticas""
71 |
Modélisation de la propagation atmosphérique d'ondes électromagnétiques en 2D et 3D à partir de transformées de Fourier et en ondelettes / Modeling the atmospheric propagation of electromagnetic waves in 2D and 3D using fourier and wavelet transformsZhou, Hang 06 April 2018 (has links)
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.
|
72 |
A Spatially-filtered Finite-difference Time-domain Method with Controllable Stability Beyond the Courant LimitChang, Chun 19 July 2012 (has links)
This thesis introduces spatial filtering, which is a technique to extend the time step size beyond the conventional stability limit for the Finite-Difference Time-Domain (FDTD) method, at the expense of transforming field nodes between the spatial domain and the discrete spatial-frequency domain and removing undesired spatial-frequency components at every FDTD update cycle. The spatially-filtered FDTD method is demonstrated to be almost as accurate as and more efficient than the conventional FDTD method via theories and numerical examples. Then, this thesis combines spatial filtering and an existing subgridding scheme to form the spatially-filtered subgridding scheme. The spatially-filtered subgridding scheme is more efficient than existing subgridding schemes because the former allows the time step size used in the dense mesh to be larger than the dense mesh CFL limit. However, trade-offs between accuracy and efficiency are required in complicated structures.
|
73 |
A Spatially-filtered Finite-difference Time-domain Method with Controllable Stability Beyond the Courant LimitChang, Chun 19 July 2012 (has links)
This thesis introduces spatial filtering, which is a technique to extend the time step size beyond the conventional stability limit for the Finite-Difference Time-Domain (FDTD) method, at the expense of transforming field nodes between the spatial domain and the discrete spatial-frequency domain and removing undesired spatial-frequency components at every FDTD update cycle. The spatially-filtered FDTD method is demonstrated to be almost as accurate as and more efficient than the conventional FDTD method via theories and numerical examples. Then, this thesis combines spatial filtering and an existing subgridding scheme to form the spatially-filtered subgridding scheme. The spatially-filtered subgridding scheme is more efficient than existing subgridding schemes because the former allows the time step size used in the dense mesh to be larger than the dense mesh CFL limit. However, trade-offs between accuracy and efficiency are required in complicated structures.
|
74 |
Fast Solvers for Integtral-Equation based Electromagnetic SimulationsDas, Arkaprovo January 2016 (has links) (PDF)
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.
|
75 |
Study of RCS from Aerodynamic Flow using Parallel Volume-Surface Integral EquationPadhy, Venkat Prasad January 2016 (has links) (PDF)
Estimation of the Radar Cross Section of large inhomogeneous scattering objects such as composite aircrafts, ships and biological bodies at high frequencies has posed large computational challenge. The detection of scattering from wake vortex leading to detection and possible identification of low observable aircrafts also demand the development of computationally efficient and rigorous numerical techniques. Amongst the various methods deployed in Computational Electromagnetics, the Method of Moments predicts the electromagnetic characteristics accurately. Method of Moments is a rigorous method, combined with an array of modeling techniques such as triangular patch, cubical cell and tetrahedral modeling. Method of Moments has become an accurate technique for solving electromagnetic problems from complex shaped homogeneous and inhomogeneous objects. One of the drawbacks of Method of Moments is the fact that it results into a dense matrix, the inversion of which is a computationally complex both in terms of physical memory and compute power. This has been the prime reason for the Method of Moments hitherto remaining as a low frequency method. With recent advances in supercomputing, it is possible to extend the range of Method of Moments for Radar Cross Section computation of aircraft like structures and radiation characteristic of antennas mounted on complex shaped bodies at realistic frequencies of practical interest. This thesis is a contribution in this direction.
The main focus of this thesis is development of parallel Method of Moments solvers, applied to solve real world electromagnetic wave scattering and radiation problems from inhomogeneous objects. While the methods developed in this thesis are applicable to a variety of problems in Computational Electromagnetics as shown by illustrative examples, in specific, it has been applied to compute the Radar Cross Section enhancement due to acoustic disturbances and flow inhomogeneities from the wake vortex of an aircraft, thus exploring the possibility of detecting stealth aircraft. Illustrative examples also include the analysis of antenna mounted on an aircraft.
In this thesis, first the RWG basis functions have been used in Method of Moments procedure, for solving scattering problems from complex conducting structures such as aircraft and antenna(s) mounted on airborne vehicles, of electrically large size of about 45 and 0.76 million unknowns.
Next, the solver using SWG basis functions with tetrahedral and pulse basis functions with cubical modeling have been developed to solve scattering from 3D inhomogeneous bodies. The developed codes are validated by computing the Radar Cross Section of spherical homogeneous and inhomogeneous layered scatterers, lossy dielectric cylinder with region wise inhomogeneity and high contrast dielectric objects.
Aerodynamic flow solver ANSYS FLUENT, based on Finite Volume Method is used to solve inviscid compressible flow problem around the aircraft. The gradients of pressure/density are converted to dielectric constant variation in the wake region by using empirical relation and interpolation techniques. Then the Radar Cross Section is computed from the flow inhomogeneities in the vicinity of a model aircraft and beyond (wake zone) using the developed parallel Volume Surface Integral Equation using Method of Moments and investigated more rigorously. Radar Cross Section enhancement is demonstrated in the presence of the flow inhomogeneities and detectability is discussed. The Bragg scattering that occurs when electromagnetic and acoustic waves interact is also discussed and the results are interpreted in this light. The possibility of using the scattering from wake vortex to detect low visible aircraft is discussed.
This thesis also explores the possibility of observing the Bragg scattering phenomenon from the acoustic disturbances, caused by the wake vortex. The latter sets the direction for use of radars for target identification and beyond target detection.
The codes are parallelized using the ScaLAPACK and BiCG iterative method on shared and distributed memory machines, and tested on variety of High Performance Computing platforms such as Blue Gene/L (22.4TF), Tyrone cluster, CSIR-4PI HP Proliant 3000 BL460c (360TF) and CRAY XC40 machines. The parallelization speedup and efficiency of all the codes has also been shown.
|
76 |
Modélisation Multi-échelles : de l'Electromagnétisme à la Grille / Multi-scale Modeling : from Electromagnetism to GridKhalil, Fadi 14 December 2009 (has links)
Les performances des outils numériques de simulation électromagnétique de structures complexes, i.e., échelles multiples, sont souvent limitées par les ressources informatiques disponibles. De nombreux méso-centres, fermes et grilles de calcul, se créent actuellement sur les campus universitaires. Utilisant ces ressources informatiques mutualisées, ce travail de thèse s'attache à évaluer les potentialités du concept de grille de calcul (Grid Computing) pour la simulation électromagnétique de structures multi-échelles. Les outils numériques de simulation électromagnétique n'étant pas conçus pour être utilisés dans un environnement distribué, la première étape consistait donc à les modifier afin de les déployer sur une grille de calcul. Une analyse approfondie a ensuite été menée pour évaluer les performances des outils de simulation ainsi déployés sur l'infrastructure informatique. Des nouvelles approches pour le calcul électromagnétique distribué avec ces outils sont présentées et validées. En particulier, ces approches permettent la réalisation de simulation électromagnétique de structures à échelles multiples en un temps record et avec une souplesse d'utilisation. / The numerical electromagnetic tools for complex structures simulation, i.e. multi-scale, are often limited by available computation resources. Nowadays, Grid computing has emerged as an important new field, based on shared distributed computing resources of Universities and laboratories. Using these shared resources, this study is focusing on grid computing potential for electromagnetic simulation of multi-scale structure. Since the numerical simulations tools codes are not initially written for distributed environment, the first step consists to adapt and deploy them in Grid computing environment. A performance study is then realized in order to evaluate the efficiency of execution on the test-bed infrastructure. New approaches for distributing the electromagnetic computations on the grid are presented and validated. These approaches allow a very remarkable simulation time reduction for multi-scale structures and friendly-user interfaces
|
77 |
[en] INVESTIGATION OF ELECTROMAGNETIC PROPAGATION IN PLASMA STRUCTURES THROUGH EIGENFUNCTION EXPANSIONS AND FDTD TECHNIQUES / [pt] INVESTIGAÇÃO DE PROPAGAÇÃO ELETROMAGNÉTICA EM ESTRUTURAS DE PLASMA ATRAVÉS DE EXPANSÕES EM AUTOFUNÇÕES E TÉCNICAS FDTDJULIO DE LIMA NICOLINI 18 July 2017 (has links)
[pt] Plasma é um dos quatro estados fundamentais da matéria, presente em forma natural na Terra na ionosfera, em relâmpagos e nas chamas resultantes de combustão, assim como em forma artificial em lâmpadas de neônio, lâmpadas fluorescentes e processos industriais. O comportamento de plasmas é extraordinariamente complexo e variado, como por exemplo a formação espontânea de características espaciais interessantes em variadas escalas diferentes de comprimento. Uma antena de plasma, por sua vez, é uma estrutura radiante baseada em um elemento de plasma em vez de um condutor metálico, o que gera diversas vantagens e características úteis de um ponto de vista tecnológico. Nesse presente trabalho, uma investigação da propagação eletromagnética dentro de estruturas de plasma é realizada através de métodos teóricos e numéricos como um primeiro passo em direção ao desenvolvimento de modelos apropriados para o estudo de antenas de plasma. / [en] Plasma is one of the four fundamental states of matter, present on Earth in natural form at the ionosphere, in lightning strikes and in the flames resulting from combustion, as well as in artificial form in neon signs, fluorescent light bulbs and industrial processes. Plasma behaviour is extraordinarily complex and varied, e.g. the spontaneous formation of interesting spatial features over a wide range of length scales. A plasma antenna, on the other hand, is a radiating structure based in a plasma element instead of a metallic conductor, which creates several technological advantages and useful characteristics. In this present work, an investigation of electromagnetic propagation inside of plasma structures is performed through both theoretical and numerical means as a first step towards constructing appropriate models for the study of plasma antennas.
|
78 |
Investigation of Near-Field Contribution in SBR for Installed Antenna PerformanceHultin, Harald January 2019 (has links)
To investigate near-eld contributions for installed antennas, an in-house code iswritten to incorporate near-eld terms in Shooting and Bouncing Rays (SBR). SBRis a method where rays are launched toward an object and scatter using GeometricalOptics (GO). These rays induce currents on the object, from which the totalscattered eld can be found.To gauge the eect of near-eld terms, the in-house code can be set to excludenear-eld terms. Due to this characteristic, the method is named SBR Includingor Excluding Near-eld Terms (SIENT). The SIENT implementation is thoroughlydescribed. To make SIENT more exible, the code works with triangulated meshesof objects. Antennas are represented as near-eld sources, allowing complex antennasto be represented by simple surface currents. Further, some implementedoptimizations of SIENT are shown.To test the implemented method, SIENT is compared to a reference solution andcomparable commercial SBR solvers. It is shown that SIENT compares well to thecommercial options. Further, it is shown that the inclusion of near-eld terms actsas a small correction to the far-eld of the installed antenna. / För att undersöka närfältsbidrag för installerade antenner, har en kod skrivits föratt ta med närfältstermer i Shooting Bouncing Rays (SBR). SBR är en metod där strålar (”rays”) skjuts mot ett object och sprids via Geometrisk Optik (GO). Dessa strålar inducerar strömmar på objectet, från vilka det totala sprida fältet kan hittas. För att undersöka bidraget från närfältstermer, så kan koden exkludera dessa. På grund av denna karaktär, kallas koden SBR Including or Excluding Near-field Terms (SIENT). Implementationen av SIENT beskrivs utförligt. För att göra SIENT mer flexibel, arbetar SIENT med triangulerade nät av objekt. Antenner representeras av närfältskällor, vilket låter komplexa antenner representeras med enkla yt-strömmar.Implementerade optimeringar av SIENT visas också.För att testa den implementerade metoden, jämförs SIENT med en referenslösning och jämförbara kommerciella SBR-lösare. Det visas att SIENT överensstämmer bra med kommerciella alternativ. Det visas också att närfältstermer agerar som enmindre korrektion till fjärrfältet av den installerade antennen.
|
79 |
Enhancing numerical modelling efficiency for electromagnetic simulation of physical layer componentsSasse, Hugh Granville January 2010 (has links)
The purpose of this thesis is to present solutions to overcome several key difficulties that limit the application of numerical modelling in communication cable design and analysis. In particular, specific limiting factors are that simulations are time consuming, and the process of comparison requires skill and is poorly defined and understood. When much of the process of design consists of optimisation of performance within a well defined domain, the use of artificial intelligence techniques may reduce or remove the need for human interaction in the design process. The automation of human processes allows round-the-clock operation at a faster throughput. Achieving a speedup would permit greater exploration of the possible designs, improving understanding of the domain. This thesis presents work that relates to three facets of the efficiency of numerical modelling: minimizing simulation execution time, controlling optimization processes and quantifying comparisons of results. These topics are of interest because simulation times for most problems of interest run into tens of hours. The design process for most systems being modelled may be considered an optimisation process in so far as the design is improved based upon a comparison of the test results with a specification. Development of software to automate this process permits the improvements to continue outside working hours, and produces decisions unaffected by the psychological state of a human operator. Improved performance of simulation tools would facilitate exploration of more variations on a design, which would improve understanding of the problem domain, promoting a virtuous circle of design. The minimization of execution time was achieved through the development of a Parallel TLM Solver which did not use specialized hardware or a dedicated network. Its design was novel because it was intended to operate on a network of heterogeneous machines in a manner which was fault tolerant, and included a means to reduce vulnerability of simulated data without encryption. Optimisation processes were controlled by genetic algorithms and particle swarm optimisation which were novel applications in communication cable design. The work extended the range of cable parameters, reducing conductor diameters for twisted pair cables, and reducing optical coverage of screens for a given shielding effectiveness. Work on the comparison of results introduced ―Colour maps‖ as a way of displaying three scalar variables over a two-dimensional surface, and comparisons were quantified by extending 1D Feature Selective Validation (FSV) to two dimensions, using an ellipse shaped filter, in such a way that it could be extended to higher dimensions. In so doing, some problems with FSV were detected, and suggestions for overcoming these presented: such as the special case of zero valued DC signals. A re-description of Feature Selective Validation, using Jacobians and tensors is proposed, in order to facilitate its implementation in higher dimensional spaces.
|
80 |
Fast algorithms for frequency domain wave propagationTsuji, Paul Hikaru 22 February 2013 (has links)
High-frequency wave phenomena is observed in many physical settings, most notably in acoustics, electromagnetics, and elasticity. In all of these fields, numerical simulation and modeling of the forward propagation problem is important to the design and analysis of many systems; a few examples which rely on these computations are the development of
metamaterial technologies and geophysical prospecting for natural resources. There are two modes of modeling the forward problem: the frequency domain and the time domain. As the title states, this work is concerned with the former regime.
The difficulties of solving the high-frequency wave propagation problem accurately lies in the large number of degrees of freedom required. Conventional wisdom in the computational electromagnetics commmunity suggests that about 10 degrees of freedom per wavelength be used in each coordinate direction to resolve each oscillation. If K is the width of the domain in wavelengths, the number of unknowns N grows at least by O(K^2) for surface discretizations and O(K^3) for volume discretizations in 3D. The memory requirements and asymptotic complexity estimates of direct algorithms such as the multifrontal method are too costly for such problems. Thus, iterative solvers must be used. In this dissertation, I will present fast algorithms which, in conjunction with GMRES, allow the solution of the forward problem in O(N) or O(N log N) time. / text
|
Page generated in 0.4692 seconds