Modeling acoustic transducer surface waves by Transmission Line Matrix method
|
|
- Lisa Hardy
- 5 years ago
- Views:
Transcription
1 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 O Connor Unversty College Dubln Summary A model for the generaton and propagaton of acoustc and surface waves s mplemented employng the Transmsson Lne Matrx (TLM) method wth a locally reactng surface restorng boundary condton, whch s a model for an array of capactve ultrasound transducers. Frst test calculatons clearly show the ablty of the proposed model to capture acoustc and surface wave effects. The comparson wth a reference soluton derved wth FE methods reveals that further refnements may be necessary although general agreement s acheved. Keywords TLM, surface waves, ultra sound, transducer array on FEM Technology 1
2 1 Introducton Mcro-mechanc ultrasound transducer arrays are usually bult on the bass of pezo electrc bulk resonators. Alternatve transducer prncples, especally capactve transducers are not employed because bulk resonators offer hgher sound pressure levels, however wth the drawback of a small bandwdth. In the last few years technologes were developed that allow to produce capactve transducers wth extremely small gap szes (gap sze means dstance between the membrane and the counter electrode) on the bass of the standard BCMOS process. Ths small gap enables hgh sound pressure levels and a large bandwdth at reasonable load voltages (about 10 V). Arrays of capactve transducers consst of a large number of membranes next to each other whch are set n moton by an alternatng current. Experments have shown that apart from sound waves n the plane of the transducer array surface waves may be generated and propagate, much lke the gravty waves at the surface of a lake. Ths surface waves have long rng down tmes and may cause serous errors f the transducer array s used as a mcrophone. Therefor a method s needed to calculate the generaton and propagaton of these surface waves n order to develop and test strateges to suppress the unwanted surface effects. The Fnte Element Method (FEM) s a well-establshed tool for the numercal soluton of partal dfferental equatons. One dsadvantage of the method results from the fact that the numercal effort tends to grow as the square of the number of nodes whch dscretze the problem. Thus the demand of a three-dmensonal calculaton grows wth the power of sx of a typcal length. Unfortunately to resolve propagatng waves the node densty must not fall below a certan lmt whch prevents an enlargement of the element sze far from the orgn. Ths constrant makes the calculaton of larger acoustc felds mpossble. On the other hand, for wave propagaton other numercal schemes exst whch requre computatonal resources drectly proportonal to the number of nodes. In ths case the total cost of the computaton s proportonal to the power of three of a typcal length of the problem n three dmensons. One such scheme s the Transmsson Lne Matrx (TLM) method. The theoretcal foundatons of the TLM method date back to the 40 s when Whnnery et al. and Kron developed equvalent crcut models made up of dscrete elements for Maxwell s equaton to study wave gude models usng network analyzers [1,, 3]. Wth the advent of more powerful computers the TLMmethod was made popular by Johns and Beurle n 70s. Snce then the method has been appled successfully to many problems whch are descrbed by wave or dffuson equatons. The man dfference between the TLM ansatz and standard schemes such as Fnte Dfference or FE methods s as follows: Usually a physcal problem s descrbed precsely by a partal dfferental equaton whch subsequently s dscretzed. The soluton to ths equaton s then approxmated numercally. The TLM approach s to approxmate the physcal process (e.g. wave propagaton) by dscrete physcal models, whch can then be calculated exactly. Ths approach has several advantages. The physcal models may be programmed ntutvely, all dscretzaton approxmatons are apparent, and the calculaton may be parallelzed easly by doman decomposton. The dsadvantages are that TLM s best suted to equdstant nodes and fxed tme steps, whch are coupled to the spatal dstance of the nodes. However, n recent years TLM was adapted for non-equdstant and curvlnear grds. For complex techncal devces such as acoustc transducer arrays t s requred to compute several coupled physcal effects (electro-magnetcs, structural mechancs, flud structure nteracton) whle large sound felds must be mapped. One approach s to couple TLM and FEM usng FEM to calculate the structural and electromagnetc effects and usng TLM to model the acoustc flud. For ths study a major problem was to smulate surface waves whch are generated at the boundares of the actve area of a transducer array and whch then propagate along the surface of the transducer array. These surface waves decay slowly and may cause sgnfcant deflectons of the transducer membranes whch are not wanted as they are not drectly related to the acoustc feld. As a frst step the ablty of the TLM-model to on FEM Technology
3 smulate surfaces waves on a transducer array coupled to a flud s studed. The next secton descrbes the TLM theory and algorthm n more detal. Secton 3 deals wth the couplng of TLM and a smple model for the transducer array, whch s to be replaced by FEM models n future. In the followng secton the current status of the work s reported by gvng an example of a flud/structure couplng. TLM-Theory The TLM scheme may be derved n two dfferent ways. One way s to derve t as a specal and smplfed case of the more general Lattce Boltzmann schemes for the soluton of the Naver-Stokes equatons as was demonstrated by Chopard et al. [4]. Though ths approach opens the way to ncluson of many nonlnear effects of sound propagaton t s not very ntutve. Therefore another dervaton shall be presented here whch reles on a physcal model of the feld effects and was proposed by Kagawa for the acoustc feld [5, 6]. The basc dea s to follow Huygens prncple that states, that each pont on a wavefront can be consdered as a source of secondary wavelets, so that the propagaton can be consdered as a contnuous process of scatterng of elementary component waves. The superposton of all scattered (elementary) waves makes up the wave. Ths mechansm may be dscretzed as follows: The flud s replaced by a network of ppes wth equal dameters whch are jont at fxed dstances x (see fgure 1). The junctons are called nodes. Wthn the ppes acoustc pressure pulses may travel undsturbed at the speed of sound,.e. no dampng or boundary effects are consdered. The dameter of the ppes s assumed to be much smaller than the shortest wavelength of nterest so hgher order ppe modes are neglected. From ths follows, that f a set of acoustc pressure pulses s released at the same tme at arbtrary nodes, these pulses travel through the ppes and reach ther neghborng nodes after a tme nterval of t = x/c 0 wth c 0 beng the speed of sound n the ppe. For dmensons, the juncton may be seen as a meetng of the ppe carryng the ntal, ncomng pressure pulse wth three outgong ppes. Thus the acoustc pressure pulse arrvng at the node sees an apparent enlargement of the cross secton by a factor of three, whch means a correspondng change of the mpedance. At ths mpedance msmatch the ncomng pressure pulse s reflected wth half of the ntal ampltude and a reversal of sgn, whle t s transmtted nto the other three ppes wth half the ntal ampltude and no change of sgn (see fgure 1). From ths pont on the process of travelng n the ppes and scatterng at the nodes s repeated, wth four pulses arrvng smultaneoulsy at each node and each scatterng nto the other three ppes. As the process of scatterng and travelng s lnear, superposton apples and the net result s the sum of the component effects. The algorthm descrbed above can formally wrtten as follows: f ( x + x, t + t) = S j f ( x, t) (1) wth f ( x + x, t + t) representng the ampltude of the pressure pulse that travels from the node at poston x to the node at poston x + x between tmes t and t + t. x s the vector pontng to the nearest neghbor n the -th drecton, and S j s the scatterng matrx S j = Because of the smple structure of the scatterng matrx equaton (1) may also be wrtten as f ( x + x, t + t) = 1 f j ( x, t) f ( x, t) () j on FEM Technology 3
4 % % % & & & % % % & & & % % % & & & x ' ' ' ( ( ( ' ' ' ( ( ( ' ' ' ( ( ( : : : : : : : : : / / / / / / / / / !!! " " "!!! " " "!!! " ; ; ; ; < < < < ; ; ; ; < < < < ; ; ; ; < < < < before scatterng p after scatterng p/ p/ p/ p/ Fgure 1: The acoustc medum s replaced by a network of ppes, n whch pressure waves travel. The junctons are called nodes, the vectors x pont from one node to the nearest neghbors. The scatterng process s depcted n the lower part of the fgure. where s defned such that x = x. The acoustc pressure p ak and the velocty v are calculated from f ( x, t) applyng the followng formulae: p ak ( x, t) = 1 f ( x, t) (3) v( x, t) = 1 ρc x t f ( x, t) (4) Ths scheme s numercally equvalent to a second order Fnte Dfference scheme for the wave equaton, whch s shown as follows: The acoustc pressure at a poston x and tme t + t s gven as Repeated use of eq.(3) yelds p ak ( x, t + t) = 1 f ( x, t + t) p ak ( x, t + t) = 1 p ak ( x + x, t) f ( x + x, t) = 1 p ak ( x + x, t) p ak ( x, t t) + f ( x, t t) Ths may be rearranged to gve p ak ( x, t + t) p ak ( x, t) + p ak ( x, t t) = 1 p ak ( x + x, t) p ak ( x, t) whch equvalent to the correspondng Fnte Dfference equaton f x t = c on FEM Technology 4
5 speed of sound 1500 m/s flud densty 1000 kg/m 3 sze of calculaton area 4 4 mm number of nodes = Spatal resoluton x = 5 µm temporal resoluton t = s membrane complance Pa/m membrane mass kg/m frequency of exctaton 1 MHz Table 1: Parameters of the test calculaton of surface waves. Ths last result means that the speed of wave propagaton c s gven by the rato of the spatal and temporal resoluton. Note that the overall wave speed c s lower by the factor 1/ than the speed of component pressure pulses travelng n the ppes, whch s a feature of the method that s surprsng when frst encountered. 3 The ultrasound transducer array model The ultrasound transducers to be modeled here are capactve transducers made up of a flexble membranes above a fxed counter electrode. If a voltage s appled the membrane s attracted to the counter electrode untl the electrostatc force s balanced by the restorng force of the membrane stffness. Varatons of the magntude of the appled voltage lead to movements of the membrane whch then acts as a sound generator for the flud outsde the transducer. The transducer s modeled as a smple masssprng system where the sprng wth sprng constant D provdes the restorng force and the mass m corresponds to the membrane mass of the real transducer. From the flud sde the acoustc pressure p ak calculated usng eq.(3) acts on the membrane wth area A, also an electrostatc force F el may be appled. The equaton of moton for the membrane s m s = p ak A Ds + F el where s s the deflecton of the membrane. Ths equaton of moton s then solved employng the Newmark tme ntegraton scheme for every membrane and the resultng membrane velocty s appled as a boundary condton to the flud by evaluaton of eq.(4). 4 Results As a frst test a two dmensonal setup s consdered, whch conssts of water as the acoustc medum. The calculaton area s 4 4 mm, three of the surroundng boundares are rgd walls. The membrane array s placed along the bottom boundary. In ths case one membrane model per TLM-node s appled and the membrane mass s chosen to be very small whch corresponds to a free surface wth locally reactng restorng force. Further detals are gven n table 1. The frst 80 membranes on the left are excted wth a snusodal force wth 1 MHz whle all other membranes are allowed to move freely. After turnng on the exctng force two effects clearly are observed: a sound wave s generated by the excted transducer and propagates nto the flud a surface wave s generated that spreads along the boundary of the flud wth lower propagaton velocty than the sound wave on FEM Technology 5
6 Fgure : Pressure feld at tme t = s after turnng on the exctaton force ( MHz). Whte and black colors correspond to hgh and low pressure, respectvely. The surface wave that propagates along the bottom boundary can clearly be seen to be slower (have shorter wavelength) than the sound waves, that spread nto the medum. Repeated experments wth other exctaton frequences show that the typcal dsperson relatons for sound waves c = λf and surface waves c λ are satsfed by ether wave type. Fgure shows the pressure feld at the tme t = s after turnng on the exctaton force, here wth a frequency of MHz. In order to compare the results quanttatvely the same setup s also smulated usng ANSYS. The spatal and temporal resoluton s changed to x = 0µm and t = s to save calculaton tme. Fg. 3 shows the membrane deflectons as a functon of poston calculated wth ANSYS and TLM at tme t = s after turnng on the exctaton force. Although the general shape agrees well for both calculatons the ampltudes show dfferences of up to 0%. The reason for ths dfferences s not yet understood. Although ts too early to compare calculaton tmes t s noteworthy that the TLM-program takes about mnutes on 800x800 grd whle ANSYS needs 0 mnutes to complete on a 100x100 grd. 5 Summary and Conclusons A model for the calculaton of acoustc and surface waves n fluds s establshed on the bass of the Transmsson Lne Matrx (TLM) method. The free surface boundary condton s modeled by masssprng-systems whch are coupled to acoustc feld by the acoustc pressure and velocty. Frst results clearly demonstrate the ablty of ths scheme to capture acoustc and surface wave propagaton effects. Dfferences between the TLM soluton and a reference soluton calculated wth the FE method ndcate, that further checks or refnements may be necessary. The hgh computatonal speed of the TLM-scheme allows 3d problems to be modeled. A 3d verson of the TLM scheme was coded recently but the results have not yet been verfed. Frst tests show that t s possble to handle problems wth up to 4 mllon nodes on a contemporary PC n reasonable tme ( 1 s CPU tme per model tme step). on FEM Technology 6
7 1e 11 FEM TLM 5e 1 0 Deflecton n m 5e 1 1e e 11 e 11.5e Poston n m Fgure 3: Deflecton of the membranes as a functon of poston calculated wth ANSYS (sold lne) and TLM (dashed lne) at t = s after turnng on the exctaton force. References [1] G. Kron, Equvalent crcut of the feld equatons of Maxwell, Proc IRE, vol. 3, pp , [] J. Whnnery and S. Ramo, A new approach to the soluton of hgh-frequency problems, Proc IRE, vol. 3, pp , [3] J. Whnnery and S. Ramo, Network analyser studes of electromagnetc cavty resonators, Proc IRE, vol. 3, pp , [4] B. Chopard and M. Droz, Cellular Automata Modelng of Physcal Systems. Cambrdge Unversty Press, [5] Y. Kagawa and T. Yamabuch, Fnte-element equvalent crcuts for acoustc feld, Journal of the Acoustcal Socety of Amerca, vol. 64, pp , [6] Y. Kagawa, T. Tsuchya, B. Fuj, and K. Fuchoka, Dscrete Huygens model approach to sound wave propagaton, Journal of Sound and Vbraton, vol. 18, no. 3, pp , on FEM Technology 7
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 informationNumerical 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 informationNON-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 informationSIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract
SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME PACS REFERENCE: 43.20.Mv Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Zeunerstr.
More informationModule 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 informationInductance 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 informationELASTIC 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 informationAppendix 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 informationLecture 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 informationStatistical 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 informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationAmplification 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 informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationTHE 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 informationNumerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation
BULLETIN of the Malaysan Mathematcal Scences Socety http://math.usm.my/bulletn Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), 163 172 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationSUPPLEMENTARY INFORMATION
do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationDUE: 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 information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationTitle: 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 informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationCOMPOSITE 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 informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationOne-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 informationHomework 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 informationChapter - 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 informationCFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL
CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL 1. Introducton Chrstophe Vallée and Thomas Höhne In dfferent scenaros of small break Loss of Coolant Accdent (SB-LOCA), stratfed twophase
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationFrequency 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 informationComputational Electromagnetics in Antenna Analysis and Design
Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationNUMERICAL 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 informationThree-dimensional eddy current analysis by the boundary element method using vector potential
Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty
More informationFinite 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χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationIrregular 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 informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More information(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 informationIndeterminate 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 informationWeek3, 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 informationTemperature. 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 informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationThe 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 informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationMultiple Sound Source Location in 3D Space with a Synchronized Neural System
Multple Sound Source Locaton n D Space wth a Synchronzed Neural System Yum Takzawa and Atsush Fukasawa Insttute of Statstcal Mathematcs Research Organzaton of Informaton and Systems 0- Mdor-cho, Tachkawa,
More informationSupplementary 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 informationGravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)
Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationA 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 informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationComputational Modelling of the Unbalanced Magnetic Pull by Finite Element Method
Avalable onlne at www.scencedrect.com Proceda Engneerng 48 (2012 ) 83 89 MMaMS 2012 Computatonal Modellng of the Unbalanced Magnetc Pull by Fnte Element Method Martn Donát a * a Brno Unversty of Technology,
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationV.C The Niemeijer van Leeuwen Cumulant Approximation
V.C The Nemejer van Leeuwen Cumulant Approxmaton Unfortunately, the decmaton procedure cannot be performed exactly n hgher dmensons. For example, the square lattce can be dvded nto two sublattces. For
More informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationMAGNETISM MAGNETIC DIPOLES
MAGNETISM We now turn to magnetsm. Ths has actually been used for longer than electrcty. People were usng compasses to sal around the Medterranean Sea several hundred years BC. However t was not understood
More informationBoundaries, 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 informationConstitutive Modelling of Superplastic AA-5083
TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy
More informationAnalytical 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 informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationDynamics of a Superconducting Qubit Coupled to an LC Resonator
Dynamcs of a Superconductng Qubt Coupled to an LC Resonator Y Yang Abstract: We nvestgate the dynamcs of a current-based Josephson juncton quantum bt or qubt coupled to an LC resonator. The Hamltonan of
More informationElectrical 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), it produces a response (output function g (x)
Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationarxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationNumerical Simulation of Wave Propagation Using the Shallow Water Equations
umercal Smulaton of Wave Propagaton Usng the Shallow Water Equatons Junbo Par Harve udd College 6th Aprl 007 Abstract The shallow water equatons SWE were used to model water wave propagaton n one dmenson
More informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationA constant recursive convolution technique for frequency dependent scalar wave equation based FDTD algorithm
J Comput Electron (213) 12:752 756 DOI 1.17/s1825-13-479-2 A constant recursve convoluton technque for frequency dependent scalar wave equaton bed FDTD algorthm M. Burak Özakın Serkan Aksoy Publshed onlne:
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationWAVE PROPAGATION, REFLECTION AND TRANSMISSION IN CURVED BEAMS
ICSV4 Carns Australa 9- July, 7 WAVE PROPAGATION, REFECTION AND TRANSMISSION IN CURVED BEAMS Seung-Kyu ee, Bran Mace and Mchael Brennan NVH Team, R&D Centre, Hankook Tre Co., td. -, Jang-Dong, Yuseong-Gu,
More informationThe optimal delay of the second test is therefore approximately 210 hours earlier than =2.
THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple
More informationGlobal 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 informationDifference 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 informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationStudy on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI
2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationAPPLICATION OF EDDY CURRENT PRINCIPLES FOR MEASUREMENT OF TUBE CENTERLINE
APPLICATION OF EDDY CURRENT PRINCIPLES FOR MEASUREMENT OF TUBE CENTERLINE DEFLECTION E. J. Chern Martn Maretta Laboratores 1450 South Rollng Road Baltmore, MD 21227 INTRODUCTION Tubes are a vtal component
More informationStudy on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component
Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc
More informationFlow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved
Smulaton of nose generaton and propagaton caused by the turbulent flow around bluff bodes Zamotn Krll e-mal: krart@gmal.com, cq: 958886 Summary Accurate predctons of nose generaton and spread n turbulent
More information