OU blog

Personal Blogs

Valentin Fadeev

Pfaff equation

Visible to anyone in the world
Edited by Valentin Fadeev, Wednesday, 15 Sep 2010, 23:37

Moving on with Stepanov's book I have reached the subject equations which have the following form (3 variables):

sum with, 3 , summands cap p times d times x plus cap q times d times y plus cap r times d times z equals zero

Where P, Q, R are sufficiently differentiable functions of x,y,z.

Excercise 205:

open y times z minus z squared close times d times x minus x times z times d times y plus x times y times d times z equals zero

z times open y minus z close times d times x plus x times open y times d times z minus z times d times y close equals zero

Integrating factor: mu equals one divided by y squared

z divided by y times open one minus z divided by y close times d times x plus x times y times d times z minus z times d times y divided by y squared equals zero

z divided by y times open one minus z divided by y close times d times x plus x times d of z divided by y equals zero

d times x divided by x plus d of z divided by y divided by z divided by y times open one minus z divided by y close equals zero

z divided by y equals u

d times x divided by x plus d times u divided by u minus d times u divided by u minus one equals zero

x times u divided by u minus one equals cap c

x times z divided by z minus y equals cap c

Will we always be lucky to have an appropriate factor to cast the equation into a full differential form? The book gives a negative answer setting out a very specific condition on the coefficients.

Assume the equation does have a solution and this solution is a 2-dimensional manifold, i.e has the form:

normal cap phi of x comma y comma z equals cap c

or (locally, at least):

z equals phi of x comma y

Then

equation left hand side d times z equals right hand side normal partial differential times z divided by normal partial differential times x times d times x plus normal partial differential times z divided by normal partial differential times y times d times y

On the other hand, by virtue of the equation (assuming R does not vanish identically):

equation left hand side d times z equals right hand side negative cap p divided by cap r times d times x minus cap q divided by cap r times d times y

Comparing the coefficients:

equation sequence normal partial differential times z divided by normal partial differential times x equals negative cap p divided by cap r equals cap a

equation sequence normal partial differential times z divided by normal partial differential times y equals negative cap q divided by cap r equals cap b

This is an overdetermined system: one function, two equations which generally does not have a solution The integrability condition can be obtained by equation mixed second derivatives, however I will quote the geometrical argument which may also shed light on some fact presented below.

Consider an infinitesimal shift along the manifold from open x comma y comma z close to open x plus d times x comma y close . Then z will take value:

equation left hand side z plus normal partial differential times z divided by normal partial differential times x times d times x equals right hand side z plus cap a times d times x

From here we move to the point with coordinates open x plus d times x comma y plus d times y close without leaving the manifold. New value of z:

equation left hand side sum with, 3 , summands z plus cap a times d times x plus normal partial differential divided by normal partial differential times y times open z plus cap a times d times x close times d times y equals right hand side sum with, 4 , summands z plus cap a times d times x plus cap b times d times y plus open cap a sub y plus cap a sub z times cap b close times d times x times d times y

Similarly, if we first move along y and then along x ,  we arrive at the following point:

sum with, 4 , summands z plus cap a times d times x plus cap b times d times y plus open cap b sub x plus cap b sub z times cap a close times d times x times d times y

Now we require that whatever route is chosen it leads to the same point on the manifold (up to the terms of the second order). This leads to the following equation:

cap a sub y plus cap a sub z times cap b minus cap b sub x minus cap b sub z times cap a equals zero

negative normal partial differential divided by normal partial differential times y times open cap p divided by cap r close plus normal partial differential divided by normal partial differential times z times open cap p divided by cap r close times cap q divided by cap r minus cap q divided by cap r minus normal partial differential divided by normal partial differential times z times open cap q divided by cap r close times cap p divided by cap r equals zero

equation left hand side sum with, 3 , summands cap p times open normal partial differential times cap q divided by normal partial differential times z minus normal partial differential times cap r divided by normal partial differential times z close plus cap q times open normal partial differential times cap r divided by normal partial differential times x minus normal partial differential times cap p divided by normal partial differential times z close plus cap r times open normal partial differential times cap p divided by normal partial differential times y minus normal partial differential times cap q divided by normal partial differential times x close equals right hand side zero times open one close

Now that was the book, here are some thoughts about this theory.

1) First, equation open asterisk operator close   definitely points to some categories of vector analysis. Indeed, the factors of P, Q and R are the components of the rotor of the vector field cap v right arrow of cap p comma cap q comma cap r . Hence, the condition can be rewritten in a more compact form:

cap v right arrow times r times o times t times cap v right arrow equals zero

At first sight this should hold trivially for any cap v right arrow , for the rotor is by definition perpendicular to the plane defined by cap v right arrow and tangent normal partial differential times cap v right arrow . However, this would only be true, if the solution were indeed an 2-dimensional manifold. If there is no such solution, then the whole derivation becomes invalid.

2) There is another reason why I prefer the geometric argument over comparing the mixed derivatives. The logic is very similar to that used to derive Cauchy-Riemann conditions for the analytic function. Remarkably enough, we can also apply complex formalism to the above problem. Consider the following operator:

equation left hand side diamond operator equals right hand side normal partial differential divided by normal partial differential times x plus i times normal partial differential divided by normal partial differential times y

where equation left hand side i squared equals right hand side negative one .

Assuming again that the solution exists in the form z equals phi of x comma y and using the above shortcuts for partial derivatives we obtain:

equation left hand side prefix diamond operator of z equals right hand side cap a plus i times cap b

Now apply diamond operator super asterisk operator to both parts, where asterisk operator is complex conjugate:

equation sequence diamond operator super asterisk operator diamond operator z equals open normal partial differential divided by normal partial differential times x minus i times normal partial differential divided by normal partial differential times y close times open cap a plus i times cap b close equals open normal partial differential divided by normal partial differential times x times open cap a close plus normal partial differential divided by normal partial differential times y times open cap b close close minus i times open cap a sub y plus cap a sub z times cap b minus cap b sub x minus cap b sub z times cap a close

normal partial differential divided by normal partial differential times x times open cap a close stands for "full partial derivative" where dependance of z on x is taken into account. Replacing cap a and cap b with their values, we obtain:

equation left hand side diamond operator super asterisk operator diamond operator z equals right hand side normal cap delta times z minus i times open cap a sub y plus cap a sub z times cap b minus cap b sub x minus cap b sub z times cap a close

where normal cap delta is the Laplacian, On the other hand:

equation sequence diamond operator super asterisk operator postfix diamond operator equals open normal partial differential divided by normal partial differential times x plus i times normal partial differential divided by normal partial differential times y close times open normal partial differential divided by normal partial differential times x minus i times normal partial differential divided by normal partial differential times y close equals normal cap delta

Hence

equation sequence negative normal black letter cap i times diamond operator super asterisk operator diamond operator z equals cap a sub y plus cap a sub z times cap b minus cap b sub x minus cap b sub z times cap a equals zero

which gives the above integrability condition.

So this is another example of how recourse to complex values can reveal deep facts behind the otherwise unfamiliar looking expressions. And formulate them in a nice compact form as well.

In conclusion here is an example where integrability condition does not hold:

sum with, 3 , summands d times x plus open y plus z close times d times y plus z times d times z

To solve it we rewrite it as follows:

sum with, 4 , summands d times x plus y times d times y plus z times d times z plus z times d times y equals zero

d times open x plus y squared plus z squared divided by two close plus z times d times y equals zero

Now let

x plus y squared plus z squared divided by two equals phi of y

where phi is an arbitrary function. Then

z equals negative phi super prime times open y close

These 2 relations give the general solution.

Permalink Add your comment
Share post
Valentin Fadeev

Gunter's example

Visible to anyone in the world
Edited by Valentin Fadeev, Tuesday, 6 Jul 2010, 23:04

This one gave me some hard time:

z times open x plus z close times normal partial differential times z divided by normal partial differential times x minus y times open y plus z close times normal partial differential times z divided by normal partial differential times y equals zero

where z equals Square root of y when x equals one

The difference with "ordinary" linear equations is that coefficients here depend on z . At the same time right hand side is identically zero, so this is not strictly an inhomogeneous equations.

I made false starts trying to find some nice substitute to absorb z or alternatively pull out the missing term to construct an inhomogeneous equations. Finally, I found a hint in the book N.M. Günter, Integration of First-Order Partial Differential Equations, ONTI/GTTI, Leningrad/Moscow (1934). In the solution for an equation of a similar structure it was suggested to use the standard method and treat z as a constant when integrating the associated system.

So here we go. Searching for a solution in an implicit form:

cap v of z comma x comma y equals zero

z times open x plus z close times normal partial differential times cap v divided by normal partial differential times x minus y times open y plus z close times normal partial differential times cap v divided by normal partial differential times y equals zero

The associated system:

equation sequence d times x divided by z times open x plus z close equals negative d times y divided by y plus z equals d times z divided by zero equals d times t

In the same book it is hinted that one particular integral of this system is z equals cap c sub one which follows from the last identity. Another integral can be found using the first identity:

d times x divided by z times open x plus z close plus d times y divided by y plus z equals zero

Treating z as constant we can simplify the expressions:

d of x divided by z divided by one plus x divided by z plus d of y divided by z divided by y divided by z times open one plus y divided by z close equals zero

equation sequence x divided by z equals u times y divided by z equals v

d times u divided by one plus u plus d times v divided by v minus d times v divided by one plus v equals zero

equation left hand side u times open v plus one close divided by v plus one equals right hand side cap c sub two

equation left hand side open x plus z close times y divided by y plus z equals right hand side cap c sub two

Therefore the general solution can be written in the form:

equation sequence cap v of z comma x comma y equals normal cap phi of z comma open x plus z close times y divided by y plus z equals zero

It is already within reach of sheer guess to let normal cap phi of u comma v equals u minus v to establish the result, however we proceed with a more lengthy, yet rigorous way.

The system of the first integrals is written as follows:

psi sub one equals z

equation left hand side psi sub two equals right hand side open x plus z close times y divided by y plus z

Using the initial condition x equals one :

psi one macron equals z

equation left hand side psi two macron equals right hand side open one plus z close times y divided by y plus z

Solving for z and y :

z equals psi one macron

equation sequence y equals z times psi two macron divided by one plus z minus psi two macron equals psi one macron times psi two macron divided by one plus psi one macron minus psi two macron

Now following the standard method already mentioned below:

MathJax failure: TeX parse error: Extra open brace or missing close brace

equation left hand side psi sub one equals right hand side psi sub two divided by one plus psi sub one minus psi sub two

z equals open x plus z close times y divided by y plus z divided by one plus z minus open x plus z close times y divided by y plus z

equation left hand side sum with, 3 , summands z times y plus z squared plus open z squared minus x times y close times z equals right hand side x times y plus z times y

open z squared minus x times y close times open z plus one close equals zero

Suppressing the solution z equals negative one which does not satisfy intial conditions, finally we obtain:

z equals Square root of x times y

 

Permalink
Share post
Valentin Fadeev

Going further

Visible to anyone in the world
Edited by Valentin Fadeev, Monday, 21 Jun 2010, 19:59

Strolling struggling on:

equation left hand side open m times z minus n times y close times normal partial differential times z divided by normal partial differential times x plus open n times x minus l times z close times normal partial differential times z divided by normal partial differential times y equals right hand side l times y minus m times x

This is an inhomogeneous equation. Following the theory we try to
find the general solution in an implicit form:

cap v of z comma x comma y equals zero

It is proven that the solution found in this form is indeed
general, i.e. we are not losing any solutions on the way.

equation left hand side normal partial differential times z divided by normal partial differential times x equals right hand side negative normal partial differential times cap v divided by normal partial differential times x divided by normal partial differential times cap v divided by normal partial differential times z

equation left hand side normal partial differential times z divided by normal partial differential times y equals right hand side negative normal partial differential times cap v divided by normal partial differential times y divided by normal partial differential times cap v divided by normal partial differential times z

sum with, 3 , summands open m times z minus n times y close times normal partial differential times cap v divided by normal partial differential times x plus open n times x minus l times z close times normal partial differential times cap v divided by normal partial differential times y plus open l times y minus m times x close times normal partial differential times cap v divided by normal partial differential times z equals zero

Now we can write the associated system (in the symmetrical form):

equation sequence d times x divided by m times z minus n times y equals d times y divided by n times x minus l times z equals d times z divided by l times y minus m times x equals d times t

Or more conveniently for this example, in the canonical form:

equation left hand side d times x divided by d times t equals right hand side m times z minus n times y

equation left hand side d times y divided by d times t equals right hand side n times x minus l times z

equation left hand side d times z divided by d times t equals right hand side l times y minus m times x

Multiplying these equations by l, m and n respectively and summing we get:

d times open sum with, 3 , summands l times x plus m times y plus n times z close divided by d times t equals zero

equation left hand side sum with, 3 , summands l times x plus m times y plus n times z equals right hand side cap c sub one

This is one of the first integrals of the system. Now multiplying the equations by x, y and z respectively amd summing we obtain:

one divided by two times d times open sum with, 3 , summands x squared plus y squared plus z squared close divided by d times t equals zero

equation left hand side sum with, 3 , summands x squared plus y squared plus z squared equals right hand side cap c sub two

Therefore, the general solution has the following form:

normal cap phi of sum with, 3 , summands l times x plus m times y plus n times z comma sum with, 3 , summands x squared plus y squared plus z squared equals zero

Geometrically the first solution represents a plane in a 3d space with angular coefficients of the normal vector cosine of alpha colon cosine of beta colon cosine of gamma equals l colon m colon n . The second integral represents a sphere centered at the origin. Therefore, the characteristics of the equation (the curves,

resulting from intersection of these surfaces) are the circles centered on the line passing through the origin with the above mentioned angular coefficients.

Indeed, another way to look at it is rewrite the equation in the following form:

$$\left|\begin{array}{ccc}

\frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} & \frac{\partial u}{\partial z}\\

l & m & n\\

x & y & z\end{array}\right|=0$$

Or:

tau right arrow dot operator open n right arrow multiplication r right arrow close equals zero

Where tau right arrow is a tangential vector to the surface u of x comma y comma z equals zero , n right arrow is the vector of the axis of revolution and r right arrow is the radius vector of an arbitrary point on the surface. It means that for every point on the surface the tangent vector must lie in the plane

passing through the axis of revolution. This is natural, for the surface is ontained by rotating a plane curve against the axis.

Permalink
Share post
Valentin Fadeev

Going partial

Visible to anyone in the world
Edited by Valentin Fadeev, Friday, 18 Oct 2019, 07:52

Struggling through the early chapters of partial differential equations (from the all-time-classic Stepanov's book) I often run across the moves that are not really technically demanding, but the logic behind takes time to sink in.

One of them is the trick used to find a particular solution of the first order linear PDE satisfying initial condition:

f of x sub one comma x sub two comma full stop full stop full stop comma x sub n minus one comma x sub n super zero equals phi of x sub one comma x sub two comma full stop full stop full stop comma x sub n minus one

Here is an excercise: find the solution of the following equation

sum with, 3 , summands Square root of x times normal partial differential times f divided by normal partial differential times x plus Square root of y times normal partial differential times f divided by normal partial differential times x plus Square root of z times normal partial differential times f divided by normal partial differential times z equals zero

         satisfying the condition: f equals y minus z when x equals one

We start by writing out the associated system of ODEs:

equation sequence d times x divided by Square root of x equals d times y divided by Square root of y equals d times z divided by Square root of z

Use the first identity to find a "first integral" of the system:

equation left hand side d times x divided by Square root of x equals right hand side d times y divided by Square root of y

equation sequence psi sub one equals Square root of x minus Square root of y equals cap c sub one

Then equate the first and the third quotients to obtain the second "first integral":

equation left hand side d times x divided by Square root of x equals right hand side d times z divided by Square root of z

equation sequence psi sub two equals Square root of x minus Square root of z equals cap c sub two

Obviously there are no more independent first integrals

According to the theory, the general solution is given as an arbitrary function of these two integrals:

equation sequence f equals normal cap phi of psi sub one comma psi sub two equals normal cap phi of Square root of x minus Square root of y comma Square root of x minus Square root of z

Now we shall find the particular solution. Let x=1. Following the book, we introduce new functions psi one macron and psi two macron to which psi sub one and psi sub two turn when we set the value of x :

equation left hand side psi one macron equals right hand side one minus Square root of y

equation left hand side psi two macron equals right hand side one minus Square root of z

Now we solve these equations with respect to y and z :

y equals open one minus psi macron sub one close squared times open one close

z equals open one minus psi two macron close squared

And now comes the difficult part. In order to get the final result we need to substitute the above results in the expression of f replacing psi macron with psi :

equation sequence f equals y minus z equals open one minus psi sub one close squared minus open one minus psi sub two close squared equals psi sub one squared minus psi sub two squared plus two times open psi sub two minus psi sub one close

equation sequence f equals x minus two times Square root of x times y plus y minus x plus two times Square root of x times z minus z plus two times open Square root of y minus Square root of z close equals y minus z plus two times open Square root of x plus one close times open Square root of z minus Square root of y close

It took me some time to accept that I have to substitute psi which is a more general expression than  psi macron in the equation which was derived after giving x ist particular value.

Finally:

f equals y minus z plus two times open Square root of x plus one close times open Square root of z minus Square root of y close

The reasoning is clear, once it is explained.  psi sub one and psi sub two are both solutions, hence the function of them is also a solution. Then, if we let x equals one they turn to  psi one macron and psi two macron . Then by virtue of the equations (1) we get the variables y and z and consequently the expression of phi of y comma z as required.

In the book the whole argument is presented in the most general form. However, it takes several excercises and hours of thinking to get a hands-on experience with this method.

Permalink
Share post

This blog might contain posts that are only visible to logged-in users, or where only logged-in users can comment. If you have an account on the system, please log in for full access.

Total visits to this blog: 54107