Dear CFD experts:

Can someone suggestion a reference site where I can find the formula on how to compute the following:

grad(H^2) where H is a vector field.

Thanks!

phsieh2005

If I am not mistaken, it should be

grad(H^2) = 2 * H dot grad(H)

This typically shows up for ferromagnetic force, or dielectrophoretic forces.

Opaque

Dear Mr. Opaque:

You are GOOD!

Yes, I am trying to compute the force on paramagnetic particles in the flow field exposed to permanent magnets.

Thanks a lot!

phsieh2005

Sorry, I don't believe this is correct:

1/2 GRAD(h^2) = hj hj,i ei = [GRAD h]^T [h]

this derives by the following:

hj,i ei = [GRAD h] ^ T (transpose)

= hj,i (ei Tensor ej)

Applying the tensor product:

hj,i (ei Tensor ej) * hk ek = delta_jk hk hj,i ei

collapsing jk using delta function:

hj,i (ei Tensor ej) * hk ek = hj hj,i ei

In other words: GRAD(h^2) = 2 hj hj,i ei (repeating both i and j)

or

1/2 GRAD(h^2) = h1(h1,1 e1 + h1,2 e2 + h1,3 e3) + ...

h2(h2,1 e1 + h2,2 e2 + h2,3 e3) + ...

h3(h3,1 e1 + h3,2 e2 + h3,3 e3)