Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Everand Logo

Welcome to Everand!

  • Discover millions of ebooks, audiobooks, and so much more with a free trial

    From $11.99/month after trial. Cancel anytime.

    Boundary Value Problems and Fourier Expansions
    Boundary Value Problems and Fourier Expansions
    Boundary Value Problems and Fourier Expansions
    Ebook674 pages3 hours

    Boundary Value Problems and Fourier Expansions

    By Charles R. MacCluer

    ()

    Read preview

    About this ebook

    Based on modern Sobolev methods, this text for advanced undergraduates and graduate students is highly physical in its orientation. It integrates numerical methods and symbolic manipulation into an elegant viewpoint that is consonant with implementation by digital computer. The first five sections form an informal introduction that develops students' physical and mathematical intuition. The following section introduces Hilbert space in its natural environment, and the next six sections pose and solve the standard problems. The final seven sections feature concise introductions to selected topics, including Sturm-Liouville problems, Fourier integrals, Galerkin's method, and Sobolev methods. 1994 edition. 64 figures. Exercises.
    LanguageEnglish
    PublisherDover Publications
    Release dateJan 18, 2013
    ISBN9780486153179
    Boundary Value Problems and Fourier Expansions

    Read more from Charles R. Mac Cluer

    Related to Boundary Value Problems and Fourier Expansions

    Titles in the series (100)

    View More
    View More

    Related ebooks

    Mathematics For You

    View More
    View More

    Related podcast episodes

    Related articles

    • UNLIMITED

      The Math That Takes Newton Into the Quantum World

      Feb 28, 2019

      In my 50s, too old to become a real expert, I have finally fallen in love with algebraic geometry. As the name suggests, this is the study of geometry using algebra. Around 1637, René Descartes laid the groundwork for this subject by taking a plane,

      9 min read

    • UNLIMITED

      Physicists Attack Math’s $1,000,000 Question

      Apr 4, 2017

      Physicists are attempting to map the distribution of the prime numbers to the energy levels of a particular quantum system.

      4 min read

    • UNLIMITED

      Fanning The Flame

      Sep 13, 2023

      Fanning The Flame
      3 min read

    • UNLIMITED

      Is The Law Of Conservation Of Energy Cancelled?

      Dec 12, 2019

      Physics is often baffling, but one principle seems rock-solid: the law of conservation of energy. The world contains this thing called “energy” whose amount never changes. It can change its form or go from one body to another, but its total amount re

      6 min read

    • UNLIMITED

      Are Neural Networks About to Reinvent Physics?

      Nov 21, 2019

      Can AI teach itself the laws of physics? Will classical computers soon be replaced by deep neural networks? Sure looks like it, if you’ve been following the news, which lately has been filled with headlines like, “A neural net solves the three-body p

      9 min read

    • UNLIMITED

      Is The Standard Model Broken?

      Jun 2, 2021

      Is The Standard Model Broken?
      3 min read

    • UNLIMITED

      The Universe Began With a Big Melt, Not a Big Bang

      Oct 5, 2017

      There are two tantalizing mysteries about our universe, one dealing with its final fate and the other with its beginning, that have intrigued cosmologists for decades. The community has always believed these to be independent problems—but what if the

      10 min read

    • UNLIMITED

      Babbage And Running His Many Engines

      Nov 16, 2021

      Mike Bedford Like most of us, Mike tended to think of the likes of Alan Turing and John von Neumann as the pioneers of computing, so delving further into Charles Babbage’s creations was truly inspiring. Discussion about the first ever computer often

      11 min read

    • UNLIMITED

      Do The Maths

      Sep 2, 2024

      by Robyn Arianrhod (UNSW Press, $49.99) Imagine an arrow on a computer screen –in PowerPoint, say. You can drag it around the screen and it’s the same arrow. You can make different arrows by stretching it or rotating it to point in a different direc

      2 min read

    • UNLIMITED

      Thought Experiment Paradox Divides Quantum Experts

      Sep 27, 2018

      A new thought experiment may have exposed a paradox within the theory of quantum mechanics, according to new research. There is likely no other scientific theory that is as well supported as quantum mechanics. For nearly 100 years now, it has repeate

      3 min read

    • UNLIMITED

      How Physics Is Like Three-Chord Rock: Like a set of common chords, the same math appears in diverse fields of physics.

      Jun 12, 2014

      A few years back, 12 million of us clicked over to watch the “Pachelbel Rant” on YouTube. You might remember it. Strumming repetitive chords on his guitar, comedian Rob Paravonian confessed that when he was a cellist, he couldn’t stand the Pachelbel

      6 min read

    • UNLIMITED

      How Spooky Science Helps Us Peer Inside The Planets

      Dec 3, 2020

      An assistant professor of computational science at the EPFL research centre in Lausanne, Switzerland, involved in the current research on metallic hydrogen. Could you explain how the machine-learning techniques used in your research work? Why were th

      1 min read

    • UNLIMITED

      A Supermassive Lens On The Constants Of Nature

      Nov 25, 2020

      The 2020 Nobel Prize in Physics went to three researchers who confirmed that Einstein’s general relativity predicts black holes, and established that the center of our own galaxy houses a supermassive black hole with the equivalent of 4 million suns

      9 min read

    • UNLIMITED

      What Is The Theory Of Everything?

      Jun 20, 2019

      The theory will have to bring together a lot of different aspects of physics, from Einstein’s general theory of relativity that governs our understanding of the largest structures in the universe to the Standard Model that helps us understand the sma

      2 min read

    • UNLIMITED

      What Is Cfd

      Jun 21, 2022

      What Is Cfd
      3 min read

    Related categories

    Reviews for Boundary Value Problems and Fourier Expansions

    0 ratings

    0 ratings0 reviews

    What did you think?

    Tap to rate

    Review must be at least 10 words

      Book preview

      Boundary Value Problems and Fourier Expansions - Charles R. MacCluer

      Universities.

      Chapter 1

      Preliminaries

      We begin with a review of partial derivatives and their two chain rules. Several examples of partial differential equations (PDEs) are solved. The relaxing temperatures within a block are obtained by a numerical model which reveals the PDE governing heat conduction. Finally, the important V operator is introduced and the ‘big six’ PDEs are displayed.

      §1.1 Partial Derivatives

      The most interesting problems involve multiple degrees of freedom: vibrating structures, changing temperatures within solids, voltage potentials within regions, orbiting electrons, growing economies or populations, etc. Astonishingly accurate models of these complicated phenomena have been given using partial differential equations — as relations among the various partial derivatives of measured quantities.

      Recall that the partial derivative at (x⁰, y⁰) of a function f of the two real variables x, y with respect to (say) the first variable x is the limit

      (1.1)

      when it exists. The symbol fx is a common alternate notation for the partial derivative ∂f/∂x. So for instance, if f(x, y) = x³ + y² + x⁷y⁵ + 1, then fx = 3x² + 7x⁶y⁵, while fy = 2y + 5x⁷y⁴. The partial fx represents the sensitivity of f to changes in x while holding y fixed.

      More generally, the directional derivative of f at (x⁰, y⁰) in the direction of the unit vector v = (a, b) is the limit

      (1.2)

      when it exists. The existence of the directional derivative in all directions v from (x⁰, y⁰) is guaranteed when f is differentiable at (x⁰, y⁰, i.e., when f is well approximated by a plane nearby (x⁰, y⁰) — see [Apostol]. It is easy to derive the following formula for the directional derivative:

      (1.3)

      (Exercise 1.1). This formula for the directional derivative generalizes to n variables x = (x1, x2, . . . , xn):

      (1.4)

      where ∇f is the gradient of f, i.e, the vector field

      and where · is dot product. See §1.5. Because

      (1.5)

      the gradient points in the direction of the maximal increase of f.

      The formula for the directional derivative is, in turn, a special case of the first of two chain rules for partial derivatives.

      First Chain Rule. Suppose x = x(t) is a curve in n-space that is differentiable at t = t⁰, and that f(x) is a function of n variables differentiable at x⁰ = x(t⁰). Then

      (1.6)

      In more familiar notation,

      Consequently,

      the gradient is normal to all contour surfaces f = c.

      The first chain rule is itself a special case.

      Second Chain Rule. Suppose f(u) is a differentiable function of the n variables u = (u1, u2, . . . , un), while each ui is itself a differentiable function of the r variables x = (x1, x2, . . . , xr). Then

      (1.7)

      Review these chain rules by working Exercises 1.2-1.6 and 1.19-1.22.

      §1.2 Several Example PDEs

      Let us work through several simple examples of partial differential equations, each of some physical importance.

      Example 1. Consider the simple second order PDE

      (1.8)

      Since (ux)y = 0, ux cannot depend on y, and so ux = h(x). But then integrating with respect to x yields that u must be of the form

      (1.9)

      Conversely note that any function of the form (1.9) with f and g differentiable is indeed a solution of (1.8).

      If we happen to know the values of u(x, y0) and u(x0, y) on the two lines x = x0 and y = y0, then we can recover the functions f and g uniquely within constants that add to 0 (Exercise 1.7). We will see this simple PDE again when we solve the all-important wave equation in Chapter 4.

      Example 2. The first order PDE

      (1.10)

      where a and b are constants, is called the transport equation for reasons explained below. Observe that this equation is a geometric statement: The directional derivative

      ∇u · (a, b) = 0

      of u is zero in the direction (a, b), i.e., u = u(x, y) is constant along curves with tangent vector (a, b). This means u is constant along each line ay − bx = c. But then u = u(x, y) is determined solely by the value c, i.e.,

      (1.11)

      Conversely, note that any u of the form (1.11) clearly satisfies the PDE (1.10) as long as f is differentiable. These lines ay − bx = c are called the characteristic lines of the PDE (1.10).

      Physical realization of the transport equation. Think of fluid flowing through a pipe with velocity v = v(x, t) at location x at time t. This fluid is carrying immiscible particles of density p = p(x, t) per unit length at location x at time t. The total particle mass m within the tube between x = a and x = b is therefore

      (1.12)

      By Leibniz’s principle (§3.2) we may differentiate past the integral to find the instantaneous rate of change of the mass within this length of pipe:

      (1.13)

      On the other hand, this gain of mass is the net flow in from the left less the flow out from the right, i.e.,

      (1.14)

      Equating and dividing by b − a gives that

      (1.15)

      In the limit, as b → a, we have conservation of mass:

      (1.16)

      or in the case velocity v is constant, we have transport:

      (1.17)

      Example 3. Consider a generalization of transport:

      (1.18)

      This can be read as the directional derivative in the direction (a, b) is 0, i.e., u is constant along a curve with tangent vector (a, b). Such curves have normal (−b, a), so locally satisfy

      (1.19)

      Solutions φ(x, y) = c of the differential form (1.19) are called the characteristic curves of the PDE (1.18), one per integration constant c. Then

      (1.20)

      will at least locally solve (1.18) since

      aux + buy = af′(φ)φx + bf′(φ)φy = f′(φ)(−ba + ba) = 0.

      The solution u = u(x, t) is constant along the characteristic curves. Practice this method of characteristics by working Exercises 1.8 - 1. 15.

      §1.3 Transient Heat Flow in a Block

      Let us investigate the flow of heat within a rectangular block of homogeneous isotropic material as shown in Figure 1.1. The top and bottom faces are insulated, thus preventing flow through these faces. Suppose initially the block is uniformly at temperature 1 when suddenly the remaining four side faces are placed and held in contact with heat sinks at temperature 0. Again refer to Figure 1.1.

      Intuitively it is clear that the temperature u within the block will eventually relax to 0 as heat is forced through the four faces x = 0, x = L, y = 0, and y = L. What is not clear is the temperature regime during this relaxation. For instance, can overshoot occur — will portions of the block ever be at a negative temperature? Does this relaxation require an infinite amount of time? Must temperatures be symmetric in x and y? Is there but one trajectory that the regime of temperatures must follow during this relaxation to 0?

      Figure 1.1 A rectangular solid homogeneous block, initially at temperature 1, with top and bottom faces z = 0, 1 insulated. Suddenly at t = 0, the remaining four faces are put in contact with sinks at temperature 0.

      A numerical solution

      A practical method for resolving these questions is relaxation, nowadays called finite differences [Ames; Golub and Ortega; Richmeyer and Morton]. Because the top and bottom faces are insulated, it is reasonable to assume no flow occurs in the vertical direction. Thus temperature within the block is of the form

      (1.21)

      where

      i) u(x, y, 0) = 1 for 0 < x < L, 0 < y < L,

      ii) u(0, y, t) = 0 = u(L, y, t) for 0 ≤ y ≤ L, t ≥ 0,

      iii) u(x, 0, t) = 0 = u(x, L, t) for 0 ≤ x ≤ L, t ≥ 0.

      The temporal condition i) is called an initial condition while the spatial conditions ii) and iii) are called boundary conditions.

      Partition the face 0 ≤ x, y ≤ L in the usual way into nx y = L/n indexed by their upper right-hand vertex (xi, yj) for i, j = 1, 2, . . . , n.

      Make the simplifying assumption that characterizes the method of finite differences: during the kt of time the temperature of the (i, j. (Traditionally time steps are indexed as superscripts while spatial indices are subscripts.)

      Figure 1.2 The five point rule for finite differences. During the kth interval of time, the center block loses heat to its four neighboring blocks; each block is imagined to be at a constant temperature throughout.

      that heat is lost from the (i, j)-th block to its four neighbors is proportional to the temperature differences across the four faces (see Figure 1.2); the constant of this proportionality is the conductivity κ scaled inversely by the distance traveled by the heat and directly by the surface area of the membrane through which it passes:

      Conductivity κ is in units of thermal power by thickness per degree temperature gradient per area.

      On the other hand, this total heat lost by the (i, j)-th block must be made up by a decrease in temperature during the k-th interval of time:

      where c is the specific heat of the material, the heat released by the material per unit mass per degree of temperature drop, and where ρ is the density of the material (mass per unit volume).

      Solving for the subsequent temperature of the subblock yields the explicit five point rule:

      (1.22)

      where α = κ/cρ is the diffusivity of the material (in units of area/time),

      The initial condition i) is imposed by requiring that

      = 1 for all 1 ≤ i, j ≤ n.

      and requiring

      for all j and k

      and

      for all i and k.

      The relation (1.22) together with conditions i)′, ii)′, and iii)′ are now implemented on a digital computer.

      During World War II such relaxation problems were solved by placing a large number of people with good computing skills at desks, each desk representing a cell. Upon a command, each desk would perform the right hand side of (1.22), then pass copies of the result to the four neighboring desks. Nowadays of course a digital computer is employed. See Figure 1.3.

      Figure 1.3 Snapshots of the relaxing temperatures of the block of Figure 1.1 with L = α = 1 (via the symbolic manipulator Mathematica). See Chapter 15.

      Warning. The explicit formula (1.22) suffers from instabilities unless the quantities α tx² and α ty² are small. Later we will discuss a stable implicit method known as the Crank-Nicolson Method that is now universally employed in such problems.

      The analytic model

      The equation (1.22) of relaxation can be rearranged into the form

      (1.23)

      xy → 0, it is reasonable that temperature u must satisfy in the interior of the square the partial differential equation (PDE)

      (1.24)

      an equation known as the heat equation or equation of diffusion. We will derive this famous equation by analytic means in §3.2.

      §1.4 The ∇ Operator

      Recall that the del operator

      (1.25)

      when applied to a scalar function f yields the gradient

      (1.26)

      a vector field, while in contrast the divergence transforms a vector field

      F = Pi + Qj + Rk

      into the scalar field

      (1.27)

      The curl

      (1.28)

      transforms vector fields to vector fields.

      As will be seen from the list of the ‘big six’ PDEs in §1.5, the most physically important operator is the Laplacian, the divergence of the gradient

      (1.29)

      transforming scalar fields to scalar fields. In particular, when f = f(x, y) and g = g(x),

      (1.30)

      The Laplacian in cylindrical coordinates (Exercise 1.19) is

      (1.31)

      The Laplacian in spherical coordinates (Exercise 1.20) is

      (1.32)

      §1.5 The Big Six PDEs

      The most frequently encountered distributed problems are modeled by variations on one of the Big Six PDEs:

      The heat equation (or equation of diffusion)

      (1.33)

      The wave equation

      (1.34)

      Laplace’s equation

      (1.35)

      Poisson’s equation

      (1.36)

      Schrödinger’s equation

      (1.37)

      The plate (beam) equation

      (1.38)

      Some of the above have been brought to their nondimensional form by an appropriate rescaling of variables.

      Exercises

      A function is harmonic on a domain Ω if it satisfies Laplace’s equation ∇²u = 0 on Ω.

      1.1 Prove the formula (1.3) for the directional derivative.

      1.2 Show that u = log r is harmonic on the punctured plane r > 0.

      1.3 Show that u = 1/ρ is harmonic on p > 0.

      Is v = 1/r harmonic on r > 0?

      1.4 = 0 within the unit disk x² + y² < 1. Hint: showing that ∇²u = 0 directly is difficult. I recommend either

      the use of a symbolic manipulator, or

      an algebraic reduction to a simpler problem, or

      noting that u is the real part of the analytic function f(z) = (1 + z)/(1 − z).

      1.5 Let v(x, y= u(ax + by, ay − bx) for a² + b² ≠ 0. This change of variables

      is a rotation and a dilation. Show using the chain rule that wherever sensible,

      (More generally, harmonicity is conformally invariant for functions in the plane.)

      1.6 Using the second chain rule, show that if α = x − ct and β = x + ct, then utt = c²uxx if and only if uαβ = 0.

      1.7 Show that the functions f and g of (1.9) are uniquely determined within constants that sum to zero.

      1.8 Solve ut + xux = 0 via the method of characteristics.

      1.9 What is the physical interpretation in (1.17) of its characteristic lines?

      1.10 Solve ux + uxy = 0.

      1.11 Solve ux + 2ut = 0 given that u(x, 0) = sin x.

      1.12 = ax + by = bx − ay, transport aux + buy = 0.

      1.13 Solve yux + xuy = 0.

      1.14 Solve uxx = 0 when u = u(x, y).

      1.15 Solve uxx = y²u when u = u(x, y).

      1.16 Argue that if ρ = ρ(x, t) is the density of freeway traffic flow at time t past point x and υ = v(x, t) its velocity, then ρt + (υρ)x = 0.

      1.17* Let u = u(x, t) be the population density of a country, i.e.,

      is the population at time t between the ages of a and b. Let q(x) be the probability of death less net immigration at age x per year. Argue for von Foerster’s model

      q(y) dy]. Argue that the left boundary condition is of the form

      and experiment with various initial populations, birth, and death rates. Is this model reasonable?

      1.18 Generically let f be a scalar field, F a vector field. Let ∇, ∇·, and ∇x denote the gradient, divergence, and curl respectively [Schey]. Of the 18 possible triple combinations of symbols ∇ * ∇ * * formed from the table

      only five make sense. Which five? Of the five, two are identically 0. Which two?

      1.19 Establish in two ways the Laplacian in polar coordinates

      a) first by applying the chain rule to the change of variables x = r cos θ and y = r sin θ, and then

      b) with physical intuition — rederive the heat equation by accounting for the flux leaving the four edges of r0 < r < rr and θ0 < θ < θθ. Take the limit.

      Corollary. (Laplacian in Cylindrical Coordinates)

      1.20 Establish that the Laplacian in Spherical Coordinates is

      in two ways:

      a) first via the chain rule on the change of variables x = ρ sin ϕ cos θ, y = ρ sin ϕ sin θ, z = ρ cos ϕ, then again

      b) via physical intuition — rederive the heat equation by accounting for the flux leaving the six faces of the box determined by spherical infinitesimals.

      1.21 Show using the form of the Laplacian in polar coordinates that u = θ is harmonic, i.e., ∇²u = 0 within any sector of opening less than 2π radians.

      1.22 Redo Exercise 1.3 using Exercises 1.19-1.20.

      1.23 Show that the units of diffusivity α must be area/time by two methods: first from examining α = κ/ρc and second by balancing units in the heat equation (1.23) itself.

      1.24 Find a vector field F with no potential, i.e., for which there is no scalar function ϕ such that ∇ϕ = −F.

      Hint: Give an infamous example like ‘dQ’ = cVdT + nRTdV/V or ‘∇θ.’

      1.25 Clausius’s version of the Second Law of Thermodynamics is that heat cannot of itself pass from a cold to a hot body. [Planck, §112]. Attempt to deduce that no portion of the block of §1.4 can be at a negative temperature during the above relaxation. More strongly, attempt to show that temperature at each location is a decreasing function of time.

      1.26 Argue that the dispersion of a spilled pollutant may in some circumstances be modeled by the three dimensional equation of diffusion

      1.27 Using a spreadsheet or MatLab, find the approximate temperatures of the block of Figure 1.1 at selected times 0 < t < 1 for the case α = L = 1 [Orvis].

      1.28 Argue that the thickness of the block of Figure 1.1 is irrelevant, that heat transfer within the block is modeled by (1.3) for any thickness.

      1.29 Consider a cylindrical rod of length L = 1 of homogeneous material of diffusivity α that is insulated everywhere but at the end faces. Initially the rod is uniformly at temperature 0, when suddenly the right end is put in contact with a source at temperature 1 while the left end is held at temperature 0. Derive Finite Difference and PDE models of the transient temperatures within the rod. Using a spreadsheet or MatLab, experimentally verify the intuitively clear eventual steady state. Experimentally discover for what value of h = α tx² your numerical model becomes unstable.

      1.30 Argue that the results of the previous Exercise could equally well hold for an infinite slab of thickness L.

      1.31 In the early days of the telegraph, undersea cables experienced what is today called intersymbol distortion: transmitted pulses would smear spatially as they traveled along the cable, combining with previous and subsequent pulses. William Thomson in 1854 postulated that the cable can be modeled as cascaded infinitesimal sections of distributed series resistance R (Ohms per meter) together with distributed shunt capacitance C (farads per meter), as shown in Figure 1.4. Recalling that voltage v and current i are related by the rules v = Ri and Cdv/dt = i for a resistor and capacitor respectively, reconstruct Kelvin’s deduction that the voltage V and current I on the cable must satisfy the system of coupled PDEs

      giving that the voltage V on the undersea cable is modeled by diffusion

      with diffusivity

      (The above problem is seminal for the field of Electrical

      Enjoying the preview?
      Page 1 of 1