A semi-analytic technique to determine the propagation constant of periodically segmented Ti:LiNbO 3 waveguide

Size: px
Start display at page:

Download "A semi-analytic technique to determine the propagation constant of periodically segmented Ti:LiNbO 3 waveguide"

Transcription

1 Avalable onlne at Pelaga Research Lbrary Advances n Appled Scence Research, 011, (1): ISSN: CODEN (USA): AASRFC A sem-analytc technque to determne the propagaton constant of perodcally segmented T:LNbO wavegude Pranabendu Ganguly Advanced Technology Development Centre, Indan Insttute of Technology, Kharagpur, Inda ABSTRACT A sem-analytcal technque to determne the propagaton constant of perodcally segmented T:LNbO wavegude (PSW) s descrbed. The -D refractve ndex profles of the wavegude segments are computed from ts fabrcaton parameters, such as, T-layer thckness and wdth, and the dffuson parameters. WKB method s used to transform -D refractve ndex profle to 1- D lateral effectve ndex profle, whch s then converted to equvalent refractve ndex profle of the PSW. Fnally transfer matrx method s appled to compute the propagaton constant of the wavegude. Key words: T:LNbO, Perodcally segmented wavegude, Effectve-ndex, Equvalent refractve ndex, Matrx method. INTRODUCTION There has been an ncreasng amount of nterest n the applcaton of perodcally segmented wavegudes (PSW s) n ntegrated optcs. To date these devces have been reported n a number of dfferent materal systems ncludng proton-exchanged LNbO [1], KTP [], InP [], annealed proton exchanged LNbO [4-6] and T:LNbO [7-9]. Applcaton of PSWs ranges from nonlnear devces, whch uses quas phase matched second harmonc generaton [], to lnear devces lke a -D mode taper [10], asymmetrc drectonal coupler flter [11], and wavelength demultplexer [9]. In a PSW the ncrease n refractve ndex ( n) s modulated perodcally durng fabrcaton, as shown n Fg.1. Pelaga Research Lbrary 16

2 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): Fgure 1 Perodcally segmented wavegude (PSW). As a consequence of the segmentaton, the loss n the gude s ncreased and the effectve refractve ndex s reduced when compared to a contnuous wavegude. A PSW s characterzed by ts perod, Λ, and duty-cycle, η (the rato of the length of a segment and the perod of the gude). It has been demonstrated that a PSW can be represented by an equvalent contnuous wavegude wth same depth and wdth, n whch the average ndex dfference ( n eq ) s taken to be the weghted average of the ndex along the propagaton drecton. Ths s represented by equaton (1), [1] n eq = η n (1) By choosng the duty-cycle, the average ndex can be spatally modfed along the wavegude. The refractve ndex change determnes the mode sze, propagaton constant and cut-off wavelength of the PSW. Perodcally segmented wavegudes formed by T-ndffuson n a LNbO substrate wll have a graded concentraton varaton along depth and lateral drectons. In ths work, the propagaton constants of perodcally segmented T:LNbO wavegudes have been computed by applyng the effectve-ndex-based matrx method (EIMM) along wth equvalent wavegude concept. In the frst step, the depth and lateral refractve ndex profles of each T:LNbO segment has been computed drectly from ts fabrcaton parameters, such as, T-layer thckness and wdth, and the dffuson parameters. Effectve ndex method has been appled to transfer the -dmentsonal ndex dstrbuton to 1-dmensonal lateral effectve-ndex dstrbuton. In the next step, the average refractve ndex has been computed to transfer the PSW nto an equvalent contnuous wavegude. Fnally, the transfer matrx method s used to compute the propagaton constant of the equvalent wavegude. Theoretcal Approach Determnaton of effectve refractve ndex of T:LNbO PSW In the case of T:LNbO wavegudes the T-concentraton profle can be represented by [1] 1 τ z τ z C ( x, z) = C o erf erf 4 d z d z Pelaga Research Lbrary 17

3 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): W x W x erf erf () d x d x where C o s the sold solublty of T nto LNbO, W and τ are the deposted T strp wdth and thckness, and d x and d z are the dffuson lengths along x and z axes of the crystal The refractve ndces of ordnary and extraordnary ray, n o and n e, of the congruently grown LNbO crystal, at room temperature (5 o C), can be obtaned usng the modfed Sellmeer equatons [14] n = λ λ o () ( ) λ n e = λ (4) ( ) where λ s the wavelength n µm. The refractve ndex change nduced by T-ndffuson for both ordnary and extraordnary rays are related to the ttanum concentraton for LNbO as [1] α o, n x z, λ = A ( λ) C ( x, z) (5) o, e ( ) [ ] e, o, e where A o,e are dependent on λ, and α o,e are 0.5 and 0.85 for ordnary and extraordnary rays. Detals of the λ dependence of A o,e are gven n ref.[1], whch are vald wthn the wavelength range 0.6 λ (µm) 1.6. So by usng equatons (-5) one can determne the -dmensonal refractve ndex dstrbuton of the segment of a T:LNbO PSW drectly from ts fabrcaton parameters. Now a wavegude mode exsts only f the total transverse phase shft (along Z) for one round-trp across the gude equals an ntegral multple of п,.e., [15] Zb 1/ [ n ( z, x) n ( x) ] dz φ φ = mπ 4π eff t b λ 0 m = 0,1,,K (6) Ths refers to the wave confnement over the YZ plane at a partcular lateral poston (x). The frst term s the total phase change of the wave as t travels between the crystal surface z = 0 and the return pont z = Z b n the bulk and back. ф t and ф b are Goos-Hanchen phase shfts correspondng to total nternal reflecton at z = 0 and z = Z b, respectvely. For the LNbO -ar nterface, ф t s about -0.9п and -0.98п for TE and TM modes, respectvely, and ф b approaches (-п/) for both polarzatons [15]. Equaton (6) s known as the WKB quantzaton condton. Ths leads to a set of dscrete angles of propagaton correspondng to dfferent guded modes. So by solvng equaton (6) numercally for the fundamental mode (m = 0) one can determne the effectve refractve ndex, n eff (x), of the sngle-mode T:LNbO wavegude. For the computaton of refractve ndex of equvalent wavegude for a PSD wth duty-cycle η one has to use equaton (1),.e., n eq ( x) = η n ( x) (7) eff In the next step, ths n eg (x) s approxmated by a starcase-type step-ndex profle and matrx method s appled to that layered structure. Pelaga Research Lbrary 18

4 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): Determnaton of propagaton constant The propagaton constants of the T:LNbO PSW s are here determned by matrx method. If the medum s consdered to be made of a number of layers (each of thckness d and refractve ndex n ), then for an ncdent plane wave n the frst layer, the electrc feld assocated wth each layer can be represented as [15] E = u E exp( j ) exp j ωt ( K cosθ ) x β y where E and u E [ { }] j )exp[ j{ ωt ( K cosθ ) x β y} ] exp( (8) 1 = = 0 ; = K d cosθ ; 4 = K 4 ( d d ) cosθ 4 ; K ( d d d )cos π = L 1 θ ; K = K o n = n ; λ β = K snθ. (9) E are the magntude of electrc feld vectors assocated wth the downward and upward-propagatng waves, respectvely, and u and u are ther correspondng unt vectors, and β s the propagaton constant along the layers and s the same for all the layers, whch s consstent wth the law of refracton. Applyng the boundary condtons at the nterfaces, one obtans E 1 E E N = = = S 1 L S1 S LS N 1 (10) E1 E E N where 1 exp( jδ ) r exp( jδ ) S 1 = (11) t r exp( jδ ) exp( jδ ) π δ = n d cosθ λ, =,,4,,N. o S -1 s the transmsson matrx from the (-1)-th and -th layer; r and t represent the ampltude reflecton and transmsson coeffcents, respectvely, from the (-1)-th layer to the -th layer. Usng equatons (8-11) one can determne the feld ampltude n any layer n terms of the ncdent feld ampltude E. 1 The above method when appled to equvalent refractve ndex profle of equaton (7) can be used to determne the propagaton constant of T:LNbO perodcally segmented wavegudes. Prsm couplng approach has been appled for the purpose. Intally a hgher refractve ndex layer s consdered as layer 1 and exctaton effcency of the gudng layer s computed by matrx method for dfferent ncdent angles. Out of the dfferent ncdent angles only a sngle wll excte the guded mode n the sngle-mode wavegude refractve ndex structure. Exctaton effcency versus ncdent angle or propagaton constant characterstcs wll show only one peak for the guded mode, from whch propagaton constant of the PSW can be determned. The detals of the matrx method are gven n ref. [16]. Pelaga Research Lbrary 19

5 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): Computatonal Procedure and Results Computatonal Procedure The computer program wrtten for calculatng propagaton constant of T:LNbO PSW requres the followng parameters as nputs: 1. T flm parameters: thckness, strp wdth, gap between the segments, and length of segments,. T dffuson parameters: dffuson temperature, tme, and ambent,. Beam parameters: wavelength and polarzaton (TE / TM), 4. Computatonal parameters: samplng area (x range, z range) and the layer thcknesses of the dscretsed n eq (x) profle. The frst three sets of parameters are taken as the fabrcaton parameters of T:LNbO PSW, whle the x range and z range are taken as -15 to 15 µm and 0 to 15 µm, respectvely, for most computatons. A unform layer thckness of 0.01 µm has been used. However, n order to reduce the computaton tme wthout loss of accuracy, nonunform parttonng of n eq (x) nvolvng a fewer number of layers may also be consdered. The computatonal steps to determne the propagaton constants of T:LNbO PSWs are as follows: 1. Computaton of two-dmensonal refractve ndex profles, n(x, z), of T:LNbO wavegude from ts fabrcaton parameters.. Computaton of one-dmensonal lateral effectve ndex profle, n eff (x), from the twodmensonal refractve ndex dstrbuton by WKB method.. Determnaton of equvalent refractve ndex profle, n eq (x), for the PSW from n eff (x). 4. Applcaton of matrx method onto n eq (x) to determne the propagaton constant of T:LNbO PSW. Followng the matrx method outlned n prevous secton, the exctaton effcency s computed for varous values of ncdence angle θ 1 from frst hgh refractve ndex medum (added) for θ c θ 1 90 o to fnd the values of θ 1 and hence β [usng equaton (9)] for whch exctaton effcency shows sharp resonance peaks. The gap between the added medum and the boundary of the wavegude s ncreased untl the lmtng values of β are obtaned. The flow chart of the entre computatonal procedure s gven n Fg.. Pelaga Research Lbrary 140

6 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): User Screen START Computaton of r.. of the wavegude from ts fabrcaton parameters Converson of D r.. profle to 1D effectve ndex profle by WKB method Computaton of n eq (x) from the effectve ndex profle INPUTS x, n 1, d Determnaton of θ C for the structure θ = θ C Add ncrement to θ Generaton of exctaton effcency versus β data No Is θ 90 o? Yes Determnaton of exctaton effcency versus β values for the resonance peak OUTPUT β value of PSW Fgure Flow-chart for computaton of the propagaton constant of T:LNbO PSW. Pelaga Research Lbrary 141

7 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): RESULTS AND DISCUSSION The computed propagaton constant versus gap between the segments of T:LNbO PSW s shown n Fg. for TM mode at 1.1 µm transmttng wavelength. The Propagaton Constant (µm -1 ) Length of each segment = 10 µm T-thckness = µm T-wdth = 6. µm Dffuson temperatute = 1050 o C Dffuson tme = 6.0 hours Wavelength = 1.1 µm, TM mode Gap between the segments (µm) Fgure Computed propagaton constant versus gap between the segments of a T:LNbO PSW. Pelaga Research Lbrary 14

8 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): Segment length = 10.0 µm Gap between the segments = 9.5 µm Exctaton Effcency Propagaton Constant (µm -1 ) Fgure 4 Typcal resonance peak of exctaton effcency versus propagaton constant plot of T:LNbO PSW. length of the segments s taken as 10 µm. It may be noted from the fgure that as the gap between the segments ncreases the propagaton constant decreases. Ths nature of the characterstc s obvous from the fact that as gap between the segments ncreases, the n eq of PSW decreases and hence propagaton constant of the wavegude decreases. Typcal exctaton effcency versus propagaton constant plot for gap 9.5 µm s shown n Fg.4, whch shows a sharp peak ndcatng the guded mode propagaton constant of the PSW. It may be noted that the present analytcal model s usable for both TE and TM polarzed lghts for dfferent sets of PSW fabrcaton parameters wthn a wavelength range from 0.6 to 1.6 µm. The method nvolves only straghtforward multplcaton of x transfer matrces of the layers. It does not requre the soluton of any transcendental or dfferental equaton. Iteratons are used only once to evaluate the effectve refractve ndex of the wavegude. Ths makes the present model very fast and smple to use wth a PC. The computed propagaton constant wll be an essental nput for desgn consderaton for several PSW components of optcal ntegrated crcuts, such as asymmetrc drectonal couplers and Y-junctons. Pelaga Research Lbrary 14

9 Pranabendu Ganguly Adv. Appl. Sc. Res., 011, (1): CONCLUSION Propagaton constants of T:LNbO perodcally segmented wavegudes (PSW s) are determned from ts fabrcaton parameters by effectve ndex based matrx method (EIMM) along wth equvalent contnuous wavegude model. The -D refractve ndex profles of the wavegude segments are computed from ts fabrcaton parameters. WKB method s used to transform -D refractve ndex profle to 1-D lateral effectve ndex profle, whch s then converted to equvalent refractve ndex profle of the PSW. Fnally transfer matrx method s appled to compute the propagaton constant of the wavegude. Ths sem-analytcal technque s computatonally fast and smple to use wth a PC. REFERENCES [1] E.B. Pun, K.K. Wong, I. Andonovc, P.J.R. Laybourn, and R. De La Rue, Electron. Lett., 18, pp , 198. [] J.D. Berlen, D.B. Laubacher, and J.B. Brown, Appl. Phys. Lett., 56, pp , [] F. Dorgeulle, B. Mersal, S. Francos, G. Herve-Gruyer, and M. FlocheOpt. Lett., 0, pp , [4] M.H. Chou, M.A. Arbone, and M.M. Fejer, Opt. Lett., 1, pp , [5] K. Thyagarajan, C.W. Chen, R.V. Ramaswamy, H.S. Km, and H.C. Cheng, Opt. Lett., 19, pp , [6] D. Nr, S. Ruschn, A. Hardy, and D. Brooks, Electron. Lett., 1, pp , [7] D. Nr, Z. Wessman, S. Ruschn, and A. Hardy, Opt. Lett., 19, pp , [8] J. Webjorn, F. Laurell, and G. Arvdsson, J. Lghtwave Technol., 7, pp , [9] Z. Wessman, D. Nr, S. Ruschn, and A. Hardy, Appl. Phys. Lett., 67, pp.0-04, [10] Z. Wessman and I. Hendel, J. Lghtwave Technol., 1, pp , [11] Z. Wessman, F. Sant-Andre, and A. Kevorkan, Proc. ECIO 97, pp.5-55, (Aprl 1997). [1] D. Ortega, R.M. De La Rue, and J.S. Atchson, J. Lghtwave Technol., 16, pp.84-91, [1] P. Ganguly, D.C. Sen, S. Dutt, J.C. Bswas, and S.K. Lahr, Fber and Integrated Optcs, 15, pp , [14] D.S. Smth and H.D. Rccus, Opt. Commun., 17, pp.-5, [15] P. Ganguly, J.C. Bswas, and S.K. Lahr, Fber and Integrated Optcs, 17, pp , [16] A.K. Ghatak, K. Thyagarajan, and M.R. Shenoy, J. Lghtwave Technol., 5, pp , Pelaga Research Lbrary 144

Boundaries, Near-field Optics

Boundaries, Near-field Optics Boundares, Near-feld Optcs Fve boundary condtons at an nterface Fresnel Equatons : Transmsson and Reflecton Coeffcents Transmttance and Reflectance Brewster s condton a consequence of Impedance matchng

More information

CHAPTER II THEORETICAL BACKGROUND

CHAPTER II THEORETICAL BACKGROUND 3 CHAPTER II THEORETICAL BACKGROUND.1. Lght Propagaton nsde the Photonc Crystal The frst person that studes the one dmenson photonc crystal s Lord Raylegh n 1887. He showed that the lght propagaton depend

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and

More information

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011 Interface Energy Couplng between -tungsten Nanoflm and Few-layered Graphene Meng Han a, Pengyu Yuan a, Jng Lu a, Shuyao S b, Xaolong Zhao b, Yanan Yue c, Xnwe Wang a,*, Xangheng Xao b,* a 2010 Black Engneerng

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Summer Term 01, Abbe School of Photoncs, FSU Jena, Prof. Thomas Pertsch Computatonal Photoncs Semnar 03, 7 May 01 Implementaton of the Matr Method calculaton of the transfer matr

More information

Effect of Losses in a Layered Structure Containing DPS and DNG Media

Effect of Losses in a Layered Structure Containing DPS and DNG Media PIERS ONLINE, VOL. 4, NO. 5, 8 546 Effect of Losses n a Layered Structure Contanng DPS and DNG Meda J. R. Canto, S. A. Matos, C. R. Pava, and A. M. Barbosa Insttuto de Telecomuncações and Department of

More information

ECE 107: Electromagnetism

ECE 107: Electromagnetism ECE 107: Electromagnetsm Set 8: Plane waves Instructor: Prof. Vtaly Lomakn Department of Electrcal and Computer Engneerng Unversty of Calforna, San Dego, CA 92093 1 Wave equaton Source-free lossless Maxwell

More information

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices Supplementary Informaton for Observaton of Party-Tme Symmetry n Optcally Induced Atomc attces Zhaoyang Zhang 1,, Yq Zhang, Jteng Sheng 3, u Yang 1, 4, Mohammad-Al Mr 5, Demetros N. Chrstodouldes 5, Bng

More information

16 Reflection and transmission, TE mode

16 Reflection and transmission, TE mode 16 Reflecton transmsson TE mode Last lecture we learned how to represent plane-tem waves propagatng n a drecton ˆ n terms of feld phasors such that η = Ẽ = E o e j r H = ˆ Ẽ η µ ɛ = ˆ = ω µɛ E o =0. Such

More information

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

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

( ) + + REFLECTION FROM A METALLIC SURFACE

( ) + + REFLECTION FROM A METALLIC SURFACE REFLECTION FROM A METALLIC SURFACE For a metallc medum the delectrc functon and the ndex of refracton are complex valued functons. Ths s also the case for semconductors and nsulators n certan frequency

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

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

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

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

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

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

Homework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich Homework 4 Contact: frmmerm@ethz.ch Due date: December 04, 015 Nano Optcs, Fall Semester 015 Photoncs Laboratory, ETH Zürch www.photoncs.ethz.ch The goal of ths problem set s to understand how surface

More information

Fresnel's Equations for Reflection and Refraction

Fresnel's Equations for Reflection and Refraction Fresnel's Equatons for Reflecton and Refracton Incdent, transmtted, and reflected beams at nterfaces Reflecton and transmsson coeffcents The Fresnel Equatons Brewster's Angle Total nternal reflecton Power

More information

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton

More information

Supporting Information

Supporting Information Supportng Informaton Water structure at the ar-aqueous nterface of dvalent caton and ntrate solutons Man Xu, Rck Spnney, Heather C. Allen* allen@chemstry.oho-state.edu Fresnel factors and spectra normalzaton

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

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

Note: Please use the actual date you accessed this material in your citation. MIT OpenCourseWare http://ocw.mt.edu 6.13/ESD.13J Electromagnetcs and Applcatons, Fall 5 Please use the followng ctaton format: Markus Zahn, Erch Ippen, and Davd Staeln, 6.13/ESD.13J Electromagnetcs and

More information

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

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

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

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

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

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal

More information

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

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Grover s Algorithm + Quantum Zeno Effect + Vaidman Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the

More information

Problem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that

Problem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that PHY 387 K. Solutons for problem set #7. Problem 1: To prove that under the assumptons at hand, the group velocty of an EM wave s less than c, I am gong to show that (a) v group < v phase, and (b) v group

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

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

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Interconnect Modeling

Interconnect Modeling Interconnect Modelng Modelng of Interconnects Interconnect R, C and computaton Interconnect models umped RC model Dstrbuted crcut models Hgher-order waveform n dstrbuted RC trees Accuracy and fdelty Prepared

More information

Introductory Optomechanical Engineering. 2) First order optics

Introductory Optomechanical Engineering. 2) First order optics Introductory Optomechancal Engneerng 2) Frst order optcs Moton of optcal elements affects the optcal performance? 1. by movng the mage 2. hgher order thngs (aberratons) The frst order effects are most

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Theoretical and experimental investigations of macro-bend Losses for standard single mode fibers

Theoretical and experimental investigations of macro-bend Losses for standard single mode fibers Theoretcal and epermental nvestgatons of macro-bend Losses for standard sngle mode fbers Qan Wang Gerald Farrell and Thomas Frer ppled Optoelectroncs Center ubln Insttute of Technology Kevn Street ubln

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Modeling acoustic transducer surface waves by Transmission Line Matrix method

Modeling acoustic transducer surface waves by Transmission Line Matrix method Modelng acoustc transducer surface waves by Transmsson Lne Matrx method Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Peter-Chrstan Eccardt Semens AG, CT MS, München Wllam

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Title: Radiative transitions and spectral broadening

Title: Radiative transitions and spectral broadening Lecture 6 Ttle: Radatve transtons and spectral broadenng Objectves The spectral lnes emtted by atomc vapors at moderate temperature and pressure show the wavelength spread around the central frequency.

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

More information

Color Rendering Uncertainty

Color Rendering Uncertainty Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:

More information

A Fast Computer Aided Design Method for Filters

A Fast Computer Aided Design Method for Filters 2017 Asa-Pacfc Engneerng and Technology Conference (APETC 2017) ISBN: 978-1-60595-443-1 A Fast Computer Aded Desgn Method for Flters Gang L ABSTRACT *Ths paper presents a fast computer aded desgn method

More information

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

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Polynomial Regression Models

Polynomial Regression Models LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

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

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

More information

Chapter 4. Velocity analysis

Chapter 4. Velocity analysis 1 Chapter 4 Velocty analyss Introducton The objectve of velocty analyss s to determne the sesmc veloctes of layers n the subsurface. Sesmc veloctes are used n many processng and nterpretaton stages such

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Multiline Overview. Karl Bois, PhD Lead SI Engineer Servers & Blades. December 5th, 2016

Multiline Overview. Karl Bois, PhD Lead SI Engineer Servers & Blades. December 5th, 2016 Multlne Overvew Karl Bos, PhD Lead SI Engneer Servers & Blades December 5th, 216 Why are we here? Orgnal Equpment Manufacturers (OEMs) and Orgnal Desgn Manufacturers (ODMs) have desgns that requre that

More information

Electron-Impact Double Ionization of the H 2

Electron-Impact Double Ionization of the H 2 I R A P 6(), Dec. 5, pp. 9- Electron-Impact Double Ionzaton of the H olecule Internatonal Scence Press ISSN: 9-59 Electron-Impact Double Ionzaton of the H olecule. S. PINDZOLA AND J. COLGAN Department

More information

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES Abdelkader Benchou, PhD Canddate Nasreddne Benmoussa, PhD Kherreddne Ghaffour, PhD Unversty of Tlemcen/Unt of Materals

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Rate of Absorption and Stimulated Emission

Rate of Absorption and Stimulated Emission MIT Department of Chemstry 5.74, Sprng 005: Introductory Quantum Mechancs II Instructor: Professor Andre Tokmakoff p. 81 Rate of Absorpton and Stmulated Emsson The rate of absorpton nduced by the feld

More information

( ) = ( ) + ( 0) ) ( )

( ) = ( ) + ( 0) ) ( ) EETOMAGNETI OMPATIBIITY HANDBOOK 1 hapter 9: Transent Behavor n the Tme Doman 9.1 Desgn a crcut usng reasonable values for the components that s capable of provdng a tme delay of 100 ms to a dgtal sgnal.

More information

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features Analyss of the Magnetomotve Force of a Three-Phase Wndng wth Concentrated Cols and Dfferent Symmetry Features Deter Gerlng Unversty of Federal Defense Munch, Neubberg, 85579, Germany Emal: Deter.Gerlng@unbw.de

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Over-Temperature protection for IGBT modules

Over-Temperature protection for IGBT modules Over-Temperature protecton for IGBT modules Ke Wang 1, Yongjun Lao 2, Gaosheng Song 1, Xanku Ma 1 1 Mtsubsh Electrc & Electroncs (Shangha) Co., Ltd., Chna Room2202, Tower 3, Kerry Plaza, No.1-1 Zhongxns

More information

1.4 Small-signal models of BJT

1.4 Small-signal models of BJT 1.4 Small-sgnal models of J Analog crcuts often operate wth sgnal levels that are small compared to the bas currents and voltages n the crcut. Under ths condton, ncremental or small-sgnal models can be

More information

Frequency dependence of the permittivity

Frequency dependence of the permittivity Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM Smulaton of 2D Elastc Bodes wth Randomly Dstrbuted Crcular Inclusons Usng the BEM Zhenhan Yao, Fanzhong Kong 2, Xaopng Zheng Department of Engneerng Mechancs 2 State Key Lab of Automotve Safety and Energy

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

AGC Introduction

AGC Introduction . Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Lecture 6. P ω ω ε χ ω ω ω ω E ω E ω (2) χ ω ω ω χ ω ω ω χ ω ω ω (2) (2) (2) (,, ) (,, ) (,, ) (2) (2) (2)

Lecture 6. P ω ω ε χ ω ω ω ω E ω E ω (2) χ ω ω ω χ ω ω ω χ ω ω ω (2) (2) (2) (,, ) (,, ) (,, ) (2) (2) (2) Lecture 6 Symmetry Propertes of the Nonlnear Susceptblty Consder mutual nteracton of three waves: ω, ω, ω = ω + ω 3 ω = ω ω ; ω = ω ω 3 3 P ω ω ε ω ω ω ω E ω E ω n + m = 0 jk m + n, n, m j n k m jk nm

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

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

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA 14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton

More information

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL The Synchronous 8th-Order Dfferental Attack on 12 Rounds of the Block Cpher HyRAL Yasutaka Igarash, Sej Fukushma, and Tomohro Hachno Kagoshma Unversty, Kagoshma, Japan Emal: {garash, fukushma, hachno}@eee.kagoshma-u.ac.jp

More information

Supplementary Materials for

Supplementary Materials for advances.scencemag.org/cg/content/full/2/7/e1600304/dc1 Supplementary Materals for Interface-drven topologcal Hall effect n SrRuO3-SrIrO3 blayer Jobu Matsuno, Naok Ogawa, Kenj Yasuda, Fumtaka Kagawa, Wataru

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information