4
$\begingroup$

This comes up in the context of a homework assignment. We're given the Euler equations for invisicid fluid flow. The variables at play:

  • $p=p(x,y,z,t)$ is pressure
  • $\rho = \rho(x,y,z,t)$ is mass density
  • $\vec v = \vec v(x,y,z,t)$ is velocity of the fluid
  • $\vec f$ is external force per volume

Then we're given that the Euler equations are:

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec v) = 0 \qquad \frac{\partial \rho \vec v}{\partial t} + \vec v \cdot \nabla (\rho \vec v) = \vec f - \nabla p$$

for continuity and momentum, respectively.


My concerns lie with the order of operations dot products and the $\nabla$ operator, and in turn the consistency of the dimensions of these equations. (Which might sound a bit silly, but my physics knowledge isn't the best.)

More explicitly, consider the expression $\vec x \cdot \nabla \vec y$ as in the question title. Which order of operations would be correct: finding $\nabla \vec y$ and then taking the dot product, or taking the adjunction $\vec x \cdot \nabla$ and multiplying that by $\vec y$? Or, symbolically, which does the title mean:

$$\vec x \cdot \Big( \nabla \vec y \Big) \qquad \text{or} \qquad \Big( \vec x \cdot \nabla \Big) \vec y$$

The reason this is a concern lies with the momentum equation given earlier. $\partial_t (\rho \vec v)$ should be a vector, and the right-hand side is the difference of vectors. However, my intuition suggests that of the two above conventions, we should be using the left one, which would render $\vec v \cdot \nabla (\rho \vec v)$ a scalar, which doesn't make sense.

So I was wondering if it happens to instead be $(\vec v \cdot \nabla) (\rho \vec v)$, which, to my understanding, would indeed be a vector. Or is there something else I'm missing?

$\endgroup$
3
  • $\begingroup$ If you aren’t comfortable with tensors, take the dot product first and then apply the resulting “directional derivative” operator. $\endgroup$
    – G. Smith
    Commented Aug 30, 2020 at 21:29
  • 4
    $\begingroup$ If you work out the product $\nabla \vec{y}$ you will see that it is in fact a 2D array and not just a vector. The order doesn't matter in the final expression, it will all produce a vector in the end. $\endgroup$
    – Triatticus
    Commented Aug 30, 2020 at 22:23
  • 2
    $\begingroup$ Very related: physics.stackexchange.com/questions/334509/… $\endgroup$
    – user258881
    Commented Aug 31, 2020 at 5:06

1 Answer 1

9
$\begingroup$

They're all the same.

$$\vec a \cdot \nabla \vec b \equiv (\vec a \cdot \nabla) \vec b \equiv \vec a \cdot (\nabla \vec b)$$

In cartesian coordinates, the $\nabla$ operator is defined as:

$$ \nabla = \left( \begin{array}{ccc}\frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \end{array}\right) $$

When this operates on a vector $\vec b = \left(\begin{array}{ccc}b_x & b_y & b_z\end{array}\right)^T$, we get a tensor of the form:

$$ \nabla \vec b = \left[ \begin{array}{ccc} \frac{\partial b_x}{\partial x} & \frac{\partial b_x}{\partial y} & \frac{\partial b_x}{\partial z} \\ \frac{\partial b_y}{\partial x} & \frac{\partial b_y}{\partial y} & \frac{\partial b_y}{\partial z} \\ \frac{\partial b_z}{\partial x} & \frac{\partial b_z}{\partial y} & \frac{\partial b_z}{\partial z} \\ \end{array} \right] $$

And when we take the dot product of this tensor with the vector $\vec a$, the result is a vector:

\begin{align} \vec a \cdot ( \nabla \vec b ) &= \left( \begin{array}{c} a_x \frac{\partial b_x}{\partial x} + a_y \frac{\partial b_x}{\partial y} + a_z \frac{\partial b_x}{\partial z} \\ a_x \frac{\partial b_y}{\partial x} + a_y \frac{\partial b_y}{\partial y} + a_z \frac{\partial b_y}{\partial z} \\ a_x \frac{\partial b_z}{\partial x} + a_y \frac{\partial b_z}{\partial y} + a_z \frac{\partial b_z}{\partial z} \\ \end{array} \right) \tag 1 \end{align}

Now because $\vec a$ and $\nabla$ are both vectors, $\vec a \cdot \nabla$ is a scalar operation given by $$ \vec a \cdot \nabla = a_x \frac{\partial}{\partial x} + a_y \frac{\partial}{\partial y} + a_z \frac{\partial}{\partial z} $$

And when this scalar operation is applied to the vector $\vec b$, we get a vector:

\begin{align} (\vec a \cdot \nabla) \vec b &= \left( a_x \frac{\partial}{\partial x} + a_y \frac{\partial}{\partial y} + a_z \frac{\partial}{\partial z} \right) \cdot \left( \begin{array}{c} b_x \\ b_y \\ b_z \end{array} \right) \\ &= \left( \begin{array}{c} a_x \frac{\partial b_x}{\partial x} + a_y \frac{\partial b_x}{\partial y} + a_z \frac{\partial b_x}{\partial z} \\ a_x \frac{\partial b_y}{\partial x} + a_y \frac{\partial b_y}{\partial y} + a_z \frac{\partial b_y}{\partial z} \\ a_x \frac{\partial b_z}{\partial x} + a_y \frac{\partial b_z}{\partial y} + a_z \frac{\partial b_z}{\partial z} \\ \end{array} \right) \tag 2 \end{align}

Hey, look at that! $(1) = (2)$

$\endgroup$
2
  • $\begingroup$ Ah, I see, thanks for the help - I guess I was mixing things up with some other notation or something I found along the way. Coincidentally, I noticed that the matrix for the tensor $\vec b$ looks quite like a Jacobian, too. That in turn led me to the Wikipedia page where, it turns out, that notation is sometimes used in that manner as well! So thanks further for helping to connect even those two dots for me! $\endgroup$ Commented Sep 1, 2020 at 21:24
  • $\begingroup$ Yes! The del operator applied to a vector field is the Jacobian (or its transpose, can never remember which) of the vector field, because it gives the gradient of the vector field in all directions. $\endgroup$
    – pho
    Commented Sep 1, 2020 at 21:43

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.