Partial differential equations occupy a wide swath of physics and mathematical physics areas, from quantum mechanics, to plasma physics to the Lagrangian and Hamiltonian formulations of planetary motion and non-holonomic systems to magneto-hydrodynamics.
Perhaps the simplest form of such an equation is the 1st order partial:
¶ z/ ¶ y = (¶ z/ ¶ x
- 2) 2 + 1
For which we can use substitution and let:
¶ z/ ¶ x = p
= (a + 2)
And:
¶ z/ ¶ y = q
= a 2 + 1
Then general solution can be written:
z = (p) x +
(q) y + g
Or: z = (a + 2) x
+ (a 2 + 1) y +
f
¶ 2 u/ ¶ y ¶ x
= x 3 - y
Which we can rewrite in the form:
¶ y (¶ u/ ¶ x)
= x 3 - y
Then integrate partially with respect to y and holding x constant to get:
¶ u/ ¶ x = x 3 y -
y 2 /2 + j(x)
Next, integrate 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.
I point out at this point that while ordinary differential equations when integrated yield results with arbitrary constants, PDEs yield results with arbitrary functions.Homogeneous linear partial differential equations of 2nd order are also relatively tame to solve, and have the general form:
a(¶ 2 u/ ¶ x2 ) + b ( ¶ 2 u/ ¶ x ¶ y ) + c (¶ 2 u/ ¶ y2 ) = 0 (Original eqn.)
Where a, b, c are constants. And solutions will be of the form: f (y + mx)
Where f is an arbitrary function and m is a constant.
Then rewrite.:
¶ 2 u/ ¶ x2 = m2 f" (y + mx)
¶ 2 u/ ¶ x ¶ y = m f" (y + mx)
¶ 2 u/ ¶ y2 = f" (y + mx)
Substituting these into the original eqn.:
am2 f" (y + mx) + b m f" (y + mx) + cf" (y + mx) = 0
Or: f" (y + mx) [am2 + bm + c] = 0
Then f(y + mx) is a solution of the original equation if m satisfies the quadratic:
am2 + bm + c = 0
Note here that four conditions attach to the possible solutions:
i) a ≠ 0 and roots of quadratic are distinct
ii) a ≠ 0 and roots of quadratic are equal
iii) a = 0, b ≠ 0
iv) a = 0, b = 0, c ≠ 0
Another 2nd order homogeneous example but with the constant coefficient a:
¶ 2 z/ ¶ x2 -
a [ ¶ 2 z/ ¶ x ¶ y]
- 6a 2 [ ¶ 2 z/ ¶ y2]
= 0
To solve, first write the auxiliary equation, i.e. bearing the quadratic:
m2 - a m – 6a = 0
This will be found to have the roots: m = 3a, m = -2a
The general solution will therefore be:
z = j(y + 3 ax) + x y (y - 2 ax)
Suggested Problems For The Math Whiz:
1) Solve: ¶ z/ ¶ x = ax + y
2) For the 2nd order PDE:
am2 f" (y + mx) + b m f" (y + mx) + cf" (y + mx) = 0
Let the distinct roots be m1 and m2. What would the solutions be?
3) Same PDE as (2) but now let there be a double root, e.g. m o.
What would the solutions be?
4) Apply condition (iii) to the same equation what would the solution be?
5) Solve: ¶ 2 u/ ¶ x2 - 4 [ ¶ 2 u / ¶ x ¶ y] + 4[ ¶ 2 u/ ¶ y 2 ] = 0
No comments:
Post a Comment