0

My question is similar to this one about convective derivative. But how do I adapt it to the cylindrical coordinate? Simply changing Grad to cylindrical coordinate does not produce correct results, which can be found here, because we need to include the terms with the derivative of the unit vector too.

Here is the adapted code,

uVec = Through[{Subscript[u, r], Subscript[u, \[Theta]], 
     Subscript[u, z]}[r, \[Theta], z, t]];
dConvect[vec_] := 
 D[vec, t] + ((vec.Grad[#, {r, \[Theta], z}, "Cylindrical"]) & /@ 
    vec)
dConvect[uVec] // MatrixForm // TraditionalForm

which produces the result below.

enter image description here

iplasma
  • 1
  • 1

0 Answers0