Operation that will undo a curl operation?

[Brad_Ad23]

Does anyone know if there is an operation that will undo a curl operation?

 

[matt grime]

No, curl isn't bijective.

 

[Brad_Ad23]

That isn't good news. Hmm, if it isn't bijective, what is it? I'm going to figure this would have to do with the fact that the matrix operator for curl is a singular matrix?

 

[matt grime]

Why ought it to be bijective? A simple explanation from a decent differential geometer would be handy right now, but I'm not one of those.

It is not bijective as there are plenty of things with curl equal to zero (irrotational fields)

As the matrix operator is not a linear map on R^3, and is just a matrix of partial derivatives, I'm not sure how you'd even begin to define singularity.

 

[Organic]

http://mathworld.wolfram.com/Curl.html

 

[matt grime]

Originally posted by Organic
http://mathworld.wolfram.com/Curl.html

And the point of that was...?

 

[Hurkyl]

The other problem is that the curl of a vector field is always divergence free, so unless your field is divergence free, you cannot find an anticurl.


Nonuniqueness of anticurl isn't a problem; just like with ordinary antiderivatives you return a class of solutions. In this case, you'd have a base solution plus an arbitrary irrotational field. Sure, it's a lot messier than just adding an arbitrary constant, but it can be done.


Anyways, there are almost certainly better ways to do it, but you can find an anti curl just by looking at the partial differential equations. The procedure is pretty much in the same spirit as when you're looking for the antigradient of an exact vector field.

(NOTE: I'm assuming that everything involved is well-behaved)

Suppose that you're given this vector field:

F = <0, y e^z - z \cos (xz), 2x - e^z>

(exercise: show that this is divergence free)

You want to find A = <f, g, h> such that \nabla \times A = F.

First, note that f can be anything, because our answer is only unique up to adding an arbitrary irrotational field. I'm going to hunt for the most generic answer, but if you only need one answer, you could simplify life by setting f = 0.

From the definition of curl, we have:

<br /> \begin{equation*}\begin{split}<br /> h_y - g_z &amp;= 0 \\<br /> f_z - h_x &amp;= y e^z - z \cos (xz) \\<br /> g_x - f_y &amp;= 2x - e^z<br /> \end{split}\end{equation*}<br />

So now solve.

<br /> \begin{equation*}\begin{split}<br /> g_x &amp;= f_y + 2x - e^z \\<br /> g &amp;= x^2 - x e^z + \int f_y \, dx \\<br /> h_x &amp;= f_z - y e^z + z \cos (xz) \\<br /> h &amp;= \sin (xz) - y e^z + \int f_z \, dx<br /> \end{split}\end{equation*}<br />

Now, remember that the "constants" of integration here are arbitrary functions of both y and z. You now use the third equation, h_y - g_z &amp;= 0 to fix what those "constants" may be.

One particular solution is A = &lt;0, x^2 - x e^z, \sin (xz) - x y e^z&gt;, so the general solution can be found by adding an arbitrary irrotational field to this.


(incidentally, the solution I used to generate the initial vector field was A = &lt;y e^z, x^2, \sin (xz)&gt;)

 

[matt grime]

Strictly speaking it is just the non-injectivity that matters: undo is a vague term. For a strict inverse one would require bijectivity, however the word undo could be interpreted to mean F undoes G if

FG(x) =x,

without requiring F to be defined on the whole of the set in which G maps to, merely defined on the image of G.

Anyway, curl(f) =0 for all constant fields (and a whole lot more besides) so even that relaxation won't do.

edit, this was written as Hurkyl's post appeared. the phrase 'so it's only injectivity that counts' is not a contradiction of his post, merely my presonal opinion about interpreting the word 'undo', which means inverse in my mind.

 

[Brad_Ad23]

Hmm, very intriguing. In the case I am looking at, the vector field is indeed divergence free. And I think in this case the irrotational fields will probably wind up canceling out, or not having much of a contribution at any rate, but that does indeed help a bit.

My next question is, suppose you had a vector equation of the form

u_t = curlu. If I were to find the general solution to the case where \nabla \times A = u, would then the property of u_t \times A = u hold?

Note: the vector u is such that each component of u is defined by the form u_1(x,y,z,t) and so on for u_2 and u_3.

 

[Hurkyl]

FYI, the (generally very large) field

<br /> F = &lt;e^{x+y+z}, e^{x+y+z}, e^{x+y+z}&gt;<br />

is irrotational, as is the field

<br /> F = &lt;10^{10^{10}}, 10^{10^{10}}, -10^{10^{10}} &gt;<br />

You can't assume if you get one solution that it will be close to the desired solution.

 

[matt grime]

ignore the t thing

you want to know what

(\nabla\wedge u) \wedge A is when

\nabla\wedge A = u

how are you with summation convention?

(\nabla\wedge u) \wedge A)_r = \epsilon_{riq}\epsilon_{ijk}\partial_ju_kA_q = (\delta_{qj}\delta_{rk}-\delta_{qk}\delta_{rj})\partial_ju_kA_q

which simplifes to

\partial_qu_rA_q - \partial_ru_qA_q

so you want to know what

u(\nabla.A) - \nabla(u.A)

is if it is also true that \nabla \wedge A = u

 

[matt grime]

Ok, it's coming back to me now.

u.A = (\nabla\wedge A).A=0 cos its the product of symmetric and anti symmetric things,

hence we only have to examine \nabla . A I see no reason why this should be 1.

 

[Hurkyl]

My next question is, suppose you had a vector equation of the form...

Have you tried any examples?

 

[Brad_Ad23]

Touche'. I'll try working on it some more now that classes are done for the day.

 

[Kalimaa23]

One thing of interest might be Helmholtz's theorem:

If the divergence and the rotation of a vector field, which varies asymptotically as at least r^(-2) and is sufficiently smooth, is given, then there exists a unique vector field (barring an additive constant) with that particular curl and divergence.

You can actually write it down, but I'm I don't have time to type it up in TeX right now. Watch this space!

 

[Hurkyl]

Don't forget that it requires pesky things like vanishing at infinity.

http://mathworld.wolfram.com/HelmholtzsTheorem.html

and the answer still isn't unique, since there exist nonconstant vector fields with zero divergence and zero curl.

 

[Kalimaa23]

Originally posted by Hurkyl
Don't forget that it requires pesky things like vanishing at infinity.

That's what I meant with the correct asymptotic behaviour. And a really important pesky thing it is, if you can write down a gazzilion forms for a constant B-field that doens't vanish at infinity.