Because the metal enclosure with apertures is the structure of the cavity form, FDTD method belonging to full-wave analytical algorithm isn¡¦t efficient for analyzing such case. Thereby the ETL method belonging to analytical formulation provides another way to analyze it. Although the structure which it can apply to isn¡¦t as complex as that in full-wave analytical algorithm, the computing time is shorter than FDTD. In contrast, the method can¡¦t get good agreement. Thus two improved ways are introduced in this paper. One is to make the suitable structure wider and the other is to consider the higher mode original ETL method isn¡¦t included to improve the accuracy. Further these two ways make the suitable range of the method wider and the method has good practicability.

Σκοπός της παρούσας εργασίας είναι να παράσχει χρήσιμες πληροφορίες για την πρακτική κατασκευή συσκευών αορατότητας. Αρχικά, παρουσιάζονται οι διαφορετικοί τύποι των δισδιάστατων μέσων μετασχηματισμού (συσκευών αορατότητας) και διαπιστώνεται η υπεροχή του δισδιάστατου μανδύα αορατότητας επιπέδου από διηλεκτρικό. Προκειμένου να διερευνήσουμε αριθμητικά την οπτική συμπεριφορά αυτού του τύπου μανδύα, χρησιμοποιούμε τις εξισώσεις που εξάγονται από τη θεωρία της οπτικής μετασχηματισμών, και μέσω της μεθόδου F.D.T.D υπολογίζουμε το ηλεκτρομαγνητικό πεδίο σε όλα τα σημεία του υπολογιστικού πλέγματος. Στα όρια του πλέγματος αυτού τοποθετείται ένα στρώμα τέλειας προσαρμογής , έτσι ώστε να ελαχιστοποιούνται οι ανακλάσεις στα όρια αυτά, σε αναλογία με το πραγματικό πείραμα. Υλοποιώντας έναν κώδικα σε Μatlab, υπολογίζουμε την ένταση της σκεδαζόμενης ακτινοβολίας, για εγκάρσια μαγνητικά επίπεδα κύματα που προσπίπτουν υπό γωνία σε ένα αντικείμενο που περιβάλλεται από έναν μανδύα από διηλεκτρικό και είναι τοποθετημένο σε επίπεδο. Συγκεκριμένα, μελετάμε: την επίδοση της συσκευής συναρτήσει του βαθμού διακριτοποίησης του αριθμητικού πλέγματος, την επίδραση της απορρόφησης και της τυχαίας διαταραχής της επιτρεπτότητας των δομικών υλικών του μανδύα, καθώς και την επίδοση για πρόσπτωση υπό διαφορετικές γωνίες. Τέλος, εκτιμάμε τις ιδιότητες των στρώσεων από διηλεκτρικά που θα μπορούσαν να περιβάλουν τον μανδύα και να βελτιώσουν την αποδοτικότητά του στο ορατό φάσμα. / This work aims to provide useful information for the fabrication of practical cloaking devices. Initially, different types of two dimensional transformation media (invisibility devices) are presented and it is concluded that the all dielectric, two dimensional electromagnetic ground plane cloak dominates. In order to perform a numerical study of the sensitivity of this type of cloak, we use the equations deduced from the transformation optics via finite difference time domain method, for the electromagnetic wave, at every single point of the numerical domain. In the boundaries of this domain, we place a perfectly matched layer in order to minimize the reflections in those boundaries. That would replicate the real experiment with our simulations. By the application of a code written in Matlab, we calculate the scattering signature of an object surrounded by an all dielectric cloak and placed on a ground plane. Transverse magnetic plane waves are launched at an angle towards the object. In particular, we study the performance of the system as a function of the number of distinct components the cloak is divided into, the effects of lossy elements, the angle of incidence, as well as typical random variations of the permittivity of the building materials. Finally, we evaluate impedance matching layers that can surround the cloak and improve its effectiveness in the visible spectrum.

This thesis has proposed a simplified model of corona discharge from an overhead wire struck by lightning for surge computations using the FDTD method. In the corona model, the progression of corona streamers from the wire is represented as the radial expansion of cylindrical conducting region around the wire. The validity of this corona model has been tested against experimental data. Then, its applications to lightning electromagnetic pulse computations have been reviewed. / 博士(工学) / Doctor of Philosophy in Engineering / 同志社大学 / Doshisha University

This thesis covers recent advances in the adjoint variable method for the sensitivity estimations through time-domain electromagnetic simulations. It considers both frequency-independent and frequency-dependent response functions, and at the same time, provides a novel adjoint treatment for addressing dispersive sensitivity parameters in the material constitutive relation. With this proposed adjoint technique, response sensitivities with respect to all N sensitivity parameters can be computed through at most one extra simulations regardless of the value of N. This thesis also extends the existing adjoint technique to estimate all N^2 second-order sensitivity entries in the response Hessian matrix through N additional simulations. All adjoint sensitivity techniques presented in this thesis are numerically validated through various practical examples. Comparison shows that our produced adjoint results agree with those produced through central finite-difference approximations or through exact analytical approaches. / Dissertation / Doctor of Engineering (DEng)

風力発電システムは、現在世界中で普及しているが、その地上高と立地条件からしばしば落雷の被害を受け、稼働率の低下が問題視されている。本論文は、垂直導体である風力発電タワーにおける雷サージ解析の高精度化、高速化を目的とした垂直導体回路解析モデルの提案、およびその有用性のについて述べるとともに、風力発電所や洋上風車への拡張性についてまとめたものである。 / 博士(工学) / Doctor of Philosophy in Engineering / 同志社大学 / Doshisha University

Ground-Penetrating Radar (GPR) is a non-destructive electromagnetic investigative tool used in many applications across the fields of engineering and geophysics. The propagation of electromagnetic waves in lossy materials is complex and over the past 20 years, the computational modelling of GPR has developed to improve our understanding of this phenomenon. This research focuses on the development of accurate numerical models of widely-used, high-frequency commercial GPR antennas. High-frequency, highresolution GPR antennas are mainly used in civil engineering for the evaluation of structural features in concrete i. e., the location of rebars, conduits, voids and cracking. These types of target are typically located close to the surface and their responses can be coupled with the direct wave of the antenna. Most numerical simulations of GPR only include a simple excitation model, such as an infinitesimal dipole, which does not represent the actual antenna. By omitting the real antenna from the model, simulations cannot accurately replicate the amplitudes and waveshapes of real GPR responses. Numerical models of a 1.5 GHz Geophysical Survey Systems, Inc. (GSSI) antenna and a 1.2 GHz MALÅ GeoScience (MALÅ) antenna have been developed. The geometry of antennas is often complex with many fine features that must be captured in the numerical models. To visualise this level of detail in 3d, software was developed to link Paraview—an open source visualisation application which uses the Visualisation Toolkit (VTK)—with GprMax3D—electromagnetic simulation software based on the Finite-Difference Time-Domain (FDTD) method. Certain component values from the real antennas that were required for the models could not be readily determined due to commercial sensitivity. Values for these unknown parameters were found by implementing an optimisation technique known as Taguchi’s method. The metric used to initially assess the accuracy of the antenna models was a cross-corellation of the crosstalk responses from the models with the crosstalk responses measured from the real antennas. A 98 % match between modelled and real crosstalk responses was achieved. Further validation of the antenna models was undertaken using a series of laboratory experiments where oil-in-water (O/W) emulsions were created to simulate the electrical properties of real materials. The emulsions provided homogeneous liquids with controllable permittivity and conductivity and enabled different types of targets, typically encountered with GPR, to be tested. The laboratory setup was replicated in simulations which included the antenna models and an excellent agreement was shown between the measured and modelled data. The models reproduced both the amplitude and waveshape of the real responses whilst B-scans showed that the models were also accurately capturing effects, such as masking, present in the real data. It was shown that to achieve this accuracy, the real permittivity and conductivity profiles of materials must be correctly modelled. The validated antenna models were then used to investigate the radiation dynamics of GPR antennas. It was found that the shape and directivity of theoretically predicted far-field radiation patterns differ significantly from real antenna patterns. Being able to understand and visualise in 3d the antenna patterns of real GPR antennas, over realistic materials containing typical targets, is extremely important for antenna design and also from a practical user perspective.

First, one of the purposes of this thesis is to estimate the complete crosstalk effects including the package and the pads on the surface acoustic wave (SAW) substrate. A new approach based on finite-difference time-domain (FDTD) with equivalent current source method is applied. Two kinds of patterns of one-port SAW resonators with the same package structure and inter-digital transducer (IDT) design are studied. Verification with the measurement results shows that our method is able to obtain good agreement and be used to observe the influence from the SAW pattern. Second, the equivalent current source method is extended to model the excitation of human-body¡¦s electrostatic discharge (ESD) situations. The efficiencies of sacrificial electrodes are also discussed. Finally, a novel sacrificial electrode with fractal to protect SAW devices from ESD break is proposed. Comparing with traditional electrode, the simulation results show that fractal can improve the protective efficiency greatly. Finally, a novel analysis model that can be used to analyze and optimize the impedance of an RFID transponder integrated circuit (IC) which uses backscatter encoding based on simultaneously maintaining the BER of the reader and maximizing the received power of the transponder IC is proposed. The analysis method utilizes mapping from signal constellation of the backscattered signal to the Smith chart to relate the two parameters. Given the system specification and characteristics of the reader and transponder antennas, the optimum impedances of transponder IC for binary communication system can be easily designed by using this model.

Yes / This work presents the development of a new approach of modelling the source excitation and the penetration of structures by continuous propagating electromagnetic (EM) plane waves. The technique incorporates the solution of time-dependent Maxwell's equations and the initial value problem as the structures are illuminated by the plane waves. The propagation of waves from source excitation is simulated by solving a finite-difference Maxwell's equation in the time domain. Subgridding method is used to condense the lattice at the point of interest locally for observing field distribution in high resolution. The computational burden due to huge number of time steps has been eased by employing quasi-static approach. An example of induced EM fields near an underground pipeline runs parallel to a 132 kV overhead power transmission line (OHTL) has been presented which paves the way in the development of new approach of EM fields interaction modelling. / MSCRC

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.

There is an increasing need for accurate models describing the electrical behaviour of individual biological cells exposed to electromagnetic fields. In this area of solving linear problem, the most frequently used technique for computing the EM field is the Finite-Difference Time-Domain (FDTD) method. When modelling objects that are small compared with the wavelength, for example biological cells at radio frequencies, the standard Finite-Difference Time-Domain (FDTD) method requires extremely small time-step sizes, which may lead to excessive computation times. The problem can be overcome by implementing a quasi-static approximate version of FDTD, based on transferring the working frequency to a higher frequency and scaling back to the frequency of interest after the field has been computed. An approach to modeling and analysis of biological cells, incorporating the Hodgkin and Huxley membrane model, is presented here. Since the external medium of the biological cell is lossy material, a modified Berenger absorbing boundary condition is used to truncate the computation grid. Linear assemblages of cells are investigated and then Floquet periodic boundary conditions are imposed to imitate the effect of periodic replication of the assemblages. Thus, the analysis of a large structure of cells is made more computationally efficient than the modeling of the entire structure. The total fields of the simulated structures are shown to give reasonable and stable results at 900MHz, 1800MHz and 2450MHz. This method will facilitate deeper investigation of the phenomena in the interaction between EM fields and biological systems. Moreover, the nonlinear response of biological cell exposed to a 0.9GHz signal was discussed on observing the second harmonic at 1.8GHz. In this, an electrical circuit model has been proposed to calibrate the performance of nonlinear RF energy conversion inside a high quality factor resonant cavity with known nonlinear device. Meanwhile, the first and second harmonic responses of the cavity due to the loading of the cavity with the lossy material will also be demonstrated. The results from proposed mathematical model, give good indication of the input power required to detect the weakly effects of the second harmonic signal prior to perform the measurement. Hence, this proposed mathematical model will assist to determine how sensitivity of the second harmonic signal can be detected by placing the required specific input power.