METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
|
|
- Marilyn Higgins
- 5 years ago
- Views:
Transcription
1 Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009, Volume 8, Issue, pages The website: Scietific Research of the Istitute of Mathematics ad Computer Sciece METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS Staisław Kula Istitute of Mathematics ad Computer Sciece, Czestochowa Uiversity of Techology, Polad ula@imi.pcz.pl Abstract. I this paper, the method of fudametal solutios for Helmholtz eigeproblems i a elliptical domai is preseted. To fid the approximate solutio of the problem, the Hael fuctio of the first id ad zero order as the fudametal solutio of the Helmholtz equatio i ubouded domai o the plae was used. Numerical examples illustratig the accuracy of the preset method are give. Itroductio The trasverse vibratio of a membrae which occupy the domai S, is govered by the followig equatio [] where is the Laplace operator, w w = () c t c = T ρ is the speed of soud, T is the uiform tesio per uit legth of the edge S, w x, y, t is the displacemet of the plate poit (x,y) at time t. Whe the free vibratio of the i t membrae is cosidered, a harmoic time depedece of the form e ω is assumed, t i.e. (,, ) (, ) e i ω w x y t = W x y. Substitutig the form of the displacemet ito (), the Helmholtz equatio is obtaied where Ω = ω ρ T assume here the Dirichlet coditio + W = ρ is the desity per uit area, ( ) W Ω, ( x y) S 0, (). The equatio () is completed by boudary coditios. We W ( x y) = 0, ( x, y) S, (3) The differetial equatio () ad boudary coditio (3) form the Helmholtz eigevalue problem. The problem for elliptic domai S is the subject of this paper.
2 86 S. Kula. The method of fudametal solutios We solve the eigevalue problem ()-(3) by usig the method of fudametal solutios (MFS). The MFS is a boudary method which does ot ivolve discretizatio ad itegratio. The method is applicable to problems for which a fudametal solutio of the goverig equatio is ow. The basic idea of the method is the usage of a liear combiatio of fudametal solutios with sources located at fictitious poits outside the domai of the cosidered problem. The fudametal solutio of a differetial equatio is ot uique. For the Helmholtz equatio () we tae the Hael fuctio [, 3] where i ( ) ( ) G P, Q = H 0 ( Ω d( P, Q) ) (4) 4 ( H ) 0 ( ) is the Hael fuctio of the first id ad zero order, ( P Q) i =. The fuctios G ( x y; ) distace betwee poits P ad Q, Q d, is the, satisfy the Helmholtz equatio i the domai S for each source poits Q located outside S. I the MFS, we approximate the solutio of the problem by a fuctio of the form w ( P) = c G( P Q ) =, (5) The approximate solutio w satisfies the differetial equatio (), ad it does ot satisfy the boudary coditio (3). The coditio ca be satisfied approximately by a suitable determiatio of the coefficiets c, =,,,. For this purpose we use the last square method. First we choose the poits P i, i =,,,, located o boudary S of the domai S. We determie the coefficiets c so that the fuctio f (( c c,..., c )) = [ w ( Pi ; c, c,..., c )] i=, (6) has a miimum. This leads to a homogeeous liear system of equatios f c j = c = i= (, Q ) G( P, Q ) = 0, G P i j i j =,,, (7) which ca be writte i the matrix form as A c = 0 (8) where A = [ a ], j a ( ) ( ) j, j G Pi, Q j G Pi, Q =, = [ c c ] T i= c... c
3 Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais 87 For a o-trivial solutio of equatio (8), the determiat of the coefficiet matrix A is set equal to zero, yieldig the eigevalue equatio det A(Ω ) = 0 (9) Equatio (9) with the uow Ω, is solved umerically. The eigefuctios for the determied eigevalues Ω m, m =,,, are give by (5) where the coefficiets c, =,,, are derived i depedece of c from - equatios of the system (7). For the cosidered Helmholtz eigevalue problem the fuctio w defied by (5), after usig (4) has the form w ( ) ( ) x, y = c H Ω ( x ξ ) + ( y η ), ( x y) S = 0 where the source poits Q ( ξ, η ) coefficiets of the matrix A are a j = i= H, (0), =,,,, are over the domai S. But the Ω x i ξ j + yi η j H0 Ω xi ξ + yi η () ( ) ( ) ( ) ( ) ( ) ( ) 0 where the collocatio poits ( ) S of the domai S. P x i, y i, i =,,,, are located o the boudary. Numerical examples We assume that the boudary of the cosidered membrae domai is a ellipse with major ad mior semi-diameters equal a ad b, respectively. The collocatio poits are uiformly distributed o this ellipse ad the source poits are uiformly distributed o a circle with the radius r = α max{ a, b}, α >. Geometry cofiguratio of a elliptical domai with the collocatio poits o the boudary ad a circle with the source poits is preseted i Figure. The Hael fuctio ( H ) 0 occurrig i equatios (0, ) assumes complex values ad that way the elemets of the matrix A are complex. Our goal is to evaluate real roots of the complex equatio (9). For this purpose the ew fuctio F is defied as where the symbol ( Ω ) = det A( Ω ) F () deotes a modulus of complex umber. It is clear that the roots of equatio (9) are determied by miima of the fuctio F. I Figure a
4 88 S. Kula exemplary graph of the fuctio i logarithmic scale for a elliptical domai with semi-diameters ratio a/b =.0 ad = 4 is preseted. Locatios of the distict drops idicate the roots of equatio (9). y Q j = (ξ j,η j ) b P i = (x i, y i ) a x Fig.. Geometry cofiguratio of the cosidered domai with collocatio poits P i o the ellipse ad source poits o Q j the circle Ω 70 log det A( Ω ) Fig.. Logarithmic values of modulus of the determiat deta(ω ) as a fuctio of Ω for the Helmholtz operator i the elliptical domai with semi-diameters ratio a/b =.0
5 Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais 89 To estimate the MFS accuracy, a example of the Helmholtz eigevalue problem i a circular domai will be solved. I this case the eigevalues are the roots of the equatio [3] m ( Ω ) = 0 J, m =,, (3) A sequece of the eigevalues for each m ca be obtaied. First members of the sequeces are umerically determied by fidig the local miima of the fuctio F Ω F Ω is give by equatio (). The calculatios were performed log ( ), where ( ) for various umbers of source poits ( =6; 0; 4). The compariso of the results obtaied by the MFS with exact values obtaied as roots of equatio (3) is preseted i Table. Note that the umerical calculatio by MFS with small umber of the source poits gives the results with a larger error ad ca lead to omissio of eigevalues. The compariso of the results obtaied by MFS ad the exact eigevalues show that the absolute error for = 4 is ot larger as 0 5. The eigevalues of the Helmholtz operator are the frequecy parameters of free vibratio of a membrae. The first seve frequecy parameter values of the elliptic membrae with clamped edge are preseted i Table. The calculatios are performed for various values of semi-diameters ratio a/b of the ellipse. The frequecy parameters obtaied by usig MFS were compared with the eigevalues determied by the fiite elemet method i the paper [] by Buchaa ad Peddieso. For assumed umber of sources = 4 i MFS a small differeces of the results calculated by usig both methods are observed. Table. Eigevalues Ω ( =,,0) of the Helmholtz operator i a circular domai obtaied by the MFS for = 6; 0; 4 ad exact eigevalues obtaied as roots of equatio () MFS = 6 MFS = 0 MFS = 4 Exact eigevalues Ω Ω Ω Ω Ω Ω Ω Ω Ω Ω
6 90 S. Kula Table. Frequecy parameters Ω ( =,,7) of the elliptic membrae for various values of semi-diameters ratio b/a of the ellipse a/b =.5 a/b =.0 a/b = 3.0 MFS FEM * MFS FEM * MFS FEM * Ω Ω Ω Ω Ω Ω Ω the results by FEM are give i the paper [4] Coclusios I the paper, the applicatio of the method of fudametal solutios to the Helmholtz eigevalue problem i the elliptical domai has bee preseted. As the fudametal solutio of the Helmholtz equatio i the plae, the Hael fuctio of the first id ad zero order was used. I order to determie the eigevalues, the miimum of a real fuctio which is the logarithm of modulus of a complex fuctio was foud. The source poits occurrig i the approximate formula of the solutio ca be selected o a circle outside of the cosidered elliptical domai. The compariso of umerical results shows that high accuracy of the calculatio is achieved for 4 sources. Refereces [] Pietrowsi L., Rówaia róŝiczowe cząstowe, PWN, Warszawa 955. [] Che J.T., Che I.L., Lee Y.T., Eigesolutios of multiply coected membraes usig the method of fudametal solutios, Boudary Aalysis with Boudary Elemets 005, 9, [3] Reutsiy S.Yu., The method of fudametal solutios for Helmholtz eigevalue problems i simply ad multiply coected domais, Boudary Aalysis with Boudary Elemets 006, 30, [4] Buchaa G.R., Peddieso J., A fiite elemet i eliptic coordiates with applicatio to membrae vibratio, Thi-Walled Structures 005, 43, [5] Duffy D.G., Gree s fuctios with applicatios, Chapma & Hall/CRC, Boca Rato, Washigto DC 00.
Damped Vibration of a Non-prismatic Beam with a Rotational Spring
Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationAnalysis of composites with multiple rigid-line reinforcements by the BEM
Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More informationNumerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials
IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationUNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA
UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA Naval Uiversity of Gdyia 81-13 Gdyia, Śmidowicza 69, Polad Gdańsk Uiversity of Techology 8-95 Gdańsk, Narutowicza 11/1,
More informationANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION
Molecular ad Quatum Acoustics vol. 7, (6) 79 ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION Jerzy FILIPIAK 1, Lech SOLARZ, Korad ZUBKO 1 Istitute of Electroic ad Cotrol Systems, Techical Uiversity of Czestochowa,
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationNotes 18 Green s Functions
ECE 638 Fall 017 David R. Jackso Notes 18 Gree s Fuctios Notes are from D. R. Wilto, Dept. of ECE 1 Gree s Fuctios The Gree's fuctio method is a powerful ad systematic method for determiig a solutio to
More informationAvoidance of numerical singularities in free vibration analysis of Euler-Bernoulli beams using Green functions
WSEAS TRASACTIOS o APPLIED ad THEORETICAL MECHAICS Goraka Štimac Ročević, Braimir Ročević, Ate Skoblar, Saji Braut Avoidace of umerical sigularities i free vibratio aalysis of Euler-Beroulli beams usig
More informationPoincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains
Advaces i Pure Mathematics 23 3 72-77 http://dxdoiorg/4236/apm233a24 Published Olie Jauary 23 (http://wwwscirporg/oural/apm) Poicaré Problem for Noliear Elliptic Equatios of Secod Order i Ubouded Domais
More informationDynamical stability and parametrical vibrations of the laminated plates with complex shape
(3) 75 88 Dyamical stability ad parametrical vibratios of the lamiated plates with comple shape Abstract The problem of oliear vibratios ad stability aalysis for the symmetric lamiated plates with comple
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationFREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES
Abdullah : FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS USING FOURIER SERIES SAWA MUBARAK ABDUAH Assistat ecturer Uiversity of Mosul Abstract Fourier series will be
More informationMath 220B Final Exam Solutions March 18, 2002
Math 0B Fial Exam Solutios March 18, 00 1. (1 poits) (a) (6 poits) Fid the Gree s fuctio for the tilted half-plae {(x 1, x ) R : x 1 + x > 0}. For x (x 1, x ), y (y 1, y ), express your Gree s fuctio G(x,
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationInverse Nodal Problems for Differential Equation on the Half-line
Australia Joural of Basic ad Applied Scieces, 3(4): 4498-4502, 2009 ISSN 1991-8178 Iverse Nodal Problems for Differetial Equatio o the Half-lie 1 2 3 A. Dabbaghia, A. Nematy ad Sh. Akbarpoor 1 Islamic
More informationChapter 7 z-transform
Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time
More informationSolution of EECS 315 Final Examination F09
Solutio of EECS 315 Fial Examiatio F9 1. Fid the umerical value of δ ( t + 4ramp( tdt. δ ( t + 4ramp( tdt. Fid the umerical sigal eergy of x E x = x[ ] = δ 3 = 11 = ( = ramp( ( 4 = ramp( 8 = 8 [ ] = (
More informationUnit 4: Polynomial and Rational Functions
48 Uit 4: Polyomial ad Ratioal Fuctios Polyomial Fuctios A polyomial fuctio y px ( ) is a fuctio of the form p( x) ax + a x + a x +... + ax + ax+ a 1 1 1 0 where a, a 1,..., a, a1, a0are real costats ad
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationLecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables
Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio
More information: Transforms and Partial Differential Equations
Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
More informationRay Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET
Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationCh3 Discrete Time Fourier Transform
Ch3 Discrete Time Fourier Trasform 3. Show that the DTFT of [] is give by ( k). e k 3. Determie the DTFT of the two sided sigal y [ ],. 3.3 Determie the DTFT of the causal sequece x[ ] A cos( 0 ) [ ],
More information8. Applications To Linear Differential Equations
8. Applicatios To Liear Differetial Equatios 8.. Itroductio 8.. Review Of Results Cocerig Liear Differetial Equatios Of First Ad Secod Orders 8.3. Eercises 8.4. Liear Differetial Equatios Of Order N 8.5.
More informationChapter 5 Vibrational Motion
Fall 4 Chapter 5 Vibratioal Motio... 65 Potetial Eergy Surfaces, Rotatios ad Vibratios... 65 Harmoic Oscillator... 67 Geeral Solutio for H.O.: Operator Techique... 68 Vibratioal Selectio Rules... 7 Polyatomic
More informationGeneralizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations
Geeraliig the DTFT The Trasform M. J. Roberts - All Rights Reserved. Edited by Dr. Robert Akl 1 The forward DTFT is defied by X e jω = x e jω i which = Ω is discrete-time radia frequecy, a real variable.
More informationChapter 7 Isoperimetric problem
Chapter 7 Isoperimetric problem Recall that the isoperimetric problem (see the itroductio its coectio with ido s proble) is oe of the most classical problem of a shape optimizatio. It ca be formulated
More informationNUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS
THERMAL SCIENCE, Year 07, Vol., No. 4, pp. 595-599 595 NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS by Yula WANG *, Da TIAN, ad Zhiyua LI Departmet of Mathematics,
More informationAPPENDIX F Complex Numbers
APPENDIX F Complex Numbers Operatios with Complex Numbers Complex Solutios of Quadratic Equatios Polar Form of a Complex Number Powers ad Roots of Complex Numbers Operatios with Complex Numbers Some equatios
More informationFINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,
More informationNBHM QUESTION 2007 Section 1 : Algebra Q1. Let G be a group of order n. Which of the following conditions imply that G is abelian?
NBHM QUESTION 7 NBHM QUESTION 7 NBHM QUESTION 7 Sectio : Algebra Q Let G be a group of order Which of the followig coditios imply that G is abelia? 5 36 Q Which of the followig subgroups are ecesarily
More informationProject 10.3 Vibrating Beams and Diving Boards
Project 10.3 Vibratig Beams ad Divig Boards I this project you are to ivestigate further the vibratios of a elastic bar or beam of legth L whose positio fuctio y(x,t) satisfies the partial differetial
More informationFree Vibration Analysis of Beams Considering Different Geometric Characteristics and Boundary Conditions
Iteratioal Joural of Mechaics ad Applicatios 1, (3): 9-1 DOI: 1.593/j.mechaics.13.3 Free Vibratio Aalysis of Beams Cosiderig Differet Geometric Characteristics ad Boudary Coditios Mehmet Avcar Civil Egieerig
More informationComputation of Error Bounds for P-matrix Linear Complementarity Problems
Mathematical Programmig mauscript No. (will be iserted by the editor) Xiaoju Che Shuhuag Xiag Computatio of Error Bouds for P-matrix Liear Complemetarity Problems Received: date / Accepted: date Abstract
More informationDYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS
DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS Ivaa Štimac 1, Ivica Kožar 1 M.Sc,Assistat, Ph.D. Professor 1, Faculty of Civil Egieerig, Uiverity of Rieka, Croatia INTRODUCTION The vehicle-iduced
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationMonte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem
Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic
More informationMath 5311 Problem Set #5 Solutions
Math 5311 Problem Set #5 Solutios March 9, 009 Problem 1 O&S 11.1.3 Part (a) Solve with boudary coditios u = 1 0 x < L/ 1 L/ < x L u (0) = u (L) = 0. Let s refer to [0, L/) as regio 1 ad (L/, L] as regio.
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationA Lattice Green Function Introduction. Abstract
August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:
More informationIIT JAM Mathematical Statistics (MS) 2006 SECTION A
IIT JAM Mathematical Statistics (MS) 6 SECTION A. If a > for ad lim a / L >, the which of the followig series is ot coverget? (a) (b) (c) (d) (d) = = a = a = a a + / a lim a a / + = lim a / a / + = lim
More informationGROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN
The th World Coferece o Earthquake Egieerig October -7, 8, Beiig, Chia GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE Hui QI, Yog SHI, Jigfu NAN ABSTRACT : Professor, Dept. of
More informationThe Numerical Solution of Singular Fredholm Integral Equations of the Second Kind
WDS' Proceedigs of Cotributed Papers, Part I, 57 64, 2. ISBN 978-8-7378-39-2 MATFYZPRESS The Numerical Solutio of Sigular Fredholm Itegral Equatios of the Secod Kid J. Rak Charles Uiversity, Faculty of
More informationOn the Weak Localization Principle of the Eigenfunction Expansions of the Laplace-Beltrami Operator by Riesz Method ABSTRACT 1.
Malaysia Joural of Mathematical Scieces 9(): 337-348 (05) MALAYSIA JOURAL OF MATHEMATICAL SCIECES Joural homepage: http://eispemupmedumy/joural O the Weak Localizatio Priciple of the Eigefuctio Expasios
More informationA) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set.
M.A./M.Sc. (Mathematics) Etrace Examiatio 016-17 Max Time: hours Max Marks: 150 Istructios: There are 50 questios. Every questio has four choices of which exactly oe is correct. For correct aswer, 3 marks
More informationWarped, Chirp Z-Transform: Radar Signal Processing
arped, Chirp Z-Trasform: Radar Sigal Processig by Garimella Ramamurthy Report o: IIIT/TR// Cetre for Commuicatios Iteratioal Istitute of Iformatio Techology Hyderabad - 5 3, IDIA Jauary ARPED, CHIRP Z
More informationChapter 10 Partial Differential Equations and Fourier Series
Math-33 Chapter Partial Differetial Equatios November 6, 7 Chapter Partial Differetial Equatios ad Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7. Boudary Value Problems for
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationMTH Assignment 1 : Real Numbers, Sequences
MTH -26 Assigmet : Real Numbers, Sequeces. Fid the supremum of the set { m m+ : N, m Z}. 2. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a
More informationEstimation of Backward Perturbation Bounds For Linear Least Squares Problem
dvaced Sciece ad Techology Letters Vol.53 (ITS 4), pp.47-476 http://dx.doi.org/.457/astl.4.53.96 Estimatio of Bacward Perturbatio Bouds For Liear Least Squares Problem Xixiu Li School of Natural Scieces,
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS
July 3. Vol. 4, No. Iteratioal Joural of Egieerig ad Applied Scieces EAAS & ARF. All rights reserved www.eaas-joural.org ISSN35-869 ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S
More informationON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS
July 3 Vol 4, No Iteratioal Joural of Egieerig ad Applied Scieces EAAS & ARF All rights reserved wwweaas-jouralorg ISSN35-869 ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS
More informationPHY4905: Nearly-Free Electron Model (NFE)
PHY4905: Nearly-Free Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. No-degeerate eigestates
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More informationSolutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.
Quiz November 4th, 23 Sigals & Systems (5-575-) P. Reist & Prof. R. D Adrea Solutios Exam Duratio: 4 miutes Number of Problems: 4 Permitted aids: Noe. Use oly the prepared sheets for your solutios. Additioal
More informationREAL ANALYSIS II: PROBLEM SET 1 - SOLUTIONS
REAL ANALYSIS II: PROBLEM SET 1 - SOLUTIONS 18th Feb, 016 Defiitio (Lipschitz fuctio). A fuctio f : R R is said to be Lipschitz if there exists a positive real umber c such that for ay x, y i the domai
More informationa b c d e f g h Supplementary Information
Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;
More informationA COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS
A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS Y. S. PARK a adw.t.ang b Computatioal ad Mathematical Scieces Group Faculty of Iformatio
More informationSolution of Differential Equation from the Transform Technique
Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from
More informationMATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and
MATH01 Real Aalysis (2008 Fall) Tutorial Note #7 Sequece ad Series of fuctio 1: Poitwise Covergece ad Uiform Covergece Part I: Poitwise Covergece Defiitio of poitwise covergece: A sequece of fuctios f
More informationOn forward improvement iteration for stopping problems
O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal
More informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More informationIN many scientific and engineering applications, one often
INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS, VOL 3, NO, FEBRUARY 07 5 Secod Degree Refiemet Jacobi Iteratio Method for Solvig System of Liear Equatio Tesfaye Kebede Abstract Several
More information1. Hydrogen Atom: 3p State
7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More informationReview Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =
Review Problems ICME ad MS&E Refresher Course September 9, 0 Warm-up problems. For the followig matrices A = 0 B = C = AB = 0 fid all powers A,A 3,(which is A times A),... ad B,B 3,... ad C,C 3,... Solutio:
More informationSynopsis of Euler s paper. E Memoire sur la plus grande equation des planetes. (Memoir on the Maximum value of an Equation of the Planets)
1 Syopsis of Euler s paper E105 -- Memoire sur la plus grade equatio des plaetes (Memoir o the Maximum value of a Equatio of the Plaets) Compiled by Thomas J Osler ad Jase Adrew Scaramazza Mathematics
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationAssignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1
Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate
More information2C09 Design for seismic and climate changes
2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationAnalytical solutions for multi-wave transfer matrices in layered structures
Joural of Physics: Coferece Series PAPER OPEN ACCESS Aalytical solutios for multi-wave trasfer matrices i layered structures To cite this article: Yu N Belyayev 018 J Phys: Cof Ser 109 01008 View the article
More information