Instant download of Separation of Variables for Partial Differential Equations An Eigenfunction Approach 1st Edition George Cain ebook PDF, every chapter
Instant download of Separation of Variables for Partial Differential Equations An Eigenfunction Approach 1st Edition George Cain ebook PDF, every chapter
https://ebookultra.com/download/partial-differential-equations-an-
introduction-2nd-edition-strauss-w-a/
https://ebookultra.com/download/handbook-of-differential-equations-
stationary-partial-differential-equations-volume-6-1st-edition-michel-
chipot/
https://ebookultra.com/download/partial-differential-equations-a-
unified-hilbert-space-approach-1st-edition-rainer-picard/
https://ebookultra.com/download/applied-partial-differential-
equations-john-ockendon/
Stochastic Partial Differential Equations 2nd Edition Chow
https://ebookultra.com/download/stochastic-partial-differential-
equations-2nd-edition-chow/
https://ebookultra.com/download/solution-techniques-for-elementary-
partial-differential-equations-third-edition-constanda/
https://ebookultra.com/download/introduction-to-partial-differential-
equations-rao-k-s/
https://ebookultra.com/download/control-theory-of-partial-
differential-equations-1st-edition-guenter-leugering/
https://ebookultra.com/download/partial-differential-equations-2nd-
edition-lawrence-c-evans/
Separation of Variables for Partial Differential Equations
An Eigenfunction Approach 1st Edition George Cain
Digital Instant Download
Author(s): George Cain, Gunter H. Meyer
ISBN(s): 9781584884200, 1584884207
Edition: 1
File Details: PDF, 9.72 MB
Year: 2005
Language: english
Separation
of Variables
for Partial
Differential
Equations
An Eigenfunction
Approach
George Cain
Georgia Institute of Technology
Atlanta, Georgia, USA
Gunter H. Meyer
Georgia Institute of Technology
A.r!anta, Georgia, USA
This book contains information obtained from authentic and highly regarded sources. Reprinted material is
quoted with permission, and sources are indicated. A wide variety of references are listed. Reasonable efforts
have been made to publish reliable data and infonnation, but the author and the publisher cannot assume
responsibility for the validity of all materials or for the consequences of their use.
No part of this book may be reprinted, reproduced, transmitted, or utilized in any form by any electronic.
mechanical, or other means, now known or hereafter invented, including photocopying, microfilming, and
recording, or in any information storage or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, please access www.copyright.com
(hltp://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC) 222 Rosewood Drive,
Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and registration
for a variety of users. For organizations that have been granted a photocopy license by the CCC, a separate
system of payment has been arranged.
Trademark Notice: Product or corporate names may be trademarks or registered trademarks. and are used only
for identification and explanation without intent to infringe.
Cain, George L.
Separation of variables for partial differential equations : an eigenfunction approach I George Cain,
Gunter H. Meyer.
p. cm. -- (Studies in advanced mathematics)
Includes bibliographical references and index.
ISBN 1-58488-420-7 (alk. paper)
l. Separation of variables. 2. Eigenfunctions. L Meyer, Gunter H. IL Title. III. Series.
QA377.C247 2005
5 l 5'.353--dc22 2005051950
\Ve would like to thank our editor Sunil Nair for welcoming the project and
for his willingness to stay with it as it changed its scope and missed promised
deadlines. ·
We also wish to express our gratitude to Ms. Annette Rohrs of the School
of Mathematics of Georgia Tech who transformed decidedly low-tech scribbles
into a polished manuscript. Without· her talents, and patience, we would· not
have completed the book.
Preface
vii
viii PREFACE
Chapter 8 deals with the potential equation. It describes how one can pr~
condition the data of problems with smooth solutions in order not to introdi.:ce
artificial discontinuities into the separation of variables solution. We solve P'=-
tential problems with various boundary conditions and conclude with a calcu-
lation of eigenfunctions for the two-dimensional Laplacian.
Chapter 9 uses the eigenfunctions of the preceding chapter to find eigenfunc-
tion expansion solutions of two- and three-dimensional heat, wave, and potenti.a::
equations.
This text is written for advanced undergraduate and graduate students :.::.
science and engineering with previous exposure to a course in engineering math-
ematics, but not necessarily separation of variables. Basic prerequisites beyond
calculus are familiarity with linear algebra, the concept of vector spaces of func-
tions, norms and inner products, the ability to solve linear inhomogeneous first
and second order ordinary differential equations, and some contact with practi-
cal applications of partial differential equations.
The book contains more material than can (and should) be taught in a cour~
on separation of variables. We have introduced the eigenfunction approach to
our own students based on an early version of this text. We covered parts of
Chapters 2-4 to lay the groundwork for an extensive discussion of Chapter 5.
The remainder of the term was filled by working through selected examples
involving the heat, wave, and potential equation. We believe that by term ·s
end the students had an appreciation that they could solve realistic problems.
Since we view Chapters 2-5 as suitable for teaching separation of variables,
we have included exercises to help deepen the reader's understanding of the
eigenfunction approach. The examples of Chapters 6-8 and their exercise sets
generally lend themselves for project assignments.
This text will put a bigger burden on the instructor to choose topics and
guide students than more elementary texts on separation of variables that start
with product solutions. The instructor who subscribes to the view put forth
in Chapter 5 should find this text workable. The more advanced applications,
such as interface, inverse, and multidimensional problems, as well as the the
more theoretical topics require more mathematical sophistication and may be
skipped without breaking continuity.
The book is also meant to serve as a reference text for the method of separa-
tion of variables. We hope the many examples will guide the reader in deciding
whether and how to apply the method to any given problem. The examples
should help in interpreting computed solutions, and should give insight into
those cases in which formal answers are useless because of lack of convergence
or unacceptable oscillations. Chapters 1 and 9 are included to support the
reference function. They do not include exercises.
We hasten to add that this text is not a complete reference book. We
do not attempt to characterize the equations and coordinate systems where
a separation of variables is applicable. We do not even mention the various
coordinate systems (beyond cartesian, polar, cylindrical, and spherical) in which
the Laplacian is separable. We _have not scoured the literature for new and
innovative applications of separation of variables. Moreover, the examples we do
PREFACE xi
include are often meant to show structure rather than represent reality because
in general little attention is given to the proper scaling of the equations.
There does not appear to exist any other source that could serve as a prac-
tical reference book for the practicing engineer or scientist. We hope this book
will alert the reader that separation of variables has more to offer than may be
apparent from elementary texts.
Finally, this text does not mention the implementation of our formulas and
calculations on the computer, or do we provide numerical algorithms or pro-
grams. Yet the text, and in particular our numerical examples, could not
have been presented without access to symbolic and numerical packages such as
!viaple, Mathematica, and Matlab. We consider our calculations and the graph-
ical representation of their results routine and well within the competence of
today's students and practitioners of science and engineering.
Contents
Acknowledgments v
Preface vii
3 Sturm-Liouville Problems 45
3.1 Sturm-Liouville problems for q/' = µ¢ 45
3.2 Sturm-Liouville problems for £¢ = µ¢ . . . 53
3.3 A Sturm-Liouville problem with an interface 59
4 Fourier Series 67
4.1 Introduction . 67
4.2 Convergence . 68
4.3 Convergence of Fourier series 74
4.4 Cosine and sine series . . . . 77
4.5 Operations on Fourier series . 80
4.6 Partial sums of the Fourier series and the Gibbs phenomenon 84
Bibliography 277
Index 279
Chapter 1
This chapter provides a quick look into the vast field of partial differential
equations. The main goal is to extract some qualitative results on the three
dominant equations of mathematical physics, the potential, heat, and wave
equation on which our attention will be focused throughout this text.
1.1 Overview
When processes that change smoothly with two or more independent variables
are modeled mathematically, then partial differential equations arise. Most
common are second order equations of the general form
M M
Lu= L aijUxix, +L biuxi +cu = F, xE DC IRM (1.1)
i,j=l i=l
where the coefficients and the source term may depend on the independent vari-
ables {x 1, ... , x M}, on u, and on its derivatives. D is a given set in !RM (whose
boundary will be denoted by [) D). The equation may reflect conservation and
balance laws, empirical relationships, or may be purely phenomenological. Its
solution is used to explain, predict, and control processes in a bewildering array
of applications ranging from heat, mass, and fluid flow, migration of biological
species, electrostatics, and molecular vibration to mortgage banking.
In (1.1) C. is known as a partial differential operator that maps a smooth
function u to the function C.u. Throughout this text a smooth function denotes
a function with as many continuous derivatives as are necessary to carry out
the operations to which it is subjected. lu = F is the equation to be solved.
Given a partial differential equation and side constraints on its solution,
typically initial and boundary conditions, it becomes a question of mathematical
analysis to establish whether the problem has a solution, whether the solution
2 CHAPTER 1. POTENTIAL, HEAT, AND WAVE EQUATION
is unique, and whether the solution changes continuously with the data of the
problem. If that is the case, then the given problem for (1.1) is said to be well
posed; if not then it is ill posed. We note here that the data of the problem are
the coefficients of£, the source term F, any side conditions imposed on u, and
the shape of D. However, dependence on the coefficients and on the shape of D
will be ignored. Only continuous dependence with respect to the source term
and the side conditions will define well posedness for our purposes.
The technical aspects of in what sense a function u solves the problem and
in what sense it changes with the data of the problem tend to be abstract and
complex and constitute the mathematical theory of partial differential equations
(e.g., [5]). Such theoretical studies are essential to establish that equation (1.1)
and its side conditions are a consistent description of the processes under consid-
eration and to characterize the behavior of its solution. Outside mathematics
the validity of a mathematical model is often taken on faith and its solution
is assumed to exist on "physical grounds." There the emphasis is entirely on
solving the equation, analytically if possible, or approximately and numerically
otherwise. Approximate solutions are the subject of this text.
so that
A= G ~)
We can now introduce three broad classes of differential equations.
is
i) Elliptic at x = (x 1 , ••. , XM) if all eigenvalues of the symmetric matrix A are
nonzero and have the same algebraic sign,
1.3. LAPLACE'S AND POISSON'S EQUATION 3
ii) Hyperbolic at i if all eigenvalues of A are nonzero and one has a different
algebraic sign from all others,
ili) Parabolic at x if A has a zero eigenvalue.
£u = \1 ·\Ju= \7 2 u = 6.u.
The last form is common in the mathematical literature and will be used con-
sistently throughout this text. The Laplacian in cartesian coordinates
M
~u. = LUxixi
i=l
f).u=F (1.2)
for a given source term Fis known as Poisson's equation. It will be the dominant
elliptic equation in this text.
In general equation (1.2) is to be solved for i! E D where D is an open set
in JRM. For the applications in this text D will usually be a bounded set with a
sufficiently smooth boundary 8D. On 8D the solution u has to satisfy boundary
conditions. We distinguish between three classes of boundary data for (1.2) and
its generalizations.
i) The Dirichlet problem (also known as a problem of the first kind)
u = g(i!), i! E 8D.
ii) The Neumann problem (also known as a problem of the second kind)
au
an=gx,
(-) i! E 8D
au + a2u = g (-)
°''on x ' i! E 8D,
where
b.u1 = 0, xED
U1 g(x),
= xE8D
b.u2 = F(x), xED
Uz = 0, xE8D
because different mathematical tools are available for Laplace's equation with
nonzero boundary data and for Poisson's equation with zero boundary data
(which, in an abstract sense, has a good deal in common with the matrix problem
Au= b). Splittings of this type will be used routinely in Chapter 8. Of course,
if g is defined and continuous on all of fJ and twice continuously differentiable
in the open set D, then it is usually advantageous to introduce the new function
w=u-g
b.w = F - b.g, xE D
w =0, xE8D
without splitting.
Given a classical solution we now wish to show that it depends continuously
on F and g. To give meaning to this phrase we need to be able to measure
change in the functions F and g. Here this will be done with respect to the so-
called supremum norm. We recall from analysis that for any function G defined
on a set D c RM
sup [G(i7)[ =least upper bound on the set of values {[G(i7)[ : if ED}.
fjED
A common notation is
llGI[ =sup [G(i7)[
fjED
which is called the supremum norm of G and which is just one example of the
concept of a norm discussed in Chapter 2.
If fJ is a closed and bounded set and G is continuous in fJ, then G must
take on its maximum and minimum on fJ so that
whenever
llFll + llgll '.": J.
Here
!lull =ma!' iu(x)I, llFll =sup IF(x)I, and llgll = xrp.Eaaxo lg(x)l-
XED :lED
Continuous dependence on the data in this sense, and uniqueness follow from
foe maximum principle for elliptic equations. Since it is used later on and always
provides a quick check on computed or approximate solutions of the Dirichlet
problem, and since it is basically just the second derivative test of elementary
calculus, we shall briefly discuss it here.
Assume that
F >0 for all x ED.
Then u(x) cannot assume a relative maximum in D.
Proof. Suppose that D is such that a < .xk < b for some k, 1 < k < M where
a and b denote finite lower and upper bounds on the kth coordinate of E D. x
Define the function
¢(£) = (llFll + <) (xk; a) 2
8 CHAPTER l. POTENTIAL, HEAT, AND WAVE EQUATION
where E > 0 is arbitrary and where llFll = supxED JF(x)I and 11911 = maxxEi:W l9(i
are assumed to be finite. Then
Proof. The difference between two solutions satisfies the Dirichlet problem
b..u = 0 in D
u = 0 on 8D
Corollary 1.4 Let F ~ 0. Then the solution of the Dirichlet problem assumes
its maximum on 8D.
(b - a) 2
lu,(x) - u(x)I :SE - 2- .
b..u=O, xED
1.3. LAPLACE'S AND POISSON'S EQUATION 9
xE D (1.5)
x E EJD.
The existence of UN is given because it will be found explicitly. Uniqueness of
rhe solution guarantees that no other solution of (1.5) exists. It only remains
to establish in what sense UN approximates the analytic solution u. But it is
clear from Theorem 1.2 that for all x ED
(1.6)
wher~ the constant K depends only on the geometry of D. Thus the error in
the approximation depends on how well FN and 9N approximate the given data
F and g. These issues are discussed in Chapters 3 and 4.
When the Dirichlet problem does not have a classical solution because the
data are not smooth, then continuous dependence for weak solutions must be
established. It generally is possible to show that if the data tend to zero in a
mean square sense, then the weak solution of the Dirichlet problem tends to
zero in a mean square sense. This translates into mean square convergence of
UN to u. The analysis of such problems becomes demanding and we refer to [6]
for details. A related result for the heat equation is discussed in Section 6.2.
Neumann problem: In contrast to the Dirichlet problem, the Neumann
problem for Poisson's equation is not well posed because if u is a solution then
u + c for any constant c is also a solution, hence a solution is not unique. But
there may not be a solution at all if the data are inconsistent. Suppose that
u is a solution of the Neumann problem; then it follows from the divergence
theorem that
= r
./!'JD
'llu(S). n(S)ds= r
./8D
g(S)ds,
where ii is the outward unit normal on EJD. Hence a necessary condition for the
existence of a solution is the compatibility condition
(1.7) implies that the solution does not change continuously with the data since.
F and g cannot be changed independently. If, however, (1.7) does hold, then the
Neumann problem is known to have a classical solution which is unique up to ,
an additive constant. In practice the solution is normalized by assigning a value
to the additive constant, e.g., by requiring u(x0 ) = 1 for some fixed x 0 E D.
The Neumann problems arise frequently in applications and can be solved with
separation of variables. This requires care in formulating the approximating
problem because it, too, must satisfy the compatibility condition (1.7) We shall
address these issues in Chapter 8.
We note that problems of the third kind formally include the Dirichlet and
Neumann problem. The examples considered later on are simply assumed (on
physical grounds) to be well posed. Uniqueness, however, is easy to show with
the maximum principle, provided
W = U1 - Uz
satisfies
~w =0, xE D
ow
0<1 on + 0<2W = 0, x E oD.
From the maximum principle we know that w must assume its maximum and
minimum.on oD. Suppose that w has a positive maximum at io E oD; then
the boundary condition implies that
ow
- (-
xo ) = --
°'2 w (-
xo ) < 0
on °'1
so that w has a strictly positive directional derivative along the inward unit
normal -ii. This contradicts that w(£0 ) is a maximum of w on D. Hence w
cannot have a positive maximum on D. An analogous argument rules out a
negative minimum so that w = 0 is the only possibility.
Finally, let us illustrate the danger of imposing the wrong kind of boundary
conditions on Laplace's equation.
For any positive integer k let us set
Ek = (2k + l)7r
u =0 on x = 0, 1
IA. THE HEAT EQUATION 11
and
u(x,O) = 0 for x E (0, 1)
while
lu G, I 1) = E~ sinh ~ _, oo ask-. oo.
Hence the boundary data tend to zero uniformly while the solution blows up.
Thus the problem cannot be well posed.
In general it is very dangerous to impose simultaneously Dirichlet and Neu-
mann data (called Cauchy data) on the solution of an elliptic problem on a
portion of 8D even if the application does furnish such data. The resulting
problem, even if formally solvable, tends to hav: an unstable solution.
for a1, 0<2 ;::: 0 and a1 + a2 > 0. Let us define the set
Qr={(x,t):xED, O<t'.'OT},
where T > 0 is an arbitrary but fixed final time. Then a solution of (1.9) has
to satisfy (in some sense) (1.9a) in Qr and the boundary and initial conditions
on &Qr.
A classical solution of (1.9) is a function which is smooth in Qr, which is
continuous on Qr U &Qr (together with its spacial derivatives if a 1 f= 0) and
which satisfies the given data at every point of &Qr.
It is common for diffusion problems that the initial and boundary condi-
tions are not continuous at all points of the parabolic boundary. In this case
u cannot be continuous on &Qr and one again has to accept suitably defined
weak solutions which only are required to solve the diffusion equation and ini-
tial/boundary conditions in an integral equation sense.
First and foremost in the discussion of well posedness for the heat equation
is the question of existence of a solution. For the one-dimensional heat equation
one can sometimes exhibit and analyze a solution in terms of an exponential
integral - see the discussion of (1.13), (1.14) below - but in generirl~is
question is resolved with fairly abstract classical and functional analysis. As in
the case of Poisson's equation it is possible to split the problem by writing
where
(x,t) E Qr
1.4. THE HEAT EQUATION 13
°'1
8u1 + 0<2U1 =
Bn (- t ) ,
g X, xEaD, 0<t:5.T
u1(x,O) = 0, x E [J
and
.Cu2 = F(x, t), (x,t) E Qr
xECID, 0<t:5.T
u2(x, O) = uo(x),
and employ special techniques to establish the existence of u 1 and u 2 . In par-
ticular, the problem for Uz in an abstract sense has a lot in common with an
n-dimensional first order system
du
dt - A(t)u = G(t), u(O) = uo
.Cu= F, (x,t) E Qr .
Proof. If u had a maximum at some point (£0 , t 0 ) E Qr, i.e., £ 0 lies in the
open set D, then necessarily
ux,x,(£0 , t 0 ) :5. 0 and u,(£0 , to)= 0 if to< Tor u,(xo, to) 2: 0 if to= T.
14 CHAPTER 1. POTENTIAL, HEAT, AND WAVE EQUATION
fluff :Sf
ff Pff =sup fF(x, t)f, ffgff = ~~x fg(x, t)f, ffuoff = m11x fuo(x)f.
QT D
O:St~T
Theorem 1.6 The solution of the Dirichlet problem for (1.8) depends contin-
uously on the data F, g, and u 0 .
Proof. We again assume that the kth coordinate Xk satisfies a :S Xk :S b for all
x ED. For arbitrary f > 0 define
(b a) 2
fluff :S ffFff ~ + max{ffgff, ffuoff}. (1.10)
~imilarly, we can conclude that the solution of the Dirichlet problem for the
heat equation is unique and that the solution of the heat equation with F 0 =
(the homogeneous heat equation) must take on-its maximum and minimum on
iJQy.
In addition to boundary value problems we also can consider a pure initial
value problem for the heat equation
u(x, O) = uo(x).
It can be verified that the problem is solved by the formula
u(x, t) = J.
IRM
s(x - fj, t)u0 (fj)dfj (1.11)
where
- 1 _jg!
s(x, t) = (47rt)M/2 e " -
Here (x,x) denotes the dot product for vectors in IRM- s(x, t) is known as the
fundamental solution of the heat equation. A simple calculation shows that s is
infinitely differentiable with respect to each component X; and t for t > 0 and
that
!is(x, t) - s,(x, t) = 0 fort> 0.
It follows from (Lll) that u(x, t) is infinitely differentiable with respect to all
variables for t > 0 provided only that the resulting integrals remain defined
and bounded. In particular, if u0 is a bounded piecewise continuous function
defined on IRM, then the solution to the initial value problem exists and is
infinitely differentiable for all X; and all t > 0. It is harder to show that u(x, t)
is continuous at (£0 , 0) at all points x 0 where u 0 is continuous and that u(x, t)
assumes the initial value u 0 (x0 ) as (x, t) _, (x0 , 0). We refer to [5) for a proof
of these results. Note that for discontinuous u 0 the expression (1.ll) is only a
weak solution because u(x, t) is not continuous at t = 0.
We see from (1.11) that if for any E > 0
( -) _ {1 for \\x\\::; E
uo x - 0 otherwise,
then u(x, t) > 0 for t > 0 at all x E !RM- In other words, the initial condition
spreads throughout space infinitely fast. This property is a consequence of the
mathematical model and contradicts the observation that heat does not flow
infinitely fast. But in fact, the change in the solution (1.11) at \\x\\ » E remains
unmeasurably small for a certain time interval before a detectable heat wave
arrives so that defacto the wave speed is finite. We shall examine this issue at
length in Example 6.3 where the speed of an isotherm is found numerically.
The setting of diffusion in all of IRM would seem to preclude the application of
(1.ll) to practical problems such as heat flow in a slab or bar. But (1.11) is not
as restrictive as it might appear. This is easily demonstrated if M = 1. Suppose
that u 0 is odd with respect to a given point x 0 • i.e., u 0 (x 0 + x) = -u 0 (x 0 - x);
u(x 0 1:
then with the obvious changes of variables
= -1
00
s(x - z, t)uo(xo - z)dz = 1-oo s(x - x 0 + y)u 0 (y)dy
=- 1:
-00
- x, t)
16 CHAPTER 1. POTENTIAL, HEAT, AND WAVE EQUATION,
we see that u(x, t) is odd in x with respect to x 0 for all t. Since u(x, t) is
smooth for t > 0, this implies that u(x 0 , t) = 0. Furthermore, if uo is periodic
with period w, then a similar change of variables techniques establishes that
u(x, t) is periodic in x with period w. Hence if for an integer n
. mrx
u 0 (x ) =smy,
for scaJars {&n}, then uN(x, t) given by (1.ll) is the unique classical solution
of the initial/boundary value problem
LU = Uxx - Ut = 0 - (1.13)
u(O, t) = u(L, t) = 0
u(x, 0) = ua(x).
= U(L, t, r) = 0
U(O, t, r)
U(x,r,r) = -F(x,r),
where r is a nonnegative parameter, then
£u = b.u + u, = 0.
u(x,O) =u 0 (£).
Suppose we wish to find u(:i', T) for T > 0. If we set r = T - t and w(x, r) =
u(x, T - r), then the problem is equivalent to finding w(x, 0) of the problem
£w =bow - w, = 0
w(x, T) = uo(x).
In a thermal setting this implies that from_ knowledge of the temperature at
some future time T we wish to find the temperature today. Intuition tells us
that if T is large and u 0 is near a steady-state temperature, then all initial
temperatures w(x, 0) will decay to near u 0 . In other words, small changes in u 0
could be consistent with large changes in w(£, 0), suggesting that the problem
is not well posed when the heat equation is integrated backward in time (or the
backward heat equation is integrated forward in time).
18 CHAPTER l. POTENTIAL, HEAT, AND WAVE EQUATION
The equation is usually associated with oscillatory phenomena and shows mark-
edly different properties compared to Poisson's and the heat equation. Equation
(1.16) is an example of a hyperbolic equation. Here the (M + 1) x (M + 1) matrix
A has the form
A= c~ !1)
where IM is the M-dimensional i<lentity matrix.
It is easy to show that (1.16) allows wave-like solutions. For example, let f
be an arbitrary twice continuously differentiable function of a scalar variable y.
Let ii be a (Euclidean) unit vector in !RM and define
y = (ii,x) - ct
where (ii, x) is the dot product of n and x. Then differentiation shows that
is a solution of (1.16). Suppose that c,t > 0, then the set {x : (n,x) -
ct = constant} is a plane in !RM traveling in the direction of ii with speed
J.5. THE WAVE EQUATION 19
c, and f( (ii., x) - ct) describes a wave with constant value on this plane. For
example, if
j(y) = eiY,
then
/((ii., x) - ct) = ei((n,x)-ct]
is known as a plane wave. Similarly, g( (ii., x) +ct) describes a wave traveling in
the direction of -n with speed c.
Our aim is to discuss again what constitutes well posed problems for (1.16).
We begin by exhibiting a solution which is somewhat analogous to the solution
of the one-dimensional heat equation discussed at the end of Section 1.4.
Ifwe set x = (ii.,x), then the solutions f(x-ct) and g(x+ct) of (1.16) solve
the one-dimensional wave equation
1
£u '=' Uxx - ;;2 Utt = 0 (1.17)
which, for example, describes the motion of a vibrating uniform string. Here
u(x, t) is the vertical displacement of the string from its equilibrium position.
We show next that any smooth solution of _(1.17) must be of the form
i.e., the superposition of a right and left traveling wave. This observation follows
if we introduce new variables
E=x-ct
1)=x+ct
and express the wave equation in the new variables. The chain rule shows that
and determine f and g so that u satisfies the initial conditions. Hence we need
-11"
c xo
ui(s)ds=-f(x)+g(x)+K
where xo is some arbitrary but fixed point in (-oo, oo) and f( = f(xo) - g(xo).
When we solve algebraically for f(x) and g(x), we obtain
f(x) = ~ [uo(x) -
.
~
C
{" u1(s)ds +
lxo
K]
Hence
which simplifies to
u(x, t) =
1
2 [uo(x - ct)+ uo(x +ct)]+
1
2C
;·x+ct u (s)ds.
1 (1.19)
x-ct
The expression (1.19) is known as d'Alembert's solution for the initial value
problem of the one-dimensional wave equation. If u 0 is twice continuously dif-
ferentiable and u 1 is once continuonsly differentiable on ( -oo, oo), then the
d' Alembert solution is a classical solution of the initial value problem for all fi-
nite t and x. Moreover, it is unique because u has to be the superposition of two
traveling waves and the d'Alembert construction determines f and g uniquely.
Moreover, if we set
then
lu(x, t)i ~ lluoll + tllu1 I
which implies continuous dependence for all t ~ T where T is an arbitrary but
fixed time. Hence the initial value problem (1.17), (l.18) is well posed.
We observe from (1.19) that the value of u(x 0 , t 0 ) at a given point (x 0 , to)
depends only on the initial value u 0 at x 0 - ct 0 and x 0 + ct 0 and on the initial
value u 1 over the interval [xa - eta, xa +eta]. Thus, if
l.5. THE WAVE EQUATION 21
and xo > t, then regardless of the form of the data on the set lxl < e we have
u(x 0 , t) = {o ,
fort<=
c
..!..
2c J
-£
u 1 (s)ds for t > ~
c .
Hence these initial conditions travel with speed c to the point Xo but in general
~(xo, t) will not decay to zero as t ..... oo. This is a peculiarity of the M-
.iimensional wave equation for M = 1 and all even M [5]. If uo and u 1 do
not have the required derivatives but (1.19) remains well defined, then (1.19)
represents a weak solution of (1.17), (1.18). We shall comment on this aspect
when discussing a plucked string in Example 7.1.
As in the case of the heat kernel solution we can exploit symmetry properties
of the initial conditions to solve certain initial/boundary value problems for
the one-dimensional wave equation with d'Alembert's solution. For example,
suppose that u 0 and u 1 are smooth and odd with respect to the point x 0 ; then
the d'Alembert solution is odd with respect to x 0 •
To see this suppose that
and that u 1 is odd with respect to the point x 0 , i.e., u 1 (x 0 +x) = -u 1 (x 0 -x).
Then with y = xo - s and r = Xo + y we obtain
Since by hypothesis
we conclude that the d' Alembert solution is odd with respect to the point x 0 and
hence equal to zero at x 0 for all t. It follows that the boundary value problem
I
lu '= Uxx - ;::2" Utt = 0
u(O, t) = u(L, t) = 0
N
u(x, 0) = uo(x) = L &n sin n~x
n=l
N
u,(x, 0) = u1(x) = "'\"""' • mr
~ f3n sin L
n=l
22 CHAPTER 1. POTENTIAL, HEAT, AND WAVE EQUATION
for constant {&n} and {~n} is solved by the d'Alembert solution (1.19) because
u 0 an:d u 1 a.re defined on ( -oo, oo) and odd with respect to x = 0 and x = L.
As we remarked in Section 1.4, inhomogeneous boundary conditions can
often be made homogeneous at the expense of adding a source term to the
differential equation. This leads to problems of the type
1
LU= Uxx - - Utt = F(x, t)
C2
u(O, t) = u(L, t) = 0
u(x, 0) = u,(x, 0) = 0.
Now a Duhamel superposition principle can be applied. It is straightforward to
l
show that the function
u(x,t) = U(x,t,T)dT
1
[,U = Uxx - ;::2 Utt = 0
U(O, t, T) = U(L, t, T) = 0
U(x,T,T)=O
U,(x,T,T) = -c2 F(x,T)
where T is a parameter. It follows that if F is of the form
U(x, t, T) = - -
C
2
1 x+c(t-T)
:L-C(l-T)
FN(s, T)ds.
All integrations can be carried out analytically and the resultant solution UN (x, t)
can be shown to be identical to the separation of variables solution found in
Chapter 7 when an arbitrary source term F is approximated by a trigonometric
sumFN.
Let us now turn to the general initial/boundary value problem of the form
u(i, 0) = uo(i)
u,(i,0) = u 1 (i)
J.5. THE WAVE EQUATION 23
au + 0<2U = g (-X, t ),
0<1 [}n x E BD, t>0
where D is a given domain in !RM. For convenience we have set c = 1 which can
always be achieved by scaling time. We point out that if D = !RM and we have
a pure initial value problem, then it again is possible to give a formula for u(x, t)
analogous to the d'Alembert solution of the one-dimensional problem (see [5,
Chapter 2]), but for a true initial/boundary value problem the existence of a
solution will generally be based on abstract theory. (We note in this context
=
that if g 0 on [JD, then, again in a very general sense, the problem has a lot
in common with the ordinary differential equation
cf2u
dt 2 - A(t)u = G(t), u(O) = u0 , u'(O) = ui-)
We shall henceforth assume that the existence theory of [5] applies so that we
can concentrate on uniqueness of the solution and on its continuous dependence
on the data of the problem.
Uniqueness and continuous dependence follow from a so-called energy esti-
mate. If u is a smooth solution of (1.20) with g 0 on BD, then =
u,.Cu = Ut/:;.U - UtUtt = u,F(x, t)
so that
l [V' · u, V'u - V'u, · '\i'u - u,u,,] dx = l u,F di.
imply
so that
+ n2)u, Bu -0<2u,u - Bu (Bu) .
(n1
an = °'1
an an t
Since IUtl' I :S ¥- + t;-, we obtain from (1.21) the estimate
(1.22)
where
24 CHAPTER 1. POTENTIAL, HEAT, AND WAVE EQUATION
and
llF(.,t)ll = (L F (x,t)dxf
2
2
where
F = g = uo = u1 = 0.
This implies that E(O) = 0 so that (1.23) assures that E(t) = 0 for all t. Then
by Schwarz's inequality (see Theorem 2.4)
Basic Approximation
Theory
This chapter will review the abstract ideas of approximation that will be used
in the sequel. Let X be a linear space (sometimes called a vector space) over the
field S of real or complex numbers. The elements of X ~re called vectors and
those of S are called scalars. The vector spaces appearing in this book are Fl,,
and Cn with elements x = (x1, ... , Xn), jj, etc., or spaces of real- or complex-
valued functions/, g, etc. d-efined on a real interval or, more generally, a subset
of Euclidean n-space. In all spaces 0 denotes the zero vector. The scalars of S
are denoted by a., /3, or a, b, etc.
We now recall a few definitions from linear algebra which are central in our
discussion of the approximation of functions.
25
26 CHAPTER 2. BASIC APPROXIMATION THEORY
Example 2.3 a) In real Euclidean n-space Rn, the usual dot product
where x = (x1, x2, ... , Xn) and iJ = (Y1, y2, ... , Yn) is an inner product.
b) In R2, for x = (x 1,x2) and iJ = (y 1,yz) define (x,fj) by (x,fj) =Ax· i],
where A is the matrix
A=(2 1)
. 1 2
and i1 · iJ is the usual dot product of i1 and ii. Then (x, fj) is an inner product.
First consider
It is clear that (x, x) :'.'.: 0 and (x, x) = 0 if and only if x = (0, 0).
To see that (x, flJ = (y, x), simply compute both inner products. The re-
maining two properties are evident.
c) On the space of all continuous functions (real- or complex-valued) defined
on the reals having period 2£, it is easy to verify that
L
(f,g) = J f(t)g(t)dt
-L
is an inner product.
Now
(f + ag, f + ag) = (!, f) +Ci(!, g) + a(g, f) + lnl 2 (g, g).
Next, let a= t(f,g), where tis any real number. Then-·
This expression is quadratic in t and so the fact that it is never negative means
that
4 l(f, g)l 4 - 4 l(f,g)l 2 (!, !) (g,g) ::; 0.
In other words,
l(f,g)l 2 ::; (!,!) (g,g),
which completes the proof.
The inequality
l(f,g)I < J(f,f) ~
is known as Schwarz's inequality.
Proof. The proofs that F(f) :'.'.: 0 and F( a!) = lnl F(f) are simple and omitted.
We prove the triangle inequality.
And now happened the event for which Archie was looking and
for which he was prepared. The other horses crowded in upon the
bay—his own would have followed had he not restrained him—and
in less time than it takes to tell it the spring was full of rearing,
kicking, biting animals. But they were not permitted to remain there
long. The leader began to assert his authority. Raising his head with
an angry snort, he sent his heels right and left to such good purpose
that he quickly cleared a space about him. The instant that his head
was elevated was enough for Archie. The lariat left his hand, true to
its aim this time, and settling down over the bay’s neck was quickly
drawn tight. Finding himself fast, the horse started to run, but the
other end of the lasso was securely fastened to Archie’s saddle, and
his own horse bracing himself to meet the shock, the captured steed
was thrown flat on his side. The others snorted with terror and took
to their heels in short order.
When Fred and Eugene arrived upon the ground they saw
something exciting.
CHAPTER VIII.
THE EMIGRANT TRAIN.
“I HEARD Dick talking about Pikes the other day,” said Eugene,
who had followed Archie in order to lend his assistance in case
the bay had succeeded in freeing himself from his halter or lariat.
“Did he refer to such people as our friends here?”
“He did,” answered Archie. “As I understand it, these Pikes are
one of the results of our late war. A good many of the people in the
border states were completely broken up, their houses being
burned, their crops destroyed and their stock killed or driven off. It
takes time and patience to accumulate property by farming, and it is
hard for a man to begin over again where he began twenty or forty
years ago. There are not many who have the courage to do it. These
people had heard wonderful stories of the mines and the sudden
fortunes sometimes made there, and believing that since they had to
dig in the ground for a living, it would be easier to dig gold and
silver than potatoes, they emigrated by regiments. But where one
family made a permanent settlement in California, a dozen turned
back. But still some of them were not satisfied, and after staying a
year or two in their old homes, they would bundle up again and go
back to California; and from roaming about so much, they fell into
the habit of leading a wandering, gipsy life. They are not satisfied
with a permanent settlement anywhere. How the name they bear
came to be applied to them I don’t know, unless it was because the
first of them came from Pike county, Missouri. This man is a genuine
Pike. There is no telling how many times he has wandered back and
forth over these prairies, but he hasn’t learned much during his
wanderings. I have crossed the prairie only once, and I know that he
is a hundred miles north of where he belongs.”
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
ebookultra.com