Register to reply

Curl of gradient

by TrickyDicky
Tags: curl, gradient
Share this thread:
TrickyDicky
#1
Apr5-12, 06:16 AM
P: 3,012
Under what circumstances can a conservative field NOT be irrotational?
Phys.Org News Partner Science news on Phys.org
Bees able to spot which flowers offer best rewards before landing
Classic Lewis Carroll character inspires new ecological model
When cooperation counts: Researchers find sperm benefit from grouping together in mice
NasuSama
#2
Sep17-12, 03:09 PM
P: 326
Quote Quote by TrickyDicky View Post
Under what circumstances can a conservative field NOT be irrotational?
According to Wikipedia:

"Conservative vector fields are also irrotational, meaning that (in three-dimensions) they have vanishing curl. In fact, an irrotational vector field is necessarily conservative provided that a certain condition on the geometry of the domain holds: it must be simply connected.
An irrotational vector field which is also solenoidal is called a Laplacian vector field because it is the gradient of a solution of Laplace's equation."

We say that the field is irrotational if its curl is zero.

See more: http://en.wikipedia.org/wiki/Conservative_vector_field

Good luck!
Simon Bridge
#3
Oct13-12, 08:38 PM
Homework
Sci Advisor
HW Helper
Thanks
Simon Bridge's Avatar
P: 12,475
eg.
there are no conditions in which a conservative field is irrotational.
Did you have a particular situation in mind?

A Dhingra
#4
Oct18-12, 08:30 AM
P: 196
Curl of gradient

TrickyDicky,
Do you have any situations in your mind for existence of the curl of a gradient or conservative field being rotational? Honestly I am looking forward to hear that..
(and that's why after waiting for a while I wrote this.. )
Zelyucha
#5
Feb21-13, 05:27 PM
P: 25
Quote Quote by TrickyDicky View Post
Under what circumstances can a conservative field NOT be irrotational?
When the vector field is on a surface that is not simply connected. How is this so?
Have a look at the curl operator in matrix form:



If you look closely you can see that the curl operator, which I will refer to from now on as is a skew-symmetric matrix. It turns out that the exponential map of a skew-symmetric matrix is an orthogonal matrix(its transpose is equal to its inverse). This can be shown by defining the matrix S = e and ST=e-∇. Now we take the matrix product SST = ee-∇ = e(∇)-(∇) = e[0] = I3 = e-∇+∇ = e-∇e = STS = I3 where [0] is the zero matrix(a 33 matrix of all zero entries) and I3 is the square identity matrix of dimension 3.

Now that we have shown that the exponential of is an orthogonal matrix, there's an even more detailed proof here showing that the exponential of a skew symmetric matrix like the Curl operator is a rotation matrix.

So if G is the gradient of a 3 dimensional scalar field and ∇G = 0, then e∇G = I3 which is the identity element of SO(3)( the group of rotations about the origin of vectors in Euclidean 3 space R3) which is itself not simply connected. So the map of all simply connected scalar fields in R3 maps to the Kernel{SO(3)} as is not 1-to-1. Whereas scalar fields in
R3 that are not simply connected map into SO(3)~I3. This map is injective but not surjective. So the curl of the gradient of a scalar field that is not simply connected corresponds to a rotation operator.
micromass
#6
Feb21-13, 06:04 PM
Mentor
micromass's Avatar
P: 18,062
Quote Quote by Zelyucha View Post
When the vector field is on a surface that is not simply connected. How is this so?
Have a look at the curl operator in matrix form:



If you look closely you can see that the curl operator, which I will refer to from now on as is a skew-symmetric matrix. It turns out that the exponential map of a skew-symmetric matrix is an orthogonal matrix(its transpose is equal to its inverse). This can be shown by defining the matrix S = e and ST=e-∇. Now we take the matrix product SST = ee-∇ = e(∇)-(∇) = e[0] = I3 = e-∇+∇ = e-∇e = STS = I3 where [0] is the zero matrix(a 33 matrix of all zero entries) and I3 is the square identity matrix of dimension 3.

Now that we have shown that the exponential of is an orthogonal matrix, there's an even more detailed proof here showing that the exponential of a skew symmetric matrix like the Curl operator is a rotation matrix.

So if G is the gradient of a 3 dimensional scalar field and ∇G = 0, then e∇G = I3 which is the identity element of SO(3)( the group of rotations about the origin of vectors in Euclidean 3 space R3) which is itself not simply connected. So the map of all simply connected scalar fields in R3 maps to the Kernel{SO(3)} as is not 1-to-1. Whereas scalar fields in
R3 that are not simply connected map into SO(3)~I3. This map is injective but not surjective. So the curl of the gradient of a scalar field that is not simply connected corresponds to a rotation operator.

You're mixing things up. The curl of a gradient is always zero. No matter what the domain of the vector field is. The "converse" is not generally true. That is: if the curl of a vector field is 0, then it is the gradient of something. That requires the domain to be simply connected. This last statement is essentially the Poincar lemma.
WannabeNewton
#7
Feb21-13, 06:04 PM
C. Spirit
Sci Advisor
Thanks
WannabeNewton's Avatar
P: 5,435
It is much more appropriate to talk about conservative co - vector fields NOT vector fields because conservative has line integral definitions for which you need co - vector fields but if you care only about euclidean 3 - space then it's not really a big deal. A conservative vector field corresponds to an exact one form and an irrotational vector field corresponds to a closed one - form in the context of euclidean 3 - space. Every exact form is closed. On the other hand, a closed form is exact on a contractible domain. All these things are MUCH MUCH more elegant and natural in the more general context of differential forms.
mathwonk
#8
Feb21-13, 09:11 PM
Sci Advisor
HW Helper
mathwonk's Avatar
P: 9,453
there might be a continuity assumption. otherwise, always.
Zelyucha
#9
Feb21-13, 10:21 PM
P: 25
Quote Quote by micromass View Post
You're mixing things up. The curl of a gradient is always zero. No matter what the domain of the vector field is. The "converse" is not generally true. That is: if the curl of a vector field is 0, then it is the gradient of something. That requires the domain to be simply connected. This last statement is essentially the Poincar lemma.

UGH.....Ya got me! iStand corrected.

But I'd like to see the proof of this if you will.
WannabeNewton
#10
Feb21-13, 10:35 PM
C. Spirit
Sci Advisor
Thanks
WannabeNewton's Avatar
P: 5,435
Quote Quote by Zelyucha View Post
UGH.....Ya got me! iStand corrected.

But I'd like to see the proof of this if you will.
[itex]\omega = df[/itex] hence [itex]d\omega = d(df) = d^{2}f = 0[/itex]. See how pretty it is with differential forms =D. Or you know, if you're a masochist, go ahead and work out the curl of a gradient in coordinates xD.
Zelyucha
#11
Feb23-13, 12:57 PM
P: 25
micromass, upon looking at some old notes of mine attempting to disprove that ∇∇f = 0
for any scalar function f mapping R2 into R3 which is continuous and is differentiable up to at least order 2, I found a page which actually proves it to be correct.

One can write δ2f/δyδx = lim(Δy→0) {(f'(x,y+Δy,z) -f'(x,y,z))/Δy} and the quantity in the brackets can be re-written as lim(Δx→0){ (f(x+Δx,y,z) -f(x,y,z))/Δx}. Substituting we get:

δ2f/δyδx = lim(ΔyΔx→0) { (f(x+Δx,y+Δy,z) - f(x+Δx,y,z) + f(x,y,z))/ΔyΔx}.

and

δ2f/δxδy = lim(ΔxΔy→0) { (f(x+Δx,y+Δy,z) - f(x,y+Δy,z) + f(x,y,z))/ΔxΔy}.

Since multiplication of real numbers is commutative, ΔyΔx = ΔxΔy and if we subtract the 2 quantities in brackets we now have:

2f/δyδx -δ2f/δxδy] = lim(ΔxΔy→0){ (f(x,y+Δy,z)-f(x+Δx,y,z))/ΔxΔy}

By the continuity of f: R2→R3, there is a point a = (x,y) in R2 whose image is f(x,y,z) and any open δ-neighborhood has an ε-neighborhood on G centered at Img(a) in G. So if we consider the points ax =(x+Δx,y) and ay=(x,y+Δy) in R2 and the vector X connecting them, then |X| = √(Δx2+Δy2). Now I claim that
lim(ΔxΔy→0){ (f(x,y+Δy,z)-f(x+Δx,y,z))/ΔxΔy} = 0. Which implies that for any ε > 0, there exists a δ > 0 such that:

|(f(x,y+Δy,z)-f(x+Δx,y,z))|/|ΔxΔy| < ε whenever the distance d(ax,ay) < δ. So what is the delta? Well since the vectors from a→ax and a→ay are orthogonal(to see this just set (x,y)=0), if you were to move one of those points towards the other along the X which connects them such that d(ax,ay) < √(Δx2+Δy2) = |X|, then we can choose δ = |X| and place our δ-disk at the midpoint of X between ax and y. So shrinking ε means we shrink the magnitude of Δy and Δx to move the 2 a-points even closer together along X to re-size our δ-disk. And as ΔxΔy→0(along with |X|), then the |Img(X)|→0 and the proof is complete. [δxiδxj]δf = 0 whenever i≠j.
micromass
#12
Feb23-13, 01:14 PM
Mentor
micromass's Avatar
P: 18,062
Quote Quote by Zelyucha View Post
micromass, upon looking at some old notes of mine attempting to disprove that ∇∇f = 0
for any scalar function f mapping R2 into R3 which is continuous and is differentiable up to at least order 2, I found a page which actually proves it to be correct.

One can write δ2f/δyδx = lim(Δy→0) {(f'(x,y+Δy,z) -f'(x,y,z))/Δy} and the quantity in the brackets can be re-written as lim(Δx→0){ (f(x+Δx,y,z) -f(x,y,z))/Δx}. Substituting we get:

δ2f/δyδx = lim(ΔyΔx→0) { (f(x+Δx,y+Δy,z) - f(x+Δx,y,z) + f(x,y,z))/ΔyΔx}.

and

δ2f/δyδx = lim(ΔxΔy→0) { (f(x+Δx,y+Δy,z) - f(x,y+Δy,z) + f(x,y,z))/ΔxΔy}.

Since multiplication of real numbers is commutative, ΔyΔx = ΔxΔy and if we subtract the 2 quantities in brackets we now have:

2f/δyδx -δ2f/δyδx] = lim(ΔxΔy→0){ (f(x,y+Δy,z)-f(x+Δx,y,z))/ΔxΔy}

By the continuity of f: R2→R3, there is a point a = (x,y) in R2 whose image is f(x,y,z) and any open δ-neighborhood has an ε-neighborhood on G centered at Img(a) in G. So if we consider the points ax =(x+Δx,y) and ay=(x,y+Δy) in R2 and the vector X connecting them, then |X| = √(Δx2+Δy2). Now I claim that
lim(ΔxΔy→0){ (f(x,y+Δy,z)-f(x+Δx,y,z))/ΔxΔy} = 0. Which implies that for any ε > 0, there exists a δ > 0 such that:

|(f(x,y+Δy,z)-f(x+Δx,y,z))|/|ΔxΔy| < ε whenever the distance d(ax,ay) < δ. So what is the delta? Well since the vectors from a→ax and a→ay are orthogonal(to see this just set (x,y)=0), if you were to move one of those points towards the other along the X which connects them such that d(ax,ay) < √(Δx2+Δy2) = |X|, then we can choose δ = |X| and place our δ-disk at the midpoint of X between ax and y. So shrinking ε means we move the 2 a points even closer together along X to re-size our δ-disk. And as ΔxΔy→0(along with |X|), then the |Img(X)|→0 and the proof is complete. [δxiδxj]δf = 0 whenever i≠j.
Your proof looks more like a proof for the theorem of Clairaut-Schwarz... http://en.wikipedia.org/wiki/Symmetr...nd_derivatives
Zelyucha
#13
Feb23-13, 05:33 PM
P: 25
Quote Quote by micromass View Post
Your proof looks more like a proof for the theorem of Clairaut-Schwarz... http://en.wikipedia.org/wiki/Symmetr...nd_derivatives
Well I'm not surprised that someone didn't construct a theorem about mixed partial derivatives(and prove it too). My proof was a specific case of n=3. And I can see that the Clairaut-Schwartz theorem is an induction of this concept to n dimensions. The commutativity of mixed partials for a continuous function f of differential class C2 that defines a (2-D) surface in Euclidean 3 space is a consequence of continuity as well as the Cauchy-Schwartz inequality used for defining ε-δ neighborhoods.
Zelyucha
#14
Feb25-13, 12:04 AM
P: 25
Might I add, the fact that the Curl of the Gradient of a scalar function is always zero is a direct consequence of Clairaut's theorem since the theorem states that mixed partials commute.

But what about a continuous well defined function of differential class C2+ on real projective space? Or more specifically a scalar function f in R3 that defines a non-orientable surface?
micromass
#15
Feb25-13, 03:49 AM
Mentor
micromass's Avatar
P: 18,062
Quote Quote by Zelyucha View Post
Might I add, the fact that the Curl of the Gradient of a scalar function is always zero is a direct consequence of Clairaut's theorem since the theorem states that mixed partials commute.

But what about a continuous well defined function of differential class C2+ on real projective space? Or more specifically a scalar function f in R3 that defines a non-orientable surface?
Of course you can calculate partial derivatives on weird manifolds. Then the partial derivatives don't necessarily commute since you can have a weird chart. You can define a Lie-bracket to see how much the partial derivatives commute.

But curl and gradient are defined completely different on manifolds. They are not longer vector fields but differential forms. The curl and gradient correspond to the exterior derivative and you can show that applying the exterior derivative twice gives 0.
Zelyucha
#16
Feb26-13, 01:11 PM
P: 25
Quote Quote by micromass View Post
Of course you can calculate partial derivatives on weird manifolds. Then the partial derivatives don't necessarily commute since you can have a weird chart. You can define a Lie-bracket to see how much the partial derivatives commute.

But curl and gradient are defined completely different on manifolds. They are not longer vector fields but differential forms. The curl and gradient correspond to the exterior derivative and you can show that applying the exterior derivative twice gives 0.
A conservative vector field for this purpose is simply the gradient of a function that defines a 2-dimensional manifold embedded in R3. So when you speak of manifolds you are speaking of a completely separate topological space whereas what I'm talking about is restricted to (2D)surfaces in Euclidean 3D space. The CURL operator is a differential 2-form on such an object.

Now consider the scalar function ψ: (x,y,z) → (x/√(x2+y2+z2),y/√(x2+y2+z2),z/√(x2+y2+z2)) = (ψ(x),ψ(z),ψ(z))

This is a mapping of R3 into Real projective space(of dimension 2). If you take ∇ψ = (δψ(x)/δx, δψ(y)/δy,δψ(z)/δz) and then apply the CURL operator, you will find that ∇∇ψ ≠ 0 ! So any scalar field f(x,y,z) that is defined by a continuous mapping of the Real Projective Plane will have asymmetric mixed partial derivatives and a non-vanishing CURL. And while I haven't proved it just yet; I'm going to conjecture that:

Let f is a continuous scalar function of 3 real variables and differential class C2+ that defines a surface in R3. If f is orientable then ∇f is a conservative vector field.
micromass
#17
Feb26-13, 06:25 PM
Mentor
micromass's Avatar
P: 18,062
I am very sorry to say that I could make very little sense of your post.

Quote Quote by Zelyucha View Post
A conservative vector field for this purpose is simply the gradient of a function that defines a 2-dimensional manifold embedded in R3.
I'm not following. If you define your manifold as a level set of [itex]F:\mathbb{R}^3\rightarrow \mathbb{R}[/itex], then the gradient is a normal field on F. I don't get why you call this a conservative vector field at all.

So when you speak of manifolds you are speaking of a completely separate topological space whereas what I'm talking about is restricted to (2D)surfaces in Euclidean 3D space. The CURL operator is a differential 2-form on such an object.
The curl operator is not a 2-form.

Now consider the scalar function ψ: (x,y,z) → (x/√(x2+y2+z2),y/√(x2+y2+z2),z/√(x2+y2+z2)) = (ψ(x),ψ(z),ψ(z))
You say that

[tex]\psi(x)=\frac{x}{\sqrt{x^2 + y^2 + z^2}}[/tex]

But the right hand side depends clearly on x, y and z. And the left-hand side depends on x. So the notation is not very legal.

This is a mapping of R3 into Real projective space(of dimension 2). If you take ∇ψ = (δψ(x)/δx, δψ(y)/δy,δψ(z)/δz)
In my knowledge, a gradient is only defined for a scalar field, that is for a sufficiently smooth function [itex]f:\mathbb{R}^3\rightarrow \mathbb{R}[/itex]. You seem to be taking a gradient for something whose codomain is not [itex]\mathbb{R}[/itex]. I don't know how that is defined.
lavinia
#18
Feb26-13, 07:20 PM
Sci Advisor
P: 1,716
In the planar case one can think of the curl of a vector field as the divergence of its rotation by -90 degrees.

∇xV = ∇.A(V) where a is the -90 degree rotation of V(x) around x.

For if V = (u,w) then A(v) = (w,-u) so ∇.A(V) = w[itex]_{x}[/itex] - u[itex]_{y}[/itex]

The divergence gives the infinitesimal flux density of the vector field away or towards the point. It is a simple limiting arguement to show this.

So one can think of the curl of the vector field as the infinitestimal rotation density either clockwise or counter clockwise around the point.


Register to reply

Related Discussions
Zero curl and gradient of some scalar potential Calculus 1
Gradient,divergence and curl. Calculus 36
Purpose of each of the operators , divergence, gradient and curl? Calculus 12
Use of curl of gradient of scalar Differential Equations 10
Vecctor analysis and got the mathematical formulae for gradient General Physics 4