Fourier-Bessel Expansions with Arbitrary Radial Boundaries

Similar documents
Review of Mathematical Concepts

9.4 The response of equilibrium to temperature (continued)

Integrals and Polygamma Representations for Binomial Sums

Physics 505 Fall 2005 Midterm Solutions. This midterm is a two hour open book, open notes exam. Do all three problems.

Available online at ScienceDirect. Procedia Engineering 91 (2014 ) 32 36

Discrete Model Parametrization

( ) D x ( s) if r s (3) ( ) (6) ( r) = d dr D x

10 Statistical Distributions Solutions

About Some Inequalities for Isotonic Linear Functionals and Applications

Two dimensional polar coordinate system in airy stress functions

FI 2201 Electromagnetism

On the Eötvös effect

U>, and is negative. Electric Potential Energy

Important design issues and engineering applications of SDOF system Frequency response Functions

π,π is the angle FROM a! TO b

Homework: Study 6.2 #1, 3, 5, 7, 11, 15, 55, 57

Friedmannien equations

Physics 604 Problem Set 1 Due Sept 16, 2010

Homework 3 MAE 118C Problems 2, 5, 7, 10, 14, 15, 18, 23, 30, 31 from Chapter 5, Lamarsh & Baratta. The flux for a point source is:

Chapter 6 Thermoelasticity

Electric Potential. and Equipotentials

Week 8. Topic 2 Properties of Logarithms

Radial geodesics in Schwarzschild spacetime

3.1 Magnetic Fields. Oersted and Ampere

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

MATHEMATICS IV 2 MARKS. 5 2 = e 3, 4

Chapter 7. Kleene s Theorem. 7.1 Kleene s Theorem. The following theorem is the most important and fundamental result in the theory of FA s:

Michael Rotkowitz 1,2

A COMPARISON OF MEMBRANE SHELL THEORIES OF HYBRID ANISOTROPIC MATERIALS ABSTRACT

Mark Scheme (Results) January 2008

EECE 260 Electrical Circuits Prof. Mark Fowler

dx was area under f ( x ) if ( ) 0

Section 35 SHM and Circular Motion

Physics 1502: Lecture 2 Today s Agenda

Electricity & Magnetism Lecture 6: Electric Potential

Lectures # He-like systems. October 31 November 4,6

Qualitative Analysis for Solutions of a Class of. Nonlinear Ordinary Differential Equations

u(r, θ) = 1 + 3a r n=1

Continuous Charge Distributions

(a) Counter-Clockwise (b) Clockwise ()N (c) No rotation (d) Not enough information

ELECTRO - MAGNETIC INDUCTION

This immediately suggests an inverse-square law for a "piece" of current along the line.

Introductions to ArithmeticGeometricMean

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 20

1 Using Integration to Find Arc Lengths and Surface Areas

Lecture 11: Potential Gradient and Capacitor Review:

General Physics II. number of field lines/area. for whole surface: for continuous surface is a whole surface

Study of Electromagnetic Wave Propagation in Periodic Dielectric Structure; MathCAD Analysis

Right-indefinite half-linear Sturm Liouville problems

RELATIVE KINEMATICS. q 2 R 12. u 1 O 2 S 2 S 1. r 1 O 1. Figure 1

Class Summary. be functions and f( D) , we define the composition of f with g, denoted g f by

PX3008 Problem Sheet 1

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING FLUID MECHANICS III Solutions to Problem Sheet 3

Review of Gaussian Quadrature method

r = (0.250 m) + (0.250 m) r = m = = ( N m / C )

Investigations of Boundary Treatments in Incompressible Smoothed Particle Hydrodynamics for Fluid-Structural Interactions

STUDY OF THE UNIFORM MAGNETIC FIELD DOMAINS (3D) IN THE CASE OF THE HELMHOLTZ COILS

Mathematical formulation of the F 0 motor model

STD: XI MATHEMATICS Total Marks: 90. I Choose the correct answer: ( 20 x 1 = 20 ) a) x = 1 b) x =2 c) x = 3 d) x = 0

7.5-Determinants in Two Variables

Optimization. x = 22 corresponds to local maximum by second derivative test

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

Unit 6. Magnetic forces

CHAPTER 2 ELECTROSTATIC POTENTIAL

Previously. Extensions to backstepping controller designs. Tracking using backstepping Suppose we consider the general system

PROGRESSION AND SERIES

2 mv2 qv (0) = 0 v = 2qV (0)/m. Express q. . Substitute for V (0) and simplify to obtain: v = q

Topics for Review for Final Exam in Calculus 16A

Lecture 10. Solution of Nonlinear Equations - II

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

r a + r b a + ( r b + r c)

6. Numbers. The line of numbers: Important subsets of IR:

SPA7010U/SPA7010P: THE GALAXY. Solutions for Coursework 1. Questions distributed on: 25 January 2018.

Chapter 23 Electrical Potential

International Journal of Pure and Applied Sciences and Technology

Refutation of comment by Jadczyk et al

10.3 The Quadratic Formula

Polynomials and Division Theory

An Alternative Approach to Estimating the Bounds of the Denominators of Egyptian Fractions

Summary: Binomial Expansion...! r. where

CHAPTER 18: ELECTRIC CHARGE AND ELECTRIC FIELD

Language Processors F29LP2, Lecture 5

defined on a domain can be expanded into the Taylor series around a point a except a singular point. Also, f( z)

10 m, so the distance from the Sun to the Moon during a solar eclipse is. The mass of the Sun, Earth, and Moon are = =

AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir

Electronic Supplementary Material

Comparative Studies of Law of Gravity and General Relativity. No.1 of Comparative Physics Series Papers

Elastic limit angular speed of solid and annular disks under thermomechanical

Ch 26 - Capacitance! What s Next! Review! Lab this week!

The Formulas of Vector Calculus John Cullinan

Prerna Tower, Road No 2, Contractors Area, Bistupur, Jamshedpur , Tel (0657) ,

Tests for Correlation on Bivariate Non-Normal Data

Algebra Based Physics. Gravitational Force. PSI Honors universal gravitation presentation Update Fall 2016.notebookNovember 10, 2016

1B40 Practical Skills

Elastic scattering of 4 He atoms at the surface of liquid helium

Fluids & Bernoulli s Equation. Group Problems 9

Dynamically Equivalent Systems. Dynamically Equivalent Systems. Dynamically Equivalent Systems. ME 201 Mechanics of Machines

igid nd non-leky two-comptment building. Yu et l [8] developed non-line govening equtions by consideing the effect of bckgound lekge. Howeve, thee e n

Wave Generation by Oscillating Wall in Static Media

Effect of Heat Generation on Quasi- Static Thermal Stresses in a Solid Sphere

Transcription:

Applied Mthemtics,,, - doi:./m.. Pulished Online My (http://www.scirp.og/jounl/m) Astct Fouie-Bessel Expnsions with Aity Rdil Boundies Muhmmd A. Mushef P. O. Box, Jeddh, Sudi Ai E-mil: mmushef@yhoo.co.uk Received Jnuy, ; evised Feuy, ; ccepted Feuy, Seies expnsion of single vile functions is epesented in Fouie-Bessel fom with unknown coefficients. The poposed seies expnsions e deived fo ity dil oundies in polems of cicul domin. Zeos of the geneted tnscendentl eqution nd the eltionship of othogonlity e employed to find the unknown coefficients. Sevel numeicl nd gphicl exmples e explined nd discussed. Keywods: Fouie-Bessel Anlysis, Boundy Vlue Polems, Othogonlity of Bessel Functions. Intoduction Sevel oundy vlue polems in the pplied sciences e fequently solved y expnsions in cylindicl hmonics with infinite tems. Polems of cicul domin with ounded sufces often genete infinite seies of Bessel functions of the fist nd second types with unknown coefficients. In this cse, the intention is to find the seies coefficients which should stisfy the oundy conditions. The suject of Fouie-Bessel seies expnsions ws investigted nd exmined in mny texts [-]. Nely ll of them hs deived cylindicl hmonics expnsions in J () fo the intevl [, ] only, whee J () is the Bessel function of the fist kind with ode zeo nd gument []. The existence of the oigin point excludes Y (), Bessel function of the second kind with ode zeo nd gument, ecuse it goes to negtive infinity s ppoches zeo []. Both J () nd Y () e shown plotted in Figue. In mny othe polems in the pplied sciences, the intevl of expnsion is found to e [, ] such tht, R. An exmple of this could e hollow cylinde in het conduction polems o cicul nd in vitions nlysis solved in the cylindicl coodinte system. In this cse, cylindicl hmonics expnsions in oth J () nd Y () e necessy. In this ppe, the deivtion of cylindicl hmonics expnsion of single vile function in [, ] in oth J () nd Y () is solved. In ccodnce with the oundies t = nd =, zeos of the otined tnscendentl eqution e fist clculted. As shown in Figue, the solution egion is fo whee the desied seies expnsions e foced to e zeo t = nd = espectively. Unknown coefficients e then found nd the complete seies expnsion cn e chieved. J() nd Y() Figue. Eqution (), J (), Y (). Solution Region Figue. The solution egion in dil oundies. Copyight SciRes.

M. A. MUSHREF. Fomultion nd Solution The Bessel diffeentil eqution of ode zeo is well known s [, ]: d d f ( ) + f ( ) + α f ( ) = () d d α nd R nd. The genel solution to Eqution () fo el vlues of αα is known to e [, ]: f ( ) = n= An J ( α ) + BnY ( α) () As in Eqution (), the ssumed oundy conditions t = nd = e of Diichlet type s f() = nd f() = espectively. Both A n nd B n e then elted s: Y ( α) An = B n () J ( α) A Y ( α) = () n B n J ( α) Going fte the elimintion method, the tnscendentl eqution cn e otined s: J α ) Y ( α) J ( α) Y ( α) () ( = In ode fo Eqution () to e stisfied, thee exist mny zeos o vlues of αα to e clculted. Thus, in ll fome nd coming equtions αα cn e eplced y αα n which e the zeos otined fom the tnscendentl eqution n I. Tht is: J α ) Y ( α ) J ( α ) Y ( α ) () ( n n n n = The othogonlity fetue of Bessel functions cn e pplied to Eqution () y multiplying oth sides y [ Am J( α m) + BmY ( α m) ] nd integting it ove ll possile vlues of fom to s: whee, n ( ) Cm ( ) d = n= C C C ( ) f ( ) d m ( m m m m n ( ) An J ( α n) + BnY ( α n) C m () ) = A J ( α ) + B Y ( α ) () = () The tems unde the summtion in the left side of Eqution () e zeos fo ll vlues of m n [,, ]. Hence, Eqution () cn e simplified to: n ) { Cn ( ) f ( ) } d = C ( () Eithe Eqution () o () cn help. Using Eqution () we cn otin the B n coefficients s: B n = S ( α ) f ( ) d n [ S ( α ) ] n d () whee, S (αα n ) is given y: Y ( αn) S( αn) = Y ( αn) J( αn) () J ( α ) By Eqution () o (), the A n coefficients cn lso e found. Once the coefficients A n nd B n e clculted, the function f() cn e expnded s in Eqution ().. Numeicl Exmples The tnscendentl expession in Eqution () shows gdul decy s αα inceses which men smll mgnitudes etween high zeos. This leds to the convegence of the seies in Eqution () ove s n inceses. As consequence, finite nume of tems in Eqution () cn e sufficient fo numeicl ppoximtions. The zeos e fist evluted using the tnscendentl coss poduct Bessel functions eqution fo the intevl [, ]. A gph of Eqution () is shown in Figue fo the solution egions [.,.] nd [., ]. Tle shows the fist zeos of Eqution () fo =. nd =.. Zeos otined fom the tnscendentl eqution chnges ccoding to the vlues of nd ssumed fo the solution egion. The dt pesented in Tle indictes tht the clculted zeos e not peiodic nd should e clculted using pope numeicl technique. Let s ssume tht the function f() to e expnded s in Eqution () is sin() with dil solution egion in [.,.]. The coefficients B n cn e evluted fom Eqution () nd the A n coefficients e then otined y Eqution (). Both coefficients e shown in Tles nd espectively fo n = to. J(αn)Y(αn) J(αn)Y(αn) Figue. Eqution (), [.,.], [., ]. α n Copyight SciRes.

M. A. MUSHREF Mny vitions cn e noticed fo the numeicl vlues of A n nd B n with genel solute scle of < except fo B =.. Some coefficients e in the ode of - mening tht thei ssocited tems e vey smll such s B nd A in Tles nd espectively. The function sin() nd its ppoximte expnsions e plotted in Figue. Summtion ove the fist tems poduced n cceptle estimtion in the intevl [.,.] with some ppent oscilltions ound the exct function. An impoved ppoximte expnsion is lso plotted fo n = to with less fluctutions in the sme dil domin. Tle. Fist fifty zeos of Eqution () in [.,.]. n α n n α n n α n n α n n α n.................................................. Tle. Fist fifty B n fo f() = sin() in [.,.]. n B n n B n n B n n B n n B n.....e-.................e-....................e-........ Tle. Fist fifty A n fo f() = sin() in [.,.].................................E-.................. Tle. Fist fifty B n fo f() = cos() in [.,.]. n B n n B n n B n n B n n B n.....e-.....................................e-........ In ddition, f() = cos() is expnded s in Eqution () nd the fist fifty coefficients e listed in Tles nd fo the B n nd A n espectively. Simil to the sin(), the cos() coefficients go though sevel vitions with genel solute scle of < except A =.. Also, only fou coefficients e in the ode of - implying tht thei elted tems in the seies e extemely smll such s B nd A in Tles nd espectively. Copyight SciRes.

M. A. MUSHREF Tle. Fist fifty A n fo f() = cos() in [.,.]................E-.................E-.................. sin() Figue. sin(), Eqution () with n = to, Eqution () with n = to. cos() with n solute level of > o <. The getest vlues in Tles nd e found s B =. nd A =.. In ddition, no coefficients e clculted in the ode of - implying tht ll coefficients e to e included in the seies expnsion. The function exp() nd its estimted expnsions e shown plotted in Figue in [.,.]. A stisfctoy estimtion of finite summtion ove the fist tems e geneted with sevel oscilltions close to the exct function. A good ppoximted expnsion is lso plotted fo n = to with fewe vitions in the sme solution egion. The lst numeicl exmple to e discussed is the sque function expessed s:.. f ( ) = () othewise The clculted B n nd A n coefficients fo this function e shown in Tles nd espectively. Simil to fome expnsions, oth coefficients vy out the exct vlues of Eqution (). The B n coefficients hve genel solute level of < except B, B, B, B nd B tht hve n solute scle of >. Futhemoe, the A n coefficients show n solute level of < except the solute vlues of A, A, A nd A tht e >. Some B n nd A n coefficients e clculted in the ode of - like A o Figue. cos(), Eqution () with n = to, Eqution () with n = to. The function cos() nd its estimted expnsions e shown plotted in Figue. Finite summtion ove the fist tems geneted stisfctoy estimtion in the intevl [.,.] with sevel ovious oscilltions close to the exct function. A ette ppoximte expnsion is lso plotted fo n = to with less fluctutions in the sme solution egion. The clculted coefficients fo the function e e lso shown in Tles nd fo B n nd A n espectively. Appently, the coefficients swing ound the exct vlues exp() Figue. exp(), Eqution () with n = to, Eqution () with n = to. Copyight SciRes.

M. A. MUSHREF Tle. Fist fifty B n fo f() = exp() in [.,.]. n B n N B n n B n n B n n B n.................................................. Tle. Fist fifty A n fo f() = exp() in [.,.]................................................... Tle. Fist fifty B n fo Eqution () in [.,.]. n B n n B n n B n n B n n B n.... E-......E-......E-.... E-.. E-.. E-..... E-. E-..... E-..... E-.. Tle. Fist fifty A n fo Eqution () in [.,.]. E-..........E-.E-.................E-......E-.............E- E-. in the ode of - such s B indicting tht thei ssocited tems in the seies e vey smll. The function expessed y Eqution () nd its ppoximte expnsions e plotted in Figue. Summtion ove the fist tems poduced n cceptle estimtion in the intevl [.,.] with some noticele oscilltions ound the exct function. A ette ppoximte expnsion is lso plotted fo n = to with less fluctutions in the sme dil domin. In ll gphicl plots peviously shown, the cuves etun to zeo t the ssumed oundies =. nd =.. In ddition, ccucy of the expnded cuves my ppe ette s n inceses due to lge nume of tems involved in the seies nd less fluctutions seen ound the exct vlues.. Conclusions Functions wee expnded s Fouie-Bessel seies summtion in oth J () nd Y (). A finite seies expn- Copyight SciRes.

M. A. MUSHREF f() Figue. Eqution (), Eqution () with n = to, Eqution () with n = to. sion ws otined fo ity dil oundies in [, ]. Coefficients wee found y clculting the zeos of the tnscendentl eqution nd y employing the eltionship of othogonlity. A nume of exmples wee numeiclly nd gphiclly discussed.. Refeences tions, Cmidge Univesity Pess, United Kingdom,. [] H. F. Dvis, Fouie Seies nd Othogonl Functions, Dove Pulictions, Inc., New Yok,. [] W. E. Byely, An Elementy Tetise on Fouie s Seies, Dove Pulictions, Inc., New Yok,. [] F. Bowmn, Intoduction to Bessel Functions, Dove Pulictions, Inc., New Yok,. [] N. N. Leedev, Specil Functions nd Thei Applictions, Dove Pulictions, Inc., New Yok,. [] E. A. Gonzlez Velsco, Fouie Anlysis nd Boundy Vlue Polems, Acdemic Pess, Sn Diego,. [] A. Bomn, Intoduction to Ptil Diffeentil Equtions: Fom Fouie Seies to Boundy Vlue Polems, Dove Pulictions, Inc., New Yok,. [] P. V. O Neil, Advnced Engineeing Mthemtics, Wdswoth Pulishing Compny, Clifoni,. [] H. Sgn, Boundy nd Eigenvlue Polems in Mthemticl Physics, Dove Pulictions, Inc., New Yok,. [] N. N. Leedev, I. P. Sklsky nd Y. S. Uflynd, Woked Polems in Applied Mthemtics, Dove Pulictions, Inc., New Yok,. [] G. N. Wtson, A Tetise on the Theoy of Bessel Func- Copyight SciRes.