Computational Electromagnetics in Antenna Analysis and Design

Similar documents
Numerical Heat and Mass Transfer

Appendix B. The Finite Difference Scheme

Formal solvers of the RT equation

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

The Finite Element Method

NUMERICAL DIFFERENTIATION

Interconnect Modeling

DUE: WEDS FEB 21ST 2018

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures

A constant recursive convolution technique for frequency dependent scalar wave equation based FDTD algorithm

Hybrid FDTD/FETD Technique Using Parametric Quadratic Programming for Nonlinear Maxwell s Equations

Lecture 12: Discrete Laplacian

CHAPTER 14 GENERAL PERTURBATION THEORY

Integrals and Invariants of Euler-Lagrange Equations

One-sided finite-difference approximations suitable for use with Richardson extrapolation

B. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea

Consistency & Convergence

2 Finite difference basics

A new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Global Sensitivity. Tuesday 20 th February, 2018

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Modeling acoustic transducer surface waves by Transmission Line Matrix method

Research Article Green s Theorem for Sign Data

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Lecture 21: Numerical methods for pricing American type derivatives

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

PART 8. Partial Differential Equations PDEs

New Method for Solving Poisson Equation. on Irregular Domains

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

The Finite Element Method: A Short Introduction

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

A Hybrid Variational Iteration Method for Blasius Equation

CIE4801 Transportation and spatial modelling Trip distribution

Integrals and Invariants of

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Three-dimensional eddy current analysis by the boundary element method using vector potential

Field computation with finite element method applied for diagnosis eccentricity fault in induction machine

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

Kernel Methods and SVMs Extension

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

4DVAR, according to the name, is a four-dimensional variational method.

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

Optimal Control of Temperature in Fluid Flow

COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC

ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY

Polynomial Regression Models

International Journal of Pure and Applied Sciences and Technology

Homework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Inductance Calculation for Conductors of Arbitrary Shape

Chapter 3 Differentiation and Integration

Difference Equations

Tensor Smooth Length for SPH Modelling of High Speed Impact

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

Monte Carlo Integration Technique for Method of Moments Solution of EFIE in Scattering Problems

Imposition of the essential boundary conditions in transient heat conduction problem based on Isogeometric analysis

A new Approach for Solving Linear Ordinary Differential Equations

Frequency dependence of the permittivity

VQ widely used in coding speech, image, and video

18.1 Introduction and Recap

Generalized Linear Methods

Aging model for a 40 V Nch MOS, based on an innovative approach F. Alagi, R. Stella, E. Viganò

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Waveguides and resonant cavities

2.29 Numerical Fluid Mechanics

Numerical Transient Heat Conduction Experiment

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

Grid Generation around a Cylinder by Complex Potential Functions

DFT with Planewaves pseudopotential accuracy (LDA, PBE) Fast time to solution 1 step in minutes (not hours!!!) to be useful for MD

The Discretization Process

Single-Facility Scheduling over Long Time Horizons by Logic-based Benders Decomposition

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Chapter 4 The Wave Equation

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Numerical Solution of Ordinary Differential Equations

Solution for singularly perturbed problems via cubic spline in tension

NUMERICAL MODELING OF ACTIVE DEVICES CHARACTERIZED BY MEASURED S-PARAMETERS IN FDTD

Module 3: Element Properties Lecture 1: Natural Coordinates

Constitutive Modelling of Superplastic AA-5083

Feb 14: Spatial analysis of data fields

Linear Approximation with Regularization and Moving Least Squares

Relaxation Methods for Iterative Solution to Linear Systems of Equations

Note: Please use the actual date you accessed this material in your citation.

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

16 Reflection and transmission, TE mode

Transcription:

Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures. Classcal approaches may fal f: the materal s not lnear and cannot be lnearzed wthout serously affectng the result the soluton regon s complex (.e. the varous boundares do not concde wth any well descrbed coordnate system). the boundary condtons are tme-dependent the medum s nhomogeneous or ansotropc Whenever a problem wth such complexty arses numercal solutons must be employed. Fortunately there are a large number of very good commercal programs avalable for solvng antenna problems.

Computatonal Electromagnetcs computatonal electromagnetcs rgorous methods Hgh frequency IE DE VM TD FD TD FD feld based current based MoM FEM FDTD TLM GO/GTD PO/PTD

Comparson of Methods Doman Generalty Accuracy Memory N= number of elements MoM Frequency Homogeneous or dscretely homogeneous regons FDTD Tme (all frequences n one run) Very general. nhomogeneous dspersve ansotropc FEM Frequency Very general. nhomogeneous dspersve ansotropc Hgh Frequency Methods Frequency Only good for structures much larger than the wavelength Antenna Types Very accurate @(N ) All but harder for large reflector antennas Moderately accurate @(N) All but harder for large reflector antennas Very accurate @(N log N) All but harder for large reflector antennas Only accurate for large structures @(N) Only good for large antennas (mostly used for reflectors)

Goal of all of these methods Approxmate these B D E = M H = t t D = ρ B = ρ Wth ths m J [ A ] [ x] = [ b]

Fnte Dfference Tme Doman FDTD

Reason for nterest n FDTD In the tme doman Maxwell s equatons gve rse to PDEs nvolvng tme and spatal dervatves. Some good reasons for dealng wth PDE s are: Complex-value materals easly accommodated. Computer resources are adequate. PDE solutons are robust. Tme doman PDE methods usually have no matrces. Geometres to be solved can be more vared.

Some bg advantages Broadband response wth a sngle exctaton. 3D models easly. Memory requrement scales lnearly wth problem sze Frequency dependent materals accommodated. Most parameters can be generated e.g. Scattered felds antenna patters RCS S-parameters etc..

How does t wor? Based on the two Maxwell curl equatons n dervatve form. These are lnearzed by central fnte dfferencng. We only consder nearest neghbor nteractons because all the felds are advanced temporally n dscrete tme steps over spatal cells. e we sample n space & tme embeddng of an antenna n a FDTD space lattce (note that the whole volume s meshed!)

Dscretze Obects n Space usng Cartesan Grd D Dscretzaton 3D Dscretzaton z x D Dscretzaton Z z = 0 Ex () zt z = Z

Defne Locatons of Feld Components: FDTD Cell called Yee Cell Fnte-Dfference Space s dvded nto small cells One Cell: (dx)(dy)(dz) E and H components are dstrbuted n space around the Yee cell (note: feld components are not collocated) FDTD: Yee K. S.: Numercal soluton of ntal boundary value problems nvolvng Maxwell's equatons n sotropc meda. IEEE Transactons on Antennas Propagaton Vol. AP-4 pp. 30-307 966.

3D formulaton y x z x y x z y y x z z E H E H t z y H E E H t x z E E H H t y x ρ µ ρ µ ρ µ = = = y x z x y x z y y x z z H E H E t y z E H H E t z x H H E E t x y σ ε σ ε σ ε = = = ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) = = z E E x E E H H y E E z E E H H n x n x n z n z r n x n x n z n z n y n y r n x n x 0 0 µ µ µ µ ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) = = 0 0 E E z H H y H H t E E x E E y E E t H H n x n x n y n y n z n z r n x n x n y n y n x n x r n z n z σ ε ε µ µ Convert equatons le these To ones le these

QucWave (http://www.qwed.eu/) Commercal FDTD pacage wth CAD nterface Uses conformal FDTD mesh Many specal features for antenna problems Wrtten and supported by QWED Poland Runs under 3/x64 bt Wndows platforms and Lnux

QucWave (http://www.qwed.eu/) Examples:

QucWave (http://www.qwed.eu/) Examples:

XFdtd (http://www.remcom.com /) Commercal FDTD pacage wth CAD nterface Probably the most popular FDTD pacage for antenna problems Many specal features for antenna problems ncludng full human body mesh Wrtten and supported by Remcom USA Runs under 3/x64 bt Wndows Mac OS X and Lnux

XFdtd (http://www.remcom.edu /)

Lumercal (http://www.lumercal.com /) Commercal FDTD pacage wth CAD nterface Popular FDTD pacage for the optcs fols (.e. ntegrated optcs lght scatterng plasmoncs) Large lbrary of materals at optcal wavelengths Has very nce scrptng language for defnng large complcated problems. Has nce bult n optmzaton Runs under 3/x64 bt Wndows Mac OS X and Lnux

Lumercal (http://www.lumercal.com /)

Fnte Element Method

Fnte Element Method

Varatonal Approach In solvng problems arsng n physcs and engneerng t s often possble to replace the problem of ntegratng a dfferental equaton by the equvalent problem of seeng a functon that gves a mnmum value of some ntegral. Problems of ths type are called varatonal problems. The methods that allow us to reduce the problem of ntegratng a dfferental equaton to the equvalent varatonal problem are usually called varatonal methods.

Varatonal Approach Name of equatons PDE Varatonal prncple Homogeneous wave equaton wth sources Homogeneous wave equaton wthout sources Dffuson equaton Posson s equaton Homogenous Laplace s equaton Φ Φ = Φ Φ Φ = Φ t Φ = Φ = g 0 = g 0 0 [ I( Φ) = Φ Φ gφ]dv v I( Φ) = v [ Φ Φ ]dv Φ I( Φ) = Φ Φ dvdt t t v [ I( Φ) = Φ gφ]dv v I( Φ) = v [ Φ ]dv

Fnte Element Method The fnte element method (FEM) has ts orgn n the feld of structural analyss. However snce then the method has been employed n nearly all areas of computatonal physcs and engneerng. The FEM method whle more dffcult to program than ether the fnte dfference (FD) or method of moments (MOM) s a more powerful and versatle numercal technque for handlng problems nvolvng complex geometres and nhomogeneous meda.

Basc concept Although the behavour may be complex when vewed over a large regon a smple approxmaton may suffce over a small subregon. The regon s dvded up nto fnte elements. (usually trangles or squares but can be more complcated) Regardless of the shape the feld s approxmated by a dfferent expresson over each element mantanng contnuty at adonng elements.

Soluton Strategy: Varatonal Approach The equatons to be solved are usually stated not n terms of feld the varables but n terms of an ntegral-type functonal such as energy. The functonal s chosen such that the feld soluton maes the functonal statonary The total functonal s the sum of the ntegral over each element

Fnte Element Method The fnte element method (FEM) nvolves bascally four steps: () Dscretze the soluton regon nto a fnte number of subregons or elements () Derve the governng equatons for each element based on ether a varatonal approach or Galern s method (3) Assemble all the elements together n the soluton space. (4) Solve the resultng system of equatons

HFSS (http://www.ansoft.com/) Commercal FEM pacage wth CAD nterface Uses adaptve meshng Probably the most popular commercal pacage for antenna applcatons. Wrtten and supported by Ansoft USA Runs under 3/x64 bt Wndows platforms Redhat Lnux Solars (Sun worstatons). Has ntegrated hybrd fnte element / boundary ntegral methods (MoM) Knd of expensve! Optonal optmzaton pacage (optmetrcs)

HFSS (http://www.ansoft.com/)

CST Mcrowave Studo (http://www.cst.com/) Commercal FEM MoM and TLM pacage wth CAD nterface Mature and easy to use nterface Popular program for mcrowave crcut applcatons but also very useful for antennas. Wrtten and supported by CST Internatonal Runs under 3/x64 bt Wndows platforms and Redhat Lnux Knd of expensve! Has a specfc antenna desgn opton called Magus

CST Mcrowave Studo (http://www.cst.com/)

Comsol Multphyscs(http://www.comsol.com/) Commercal FEM pacage wth CAD nterface Is best nown for ts ablty to solve multphyscs problems Becomng a very popular program. Lns ncely wth Matlab Easy to learn nterface Runs under 3/x64 bt Wndows platforms Moderately prces

Comsol Multphyscs(http://www.comsol.com/)

References:. Clemson ste lsts all the free EM modelng tools http://www.clemson.edu/ces/cvel/modelng/emag/free-codes.html. Clemson ste lsts all the commercal EM modelng tools http://www.clemson.edu/ces/cvel/modelng/emag/csoft.html