No description available.

Yes / Potential demodulating ability in biological tissue exposed to Radio Frequency (RF) signals intrinsically requires an unsymmetrical diode-like nonlinear response in tissue samples. This may be investigated by observing possible generation of the second harmonic in a cavity resonator designed to have fundamental and second harmonic resonant frequencies with collocated antinodes. Such a response would be of interest as being a mechanism that could enable demodulation of information-carrying waveforms having modulating frequencies in ranges that could interfere with cellular processes. Previous work has developed an experimental system to test for such responses: the present work reports an electric circuit model devised to facilitate calibration of any putative nonlinear RF energy conversion occurring within a nonlinear test-piece inside the cavity. The method is validated computationally and experimentally using a well-characterised nonlinear device. The variations of the reflection coefficients of the fundamental and second harmonic responses of the cavity due to adding nonlinear and lossy material are also discussed. The proposed model demonstrates that the sensitivity of the measurement equipment plays a vital role in deciding the required input power to detect any second harmonic signal, which is expected to be very weak. The model developed here enables the establishment of a lookup table giving the level of the second harmonic signal in the detector as a function of the specific input power applied in a measurement. Experimental results are in good agreement with the simulated results. / Engineering and Physical Science Research Council through Grant EP/E022936A

Higher-order basis functions are widely used to model currents and fields in numerical simulations of electromagnetics problems because of the greater accuracy and computational efficiency they can provide. Different problem formulations, such as method of moments (MoM) and the finite element method (FEM) require different constraints on basis functions for optimal performance, such as normal or tangential continuity between cells. In this thesis, a method of automatically generating bases that satisfy the desired basis constraints is applied to a MoM formulation for scattering problems using surface integral equations. Numerical results demonstrate the accuracy of this approach, and show good system matrix conditioning when compared to other higher-order bases.

[pt] O objetivo principal deste trabalho é o desenvolvimento de modelos computacionais para analisar a resposta eletromagnética de ferramentas de perfilagem LWD/MWD em formações geofísicas arbitrárias. Essa modelagem envolve a determinação precisa de campos eletromagnéticos em regiões tridimensionais (3D) complexas e, conseqüentemente, a solução de sistemas lineares não-hermitianos de larga escala. A modelagem numérica é realizada através da aplicação do método dos volumes finitos (FVM) no domínio da freqüência. Desenvolvem-se dois modelos computacionais, o primeiro válido em regiões isotrópicas e o segundo considerando a presença de anisotropias no meio. As equações de Maxwell são resolvidas através de duas formulações distintas: formulação por campos e formulação por potenciais vetor e escalar. A discretização por volumes finitos utiliza um esquema de grades entrelaçadas em coordenadas cilíndricas para evitar erros de aproximação de escada da geometria da ferramenta. Os modelos desenvolvidos incorporam quatro técnicas numéricas para aumentar a eficiência computacional e a precisão do método. As formulações por campos e por potenciais vetor e escalar são comparadas em termos da taxa de convergência e do tempo de processamento em cenários tridimensionais. Os modelos foram validados e testados em cenários tridimensionais complexos, tais como: (i) poços horizontais ou direcionais; (ii) formações não homogêneas com invasões de fluído de perfuração; (iii) formações anisotrópicas e (iv) poços excêntricos. Motivado pela flexibilidade dos modelos e pelos resultados numéricos obtidos em diferentes cenários tridimensionais, estende-se a metodologia para analisar a resposta de ferramentas LWD que empregam antenas inclinadas em relação ao eixo da ferramenta. Tais ferramentas podem prover dados com sensibilidade azimutal, assim como estimativas da anisotropia da formação, auxiliando o geodirecionamento de poços direcionais e horizontais. / [en] The main objective of this work is to develop computational models to analyze electromagnetic logging-while-drilling tool response in arbitrary geophysical formations. This modeling requires the determination of electromagnetic fields in three- dimensional (3-D) complex regions and consequently, the solution of large scale non-hermitian systems. The numerical modeling is done by using Finite Volume Methods (FVM) in the frequency domain. Both isotropic and anisotropic models are developed. Maxwell's equations are solved by using both the field formulation and the coupled vector-scalar potentials formulation. The proposed FVM technique utilizes an edge-based staggered-grid scheme in cylindrical coordinates to avoid staircasing errors on the tool geometry. Four numerical techniques are incorporated in the models in order to increase the computational efficiency and the accuracy of the method. The field formulation and the coupled vector-scalar potentials formulation are compared in terms of their accuracy, convergence rate, and CPU time for three-dimensional environments. The models were validated and tested in 3-D complex environments, such as:(i) horizontal and directional boreholes; (ii) multilayered geophysical formations including mud-filtrate invasions; (iii) anisotropic formations and (iv)eccentric boreholes. The methodology is extended to analyze LWD tools that are constructed with the transmitters and/or receivers tilted with respect to the axis of the drill collar. Such tools can provide improved anisotropy measurements and azimuthal sensitivity to benefit geosteering.

The phenomenon of polarization interaction in certain nonlinear materials is presented, and the design of an all-optical logic device based on this concept is described. An efficient two-dimensional finite-difference time-domain code for studying third-order nonlinear optical phenomena is discussed, in which both the slowly varying and the rapidly varying components of the electromagnetic fields are considered. The algorithm solves the vector form Maxwell s equations for all field components and uses the nonlinear constitutive relation in matrix form as the equations required to describe the nonlinear system. The stability of the code is discussed and its accuracy demonstrated through the simulation of the self-phase modulation effect observed in Kerr media. Finally, the code is used to simulate polarization mixing in photonic crystal-based line defect and coupled resonator optical waveguides.

In the thesis, we propose a newly-developed method called the method of Connected Local Fields (CLF) to analyze opto-electromagnetic passive devices. The method of CLF somewhat resembles a hybrid between the finite difference and pseudo-spectral methods. For opto-electromagnetic passive devices, our primary concern is their steady state behavior, or narrow-band characteristics, so we use a frequency-domain method, in which the system is governed by the Helmholtz equation. The essence of CLF is to use the intrinsic general solution of the Helmholtz equation to expand the local fields on the compact stencil. The original equation can then be transformed into the discretized form called LFE-9 (in 2-D case), and the intrinsic reconstruction formulae describing each overlapping local region can be obtained. Further, we present rigorous analysis of the numerical dispersion equation of LFE-9, by means of first-order approximation, and acquire the closed-form formula of the relative numerical dispersion error. We are thereby able to grasp the tangible influences brought both by the sampling density as well as the propagation direction of plane wave on dispersion error. In our dispersion analysis, we find that the LFE-9 formulation achieves the sixth-order accuracy: the theoretical highest order for discretizing elliptic partial differential equations on a compact nine-point stencil. Additionally, the relative dispersion error of LFE-9 is less than 1%, given that sampling density greater than 2.1 points per wavelength. At this point, the sampling density is nearing that of the Nyquist-Shannon sampling limit, and therefore computational efforts can be significantly reduced.

Optical whispering-gallery mode (WGM) cavities, which exhibit extraordinary spatial and temporal confinement of light, are one of the leading transducers for examining molecular recognition at low particle counts. With the advent of hybrid photonic-plasmonic and increasingly sophisticated forms of these resonators, the importance of supporting numerical methods has correspondingly become evident. In response, we adopt a full-vector finite difference approximation in order to solve for WGM's in terms of their field distributions, resonant wavelengths, and quality factors in the context of naturally discontinuous permittivity structure. A segmented Taylor series and alignment/rotation operator are utilized at such singularities in conjunction with arbitrarily spaced grid points. Simulations for microtoroids, with and without dielectric nanobeads, and plasmonic microdisks are demonstrated for short computation times and shown to be in agreement with data in the literature. Constricted surface plasmon polariton (SPP) WGM's are also featured within this document. The module of this thesis is devised as a keystone for composite WGM models that may guide experiments in the field. / Graduate

The Method of Moments (MoM) has been widely used for the full-wave electromagnetic analysis of planar multilayered media. However, the MoM is a computationally intensive process and requires considerable computer resources to perform the analysis. Thus, several efficient numerical techniques both in the spectral domain and spatial domain are investigated and further developed in this research. Two fitting procedures, i.e., the Rational Function Fitting Method (RFFM) and the Discrete Complex Image Method (DCIM), are investigated and developed in order to obtain closed-form spatial-domain Green's functions (GFs). Because the subtraction of the surface-wave pole contribution plays an important role for the accurate estimation of the spatial-domain GFs via DCIM, an efficient and accurate surface-wave pole location method is developed to find all the surface-wave poles for general multilayered media. The RFFM can be realized through either the Total Least Square Algorithm (TLSA) or the Vector Fitting (VECTFIT) method. Both the RFFM using VECTFIT and DCIM are detailed in step by step procedures. An efficient and low cost algorithm combining the advantages of DCIM and TLSA is also developed to evaluate the closed-form Green's functions for general multilayered media. A prototype version of the Full-Wave Layered-Interconnect Simulator (UA-FWLIS) was developed by analytically calculating the MoM reaction elements via Cauchy's residue theorem and the Complementary Incomplete Lipschitz-Hankel Integrals in stripline structures. After applying RFFM via VECTFIT to the G-functions, which are directly related to the spectral-domain GFs for microstrip structures, a procedure that is similar to the one used in the previously developed UA-FWLIS for stripline structures can be applied to calculate the MoM reaction elements analytically when the two reaction cells are close (< 0.1λ₀) in the spectral domain via the Electrical Field Integral Equation (EFIE). When the two reaction cells are far enough away (> 0.10.1λ₀), a simple expression for the reaction element can be obtained in the spatial domain via the Mixed Potential Integral Equation (MPIE) by a summation of a few complex image terms for the evaluation of the vector and scalar GFs. An efficient hybrid spectral-spatial method is thus developed to extend UA-FWLIS to microstrip structures. The method is validated by several numerical examples.

Maxwell's equations are a system of partial differential equations of vector fields. Imposing the constitutive relations for material properties yields equations for the curl and divergence of the electric and magnetic fields. The curl and divergence equations must be solved simultaneously, which is not the same as solving three separate scalar problems in each component of the vector field. This thesis describes a new method for solving partial differential equations of vector fields using the finite element method. New basis functions are used to solve the curl equation while allowing the divergence to be set as a constraint. The basis functions are defined on a mesh of bricks and the method is applicable for geometries that conform to a Cartesian coordinate system. The basis functions are a combination of cubic Hermite splines and second order Lagrange interpolation polynomials. The method yields a linearly independent set of constraints for the divergence, which is modelled to second order accuracy within each brick. Mesh refinement is accomplished by dividing selected bricks into $2\times 2\times 2$ smaller bricks of equal size. The change in the node pattern at an interface where mesh refinement occurs necessitates a modified implementation of the divergence constraints as well as additional constraints for hanging nodes. The mesh can be refined to an arbitrary number of levels. The basis functions can exactly model the discontinuity in the normal component of the field at a planar interface. The method is modified to solve problems with singularities at material boundaries that form $90^{\circ}$ edges and corners. The primary test problem of the new basis functions is to obtain the resonant frequencies and fields of three-dimensional cavities. The new basis functions can resolve physical solutions and non-physical, spurious modes. The eigenvalues obtained with the new method are in good agreement with exact solutions and experimental values in cases where they exist. There is also good agreement with results from second-order edge elements that are obtained with the software package HFSS. Finally, the method is modified to solve problems in cylindrical coordinates provided the domain does not contain the coordinate axis.

Maxwell's equations are a system of partial differential equations of vector fields. Imposing the constitutive relations for material properties yields equations for the curl and divergence of the electric and magnetic fields. The curl and divergence equations must be solved simultaneously, which is not the same as solving three separate scalar problems in each component of the vector field. This thesis describes a new method for solving partial differential equations of vector fields using the finite element method. New basis functions are used to solve the curl equation while allowing the divergence to be set as a constraint. The basis functions are defined on a mesh of bricks and the method is applicable for geometries that conform to a Cartesian coordinate system. The basis functions are a combination of cubic Hermite splines and second order Lagrange interpolation polynomials. The method yields a linearly independent set of constraints for the divergence, which is modelled to second order accuracy within each brick. Mesh refinement is accomplished by dividing selected bricks into $2\times 2\times 2$ smaller bricks of equal size. The change in the node pattern at an interface where mesh refinement occurs necessitates a modified implementation of the divergence constraints as well as additional constraints for hanging nodes. The mesh can be refined to an arbitrary number of levels. The basis functions can exactly model the discontinuity in the normal component of the field at a planar interface. The method is modified to solve problems with singularities at material boundaries that form $90^{\circ}$ edges and corners. The primary test problem of the new basis functions is to obtain the resonant frequencies and fields of three-dimensional cavities. The new basis functions can resolve physical solutions and non-physical, spurious modes. The eigenvalues obtained with the new method are in good agreement with exact solutions and experimental values in cases where they exist. There is also good agreement with results from second-order edge elements that are obtained with the software package HFSS. Finally, the method is modified to solve problems in cylindrical coordinates provided the domain does not contain the coordinate axis.