A Survey of Azimuthal Angle and Eigenvalues of the Laplace Equation

Similar documents
THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

On Polynomials Construction

15 Solving the Laplace equation by Fourier method

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

SOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4

DonnishJournals

8 Separation of Variables in Other Coordinate Systems

On the Quasi-inverse of a Non-square Matrix: An Infinite Solution

Measure Estimates of Nodal Sets of Polyharmonic Functions

POISSON S EQUATION 2 V 0

Do not turn over until you are told to do so by the Invigilator.

Green s Identities and Green s Functions

B da = 0. Q E da = ε. E da = E dv

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4

FREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS


SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES

arxiv: v1 [physics.pop-ph] 3 Jun 2013

=0, (x, y) Ω (10.1) Depending on the nature of these boundary conditions, forced, natural or mixed type, the elliptic problems are classified as

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

Liquid gas interface under hydrostatic pressure

1D2G - Numerical solution of the neutron diffusion equation

The Divergence Theorem

FOURIER-BESSEL SERIES AND BOUNDARY VALUE PROBLEMS IN CYLINDRICAL COORDINATES

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

Question 1: The dipole

Math 2263 Solutions for Spring 2003 Final Exam

The first nontrivial curve in the fučĺk spectrum of the dirichlet laplacian on the ball consists of nonradial eigenvalues

Vectors, Vector Calculus, and Coordinate Systems

Exceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs

I. CONSTRUCTION OF THE GREEN S FUNCTION

Numerical Integration

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2

2.25 Advanced Fluid Mechanics

Lecture 7: Angular Momentum, Hydrogen Atom

Physics 506 Winter 2006 Homework Assignment #9 Solutions

1 Spherical multipole moments

Dymore User s Manual Two- and three dimensional dynamic inflow models

Question Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

Chapter 3 Optical Systems with Annular Pupils

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

An Exact Solution of Navier Stokes Equation

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

Doubling property for the Laplacian and its applications (Course Chengdu 2007)

EN40: Dynamics and Vibrations. Midterm Examination Thursday March

CHAPTER 25 ELECTRIC POTENTIAL

Stress, Cauchy s equation and the Navier-Stokes equations

Math 124B February 02, 2012

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

A Lower Bound for the First Steklov Eigenvalue on a Domain

Regularity for Fully Nonlinear Elliptic Equations with Neumann Boundary Data

A Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation

2 Governing Equations

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

Moment-free numerical approximation of highly oscillatory integrals with stationary points

Scattering in Three Dimensions

A Method of Solving Compressible Navier Stokes Equations in Cylindrical Coordinates Using Geometric Algebra

AXIS-SYMMETRIC FRACTIONAL DIFFUSION-WAVE PROBLEM: PART I-ANALYSIS

Nuclear size corrections to the energy levels of single-electron atoms

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Identification of the degradation of railway ballast under a concrete sleeper

Supplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

Classical Mechanics Homework set 7, due Nov 8th: Solutions

COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT

EM Boundary Value Problems

Right-handed screw dislocation in an isotropic solid

Precessing Ball Solitons as Self-Organizing Systems during a Phase Transition in a Ferromagnet

Gradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function

Lecture 23. Representation of the Dirac delta function in other coordinate systems

Multiple Experts with Binary Features

PES 3950/PHYS 6950: Homework Assignment 6

In the previous section we considered problems where the

On the Poisson Approximation to the Negative Hypergeometric Distribution

Results on the Commutative Neutrix Convolution Product Involving the Logarithmic Integral li(

Review: Electrostatics and Magnetostatics

New problems in universal algebraic geometry illustrated by boolean equations

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

A new class of exact solutions of the Navier Stokes equations for swirling flows in porous and rotating pipes

Homework Set 3 Physics 319 Classical Mechanics

Lecture 10. Vertical coordinates General vertical coordinate

Numerical Methods for Two-Fluid Dispersive Fast MHD Phenomena

A STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction

Surveillance Points in High Dimensional Spaces

Multipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source

On the global uniform asymptotic stability of time-varying dynamical systems

arxiv: v1 [physics.gen-ph] 18 Aug 2018

Chapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in

A Power Method for Computing Square Roots of Complex Matrices

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

Safety variations in steel designed using Eurocode 3

(n 1)n(n + 1)(n + 2) + 1 = (n 1)(n + 2)n(n + 1) + 1 = ( (n 2 + n 1) 1 )( (n 2 + n 1) + 1 ) + 1 = (n 2 + n 1) 2.

Transcription:

Contempoay Engineeing Sciences, Vol., 08, no. 95, 4743-4749 HIKAI Ltd, www.m-hikai.com https://doi.og/0.988/ces.08.8950 A Suvey of Azimuthal Angle Eigenvalues of the Laplace Equation Luz aía ojas Duque eseach in Statistics Epidemiology GIEE) Fundación Univesitaia del Áea Andina, Peeia, Colombia Peeia, Colombia Pedo Pablo Cádenas Alzate Depatment of athematics GEDNOL Univesidad Tecnológica de Peeia Peeia, Colombia José Geado Cadona Too Depatment of athematics GIEE Univesidad Tecnológica de Peeia Peeia, Colombia Copyight c 08 Luz aía ojas Duque et al. This aticle is distibuted unde the Ceative Commons Attibution License, which pemits unesticted use, distibution, epoduction in any medium, povided the oiginal wok is popely cited. Abstact In this wok we show that the eigenvalues of the Laplacian with homogeneous Neumann B.C. cannot be positive we solve the Diichlet BVP fo the Laplace equation u = 0 in the egion between two concentic sphees = = ). Finally, we show the simplifications in the solution if Diichlet bounday data ae independent of the Azimuthal angle. Keywods: Laplace equation, Azimuthal angle, spheical hamonics

4744 Luz aía ojas Duque et al. Intoduction Let us denote by ω = ω, ω,..., ω n ) the unit nomal vecto outwad) to the bounday. Note that ω is a vecto field that is defined at each point of the bounday. That u is an eigenfunction of the Laplacian with the homogeneous Neumann condition means that thee is a numbe λ called the eigenvalue coesponding to u, such that [] u = λu, ) in ω u=0 on, whee ω u is the nomal deivative ω u = ω u x ω n u xn We equie u to be not identically zeo, since the zeo function u = 0 satisfies Eq.) fo any λ. ecall Geen s fist identity ω u = ω ω u ω u, which can be poven by applying the divegence theoem to the vecto field ω u. Now, multiplying by u integating ove we get u u = λ u. Theefoe, by Geen s identity we obtain u ω u u = λ which implies that λ u = u, ) then due to the homogeneous Neumann condition, the bounday integal tem is null. Now, si λ 0, then u must be nonzeo in ode to be an eigenfunction []. Diichlet case Fo the Diichlet case, by a simila easoning it is concluded that λ 0. Now, if thee is an eigenfunction u with λ = 0, then in u = 0 on. u = 0, u,

Laplace equation 4745 Hee, the identically-zeo function u satisfies this equation it is the only solution). Finally, thee is no nontivial function u = that satisfies u = 0, concluding that all eigenvalues ae nonzeo, i.e., λ < 0 [3].. Diichlet B.V.P fo the Laplace equation Hee, employing spheical coodinates, the equation we need to solve is [4] { u = 0, {, θ, φ) : < < } u, θ, φ) = fθ, φ), u, θ, φ) = gθ, φ). The 3-dimensional Laplacian in spheical coodinates is u = u ) u θ sin θ) θ sin θ u φφ sin θ, o u = u u u θθ cot θ u θ sin θ u φφ. Now, substituting u, θ, φ) = )Nθ, φ) into u = 0 we obtain If, we multiply by N N N N θθ cot θ N θ sin θ N φφ we get = N θθ N N θ cot θ N φφ N N sin θ = λ Hee, the equation fo the spheical hamonics, we have the solution fo n = 0,,... m = ±,..., ±n 3) with λ = n )n. Nθ, φ) = Y n,m θ, φ), whee Theefoe, the spheical hamonics ae explicitly given by Y n,m θ, φ) = sin m θ)p n m ) cos θ)ϕ m φ), ϕ m φ) = { cosmφ), m 0 sinmφ), m < 0. Hee, P n is the Legende polynomial of degee n [6]. Taking into account that λ = n )n, the equation fo is = nn ).

u = 0, in {,, ') : < < }, u,, ') = f, '), u,, ') = g, '). ecall that the 3-dimensional Laplacian in spheical coodinates is u'' u ) u sin ) sin sin Luz aı a cot Duque et al. ojas = u u u u u''. sin 7 Let us look fo a solution that can be witten as u,, ') = ), '). Substituting this y 00 into the equation u = 0, we get u= 4746 cot 00 0 '' = 0, sin multiplying by y 0 y eaanging, we have y '' 00 0 cot = =. SOLUTIONS TO SELECTED POBLES FO ASSIGNENTS 3, 4 sin A vaiation of Poblem fom Assignment 4 We 4ecognise the equation given by the second equality sign as the equation fo the spheical hamonics, so the fo Laplace each = 0,,y..u.,= m = 0, ±,..., ±n, we have the solution y value y fo y 0 nequation Statement. Solve the Diichlet bounday poblem 0 infoy each the egion between two concentic sphees of adii. What simplifications do we get = Yn,mangle, '), in the solution pocess if the Diichlet bounday data ae independent of, the') azimuthal longitude) '? with ecall that the spheical hamonics ae explicitly given by Solution. Employing spheical coodinates, the equation we need to solve is u = 0, in {,, ') : < Figue B.: Plots of the eal-valued spheical hamonic basis functions. Geen indicates positive values m m ) values., ') = f, '), ed<indicates }, negative n,m u, u,, ') = g, ').m Y, ') = sin ) Pn cos ) '), whee FigueLaplacian : eal-valued spheical hamonics basis functions. ecall that the 3-dimensional in spheical coodinates is = nn ). cosm ) fo m 0, m ') = sinm ) fo m < 0, only the fequencies of the function up to some theshold, obtaining an n th ode b-limited Pn is the Legende polynomial of degee n. u'' function u )appoximation u sin ) whee f: u= f of the oiginal sin sin cot n l = u u u u f ~!) = X ux ''.y m ~!) f lm. B.) sin l =0 m= l l Let us look fo a solution that can be witten as u,, ') = ), '). Substituting this into the equation u = 0, we getlow-fequency functions can be well appoximated using only a few bs, as the numbe of coefficients inceases, cot highe fequency signals can be appoximated moe accuately. 00 0 '' = 0, It is often convenient toefomulate sin the indexing scheme to use a single paamete multiplying by i = l l ) m. have With this convension it is easy to see that an n th ode appoximation can be eaanging, we econstucted coefficients, 0 '' using n cot = =. sin nx!) = equation y i ~!) f i. fo the spheical We ecognise the equation given by the second equality sign asf ~the i =0 hamonics, so fo each n = 0,,..., fo each m = 0, ±,..., ±n, we have the solution 00 B.), ') = Yn,m, '), a) with ). n =,,=3, 4nn fom top to bottom, b) n = 3, 4, 5, 6 fom top to bottom, m =,..., n m =given 0,... by, n fom left to ight fom left to ight ecall that the spheical hamonics ae explicitly n,m cosm ) fo m 0, Yn,m, ') = sin m ) Pn m ) cos ) m '), whee m ') = sinm )Figue fo m <. 0, Spheical hamonics Ynm Figue : eal-valued spheical hamonics Y 0,..., n. fo n =,, 3, 4 m = whee Pn is the Legende polynomial of degee n. a) n =,, 3, 4 fom top to bottom, m = 0,..., n fom left to ight b) n = 3, 4, 5, 6 fom top to bottom, m =,..., n fom left to ight Figue 3: eal-valued spheical hamonics Yn,m fo n = 3, 4, 5, 6 m = Figue. Spheical hamonics Ynm,..., n. Hee, taking ) = α gives αα ) α = nn ),

Laplace equation 4747 which has the solutions α = n α = n ). We conclude that fo each n N 0, fo each m Z, we have two independent solutions u, θ, φ) = n Y n,m θ, φ) u, θ, φ) = n) Y n,m θ, φ), so the geneal solution is given by a linea combination of all u, θ, φ) = n n=0 m= n A n,m n B ) n,m Y n n,m θ, φ) The coefficients A n,m B n,m ae to be found the bounday conditions. Then, we have, n fθ, φ) = fn,my n,m θ, φ) n=0 m= n n gθ, φ) = gn,my n,m θ, φ), n=0 m= n Finally A n,m B n,m A n,m n B n,m n = g n,m A n,m = n g n,m f m,n n B n,m = n fn,m n gm,n. n Now, when the bounday data ae independent of φ. The spheical hamonics Y n,m that do not depend on φ ae the ones with m = 0, i.e, the zonal hamonics hence the solution is simply u, θ, φ) = Y n,0 θ, φ) = Z n θ) = P n cos θ), n=0 A n n B ) n P n n cos θ).

4748 Luz aía ojas Duque et al. Theefoe, so we look fo a solution u, θ) = )ϕθ) athe than u, θ, φ) = )Nθ, φ) we get = ϕ ϕ ϕ cot θ ϕ = λ. We ecognise the equation fo ϕ as the equation leading to the Legende equation x )y xy = λy, whee yx) = ϕθ) with x = cos θ, then we have a solution ϕθ) = P n cos θ) with λ = nn ) fo each n = 0,,... [7]. 3 Conclusion In this pape we demonstate that the eigen-values of the Laplacian on with homogeneous Neumann bounday conditions cannot be positive. Also, we have poven that the eigenvalues of the Laplacian on with homogeneous Diichlet bounday condition) ae stictly negative. Theefoe, fo the Diichlet case, the Laplace equation has a unique solution with the homogeneous Diichlet bounday condition. We have solve the Diichlet bounday value poblem fo the Laplace equation u = 0 in the egion between two concentic sphees of = =. Acknowledgements. We would like to thank the efeee fo his valuable suggestions that impoved the pesentation of this pape ou gatitude to the Depatment of athematics of the Univesidad Tecnológica de Peeia Colombia), the goup GEDNOL eseach in Statistics Epidemiology GIEE). efeences [] P. Cádenas Alzate, A suvey of the implementation of numeical schemes fo linea advection equation,advances in Pue athematics, 4 04), no. 8, 467-479. https://doi.og/0.436/apm.04.4805 [] O. Agawal, Genealized Eule-Lagange equations tansvesality conditions fo FVPs in tems of the Caputo deivative, Jounal of Vibation Contol., 3 007), 7-37. https://doi.og/0.77/07754630707747

Laplace equation 4749 [3] Dinh Nho Hao, B.T. Johansson, D. Lesnic, Pham inh Hien, A vaiational method appoximations of a Cauchy poblem fo elliptic equations, J. Algoithms Comput. Technol., 4 00), 89-9. https://doi.og/0.60/748-308.4..89 [4]. Fontelos, Fundamentos atemáticos de la Ingenieía, Libeía- Editoial Dykinson, Vol. 0, 007. [5] K. Atkinson, The Numeical Solution of Integal Equations of the Second Kind, Cambidge: Cambidge Univesity Pess, 997. https://doi.og/0.07/cbo9780566340 [6] C. uñoz, Pedo Pablo Cadenas Alzate, Anibal unoz Loaiza, An iteative method fo solving two special cases of nonlinea PDEs, Contempoay Egineeing Sciences, 0 07), no., 545-553. https://doi.og/0.988/ces.07.765 [7] P. Cádenas Alzate, Jose Geado Cadona, Luz aia ojas A special case on the stability accuacy fo the D heat equation using 3-Level θ-schemes, Applied athematics, 6 05), no. 3, 476-483. https://doi.og/0.436/am.05.63045 eceived: Octobe 7, 08; Published: Novembe, 08