This post will deal with the solutions of some partial differential equations starting with relatively simple ones.
Examples:
Ex. (1) Solving a first order partial:
¶ z/ ¶ y = (¶ z/ ¶ x - 2) 2 + 1
Let: ¶ z/ ¶ x = p = (a + 2)
Let: ¶ z/ ¶ y = q = a 2 + 1
So: q = (p - 2) 2 + 1
Then general solution can be written:
z = (p) x + (q) y + g
Or: z = (a + 2) x + (a 2 + 1) y + f
Ex. (2) Solving a 2nd order partial
¶ 2 u/ ¶ y ¶ x = x 3 - y
First, write in the form:
¶ y (¶ u/ ¶ x) = x 3 - y
Next integrate partially with respect to y, holding x constant to obtain:
¶ u/ ¶ x = x 3 y - y 2 /2 + j (x)
Next integrate the result partially with respect to x, holding y constant to obtain the solution in the form:
u = x 4 y /4 - x 2 y/ 2 + f(x) + g(y)
Where f(x) = ò j (x) dx
is an arbitrary function of x and g is an arbitrary function of y. We note here that while ordinary DEs have solutions which contain arbitrary constants, partial differential DE solutions have arbitrary functions.
Ex. (3) 2nd order (homogenous with constant coefficients)
¶ 2 z/ ¶ x2 - a [ ¶ 2 z/ ¶ x ¶ y] - 6a 2 [ ¶ 2 z/ ¶ y2] = 0
First write the auxiliary equation:
m2 - a m – 6a = 0
Which has roots: m = 3a, m = -2a
The general solution will therefore be:
z = j (y + 3 ax) + x y (y - 2 ax)
Ex. (4): Application in three dimensions, using separation of variables.
A particle confined in a 3-dimensional Box
We want to find the applicable wave function solution ( y (x,y,z)), and the energy applicable
We can set out the rectangular regions (limits) according to:
0 < x < a
0 < y < b
0 < z < c
The partial differential equation to solve then becomes:
¶2 y/ ¶ x2 + ¶2 y/ ¶ y2 + ¶2 y/ ¶ z2 + 2mE/ ħ2 y = 0
This is a 2nd order partial differential equation easily solved by the separation of variables:
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
Then, 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, 3 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:
n2 x p2 /a2 - n2 y p2 /b2 - n2 z p2 /c2 + 2mE/ ħ2 = 0
And:
E = p2 ħ2 / 2m [n2x /a2 + n2 y /b2 + n2 z /c2 )
If the box is a cube, so a = b = c:
E = p2 ħ2 / 2m a2 (n2x + n2 y + n2 z )
Problems for the Math Whiz:
1) Solve: ¶ z/ ¶ x = ax + y
2) Solve: ¶ 2 u/ ¶ x2 - 4 [ ¶ 2 u / ¶ x ¶ y] + 4[ ¶ 2 u/ ¶ y 2 ] = 0
3) Solve the following higher order partial DE:
¶ 4 z/ ¶ x4 - 6 [ ¶ 4 z/ ¶ x 3 ¶ y] + 14 [¶ 4 z/ ¶ x2 ¶ y2 ] -
16 [ ¶ 4 z/ ¶ x ¶ y 3 ] + 8 [ ¶ 4 z/ ¶ y 4 ] = 0
3) Consider a wire of length ℓ :
]-------------- ℓ ---------------------->[
for which the relevant wave equation is:
¶ 2 u/ ¶ t 2 = c 2 [¶ 2 u/ ¶ x2 ]
If 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 write two different general solutions for X( ℓ ) and state why one is unacceptable
Hints: i) u(0, t) = X(0) T(t) and
(ii) u(ℓ , t ) = X( ℓ ) T(t)
No comments:
Post a Comment