DOCSLIB.ORG

The Schrödinger Wave Equation

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
The Schrödinger Wave Equation Chapter 12 The SchrÄodinger Wave Equation Topics The double-slit experiment. Representing particles by waves. Heisenberg's Uncer- tainty Principle. SchrÄodinger'swave equation. Stationary states. Interpretation of the wave-function. One dimensional solutions for a particle in an in¯nite square potential well. General features of solutions of SchrÄodinger'swave equation. The correspondence principle. This chapter is the centrepiece of our development of quantum mechanics and the `derivation' of the SchrÄodingerwave equation. Let us build up the equation, reviewing the clues provided in the pre- vious chapters. We have established that the wave- particle duality works in both directions: ² waves have particle properties, most vividly demonstrated by the photoelectric e®ect and the derivation of the Planck distribution. For photons, hº ¹h! E = hº = ¹h! p = = = ¹hk: c c ² particles have wave properties as illustrated by electron di®raction experiments. h p = ¹hk = : ¸ 1 The SchrÄodingerWave Equation 2 12.1 The Double-slit Experiment The famous double-slit experiment was ¯rst car- ried out by Thomas Young in Emmanuel College in 1802 { the importance of the experiment was that it demonstrated beyond any reasonable doubt the correctness of the wave theory of light. A pair of narrow slits separated by distance d is illuminated by a uniform plane wave and the di®racted rays observed on a screen at some distance R from the slits. A series of light and dark bands is observed, corresponding to constructive and destructive in- terference between the light rays from the two slits. The geometry of the experiment is shown in Figure 12.1. Let us revise the calculation of the observed inten- Figure 12.1. The geometry of Young's sity, but now using complex wave notation double-slit experiment. à = <[A ei(kx¡!t)] = < fA exp [i(kx ¡ !t)]g : where < means `take the real part of'. It will turn out that this is not just a mathematical conve- nience, but rather goes right to the heart of quan- tum mechanics. Complex numbers are the natural language of quantum mechanics. For convenience, we will drop the explicit use of <. From Figure 12.1, provided R À d, the paths travelled from the two slits to the observation point P are: d d r ¼ r ¡ sin θ and r ¼ r + sin θ: 1 2 2 2 We use Huyghens' construction according to which we can replace the slits by sources of light which are in phase and have the same amplitude. Therefore, the amplitude of the wave at P is à = A exp [i(kr1 ¡ !t)] + A exp [i(kr2 ¡ !t)] d d ¼ A exp i[k(r ¡ 2 sin θ) ¡ !t] + A exp i[k(r + 2 sin θ) ¡ !t] = A exp i(kr ¡ !t) fexp [ik(d=2) sin θ] + exp [¡ik(d=2) sin θ]g = 2A exp i(kr ¡ !t) cos[k(d=2) sin θ]: The intensity is proportional to the square of the The SchrÄodingerWave Equation 3 amplitude and so I / 4A2 cos2[k(d=2) sin θ] (12.1) (Figure 12.2). For small values of θ, sin θ ¼ θ and so the intensity distribution becomes µ ¶ kθd I / 4A2 cos2 : (12.2) 2 There is an important piece of mathematics in what Figure 12.2 Di®raction pattern from double we have just done. In the complex representation, slits. the intensity can be found by taking the square of the modulus of the complex function Ã, that is, I = Ãä = jÃj2; where ä is the complex conjugate of Ã. We do not need to take the real part at all. Let us repeat the double-slit experiment using a photon counting detector. In this case, the arrival The experiment with light was ¯rst performed by of each photon on the screen is detected separately. G.I. Taylor in 1909, and the electron di®raction When the light is of high intensity, we observe the experiment was performed by JÄonnsonin 1961 us- 2 ing 50 keV electrons and a slit separation of only usual cos (kθd=2) variation of light intensity on the 1 ¹m. screen (Figure 12.3(c)). As the light intensity is de- creased, we begin to see the pattern `shimmering' as we begin to see the individual photons arriving (Figure 12.3(b)). At the very lowest light intensi- ties, we observe individual photons arriving at the screen (Figure 12.3(a)). The extraordinary thing about this experiment is that the individual par- ticles of light know where to appear on the screen to produce the standard di®raction pattern, despite the fact that they only arrive one at a time. To put it another way, the physics of light quanta must be such that there is zero probability that the individ- ual photons will land on the screen at the nulls in the di®raction pattern. It is interesting to quote the words of Richard Feynman about the signif- icance of the double-slit experiment for quantum mechanics. `We choose to examine a phenomenon which is impossible, absolutely impossi- ble, to explain in any classical way, and which has in it the heart of quantum mechanics. In reality, it contains the Figure 12.3. Young's double-slit experiment only mystery.' carried out using a photon counting detector. The SchrÄodingerWave Equation 4 These results provide the clue as to how we should interpret the cos2(kθd=2) intensity distribution. This function describes the probability that a single par- ticle arrives at any particular location on the detec- tor. Therefore, the quantum mechanical wave func- tion Ã(x) will be used to describe a particle such as The Wave Function an electron with the property that the probability The wave function is de¯ned so that the of ¯nding the particle in the interval x to x + dx is probability of ¯nding the particle in the interval x to x + dx is p(x) dx = jÃj2 dx = Ãä dx; (12.3) p(x) dx = jÃj2 dx = Ãä dx which has the form of a probability density func- tion. This represents a huge change from classical me- chanics. According to classical physics we can spec- ify precisely the position and momentum of a par- ticle. This has been replaced by a probabilistic de- scription of phenomena at the atomic level. This di®erent perspective is forced upon us as a result of the need to describe particles by waves and vice versa. The consequences are very profound as we shall see. Our next task is to determine the di®erential equa- tion which determines the wave function à { that will turn out to be the SchrÄodingerwave equation. First, however, let us investigate how we go about representing particles by waves. 12.2 The Representation of Particles by Waves We need a means of localising the particle in space { in other words, of ¯nding a probability density which describes where the particle is located. This is done by developing the idea of a wave packet which we construct out of a superposition of waves. This leads to a consideration of the properties of Fourier series and Fourier transforms, which come only late in the mathematics course. We need only one key idea from these topics and that is that we can use sums of sine and cosine functions to repre- sent any function f(x) within some range x1 · x · x2. The reason this works is that, mathematically, The SchrÄodingerWave Equation 5 the sequences of sine waves sin k0x; sin 2k0x; sin 3k0x; : : : plus the series of cosine waves cos k0x; cos 2k0x; cos 3k0x; : : : form a complete, orthogonal set of functions, mean- ing that we can synthesise any pattern we like by selecting appropriately the constants A0;A1, A2, ::: and B1, B2, ::: in the series expansion X1 X1 Ã(x) = A0 + An cos nk0x + Bn sin nk0x: n=1 n=1 (12.4) This sum is the Fourier series for the function Ã(x) de¯ned in the interval x1 · x · x2. These ideas suggest a way in which we can rep- resent particles by waves. As an example, con- sider building up a gaussian distribution in the x- direction, which is to represent the probability den- sity distribution of the particle in that coordinate. Suppose the particle has momentum p, so that, according to the de Broglie hypothesis, the wave- length associated with it is ¸0 = h=p. The corre- sponding wavevector is 2¼ 2¼p k0 = = : (12.5) ¸0 h Figure 12.4. The distribution of terms The waves associated with k are sine or cosine 0 in the Fourier Series with ¢k = 0:05. waves, A sin k x or A cos k x, and they are not lo- 0 0 0 We are bearing in mind that in the calised at all { the extend all the way from ¡1 to limit of a large number of terms, the +1 and this is not very helpful. binomial series tends to a gaussian Now, let us start with A cos k0x and add some other distribution. cosine waves and, for illustrative purposes, let us add together a binomial sequence of terms, cen- . tred on k0, but separated from it by §n ¢k0, where . ....... .. p(k) . .. .. .. ¢k ¿ k . Figure 12.4 shows the binomial distri- . .. 0 0 . .. .. ................ ...... bution of amplitudes of cosine waves, which are the . ...... ...... ...... ...... 8 . ...... ...... coe±cients of the expansion of (1 + x) , that is, . ...... ...... ...... ...... ...... ...... ...... 8 2 3 4 . ...... .................... (1 + x) = 1 + 8x + 28x + 56x + 70x . .......... .......... .......... .......... 5 6 7 8 . .......... .......... + 56x + 28x + 8x + x ; . .......... .......... .......... .......... ............................ .............. .............. .............. where we have chosen k0 = 1 and ¢k0 = 0:05. .................................................................................................................................................................................................................................................................................................................................................................................................................... .......................... .... .... .... .......... 0 k0 k The SchrÄodingerWave Equation 6 The result of adding together this series of cosine waves is shown in Figure 12.5 by the oscillating line which has maximum at the origin, but tends rapidly to zero in a ¯nite distance rather than os- cillating to in¯nity.
Load more

Recommended publications
  • Glossary Physics (I-Introduction)
    1 Glossary Physics (I-introduction) - Efficiency: The percent of the work put into a machine that is converted into useful work output; = work done / energy used [-]. = eta In machines: The work output of any machine cannot exceed the work input (<=100%); in an ideal machine, where no energy is transformed into heat: work(input) = work(output), =100%. Energy: The property of a system that enables it to do work. Conservation o. E.: Energy cannot be created or destroyed; it may be transformed from one form into another, but the total amount of energy never changes. Equilibrium: The state of an object when not acted upon by a net force or net torque; an object in equilibrium may be at rest or moving at uniform velocity - not accelerating. Mechanical E.: The state of an object or system of objects for which any impressed forces cancels to zero and no acceleration occurs. Dynamic E.: Object is moving without experiencing acceleration. Static E.: Object is at rest.F Force: The influence that can cause an object to be accelerated or retarded; is always in the direction of the net force, hence a vector quantity; the four elementary forces are: Electromagnetic F.: Is an attraction or repulsion G, gravit. const.6.672E-11[Nm2/kg2] between electric charges: d, distance [m] 2 2 2 2 F = 1/(40) (q1q2/d ) [(CC/m )(Nm /C )] = [N] m,M, mass [kg] Gravitational F.: Is a mutual attraction between all masses: q, charge [As] [C] 2 2 2 2 F = GmM/d [Nm /kg kg 1/m ] = [N] 0, dielectric constant Strong F.: (nuclear force) Acts within the nuclei of atoms: 8.854E-12 [C2/Nm2] [F/m] 2 2 2 2 2 F = 1/(40) (e /d ) [(CC/m )(Nm /C )] = [N] , 3.14 [-] Weak F.: Manifests itself in special reactions among elementary e, 1.60210 E-19 [As] [C] particles, such as the reaction that occur in radioactive decay.
    [Show full text]
  • Linear Elastodynamics and Waves
    Linear Elastodynamics and Waves M. Destradea, G. Saccomandib aSchool of Electrical, Electronic, and Mechanical Engineering, University College Dublin, Belfield, Dublin 4, Ireland; bDipartimento di Ingegneria Industriale, Universit`adegli Studi di Perugia, 06125 Perugia, Italy. Contents 1 Introduction 3 2 Bulk waves 6 2.1 Homogeneous waves in isotropic solids . 8 2.2 Homogeneous waves in anisotropic solids . 10 2.3 Slowness surface and wavefronts . 12 2.4 Inhomogeneous waves in isotropic solids . 13 3 Surface waves 15 3.1 Shear horizontal homogeneous surface waves . 15 3.2 Rayleigh waves . 17 3.3 Love waves . 22 3.4 Other surface waves . 25 3.5 Surface waves in anisotropic media . 25 4 Interface waves 27 4.1 Stoneley waves . 27 4.2 Slip waves . 29 4.3 Scholte waves . 32 4.4 Interface waves in anisotropic solids . 32 5 Concluding remarks 33 1 Abstract We provide a simple introduction to wave propagation in the frame- work of linear elastodynamics. We discuss bulk waves in isotropic and anisotropic linear elastic materials and we survey several families of surface and interface waves. We conclude by suggesting a list of books for a more detailed study of the topic. Keywords: linear elastodynamics, anisotropy, plane homogeneous waves, bulk waves, interface waves, surface waves. 1 Introduction In elastostatics we study the equilibria of elastic solids; when its equilib- rium is disturbed, a solid is set into motion, which constitutes the subject of elastodynamics. Early efforts in the study of elastodynamics were mainly aimed at modeling seismic wave propagation. With the advent of electron- ics, many applications have been found in the industrial world.
    [Show full text]
  • 1 Fundamental Solutions to the Wave Equation 2 the Pulsating Sphere
    1 Fundamental Solutions to the Wave Equation Physical insight in the sound generation mechanism can be gained by considering simple analytical solutions to the wave equation. One example is to consider acoustic radiation with spherical symmetry about a point ~y = fyig, which without loss of generality can be taken as the origin of coordinates. If t stands for time and ~x = fxig represent the observation point, such solutions of the wave equation, @2 ( − c2r2)φ = 0; (1) @t2 o will depend only on the r = j~x − ~yj. It is readily shown that in this case (1) can be cast in the form of a one-dimensional wave equation @2 @2 ( − c2 )(rφ) = 0: (2) @t2 o @r2 The general solution to (2) can be written as f(t − r ) g(t + r ) φ = co + co : (3) r r The functions f and g are arbitrary functions of the single variables τ = t± r , respectively. ± co They determine the pattern or the phase variation of the wave, while the factor 1=r affects only the wave magnitude and represents the spreading of the wave energy over larger surface as it propagates away from the source. The function f(t − r ) represents an outwardly co going wave propagating with the speed c . The function g(t + r ) represents an inwardly o co propagating wave propagating with the speed co. 2 The Pulsating Sphere Consider a sphere centered at the origin and having a small pulsating motion so that the equation of its surface is r = a(t) = a0 + a1(t); (4) where ja1(t)j << a0.
    [Show full text]
  • Phys 234 Vector Calculus and Maxwell's Equations Prof. Alex
    Prof. Alex Small Phys 234 Vector Calculus and Maxwell's [email protected] This document starts with a summary of useful facts from vector calculus, and then uses them to derive Maxwell's equations. First, definitions of vector operators. 1. Gradient Operator: The gradient operator is something that acts on a function f and produces a vector whose components are equal to derivatives of the function. It is defined by: @f @f @f rf = ^x + ^y + ^z (1) @x @y @z 2. Divergence: We can apply the gradient operator to a vector field to get a scalar function, by taking the dot product of the gradient operator and the vector function. @V @V @V r ⋅ V~ (x; y; z) = x + y + z (2) @x @y @z 3. Curl: We can take the cross product of the gradient operator with a vector field, and get another vector that involves derivatives of the vector field. @V @V @V @V @V @V r × V~ (x; y; z) = ^x z − y + ^y x − z + ^z y − x (3) @y @z @z @x @x @y 4. We can also define a second derivative of a vector function, called the Laplacian. The Laplacian can be applied to any function of x, y, and z, whether the function is a vector or a scalar. The Laplacian is defined as: @2f @2f @2f r2f(x; y; z) = r ⋅ rf = + + (4) @x2 @y2 @z2 Notes that the Laplacian of a scalar function is equal to the divergence of the gradient. Second, some useful facts about vector fields: 1. If a vector field V~ (x; y; z) has zero curl, then it can be written as the gradient of a scalar function f(x; y; z).
    [Show full text]
  • Fundamentals of Acoustics Introductory Course on Multiphysics Modelling
    Introduction Acoustic wave equation Sound levels Absorption of sound waves Fundamentals of Acoustics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELINSKI´ bluebox.ippt.pan.pl/˜tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw • Poland 3 Sound levels Sound intensity and power Decibel scales Sound pressure level 2 Acoustic wave equation Equal-loudness contours Assumptions Equation of state Continuity equation Equilibrium equation 4 Absorption of sound waves Linear wave equation Mechanisms of the The speed of sound acoustic energy dissipation Inhomogeneous wave A phenomenological equation approach to absorption Acoustic impedance The classical absorption Boundary conditions coefficient Introduction Acoustic wave equation Sound levels Absorption of sound waves Outline 1 Introduction Sound waves Acoustic variables 3 Sound levels Sound intensity and power Decibel scales Sound pressure level Equal-loudness contours 4 Absorption of sound waves Mechanisms of the acoustic energy dissipation A phenomenological approach to absorption The classical absorption coefficient Introduction Acoustic wave equation Sound levels Absorption of sound waves Outline 1 Introduction Sound waves Acoustic variables 2 Acoustic wave equation Assumptions Equation of state Continuity equation Equilibrium equation Linear wave equation The speed of sound Inhomogeneous wave equation Acoustic impedance Boundary conditions 4 Absorption of sound waves Mechanisms of the acoustic energy dissipation A phenomenological approach
    [Show full text]
  • Elastic Wave Equation
    Introduction to elastic wave equation Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012 Outline • Motivation • Elastic wave equation – Equation of motion, Definitions and The linear Stress- Strain relationship • The Seismic Wave Equation in Isotropic Media • Seismic wave equation in homogeneous media • Acoustic wave equation • Short Summary Motivation • Elastic wave equation has been widely used to describe wave propagation in an elastic medium, such as seismic waves in Earth and ultrasonic waves in human body. • Seismic waves are waves of energy that travel through the earth, and are a result of an earthquake, explosion, or a volcano. Elastic wave equation • The standard form for seismic elastic wave equation in homogeneous media is : ρu ( 2).u u ρ : is the density u : is the displacement , : Lame parameters Equation of Motion • We will depend on Newton’s second law F=ma m : mass ρdx1dx2dx3 2u a : acceleration t 2 • The total force from stress field: body F Fi Fi body Fi fidx1dx2dx3 ij F dx dx dx dx dx dx i x j 1 2 3 j ij 1 2 3 Equation of Motion • Combining these information together we get the Momentum equation (Equation of Motion) 2u 2 j f t ij i • where is the density, u is the displacement, and is the stress tensor. Definitions • Stress : A measure of the internal forces acting within a deformable body. (The force acting on a solid to deform it) The stress at any point in an object, assumed to behave as a continuum, is completely defined by nine component stresses: three orthogonal normal stresses and six orthogonal shear stresses.
    [Show full text]
  • The Wave Equation (30)
    GG 454 March 31, 2002 1 THE WAVE EQUATION (30) I Main Topics A Assumptions and boundary conditions used in 2-D small wave theory B The Laplace equation and fluid potential C Solution of the wave equation D Energy in a wavelength E Shoaling of waves II Assumptions and boundary conditions used in 2-D small wave theory Small amplitude surface wave y L H = 2A = wave height η x Still water u level ε v ζ Depth above Water depth = +d bottom = d+y Particle Particle orbit velocity y = -d Modified from Sorensen, 1978 A No geometry changes parallel to wave crest (2-D assumption) B Wave amplitude is small relative to wave length and water depth; it will follow that particle velocities are small relative to wave speed C Water is homogeneous, incompressible, and surface tension is nil. D The bottom is not moving, is impermeable, and is horizontal E Pressure along air-sea interface is constant F The water surface has the form of a cosine wave Hx t ηπ=−cos2 2 L T Stephen Martel 30-1 University of Hawaii GG 454 March 31, 2002 2 III The Laplace equation and fluid potential A Conservation of mass Consider a box the shape of a cube, with fluid flowing in and out or it. ∆z v+∆v (out) u+∆u (out) ∆ u y (in) ∆x v (in) The mass flow rate in the left side of the box is: ∆m ∆ρV ∆V ∆∆xyz∆ ∆x 1a 1 = = ρρ= =∆ ρyz∆ =∆∆()( ρyzu )() ∆t ∆t ∆t ∆t ∆t The mass flow rate out the right side is ∆m 1b 2 =∆∆+()(ρ yzu )(∆ u ).
    [Show full text]
  • The Acoustic Wave Equation and Simple Solutions
    Chapter 5 THE ACOUSTIC WAVE EQUATION AND SIMPLE SOLUTIONS 5.1INTRODUCTION Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluido In addition to the audible pressure fields of modera te intensity, the most familiar, there are also ultrasonic and infrasonic waves whose frequencies lie beyond the limits of hearing, high-intensity waves (such as those near jet engines and missiles) that may produce a sensation of pain rather than sound, nonlinear waves of still higher intensities, and shock waves generated by explosions and supersonic aircraft. lnviscid fluids exhibit fewer constraints to deformations than do solids. The restoring forces responsible for propagating a wave are the pressure changes that oc­ cur when the fluid is compressed or expanded. Individual elements of the fluid move back and forth in the direction of the forces, producing adjacent regions of com­ pression and rarefaction similar to those produced by longitudinal waves in a bar. The following terminology and symbols will be used: r = equilibrium position of a fluid element r = xx + yy + zz (5.1.1) (x, y, and z are the unit vectors in the x, y, and z directions, respectively) g = particle displacement of a fluid element from its equilibrium position (5.1.2) ü = particle velocity of a fluid element (5.1.3) p = instantaneous density at (x, y, z) po = equilibrium density at (x, y, z) s = condensation at (x, y, z) 113 114 CHAPTER 5 THE ACOUSTIC WAVE EQUATION ANO SIMPLE SOLUTIONS s = (p - pO)/ pO (5.1.4) p - PO = POS = acoustic density at (x, y, Z) i1f = instantaneous pressure at (x, y, Z) i1fO = equilibrium pressure at (x, y, Z) P = acoustic pressure at (x, y, Z) (5.1.5) c = thermodynamic speed Of sound of the fluid <I> = velocity potential of the wave ü = V<I> .
    [Show full text]
  • Chapter 4 – the Wave Equation and Its Solution in Gases and Liquids
    Chapter 4 – The Wave Equation and its Solution in Gases and Liquids Slides to accompany lectures in Vibro-Acoustic Design in Mechanical Systems © 2012 by D. W. Herrin Department of Mechanical Engineering University of Kentucky Lexington, KY 40506-0503 Tel: 859-218-0609 [email protected] Simplifying Assumptions • The medium is homogenous and isotropic • The medium is linearly elastic • Viscous losses are negligible • Heat transfer in the medium can be ignored • Gravitational effects can be ignored • Acoustic disturbances are small Dept. of Mech. Engineering 2 University of Kentucky ME 510 Vibro-Acoustic Design Sound Pressure and Particle Velocity pt (r,t) = p0 + p(r,t) Total Pressure Undisturbed Pressure Disturbed Pressure Particle Velocity ˆ ˆ ˆ u(r,t) = uxi + uy j + uzk Dept. of Mech. Engineering 3 University of Kentucky ME 510 Vibro-Acoustic Design Density and Temperature ρt (r,t) = ρ0 + ρ (r,t) Total Density Undisturbed Density Disturbed Density Absolute Temperature T (r,t) Dept. of Mech. Engineering 4 University of Kentucky ME 510 Vibro-Acoustic Design Equation of Continuity A fluid particle is a volume element large enough to contain millions of molecules but small enough that the acoustic disturbances can be thought of as varying linearly across an element. Δy # ∂ρ &# ∂u & ρ + t Δx u + x Δx ρ u % t (% x ( x x Δz $ ∂x '$ ∂x ' Δx Dept. of Mech. Engineering 5 University of Kentucky ME 510 Vibro-Acoustic Design Equation of Continuity Δy # ∂ρt &# ∂ux & %ρt + Δx(%ux + Δx( ρ u $ x '$ x ' x x Δz ∂ ∂ Δx ∂ $ ∂ρt '$ ∂ux ' (ρtΔxΔyΔz) = ρtuxΔyΔz −&ρt + Δx)&ux + Δx)ΔyΔz ∂t % ∂x (% ∂x ( Dept.
    [Show full text]
  • Analysis and Computation of the Elastic Wave Equation with Random Coefficients
    Computers and Mathematics with Applications 70 (2015) 2454–2473 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Analysis and computation of the elastic wave equation with random coefficients Mohammad Motamed a,∗, Fabio Nobile b, Raúl Tempone c a Department of Mathematics and Statistics, The University of New Mexico, Albuquerque, NM 87131, USA b MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland c King Abdullah University of Science and Technology, Jeddah, Saudi Arabia article info a b s t r a c t Article history: We consider the stochastic initial–boundary value problem for the elastic wave equation Received 23 August 2014 with random coefficients and deterministic data. We propose a stochastic collocation Received in revised form 6 September 2015 method for computing statistical moments of the solution or statistics of some given Accepted 12 September 2015 quantities of interest. We study the convergence rate of the error in the stochastic Available online 9 October 2015 collocation method. In particular, we show that, the rate of convergence depends on the regularity of the solution or the quantity of interest in the stochastic space, which is in Keywords: turn related to the regularity of the deterministic data in the physical space and the type Uncertainty quantification Stochastic partial differential equations of the quantity of interest. We demonstrate that a fast rate of convergence is possible in Elastic wave equation two cases: for the elastic wave solutions with high regular data; and for some high regular Collocation method quantities of interest even in the presence of low regular data.
    [Show full text]
  • Simple Derivation of Electromagnetic Waves from Maxwell's Equations
    Simple Derivation of Electromagnetic Waves from Maxwell’s Equations By Lynda Williams, Santa Rosa Junior College Physics Department Assume that the electric and magnetic fields are constrained to the y and z directions, respectfully, and that they are both functions of only x and t. This will result in a linearly polarized plane wave travelling in the x direction at the speed of light c. Ext( , ) Extj ( , )ˆ and Bxt( , ) Bxtk ( , ) ˆ We will derive the wave equation from Maxwell’s Equations in free space where I and Q are both zero. EB 0 0 BE EB με tt00 iˆˆ j kˆ E Start with Faraday’s Law. Take the curl of the E field: E( x , t )ˆ j = kˆ x y z x 0E ( x , t ) 0 EB Equating magnitudes in Faraday’s Law: (1) xt This means that the spatial variation of the electric field gives rise to a time-varying magnetic field, and visa-versa. In a similar fashion we derive a second equation from Ampere-Maxwell’s Law: iˆˆ j kˆ BBE Take the curl of B: B( x , t ) kˆ = ˆ j (2) x y z x x00 t 0 0B ( x , t ) This means that the spatial variation of the magnetic field gives rise to a time-varying electric field, and visa-versa. Take the partial derivative of equation (1) with respect to x and combining the results from (2): 22EBBEE 22 = 0 0 0 0 x x t t x t t t 22EE (3) xt2200 22BB In a similar manner, we can derive a second equation for the magnetic field: (4) xt2200 Both equations (3) and (4) have the form of the general wave equation for a wave (,)xt traveling in 221 the x direction with speed v: .
    [Show full text]
  • On Fourier Transforms and Delta Functions
    CHAPTER 3 On Fourier Transforms and Delta Functions The Fourier transform of a function (for example, a function of time or space) provides a way to analyse the function in terms of its sinusoidal components of different wavelengths. The function itself is a sum of such components. The Dirac delta function is a highly localized function which is zero almost everywhere. There is a sense in which different sinusoids are orthogonal. The orthogonality can be expressed in terms of Dirac delta functions. In this chapter we review the properties of Fourier transforms, the orthogonality of sinusoids, and the properties of Dirac delta functions, in a way that draws many analogies with ordinary vectors and the orthogonality of vectors that are parallel to different coordinate axes. 3.1 Basic Analogies If A is an ordinary three-dimensional spatial vector, then the component of A in each of the three coordinate directions specified by unit vectors xˆ1, xˆ2, xˆ3 is A · xˆi for i = 1, 2, or 3. It follows that the three cartesian components (A1, A2, A3) of A are given by Ai = A · xˆi , for i = 1, 2, or 3. (3.1) We can write out the vector A as the sum of its components in each coordinate direction as follows: 3 A = Ai xˆi . (3.2) i=1 Of course the three coordinate directions are orthogonal, a property that is summarized by the equation xˆi · xˆ j = δij. (3.3) Fourier series are essentially a device to express the same basic ideas (3.1), (3.2) and (3.3), applied to a particular inner product space.
    [Show full text]