Modelling of Ridge Waveguide Bends for Sensor Applications
|
|
- Piers Ryan
- 6 years ago
- Views:
Transcription
1 Modelling of Ridge Waveguide Bends for Sensor Applications Wilfrid Pascher FernUniversität, Hagen, Germany n 1 n2 n 3 R
2 Modelling of Ridge Waveguide Bends for Sensor Applications Wilfrid Pascher FernUniversität, Hagen, Germany n 1 n2 n 3 Radiation losses Evanescent field of a rib waveguide are precisely modelled R by the Method of Lines
3 Why employ the Method of Lines? The MoL is a semianalytic approach
4 Why employ the Method of Lines? The MoL is a semianalytic approach Analytic solution in one coordinate direction (perpendicular to the layers)
5 Why employ the Method of Lines? The MoL is a semianalytic approach Analytic solution in one coordinate direction (perpendicular to the layers) Discretization in the other direction(s) (with Finite Differences) = 3D problem 2D discretization
6 Why employ the Method of Lines? The MoL is a semianalytic approach Analytic solution in one coordinate direction (perpendicular to the layers) Discretization in the other direction(s) (with Finite Differences) = 3D problem 2D discretization For reasons of technology, waveguide structures are multilayered (e.g., planar waveguides) cascaded (e.g., waveguide circuits)
7 Advantages and Disadvantages + precise modeling + low memory and computing time
8 Advantages and Disadvantages + precise modeling + low memory and computing time reduced flexibility = different geometries require new algorithms = extension to hybrid methods
9 Transition to Cylindrical Coordinates n 1 n1 n 2 n 2 R
10 Transition to Cylindrical Coordinates z Propagation ϕ r P o p a a g t on i n 1 n1 n 2 n 2 R
11 Transition to Cylindrical Coordinates z Propagation ϕ r P o p a a g t on i n 1 n1 n 2 n 2 R Propagation z ϕ exp( jβz) exp( jνϕ)
12 Discretization of Straight / Bent Waveguides
13 Discretization of Straight / Bent Waveguides
14 Discretization of Straight / Bent Waveguides Discretization x Analytic Solution r Analytic Solution y z Discretization ψ ψ e h ψ h ψ e
15 Discretization of Straight / Bent Waveguides Discretization x Analytic Solution r Analytic Solution y z Discretization ψ ψ e h ψ h ψ e Discretization x z P x P z
16 Discretization of Straight / Bent Waveguides Discretization x Analytic Solution r Analytic Solution y z Discretization ψ ψ e h ψ h ψ e Discretization x z Analytic solution y r P x P z sin(k y y) J ν ( λr)
17 The Method of Lines (MoL) for circular bends in waveguides 1. Transition cartesian cylindrical (x, y, z) (z, r, ϕ)
18 The Method of Lines (MoL) for circular bends in waveguides 1. Transition cartesian cylindrical (x, y, z) (z, r, ϕ) 2. Discretization of the wave equation 2 z 2 P z (3. Transformation to diagonal form) diag ( λ 2 k )
19 The Method of Lines (MoL) for circular bends in waveguides 1. Transition cartesian cylindrical (x, y, z) (z, r, ϕ) 2. Discretization of the wave equation 2 z 2 P z (3. Transformation to diagonal form) diag ( λ 2 k ) ) 4. Solution of the wave equation J ν ( λk r +...
20 The Method of Lines (MoL) for circular bends in waveguides 1. Transition cartesian cylindrical (x, y, z) (z, r, ϕ) 2. Discretization of the wave equation 2 z 2 P z (3. Transformation to diagonal form) diag ( λ 2 k ) ) 4. Solution of the wave equation J ν ( λk r Field computation (6. Inverse transformation)
21 The Method of Lines (MoL) for circular bends in waveguides 1. Transition cartesian cylindrical (x, y, z) (z, r, ϕ) 2. Discretization of the wave equation 2 z 2 P z (3. Transformation to diagonal form) diag ( λ 2 k ) ) 4. Solution of the wave equation J ν ( λk r Field computation (6. Inverse transformation) 7. Characteristic equation in C radiation loss L Im(n eff )
22 Wave equations in coordinate free form Vector MoL with two potentials Π e, Π h accurate fulfillment of the continuity conditions for all field components
23 Wave equations in coordinate free form Vector MoL with two potentials Π e, Π h accurate fulfillment of the continuity conditions for all field components coordinate free approach: a) Helmholtz equation for Π h { + ε r (z)k 2 0 } Π h = 0
24 Wave equations in coordinate free form Vector MoL with two potentials Π e, Π h accurate fulfillment of the continuity conditions for all field components coordinate free approach: a) Helmholtz equation for Π h { + ε r (z)k 2 0 } Π h = 0 b) Sturm-Liouville differential equation for Π e { + ε r (z)k ε r (z) grad ε r(z) div } Π e = 0
25 Wave equations in cylindrical coordinates Potentials with one component in z direction only Π e,h = k 2 0 exp( jνϕ) ψ e,h a z order proportional to effective model index ν = n eff R using normalized coordinates: e.g. R = k 0 R
26 Wave equations in cylindrical coordinates Potentials with one component in z direction only Π e,h = k 2 0 exp( jνϕ) ψ e,h a z order proportional to effective model index ν = n eff R using normalized coordinates: e.g. R = k 0 R Consideration of the radiation losses n eff, ν complex no artificial increase in the guiding
27 Discretization of the wave equation in the cartesian z direction Partial differential equations in cylindrical coordinates { ( 1 r ) } ν2 r r r r 2 + ε r(z) + 2 z 2 ψ h = 0 { 1 r r ( r ) ν2 r r 2 + ε r(z) + ε r (z) z ( 1 ε r (z) ) } ψ e = 0 z
28 Discretization of the wave equation in the cartesian z direction Partial differential equations in cylindrical coordinates { ( 1 r ) } ν2 r r r r 2 + ε r(z) + 2 z 2 ψ h = 0 { 1 r r ( r ) ν2 r r 2 + ε r(z) + ε r (z) z ( 1 ε r (z) ) } ψ e = 0 z Potentials and dielectric constants continuous discretized ψ e, ψ h Ψ e, Ψ h (column vector) ε r (z) ε e, ε h (diagonal matrix)
29 Discretization of the wave equation in the cartesian z direction Differential operators difference operators ε r (z) z ( 1 ε r (z) 2 z 2 ψ h P zh Ψ h (tridiagonal) ) ψ e Pze ε Ψ e (tridiagonal) z
30 Discretization of the wave equation in the cartesian z direction Differential operators difference operators ε r (z) z ( 1 ε r (z) 2 z 2 ψ h P zh Ψ h (tridiagonal) ) ψ e Pze ε Ψ e (tridiagonal) z Coupled ordinary differential equations 1 r d dr ( r d ) I ν2 dr r 2 I + ε P ε ze e,h P zh } {{ } (tridiagonal) Ψ e,h = 0
31 Transformation to diagonal form with T ε e 1 (tridiagonal) {}}{ (ε e P ε ze) T ε e = λ 2 e = diag ( λ 2 e,k )
32 Transformation to diagonal form with (tridiagonal) {}}{ Te ε 1 (ε e Pze) ε Te ε = λ 2 e = diag ( λ 2 e,k ) Decoupled ordinary differential equations { ( 1 d r d ) } I ν2 r dr dr r 2 I + λ 2 e Ψ e = 0 with the transformed potential Ψ e = T ε e 1 Ψ e
33 Transformation to diagonal form with (tridiagonal) {}}{ Te ε 1 (ε e Pze) ε Te ε = λ 2 e = diag ( λ 2 e,k ) Decoupled ordinary differential equations { ( 1 d r d ) } I ν2 r dr dr r 2 I + λ 2 e Ψ e = 0 with the transformed potential Ψ e = T ε e 1 Ψ e completely analogous for Ψ h = T h 1 Ψ h
34 Solution of the Bessel differential equation 1 r d dr ( r d ) dr + ( λ 2 e,k ν2 r 2 ) Ψ e,k = 0 with ν = n eff R b for one component Ψ e,k of the transformed potential
35 Solution of the Bessel differential equation
36 Solution of the Bessel differential equation Solution in three regions
37 Solution of the Bessel differential equation Solution in three regions 1 finite for r 0 Ψ e,k = A k J ν ( λ e,k r)
38 Solution of the Bessel differential equation Solution in three regions 1 finite for r 0 Ψ e,k = A k J ν ( λ e,k r) 2 Ψ e,k = B k J ν ( λ e,k r) + C k Y ν ( λ e,k r)
39 Solution of the Bessel differential equation Solution in three regions 1 finite for r 0 Ψ e,k = A k J ν ( λ e,k r) 2 Ψ e,k = B k J ν ( λ e,k r) + C k Y ν ( λ e,k r) 3 radiation condition Ψ e,k = D k H (2) ν ( λ e,k r)
40 Radial derivation matrices Γ Transmission line equation A = Γ A with r B B Γ = j λ p ν r ν W B W A q ν
41 Radial derivation matrices Γ Transmission line equation A = Γ A with r B B Γ = j λ p ν r ν W B W A q ν The radial derivation matrix Γ is computed from the cross products and the Wronskian of the Bessel functions p ν = J ν ( r A )Y ν ( r B ) J ν ( r B )Y ν ( r A ) q ν = J ν ( r A )Y ν( r B ) J ν( r B )Y ν ( r A ) r ν = J ν( r A )Y ν ( r B ) J ν ( r B )Y ν( r A ) W A,B = 2 π r A,B with r A,B = j λr A,B
42 Programming of the cylinder functions The problem 1. very high complex order: ν j with small imaginary part
43 Programming of the cylinder functions The problem 1. very high complex order: ν j with small imaginary part 2. order nearly equal to argument: ν λ 1 r for the first eigenvalue
44 Programming of the cylinder functions The problem 1. very high complex order: ν j with small imaginary part 2. order nearly equal to argument: ν λ 1 r for the first eigenvalue The solution 1, 2 Uniform asymptotic series (high argument and order)
45 Programming of the cylinder functions The problem 1. very high complex order: ν j with small imaginary part 2. order nearly equal to argument: ν λ 1 r for the first eigenvalue The solution 1, 2 Uniform asymptotic series (high argument and order) Alternative for the ring regions Multiplication formulas for cross products
46 Uniform asymptotic expansions Bessel functions, e.g. Abramowitz-Stegun (9.3.35) J ν (νy) ( 4ζ 1 y 2 ) 1/4 Ai(ν 2/3 ζ) ν 1/3 i=0 a i (ζ) ν 2i + Ai (ν 2/3 ζ) ν 5/3 i=0 b i (ζ) ν 2i ν with 2 3 ζ3/2 = log y 2 y 1 y 2
47 Uniform asymptotic expansions Bessel functions, e.g. Abramowitz-Stegun (9.3.35) J ν (νy) ν ( 4ζ 1 y 2 ) 1/4 Ai(ν 2/3 ζ) ν 1/3 i=0 a i (ζ) ν 2i + Ai (ν 2/3 ζ) ν 5/3 with 2 3 ζ3/2 = log y 2 y Coefficients a i (ζ), b i (ζ) descending and independent of ν i=0 b i (ζ) ν 2i 1 y 2
48 Uniform asymptotic expansions Bessel functions, e.g. Abramowitz-Stegun (9.3.35) J ν (νy) ν ( 4ζ 1 y 2 ) 1/4 Ai(ν 2/3 ζ) ν 1/3 i=0 a i (ζ) ν 2i + Ai (ν 2/3 ζ) ν 5/3 with 2 3 ζ3/2 = log y 2 y Coefficients a i (ζ), b i (ζ) descending and independent of ν Series terms decrease by 1 ν Series can be truncated already after i = 2 for double precision arithmetic i=0 b i (ζ) ν 2i 1 y 2
49 Uniform asymptotic expansions Bessel functions, e.g. Abramowitz-Stegun (9.3.35) J ν (νy) ν ( 4ζ 1 y 2 ) 1/4 Ai(ν 2/3 ζ) ν 1/3 i=0 a i (ζ) ν 2i + Ai (ν 2/3 ζ) ν 5/3 with 2 3 ζ3/2 = log y 2 y Coefficients a i (ζ), b i (ζ) descending and independent of ν Series terms decrease by 1 ν Series can be truncated already after i = 2 for double precision arithmetic i=0 b i (ζ) ν 2i 1 y 2 Y ν (νy), J ν(νy), Y ν(νy) computed analogously
50 Multiplication theorem for direct calculation of cross products Bessel function of outer radius r B = µ r A C ν ( r B ) = C ν (µ r A ) = µ ν i=0 (δ r A ) i i! C ν i ( r A ) with δ = (µ 2 1)/2 Abramowitz-Stegun (9.1.74)
51 Multiplication theorem for direct calculation of cross products Bessel function of outer radius r B = µ r A C ν ( r B ) = C ν (µ r A ) = µ ν i=0 (δ r A ) i i! C ν i ( r A ) with δ = (µ 2 1)/2 Abramowitz-Stegun (9.1.74) yields series for cross products p ν, q ν p ν = (δ r A ) i µ ν i! q ν i=0 P ν,i ( r A ) Q ν,i ( r A )
52 Multiplication theorem for direct calculation of cross products Bessel function of outer radius r B = µ r A C ν ( r B ) = C ν (µ r A ) = µ ν i=0 (δ r A ) i i! C ν i ( r A ) with δ = (µ 2 1)/2 Abramowitz-Stegun (9.1.74) yields series for cross products p ν, q ν p ν = (δ r A ) i µ ν i! q ν i=0 P ν,i ( r A ) Q ν,i ( r A ) P, Q are computed from a three term recurrence relation Other cross products r ν, s ν are determined from p ν, q ν
53 Multiplication theorem for direct calculation of cross products Bessel function of outer radius r B = µ r A C ν ( r B ) = C ν (µ r A ) = µ ν i=0 (δ r A ) i i! C ν i ( r A ) with δ = (µ 2 1)/2 Abramowitz-Stegun (9.1.74) yields series for cross products p ν, q ν p ν = (δ r A ) i µ ν i! q ν i=0 P ν,i ( r A ) Q ν,i ( r A ) P, Q are computed from a three term recurrence relation Other cross products r ν, s ν are determined from p ν, q ν = Two independent procedures for checking
54 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B
55 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B Transfer of the fields from cylinder r = r A to r = r B H A H B = ỹ1 ỹ 2 ỹ 2 ỹ 1+ Ẽ A ẼB
56 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B Transfer of the fields from cylinder r = r A to r = r B H A H B = ỹ1 ỹ 2 ỹ 2 ỹ 1+ Ẽ A ẼB yields HB = ỸBẼB
57 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B Transfer of the fields from cylinder r = r A to r = r B H A H B = ỹ1 ỹ 2 ỹ 2 ỹ 1+ Ẽ A ẼB yields HB = ỸBẼB with the recurrence Ỹ B = ỹ 2 (ỹ 1 ỸA) 1 ỹ2 ỹ 1+
58 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B Transfer of the fields from cylinder r = r A to r = r B H A H B = ỹ1 ỹ 2 ỹ 2 ỹ 1+ Ẽ A ẼB yields HB = ỸBẼB with the recurrence Ỹ B = ỹ 2 (ỹ 1 ỸA) 1 ỹ2 ỹ 1+ In cartesian coordinates matrices ỹ 1,2 depend only on the layer thickness
59 Field computation in every subregion separately Matching on the cylinders r = r A and r = r B Transfer of the fields from cylinder r = r A to r = r B H A H B = ỹ1 ỹ 2 ỹ 2 ỹ 1+ Ẽ A ẼB yields HB = ỸBẼB with the recurrence Ỹ B = ỹ 2 (ỹ 1 ỸA) 1 ỹ2 ỹ 1+ In cartesian coordinates matrices ỹ 1,2 depend only on the layer thickness now ỹ 1,2 depend on the radii r A, r B now ỹ 1,2 are computed using cylinder functions
60 Determination of the radiation losses Inverse transformation to spatial domain characteristic equation Z(n eff ) H A = 0
61 Determination of the radiation losses Inverse transformation to spatial domain characteristic equation Z(n eff ) H A = 0 Solution det(z(n eff )) = 0 by a search for zeros in the complex plane C
62 Determination of the radiation losses Inverse transformation to spatial domain characteristic equation Z(n eff ) H A = 0 Solution det(z(n eff )) = 0 by a search for zeros in the complex plane C Radiation losses (db/90 ): L = Im(n eff ) R π 10 ln 10
63 Geometry of a Bent Rib Waveguide Sensor gas w t Si Si 3 N 4 Si O 2 n 1 n2 n 3 h 3 h 2 R OIWS108A
64 Geometry of a Bent Rib Waveguide Sensor gas w t Si 3 N 4 Si O 2 n 1 n2 n 3 h 2 h 3 w = 3 5µm n 1 = t = µm n 2 = h 2 = 5µm Si n 3 = 3.5 R OIWS108A
65 Sensitivity Sensitivity in % w [µm] Sensitivity depending on thickness t for w = 3.0 µm
66 Distribution of the Radial Electric Field z [µm] dr [µm] Fundamental mode, w = 3.0 µm
67 Distribution of the Radial Electric Field z [µm] dr [µm] First higher order mode, w = 3.50 µm
68 Distribution of the Radial Electric Field z [µm] dr [µm] First higher order mode, w = 3.13 µm
69 Conclusion: Advantages and Disadvantages The presented model yields most accurate results for propagation constant radiation loss
70 Conclusion: Advantages and Disadvantages The presented model yields most accurate results for propagation constant sensitivity radiation loss evanescent field
71 Conclusion: Advantages and Disadvantages The presented model yields most accurate results for propagation constant sensitivity radiation loss evanescent field Comparison at the TU Delft for a polarization converter: MoL is superior to the Effective Index Method (EIM) and Finite Element Method (FEM)
72 Conclusion: Advantages and Disadvantages The presented model yields most accurate results for propagation constant sensitivity radiation loss evanescent field Comparison at the TU Delft for a polarization converter: MoL is superior to the Effective Index Method (EIM) and Finite Element Method (FEM) Disadvantage only applicable to strictly rotational structures Extension: Combination with other numerical methods
Cylindrical Dielectric Waveguides
03/02/2017 Cylindrical Dielectric Waveguides Integrated Optics Prof. Elias N. Glytsis School of Electrical & Computer Engineering National Technical University of Athens Geometry of a Single Core Layer
More informationAnalytical modal analysis of bent slot waveguides
Analytical modal analysis of bent slot waveguides Kirankumar R. Hiremath Department of Mathematics, Institute for Scientific Computing and Mathematical Modeling, University of Karlsruhe, Engesserstrasse
More informationElectromagnetic Theory for Microwaves and Optoelectronics
Keqian Zhang Dejie Li Electromagnetic Theory for Microwaves and Optoelectronics Second Edition With 280 Figures and 13 Tables 4u Springer Basic Electromagnetic Theory 1 1.1 Maxwell's Equations 1 1.1.1
More informationElectromagnetic Theory for Microwaves and Optoelectronics
Keqian Zhang Dejie Li Electromagnetic Theory for Microwaves and Optoelectronics Translated by authors With 259 Figures Springer Contents 1 Basic Electromagnetic Theory 1 1.1 Maxwell's Equations 1 1.1.1
More informationCut-on, cut-off transition of sound in slowly varying flow ducts
Cut-on, cut-off transition of sound in slowly varying flow ducts Sjoerd W. Rienstra 19-3-1 Department of Mathematics and Computing Science, Eindhoven University of Technology, The Netherlands, s.w.rienstra@tue.nl
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationCOMSOL Design Tool: Simulations of Optical Components Week 6: Waveguides and propagation S matrix
COMSOL Design Tool: Simulations of Optical Components Week 6: Waveguides and propagation S matrix Nikola Dordevic and Yannick Salamin 30.10.2017 1 Content Revision Wave Propagation Losses Wave Propagation
More informationComparative Analysis of Techniques for Source Radiation in Cylindrical EBG with and without Periodic Discontinuities
1398 Progress In Electromagnetics Research Symposium Abstracts, St Petersburg, Russia, 22 25 May 2017 Comparative Analysis of Techniques for Source Radiation in Cylindrical EBG with and without Periodic
More informationAnalysis of Characteristics of Coplanar Waveguide with Finite Ground-planes by the Method of Lines
PIERS ONLINE, VOL. 6, NO. 1, 21 46 Analysis of Characteristics of Coplanar Waveguide with Finite Ground-planes by the Method of Lines Min Wang, Bo Gao, Yu Tian, and Ling Tong College of Automation Engineering,
More information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationFundamentals of fiber waveguide modes
SMR 189 - Winter College on Fibre Optics, Fibre Lasers and Sensors 1-3 February 007 Fundamentals of fiber waveguide modes (second part) K. Thyagarajan Physics Department IIT Delhi New Delhi, India Fundamentals
More informationA RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES
Progress In Electromagnetics Research, PIER 22, 197 212, 1999 A RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES M. Akbari, M. Shahabadi, and K. Schünemann Arbeitsbereich Hochfrequenztechnik Technische
More informationDepartment of Physics Preliminary Exam January 2 5, 2013
Department of Physics Preliminary Exam January 2 5, 2013 Day 2: Electricity, Magnetism and Optics Thursday, January 3, 2013 9:00 a.m. 12:00 p.m. Instructions: 1. Write the answer to each question on a
More informationGuided and Leaky Modes of Planar Waveguides: Computation via High Order Finite Elements and Iterative Methods
PERS ONLNE, VOL. 6, NO. 7, 669 Guided and Leaky Modes of Planar Waveguides: Computation via High Order Finite Elements and terative Methods D. Stowell and J. Tausch Southern Methodist University, USA Abstract
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationA comparative study of higher order Bragg waveguide gratings using coupled mode theory and mode expansion modeling. H. Wenzel and R.
A comparative study of higher order Bragg waveguide gratings using coupled mode theory and mode expansion modeling H. Wenzel and R. Güther Motivation Widespread application of Bragg waveguide gratings
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
More informationDrawbacks in Metallic Waveguides. α f. Skin effect. Surface roughness Mono-mode operation Narrow band in metallic rectangular waveguide
Drawbacks in Metallic Waveguides Skin effect α f Surface roughness Mono-mode operation 1 Narrow band in metallic rectangular waveguide 2 α f α 3 f 2 to f 5 2 Types of Dielectric Waveguide ε r2 (a) Circular
More informationBessel function - Wikipedia, the free encyclopedia
Bessel function - Wikipedia, the free encyclopedia Bessel function Page 1 of 9 From Wikipedia, the free encyclopedia In mathematics, Bessel functions, first defined by the mathematician Daniel Bernoulli
More informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
More informationUndulator Radiation Inside a Dielectric Waveguide
Undulator Radiation Inside a Dielectric Waveguide A.S. Kotanjyan Department of Physics, Yerevan State University Yerevan, Armenia Content Motivation On features of the radiation from an electron moving
More informationON THE HYBRID FIELD PATTERNS OF HELICAL CLAD DIELECTRIC OPTICAL FIBERS
Progress In Electromagnetics Research, PIER 91, 69 84, 2009 ON THE HYBRID FIELD PATTERNS OF HELICAL CLAD DIELECTRIC OPTICAL FIBERS A. H. B. M. Safie and P. K. Choudhury Faculty of Engineering Multimedia
More informationLecture 3 Fiber Optical Communication Lecture 3, Slide 1
Lecture 3 Optical fibers as waveguides Maxwell s equations The wave equation Fiber modes Phase velocity, group velocity Dispersion Fiber Optical Communication Lecture 3, Slide 1 Maxwell s equations in
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationINGENIERÍA EN NANOTECNOLOGÍA
ETAPA DISCIPLINARIA TAREAS 385 TEORÍA ELECTROMAGNÉTICA Prof. E. Efren García G. Ensenada, B.C. México 206 Tarea. Two uniform line charges of ρ l = 4 nc/m each are parallel to the z axis at x = 0, y = ±4
More informationMathematical and numerical techniques for open periodic waveguides
Mathematical and numerical techniques for open periodic waveguides Johannes Tausch Department of Mathematics Southern Methodist University Dallas, TX 75275 tausch@smu.edu Abstract The propagation of electromagnetic
More informationAdvanced Electrodynamics Exercise 11 Guides
Advanced Electrodynamics Exercise 11 Guides Here we will calculate in a very general manner the modes of light in a waveguide with perfect conductor boundary-conditions. Our derivations are widely independent
More informationCartesian Coordinates
Cartesian Coordinates Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Cartesian Coordinates Outline Outline Separation of Variables Away from sources,
More informationMicrowave Phase Shift Using Ferrite Filled Waveguide Below Cutoff
Microwave Phase Shift Using Ferrite Filled Waveguide Below Cutoff CHARLES R. BOYD, JR. Microwave Applications Group, Santa Maria, California, U. S. A. ABSTRACT Unlike conventional waveguides, lossless
More informationSurvival Guide to Bessel Functions
Survival Guide to Bessel Functions December, 13 1 The Problem (Original by Mike Herman; edits and additions by Paul A.) For cylindrical boundary conditions, Laplace s equation is: [ 1 s Φ ] + 1 Φ s s s
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationProperties of axial surface waves along dielectrically coated conducting cylinders
Adv. Radio Sci., 10, 79 84, 2012 doi:10.5194/ars-10-79-2012 Author(s) 2012. CC Attribution 3.0 License. Advances in Radio Science Properties of axial surface waves along dielectrically coated conducting
More informationUsing quasi-guided modes for modeling the transfer behavior of bent dielectric slab waveguides
Adv. Radio Sci., 8, 19 26, 2010 doi:10.5194/ars-8-19-2010 Authors 2010. CC Attribution 3.0 License. Advances in Radio Science Using quasi-guided modes for modeling the transfer behavior of bent dielectric
More informationIntroduction to optical waveguide modes
Chap. Introduction to optical waveguide modes PHILIPPE LALANNE (IOGS nd année) Chapter Introduction to optical waveguide modes The optical waveguide is the fundamental element that interconnects the various
More informationModified Bessel functions : Iα, Kα
Modified Bessel functions : Iα, Kα The Bessel functions are valid even for complex arguments x, and an important special case is that of a purely imaginary argument. In this case, the solutions to the
More informationTASK A. TRANSMISSION LINE AND DISCONTINUITIES
TASK A. TRANSMISSION LINE AND DISCONTINUITIES Task A. Transmission Line and Discontinuities... 1 A.I. TEM Transmission Line... A.I.1. Circuit Representation of a Uniform Transmission Line... A.I.. Time
More informationECE 6341 Spring 2016 HW 2
ECE 6341 Spring 216 HW 2 Assigned problems: 1-6 9-11 13-15 1) Assume that a TEN models a layered structure where the direction (the direction perpendicular to the layers) is the direction that the transmission
More informationFiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:
Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total
More informationUniversity of Warwick institutional repository:
University of Warwick institutional repository: http://go.warwick.ac.uk/wrap This paper is made available online in accordance with publisher policies. Please scroll down to view the document itself. Please
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More information1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationAzimuthally polarized surface plasmons as effective terahertz waveguides
Azimuthally polarized surface plasmons as effective terahertz waveguides Qing Cao and Jürgen Jahns Optische Nachrichtentechnik, FernUniversität Hagen, Universitätsstrasse 27/PRG, D-58084 Hagen, Germany
More informationA STUDY OF THE SOUND PROPAGATION IN THE ROUNDED SECTION OF A PIPE
A STUDY OF THE SOUND PROPAGATION IN THE ROUNDED SECTION OF A PIPE Professor habil. dr. P. ILGAKOJIS, Lecturer dr. E. JOTAUTIENE - Lithuanian University of Agriculture Abstract: In this paper, the sound
More informationA simple introduction to the method of lines
M. N. O. Sadiku1 and C. N. Obiozor2 1Department of Electrical and Computer Engineering, Temple University, Philadelphia, U.S.A. E-mail: sadiku@astro.ocis.temple.edu 2Department of Electrical Engineering,
More informationSurface plasmon polariton propagation around bends at a metal-dielectric interface
Surface plasmon polariton propagation around bends at a metal-dielectric interface Keisuke Hasegawa, Jens U. Nöckel and Miriam Deutsch Oregon Center for Optics, 1274 University of Oregon, Eugene, OR 97403-1274
More information! #! % && ( ) ) +++,. # /0 % 1 /21/ 3 && & 44&, &&7 4/ 00
! #! % && ( ) ) +++,. # /0 % 1 /21/ 3 &&4 2 05 6. 4& 44&, &&7 4/ 00 8 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 56, NO. 2, FEBRUARY 2008 345 Moment Method Analysis of an Archimedean Spiral Printed
More informationChapter 2 Basic Optics
Chapter Basic Optics.1 Introduction In this chapter we will discuss the basic concepts associated with polarization, diffraction, and interference of a light wave. The concepts developed in this chapter
More informationRing-shaped crack propagation in a cylinder under nonsteady cooling
High Performance Structures and Materials III 5 Ring-shaped crack propagation in a cylinder under nonsteady cooling V. A. Zhornik, Yu. A. Prokopenko, A. A. Rybinskaya & P. A. Savochka Department of Theoretical
More informationHIGH VOLTAGE TECHNIQUES Basic Electrode Systems (3)
HIGH VOLTAGE TECHNIQES Basic Electrode Systems (3) Assistant Professor Suna BOLAT KRÖGER Eastern Mediterranean niversity Department of Electric & Electronic Engineering 1 Basic electrode systems Different
More informationUSPAS Accelerator Physics 2017 University of California, Davis
USPAS Accelerator Physics 2017 University of California, Davis Chapter 9: RF Cavities and RF Linear Accelerators Todd Satogata (Jefferson Lab) / satogata@jlab.org Randika Gamage (ODU) / bgama002@odu.edu
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationPlanar QED in Magnetic or Electric Solitonic Backgrounds
Planar QED in Magnetic or Electric Solitonic Backgrounds Gerald Dunne Physics Department University of Connecticut Storrs, CT 0669 USA Abstract We review the resolvent technique for computing the effective
More information4.1 Solving Maxwell s equations in linear, isotropic, homogeneous, circularly symmetric media
Maxwell s Equations in Cylindrical Coordinates Masud Mansuripur College of Optical Sciences, The University of Arizona, Tucson, Arizona 8571 Introduction In many systems of practical interest, electromagnetic
More informationSIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(0,n) MODES*
25th IEEE International Conference on Plasma Science Raleigh, North Carolina June 1-4, 1998 SIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(,n) MODES* Ron L. Kinder and Mark J.
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationSound Transmission in an Extended Tube Resonator
2016 Published in 4th International Symposium on Innovative Technologies in Engineering and Science 3-5 November 2016 (ISITES2016 Alanya/Antalya - Turkey) Sound Transmission in an Extended Tube Resonator
More informationQuasi-Optical Design and Analysis (MBI) Créidhe O Sullivan, J.Anthony Murphy, Marcin Gradziel, Neil Trappe, Tully Peacocke & graduate students
Quasi-Optical Design and Analysis (MBI) Créidhe O Sullivan, J.Anthony Murphy, Marcin Gradziel, Neil Trappe, Tully Peacocke & graduate students Outline Corrugated Horns Analysis Techniques MBI/MODAL 2 Analysis
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationELECTRICAL CALCULATIONS SOFTWARE. R.Ubic 1996
Page 1 of 5 ELECTRICAL CALCULATIONS SOFTWARE R.Ubic 1996 These programmes are designed to calculate the dielectric constant of a microwave resonator from its geometry, the geometry of the test cavity,
More informationRezonanse typu spin-orbit w meta-materiałowych elementach nanofotoniki Spin-orbit resonances in meta-material elements of nanophotonics
Rezonanse typu spin-orbit w meta-materiałowych elementach nanofotoniki Spin-orbit resonances in meta-material elements of nanophotonics Wojciech Nasalski Zespół Badawczy Nanofotoniki Instytut Podstawowych
More informationAnalysis of Metamaterial Cloaks Using Circular Split Ring Resonator Structures
Copyright 216 Tech Science Press CMC, Vol.53, No.3, pp.132-14, 216 Analysis of Metamaterial Cloaks Using Circular Split Ring Resonator Structures Susan Thomas 1 and Dr. Balamati Choudhury 2 Abstract A
More informationMUTUAL INDUCTANCE FOR AN EXPLICITLY FINITE NUMBER OF TURNS
Progress In Electromagnetics Research B, Vol. 28, 273 287, 211 MUTUAL INDUCTANCE FOR AN EXPLICITLY FINITE NUMBER OF TURNS J. T. Conway University of Agder 9 Lilletuns Vei, Grimstad 4879, Norway Abstract
More informationElectromagnetic Implosion Using a Lens
Sensor and Simulation Notes Note 516 July 2006 Electromagnetic Implosion Using a Lens Carl E. Baum University of New Mexico Department of Electrical and Computer Engineering Albuquerque New Mexico 87131
More informationA GENERAL APPROACH FOR CALCULATING COUPLING IMPEDANCES OF SMALL DISCONTINUITIES
Presented at 16th IEEE Particle Accelerator Conference (PAC 95) and Int l Conference on on High Energy Accelerators, 5/1/1995-5/5/1995, Dallas, TX, USA SLAC-PUB-9963 acc-phys/9504001 A GENERAL APPROACH
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More informationDesign of a Multi-Mode Interference Crossing Structure for Three Periodic Dielectric Waveguides
Progress In Electromagnetics Research Letters, Vol. 75, 47 52, 2018 Design of a Multi-Mode Interference Crossing Structure for Three Periodic Dielectric Waveguides Haibin Chen 1, Zhongjiao He 2,andWeiWang
More informationBehavior of light at photonic crystal interfaces
Behavior of light at photonic crystal interfaces Emanuel Istrate, Alexander A. Green, and Edward H. Sargent Department of Electrical and Computer Engineering, University of Toronto, 10 King s College Road,
More informationPARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS
PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS NAKHLE H. ASMAR University of Missouri PRENTICE HALL, Upper Saddle River, New Jersey 07458 Contents Preface vii A Preview of Applications and
More information9. Dispersion relation of metal nanorods and nanotips
9. Dispersion relation of metal nanorods and nanotips D. E. Chang, A. S. Sørensen, P. R. Hemmer, and M. D. Lukin, Strong coupling of single emitters to surface plasmons, PR B 76,035420 (2007) M. I. Stockman,
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationTitle. Author(s) 牧野, 俊太郎. Issue Date DOI. Doc URL. Type. File Information /doctoral.k
Title A Study on Nonlinear Optical Devices Based on One-Di Author(s) 牧野, 俊太郎 Issue Date 2016-03-24 DOI 10.14943/doctoral.k12194 Doc URL http://hdl.handle.net/2115/61764 Type theses (doctoral) File Information
More informationEigenvalue Analysis of Waveguides and Planar Transmission Lines Loaded with Full Tensor Anisotropic Materials
PIERS ONLINE, VOL. 5, NO. 5, 2009 471 Eigenvalue Analysis of Waveguides and Planar Transmission Lines Loaded with Full Tensor Anisotropic Materials C. S. Lavranos, D. G. Drogoudis, and G. A. Kyriacou Department
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationUNIT I ELECTROSTATIC FIELDS
UNIT I ELECTROSTATIC FIELDS 1) Define electric potential and potential difference. 2) Name few applications of gauss law in electrostatics. 3) State point form of Ohm s Law. 4) State Divergence Theorem.
More informationElectric Field Mapping
Name: Date: PC1143 Physics III Electric Field Mapping 5 Laboratory Worksheet Part A: Parallel Electrodes Distance between parallel electrodes: d = Attach your plots of equipotential lines and electric
More informationDielectric Slab Waveguide
Chapter Dielectric Slab Waveguide We will start off examining the waveguide properties of a slab of dielectric shown in Fig... d n n x z n Figure.: Cross-sectional view of a slab waveguide. { n, x < d/
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
More informationSpherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios. Nicola Tedeschi
Spherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios Nicola Tedeschi Outline Solutions of the vector wave equation Spherical wave approach Spherical reflected and transmitted
More informationNumerical Analysis of Low-order Modes in Thermally Diffused Expanded Core (TEC) Fibers
Proceedings of the 4th WSEAS Int. Conference on Electromagnetics, Wireless and Optical Communications, Venice, Italy, November 2-22, 26 26 Numerical Analysis of Low-order Modes in Thermally Diffused Expanded
More informationDielectric waveguides
Chapter 5 Dielectric waveguides Contents 5.1 Introduction......................................... 5 1 5.2 Modes of Optical Waveguides.............................. 5 5 5.3 Propagation through dielectric
More informationPhysics 505 Fall 2005 Homework Assignment #7 Solutions
Physics 505 Fall 005 Homework Assignment #7 Solutions Textbook problems: Ch. 4: 4.10 Ch. 5: 5.3, 5.6, 5.7 4.10 Two concentric conducting spheres of inner and outer radii a and b, respectively, carry charges
More informationCalculus C (ordinary differential equations)
Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear
More informationMIXED-POTENTIAL GREEN S FUNCTION OF AN AXIALLY SYMMETRIC SHEET MAGNETIC CURRENT ON A CIRCULAR CYLINDRICAL METAL SURFACE
Progress In Electromagnetics Research, PIER 6, 45 64, 6 MIXED-POTENTIAL GREEN S FUNCTION OF AN AXIALLY SYMMETRIC SHEET MAGNETIC CURRENT ON A CIRCULAR CYLINDRICAL METAL SURFACE A. Y. Svezhentsev The A.Ya.
More informationPhotonic/Plasmonic Structures from Metallic Nanoparticles in a Glass Matrix
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Photonic/Plasmonic Structures from Metallic Nanoparticles in a Glass Matrix O.Kiriyenko,1, W.Hergert 1, S.Wackerow 1, M.Beleites 1 and
More informationSPHERICAL NEAR-FIELD ANTENNA MEASUREMENTS
SPHERICAL NEAR-FIELD ANTENNA MEASUREMENTS Edited by J.E.Hansen Peter Peregrinus Ltd. on behalf of the Institution of Electrical Engineers Contents Contributing authors listed Preface v xiii 1 Introduction
More informationGreen s Tensors for the Electromagnetic Field in Cylindrical Coordinates
Chapter 25 Green s Tensors for the Electromagnetic Field in Cylindrical Coordinates The solutions of the vector Helmholtz equation in three dimensions can be expressed by a complete set of vector fields
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationElectromagnetic Waves
Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic
More informationarxiv: v2 [math-ph] 14 Apr 2008
Exact Solution for the Stokes Problem of an Infinite Cylinder in a Fluid with Harmonic Boundary Conditions at Infinity Andreas N. Vollmayr, Jan-Moritz P. Franosch, and J. Leo van Hemmen arxiv:84.23v2 math-ph]
More informationBack to basics : Maxwell equations & propagation equations
The step index planar waveguide Back to basics : Maxwell equations & propagation equations Maxwell equations Propagation medium : Notations : linear Real fields : isotropic Real inductions : non conducting
More informationOn Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material Surfaces
Vol. 114 2008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology On Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationSound radiation from the open end of pipes and ducts in the presence of mean flow
Sound radiation from the open end of pipes and ducts in the presence of mean flow Ray Kirby (1), Wenbo Duan (2) (1) Centre for Audio, Acoustics and Vibration, University of Technology Sydney, Sydney, Australia
More informationThermal lensing in high power ridge waveguide lasers. H. Wenzel, M. Dallmer and G. Erbert
Thermal lensing in high power ridge waveguide lasers H. Wenzel, M. Dallmer and G. Erbert Outline motivation and laser structure experimental results theoretical model simulation results conclusions Motivation
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationBeam Propagation Modelling of Whispering Gallery Microcavities. Mohammad Amin Cheraghi Shirazi B.Sc., University of Tehran, 2011
Beam Propagation Modelling of Whispering Gallery Microcavities by Mohammad Amin Cheraghi Shirazi B.Sc., University of Tehran, 2011 A Thesis Submitted in Partial Fulfillment of the Requirements for the
More informationOptical fibers. Weak guidance approximation in cylindrical waveguides
Optical fibers Weak guidance approximation in cylindrical waveguides Propagation equation with symmetry of revolution r With cylindrical coordinates : r ( r,θ) 1 ψ( r,θ) 1 ψ( r,θ) ψ + + r r + r θ [ ( )
More informationA FLEXIBLE SOLVER OF THE HELMHOLTZ EQUATION FOR LAYERED MEDIA
European Conference on Computational Fluid Dynamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 26 A FLEXIBLE SOLVER OF THE HELMHOLTZ EQUATION FOR LAYERED MEDIA
More information