Characteristic equation of a PDE

Probably is a particular case of method of characteristics...

Documento 331_2.jpg

The following discussion is made for first order homogeneous linear PDEs, but according to this entry of Wikipedia it can be generalized to any first order PDE (see section Fully nonlinear case).

Given a first order homogeneous linear PDE

(1)a1W/x1++anW/xn=0

where ak=ak(x1,,xn), we look for solutions

W:RnR.

We can consider a solution like "surfaces" S given implicitly by W=cte. Their tangent planes are given by the kernel of dW.

If we define the vector field A=a1x1+a2x2++anxn, a function W is solution of (1) if A(W)=0, i.e., the integral curves of A are contained in the surfaces W=cte.

The ODE

(2)x˙=A(x)

is called the characteristic system or characteristic equation, and their integral curves are called the characteristic curves of the original PDE. Solutions of the PDE are first integrals of the characteristic equation.

The ODE (2) is, for me, like a parametric expression of a line r in linear algebra (it is the case when A(x) is a constant vector field). On the contrary, equation (1) is similar to an implicit expression of all the planes containing the line r. The 1-form dW is like an implicit equation

Method of characteristic
I think that method of characteristics is usually called to what follows.
The line r above can be described by saying that it is made of vectors such that they are proportional to (a1,,an), that is,

dx1a1=dx2a2=

We can think of this expression as the symmetric form of a line in space ("expresión continua de la recta en el espacio" in Spanish).

From here, in the case of linear algebra, we recover several planes whose intersection gives the line r:

a2x1a1x2=0a3x2a2x3=0...

In the differential context, we could recover an independent set of first integral if we were able to solve

a2dx1a1dx2=0a3dx2a2dx3=0...

The solutions are independent first integrals of the ODE or functions that generate functionally the solutions of the PDE.

In the general case, i.e. ai depending on too many variables, the system cannot be solved, although is useful for numerical approximation of initial or boundary value problems.

Particular case
If in the PDE we look for solutions in the form of explicit functions

W=xnf(x1,,xn1)

we get

a1W/x1++an1W/xn1=an

which is the case for inhomogeneous linear PDEs.