DOCSLIB.ORG

Perfectly Matched Layers for T,Φ Formulation

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Perfectly Matched Layers for T,Φ Formulation Proceedings - 49 - 12th International IGTE Symposium Perfectly Matched Layers for T, Φ Formulation Gernot Matzenauer, Oszkar´ B´ıro,´ Karl Hollaus, Kurt Preis and Werner Renhart Institute for Fundamentals and Theory in Electrical Engineering, Kopernikusgasse 24/3, A-8010 Graz, Austria E-mail: [email protected] Abstract: Perfectly Matched Layers (PMLs) are used for II. PROOF OF THE PERFECT MATCHING PROPERTY reflectionless truncation of the problem boundaries in FEM applications. The basic concept behind PMLs is to create an Based on [2] it is shown below that the PML interface artificial material with a complex and diagonal anisotropic also works with any lossless material in front of the PML. permittivity and permeability. For the A,V formulation Similarly to (2), the permittivity and permeability tensors PMLs are well known. In this paper the method of perfectly of the PML are assumed to be: matched layers is extended to truncate any lossless medium and the method is implemented for the T, Φ formulation. [] = 0r[Λ], Keywords: Perfectly Matched Layer, Finite Element Method, T, Φ formulation, Perfectly Matched Anisotropic [µ] = µ0µr[Λ]. (4) Absorber, Absorbing Boundary Condition The question is how to choose the parameters a, b, c, µr and r, so that no reflection occurs on an interface between I. INTRODUCTION ∗ the PML and a medium with permeability µr and permit- ∗ In the present paper, an artificial anisotropic lossy tivity r. The wave equations for the electric and the mag- material has been applied to a 3D edge finite-element T, Φ netic field in the PML are: formulation [1] to form perfectly matched layers. This −1 −1 2 material fully absorbs the electromagnetic field impinging [Λ] ∇ × [Λ] ∇ × E − ω µ0 µr0 rE = 0, on it without reflections on the PML-air interface [2]. −1 −1 2 [Λ] ∇ × [Λ] ∇ × H − ω µ0 µr0 rH = 0. The permittivity [] and the permeability [µ] in PMLs (5) are considered as complex diagonal tensors. The following General solutions for the electric and the magnetic field condition is required to match the intrinsic impedances of can be constructed from plane waves of the form vacuum and the PML[2]: j(ωt−kr) E = E0 e µ0 −1 [1] = [µ][] (1) j(ωt−kr) 0 H = H0 e (6) where [1] is the unit tensor. Thus, the permittivity [] and where k = kx ex + ky ey + kz ez , r = x ex + y ey + the permeability [µ] can be written as z ez and E0 and H0 are constant vectors. The dispersion relation of the PML leads to an equation of the form: a 0 0 k2 k2 k2 [] = 0 b 0 = [Λ], x y z 2 0 0 + + = k0. (7) 0 0 c bc ac ab where a 0 0 2 2 k0 = ω µ0 µr 0 r (8) [µ] = µ0 0 b 0 = µ0[Λ]. (2) 0 0 c is the propagation constant. The solution of the dispersion relation is: Here a, b and c are complex numbers describing the mate- √ rial properties of the anisotropic absorber. Consequently, kx = k0 bc sin θ cos φ Maxwell’s equations in the PML are √ ky = k0 ac sin θ sin φ ∇ × H = jω0[Λ]E √ kz = k0 ab cos θ (9) ∇ × E = −jωµ0[Λ]H To find out which component of the electric field belongs ∇ · [Λ]E = 0 to which component of the magnetic field, the solutions 0 for the electric and the magnetic field (6) are inserted into Maxwell’s equations. Using Faraday’s law yields ∇ · µ0[Λ]H = 0. (3) r ! 1 r c b The Maxwell equations are not modified with this ap- Hx = Ez sin θ sin φ − Ey cos θ , (10) proach and the implementation is straightforward in exist- Z0 a a ing solvers with finite element methods (FEM). It is suffi- 1 ra rc cient for the code to permit complex anisotropic materials. Hy = Ex cos θ − Ez sin θ cos φ , In this paper the method of perfectly matched layers is ex- Z0 b b tended to truncate any lossless medium and the method is r r ! implemented for the T, Φ formulation. The results are il- 1 b a Hz = Ey sin θ cos φ − Ex sin θ sin φ . lustrated by two numerical examples. Z0 c c 12th International IGTE Symposium - 50 - Proceedings Using Ampere’s law leads to r ! b r c E = Z H cos θ − H sin θ sin φ , (11) x 0 y a z a rc ra E = Z H sin θ cos φ − H cos θ , y 0 z b x b r ! ra b E = Z H sin θ sin φ − H sin θ cos φ . z 0 x c y c where r µ0 µr Z0 = (12) 0 r Figure 1. Plane wave in TE mode incidence on a z-plane interface. is the wave impedance of the PML. Let us consider a plane wave incident on a PML as shown is the propagation vector and re and rr are the directions of in Fig. 1. The values of a, b, c, µr and r are to be deter- the exciting and the reflected wave. The transmitted wave mined so that independent of the incident angle αe no re- is given by: flection occurs on the interface. The reflection coefficient −j k·rt r and the transmission coefficient t are defined by: Ety = E0t e (16) r 1 b E −j k·rt r = 0r Htx = − E0t cos αt e Z0 a E0e r 1 b E0t −j k·rt t = (13) Htz = E0t sin αt e E0e Z0 c √ √ k · rt = k0(x bc sin αt + z ab cos αt) where E0e, E0r and E0t are the magnitudes of the excit- ing, the reflected and the transmitted electric field. The where r is the direction of the transmitted wave. Here Z∗ condition for no reflection is that the reflection coefficient t 0 is the wave impedance of the lossless linear space: vanishes r = 0 and the transmission coefficient is t = 1. The propagation vectors are restricted to the x − z-plane. s µ µ∗ For the TE mode, the excited and the reflected waves are Z∗ = 0 r . (17) 0 ∗ given by: 0 r Enforcing the tangential continuity of E on the interface at −j k·re Eey = E0e e (14) the position z = 0 yields the relation: √ 1 −jk∗ x sin α −jk∗ x sin α −jk x bc sin α −j k·re e 0 e + r e 0 r = t e 0 t . (18) Hex = − ∗ E0e cos αe e Z0 The phase matching condition of the tangential component 1 −j k·re of the electric field is given by: Hez = ∗ E0e sin αe e Z0 √ ∗ ∗ k0 sin αe = k0 sin αr = k0 bc sin αt. (19) k · re = k0(x sin αe + z cos αe) Provided −j k·rr Ery = E0r e √ bc = 1 (20) 1 −j k·rr Hrx = ∗ E0r cos αr e Z0 and 1 ∗ ∗ −j k·rr µr r = µr r, (21) Hrz = ∗ E0r sin αr e Z0 ∗ k · rr = k0 (x sin αr − z cos αr) αe = αr = αt (22) is obtained, i.e. no reflection occurs. Using (19), (18) be- where comes ∗2 2 ∗ ∗ k0 = ω µ0 µr 0 r (15) 1 + r = t. (23) Proceedings - 51 - 12th International IGTE Symposium Enforcing the tangential continuity of H and using (19) potential) realized by edge and nodal finite elements [1]. leads to Furthermore, the properties of the PML layer have been ∗ r extended to truncate any lossless linear material. Z0 b cos αe(1 − r) = t cos αt. (24) Z0 a The field quantities for the T, Φ formulation are derived The magnitude matching equations yield the reflection co- from the potentials as efficient H = T − gradΦ, ∗ q Z0 b cos αe − Z a cos αt J = curlT (32) r = 0 . (25) TE ∗ q Z0 b cos αe + cos αt Z0 a and satisfy the constitutive relations: B = [µ]H, D = []E and E = [ρ]J where [µ],[] and [ρ] are the tensors of For the TM mode the same procedure can be followed to permeability, permittivity and resistivity. The permeability determine the reflection coefficient: and conductivity may be written in complex form: ∗ q Z0 b Z a cos αt − cos αe [σ] = [σ ] + jω[], r = 0 (26) E TM ∗ q Z0 b cos αt + cos αe [σ ] Z0 a [µ] = M + [µ]. (33) jω Using the requirement (22) the reflection coefficient is in- dependent of the angle. By setting A magnetic conductivity tensor [σM ] is introduced to de- scribe the material properties of the PML. The complex a = b (27) material tensor for the electric resistivity can be obtained from and −1 ∗ [ρ] = ([σ]) . (34) Z0 = Z0 (28) the interface will be perfectly reflectionless for any fre- The tensors of permeability and resistivity have to be quency and angle of incidence. complex in the whole region to imply wave propagation with the T, Φ formulation. In summary, (20) and (27) lead to The governing differential equations for the T, Φ formu- 1 a = b = (29) lation are obtained from Faraday’s law and from the van- c ishing divergence of the magnetic induction: and the relations (21) and (28) yield ∇ × [ρ] ∇ × T + j ω [µ]T − j ω [µ]∇Φ = 0, ∗ r = r, ∇ · [µ]T − [µ]∇Φ = 0. (35) ∗ µr = µr. (30) The material properties of a perfectly matched layer can be Any arbitrarily polarized plane wave can be decomposed obtained from the complex PML tensor [Λ]: into a linear combination of TE and TM modes. So under the conditions (29) and (30) the PML will work with any e(1 − j) 0 0 polarization.
Load more

Recommended publications
  • Cvel-08-011.3
    TECHNICAL REPORT: CVEL-08-011.3 Survey of Current Computational Electromagnetics Techniques and Software T. Hubing, C. Su, H. Zeng and H. Ke Clemson University June 6, 2009 This report is an update of the 2008 report, CVEL-08-011.2. This work was supported by the ARMY RESEARCH OFFICE (W. Devereux Palmer) under the auspices of the U.S. Army Research Office Scientific Services Program administered by Battelle (Delivery Order 0271, Contract No. W911NF-07-D-0001). Table of Contents I. Introduction ............................................................................................................................. 3 A. Purpose of This Report ................................................................................................................. 3 B. Categorizing CEM Modeling Tools .............................................................................................. 4 1. Time vs. Frequency Domain ........................................................................................ 4 2. IE vs. PDE .................................................................................................................... 4 3. 2D vs. 3D ...................................................................................................................... 5 4. Static, Quasi-static, Full-wave or Asymptotic .............................................................. 5 5. Circuit-based Field Solvers .......................................................................................... 6 6. Linear vs. Higher-order Methods
    [Show full text]
  • RF Coils, Helical Resonators and Voltage Magnification by Coherent Spatial Modes
    TELSIKS 2001, University of Nis, Yugoslavia (September 19-21, 2001) and MICROWAVE REVIEW 1 RF Coils, Helical Resonators and Voltage Magnification by Coherent Spatial Modes K.L. Corum* and J.F. Corum** “Is there, I ask, can there be, a more interesting study than that of alternating currents.” Nikola Tesla, (Life Fellow, and 1892 Vice President of the AIEE)1 Abstract – By modeling a wire-wound coil as an anisotropically II. CYLINDRICAL HELICES conducting cylindrical boundary, one may start from Maxwell’s equations and deduce the structure’s resonant behavior. Not A. Problem Formulation only can the propagation factor and characteristic impedance be determined for such a helically disposed surface waveguide, but also its resonances, “self-capacitance” (so-called), and its voltage A uniform helix is described by its radius (r = a), its pitch magnification by standing waves. Further, the Tesla coil passes (or turn-to-turn wire spacing "s"), and its pitch angle ψ, to a conventional lumped element inductor as the helix is which is the angle that the tangent to the helix makes with a electrically shortened. plane perpendicular to the axis of the structure (z). Geometrically, ψ = cot-1(2πa/s). The wave equation is not Keywords – coil, helix, surface wave, resonator, inductor, Tesla. separable in helical coordinates and there exists no rigorous 3 solution of Maxwell's equations for the solenoidal helix. I. INTRODUCTION However, at radio frequencies a wire-wound helix with many turns per free-space wavelength (e.g., a Tesla coil) One of the more significant challenges in electrical may be modeled as an idealized anisotropically conducting cylindrical surface that conducts only in the helical science is that of constricting a great deal of insulated conductor into a compact spatial volume and determining direction.
    [Show full text]
  • The Self-Resonance and Self-Capacitance of Solenoid Coils: Applicable Theory, Models and Calculation Methods
    1 The self-resonance and self-capacitance of solenoid coils: applicable theory, models and calculation methods. By David W Knight1 Version2 1.00, 4th May 2016. DOI: 10.13140/RG.2.1.1472.0887 Abstract The data on which Medhurst's semi-empirical self-capacitance formula is based are re-analysed in a way that takes the permittivity of the coil-former into account. The updated formula is compared with theories attributing self-capacitance to the capacitance between adjacent turns, and also with transmission-line theories. The inter-turn capacitance approach is found to have no predictive power. Transmission-line behaviour is corroborated by measurements using an induction loop and a receiving antenna, and by visualising the electric field using a gas discharge tube. In-circuit solenoid self-capacitance determinations show long-coil asymptotic behaviour corresponding to a wave propagating along the helical conductor with a phase-velocity governed by the local refractive index (i.e., v = c if the medium is air). This is consistent with measurements of transformer phase error vs. frequency, which indicate a constant time delay. These observations are at odds with the fact that a long solenoid in free space will exhibit helical propagation with a frequency-dependent phase velocity > c. The implication is that unmodified helical-waveguide theories are not appropriate for the prediction of self-capacitance, but they remain applicable in principle to open- circuit systems, such as Tesla coils, helical resonators and loaded vertical antennas, despite poor agreement with actual measurements. A semi-empirical method is given for predicting the first self- resonance frequencies of free coils by treating the coil as a helical transmission-line terminated by its own axial-field and fringe-field capacitances.
    [Show full text]
  • Ec6503 - Transmission Lines and Waveguides Transmission Lines and Waveguides Unit I - Transmission Line Theory 1
    EC6503 - TRANSMISSION LINES AND WAVEGUIDES TRANSMISSION LINES AND WAVEGUIDES UNIT I - TRANSMISSION LINE THEORY 1. Define – Characteristic Impedance [M/J–2006, N/D–2006] Characteristic impedance is defined as the impedance of a transmission line measured at the sending end. It is given by 푍 푍0 = √ ⁄푌 where Z = R + jωL is the series impedance Y = G + jωC is the shunt admittance 2. State the line parameters of a transmission line. The line parameters of a transmission line are resistance, inductance, capacitance and conductance. Resistance (R) is defined as the loop resistance per unit length of the transmission line. Its unit is ohms/km. Inductance (L) is defined as the loop inductance per unit length of the transmission line. Its unit is Henries/km. Capacitance (C) is defined as the shunt capacitance per unit length between the two transmission lines. Its unit is Farad/km. Conductance (G) is defined as the shunt conductance per unit length between the two transmission lines. Its unit is mhos/km. 3. What are the secondary constants of a line? The secondary constants of a line are 푍 i. Characteristic impedance, 푍0 = √ ⁄푌 ii. Propagation constant, γ = α + jβ 4. Why the line parameters are called distributed elements? The line parameters R, L, C and G are distributed over the entire length of the transmission line. Hence they are called distributed parameters. They are also called primary constants. The infinite line, wavelength, velocity, propagation & Distortion line, the telephone cable 5. What is an infinite line? [M/J–2012, A/M–2004] An infinite line is a line where length is infinite.
    [Show full text]
  • Modal Propagation and Excitation on a Wire-Medium Slab
    1112 IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 56, NO. 5, MAY 2008 Modal Propagation and Excitation on a Wire-Medium Slab Paolo Burghignoli, Senior Member, IEEE, Giampiero Lovat, Member, IEEE, Filippo Capolino, Senior Member, IEEE, David R. Jackson, Fellow, IEEE, and Donald R. Wilton, Life Fellow, IEEE Abstract—A grounded wire-medium slab has recently been shown to support leaky modes with azimuthally independent propagation wavenumbers capable of radiating directive omni- directional beams. In this paper, the analysis is generalized to wire-medium slabs in air, extending the omnidirectionality prop- erties to even modes, performing a parametric analysis of leaky modes by varying the geometrical parameters of the wire-medium lattice, and showing that, in the long-wavelength regime, surface modes cannot be excited at the interface between the air and wire medium. The electric field excited at the air/wire-medium interface by a horizontal electric dipole parallel to the wires is also studied by deriving the relevant Green’s function for the homogenized slab model. When the near field is dominated by a leaky mode, it is found to be azimuthally independent and almost perfectly linearly polarized. This result, which has not been previ- ously observed in any other leaky-wave structure for a single leaky mode, is validated through full-wave moment-method simulations of an actual wire-medium slab with a finite size. Index Terms—Leaky modes, metamaterials, near field, planar Fig. 1. (a) Wire-medium slab in air with the relevant coordinate axes. (b) Transverse view of the wire-medium slab with the relevant geometrical slab, wire medium.
    [Show full text]
  • Waveguide Propagation
    NTNU Institutt for elektronikk og telekommunikasjon Januar 2006 Waveguide propagation Helge Engan Contents 1 Introduction ........................................................................................................................ 2 2 Propagation in waveguides, general relations .................................................................... 2 2.1 TEM waves ................................................................................................................ 7 2.2 TE waves .................................................................................................................... 9 2.3 TM waves ................................................................................................................. 14 3 TE modes in metallic waveguides ................................................................................... 14 3.1 TE modes in a parallel-plate waveguide .................................................................. 14 3.1.1 Mathematical analysis ...................................................................................... 15 3.1.2 Physical interpretation ..................................................................................... 17 3.1.3 Velocities ......................................................................................................... 19 3.1.4 Fields ................................................................................................................ 21 3.2 TE modes in rectangular waveguides .....................................................................
    [Show full text]
  • Introduction to RF Measurements and Instrumentation
    Introduction to RF measurements and instrumentation Daniel Valuch, CERN BE/RF, [email protected] Purpose of the course • Introduce the most common RF devices • Introduce the most commonly used RF measurement instruments • Explain typical RF measurement problems • Learn the essential RF work practices • Teach you to measure RF structures and devices properly, accurately and safely to you and to the instruments Introduction to RF measurements and instrumentation 2 Daniel Valuch CERN BE/RF ([email protected]) Purpose of the course • What are we NOT going to do… But we still need a little bit of math… Introduction to RF measurements and instrumentation 3 Daniel Valuch CERN BE/RF ([email protected]) Purpose of the course • We will rather focus on: Instruments: …and practices: Methods: Introduction to RF measurements and instrumentation 4 Daniel Valuch CERN BE/RF ([email protected]) Transmission line theory 101 • Transmission lines are defined as waveguiding structures that can support transverse electromagnetic (TEM) waves or quasi-TEM waves. • For purpose of this course: The device which transports RF power from the source to the load (and back) Introduction to RF measurements and instrumentation 5 Daniel Valuch CERN BE/RF ([email protected]) Transmission line theory 101 Transmission line Source Load Introduction to RF measurements and instrumentation 6 Daniel Valuch CERN BE/RF ([email protected]) Transmission line theory 101 • The telegrapher's equations are a pair of linear differential equations which describe the voltage (V) and current (I) on an electrical transmission line with distance and time. • The transmission line model represents the transmission line as an infinite series of two-port elementary components, each representing an infinitesimally short segment of the transmission line: Distributed resistance R of the conductors (Ohms per unit length) Distributed inductance L (Henries per unit length).
    [Show full text]
  • Waveguides Waveguides, Like Transmission Lines, Are Structures Used to Guide Electromagnetic Waves from Point to Point. However
    Waveguides Waveguides, like transmission lines, are structures used to guide electromagnetic waves from point to point. However, the fundamental characteristics of waveguide and transmission line waves (modes) are quite different. The differences in these modes result from the basic differences in geometry for a transmission line and a waveguide. Waveguides can be generally classified as either metal waveguides or dielectric waveguides. Metal waveguides normally take the form of an enclosed conducting metal pipe. The waves propagating inside the metal waveguide may be characterized by reflections from the conducting walls. The dielectric waveguide consists of dielectrics only and employs reflections from dielectric interfaces to propagate the electromagnetic wave along the waveguide. Metal Waveguides Dielectric Waveguides Comparison of Waveguide and Transmission Line Characteristics Transmission line Waveguide • Two or more conductors CMetal waveguides are typically separated by some insulating one enclosed conductor filled medium (two-wire, coaxial, with an insulating medium microstrip, etc.). (rectangular, circular) while a dielectric waveguide consists of multiple dielectrics. • Normal operating mode is the COperating modes are TE or TM TEM or quasi-TEM mode (can modes (cannot support a TEM support TE and TM modes but mode). these modes are typically undesirable). • No cutoff frequency for the TEM CMust operate the waveguide at a mode. Transmission lines can frequency above the respective transmit signals from DC up to TE or TM mode cutoff frequency high frequency. for that mode to propagate. • Significant signal attenuation at CLower signal attenuation at high high frequencies due to frequencies than transmission conductor and dielectric losses. lines. • Small cross-section transmission CMetal waveguides can transmit lines (like coaxial cables) can high power levels.
    [Show full text]
  • Simple Transmission Line Representation of Tesla Coil and Tesla's Wave Propagation Concept
    November, 2006 Microwave Review Simple Transmission Line Representation of Tesla Coil and Tesla’s Wave Propagation Concept Zoran Blažević1, Dragan Poljak2, Mario Cvetković3 Abstract – In this paper, we shall propose a simple Many researchers of today in explanation of the Tesla's transmission line model for a Tesla coil that includes distributed concept assign the term “Hertzian wave” to the term “space voltage source along secondary of Tesla transformer, which is wave” and “current wave” to “surface wave” justifying this able to explain some of the effects that cannot be fully simulated with the difference in terminology of 19th/20th century and by the model with the lumped source. Also, a transmission line nowadays. However Tesla’s surface waves, if we could be representation of “propagation of current-through-Earth wave” will be presented upon. allowed to call them that way, have some different properties from, for instance, the Zenneck surface waves. Namely, they Keywords – Transmission line modelling, Tesla’s resonator, travel with a variable velocity in range from infinity at the Magnifying effect, Tesla’s current wave. electric poles down to the light speed at the electric equator as a consequence of the charges pushed by the resonator to I. INTRODUCTION tunnel through Earth along her diameter with speed equal or close to the speed of light in vacuum. And its effectiveness Ever since Dr. Nikola Tesla started to unfold the principles does not depend on polarization [1]. Regardless of the fact of electricity through the results of his legendary experiments that science considers this theory as fallacious, which is well and propositions, there is a level talk about it that actually explained in [6], it is also true that his conclusions are based never ceased, even to these days.
    [Show full text]
  • DETERMINATION of PROPAGATION CONSTANTS and MATERIAL DATA from WAVEGUIDE MEASUREMENTS D. Sjöberg Department of Electrical and In
    Progress In Electromagnetics Research B, Vol. 12, 163–182, 2009 DETERMINATION OF PROPAGATION CONSTANTS AND MATERIAL DATA FROM WAVEGUIDE MEASUREMENTS D. Sj¨oberg Department of Electrical and Information Technology Lund University P. O. Box 118, 221 00 Lund, Sweden Abstract—This paper presents an analysis with the aim of characterizing the electromagnetic properties of an arbitrary linear, bianisotropic material inside a metallic waveguide. The result is that if the number of propagating modes is the same inside and outside the material under test, it is possible to determine the propagation constants of the modes inside the material by using scattering data from two samples with different lengths. Some information can also be obtained on the cross-sectional shape of the modes, but it remains an open question if this information can be used to characterize the material. The method is illustrated by numerical examples, determining the complex permittivity for lossy isotropic and anisotropic materials. 1. INTRODUCTION In order to obtain a well controlled environment for making measurements of electromagnetic properties, it is common to do the measurements in a metallic cavity or waveguide. The geometrical constraints of the waveguide walls impose dispersive characteristics on the propagation of electromagnetic waves, i.e., the wavelength of the propagating wave depends on frequency in a nonlinear manner. In order to correctly interpret the measurements, it is necessary to provide a suitable characterization of the waves inside the waveguide. This is well known for isotropic materials [2, 3, 12, 20], bi-isotropic (chiral) materials [9, 17], and even anisotropic materials where an optical axis is along the waveguide axis [1, 10, 15, 16], but for general bianisotropic Corresponding author: D.
    [Show full text]
  • Design of an Ultrahigh Birefringence Photonic Crystal Fiber with Large Nonlinearity Using All Circular Air Holes for a Fiber-Optic Transmission System
    hv photonics Article Design of an Ultrahigh Birefringence Photonic Crystal Fiber with Large Nonlinearity Using All Circular Air Holes for a Fiber-Optic Transmission System Shovasis Kumar Biswas 1,* ID , S. M. Rakibul Islam 1, Md. Rubayet Islam 1, Mohammad Mahmudul Alam Mia 2, Summit Sayem 1 and Feroz Ahmed 1 ID 1 Department of Electrical and Electronic Engineering, Independent University Bangladesh, Dhaka 1229, Bangladesh; [email protected] (S.M.R.I.); [email protected] (M.R.I.); [email protected] (S.S.); [email protected] (F.A.) 2 Department of Electronics & Communication Engineering, Sylhet International University, Sylhet 3100, Bangladesh; [email protected] * Correspondence: [email protected]; Tel.: +880-1717-827-665 Received: 11 August 2018; Accepted: 29 August 2018; Published: 5 September 2018 Abstract: This paper proposes a hexagonal photonic crystal fiber (H-PCF) structure with all circular air holes in order to simultaneously achieve ultrahigh birefringence and high nonlinearity. The H-PCF design consists of an asymmetric core region, where one air hole is a reduced diameter and the air hole in its opposite vertex is omitted. The light-guiding properties of the proposed H-PCF structure were studied using the full-vector finite element method (FEM) with a circular perfectly matched layer (PML). The simulation results showed that the proposed H-PCF exhibits an ultrahigh birefringence of 3.87 × 10−2, a negative dispersion coefficient of −753.2 ps/(nm km), and a nonlinear coefficient of 96.51 W−1 km−1 at an excitation wavelength of 1550 nm. The major advantage of our H-PCF design is that it provides these desirable modal properties without using any non-circular air holes in the core and cladding region, thus making the fiber fabrication process much easier.
    [Show full text]
  • Recent Advances in Perfectly Matched Layers in Finite Element Applications
    Turk J Elec Engin, VOL.16, NO.1 2008, c TUB¨ ITAK˙ Recent Advances in Perfectly Matched Layers in Finite Element Applications Ozlem¨ OZG¨ UN¨ and Mustafa KUZUOGLU˘ Department of Electrical and Electronics Engineering, Middle East Technical University 06531, Ankara-TURKEY Abstract We present a comparative evaluation of two novel and practical perfectly matched layer (PML) implementations to the problem of mesh truncation in the finite element method (FEM): locally-conformal PML,and multi-center PML techniques. The most distinguished feature of these methods is the simplicity and flexibility to design conformal PMLs over challenging geometries,especially those with curvature discontinuities,in a straightforward way without using artificial absorbers. These methods are based on specially- and locally-defined complex coordinate transformations inside the PML region. They can easily be implemented in a conventional FEM by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. After overviewing the theoretical bases of these methods,we present some numerical results in the context of two- and three-dimensional electromagnetic radiation/scattering problems. Key Words: Finite element method (FEM); perfectly matched layer (PML); locally-conformal PML; multi-center PML; complex coordinate stretching; electromagnetic scattering. 1. Introduction Perfectly matched layers (PMLs) have been popular in the finite element method (FEM) and the finite difference time domain method (FDTD) for solving the domain truncation problem in electromagnetic radiation and scattering problems.The PML approach is based on the truncation of computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence.The main advantage of the PML is its close proximity and conformity to the surface of the radiator or scatterer, implying that the white space can be minimized.
    [Show full text]