A simple introduction to the method of lines
|
|
- Gwenda Goodwin
- 6 years ago
- Views:
Transcription
1 M. N. O. Sadiku1 and C. N. Obiozor2 1Department of Electrical and Computer Engineering, Temple University, Philadelphia, U.S.A. 2Department of Electrical Engineering, University of North Florida, Jacksonville, U.S.A. Abstract The method of lines (MOL), a semianalytical procedure, is well known to experts in computational techniques in electromagnetics. The range of applications of the method has increased dramatically in the past few years; nevertheless, there is no introductory paper to initiate a beginner to the method. This paper illustrates the application of the MOL to solve Laplace s equation in rectangular and cylindrical coordinates. Two numerical examples are used to verify the procedure. The results obtained compare well with analytical solutions. Keywords numerical analysis; electromagnetics The method of lines (MOL) is a well established numerical technique (or rather a semianalytical method) for the analysis of transmission lines, waveguide structures, and scattering problems. The method, originally developed by mathematicians and used for boundary value problems in physics [see, e.g., Refs. (1 3)], was introduced into the electromagnetic (EM) community and further developed by Pregla et al.4 7 and other researchers. Although the method of lines has become one of the standard tools for solving practical, complex electromagnetic field problems, the available literature has failed to supply the reader with implementation ready equations and introductory material to a new beginner. Hence, there is yet to be an introductory monograph to initiate a beginner to the method. The only book3 on MOL is not geared toward EM community and the formulation of MOL is therefore different from the modern approach developed by Pregla et al. The only introductory chapter4 for the EM community is advanced because it is is geared towards experts. In this paper, a simple introduction that breaks MOL down into implementation ready equations is presented and illustrated with examples. General background The method of lines is regarded as a special finite difference method but more effective with respect to accuracy and computational time than the regular finite difference method. It basically involves discretising a given differential equation in one or two dimensions while using analytical solution in the remaining direction. MOL has the merits of both the finite difference method and analytical method; it does not yield spurious modes nor have the problem of relative convergence. Besides, the method of lines has the following properties that justify its use: (a) Computational efficiency: the semianalytical character of the formulation leads to a simple and compact algorithm, which yields accurate results
2 283 with less computational effort than other techniques. ( b) Numerical stability: by separating discretisation of space and time, it is easy to establish stability and convergence for a wide range of problems. (c) Reduced programming effort: by making use of the state-of-the-art well documented and reliable ordinary differential equations (ODE) solvers, programming effort can be substantially reduced. (d) Reduced computational time: since only a small amount of discretisation lines are necessary in the computation, there is no need to solve a large system of equations; hence computing time is small. To apply MOL usually involves the following five basic steps: 1 Partitioning the solution region into layers. 2 Discretisation of the differential equation in one coordinate direction. 3 Transformation to obtain decoupled ordinary differential equations. 4 Inverse transformation and introduction of the boundary conditions. 5 Solution of the equations. Although the method of lines is commonly used in the EM community for solving hyperbolic (wave equation), it can be used to solved parabolic and elliptic equation.1,8 11 In this paper, we consider the simple case of applying MOL to solve Laplace s equation involving two-dimensional rectangular and cylindrical regions. Laplace equation in rectectangular coordinates Laplace s equation in the Cartesian system is 2V x2 + 2V =0 (1) y2 Consider a two-dimensional solution shown in Fig. 1. The first step is discretisation of the x-variable. The region is divided into strips by N dividing straight lines ( hence the name method of lines) parallel to the y-axis. Since we are discretising along x, we replace the second derivative with respect to x with its finite difference equivalent. We apply the three-point central difference scheme, 2V i x2 = V i+1 2V i +V i 1 (2) h2 where h is the spacing between discretised line, i.e. h=dx= a (3) N+1 Replacing the derivative with respect to x by its finite difference equivalent, equation (1) becomes 2V i y2 + 1 h2 [V i+1 (y) 2V i (y)+v (y)]=0 (4) i 1
3 284 M. N. O. Sadiku and C. N. Obiozor Fig. 1 Illustration of discretisation in the x-direction. Thus the potential V in equation (1) can be replaced by a vector of size N, namely [V ]=[V 1,V 2,..., V N ]t (5a) where V i (y)=v(x i,y), i=1, 2,..., N (5b) and x i =idx. Substituting equations (4) and (5) into (1) yields 2[V (y)] y2 1 [P][V (y)]=0 (6) h2 where [P] is an N N tridiagonal matrix representing the discretised form of the second derivative with respect to x. t p 1 0, 0 l u N , 0 N N N [P]= P P P (7) N N N 0, N v 0, 0 1 p r w All the elements of matrix [P] are zeros except the tridiagonal terms; the elements of the first and the last row of [P] depends on the boundary conditions at x=0 and x=a. p =2 for the Dirichlet boundary condition and p =1 for l l the Neumann boundary condition. The same is true of p. r The next step is to solve the resulting equations analytically along the
4 285 y coordinate. To solve equation (6) analytically, we need to obtain a system of uncoupled ordinary differential equations from the coupled equations (6). To achieve this, we define the transformed potential [ V9 ] by letting [V]=[T][V9 ] (8) and requiring that [ T]t[P][ T] =[l2] (9) where [l2] is a diagonal matrix and [T]t is the transpose of [T]. [l2] and [ T] are eigenvalue and eigenvector matrices belonging to [P]. The transformation matrix [T] and the eigenvalue matrix [l2] depend on the boundary conditions and are given in Table 1 for various combinations of boundaries. It should be noted that the eigenvector matrix [ T] has the following properties: [T] 1=[T]t [T][T]t=[T]t[T]=[I] (10) where [I] is an identity matrix. Substituting equation (8) into equation (6) gives 2[T][V9 ] 1 [P][T][V9 ]=0 y2 h2 Multiplying through by [ T] 1=[ T]t yields A 2 y2 1 h2 [l2] [V9 ]=0 (11) B TABLE 1 Elements of transformation matrix [ T] and eigenvaluesa Left Right boundary boundary emitters susceptors T ij l i Dirichlet Dirichlet ip 2 sin S 2 ijp sin N+1 N+1,[T DD ] 2(N+1) Dirichlet Neumann 2 sin S 2 i( j 0.5)p (i 0.5)p sin N+0.5 N+0.5,[T DN ] 2N+1 Neumann Dirichlet S 2 (i 0.5)( j 0.5)p cos,[t ] N+0.5 N+0.5 ND 2 sin (i 0.5)p 2N+1 Neumann Neumann 2 sin S 2 (i 0.5)( j 1)p (i 1)p cos,j>1, [T ] N N NN 2N 1 N,j=1 a i, j=1, 2,..., N and subscripts D and N denote Dirichlet and Neumann conditions, respectively.
5 286 M. N. O. Sadiku and C. N. Obiozor This is an ordinary differential equation with solution V9 i =A i cosh a i y+b i sinh a i y (12) where a i =l i /h. Thus, Laplace s equation is solved numerically using a finite difference scheme in the x-direction and analytically in the y-direction. However, we have only demonstrated three out of the five basic steps for applying MOL. There remain two more steps to complete the solution: imposing the boundary conditions and solving the resulting equations. Imposing the boundary conditions is problem dependent and will be illustrated in Example 1. The resulting equations can be solved using the existing packages for solving ODE or developing our own codes in Fortran, MATLAB, C, or any programming language. We will take the latter approach in Example 1: Example 1 For the rectangular region in Fig. 1, let V(0, y)=v (a, y) =V (x, 0)=0, V (x, b)=100 and a=b=1. Find the potential at (0.25, 0.75), (0.5, 0.5), (0.75, 0.25). Solution In this case, we have Dirichlet boundaries at x=0 and x=1, which are already indirectly taken care of in the solution in (12). Hence, from Table 1, ip l =2 sin i 2(N+1) and T ij =S 2 ijp sin N+1 N+1 Let N=15 so that h=dx=1/16 and x=0.25, 0.5, 0.75 will correspond to i= 4, 8, 12 respectively. By combining equations (8) and ( 12), we obtain the required solution. To get constant A and B, we apply boundary conditions at y=0 and y=b to V i i and perform inverse transformation. Imposing V (x, y=0)=0 onto the combination of equations (8) and (12), we obtain CV T, T N V T T, T N e e, e V T ND=[0]=CT NNDCA 1 A ND 2 (15) e T, T A N1 N2 which implies that (13) (14) [A]=0 or A i =0 (16)
6 287 Imposing V (x, y=b) =100 yields C e 100D=[T]CB If we let 1 sinh a 1 b B sinh a b 2 2 bd (17) e B sinh a N N sinh a b 1 1 B sinh a b [C]=CB 2 2 (18) e e B sinh a bd=[t] 1C100 N N 100D then B i =C i /sinh a i b (19) With A i and B i found in equations (16) and (19), the potential V (x, y) is determined as V (y)= N T B sinh(a y) (20) i ij j j j=1 By applying equation (13) to equation (20), the MATLAB code in Fig. 2 was developed to obtain V ( 0.25, 0.75) =43.1, V (0.5, 0.5)=24.96, V (0.75, 0.25) =6.798 The result compares well with the exact solution:12 V ( 0.25, 0.75) =43.2, V (0.5, 0.5)=25.0, V (0.75, 0.25)=6.797 Notice that it is not necessary to invert the transformation matrix [T] in view of equation (10). Laplace equation in cylindrical coordinates Although MOL is not applicable to problems with complex geometry, the method can be used to analyse homogeneous and inhomogeneous cylindrical transmission structures13 17 and circular and elliptic waveguides.18 Here, we illustrate with the use of the MOL to solve Laplace s equation in cylindrical coordinates. The principal steps in applying MOL in cylindrical coordinates are the same as in Cartesian coordinates. We apply discretisation procedure in the angular direction. The resulting coupled ordinary differential equations are decoupled by matrix transformation and solved analytically.
7 288 M. N. O. Sadiku and C. N. Obiozor Fig. 2 MAT L AB code for Example 1. Assume that we are interested in finding the potential distribution in cylindrical transmission line with uniform but arbitrary cross-section. We assume that the inner conductor is grounded while the outer conductor is maintained at constant potential V o, as shown in Fig. 3. In cylindrical coordinates (r, w), Laplace s equation can be expressed as r2 2V V +r r2 r + 2V =0 (21) w2 subject to V(r)=0, rµc 1 (22a) V(r)=V o, rµc 2 (22b)
8 289 Fig. 3 Discretisation along w-direction. We discretise in the w-direction by using N radial lines, as shown in Fig. 3, such that V i (r)=v(r, w i ), i=1, 2,..., N (23) where w i =ih= 2pi N, h=dw= 2p N (24) and h is the angular spacing between the lines. We have subdivided the solution region into N subregions with boundaries at C 1 and C 2. In each subregion, V(r, w) is approximated by V i =V(r, w i ), with w i being constant. Applying the three-point central finite difference scheme yields 2[V] = [P] [V] (25) w2 h2 where [V]=[V 1,V 2,..., V N ]t (26)
9 290 M. N. O. Sadiku and C. N. Obiozor and t , 0 0 1u N , N N , N [P]= N N N e e e e, e e e N N , 1 2 1N v , w (27) Notice that [P] contains an element 1 in the lower left- and right-hand corners due to its angular periodicity. Also, notice that [P] is a quasi-threeband symmetric matrix which is independent of the arbitrariness of the crosssection as a result of the discretisation over a finite interval [0, 2p]. Introducing equation (25) into (21) leads to the following set of coupled differential equations r2 2[V] r2 +r [ V ] r [P] [V]=0 (28) h2 To decouple equation (28), we must diagonalise [P] by an orthogonal matrix [T] such that with [l2]=[ T]t[P][ T] (29) [T]t=[T]=[T] 1 (30) where [l2] is a diagonal matrix of the eigenvalues l2 n of [P]. The diagonalisation is achieved using T = cos a ij +sin a ij, l2=2(1 cos a ) (31) ij N n n where a ij =hωiω j, a n =hωn, i, j, n=1, 2,..., N (32) If we introduce the transformed potential U that satisfies [U]=[T][V] (33) Equation (28) becomes r2 2[U] r2 where +r [U] [m2][u]=0 (34) r [U]=[U 1,U 2,..., U N ]t (35)
10 291 is a vector containing the transformed potential function and m = l n n h = 2 h sin(a /2) (36) n Equation (34) is the Euler-type and has the analytical solution:19 U = n GA n +B n ln r, m n =0 A rm n n +B r m n n, m 0 n (37) This is applied to each subregion. By taking the inverse transform using equation (33), we obtain the potential V i (r) as V (r)= N T U (38) i ij j j=1 where T are the elements of matrix [T]. ij We now impose the boundary conditions in equation (22), which can be rewritten as V(r=r i )=0, r i µc 1 (39a) V(r=R i )=V o, R i µc 2 (39b) Applying these to equations (37) and (38), T ij [A j +B j ln r i ] mj =0 = N j=1 T ij [A j +B j ln R i ] mj =0 + N j=1 T ij [A j rm j i T ij [A j Rm j i +B r m j i j ] mj 0 =0, i=1, 2,..., N (40a) +B R m j i j ] mj 0 =V, i=1, 2,..., N o (40b) Equation (40) is solved to determine the unknown coefficients A i and B i. The potential distribution is finally obtained from equations (26) and (38). Example 2 Consider a coaxial cable with inner radius a and outer radius b. Let b=2a= 2 cm and V o =100 V. This simple example is selected to be able to compare the MOL solution with the exact solution. Solution From equation (36), it is evident that m n =0 only when n=n. Hence, we may write U as U = n GA n rm n +B n r m n, A +B ln r, n n n=1, 2,..., N 1 n=n Equation (40) can be written as N 1 T [A am ij j j +B a m i j j ]+T [A +B ln a]=0, i in N N j=1 i=1, 2,..., N (42a) (41)
11 292 M. N. O. Sadiku and C. N. Obiozor for r=a, and N 1 T [A bm ij j j +B b m i j j ]+T [A +B ln b]=v, i=1, 2,..., N (42b) i in N N o j=1 for r=b. These 2N equations will enable us find the 2N unknown coefficients A and B. They can cast into a matrix form as i i ta 1 u t 0 u tt A am , T A T B a m 1N N , B ln a N u NA 2 N N 0 N N e e N N e N N e N NT A am N N1 1 1, T A T B a m NN N N1 1 1, B ln an NA N N 0 N N N N N N = N N NT A bm , T A T B b m 1N N , B ln b N N NB 1 N N100N N e e N NB 2 N N100N vt A bm N1 1 1, T A T B b m NN N N1 1 1, B ln bw N e N N e N N vb w v100w N (43) This can be written as [D][C]=[F] (44) from which we obtain [C] =[D] 1[F] (45) where C j corresponds to A j when j=1, 2,..., N and C j corresponds to B j when j=n+1,..., 2N. Once A j and B j are known, we substitute them into equation (41) to find U j. We finally apply equation (38) to find V. The exact analytical solution of the problem is:12 r ln a V(r)=V o (46) ln b a For a<r<b, we obtain V for both exact and MOL solutions using the MATLAB codes in Fig. 4. The results of two solutions are shown in Fig. 5. The two solutions agree perfectly. Conclusions The method of lines is a differential-difference approach of solving partial differential equations. It involves a judicious combination of analysis and computation. Given a partial differential equation, all but one of the independent
12 293
13 294 M. N. O. Sadiku and C. N. Obiozor Fig. 4 MAT L AB code for Example 2. variables are discretised to obtain a system of ordinary differential equations (ODEs), which are then solved using ODE solvers. In this paper, a simple introduction on MOL with illustrative examples has been presented. It is hoped that with this introduction, one can begin to learn more about MOL and apply the semianalytical technique to solving EM problems.
14 295 Fig. 5 Comparison of exact and method of lines solutions. References 1 A. Zafarullah, Application of the method of lines to parabolic partial differential equations with error estimates, J. ACM, 17(2), (1970), M. B. Carver and H. W. Hinds, The method of lines and the advection equation, Simulation, 31(2) (1978), W. E. Schiesser, T he Numerical Method of L ines (Academic Press, San Diego, 1991). 4 R. Pregla and W. Pascher, The method of lines, in T. Itoh (ed.), Numerical T echniques for Microwave and Millimeter-wave Passive Structures (John Wiley, New York, 1989), pp R. Pregla, Analysis of planar microwave structures on magnetized ferrite substrate, Arch. Elek. Ubertragung, 40 (1986), R. Pregla, About the nature of the method of lines, Arch. Elek. Ubertragung, 41(6) (1987), U. Schulz and R. Pregla, A new technique for the analysis of the dispersion characteristics of planar waveguides, Arch. Elek. Ubertragung, 34 (1980), pp Also in R. Sorrentino (ed.), Numerical Methods for Passive Microwave and Millimeter Wave Structure (IEEE Press, New York, 1989), pp U. Schulz and R. Pregla, A new technique for the analysis of the dispersion characteristics of planar waveguides and its application to microstrips with tuning septums, Radio Sci., 16 (1981), D. J. Jones et al., On the numerical solution of elliptic partial differential equations by the method of lines, J. Comput. Phys., 9 (1972), B. H. Edwards, The numerical solution of elliptic partial differential equation by the method of lines, Revista Colombiana de Matematicas, 19 (1985),
15 296 M. N. O. Sadiku and C. N. Obiozor 11 A. G. Keen et al., Quasi-static analysis of electrooptic modulators by the method of lines, J. L ightwave T echnol., 8(1) (1990), M. N. O. Sadiku, Elements of Electromagnetics, 2nd edn. (Oxford University Press, New York, 1994), pp. 243, J. G. Ma and Z. Chen, Application of the method of lines to the Laplace equation, Micro. Opt. T ech. L ett., 14(6) (1997), S. Xiao et al., Analysis of cylindrical transmission lines with the method of lines, IEEE T rans. Micro. T heo. T ech., 44(7) (1996), K. Gu and Y. Wang, Analysis of dispersion characteristics of cylindrical microstrip line with method of lines, Elect. L ett., 26(11) (1990), R. Pregla, General formulas for the method of lines in cylindrical coordinates, IEEE T rans. Micro. T heo. T ech., 43(7) (1995), M. Thornburn et al., Application of the method of lines to cylindrical inhomogeneous propagation structures, Elect. L ett., 26(3) (1990), K. Wu and R. Vahldieck, The method of lines applied to planar transmission lines in circular and elliptical waveguides, IEEE T rans. Micro. T heo. T ech., 37(12) (1989), M. N. O. Sadiku, Numerical T echniques in Electromagnetics (CRC Press, Boca Raton, 1992), pp
The calculation of the field distribution is very accurate because of the relation of the MoL to the Discrete Fourier Transformation (DFT) [5].
CHAPTER 1 THE METHOD OF LINES 1.1 INTRODUCTION The Method of Lines (MoL) is one of the most efficient numerical methods for solving the partial differential equations with which we describe physical phenomena.
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 informationFast Direct Solver for Poisson Equation in a 2D Elliptical Domain
Fast Direct Solver for Poisson Equation in a 2D Elliptical Domain Ming-Chih Lai Department of Applied Mathematics National Chiao Tung University 1001, Ta Hsueh Road, Hsinchu 30050 Taiwan Received 14 October
More informationDESIGN OF ARBITRARY SHAPED PLANAR RES- ONATORS WITH FINE DETAILS USING MODIFIED SPACE SPECTRAL DOMAIN APPROACH
Progress In Electromagnetics Research B, Vol. 48, 249 269, 2013 DESIGN OF ARBITRARY SHAPED PLANAR RES- ONATORS WITH FINE DETAILS USING MODIFIED SPACE SPECTRAL DOMAIN APPROACH Essam A. Hashish 1 and Hossam
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationELEC 621 NUMERICAL TECHNIQUES IN ELECTROMAGNETICS
- 1 - ELEC 621 NUMERICAL TECHNIQUES IN ELECTROMAGNETICS Fall 2005 Wolfgang J.R. Hoefer, Professor Department of Electrical and Computer Engineering University of Victoria, Victoria, B.C. Tel. (250) 721-6030
More informationSpectral Domain Analysis of Open Planar Transmission Lines
Mikrotalasna revija Novembar 4. Spectral Domain Analysis of Open Planar Transmission Lines Ján Zehentner, Jan Mrkvica, Jan Macháč Abstract The paper presents a new code calculating the basic characteristics
More information4.4 Microstrip dipole
4.4 Microstrip dipole Basic theory Microstrip antennas are frequently used in today's wireless communication systems. Thanks to their low profile, they can be mounted to the walls of buildings, to the
More informationModelling of Ridge Waveguide Bends for Sensor Applications
Modelling of Ridge Waveguide Bends for Sensor Applications Wilfrid Pascher FernUniversität, Hagen, Germany n 1 n2 n 3 R Modelling of Ridge Waveguide Bends for Sensor Applications Wilfrid Pascher FernUniversität,
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 informationFast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media
Excerpt from te Proceedings of te OMSOL onference 2009 Boston Fast omputation of apacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media S. M. Musa
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
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 informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationEIGENVALUE ANALYSIS OF SPHERICAL RESONANT CAVITY USING RADIAL BASIS FUNCTIONS
Progress In Electromagnetics Research Letters, Vol. 24, 69 76, 2011 EIGENVALUE ANALYSIS OF SPHERICAL RESONANT CAVITY USING RADIAL BASIS FUNCTIONS S. J. Lai 1, *, B. Z. Wang 1, and Y. Duan 2 1 Institute
More informationDiagonalization by a unitary similarity transformation
Physics 116A Winter 2011 Diagonalization by a unitary similarity transformation In these notes, we will always assume that the vector space V is a complex n-dimensional space 1 Introduction A semi-simple
More informationFinite Element Method Analysis of Symmetrical Coupled Microstrip Lines
International Journal of Computing and Digital Systems ISSN (20-142X) Int. J. Com. Dig. Sys. 3, No.3 (Sep-2014) Finite Element Method Analysis of Symmetrical Coupled Microstrip Lines Sarhan M. Musa and
More informationEfficient Calculation of Surface Impedance for Rectangular Conductors. Emre Tuncer and Dean P. Neikirk
Efficient Calculation of Surface Impedance for Rectangular Conductors Emre Tuncer and Dean P. Neikirk Electrical and Computer Engineering Department University of Texas at Austin Austin, Texas 78712 Abstract
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationPhysics 6303 Lecture 8 September 25, 2017
Physics 6303 Lecture 8 September 25, 2017 LAST TIME: Finished tensors, vectors, and matrices At the beginning of the course, I wrote several partial differential equations (PDEs) that are used in many
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationPhy207 Exam III (Form1) Professor Zuo Spring Semester 15
Phy207 Exam III (Form1) Professor Zuo Spring Semester 15 On my honor, I have neither received nor given aid on this examination Signature: Name: ID number: #1 14 #15/16 Total Enter your name and Form 1
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationUpon successful completion of MATH 220, the student will be able to:
MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationGENERALIZED MATRIX MULTIPLICATION AND ITS SOME APPLICATION. 1. Introduction
FACTA UNIVERSITATIS (NIŠ) Ser Math Inform Vol, No 5 (7), 789 798 https://doiorg/9/fumi75789k GENERALIZED MATRIX MULTIPLICATION AND ITS SOME APPLICATION Osman Keçilioğlu and Halit Gündoğan Abstract In this
More informationMatrix Arithmetic. a 11 a. A + B = + a m1 a mn. + b. a 11 + b 11 a 1n + b 1n = a m1. b m1 b mn. and scalar multiplication for matrices via.
Matrix Arithmetic There is an arithmetic for matrices that can be viewed as extending the arithmetic we have developed for vectors to the more general setting of rectangular arrays: if A and B are m n
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationAnalysis of Single Mode Step Index Fibres using Finite Element Method. * 1 Courage Mudzingwa, 2 Action Nechibvute,
Analysis of Single Mode Step Index Fibres using Finite Element Method. * 1 Courage Mudzingwa, 2 Action Nechibvute, 1,2 Physics Department, Midlands State University, P/Bag 9055, Gweru, Zimbabwe Abstract
More informationReview of Linear Algebra
Review of Linear Algebra Dr Gerhard Roth COMP 40A Winter 05 Version Linear algebra Is an important area of mathematics It is the basis of computer vision Is very widely taught, and there are many resources
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
More informationCONSIDER a simply connected magnetic body of permeability
IEEE TRANSACTIONS ON MAGNETICS, VOL. 50, NO. 1, JANUARY 2014 7000306 Scalar Potential Formulations for Magnetic Fields Produced by Arbitrary Electric Current Distributions in the Presence of Ferromagnetic
More informationModule 6 : Wave Guides. Lecture 40 : Introduction of Parallel Waveguide. Objectives. In this course you will learn the following
Objectives In this course you will learn the following Introduction of Parallel Plane Waveguide. Introduction of Parallel Plane Waveguide Wave Guide is a structure which can guide Electro Magnetic Energy.
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationAnalysis and Design of the CRLH SICL Unit Cell using Effective Parameters
This article has been accepted and published on J-STAGE in advance of copyediting. Content is final as presented. IEICE Electronics Express, Vol.* No.*,*-* Analysis and Design of the CRH SIC Unit Cell
More information. (a) Express [ ] as a non-trivial linear combination of u = [ ], v = [ ] and w =[ ], if possible. Otherwise, give your comments. (b) Express +8x+9x a
TE Linear Algebra and Numerical Methods Tutorial Set : Two Hours. (a) Show that the product AA T is a symmetric matrix. (b) Show that any square matrix A can be written as the sum of a symmetric matrix
More informationNotation, Matrices, and Matrix Mathematics
Geographic Information Analysis, Second Edition. David O Sullivan and David J. Unwin. 010 John Wiley & Sons, Inc. Published 010 by John Wiley & Sons, Inc. Appendix A Notation, Matrices, and Matrix Mathematics
More informationLecture 2 INF-MAT : A boundary value problem and an eigenvalue problem; Block Multiplication; Tridiagonal Systems
Lecture 2 INF-MAT 4350 2008: A boundary value problem and an eigenvalue problem; Block Multiplication; Tridiagonal Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University
More informationFall Inverse of a matrix. Institute: UC San Diego. Authors: Alexander Knop
Fall 2017 Inverse of a matrix Authors: Alexander Knop Institute: UC San Diego Row-Column Rule If the product AB is defined, then the entry in row i and column j of AB is the sum of the products of corresponding
More informationFourier Decomposition Analysis af Anisotropic Inhomogeneous Dielectric Waveguide Structures
University of Pennsylvania ScholarlyCommons Departmental Papers (ESE) Department of Electrical & Systems Engineering August 2007 Fourier Decomposition Analysis af Anisotropic Inhomogeneous Dielectric Waveguide
More informationMathematics. EC / EE / IN / ME / CE. for
Mathematics for EC / EE / IN / ME / CE By www.thegateacademy.com Syllabus Syllabus for Mathematics Linear Algebra: Matrix Algebra, Systems of Linear Equations, Eigenvalues and Eigenvectors. Probability
More informationQuasi-static analysis of asymmetric cylindrical coplanar waveguides with nite-extent ground
INT. J. ELECTRONICS, 2000, VOL. 87, NO. 2, 185± 198 Quasi-static analysis of asymmetric cylindrical coplanar waveguides with nite-extent ground NIHAD DIB*{ and ASEM AL-ZOUBI{ Quasi-static characteristics
More informationLinear Equations in Linear Algebra
1 Linear Equations in Linear Algebra 1.7 LINEAR INDEPENDENCE LINEAR INDEPENDENCE Definition: An indexed set of vectors {v 1,, v p } in n is said to be linearly independent if the vector equation x x x
More information[3] (b) Find a reduced row-echelon matrix row-equivalent to ,1 2 2
MATH Key for sample nal exam, August 998 []. (a) Dene the term \reduced row-echelon matrix". A matrix is reduced row-echelon if the following conditions are satised. every zero row lies below every nonzero
More informationMatrices and Linear Algebra
Contents Quantitative methods for Economics and Business University of Ferrara Academic year 2017-2018 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2
More informationMATH 2030: MATRICES ,, a m1 a m2 a mn If the columns of A are the vectors a 1, a 2,...,a n ; A is represented as A 1. .
MATH 030: MATRICES Matrix Operations We have seen how matrices and the operations on them originated from our study of linear equations In this chapter we study matrices explicitely Definition 01 A matrix
More information5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers.
Linear Algebra - Test File - Spring Test # For problems - consider the following system of equations. x + y - z = x + y + 4z = x + y + 6z =.) Solve the system without using your calculator..) Find the
More informationProgress In Electromagnetics Research, Vol. 110, , 2010
Progress In Electromagnetics Research, Vol. 110, 371 382, 2010 CLASS OF ELECTROMAGNETIC SQ-MEDIA I. V. Lindell Department of Radio Science and Engineering School of Science and Technology, Aalto University
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationThe dependence of the cross-sectional shape on the hydraulic resistance of microchannels
3-weeks course report, s0973 The dependence of the cross-sectional shape on the hydraulic resistance of microchannels Hatim Azzouz a Supervisor: Niels Asger Mortensen and Henrik Bruus MIC Department of
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 informationTowards Classification of Separable Pauli Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 2, 768 773 Towards Classification of Separable Pauli Equations Alexander ZHALIJ Institute of Mathematics of NAS of Ukraine,
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 informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
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 informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 4: Changing Shapes: Linear Maps in 2D Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationIntroduction to Polarization
Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters
More informationT [2] have been developed over the last few years in the
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 43, NO. 3, MARCH 1995 647 Efficient Solution of the Differential Form of Maxwell s Equations in Rectangular Regions Luis Emilio Garcia-Castillo,
More informationComposites with tuned effective magnetic permeability
JOURNAL OF APPLIED PHYSICS 102, 014901 2007 Composites with tuned effective magnetic permeability Alireza V. Amirkhizi and Sia Nemat-Nasser a Center of Excellence for Advanced Materials, Department of
More informationBare minimum on matrix algebra. Psychology 588: Covariance structure and factor models
Bare minimum on matrix algebra Psychology 588: Covariance structure and factor models Matrix multiplication 2 Consider three notations for linear combinations y11 y1 m x11 x 1p b11 b 1m y y x x b b n1
More informationAbstract. Introduction
Analysis of BEM and FEM codes for the modelling of microwave passive dielectric devices C. Di Nallo, F. Frezza, A. Galli Department ofelectrical Engineering, "La Saprienza " University of Rome, Via Eudossiana
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationApplications of Time Domain Vector Potential Formulation to 3-D Electromagnetic Problems
Applications of Time Domain Vector Potential Formulation to 3-D Electromagnetic Problems F. De Flaviis, M. G. Noro, R. E. Diaz, G. Franceschetti and N. G. Alexopoulos Department of Electrical Engineering
More informationLinear Algebra Primer
Introduction Linear Algebra Primer Daniel S. Stutts, Ph.D. Original Edition: 2/99 Current Edition: 4//4 This primer was written to provide a brief overview of the main concepts and methods in elementary
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationRECENT ADVANCES IN SIMULATION OF MAGNETRONS AND CROSSED-FIELD AMPLIFIERS. Abstract
RECENT ADVANCES IN SIMULATION OF MAGNETRONS AND CROSSED-FIELD AMPLIFIERS George E. Dombrowski 69 Birchwood Heights Road Storrs, Connecticut 06268 (860)429-2478 Abstract Various improvements and enhancements
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
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 informationA Plane Wave Expansion of Spherical Wave Functions for Modal Analysis of Guided Wave Structures and Scatterers
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 51, NO. 10, OCTOBER 2003 2801 A Plane Wave Expansion of Spherical Wave Functions for Modal Analysis of Guided Wave Structures and Scatterers Robert H.
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More informationECE 451 Advanced Microwave Measurements. Circular and Coaxial Waveguides
ECE 451 Advanced Microwave Measurements Circular and Coaxial Waveguides Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jesa@illinois.edu ECE 451 Jose Schutt Aine 1 Circular
More informationWKB solution of the wave equation for a plane angular sector
PHYSICAL REVIEW E VOLUME 58, NUMBER JULY 998 WKB solution of the wave equation for a plane angular sector Ahmad T. Abawi* Roger F. Dashen Physics Department, University of California, San Diego, La Jolla,
More informationMEASURE THE COMPLEX PERMEABILITY OF FER- ROMAGNETIC THIN FILMS: COMPARISON SHORTED MICROSTRIP METHOD WITH MICROSTRIP TRANS- MISSION METHOD
Progress In Electromagnetics Research Letters, Vol. 11, 173 181, 2009 MEASURE THE COMPLEX PERMEABILITY OF FER- ROMAGNETIC THIN FILMS: COMPARISON SHORTED MICROSTRIP METHOD WITH MICROSTRIP TRANS- MISSION
More informationMicrowave Network Analysis
Prof. Dr. Mohammad Tariqul Islam titareq@gmail.my tariqul@ukm.edu.my Microwave Network Analysis 1 Text Book D.M. Pozar, Microwave engineering, 3 rd edition, 2005 by John-Wiley & Sons. Fawwaz T. ILABY,
More informationMicrowave Engineering 3e Author - D. Pozar
Microwave Engineering 3e Author - D. Pozar Sections 3.6 3.8 Presented by Alex Higgins 1 Outline Section 3.6 Surface Waves on a Grounded Dielectric Slab Section 3.7 Stripline Section 3.8 Microstrip An Investigation
More informationACM 104. Homework Set 4 Solutions February 14, 2001
ACM 04 Homework Set 4 Solutions February 4, 00 Franklin Chapter, Problem 4, page 55 Suppose that we feel that some observations are more important or reliable than others Redefine the function to be minimized
More informationChapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations
Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2
More informationKevin James. MTHSC 3110 Section 2.1 Matrix Operations
MTHSC 3110 Section 2.1 Matrix Operations Notation Let A be an m n matrix, that is, m rows and n columns. We ll refer to the entries of A by their row and column indices. The entry in the i th row and j
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationChapter 2: Numeric, Cell, and Structure Arrays
Chapter 2: Numeric, Cell, and Structure Arrays Topics Covered: Vectors Definition Addition Multiplication Scalar, Dot, Cross Matrices Row, Column, Square Transpose Addition Multiplication Scalar-Matrix,
More informationLecture 2 INF-MAT : , LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky
Lecture 2 INF-MAT 4350 2009: 7.1-7.6, LU, symmetric LU, Positve (semi)definite, Cholesky, Semi-Cholesky Tom Lyche and Michael Floater Centre of Mathematics for Applications, Department of Informatics,
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationPhysics 202 Midterm Exam 2 Nov 2, 2011
Physics 202 Midterm Exam 2 Nov 2, 2011 Name: Student ID: Section: TA (please circle): Daniel Crow Scott Douglas Yutao Gong Taylor Klaus Aaron Levine Andrew Loveridge Jason Milhone Hojin Yoo Instructions:
More informationThe Singular Value Decomposition (SVD) and Principal Component Analysis (PCA)
Chapter 5 The Singular Value Decomposition (SVD) and Principal Component Analysis (PCA) 5.1 Basics of SVD 5.1.1 Review of Key Concepts We review some key definitions and results about matrices that will
More informationNUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group
NUMERICAL METHODS lor CHEMICAL ENGINEERS Using Excel', VBA, and MATLAB* VICTOR J. LAW CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup,
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
More information9.1 Mean and Gaussian Curvatures of Surfaces
Chapter 9 Gauss Map II 9.1 Mean and Gaussian Curvatures of Surfaces in R 3 We ll assume that the curves are in R 3 unless otherwise noted. We start off by quoting the following useful theorem about self
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More information