Singular solution

From HandWiki

A singular solution ys(x) of an ordinary differential equation is a solution that is singular or one for which the initial value problem (also called the Cauchy problem by some authors) fails to have a unique solution at some point on the solution. The set on which a solution is singular may be as small as a single point or as large as the full real line. Solutions which are singular in the sense that the initial value problem fails to have a unique solution need not be singular functions. In some cases, the term singular solution is used to mean a solution at which there is a failure of uniqueness to the initial value problem at every point on the curve. A singular solution in this stronger sense is often given as tangent to every solution from a family of solutions. By tangent we mean that there is a point x where ys(x) = yc(x) and y's(x) = y'c(x) where yc is a solution in a family of solutions parameterized by c. This means that the singular solution is the envelope of the family of solutions.

Usually, singular solutions appear in differential equations when there is a need to divide in a term that might be equal to zero. Therefore, when one is solving a differential equation and using division one must check what happens if the term is equal to zero, and whether it leads to a singular solution. The Picard–Lindelöf theorem, which gives sufficient conditions for unique solutions to exist, can be used to rule out the existence of singular solutions. Other theorems, such as the Peano existence theorem, give sufficient conditions for solutions to exist without necessarily being unique, which can allow for the existence of singular solutions.

A divergent solution

Consider the homogeneous linear ordinary differential equation

[math]\displaystyle{ xy'(x) +2y(x)= 0 , \,\! }[/math]

where primes denote derivatives with respect to x. The general solution to this equation is

[math]\displaystyle{ y(x)= C x^{-2} . \,\! }[/math]

For a given [math]\displaystyle{ C }[/math], this solution is smooth except at [math]\displaystyle{ x=0 }[/math] where the solution is divergent. Furthermore, for a given [math]\displaystyle{ x\not=0 }[/math], this is the unique solution going through [math]\displaystyle{ (x,y(x)) }[/math].

Failure of uniqueness

Consider the differential equation

[math]\displaystyle{ y'(x)^2 = 4y(x) . \,\! }[/math]

A one-parameter family of solutions to this equation is given by

[math]\displaystyle{ y_c(x) = (x-c)^2 . \,\! }[/math]

Another solution is given by

[math]\displaystyle{ y_s(x) = 0 . \,\! }[/math]

Since the equation being studied is a first-order equation, the initial conditions are the initial x and y values. By considering the two sets of solutions above, one can see that the solution fails to be unique when [math]\displaystyle{ y=0 }[/math]. (It can be shown that for [math]\displaystyle{ y\gt 0 }[/math] if a single branch of the square root is chosen, then there is a local solution which is unique using the Picard–Lindelöf theorem.) Thus, the solutions above are all singular solutions, in the sense that solution fails to be unique in a neighbourhood of one or more points. (Commonly, we say "uniqueness fails" at these points.) For the first set of solutions, uniqueness fails at one point, [math]\displaystyle{ x=c }[/math], and for the second solution, uniqueness fails at every value of [math]\displaystyle{ x }[/math]. Thus, the solution [math]\displaystyle{ y_s }[/math] is a singular solution in the stronger sense that uniqueness fails at every value of x. However, it is not a singular function since it and all its derivatives are continuous.

In this example, the solution [math]\displaystyle{ y_s(x)=0 }[/math] is the envelope of the family of solutions [math]\displaystyle{ y_c(x)=(x-c)^2 }[/math]. The solution [math]\displaystyle{ y_s }[/math] is tangent to every curve [math]\displaystyle{ y_c(x) }[/math] at the point [math]\displaystyle{ (c,0) }[/math].

The failure of uniqueness can be used to construct more solutions. These can be found by taking two constant [math]\displaystyle{ c_1 \lt c_2 }[/math] and defining a solution [math]\displaystyle{ y(x) }[/math] to be [math]\displaystyle{ (x-c_1)^2 }[/math] when [math]\displaystyle{ x \lt c_1 }[/math], to be [math]\displaystyle{ 0 }[/math] when [math]\displaystyle{ c_1\leq x\leq c_2 }[/math], and to be [math]\displaystyle{ (x-c_2)^2 }[/math] when [math]\displaystyle{ x \gt c_2 }[/math]. Direct calculation shows that this is a solution of the differential equation at every point, including [math]\displaystyle{ x=c_1 }[/math] and [math]\displaystyle{ x=c_2 }[/math]. Uniqueness fails for these solutions on the interval [math]\displaystyle{ c_1\leq x\leq c_2 }[/math], and the solutions are singular, in the sense that the second derivative fails to exist, at [math]\displaystyle{ x=c_1 }[/math] and [math]\displaystyle{ x=c_2 }[/math].

Further example of failure of uniqueness

The previous example might give the erroneous impression that failure of uniqueness is directly related to [math]\displaystyle{ y(x)=0 }[/math]. Failure of uniqueness can also be seen in the following example of a Clairaut's equation:

[math]\displaystyle{ y(x) = x \cdot y' + (y')^2 \,\! }[/math]

We write y' = p and then

[math]\displaystyle{ y(x) = x \cdot p + (p)^2. }[/math]

Now, we shall take the differential according to x:

[math]\displaystyle{ p = y' = p + x p' + 2 p p' }[/math]

which by simple algebra yields

[math]\displaystyle{ 0 = ( 2 p + x )p'. }[/math]

This condition is solved if 2p+x=0 or if p′=0.

If p' = 0 it means that y' = p = c = constant, and the general solution of this new equation is:

[math]\displaystyle{ y_c(x) = c \cdot x + c^2 }[/math]

where c is determined by the initial value.

If x + 2p = 0 then we get that p = −½x and substituting in the ODE gives

[math]\displaystyle{ y_s(x) = -\tfrac{1}{2}x^2 + (-\tfrac{1}{2}x)^2 = -\tfrac{1}{4} x^2. }[/math]

Now we shall check when these solutions are singular solutions. If two solutions intersect each other, that is, they both go through the same point (x,y), then there is a failure of uniqueness for a first-order ordinary differential equation. Thus, there will be a failure of uniqueness if a solution of the first form intersects the second solution.

The condition of intersection is : ys(x) = yc(x). We solve

[math]\displaystyle{ c \cdot x + c^2 = y_c(x) = y_s(x) = -\tfrac{1}{4} x^2 }[/math]

to find the intersection point, which is [math]\displaystyle{ (-2c , -c^2) }[/math].

We can verify that the curves are tangent at this point y's(x) = y'c(x). We calculate the derivatives:

[math]\displaystyle{ y_c'(-2 c) = c \,\! }[/math]
[math]\displaystyle{ y_s'(-2 c) = -\tfrac{1}{2} x |_{x = -2 c} = c. \,\! }[/math]

Hence,

[math]\displaystyle{ y_s(x) = -\tfrac{1}{4} \cdot x^2 \,\! }[/math]

is tangent to every member of the one-parameter family of solutions

[math]\displaystyle{ y_c(x) = c \cdot x + c^2 \,\! }[/math]

of this Clairaut equation:

[math]\displaystyle{ y(x) = x \cdot y' + (y')^2. \,\! }[/math]

See also

Bibliography