

Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Material Type: Notes; Class: Mathematics in the Earth Sciences; Subject: Earth Sciences; University: University of California-Santa Cruz; Term: Unknown 1989;
Typology: Study notes
1 / 3
This page cannot be seen from the preview
Don't miss anything!


ES 111 Mathematical Methods in the Earth Sciences Lecture Outline 18 - Weds 26th Nov ’ Higher-Order Partial Differential Equations
The last thing we are going to mention is higher-order PDE’s. These are almost invariably found in real physical problems, and we can only touch on them here. We will simply look at one technique which can simplify their solution; taking solutions further involves techniques such as Fourier analysis not covered in this course.
Recall how we describe heat conduction. In one dimension, heat conduction is governed by the equation ∂T ∂t =^ κ
∂x^2 where κ is the thermal diffusivity and T (x, t) is a temperature field varying in x and time only. In 2 dimensions, the analogous equation may be written ∂T ∂t =^ κ
∂x^2 +^
∂z^2
) = κ∇^2 T
where now T is a function of x, z and t and ∇^2 is called the Laplacian operator, also referred to as ”del-squared”. The Laplacian can be written for 2,3... dimensions; it has a slightly different explicit form in curvilinear coordinate systems. The reason we use the Laplacian is that it is a useful shorthand, and that it works irrespective of which coordinate system we adopt. So ∂T∂t = κ∇^2 T is a true statement irrespective of whether we are using 2D, 3D, spherical, cylindrical etc. etc. Note, however, that the explicit form of the Laplacian does change depending on what coordinate system you’re in. Another place where the Laplacian arises is in (not surprisingly) Laplace’s equation: ∂^2 f ∂x^2 +^
∂^2 f ∂y^2 +^
∂^2 f ∂z^2 =^ ∇
(^2) f = 0
Here f is a scalar quantity f (x, y, z) and this equation is written in Cartesian coordinates. This equation is fundamental to the description of gravity and charge, among other phenomena.
The problem in dealing with a PDE is that the function we are interested in depends on several variables (e.g. T (x, t)). One trick that we can try, called Separation of Variables, is to assume that the two variables are independent. That is, we can rewrite T (x, t) as the product of a function that depends on x only, and another that depends on t only. This greatly simplifies matters. Example Say we want to solve the heat diffusion equation ∂T ∂t =^ κ
∂x^2 The we rewrite T (x, t) = A(x)B(t) and obtain
AdBdt = κB d
dx^2 Note that we have replaced partial derivatives with normal derivatives, because A and B are each functions of one variable only. We can divide through by AB to obtain 1 B
dB dt =^
κ A
d^2 A dx^2
Since the LHS is a function of t only and the RHS is a function of x only, then for this equality to be true, the expressions on both left and right must be constant. So we can write dB dt −^ kB^ = 0 and d^2 A dx^2 −^
k κA^ = 0 where k is an arbitrary constant. Note that it is sometimes more useful to employ −k or k^2 , depending on the exact problem being dealt with. Now we have two OLDE’s, which we can solve as usual. For example, a solution for B (remember that there are others) is B = B 0 ekt
where B 0 is another constant, while a solution to A (again, there are others) is
A = A 0 e(
√k/κ)x
and so one potential solution is that
T = T 0 ekte(
√k/κ)x