Partial differential equations solved by the variables separable method occupy perhaps 5 percent of the total constellation of such equations - but often provide much physical insight and to a wide variety of physics problems, from vibrating membranes, to electromagnetic waves and basic quantum mechanics. In this post we examine three examples of PDEs using the method and a few related problems are provided at the end.
Example 1:
Consider a wire of length ℓ :
]-------------- ℓ ---------------------->[
For which the relevant wave equation is:
¶ 2 u/ ¶ t 2 = c 2 [¶ 2 u/ ¶ x2 ]
And the displacement satisfies:
u(0, t) = u(ℓ , t ) = 0 ( t > 0)
Suppose at time t= 0 the displacement is u(x,0) = f(x)
And: ¶ u/ ¶ t ] t= 0 = g(x)
Using the technique of separation of variables provide two different general solutions for X( ℓ ) and X(x).
Hints: i) u(0, t) = X(0) T(t) and
(ii) u(ℓ , t )
= X( ℓ ) T(t)
Solution: By separation of variables we have:
u(x, t) = X(x) T(t)
Now let: dT/ dt = T ', d 2 T/ dt 2 = T"
dX/ dt = X ' d 2X/ dt 2 = X "
So that: ¶ 2 u/ ¶ t 2 = X(t) T "(t)
And: ¶ 2 u/ ¶ t 2 = X" (x) T (t)
The wave equation can then be written:
XT " = c 2 X" T
Or: T "/ c 2 T = X" / X
We set both sides equal to a constant K (separation constant):
T "/ c 2 T = K = X" / X
Two equations result:
i) X" - K X = 0
ii) T " - c 2 KT = 0
Look at two cases:
i) K = 0
Then X" - KX = 0 becomes X" = 0
For which: X(x) = c1 X + c2
And: X(0) = 0 Þ c2 = 0 so X(x) = c1 X
However, X(x) = 0 ⇒ c1 = 0 and X( ℓ ) = 0
Which is unacceptable, so K cannot be zero. So next look at K > 0
ii) K > 0
Then: Let
K = r 2
Has solution: c1 exp (r x) + c2 exp (-r x)
After eliminating c2 (as fn. of c1)
X(x) = c1 [exp ( r x) - exp (-r x) ] ⇒
Ñ 2 E = me ¶ 2 E / ¶ t 2
Solution:
In Cartesian coordinates, the H, E z-components partials can be written (with the angular frequency w , magnetic permeability m and the permittivity e ):
i) ¶ 2 E z / ¶ x2 + ¶ 2 E z / ¶ y 2 + g 2 E z = mew2 E z
ii) ¶ 2 H z / ¶ x2 + ¶ 2 H z / ¶ y 2 + g 2 H z = mew2 H z
The transverse mode (H z = 0 ) implies:
¶ 2 E z / ¶ x2 + ¶ 2 E z / ¶ y 2 + g 2 E z = - w2 me E z
We
assume: E z (x, y, z)
= E z (x, y) exp (- g z)
At time t = 0
Using a variables separable approach we try the solution:
E z (x, y, z) = X(x) Y(y) exp (- g z)
Þ Y d 2 X / d x2 +
X d 2 Y / d y 2 +
g 2 XY
=
- w2 me XY
Then divide by XY:
1/ X (d 2 X / d x2 ) + 1/Y (d 2 Y / d y 2 ) + g 2 + w2 me
And let h 2 = g 2 + w2 me
Then rewrite equation as:
1/ X (d 2 X / d x2 ) + 1/Y (d 2 Y / d y 2 ) + h 2 = 0
Þ 1/ X ( d 2 X / d x2 ) + h 2 = A 2 = 1/Y (d 2 Y / d y 2 )
X = C1 cos Bx + C2 sin Bx
Where: B 2 = A 2 + h 2
Similarly for Y: d 2 Y / d y 2 + A 2 Y = o
Where: Y = C3 cos Ay + C4 sin Ay
The
constants C1, C2, C3 and C4 are determined from boundary conditions.
Example 3:
3) Particle in a rectangular box, e.g.
For the particle shown find three independent solutions and the quantized energy given the PDE:
¶2 y/ ¶ x2 + ¶2 y/ ¶ y2 + ¶2 y/ ¶ z2 + 2mE/ ħ2 y = 0
And: 0 < x < a, 0 < y < b, 0 < z < c
This is
a 2nd order partial differential equation easily solved by the separation of
variables, e.g. :
y
=
X(x) Y(y) Z(z)
Then:
X’
= ¶X/ ¶ x Y’ = ¶Y/ ¶ y
and Z’ = ¶Z/ ¶ z
This
leads to the equation:
X”YZ
+ XY” Z + XYZ” + 2mE/ ħ2 XYZ
= 0
The
required normalization equation is then:
òoa òbo òc o ‖y‖ 2 dz dy dx = 1
Dividing
the equation by XYZ:
X”/ X +
Y”/ Y + Z”/Z + 2mE / ħ2 = 0
We let:
X”/ X
= - a2, Y”/ Y = -b2,
Z”/Z = -g2
With a, b and g constants.
Then, we have:
x2 + a2 x = 0, y2 + b2 y = 0, z2 + g2 z
So the independent solutions
will be:
x= Ö(2/a)
sin ( n xpx/a)]
y= Ö(2/b)
sin ( n y px/b)]
z= Ö(2/c)
sin ( n z px/c)]
where: n x = 1, 2, 3, 4.. etc.
Then:
y
(x, y, z) = (8/ abc)1/2 sin ( n xpx/a)
sin ( n y px/b) sin ( n zpx/c)
And further, to obtain
the quantized energy, E:
n2x p2 /a2
- n2y p2 /b2
- n2z p2 /c2
+ 2mE/ ħ2 =
0
And:
E = p2 ħ2 / 2m [n2x /a2
+ n2 y /b2 + n2 z /c2 )
Suggested Problems for the Math Whiz:
1) Find a partial differential equation whose solution is:
z = (a + 2) x + ( a2 + 1) y + b
2)Use separation of variables to solve the partial differential equation:
¶2 y/ ¶ x2 + ¶2 y/ ¶ y2 + ¶2 y/ ¶ z2 + 2mE/ ħ2 y = 0
For a particle in a cubic box. Include the quantized energy for the particle.
3) The partial differential equation for a deflected beam is given as:
¶ 2 u/ ¶ t 2 + c 2 [¶ 4 u/ ¶ x4 ] = 0
And: c 2 = EI/ r A
Where EI is the flexural rigidity, r is the density of the wood, and A is the area.
a) If u(x,t) is defined such that: t > 0 and
0 < x < ℓ
Use separation of variables to arrive at an initial general (but not optimal) solution and write the two ordinary differential equations arising from the approach.
b) Apply the characteristic equati0n (i.e. in order 4) to obtain improved general solutions X(x) and T(t) and also suggest values for the constants arising: C1, C2, C3 and C4.
No comments:
Post a Comment