A TIME-DOMAIN FINITE-DIFFERENCE METHOD WITH ATTENUATION BY A RECURSIVE ALGORITHM
|
|
- Franklin Hill
- 6 years ago
- Views:
Transcription
1 A TIME-DOMAIN FINITE-DIFFERENCE METHOD WITH ATTENUATION BY A RECURSIVE ALGORITHM Ningya Cheng Arthur 'C. H."Cheng and M. Nafi Tksz Earth Resurces Labratry Department f Earth Atmspheric and Planetary Sciences Massachusetts Institute f Technlgy Cambridge MA ABSTRACT A recursive algrithm t incrprate attenuatin int a time-dmain finite-difference calculatin is develped. First a rhelgical mdel f the generalized Maxwell bdy is chsen. The discrete relaxatin frequency and the peak strength f these Maxwell bdies are jintly determined by fitting t an arbitrary Q law in the frequency band f the interest. A cnjugate gradient technique and a randmly chsen starting mdel are used t determine ptimum fitting. Examples f cnstant and frequency dependent Q mdels are shwn. Secnd in rder t include the attenuatin int a finite-difference staggered-grid scheme the cnvlutin integral f stress and strain is evaluated directly. The cnvlutin integral can be expressed recursively. This is pssible because the timedmain viscelastic mdulus functin is expnential. The implementatin f 1-D wave prpagatin in a cnstant Q medium is shwn as a example. At the distance f 50 wavelengths and with three relaxatin frequencies the finite-difference results are in very gd agreement with the analytic slutins. INTRODUCTION Time-dmain finite-difference is a widely used methd in seismic wave prpagatin simulatin fr prblems with cmplex gemetry. One majr disadvantage f the time-dmain finite-difference methd is that it is hard t include attenuatin in the calculatin. Fr frequency dmain methds attenuatin can be incrprated int the calculatin simply by making the elastic mduli a cmplex functin f frequency. On the ther hand in the time-dmain stress is expressed as a cnvlutin integral f the histry f strain. 12-1
2 Cheng et al. Direct evaluatin f this cnvlutin integral requires immense strage as well as the number f peratins. Day and Minster (1984) presented a methd t incrprate attenuatin int timedmain calculatins. Their apprach is t transfrm the cnvlutin integral relating stress t strain histry int a sequence f differential peratrs. The Pade apprximatin is used in this transfrmatin. Fr a frequency independent Q the peratr cefficients can be btained in clsed frm in terms f Legendre plynmials. Emmerich and Krn (1987) imprved this apprach by using the rhelgic.al mdel fthe generalized Maxwell bdy which has a mdulus f the desired ratinal frm. This differential peratr apprach is als used in viscelastic finite-difference mdeling (Rbertssn et ai. 1994) and the pseudspectral methd (Carcine 1993). The purpse f this paper is t develp an algrithm t incrprate the attenuatin int time-dmain finite-difference calculatins. We chse the rhelgical mdel f the generalized Maxwell bdy. The discrete relaxatin frequencies and their peaks f strength are jintly determined by fitting a specific attenuatin law. The cnjugate gradient methd is used t d the ptimizatin. We_ directly evaluate the cnvlutin integral by a recursive algrithm. This is pssible because the time-dmain viscelastic mdulus functin is expnential. Finally a 1-D finite-difference wave prpagatin example is shwn. APPROXIMATION OF ATTENUATION LAW The relatinship between stress (J and strain E: be expressed in the frequency dmain as: in the case f linear viscelastisity can (J(W) = M(w)E:(w) (1) where M is the cmplex viscelastic mdulus. M beys a Kramers-Krnig relatin t make the relatinship between stress and strain causal. The quality factr Q is defined as: Q-l( ) = Im(M) W Re(M) (2) Here we apprximate the cmplex mdulus M by MN in which the relaxatin spectra cnsists f N single peaks f strength aj at discrete relaxatin frequencies Wj. This gives the cmplex mdulus as: (3) where M. is the unrelaxed mdulus. IlM = M. - Mr where Mr is the relaxed mdulus. The cmplex mdulus MN can als be expressed by using Mr. This allws us t interpret the mdel as a sum f Maxwell bdies with viscsity ajllm/wj and elastic 12-2
3 Time-Dmain Finite Difference Methd mdulus ajb.m plus an extra elastic element M r be determined by (Emmerich and Krn 1987). Q can (4) where Yj = ~Z' aj' Fr a given N the cefficients Yj and Wj have t be determined by fitting t a Q law in which we are interested. Seismic data suggest that in the seismic wave bandq-l is nearly independent f frequency. This can be characterized by a brad relaxatin spectrum (Liu et al. 1976). Fr a frequency independent Q Emmerich and Krn (1987) set the discrete relaxatin frequencies Wj equidistant n a lgarithmic scale and then determined Yj by least squares fitting. Blanch et al. (1995) apprximate Q by the summatin f Debye functins and slved the prblem analytically by the T methd but nce again the relaxatin time was pre-determined. In this paper we determine Yj and Wj by a purely numerical methd. The prblem is psed as minimizatin f expressin: with respect t Yj and Wj; [w a Wb] is the frequency band we are interested in. This methd is used fr this minimizatin purpse (Press et ai. 1992). The cnjugate gradient methd requires nt nly the value f functin but als its gradient. The gradient is cmputed by the finite-difference apprximatin. The gradient infrmatin is als updated when the mdel is updated. T avid being trapped in a lcal minimum the initial Yj is randmly chsen between zer and Q-l and Wj is randmly chsen between Wa and Wb. In the calculatin there are 1000 randm starting mdels tested and Yj and Wj are determined frm the minimum 0 f all randm starting mdels. As numerical examples we cnsider a frequency independent Q = 20 in the frequency range frm t Hz. The cmparisn f Q 0 1 and Q-;/ are shwn in Figures 1 and 2 fr the cases f N = 3 and 5 respectively. The discrete relaxatin spectra are als shwn n the plts. With N =5 Q"il prvided a very gd apprximatin t the cnstant Q mdel in the frequency band. Fr N :':: 4 the discrete relaxatin spectra are apprximately equally distanced n a lgarithmic scale. Hwver when N = 5 the spacing n lgarithmic scale is n lnger equal. The unevenly distributed relaxatin spectrum has a much smaller amplitude than an evenly distributed relaxatin spectra (Figure 2a). Next we cnsider a frequency dependent Q. It is given by Q = Q..fl. This kind f frequency dependence is due t fluid flw abve a critical frequency (Tksiiz et al. 1987). The result is pltted in Figure 3 fr N = 3. The spacing n a lgarithmic scale is nt unifrm fr discrete rela;{atin spectra distributins in the frequency dependent Q case. Even with three relaxatin spectra the apprximatin is already quite gd (5)
4 Cheng et al. FINITE-DIFFERENCE METHOD WITH ATTENUATION The wave equatin in a 1-D medium is: 8v 8u P8t=8x' where P is the density and v is the velcity. The time-dmain viscelastic mdulus MN can be btained by an '.nverse Furier transfrm f Equatin (3) MN(t)=Mu (O(t)- ~YjWje-Wj'H(t)) where H(t) is a unit step functin. In rder t use a velcity-stress frmulatin we need the time derivative f the stress. Thus we multiply by iw n bth sides f equatin (1) substituting in MN(w) and ding an inverse Furier transfrm t a time-dmain we btain: 8u 8t (6) (7) (8) and 8E: 8v = 8t 8x We let t = k6.t and assume ~~ Xj(k6.t) = xj in a discrete manner: xj = = k-l l(m+lll:>t 8 ( ) I: e-wj('-tly_~dt m=o ml:>t 8x (9) (10) t be cnstant ver each time interval 6.t. We evaluate ewjl:>t - 1 -kw.l:>t ~ mw L:>t8v(m6.t) ----e ' ~e '. Wj m=o 8x At first glance it seems that cmputing the cnvlutin summatin (Eq. 11) requires immense strage and a large number f peratins. But since it is a summatin f expnential functins it can be updated recursively (Luebbers et al. 1990). We derive the frmula t cmpute xj+l in terms f Xj. With sme algebra we btain (11) X k+l _ X k -w.l:>t (1 -e-wjl:>') 8v(k6.t). - e ' + J J w 8x J (12) 12-4
5 Time-Dmain Finite Difference Methd This is the basic equatin fr cmputing xj+1 recursively with initial cnditin XJ = O. T check its cnsistency with earlier wrk we derive the differential equatin which Xj satisfies. With straightfrward differentiatin f Equatin (9) we have X j OE: at + WjXj = at' (13) This is exactly equivalent t the first-rder differential equatin which the memry variables satisfied [Day and Minster 1984 Eq.(27)]. STAGGERED GRID ALGORITHM Equatins (6) (8) and (10) were used fr the finite-difference calculatins. Derivatives can be discretized using centered finite-difference n a staggered grid. The finitedifference scheme can be written as: and where k+l/2 _ k-i/2 6.t k k v i + I / 2 - v i +1/2 + p6.x (O'i+1 - O'i) (14) N k+l k M6.t (k+i/2 k+1/2 '" k+l O'i = O'i +~ v i +1/2 - v i _ I / 2 ) - M6.t L.. YjWjXji J:=1 -w.!!.t) k+l/2 _ k+i/2 Xk+1 = X k.e-wj!!.t e' vi+ 1/ 2 vi_ I / 2 ]t J't ( w'..6..x' J where i is the index fr the x axis discretizatin. T test ur recursive algrithm we cnsider a hmgeneus medium with velcity V p = V~" = 2000mjs density p = 1 gjcc and frequency-independent Q = 20. The discretized grid size is 6.x = 100m and 6.t = The surce signal is fed int the stress grid and a Kelly wavelet with a center frequency f 1 Hz is used. A very simple absrbing bundary cnditin is used at bth ends f the grid. The ptimal discrete relaxatin frequency and peak strength fr attenuatin are determined in the frequency bandwidth between t 10 5 Hz. The analytic slutin in the hmgeneus medium is btained by a substitute velcity via the rule: ( 1 Wi) V p = VI 1 + -QIn(-) - - 1r 21r 2Q where VI is the velcity at frequency 1 Hz. The cmparisns f the finite-difference slutin and the analytic slutin are pltted in Figure 4 fr the distance f 50 wavelengths. The finite-difference slutin with ne relaxatin frequency experienced a big errr. The slutin drifted int the lwer frequencies. With three relaxatin frequencies they prvide a fairly gd apprximatin t the cnstant Q. The finite-difference wavefrm matches the analytic slutin very well (15) (16) (17)
6 Cheng et al. CONCLUSIONS We develped a recursive algrithm t incrprate attenuatin int a time-dmain finitedifference methd. Unlike previus studies we directly evaluated the cnvlutin integral f strain and stress by a recursive algrithm. This methd is easy t understand and implement and cmputatinally efficient. In rder t apprximate the particular attenuatin law the discrete relaxatin frequency and strength are determined jintly by fitting in a least-squares sense. The system is slved by the cnjugate gradient methd with a randmized starting mdel t avid -pssible lcal minima. The numerical example shws very gd agreement between analytic slutin and the results f finite-difference with attenuatin ACKNOWLEDGMENTS One authr (N.C.) was partially supprted by Ls Alams Natinal Labratry as a Psdctral Assciate. We thank Dr. Steven Day and_ annymus reviewers fr valuable cmments n the manuscript. This wrk is supprted by the Brehle Acustics and Lgging and the Reservir Delineatin Cnsrtia at MIT. 12-6
7 Time-Dmain Finite Difference Methd REFERENCES Blanch J.O. Rbertssn J.O.A. and Symes W.W Mdeling f a cnstant Q: Methdlgy and algrithm fr an efficient and ptimally inexpensive viscelastic technique: Gephysics Carcine J.M Seismic mdeling in viscelastic media: Gephysics 58 1l0-l20. Day S.M and Minster J.B Numerical simulatin f attenuated wavefields using a Pade apprximatin methd: Gephys. l.r. Astr. Sc Emmerich H. and Krn M Incrpratin f attenuatin int time-dmain cmputatins f seismic wave field: Gephysics Liu H.P. Andersn D.L. and Kanamri H Velcity dispersin due t anelasticity: implicatins fr seismlgy and mantle cmpsitin: Gephys. l.r. astr. Sc Luebbers R. Hunsberger F.P. Kunz KS. Standler R.B. and Schneider M A frequency-dependent finite difference time dmain frmulatin fr dispersive materials: IEEE Trans. Electrmagn. Cmpat Press W.H. Teuklsky S.A. Vetterling W~T. and Flannery B.P Numerical Recipes: Cambridge University Press. Rbertssn J.O.A. Blanch J.O. and Symes W.W Viscelastic finite-difference mdeling: Gephysics 5g Tksiiz M.N. Mandai B. and Dainty A.M Frequency dependent attenuatin in the crust: Gephys. Res. Let
8 Cheng et al. 0.2 " (a) " " " Frequency (Hz) I (b)?! ' Frequency (Hz) Figure 1: Apprximatin f frequency independent Q-l with discrete relaxatin spectra with N=3. (a) Discrete relaxatin frequency and its peak strength. (b) Slid line fr exact Q-l and dashed line fr apprximated Q-l. 12-8
9 " Time-Dmain Finite Difference Methd (a). > 'I " Frequency (Hz) 10 I 100 (b) e' Frequency (Hz) Figure 2: Apprximatin f frequency independent Q-l with discrete relaxatin spectra with N=5. (a) Discrete relaxatin frequency and its peak strength. (b)slid line fr exact Q-l and dashed line fr apprximated Q-l. 12-9
10 Cheng et al f- - I 0 'r 'I 'I Freqnency (Hz) (b) CI Frequency (Hz) Figure 3: Apprximatin f frequency dependent Q-l with discrete relaxatin spectra with N=3. (a) Discrete relaxatin frequency and its peak strength. (b) Slid line fr exact Q-l and dashed line fr apprximated Q-l
11 Time-Dmain Finite Difference Methd N=l "' T XN (sec) N=2 ~ - - ~ T XN (sec) N=3 ~ T X/V (sec) 6 8 Figure 4: Cmparisn f analytic slutin with finite difference result at distance 50 wavelength fr Q=20. Discrete relaxatin spectra N = 123 are cnsidered
12 Cheng et al
ACOUSTIC LOGGING GUIDED WAVES IN TRANSVERSELY ISOTROPIC FORMATIONS
ACOUSTIC LOGGING GUIDED WAVES IN TRANSVERSELY ISOTROPIC FORMATIONS by K. J. Ellefsen, C. H. Cheng, and D.P. Schmitt" Earth Resurces Labratry Department f Earth, Atmspheric, and Planetary Sciences Massachusetts
More informationRevision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax
.7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical
More informationNumerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment
Presented at the COMSOL Cnference 2008 Hannver University f Parma Department f Industrial Engineering Numerical Simulatin f the Thermal Respsne Test Within the Cmsl Multiphysics Envirnment Authr : C. Crradi,
More informationFebruary 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA
February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal
More informationthe results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must
M.E. Aggune, M.J. Dambrg, M.A. El-Sharkawi, R.J. Marks II and L.E. Atlas, "Dynamic and static security assessment f pwer systems using artificial neural netwrks", Prceedings f the NSF Wrkshp n Applicatins
More informationDetermining the Accuracy of Modal Parameter Estimation Methods
Determining the Accuracy f Mdal Parameter Estimatin Methds by Michael Lee Ph.D., P.E. & Mar Richardsn Ph.D. Structural Measurement Systems Milpitas, CA Abstract The mst cmmn type f mdal testing system
More informationTime-domain lifted wavelet collocation method for modeling nonlinear wave propagation
Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August Time-dmain lifted wavelet cllcatin methd fr mdeling nnlinear wave prpagatin Kelvin Chee-Mun Lee and Wn-Seng Gan Digital Signal
More informationReview Problems 3. Four FIR Filter Types
Review Prblems 3 Fur FIR Filter Types Fur types f FIR linear phase digital filters have cefficients h(n fr 0 n M. They are defined as fllws: Type I: h(n = h(m-n and M even. Type II: h(n = h(m-n and M dd.
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationLead/Lag Compensator Frequency Domain Properties and Design Methods
Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin
More informationFall 2013 Physics 172 Recitation 3 Momentum and Springs
Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.
More informationDataflow Analysis and Abstract Interpretation
Dataflw Analysis and Abstract Interpretatin Cmputer Science and Artificial Intelligence Labratry MIT Nvember 9, 2015 Recap Last time we develped frm first principles an algrithm t derive invariants. Key
More informationBootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >
Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);
More informationOF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
More information7 TH GRADE MATH STANDARDS
ALGEBRA STANDARDS Gal 1: Students will use the language f algebra t explre, describe, represent, and analyze number expressins and relatins 7 TH GRADE MATH STANDARDS 7.M.1.1: (Cmprehensin) Select, use,
More informationEDA Engineering Design & Analysis Ltd
EDA Engineering Design & Analysis Ltd THE FINITE ELEMENT METHOD A shrt tutrial giving an verview f the histry, thery and applicatin f the finite element methd. Intrductin Value f FEM Applicatins Elements
More informationELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322
ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA by J. C. SPROTT December 4, 1969 PLP N. 3 These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated. They are fr private
More informationA Scalable Recurrent Neural Network Framework for Model-free
A Scalable Recurrent Neural Netwrk Framewrk fr Mdel-free POMDPs April 3, 2007 Zhenzhen Liu, Itamar Elhanany Machine Intelligence Lab Department f Electrical and Cmputer Engineering The University f Tennessee
More informationEvaluation of Eigenmodes of Dielectric Waveguides by a Numerical Method Based on the BPM
Memirs f the Faculty f EngineeringOkayama UniversityV.29 N.2 pp.49-55 March 995 Evaluatin f Eigenmdes f Dielectric Waveguides by a Numerical Methd Based n the BPM Anis AHMEDt Shinji TANIGUCHI+ Osami WADAt
More informationinitially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur
Cdewrd Distributin fr Frequency Sensitive Cmpetitive Learning with One Dimensinal Input Data Aristides S. Galanpuls and Stanley C. Ahalt Department f Electrical Engineering The Ohi State University Abstract
More informationTHE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS
THE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS M.S. Chng Department f Mechanical Engineering The University f Melburne Victria 3010 AUSTRALIA min@unimelb.edu.au J.P. Mnty
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More informationENSC Discrete Time Systems. Project Outline. Semester
ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding
More informationYou need to be able to define the following terms and answer basic questions about them:
CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f
More informationMain Menu. SEG Houston 2009 International Exposition and Annual Meeting. Summary
CO elcity Measurements and Mdels r Temperatures dwn t -10 C and up t 00 C and Pressures up t 100 MPa Min Sun* and De-hua Han, Rck Physics Lab, University Hustn Micheal Batzle, Clrad Schl Mines Summary
More informationA PLETHORA OF MULTI-PULSED SOLUTIONS FOR A BOUSSINESQ SYSTEM. Department of Mathematics, Penn State University University Park, PA16802, USA.
A PLETHORA OF MULTI-PULSED SOLUTIONS FOR A BOUSSINESQ SYSTEM MIN CHEN Department f Mathematics, Penn State University University Park, PA68, USA. Abstract. This paper studies traveling-wave slutins f the
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCurseWare http://cw.mit.edu 2.161 Signal Prcessing: Cntinuus and Discrete Fall 2008 Fr infrmatin abut citing these materials r ur Terms f Use, visit: http://cw.mit.edu/terms. Massachusetts Institute
More informationAn RWG Basis Function Based Near- to Far-Field Transformation for 2D Triangular-Grid Finite Difference Time Domain Method
An RWG Basis Functin Based Near- t Far-Field Transfrmatin fr D Triangular-Grid Finite Difference Time Dmain Methd Xuan Hui Wu, Student Member, IEEE, Ahmed A. Kishk, Fellw, IEEE, and Allen W. Glissn, Fellw,
More informationSections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.
Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage
More informationPerformance Bounds for Detect and Avoid Signal Sensing
Perfrmance unds fr Detect and Avid Signal Sensing Sam Reisenfeld Real-ime Infrmatin etwrks, University f echnlgy, Sydney, radway, SW 007, Australia samr@uts.edu.au Abstract Detect and Avid (DAA) is a Cgnitive
More informationInvestigation of a Single-Point Nonlinearity Indicator in One-Dimensional Propagation. 2 Theory
Investigatin f a Single-Pint Nnlinearity Indicatr in One-Dimensinal Prpagatin Lauren Falc, Kent Gee, Anthny Atchley, Victr Sparrw The Pennsylvania State University, Graduate Prgram in Acustics, University
More informationMore Tutorial at
Answer each questin in the space prvided; use back f page if extra space is needed. Answer questins s the grader can READILY understand yur wrk; nly wrk n the exam sheet will be cnsidered. Write answers,
More informationCHARACTERIZATION OF A PERIODIC SURFACE PROFILE BY POLE-ZERO PARAMETERIZATION OF ELASTODYNAMIC PULSE REFLECTIONS * R.A. Roberts and J.D.
CHARACTERZATON OF A PERODC SURFACE PROFLE BY POLE-ZERO PARAMETERZATON OF ELASTODYNAMC PULSE REFLECTONS * R.A. Rberts and J.D. Achenbach The Technlgical nstitute Nrthwestern University Evanstn, L. 60201
More informationRay tracing equations in transversely isotropic media Cosmin Macesanu and Faruq Akbar, Seimax Technologies, Inc.
Ray tracing equatins in transversely istrpic media Csmin Macesanu and Faruq Akbar, Seimax Technlgies, Inc. SUMMARY We discuss a simple, cmpact apprach t deriving ray tracing equatins in transversely istrpic
More information[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )
(Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well
More informationModeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function
www.ccsenet.rg/mer Mechanical Engineering Research Vl. 1, N. 1; December 011 Mdeling the Nnlinear Rhelgical Behavir f Materials with a Hyper-Expnential Type Functin Marc Delphin Mnsia Département de Physique,
More informationEngineering Approach to Modelling Metal THz Structures
Terahertz Science and Technlgy, ISSN 1941-7411 Vl.4, N.1, March 11 Invited Paper ngineering Apprach t Mdelling Metal THz Structures Stepan Lucyszyn * and Yun Zhu Department f, Imperial Cllege Lndn, xhibitin
More informationLeast Squares Optimal Filtering with Multirate Observations
Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical
More informationDESIGN OPTIMIZATION OF HIGH-LIFT CONFIGURATIONS USING A VISCOUS ADJOINT-BASED METHOD
DESIGN OPTIMIZATION OF HIGH-LIFT CONFIGURATIONS USING A VISCOUS ADJOINT-BASED METHOD Sangh Kim Stanfrd University Juan J. Alns Stanfrd University Antny Jamesn Stanfrd University 40th AIAA Aerspace Sciences
More informationSAMPLING DYNAMICAL SYSTEMS
SAMPLING DYNAMICAL SYSTEMS Melvin J. Hinich Applied Research Labratries The University f Texas at Austin Austin, TX 78713-8029, USA (512) 835-3278 (Vice) 835-3259 (Fax) hinich@mail.la.utexas.edu ABSTRACT
More informationMATHEMATICS SYLLABUS SECONDARY 5th YEAR
Eurpean Schls Office f the Secretary-General Pedaggical Develpment Unit Ref. : 011-01-D-8-en- Orig. : EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 6 perid/week curse APPROVED BY THE JOINT TEACHING COMMITTEE
More informationFIELD QUALITY IN ACCELERATOR MAGNETS
FIELD QUALITY IN ACCELERATOR MAGNETS S. Russenschuck CERN, 1211 Geneva 23, Switzerland Abstract The field quality in the supercnducting magnets is expressed in terms f the cefficients f the Furier series
More informationDead-beat controller design
J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable
More informationVerification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance
Verificatin f Quality Parameters f a Slar Panel and Mdificatin in Frmulae f its Series Resistance Sanika Gawhane Pune-411037-India Onkar Hule Pune-411037- India Chinmy Kulkarni Pune-411037-India Ojas Pandav
More informationEmphases in Common Core Standards for Mathematical Content Kindergarten High School
Emphases in Cmmn Cre Standards fr Mathematical Cntent Kindergarten High Schl Cntent Emphases by Cluster March 12, 2012 Describes cntent emphases in the standards at the cluster level fr each grade. These
More informationCAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank
CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal
More information7.0 Heat Transfer in an External Laminar Boundary Layer
7.0 Heat ransfer in an Eternal Laminar Bundary Layer 7. Intrductin In this chapter, we will assume: ) hat the fluid prperties are cnstant and unaffected by temperature variatins. ) he thermal & mmentum
More informationSection I5: Feedback in Operational Amplifiers
Sectin I5: eedback in Operatinal mplifiers s discussed earlier, practical p-amps hae a high gain under dc (zer frequency) cnditins and the gain decreases as frequency increases. This frequency dependence
More informationTheoretical study of third virial coefficient with Kihara potential
Theretical study f third virial cefficient with Kihara ptential Jurnal: Manuscript ID cjp-017-0705.r Manuscript Type: Article Date Submitted by the Authr: 6-Dec-017 Cmplete List f Authrs: Smuncu E.; Giresun
More informationand the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:
Algrithm fr Estimating R and R - (David Sandwell, SIO, August 4, 2006) Azimith cmpressin invlves the alignment f successive eches t be fcused n a pint target Let s be the slw time alng the satellite track
More informationROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS
Particle Acceleratrs, 1986, Vl. 19, pp. 99-105 0031-2460/86/1904-0099/$15.00/0 1986 Grdn and Breach, Science Publishers, S.A. Printed in the United States f America ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS
More informationA spectral scheme for wave propagation simulation in 3-D elastic-anisotropic media
GEOPHYSICS. VOL. 57, NO. 12 (DECEMBER 1992); P. 1593-1607, 13 FIGS. A spectral scheme fr wave prpagatin simulatin in 3-D elastic-anistrpic media Jse M. Carcine*, Dan Kslftt, Alfred Behle**, and Geza Serianis
More informationExaminer: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data
Benha University Cllege f Engineering at Banha Department f Mechanical Eng. First Year Mechanical Subject : Fluid Mechanics M111 Date:4/5/016 Questins Fr Final Crrective Examinatin Examiner: Dr. Mhamed
More informationResampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017
Resampling Methds Crss-validatin, Btstrapping Marek Petrik 2/21/2017 Sme f the figures in this presentatin are taken frm An Intrductin t Statistical Learning, with applicatins in R (Springer, 2013) with
More informationChapter 3: Cluster Analysis
Chapter 3: Cluster Analysis } 3.1 Basic Cncepts f Clustering 3.1.1 Cluster Analysis 3.1. Clustering Categries } 3. Partitining Methds 3..1 The principle 3.. K-Means Methd 3..3 K-Medids Methd 3..4 CLARA
More informationProgress In Electromagnetics Research M, Vol. 9, 9 20, 2009
Prgress In Electrmagnetics Research M, Vl. 9, 9 20, 2009 WIDE-ANGLE REFLECTION WAVE POLARIZERS USING INHOMOGENEOUS PLANAR LAYERS M. Khalaj-Amirhsseini and S. M. J. Razavi Cllege f Electrical Engineering
More informationAdministrativia. Assignment 1 due thursday 9/23/2004 BEFORE midnight. Midterm exam 10/07/2003 in class. CS 460, Sessions 8-9 1
Administrativia Assignment 1 due thursday 9/23/2004 BEFORE midnight Midterm eam 10/07/2003 in class CS 460, Sessins 8-9 1 Last time: search strategies Uninfrmed: Use nly infrmatin available in the prblem
More informationCalculating the optimum pressure and temperature for vacancy minimization from theory; Niobium is an example. Jozsef Garai
Calculating the ptimum pressure and temperature fr vacancy minimizatin frm thery; Nibium is an example Jzsef Garai Department f Mechanical and Materials Engineering, Flrida Internatinal University, Miami,
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationChapter 3 Kinematics in Two Dimensions; Vectors
Chapter 3 Kinematics in Tw Dimensins; Vectrs Vectrs and Scalars Additin f Vectrs Graphical Methds (One and Tw- Dimensin) Multiplicatin f a Vectr b a Scalar Subtractin f Vectrs Graphical Methds Adding Vectrs
More informationChapter 4. Unsteady State Conduction
Chapter 4 Unsteady State Cnductin Chapter 5 Steady State Cnductin Chee 318 1 4-1 Intrductin ransient Cnductin Many heat transfer prblems are time dependent Changes in perating cnditins in a system cause
More informationSOLUTION OF THREE-CONSTRAINT ENTROPY-BASED VELOCITY DISTRIBUTION
SOLUTION OF THREECONSTRAINT ENTROPYBASED VELOCITY DISTRIBUTION By D. E. Barbe,' J. F. Cruise, 2 and V. P. Singh, 3 Members, ASCE ABSTRACT: A twdimensinal velcity prfile based upn the principle f maximum
More informationCambridge Assessment International Education Cambridge Ordinary Level. Published
Cambridge Assessment Internatinal Educatin Cambridge Ordinary Level ADDITIONAL MATHEMATICS 4037/1 Paper 1 Octber/Nvember 017 MARK SCHEME Maximum Mark: 80 Published This mark scheme is published as an aid
More informationA Matrix Representation of Panel Data
web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins
More information^YawataR&D Laboratory, Nippon Steel Corporation, Tobata, Kitakyushu, Japan
Detectin f fatigue crack initiatin frm a ntch under a randm lad C. Makabe," S. Nishida^C. Urashima,' H. Kaneshir* "Department f Mechanical Systems Engineering, University f the Ryukyus, Nishihara, kinawa,
More informationLecture 17: Free Energy of Multi-phase Solutions at Equilibrium
Lecture 17: 11.07.05 Free Energy f Multi-phase Slutins at Equilibrium Tday: LAST TIME...2 FREE ENERGY DIAGRAMS OF MULTI-PHASE SOLUTIONS 1...3 The cmmn tangent cnstructin and the lever rule...3 Practical
More informationMath Foundations 20 Work Plan
Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant
More informationCalculation of Ramp Response of Lossy Transmission Lines Using Two-port Network Functions
Calculatin f amp espnse f Lssy Transmissin Lines Using Tw-prt Netwrk Functins Payam Heydari and Massud Pedram Dept. f EE-Systems University f Suthern Califrnia Ls Angeles, CA 989 E-mail: {payam,massud}@zugrs.usc.edu
More informationANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES
ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES CONSISTING OF LONG PARALLEL CONDUCTORS B. de Halleux, O. Lesage, C. Mertes and A. Ptchelintsev Mechanical Engineering Department Cathlic University
More informationk-nearest Neighbor How to choose k Average of k points more reliable when: Large k: noise in attributes +o o noise in class labels
Mtivating Example Memry-Based Learning Instance-Based Learning K-earest eighbr Inductive Assumptin Similar inputs map t similar utputs If nt true => learning is impssible If true => learning reduces t
More informationA Simple Set of Test Matrices for Eigenvalue Programs*
Simple Set f Test Matrices fr Eigenvalue Prgrams* By C. W. Gear** bstract. Sets f simple matrices f rder N are given, tgether with all f their eigenvalues and right eigenvectrs, and simple rules fr generating
More informationarxiv:hep-ph/ v1 2 Jun 1995
WIS-95//May-PH The rati F n /F p frm the analysis f data using a new scaling variable S. A. Gurvitz arxiv:hep-ph/95063v1 Jun 1995 Department f Particle Physics, Weizmann Institute f Science, Rehvt 76100,
More informationIntroduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem
A Generalized apprach fr cmputing the trajectries assciated with the Newtnian N Bdy Prblem AbuBar Mehmd, Syed Umer Abbas Shah and Ghulam Shabbir Faculty f Engineering Sciences, GIK Institute f Engineering
More information( ) (1) ρ c crustal density 2600 kg m -3 ρ w water density 1000 kg m -3. HEAT FLOW PARADOX (Copyright 2001, David T. Sandwell)
1 HEAT FLOW PARADOX (Cpyright 2001, David T. Sandwell) (See Special Issue f J. Gephys. Res., v.85, 1980: A) Turctte, Tag, and Cper, A Steady-State mdel fr the distributin f stress and temperature n the
More informationOTHER USES OF THE ICRH COUPL ING CO IL. November 1975
OTHER USES OF THE ICRH COUPL ING CO IL J. C. Sprtt Nvember 1975 -I,," PLP 663 Plasma Studies University f Wiscnsin These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated.
More informationComputational Methods CMSC/AMSC/MAPL 460. Ramani Duraiswami, Dept. of Computer Science
Cmputatinal Methds CMSC/AMSC/MAPL 460 Ramani Duraiswami, Dept. f Cmputer Science Curse Gals Intrductin t the use f scientific cmputing techniques t slve prblems in varius dmains Understand principles behind
More informationFast Fourier Transform-based micromechanical modeling of polycrystals with direct input from microstructural images
Fast Furier Transfrm-based micrmechanical mdeling f plycrystals with direct input frm micrstructural images Ricard A. Lebenshn (Ls Alams Natinal Labratry, USA) FFT-based frmulatin fr plycrystals: fast
More informationApplication of APW Pseudopotential Form Factor in the Calculation of Liquid Metal Resistivities.
Internatinal Jurnal f Pure and Applied Physics. ISSN 097-1776 Vlume 8, Number (01), pp. 11-117 Research India Publicatins http://www.ripublicatin.cm/pap.htm Applicatin f APW Pseudptential Frm Factr in
More informationKeysight Technologies Understanding the Kramers-Kronig Relation Using A Pictorial Proof
Keysight Technlgies Understanding the Kramers-Krnig Relatin Using A Pictrial Prf By Clin Warwick, Signal Integrity Prduct Manager, Keysight EEsf EDA White Paper Intrductin In principle, applicatin f the
More informationA. HAUGUEL (l), G. LABADIE (2) and B. LATTEUX (3)
A FINITE ELEMENT METHOD FOR THE SHALLOW WATER EQUATIONS by A. HAUGUEL (l), G. LABADIE (2) and B. LATTEUX (3) ABSTRACT The paper reprts the current prgress in develping a finite element methd fr the shallw
More informationChebyshev Pseudospeetral Method of Viscous Flows with Corner Singularities
Jurnal f Scientific Cmputing, Vl. 4, N. 1, 1989 Chebyshev Pseudspeetral Methd f Viscus Flws with Crner Singularities W. W. Schaltz, 1 N. Y. Lee, 1'2 and J. P. Byd 3 Received February 3, 1989 Chebyshev
More informationThermodynamics and Equilibrium
Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,
More informationECE 2100 Circuit Analysis
ECE 00 Circuit Analysis Lessn 6 Chapter 4 Sec 4., 4.5, 4.7 Series LC Circuit C Lw Pass Filter Daniel M. Litynski, Ph.D. http://hmepages.wmich.edu/~dlitynsk/ ECE 00 Circuit Analysis Lessn 5 Chapter 9 &
More informationA mathematical model for complete stress-strain curve prediction of permeable concrete
A mathematical mdel fr cmplete stress-strain curve predictin f permeable cncrete M. K. Hussin Y. Zhuge F. Bullen W. P. Lkuge Faculty f Engineering and Surveying, University f Suthern Queensland, Twmba,
More informationIncreasing Heat Transfer in Microchannels with Surface Acoustic Waves*
Increasing Heat Transfer in Micrchannels with Surface Acustic Waves* Shaun Berry 0/9/04 *This wrk was spnsred by the Department f the Air Frce under Air Frce Cntract #FA87-05-C-000. Opinins, interpretatins,
More informationPressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects
Pressure And Entrpy Variatins Acrss The Weak Shck Wave Due T Viscsity Effects OSTAFA A. A. AHOUD Department f athematics Faculty f Science Benha University 13518 Benha EGYPT Abstract:-The nnlinear differential
More informationSlide04 (supplemental) Haykin Chapter 4 (both 2nd and 3rd ed): Multi-Layer Perceptrons
Slide04 supplemental) Haykin Chapter 4 bth 2nd and 3rd ed): Multi-Layer Perceptrns CPSC 636-600 Instructr: Ynsuck Che Heuristic fr Making Backprp Perfrm Better 1. Sequential vs. batch update: fr large
More informationPattern Recognition 2014 Support Vector Machines
Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft
More informationON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT
17 Kragujevac J. Sci. 8 (006) 17-4. ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem Ali Attia
More informationInflow Control on Expressway Considering Traffic Equilibria
Memirs f the Schl f Engineering, Okayama University Vl. 20, N.2, February 1986 Inflw Cntrl n Expressway Cnsidering Traffic Equilibria Hirshi INOUYE* (Received February 14, 1986) SYNOPSIS When expressway
More informationChecking the resolved resonance region in EXFOR database
Checking the reslved resnance regin in EXFOR database Gttfried Bertn Sciété de Calcul Mathématique (SCM) Oscar Cabells OECD/NEA Data Bank JEFF Meetings - Sessin JEFF Experiments Nvember 0-4, 017 Bulgne-Billancurt,
More informationA Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture
Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu
More informationSupport-Vector Machines
Supprt-Vectr Machines Intrductin Supprt vectr machine is a linear machine with sme very nice prperties. Haykin chapter 6. See Alpaydin chapter 13 fr similar cntent. Nte: Part f this lecture drew material
More informationSIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST. Mark C. Otto Statistics Research Division, Bureau of the Census Washington, D.C , U.S.A.
SIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST Mark C. Ott Statistics Research Divisin, Bureau f the Census Washingtn, D.C. 20233, U.S.A. and Kenneth H. Pllck Department f Statistics, Nrth Carlina State
More informationMethods for Determination of Mean Speckle Size in Simulated Speckle Pattern
0.478/msr-04-004 MEASUREMENT SCENCE REVEW, Vlume 4, N. 3, 04 Methds fr Determinatin f Mean Speckle Size in Simulated Speckle Pattern. Hamarvá, P. Šmíd, P. Hrváth, M. Hrabvský nstitute f Physics f the Academy
More informationLearning to Control an Unstable System with Forward Modeling
324 Jrdan and Jacbs Learning t Cntrl an Unstable System with Frward Mdeling Michael I. Jrdan Brain and Cgnitive Sciences MIT Cambridge, MA 02139 Rbert A. Jacbs Cmputer and Infrmatin Sciences University
More informationMath 302 Learning Objectives
Multivariable Calculus (Part I) 13.1 Vectrs in Three-Dimensinal Space Math 302 Learning Objectives Plt pints in three-dimensinal space. Find the distance between tw pints in three-dimensinal space. Write
More informationTHE PARTITION OF ENERGY INTO WAVES AND CURRENTS
THE PARTITION OF ENERGY INTO WAVES AND CURRENTS W. Perrie, C. Tang, Y. Hu and B.M. DeTracy Fisheries & Oceans Canada, Bedfrd Institute f Oceangraphy, Dartmuth, Nva Sctia, Canada 1. INTRODUCTION Ocean mdels
More informationMaterials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion
Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals f Diffusin Diffusin: Transprt in a slid, liquid, r gas driven by a cncentratin gradient (r, in the case f mass transprt, a chemical ptential
More information