Semi-Lagrangian simulations of the diocotron instability
|
|
- Eunice Horn
- 6 years ago
- Views:
Transcription
1 Semi-Lagangian simulations of the diocoton instability Eic Madaule, Seve Adian Histoaga, Michel Mehenbege, Jéôme Péti To cite this vesion: Eic Madaule, Seve Adian Histoaga, Michel Mehenbege, Jéôme Péti. Semi-Lagangian simulations of the diocoton instability. [Reseach Repot] 213. <hal-84154> HAL Id: hal Submitted on 5 Jul 213 HAL is a multi-disciplinay open access achive fo the deposit and dissemination of scientific eseach documents, whethe they ae published o not. The documents may come fom teaching and eseach institutions in Fance o aboad, o fom public o pivate eseach centes. L achive ouvete pluidisciplinaie HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau echeche, publiés ou non, émanant des établissements d enseignement et de echeche fançais ou étanges, des laboatoies publics ou pivés.
2 Semi-Lagangian simulations of the diocoton instability Seve Histoaga Eic Madaule Michel Mehenbege Jéôme Péti July 3, 213 Abstact We conside a guiding cente simulation on an annulus, following [2, 3]. Thee, the PIC method was used. We popose hee to evisit this test case by using a classical semi-lagangian appoach [4]. Fist, we obtain the consevation of the electic enegy and mass fo some adapted bounday conditions. Then we ecall the dispesion elation fom [1] and discussions on diffeent bounday conditions ae detailed. Finally, the semi-lagangian code is validated in the linea phase against analytical gowth ates given by the dispesion elation. Also we have validated numeically the consevation of electic enegy and mass. Numeical issues/difficulties due to the change of geomety can be tackled in such a test case which thus may be viewed as a fist intemediate step between a classical guiding cente simulation in a 2D catesian mesh and a slab 4D dift kinetic simulation. 1 The guiding cente model We conside the guiding cente model in pola coodinates fo ρ = ρ(t,, θ) t ρ θφ ρ + Φ θρ =, (1.1) coupled with the Poisson equation fo the potential Φ = Φ(t,, θ) 2 Φ 1 Φ θφ = γρ, whee γ is a numbe which is fixed to eithe 1 o 1. Geneally, we take γ = 1 (see [1]); taking γ = 1 can be intepeted as looking fo solutions backwad in time. The domain is (, θ) = [ min, max ] [, 2π]. The initial and bounday conditions will be discussed late on. We assume fo the moment that we have peiodic bounday conditions in θ. Poposition 1.1. We define the electic enegy E(t) = Φ θφ 2 ddθ, (1.2) 1
3 and the mass We have M(t) = ρddθ. (1.3) (i) (ii) (iii) t E(t) = 2 E(t) = [Φ Φ] =max dθ + γ ρφddθ. [Φ t Φ] =max dθ + 2γ t M(t) = Poof. Note that the Poisson equation wites [Φ θ ρ] =max dθ [Φρ θ Φ] =max dθ. ( Φ) 1 2 θφ = γρ. (1.4) We have, by integating by pats the electic enegy (1.2), in fo the fist tem and in θ fo the second tem, 1 E(t) = [Φ Φ] =max dθ Φ ( Φ)ddθ Φ 2 θφ ddθ, and get (i), by using the Poisson equation (1.4). We then have t E(t) = I 1 + γ(i 2 + I 3 ), (1.5) with I 1 = [ t (Φ Φ)] =max dθ, I 2 = Φ t ρddθ, I 3 = ρ t Φddθ. By using the guiding cente equation (1.1) and integating as above (in fo the fist tem and in θ fo the second tem), we compute I 2 = Φ( θ Φ ρ Φ θ ρ)ddθ = [Φρ θ Φ] =max dθ (Φ θ Φ)ρddθ + θ (Φ Φ)ρddθ = [Φρ θ Φ] =max dθ, since we have the elation (Φ θ Φ) = θ (Φ Φ). Using the Poisson equation (1.4), and integating as usually, we have γi 3 = t Φ ( ( Φ) θφ ) ddθ = [ Φ t Φ] =max dθ + I 4, 2
4 whee I 4 = ( Φ t Φ ) ddθ + t θ Φ 1 θφddθ. Integating as usually and using then the Poisson equation (1.4), we get ( I 4 = [Φ t Φ ) ( ] = =max min dθ Φ t Φ ) ddθ Φ t θφ 2 1 ddθ = [Φ t ( Φ ) ] =max dθ + γi 2. Replacing the fomulae fo I 1, I 2, I 3 in (1.5) gives t E(t) = [ t (Φ Φ)] =max dθ + 2γ [ Φ t Φ] =max dθ + [Φρ θ Φ] =max dθ [Φ t ( Φ ) ] =max dθ which is (ii), using the elation t (Φ Φ) = t Φ Φ + Φ t Φ. It emains to pove the fomula (iii) fo the mass. We use the guiding cente equation (1.1) and integate as usually t M(t) = t ρddθ = = = ( θ Φ ρ + Φ θ ρ) ddθ [Φ θ ρ] =max dθ + (Φ θ ρ Φ θ ρ) ddθ [Φ θ ρ] =max dθ. Poposition 1.2. We suppose Diichlet bounday conditions at min and at max : Φ(t, min, θ) = Φ(t, max, θ) =. Then the electic enegy and the mass ae constant in time: t E(t) =, t M(t) =. Poof. We get immediately the esult fom Poposition 1.1. Poposition 1.3. We suppose the following bounday conditions: (i) Diichlet bounday condition at max Φ(t, max, θ) = (ii) Inhomogeneous Neumann bounday condition at min fo the Fouie mode in θ: whee Q is a given constant. Φ(t, min, θ)dθ = Q, 3
5 (iii) Diichlet bounday condition at min fo the othe modes, which eads Φ(t, min, θ) = 1 Φ(t, min, θ )dθ 2π Then, the electic enegy and the mass ae constant in time: t E(t) =, t M(t) =. Poof. Fom Poposition 1.2 and using fist the Diichlet bounday condition at max (i), we get t E(t) = 2 Φ(t, min, θ) ( min t Φ(t, min, θ) + γρ(t, min, θ) θ Φ(t, min, θ)) dθ. Note that, by (iii), Φ(t, min, θ) does not depend on θ and thus, we get Φ(t, min, θ) = Φ(t, min, ) and θ Φ(t, min, θ) = which leads to We thus get, fom (ii) t E(t) = 2Φ(t, min, ) min t Φ(t, min, θ)dθ. t E(t) = 2Φ(t, min, ) min t Φ(t, min, θ)dθ = 2Φ(t, min, ) min t Q =. Fo the mass, we get similaly t M(t) = Φ(t, min, ) θ ρ(t, min, θ)dθ =. Remak 1. We do not know whethe the electic enegy and the mass emain constant in time when we conside the following bounday conditions: (i) Diichlet at max : Φ(t, max, θ) =. (ii) Neumann at min : Φ(t, min, θ) =. If we impose that ρ(t, min, ) =, we can check that we get also consevation of mass and electic enegy. In the semi-lagangian code, we have not imposed this condition (we have tied such a condition but got wose consevation fo such bounday condition on Φ); if the foot of the chaacteistic is out of the domain at a point (, θ ), we take the value at ( min, θ ), if < min o at ( max, θ ), if > max. 4
6 2 Study of the electostatic eigenvalue equation and discetization The eigenvalue equation which pemits to obtain the gowth ates of instability, has been deived in Davidson [1] (equation (6.28), p 296). We ecall this deivation, in ou notations. We begin with the following poblem t ρ θφ ρ + Φ θρ = 2 2 Φ Φ 2 θ 2 Φ 2 = γρ (2.6) We conside at fist an equilibium n () that does not depend on θ and coesponding potential Φ (), which satisfies ( Φ ()) = n (). As an equilibium, the couple (ρ, Φ) = (n, γφ ) satisfies (2.6). Then, we petub this equilibium and seach an appoximate solution of (2.6) in the fom We set also ρ(t,, θ) n () + εn 1 (t,, θ) Φ(t,, θ) γφ () + εγφ 1 (t,, θ) We inject these appoximations in (2.6) to get γ t n 1 θφ 1 n + Φ θ n 1 = O(ε) 2 θ Φ 2 1 = n O(ε) 2 2 Φ 1 Φ 1 We now look fo a paticula solution of the lineaized poblem in the fom γ t n 1 θφ 1 n + Φ θ n 1 = 2 2 Φ 1 Φ 1 2 θ 2 Φ 1 2 = n 1 n 1 (t,, θ) = n 1,l () exp(ilθ) exp( iωt/γ), Φ 1 (t,, θ) = Φ 1,l () exp(ilθ) exp( iωt/γ). This type of solution can be obtained though Laplace tansfom in time and Fouie tansfom in θ. We look fo unstable equilibium, that means that we want to have a solution with I(ω)/γ >. We obtain iω n 1,l il n 2 2 Φ 1,l Φ 1,l Φ 1,l + il Φ n 1,l = + l2 2 Φ 1,l = n 1,l 5
7 We expess n 1,l in tems of Φ 1,l to get an equation on Φ 1,l : ( iω + il Φ )( 2 Φ 2 1,l Φ 1,l Φ 2 1,l ) = il n Φ 1,l, (2.7) whee i = min + i and = ( max min )/N and N N. We can poceed in the numeical esolution of the poblem by making the following discetization: Φ 1,l ( i ) = φ i, Φ1,l ( i ) = φ i+1 φ i 1, 2 2 Φ 2 1,l ( i ) = φ i+1 2φ i + φ i 1 (2.8) 2 The poblem can then be witten as follows Aφ = cbφ. We then have to find the eigenvalues c of the poblem B 1 Aφ = cφ, whee c = ω. l We look fo an eigenvalue c such that c/γ has the geatest imaginay pat which should be stictly positive. 3 Diocoton instability fo an annula electon laye + l2 We conside the following initial data, min <, ρ(,, θ) = 1 + ε cos(lθ), +,, + < max, whee ε is a small paamete. In that specific case, we can do analytical computations. The linea analysis is pefomed in [1]. We thus can find an explicit solution of (2.7): Φ 1,I (), min <, Φ 1,l () = Φ 1,II (), < +, Φ 1,III (), + < max, whee Φ 1,I () = ( /) l (B( ) l + C( ) l ) (/ min) 2l + ɛ ( / min ) 2l + ɛ, Φ 1,II () = B l + C l, Φ 1,III () = ( + /) l (B( + ) l + C( + ) l ) 1 (/ max) 2l 1 ( + / max ) 2l. Hee, we have consideed Diichlet bounday condition at max. At min, we have Diichlet bounday condition, if ɛ = 1 and Neumann bounday condition if ɛ = 1. The constants B and C ae chosen so that (see [1]) and ( Φ 1,II ( ) Φ 1,I ( Φ 1,I ( ) )) + l ω lω m ( ) =, 6
8 with and + ( Φ 1,III ( + ) Φ 1,II ( + Φ 1,II ( + ) )) + l ω lω m ( + ) =, ( ω m () = 1 ( ) ) 2 ( ) ω q, 2 ω q = 2 Q ( ). 2 This leads to a 2 2 linea system fo (B, C). We look fo a solution with (B, C) (, ). This implies that the deteminant of the coesponding matix is zeo. This is the dispesion elation. In the case of bounday conditions given in Poposition 1.3, the dispesion elation can be explicitely given (see [1]) ω 2 b l ω + c l =, whee and and b l = l(1 ( / + ) 2 + ω q (1 + ( / + ) 2 ))(1 ( min / max ) 2l ) + (1 ( / + ) 2l ) (+ / max ) 2l ( min / ) 2l 1 ( min / max ) 2l (3.9) c l (1 ( min / max ) 2l ) = l 2 ω q (1 ( / + ) 2 + (ω q )( / + ) 2 )(1 ( min / max ) 2l ) lω q (1 ( min / + ) 2l )(1 ( + / max ) 2l ) + l(1 ( / + ) 2 + ω q ( / + ) 2 )(1 ( / max ) 2l )(1 ( min / ) 2l ) (1 ( + / max ) 2l )(1 ( min / ) 2l )(1 ( / + ) 2l ) (3.1) ω q = 4 Q ( ). 2 The case of the bounday conditions of Poposition 1.2 is obtained by taking in the dispesion elation. Q = ( ) 2 ( + ) 2 + 2( ) 2 ln( max / ) + 2( + ) 2 ln( + / max ), 8 ln( min / max ) Remak 2. In [1], bounday conditions of Poposition 1.3 wee consideed. Hee, we can also conside bounday conditions of Poposition 1.2 and Neumann bounday conditions (using ɛ = 1). The computations have been done using Maple. Remak 3. The system (2.6) fo γ = 1 o γ = 1 leads to the same gowth ate in the linea phase, since the dispesion elation has eal coefficients. 7
9 Remak 4. By using the linea appoximation ρ = n + εn 1 and Φ = γ(φ + εφ 1 ), which is supposed to be valid in the linea phase, we obtain fo bounday conditions of Remak 1 (Neumann fo all the modes at min ) t E(t) 2γ 3 Φ ( min )ε 2 n 1 (t, min, θ) θ Φ 1 (t, min, θ)dθ + O(ε 3 ), which leads a pioi to a loss of enegy consevation fo this linea appoximation, which may even gow exponentially in time, with a gowth ate dictated by the dispesion elation! On the othe hand, we can check that the enegy emains constant fo the linea appoximation when using bounday conditions of Popositions 1.2 and 1.3, as we have poven it fo the initial non linea poblem. This may explain the linea gowth ate obtained by numeical appoximations fo the electic enegy in [3]. 4 Numeical esults 4.1 Analytical and discetized dispesion elation We obtain the following esults in Tables 1, 2, 3 fo the analytical dispesion elation. We take min = 1, max = 1, γ = 1. The esults ae coheent with those given in [2]. Remak 5. We have taken the same modes as in [2]. Fo mode 5 of Table 3, we had to change + in ode to get an unstable mode. We have also tested a discetization of the dispesion elation (2.8) with Diichlet bounday conditions. The esults ae shown on Table 4. The esults ae coheent with espect to the analytical case. 4.2 Semi-Lagangian method We use a classical backwad semi-lagangian method [4]. A pedicto coecto scheme is used as time scheme, togethe with a symplectic Velet algoithm fo the computation of the chaacteistics, The intepolation is done with cubic splines. Fo the Poisson equation, finite diffeence of ode two ae used in the adial diection and Fouie in the θ diection. We denote by N (esp. N θ ) the numbe of intevals in the adial (esp. θ) diection. The time step is t. We give the gowth ate obtained by the semi-lagangian method in Tables 5, 6, 7. We obtain esults in ageement with those obtained with the dispesion elation, except fo the mode 5, when Neumann bounday conditions ae used fo the mode (o fo all the modes). Fo this specific case, the obtained value is stongly dependent on the discetization paametes. By taking N = 256, N θ = 512 and t =.25, we get I(ω) =.51 on the inteval [123, 241], which is close to the theoetical value. An example of gowth ate and of density function is given on Figue 1. The time evolution of electic enegy and of elative mass eo is given on Figues 2,3,4,5, fo diffeent bounday conditions at min (we suppose Diichlet bounday condition at max ). 8
10 We take hee ε =.5 and l = 3. We obseve convegence of enegy and mass consevation fo Diichlet and Neumann mode bounday conditions at min. On the othe hand, fo Neumann bounday conditions at min, these quantities ae no moe conseved in long time. We had to take a quite lage amplitude of ε to see this phenomenon. 5 Conclusion We have studied the consevation of electic enegy and mass of the continuous model; we have enlightened adhoc bounday conditions fo having such consevations. Then, we have detailed the dispesion elations fo such bounday conditions. Finally, we have given some numeical esults with a classical semi-lagangian method and validated the linea phase and the consevation of mass and electic enegy. A compaison on numeical methods (PIC vs. semi-lagangian method) with egads to consevative quantities fo example, as well as highe dimensional poblem with velocity discetization ae envisaged. Refeences [1] R. C. Davidson, Physics of non neutal plasmas, 199 [2] J. Péti, The diocoton instability in a pulsa cylindical electosphee, Astonomy & Astophysics, Febuay 5, 28. [3] J. Péti, Non-linea evolution of the diocoton instability in a pulsa electosphee: 2D PIC simulations, Astonomy & Astophysics, May 7, 29. [4] E. Sonnendücke, J. Roche, P. Betand, A. Ghizzo, The Semi-Lagangian Method fo the Numeical Resolution of the Vlasov Equation, J. Comput. Phys. 149, (1999). 9
11 l + ω i i i i i Table 1: Bounday conditions of Poposition 1.3 with Q = (Neumann fo mode and Diichlet fo othe modes) l + ω i i i i i Table 2: Bounday conditions of Remak 1 (Neumann fo all the modes). l + ω i i i i i Table 3: Bounday conditions of Poposition 1.2 (Diichlet fo all the modes). 1
12 l + I(ω), fo N = 256, N = 512, N = ,.2993, ,.36899, ,.38498, ,.32539, ,.33783, Table 4: Numeical esolution of the dispesion elation with bounday conditions of Poposition 1.2 (Diichlet fo all the modes). l + I(ω) time inteval of validity [26, 55] [26, 72] [26, 51] [22, 62] [26, 61] Table 5: Gowth ate fo the semi-lagangian method with bounday conditions of Poposition 1.2 (Diichlet fo all the modes). N = N θ = 128, t =.5 l + I(ω) time inteval of validity [93, 138] [33, 87] [34, 69] [83, 165] [35, 53] Table 6: Gowth ate fo the semi-lagangian method with bounday conditions of Poposition 1.3 (Neumann fo mode and Diichlet fo the othe modes). N = N θ = 128, t =.5 l + I(ω) time inteval of validity [89, 151] [32, 89] [34, 69] [35, 53] Table 7: Gowth ate fo the semi-lagangian method with bounday conditions of Poposition 1.3 (Neumann fo all the modes). N = N θ = 128, t =.5 11
13 Figue 1: (Left) Squae modulus of the 7 th Fouie mode of max min Φ(t,, θ)d vs time t fo bounday condition of Poposition 1.3 with Q = (Neumann fo mode and Diichlet fo othe modes). (Right) Density ρ at t = 95. Discetization paametes ae N = 512, N θ = 256 and t =.5. Figue 2: Time evolution of electic enegy (left) and elative mass eo (ight) fo Neumann and Neumann mode bounday conditions, with diffeent discetizations (N N θ t on legend). 12
14 Figue 3: Long time evolution of electic enegy (left) and elative mass eo (ight) fo Neumann bounday conditions, with diffeent discetizations (N N θ t on legend). Figue 4: Long time evolution of electic enegy (left) and elative mass eo (ight) fo Neumann mode bounday conditions, with diffeent discetizations (N N θ t on legend). Figue 5: Long time evolution of electic enegy (left) and elative mass eo (ight) fo Diichlet bounday conditions, with diffeent discetizations (N N θ t on legend). 13
1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationA Probabilistic Approach to Susceptibility Measurement in a Reverberation Chamber
A Pobabilistic Appoach to Susceptibility Measuement in a Revebeation Chambe Emmanuel Amado, Chistophe Lemoine, Philippe Besnie To cite this vesion: Emmanuel Amado, Chistophe Lemoine, Philippe Besnie. A
More informationSupplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract
Supplementay mateial fo the pape Platonic Scatteing Cancellation fo Bending Waves on a Thin Plate M. Fahat, 1 P.-Y. Chen, 2 H. Bağcı, 1 S. Enoch, 3 S. Guenneau, 3 and A. Alù 2 1 Division of Compute, Electical,
More information8 Separation of Variables in Other Coordinate Systems
8 Sepaation of Vaiables in Othe Coodinate Systems Fo the method of sepaation of vaiables to succeed you need to be able to expess the poblem at hand in a coodinate system in which the physical boundaies
More informationCOLLISIONLESS PLASMA PHYSICS TAKE-HOME EXAM
Honou School of Mathematical and Theoetical Physics Pat C Maste of Science in Mathematical and Theoetical Physics COLLISIONLESS PLASMA PHYSICS TAKE-HOME EXAM HILARY TERM 18 TUESDAY, 13TH MARCH 18, 1noon
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationA New Approach to General Relativity
Apeion, Vol. 14, No. 3, July 7 7 A New Appoach to Geneal Relativity Ali Rıza Şahin Gaziosmanpaşa, Istanbul Tukey E-mail: aizasahin@gmail.com Hee we pesent a new point of view fo geneal elativity and/o
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More informationClassical Mechanics Homework set 7, due Nov 8th: Solutions
Classical Mechanics Homewok set 7, due Nov 8th: Solutions 1. Do deivation 8.. It has been asked what effect does a total deivative as a function of q i, t have on the Hamiltonian. Thus, lets us begin with
More informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More informationMultipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source
Multipole Radiation Febuay 29, 26 The electomagnetic field of an isolated, oscillating souce Conside a localized, oscillating souce, located in othewise empty space. We know that the solution fo the vecto
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationCompactly Supported Radial Basis Functions
Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically
More informationBut for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,
Stella elaxation Time [Chandasekha 1960, Pinciples of Stella Dynamics, Chap II] [Ostike & Davidson 1968, Ap.J., 151, 679] Do stas eve collide? Ae inteactions between stas (as opposed to the geneal system
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More informationModeling and Calculation of Optical Amplification in One Dimensional Case of Laser Medium Using Finite Difference Time Domain Method
Jounal of Physics: Confeence Seies PAPER OPEN ACCESS Modeling and Calculation of Optical Amplification in One Dimensional Case of Lase Medium Using Finite Diffeence Time Domain Method To cite this aticle:
More informationDo Managers Do Good With Other People s Money? Online Appendix
Do Manages Do Good With Othe People s Money? Online Appendix Ing-Haw Cheng Haison Hong Kelly Shue Abstact This is the Online Appendix fo Cheng, Hong and Shue 2013) containing details of the model. Datmouth
More informationEN40: Dynamics and Vibrations. Midterm Examination Thursday March
EN40: Dynamics and Vibations Midtem Examination Thusday Mach 9 2017 School of Engineeing Bown Univesity NAME: Geneal Instuctions No collaboation of any kind is pemitted on this examination. You may bing
More informationLecture 7: Angular Momentum, Hydrogen Atom
Lectue 7: Angula Momentum, Hydogen Atom Vecto Quantization of Angula Momentum and Nomalization of 3D Rigid Roto wavefunctions Conside l, so L 2 2 2. Thus, we have L 2. Thee ae thee possibilities fo L z
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationTransverse Wakefield in a Dielectric Tube with Frequency Dependent Dielectric Constant
ARDB-378 Bob Siemann & Alex Chao /4/5 Page of 8 Tansvese Wakefield in a Dielectic Tube with Fequency Dependent Dielectic Constant This note is a continuation of ARDB-368 that is now extended to the tansvese
More information-Δ 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.
The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the
More informationON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},
ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability
More informationSupplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in
Supplementay Figue 1. Cicula paallel lamellae gain size as a function of annealing time at 50 C. Eo bas epesent the σ uncetainty in the measued adii based on image pixilation and analysis uncetainty contibutions
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationHomework 7 Solutions
Homewok 7 olutions Phys 4 Octobe 3, 208. Let s talk about a space monkey. As the space monkey is oiginally obiting in a cicula obit and is massive, its tajectoy satisfies m mon 2 G m mon + L 2 2m mon 2
More informationOPTIMISATION OF A DRIVE SYSTEM AND ITS EPICYCLIC GEAR SET
OPTIISATION OF A DRIVE SYSTE AND ITS EPICYCLIC GEAR SET Nicolas Bellegade, Philippe Dessante, Piee Vidal, Jean-Claude Vannie To cite this vesion: Nicolas Bellegade, Philippe Dessante, Piee Vidal, Jean-Claude
More informationCOMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS
Pogess In Electomagnetics Reseach, PIER 73, 93 105, 2007 COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS T.-X. Song, Y.-H. Liu, and J.-M. Xiong School of Mechanical Engineeing
More informationNuclear size corrections to the energy levels of single-electron atoms
Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea
More information12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationPHYS 110B - HW #7 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 0B - HW #7 Sping 2004, Solutions by David Pace Any efeenced euations ae fom Giffiths Poblem statements ae paaphased. Poblem 0.3 fom Giffiths A point chage,, moves in a loop of adius a. At time t 0
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More information= e2. = 2e2. = 3e2. V = Ze2. where Z is the atomic numnber. Thus, we take as the Hamiltonian for a hydrogenic. H = p2 r. (19.4)
Chapte 9 Hydogen Atom I What is H int? That depends on the physical system and the accuacy with which it is descibed. A natual stating point is the fom H int = p + V, (9.) µ which descibes a two-paticle
More informationSchool of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007
School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More information=0, (x, y) Ω (10.1) Depending on the nature of these boundary conditions, forced, natural or mixed type, the elliptic problems are classified as
Chapte 1 Elliptic Equations 1.1 Intoduction The mathematical modeling of steady state o equilibium phenomena geneally esult in to elliptic equations. The best example is the steady diffusion of heat in
More informationTransformation of the Navier-Stokes Equations in Curvilinear Coordinate Systems with Maple
Global Jounal of Pue and Applied Mathematics. ISSN 0973-1768 Volume 12, Numbe 4 2016, pp. 3315 3325 Reseach India Publications http://www.ipublication.com/gjpam.htm Tansfomation of the Navie-Stokes Equations
More informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationQuantum Mechanics II
Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationPhysics 121 Hour Exam #5 Solution
Physics 2 Hou xam # Solution This exam consists of a five poblems on five pages. Point values ae given with each poblem. They add up to 99 points; you will get fee point to make a total of. In any given
More informationPreliminary Exam: Quantum Physics 1/14/2011, 9:00-3:00
Peliminay Exam: Quantum Physics /4/ 9:-: Answe a total of SIX questions of which at least TWO ae fom section A and at least THREE ae fom section B Fo you answes you can use eithe the blue books o individual
More informationS7: Classical mechanics problem set 2
J. Magoian MT 9, boowing fom J. J. Binney s 6 couse S7: Classical mechanics poblem set. Show that if the Hamiltonian is indepdent of a genealized co-odinate q, then the conjugate momentum p is a constant
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More information763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012
763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012 1. Continuous Random Walk Conside a continuous one-dimensional andom walk. Let w(s i ds i be the pobability that the length of the i th displacement
More informationExceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs
(May 14, 2011 Exceptional egula singula points of second-ode ODEs Paul Gaett gaett@math.umn.edu http://www.math.umn.edu/ gaett/ 1. Solving second-ode ODEs 2. Examples 3. Convegence Fobenius method fo solving
More informationBlack Body Radiation and Radiometric Parameters:
Black Body Radiation and Radiometic Paametes: All mateials absob and emit adiation to some extent. A blackbody is an idealization of how mateials emit and absob adiation. It can be used as a efeence fo
More informationChapter 7-8 Rotational Motion
Chapte 7-8 Rotational Motion What is a Rigid Body? Rotational Kinematics Angula Velocity ω and Acceleation α Unifom Rotational Motion: Kinematics Unifom Cicula Motion: Kinematics and Dynamics The Toque,
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationNumerical approximation to ζ(2n+1)
Illinois Wesleyan Univesity Fom the SelectedWoks of Tian-Xiao He 6 Numeical appoximation to ζ(n+1) Tian-Xiao He, Illinois Wesleyan Univesity Michael J. Dancs Available at: https://woks.bepess.com/tian_xiao_he/6/
More informationKEPLER S LAWS OF PLANETARY MOTION
EPER S AWS OF PANETARY MOTION 1. Intoduction We ae now in a position to apply what we have leaned about the coss poduct and vecto valued functions to deive eple s aws of planetay motion. These laws wee
More informationLINEAR PLATE BENDING
LINEAR PLATE BENDING 1 Linea plate bending A plate is a body of which the mateial is located in a small egion aound a suface in the thee-dimensional space. A special suface is the mid-plane. Measued fom
More informationHomework # 3 Solution Key
PHYSICS 631: Geneal Relativity Homewok # 3 Solution Key 1. You e on you hono not to do this one by hand. I ealize you can use a compute o simply look it up. Please don t. In a flat space, the metic in
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More informationGoodness-of-fit for composite hypotheses.
Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.
More informationPhysics 506 Winter 2006 Homework Assignment #9 Solutions
Physics 506 Winte 2006 Homewok Assignment #9 Solutions Textbook poblems: Ch. 12: 12.2, 12.9, 12.13, 12.14 12.2 a) Show fom Hamilton s pinciple that Lagangians that diffe only by a total time deivative
More informationPhys-272 Lecture 17. Motional Electromotive Force (emf) Induced Electric Fields Displacement Currents Maxwell s Equations
Phys-7 Lectue 17 Motional Electomotive Foce (emf) Induced Electic Fields Displacement Cuents Maxwell s Equations Fom Faaday's Law to Displacement Cuent AC geneato Magnetic Levitation Tain Review of Souces
More informationHopefully Helpful Hints for Gauss s Law
Hopefully Helpful Hints fo Gauss s Law As befoe, thee ae things you need to know about Gauss s Law. In no paticula ode, they ae: a.) In the context of Gauss s Law, at a diffeential level, the electic flux
More informationMath Notes on Kepler s first law 1. r(t) kp(t)
Math 7 - Notes on Keple s fist law Planetay motion and Keple s Laws We conside the motion of a single planet about the sun; fo simplicity, we assign coodinates in R 3 so that the position of the sun is
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More information2 Governing Equations
2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More informationCrack Propagation and the Wave Equation on time dependent domains
Cack Popagation and te Wave Equation on time dependent domains A Majo Qualifying Poject Repot submitted to te Faculty of WPI in patial fulfillment of te equiements fo te Degee of Bacelo of Science in Matematics.
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationOn the Stabilizing Effect of Convection in 3D Incompressible Flows
On the Stabiliing Effect of Convection in 3D Incompessible Flows Thomas Y. Hou Zhen Lei Febuay 5, 8 Abstact We investigate the stabiliing effect of convection in 3D incompessible Eule and Navie- Stokes
More informationCHAPTER 25 ELECTRIC POTENTIAL
CHPTE 5 ELECTIC POTENTIL Potential Diffeence and Electic Potential Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic foce on the paticle given by F=E. When
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More informationDonnishJournals
DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş
More informationarxiv: v1 [physics.pop-ph] 3 Jun 2013
A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationarxiv:gr-qc/ v2 8 Jun 2006
On Quantization of the Electical Chage Mass Dmitiy M Palatnik 1 6400 N Sheidan Rd 2605, Chicago, IL 60626 axiv:g-qc/060502v2 8 Jun 2006 Abstact Suggested a non-linea, non-gauge invaiant model of Maxwell
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationEnergy Levels Of Hydrogen Atom Using Ladder Operators. Ava Khamseh Supervisor: Dr. Brian Pendleton The University of Edinburgh August 2011
Enegy Levels Of Hydogen Atom Using Ladde Opeatos Ava Khamseh Supeviso: D. Bian Pendleton The Univesity of Edinbugh August 11 1 Abstact The aim of this pape is to fist use the Schödinge wavefunction methods
More informationSolutions. V in = ρ 0. r 2 + a r 2 + b, where a and b are constants. The potential at the center of the atom has to be finite, so a = 0. r 2 + b.
Solutions. Plum Pudding Model (a) Find the coesponding electostatic potential inside and outside the atom. Fo R The solution can be found by integating twice, 2 V in = ρ 0 ε 0. V in = ρ 0 6ε 0 2 + a 2
More informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationC/CS/Phys C191 Shor s order (period) finding algorithm and factoring 11/12/14 Fall 2014 Lecture 22
C/CS/Phys C9 Sho s ode (peiod) finding algoithm and factoing /2/4 Fall 204 Lectue 22 With a fast algoithm fo the uantum Fouie Tansfom in hand, it is clea that many useful applications should be possible.
More informationAnalytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra
Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba
More informationFOURIER-BESSEL SERIES AND BOUNDARY VALUE PROBLEMS IN CYLINDRICAL COORDINATES
FOURIER-BESSE SERIES AND BOUNDARY VAUE PROBEMS IN CYINDRICA COORDINATES The paametic Bessel s equation appeas in connection with the aplace opeato in pola coodinates. The method of sepaation of vaiables
More information