Time-domain lifted wavelet collocation method for modeling nonlinear wave propagation
|
|
- Cathleen Cunningham
- 6 years ago
- Views:
Transcription
1 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 Prcessing Labratry, Schl f Electrical and Electrnic Engineering, Nanyang Technlgical University, S-Ba-, Nanyang Ave Singapre 998 ecmlee@ntu.edu.sg and ewsgan@ntu.edu.sg Abstract: A time-dmain adaptive numerical methd fr mdeling nnlinear wave prpagatin is develped. This methd is based n a secnd-generatin wavelet cllcatin using a lifting scheme and makes use f the multilevel decmpsitin nature f the scheme t allw fr autmatic grid refinement accrding t the magnitude f wavefrm steepening. The multiplicatin in the nnlinear term is als easy due t the cllcatin nature. With threshlding, the slutin is cmpact at every level f reslutin and cmputed nly at cllcatin pints assciated with the remaining significant wavelet cefficients. The errr tlerance and cmpressin rati f the new methd are ttally cntrlled by the threshld value used. This brings substantial savings in cmputatin time when cmpared t the cnventinal finite difference scheme n a unifrmly fine grid. Acustical Sciety f America PACS numbers:..c,..l Date Received: March Date Accepted: April. Intrductin The time-dmain numerical methd f simulating nnlinear wave prpagatin in a sund beam by using an implicit backward finite difference (IBFD) scheme t slve the Khkhlv- Zabaltskaya-Kuznetsv (KZK) equatin has been described in previus literature. The accuracy f the IBFD scheme depends n the implicit integratin slver and the fineness f the cmputatinal grid used and is therefre prprtinal t the cmputatin time required t slve the prblem. The purpse f this paper is t intrduce a wavelet cllcatin apprach that slves in a cmpact manner withut the redundant cmputatins in the wellapprximated slutin regins. We use an adaptive cmputatinal grid that is able t accmmdate lcal structures thrugh dynamic grid refinement. We first study the new methd in a -D Burgers prblem and then apply it in an axisymmetric -D frm f the KZK equatin. A cmparisn between the new results and that frm a cnventinal finite difference scheme is als made t access the perfrmances f bth.. Lifted Interplating Wavelet Transfrm Wavelets frm a versatile tl fr representing general functins r data sets because they are highly capable f capturing the essence f a functin with nly a small set f cefficients. Generally, wavelets are defined as the dyadic translates and dilates f a particular mther wavelet functin. Hwever, in this sectin, a new class f wavelets termed secndgeneratin wavelets is intrduced. These birthgnal wavelets are nt necessarily translates and dilates f each ther but rather cnstructed frm a spatial dmain apprach knwn as the lifting scheme,. They allw interplatin prcedures t be defined n mre cmplicated nnunifrm cmputatinal grid structures since the filter cefficients used in these wavelets ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America
2 Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August can be psitin- and level-dependent. The simplest frm f such wavelets is the lifted interplating wavelet transfrm and is the basis set fr the wavelet cllcatin methd because f its simplicity and its in-place cmputatin withut the need f auxiliary memry. Essentially, the transfrm is cnstructed ut f the spatial dmain using basic steps: split/merge, predict (dual lifting) and update (regular lifting). The first step splits the data int subsets (dd and even); the secnd step predicts the even set based n the dd ne and finds the difference between them; and the third step updates the dd set t preserve the mean f the functin. The resulting apprximatin and detail values (r λ andγ cefficients respectively) are then btained, and the prcess is iterated using the apprximatin cefficient λ as the new input until a specified level f carsest reslutin is reached. Fig.. The frward lifted interplating wavelet transfrm (left) and its inverse (right). The multireslutin analysis prperty f wavelets is still valid fr secnd-generatin J- j j wavelets and can be expressed as V W W... W L, wherev andw are the rthgnal subspaces that crrespnd t the scaling and wavelet functins respectively, and L is the space spanned by the riginal functin. Thus, when wavelet decmpsitin is perfrmed n the riginal functin P, its representatin at the finest reslutin level J can be btained as J J j j ( ) = λφ k k ( ) + γlφl ( ). () P x x x k K j = j l L The strength f this multilevel decmpsitin is then bvius: the functin is decmpsed int a set f rthgnal cmpnents at different reslutin levels. It is crucial fr a fast way t cmpute the slutin via the cnstructin f an adaptive cmputatinal grid, which will be elabrated in the next sectin. Lcal structures f the riginal functin due t wavefrm steepening at different scales can be islated and recrded in the magnitude f theγ cefficients. We use the lifted interplating wavelet transfrm with the number f dual vanishing mments N = and the number f real vanishing mments N = thrughut ur discussin in this paper.. Numerical Methd A crdinate transfrmatin is applied t the KZK equatin fr an axisymmetric and unfcused sund beam that prpagates in the psitive z directin (see Fig. ) t arrive at the fllwing dimensinless frm in terms f the nrmalized acustic pressure P : P = P + P d ' + A P + NP P ( + σ) ρ ( ). () σ ρ ρ + σ Fr -D nnlinear wave prpagatin, a dimensinless Burgers equatin identical t Eq. () withut the diffractin term is used. The prblem is cnveniently transfrmed int a regular and dimensinless cmputatinal grid represented by the axial distanceσ and the radial distance ρ and is discretized fr finite difference calculatins as Pijk,, = P( σk, min+ i, j ρ). As we have pint values, the derivative peratrs can be implemented using cnventinal finite difference stencils whereas the integral peratr in the diffractin term ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America
3 Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August can be implemented using an adaptive Simpsn quadrature w.r.t. '. Multiplicatin in the nnlinear term is als straightfrward and a lcal peratin since we are wrking with cllcatin pints. A th rder explicit Runge-Kutta slver is als used t march frward in space, integrating Eq. () numerically with respect tσ t get the slutin at the next prpagatin step. Fig.. Radiatin gemetry frm an axisymmetric sund surce lcated in the plane z =.. Adaptive Cmputatinal Grid The multilevel decmpsitin f Eq. () allws us t cnstruct a set f nested grids with different reslutins. Assuming the slutin is nrmalized, a relative errr threshld ε can be set t split the slutin P J int regins: ne with cefficients abve the threshld and ne with cefficients belw. This can be expressed as ( σ,, ρ ) ( σ, j, ρ ) j ( σ,, ρ ) P = P + P j =,... J, J J J l l γ ε γ ε l L j. () Threshlding the latter remves a substantial amunt f cllcatin pints frm the grid because they cntain trivial details and results in lesser subsequent cmputatins. This prvides a simple way f cntrlling the cmpressin factr f the cmputatinal grid and the apprximatin errr f the slutin. The selectin f the threshld value is easy by specifying the percentage amunt f signal energy t be retained after threshlding. The finite difference and quadrature schemes are then applied lcally n each level where there are nγ cefficients n finer scales. Cnvlutin with the finite difference stencil als leads t an increase in supprt, s we need t add extra neighbring cllcatin pints next t the edge pints in the cmputatinal grid at every level s as t crrect fr this. Furthermre, the differentiated functin might als have new significant cefficients appearing in the neighbring regins. This manner f dynamically allcating cllcatin pints at every level allws the cmputatinal grid t adapt dynamically t the lcal changes as harmnics are intrduced int the slutin due t the nnlinear term. Upn inverse wavelet transfrm, the full slutin at all cllcatin pints can be btained fr the current integratin step.. Results and Discussin The results fr the prpsed numerical methd are presented. As mentined earlier, a th rder Runge-Kutta integratin slver (MATLAB s de) is used tgether with a relative errr threshld f ε = fr the first example using Burgers equatin and a relative errr threshld f ε = fr the next example using KZK equatin. Fr stability in the integratin, we use = ρ =. and an initial prpagatin step-size f σ = e, which will then be adjusted by de s variable step-size cntrl algrithm during the simulatin. Fr higher efficiency, the prgram is cmpiled as a MEX (MATLAB-EXecutable) file befre executin in MATLAB.. -D Dimensinless Burgers Equatin ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America
4 Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August The dimensinless Burgers equatin fr simulating -D nnlinear wave prpagatin is generally expressed as P P P = A NP () σ where, fr this example, A =., N =, and the initial functin atσ = is represented by P( σ =, ) = P sin( π) given.. is the range f the time windw used and the bundary cnditins are P( σ, =±.) =. The number f decmpsitin levels is fixed at and, because the functin is largely smth, we chse a large errr threshld f ε =. The slutin starts develping a sharp gradient at = and evlves int a sawtth-like wavefrm as it prpagates away frm the surce. This is a gd test f the new methd s adaptivity t refine the slutin at lcal regins where shcks are frming, i.e., t refine the cmputatinal grid arund = as the shck develps. Nte that increasingly higher reslutin is allcated t such lcal structures in the slutin. The new methd is able t refine the grid in the vicinity f the shck accrding t the magnitude f the gradient, even after much attenuatin. The number f significant cefficients in the higher levels als increases as the gradient increases until the shck is finally reslved at the finest reslutin level. The results btained are cmpared t results frm the cnventinal finite difference scheme (i.e., ε = ). The methd demnstrates sufficient rbustness in handling nnlinearities thrughut the whle cmputatin prcess. Snapsht samples f these results, taken at selected intervals σ = (,.,.,, ) are shwn belw..8.. x x x x x Fig.. Adaptivity f the lifted wavelet cllcatin methd in handling evlving shcks: refinement n a finer grid is dne as the shck gradient increases. Cmparisn with the cnventinal finite difference slutin f the Burgers equatin shws similar results within the tlerance f the given relative errr threshld ε =. ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America
5 Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August. Dimensinless Axisymmetric KZK Equatin The initial surce cnditin fr the KZK equatin in Eq. () is given as + ρ P( σ =,, ρ) = exp sin( π ( + ρ )) () π / with as the range f the time windw used and ρ as the radial distance range (Fig. ). This can als be seen as the size f the initial cmputatinal grid t be used. Initial Surce Cnditin (nrmalized).8.. Retarded Time Radial Distance ρ Fig.. Initial surce cnditin f cmputatinal grid. With the number f decmpsitin levels fixed at, the results btained frm the lifted wavelet cllcatin methd are shwn and cmpared with thse btained frm the cnventinal finite difference scheme in Figs. and. The first set f results in Fig. shws the n-axis slutins f the prpagating wave at varius prpagatin distancesσ whereas the secnd in Fig. shws the slutins alng the ρ axis at = at the same distances. The algrithm is stable with ε = thrughut the integratin range and is seen by the cnvergence f the abslute apprximatin errrs in the far-right clumns f Figs. and t within the specified ε. Hwever, it may be pssible t cntrl the errr threshld adaptively accrding t the percentage amunt f signal energy retained after threshlding s as t further reduce the cmplexity f the cmputatins. σ=, ρ=. σ=., ρ= x -. σ=, ρ= x Fig.. On-axis slutins at varius axial distances using the wavelet cllcatin methd ( N =, A =.). 8 ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America 8
6 Lee et al.: Acustics Research Letters Online [DOI./.] Published Online 8 August σ=, = σ=, = σ=., = x ρ ρ x ρ x Cnclusin Fig.. Slutins alng ρ axis at = at varius axial distances using the wavelet cllcatin methd ( N =, A =.) An adaptive numerical methd based n secnd-generatin wavelets is develped. The adaptivity f the methd is achieved by using a lifted interplating wavelet basis t apprximate the slutin and refining the cmputatinal grid in lcal regins t accmmdate develping shcks. By specifying the errr threshld value, it is pssible t actively cntrl the accuracy f the slutin and the size f the adaptive cmputatinal grid. The sparse representatin leads t significant savings in the memry strage and the number f cmputatins required t slve the Burgers and KZK equatins numerically. This methd uses an adaptive cmputatinal grid and utperfrms the cnventinal finite difference methd, which uses a unifrm fine grid thrughut. Due t the cllcatin nature, the methd can als handle general bundary cnditins and nnlinearities effrtlessly. References and links Daubechies, I., Sweldens, W., Factring Wavelet Transfrms int Lifting Steps, J. Furier Anal. Appl., (), - (998). Hlmstrm, M., Slving hyperblic PDEs using interplating wavelets, SIAM J. Sci. Cmput., (), - (999). Lee, Y. S., Hamiltn, M. F., Time Dmain Mdeling f pulsed finite-amplitude sund beams, J. Acus. Sc. Am, 9(), 9-9 (99). Sweldens, W., The lifting scheme: A cnstructin f secnd-generatin wavelets, SIAM J. Math. Anal., 9(), - (998). Vasilyev, O. V., Bwman, C., Secnd generatin wavelet cllcatin methd fr the slutin f partial differential equatins, J. Cmp. Phys,, -9 (). 9 ARLO (), Oct 9-8//()///$9. (c) Acustical Sciety f America 9
Building 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 information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationAPPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL
JP2.11 APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL Xingang Fan * and Jeffrey S. Tilley University f Alaska Fairbanks, Fairbanks,
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 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 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 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 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 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 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 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 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 informationISSUES RELATED TO THE NUMERICAL IMPLEMENTATION OF A SPARSE METHOD FOR THE SOLUTION OF VOLUME INTEGRAL EQUATIONS AT LOW FREQUENCIES
University f Kentucky UKnwledge University f Kentucky Master's Theses Graduate Schl 2010 ISSUES RELATED TO THE NUMERICAL IMPLEMENTATION OF A SPARSE METHOD FOR THE SOLUTION OF VOLUME INTEGRAL EQUATIONS
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 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 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 informationEnhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme
Enhancing Perfrmance f / Neural Classifiers via an Multivariate Data Distributin Scheme Halis Altun, Gökhan Gelen Nigde University, Electrical and Electrnics Engineering Department Nigde, Turkey haltun@nigde.edu.tr
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 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 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 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 informationECE 5318/6352 Antenna Engineering. Spring 2006 Dr. Stuart Long. Chapter 6. Part 7 Schelkunoff s Polynomial
ECE 538/635 Antenna Engineering Spring 006 Dr. Stuart Lng Chapter 6 Part 7 Schelkunff s Plynmial 7 Schelkunff s Plynmial Representatin (fr discrete arrays) AF( ψ ) N n 0 A n e jnψ N number f elements in
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
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 informationCS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007
CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is
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 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 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 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 informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins
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 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 informationChapter 2 GAUSS LAW Recommended Problems:
Chapter GAUSS LAW Recmmended Prblems: 1,4,5,6,7,9,11,13,15,18,19,1,7,9,31,35,37,39,41,43,45,47,49,51,55,57,61,6,69. LCTRIC FLUX lectric flux is a measure f the number f electric filed lines penetrating
More informationThe Creation and Propagation of Radiation: Fields Inside and Outside of Sources
Versin Date July 10, 011 1 The Creatin and Prpagatin f Radiatin: Fields Inside and Outside f Surces Stanislaw Olbert and Jhn W. Belcher Department f Physics Massachusetts Institute f Technlgy Richard H.
More informationModule 4: General Formulation of Electric Circuit Theory
Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated
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 informationELECTROMAGNETIC CHARACTERISTICS OF CON- FORMAL DIPOLE ANTENNAS OVER A PEC SPHERE
Prgress In Electrmagnetics Research M, Vl. 26, 85 1, 212 ELECTROMAGNETIC CHARACTERISTICS OF CON- FORMAL DIPOLE ANTENNAS OVER A PEC SPHERE J. S. Meiguni *, M. Kamyab, and A. Hsseinbeig Faculty f Electrical
More informationEXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE
EXPERIMENTAL STUD ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE Tmnbu Gt, Masaaki Ohba, Takashi Kurabuchi 2, Tmyuki End 3, shihik Akamine 4, and Tshihir Nnaka 2
More informationInterference is when two (or more) sets of waves meet and combine to produce a new pattern.
Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme
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 informationCHAPTER 8b Static Equilibrium Units
CHAPTER 8b Static Equilibrium Units The Cnditins fr Equilibrium Slving Statics Prblems Stability and Balance Elasticity; Stress and Strain The Cnditins fr Equilibrium An bject with frces acting n it, but
More informationBeam vibrations: Discrete mass and stiffness models
Beam vibratins: Discrete mass and stiffness mdels Ana Cláudia Susa Neves ana.neves@tecnic.ulisba.pt Institut Superir Técnic, Universidade de Lisba, Prtugal May, 2015 Abstract In the present wrk the dynamic
More informationNUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION
NUROP Chinese Pinyin T Chinese Character Cnversin NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION CHIA LI SHI 1 AND LUA KIM TENG 2 Schl f Cmputing, Natinal University f Singapre 3 Science
More informationCOASTAL ENGINEERING Chapter 2
CASTAL ENGINEERING Chapter 2 GENERALIZED WAVE DIFFRACTIN DIAGRAMS J. W. Jhnsn Assciate Prfessr f Mechanical Engineering University f Califrnia Berkeley, Califrnia INTRDUCTIN Wave diffractin is the phenmenn
More informationSection 6-2: Simplex Method: Maximization with Problem Constraints of the Form ~
Sectin 6-2: Simplex Methd: Maximizatin with Prblem Cnstraints f the Frm ~ Nte: This methd was develped by Gerge B. Dantzig in 1947 while n assignment t the U.S. Department f the Air Frce. Definitin: Standard
More information3D FE Modeling Simulation of Cold Rotary Forging with Double Symmetry Rolls X. H. Han 1, a, L. Hua 1, b, Y. M. Zhao 1, c
Materials Science Frum Online: 2009-08-31 ISSN: 1662-9752, Vls. 628-629, pp 623-628 di:10.4028/www.scientific.net/msf.628-629.623 2009 Trans Tech Publicatins, Switzerland 3D FE Mdeling Simulatin f Cld
More informationRunning head: Weighted ENO Schemes Prfs sent t: Xu-Dng Liu Curant Institute f Mathematical Sciences 5 Mercer Street New Yrk, NY 00 AMS(MOS) subject cl
Weighted Essentially Nn-Oscillatry Schemes Xu-Dng Liu Stanley Osher y Tny Chan z Department f Mathematics, UCLA, Ls Angeles, Califrnia, 9004. Research supprted by NSF grant DMS 9-043. y Department f Mathematics,
More information1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp
THE POWER AND LIMIT OF NEURAL NETWORKS T. Y. Lin Department f Mathematics and Cmputer Science San Jse State University San Jse, Califrnia 959-003 tylin@cs.ssu.edu and Bereley Initiative in Sft Cmputing*
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 information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion
.54 Neutrn Interactins and Applicatins (Spring 004) Chapter (3//04) Neutrn Diffusin References -- J. R. Lamarsh, Intrductin t Nuclear Reactr Thery (Addisn-Wesley, Reading, 966) T study neutrn diffusin
More informationPreparation work for A2 Mathematics [2018]
Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time
More informationNUMBERS, MATHEMATICS AND EQUATIONS
AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t
More informationHomology groups of disks with holes
Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.
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 information5 th grade Common Core Standards
5 th grade Cmmn Cre Standards In Grade 5, instructinal time shuld fcus n three critical areas: (1) develping fluency with additin and subtractin f fractins, and develping understanding f the multiplicatin
More informationSoliton-Effect Optical Pulse Compression in Bulk Media with χ (3) Nonlinearity. 1 Introduction
Nnlinear Analysis: Mdelling and Cntrl, Vilnius, IMI,, N 5 Lithuanian Assciatin f Nnlinear Analysts, Slitn-Effect Optical Pulse Cmpressin in Bulk Media with χ (3) Nnlinearity Received: 9.7. Accepted: 11.1.
More informationPipetting 101 Developed by BSU CityLab
Discver the Micrbes Within: The Wlbachia Prject Pipetting 101 Develped by BSU CityLab Clr Cmparisns Pipetting Exercise #1 STUDENT OBJECTIVES Students will be able t: Chse the crrect size micrpipette fr
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 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 informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More informationFree Vibrations of Catenary Risers with Internal Fluid
Prceeding Series f the Brazilian Sciety f Applied and Cmputatinal Mathematics, Vl. 4, N. 1, 216. Trabalh apresentad n DINCON, Natal - RN, 215. Prceeding Series f the Brazilian Sciety f Cmputatinal and
More information3.4 Shrinkage Methods Prostate Cancer Data Example (Continued) Ridge Regression
3.3.4 Prstate Cancer Data Example (Cntinued) 3.4 Shrinkage Methds 61 Table 3.3 shws the cefficients frm a number f different selectin and shrinkage methds. They are best-subset selectin using an all-subsets
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 informationPOLARISATION VISUAL PHYSICS ONLINE. View video on polarisation of light
VISUAL PHYSICS ONLINE MODULE 7 NATURE OF LIGHT POLARISATION View vide n plarisatin f light While all the experimental evidence s far that supprts the wave nature f light, nne f it tells us whether light
More informationABSORPTION OF GAMMA RAYS
6 Sep 11 Gamma.1 ABSORPTIO OF GAMMA RAYS Gamma rays is the name given t high energy electrmagnetic radiatin riginating frm nuclear energy level transitins. (Typical wavelength, frequency, and energy ranges
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 informationModelling of NOLM Demultiplexers Employing Optical Soliton Control Pulse
Micwave and Optical Technlgy Letters, Vl. 1, N. 3, 1999. pp. 05-08 Mdelling f NOLM Demultiplexers Emplying Optical Slitn Cntrl Pulse Z. Ghassemly, C. Y. Cheung & A. K. Ray Electrnics Research Grup, Schl
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 informationSupporting information
Electrnic Supplementary Material (ESI) fr Physical Chemistry Chemical Physics This jurnal is The wner Scieties 01 ydrgen perxide electrchemistry n platinum: twards understanding the xygen reductin reactin
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 informationLyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
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 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 informationSubject description processes
Subject representatin 6.1.2. Subject descriptin prcesses Overview Fur majr prcesses r areas f practice fr representing subjects are classificatin, subject catalging, indexing, and abstracting. The prcesses
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 informationSynchronous Motor V-Curves
Synchrnus Mtr V-Curves 1 Synchrnus Mtr V-Curves Intrductin Synchrnus mtrs are used in applicatins such as textile mills where cnstant speed peratin is critical. Mst small synchrnus mtrs cntain squirrel
More informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
More informationLegendre-Tau Spectral Elements for Incompressible Navier-Stokes Flow
Legendre-Tau Spectral Elements fr Incmpressible avier-stkes Flw Kelly Black* Abstract A spectral multi-dmain methd is intrduced and examined. After dividing the cmputatinal dmain int nnverlapping subdmains
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 informationCOMP 551 Applied Machine Learning Lecture 11: Support Vector Machines
COMP 551 Applied Machine Learning Lecture 11: Supprt Vectr Machines Instructr: (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted fr this curse
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 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 informationGENESIS Structural Optimization for ANSYS Mechanical
P3 STRUCTURAL OPTIMIZATION (Vl. II) GENESIS Structural Optimizatin fr ANSYS Mechanical An Integrated Extensin that adds Structural Optimizatin t ANSYS Envirnment New Features and Enhancements Release 2017.03
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 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 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 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 informationA New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation
III-l III. A New Evaluatin Measure J. Jiner and L. Werner Abstract The prblems f evaluatin and the needed criteria f evaluatin measures in the SMART system f infrmatin retrieval are reviewed and discussed.
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 information1 The limitations of Hartree Fock approximation
Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants
More informationBiocomputers. [edit]scientific Background
Bicmputers Frm Wikipedia, the free encyclpedia Bicmputers use systems f bilgically derived mlecules, such as DNA and prteins, t perfrm cmputatinal calculatins invlving string, retrieving, and prcessing
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 informationI. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is
Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,
More informationIntroduction to Smith Charts
Intrductin t Smith Charts Dr. Russell P. Jedlicka Klipsch Schl f Electrical and Cmputer Engineering New Mexic State University as Cruces, NM 88003 September 2002 EE521 ecture 3 08/22/02 Smith Chart Summary
More informationCHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India
CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce
More informationCurrent/voltage-mode third order quadrature oscillator employing two multiple outputs CCIIs and grounded capacitors
Indian Jurnal f Pure & Applied Physics Vl. 49 July 20 pp. 494-498 Current/vltage-mde third rder quadrature scillatr emplying tw multiple utputs CCIIs and grunded capacitrs Jiun-Wei Hrng Department f Electrnic
More informationThe Sputtering Problem James A Glackin, James V. Matheson
The Sputtering Prblem James A Glackin, James V. Mathesn I prpse t cnsider first the varius elements f the subject, next its varius parts r sectins, and finally the whle in its internal structure. In ther
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 informationA New Approach to Increase Parallelism for Dependent Loops
A New Apprach t Increase Parallelism fr Dependent Lps Yeng-Sheng Chen, Department f Electrical Engineering, Natinal Taiwan University, Taipei 06, Taiwan, Tsang-Ming Jiang, Arctic Regin Supercmputing Center,
More informationA TIME-DOMAIN FINITE-DIFFERENCE METHOD WITH ATTENUATION BY A RECURSIVE ALGORITHM
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
More information