Monday, October 7, 2019

fluid dynamics - Physical intuition for Incompressible Navier Stokes second term


I'm trying to get an intuitive understanding of the Incompressible Navier Stokes equation (as for me thats the only way I can use it effectively and avoid the rote learning method):


$\rho(\frac{\partial \bf u}{\partial t} + \sum_i u_i\frac{\partial \bf u}{\partial x_i}) = - \nabla P + \nu \nabla^2 \bf u$



So far I get the idea that its an applied Newtons Second Law, $F = ma$.


With (volume specific) mass times acceleration on the LHS and the forces from the pressure gradient field, and fluid friction forces on the RHS. (I've omitted the other force term for simplicity).


The term that I'm struggling to fit into this picture is the '$\rho\sum_i u_i\frac{\partial \bf u}{\partial x_i}$' term.


So far my best attempt to explain it is that each directional momentum component is multiplied with the change (pd) of the velocity field relative to that direction. This gives an idea of the potential change in the direction of momentum with change in positional coordinate, which in turn describes the turbulence of the system. - Personally I feel this is a poor explanation that doesn't quite ring true and I'm don't feel like I have understood it properly. Could someone offer a better explanation of the physical intuition of this term in simple language?



Answer



The term combines with the $\partial \mathbf{u} / \partial t$ term to form a convective derivative, which has a simple physical interpretation.


Suppose you're standing by a highway and a red car going 80 mph passes by. Half a second later a blue car going 10 mph passes by. Your reaction wouldn't be "wow, that guy in the red car must have really slammed the brakes, since 10 is much less than 80!", because these are not the speeds of the same object. The convective derivative fixes this mistake.


Mathematically, if you want to write Newton's second law for a particle, the acceleration is simple $d \mathbf{u} / d t$ where $\mathbf{u}(t)$ is the particle's velocity at time $t$. For a fluid, $\mathbf{u}$ is now a function $\mathbf{u}(\mathbf{r}, t)$ which describes the local velocity of the fluid at position $\mathbf{r}$ at time $t$. But $\partial \mathbf{u} / \partial t$ doesn't describe the acceleration of any physical object. If you compare $\mathbf{u}(\mathbf{r}, t)$ and $\mathbf{u}(\mathbf{r}, t + \Delta t)$ you are comparing the velocities of two different pieces of the fluid, not computing the acceleration of one piece.


To fix this, let's instead consider the velocity of a single, fixed piece of fluid, $$\mathbf{v}(t) = \mathbf{u}(\mathbf{r}(t), t)$$ where $\mathbf{r}(t)$ is the position of that piece of fluid at time $t$. The acceleration of that piece is $$\mathbf{a}_{\text{piece}} = \frac{d \mathbf{v}}{dt} = \frac{\partial \mathbf{u}}{\partial \mathbf{r}} \cdot \frac{d \mathbf{r}}{dt} + \frac{\partial \mathbf{u}}{\partial t}$$ by the chain rule. But by definition $d\mathbf{r}/dt = \mathbf{u}$, so $$\mathbf{a}_{\text{piece}} = \frac{\partial \mathbf{u}}{\partial r_i} u_i + \frac{\partial \mathbf{u}}{\partial t}$$ which is exactly what appears in the equation.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...