流体の運動を式で扱うと、最終的には「ベクトル演算子 $\nabla$ の文章」になります。
まずは勾配・発散・回転・ラプラシアンを“成分で書ける”状態にして、導出の土台を作る。
結論
この4つ＋移流演算子を“手で展開できる”ようにする。
この記事のゴールは、$\nabla$ を見た瞬間に止まらないことです。
流体の導出で最低限必要なのは、勾配 $\nabla f$／発散 $\nabla\cdot\mathbf{u}$／回転 $\nabla\times\mathbf{u}$／ラプラシアン $\nabla^2\mathbf{u}$ の4点。
どれも「成分に落とす」と同じ型になります。
さらに、導出でつまずきやすいのが $(\mathbf{u}\cdot\nabla)\mathbf{u}$（移流に対応する項）です。
これは“内積っぽい見た目”に反して、ベクトル場に対する微分演算なので、成分で書けるようになるだけで急に安心します。
そしてもう一つ、粘性や応力の話に進む直前で必ず出てくるのが「2階テンソル（行列）」です。
ここは記号として飲み込むのではなく、物理的な見え方（どの方向に動くと、どの方向の速度がどう変わるか）の表だと思うと、急に怖くなくなります。
あるある
式は読めるのに、手が動かない。
教科書で $\nabla\cdot$ や $\nabla\times$ が出ても、「まあそういうもの」と読み飛ばしてしまう。
でも導出で急に「成分で示せ」と来た瞬間に固まるやつです。ここで一度、“読む”から“展開する”に切り替えます。
シママ、流体の式ば追おうとして、どこで止まりよる？
だいたい $\nabla$ が出たところよ…。記号としては知ってるのに、成分に落とすと途端に自信がなくなるの。
よか。今日は「$\nabla$ を怖がらん日」にするばい。まずは4つだけ。勾配・発散・回転・ラプラシアン。これだけ手で開けたら勝ちたい。
4つだけって言うけど、毎回その“4つ”で転ぶのよ…。でも、今日は逃げない。ちゃんと成分まで落として。
本文
$\nabla$ は「部品」だと思えば、式が読みやすくなる。
0) まず約束：座標と記法を固定する
座標は $(x,y,z)$、ベクトル場は $\mathbf{u}=(u_x,u_y,u_z)$ とします。添字を使う場合は $u_1=u_x,\,u_2=u_y,\,u_3=u_z$ のように同一視します。
そして演算子 $\nabla$ は次の“偏微分の束”です。
$$\nabla=\left(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}\right)$$
ここで大事なのは、$\nabla$ 自体は数値ベクトルじゃなくて、作用する相手がいて初めて意味が決まることです。
1) 勾配 $\nabla f$：スカラー場の「増え方の矢印」
$f(x,y,z)$ がスカラー場（温度、圧力など）だとすると、勾配は次です。
$$\nabla f=\left(\frac{\partial f}{\partial x},\frac{\partial f}{\partial y},\frac{\partial f}{\partial z}\right)$$
方向は「最も増える向き」、大きさは「その増え方の強さ」。後で圧力のようなスカラーが出てきたとき、$\nabla p$ が“向きのある量”になる感覚が重要になります。
例題：$f(x,y,z)=x^2+yz$ の $\nabla f$ を求めよ。
解
$$\frac{\partial f}{\partial x}=2x,\ \frac{\partial f}{\partial y}=z,\ \frac{\partial f}{\partial z}=y$$
$$\therefore\ \nabla f=(2x,\ z,\ y)$$
2) 発散 $\nabla\cdot\mathbf{u}$：その場で「湧く/吸う」度合い
$\mathbf{u}=(u_x,u_y,u_z)$ の発散は次です。
$$\nabla\cdot\mathbf{u}=\frac{\partial u_x}{\partial x}+\frac{\partial u_y}{\partial y}+\frac{\partial u_z}{\partial z}$$
流体だと、これは局所的な体積膨張率の指標になります。非圧縮性の議論では次が“条件”として使われ、導出中に消える項が出てきます。
$$\nabla\cdot\mathbf{u}=0$$
例題：$\mathbf{u}=(x^2,\ yz,\ \sin z)$ の $\nabla\cdot\mathbf{u}$ を求めよ。
解
$$\frac{\partial u_x}{\partial x}=2x,\ \frac{\partial u_y}{\partial y}=z,\ \frac{\partial u_z}{\partial z}=\cos z$$
$$\therefore\ \nabla\cdot\mathbf{u}=2x+z+\cos z$$
発散が0ってのは、「そこに流体が溜まったり減ったりしない」ってことたい。条件として入れるから、式の中で消す判断ができる。
なるほど。“条件”として入れるから、式の中で消す判断ができるのね。今までは、0って書いてあるだけで置物だった。
3) 回転 $\nabla\times\mathbf{u}$：局所的な「回りやすさ」
回転は偏微分の組み合わせです。
$$\nabla\times\mathbf{u}=\left(\frac{\partial u_z}{\partial y}-\frac{\partial u_y}{\partial z},\ \frac{\partial u_x}{\partial z}-\frac{\partial u_z}{\partial x},\ \frac{\partial u_y}{\partial x}-\frac{\partial u_x}{\partial y}\right)$$
流体では渦度として使われ、速度場がどれくらい回転しているかの指標になります。さらに、粘性の整理や恒等式で $\nabla\times(\nabla\times\mathbf{u})$ の形が顔を出します。
例題：$\mathbf{u}=(y,\ 0,\ 0)$ の $\nabla\times\mathbf{u}$ を求めよ。
解
$$u_x=y,\ u_y=0,\ u_z=0$$
$$\therefore\ \nabla\times\mathbf{u}=\left(0,\ 0,\ \frac{\partial 0}{\partial x}-\frac{\partial y}{\partial y}\right)=(0,\ 0,\ -1)$$
4) ラプラシアン $\nabla^2\mathbf{u}$：ベクトルの「拡散」
ラプラシアンはスカラーに対しては次です。
$$\nabla^2 f=\frac{\partial^2 f}{\partial x^2}+\frac{\partial^2 f}{\partial y^2}+\frac{\partial^2 f}{\partial z^2}$$
ベクトル場に対しては「各成分に同じことをする」：
$$\nabla^2\mathbf{u}=(\nabla^2 u_x,\ \nabla^2 u_y,\ \nabla^2 u_z)$$
物理的には「局所の平均との差」に関係し、速度の平滑化（拡散）の形として現れます。
例題：$u_x=x^2+y^2+z^2$ の $\nabla^2 u_x$ を求めよ。
解
$$\frac{\partial^2 u_x}{\partial x^2}=2,\ \frac{\partial^2 u_x}{\partial y^2}=2,\ \frac{\partial^2 u_x}{\partial z^2}=2$$
$$\therefore\ \nabla^2 u_x=6$$
5) 2階テンソル（行列）の最小理解：$\nabla\mathbf{u}$ を“速度の地図”として見る
ここが今回の強化ポイントです。2階テンソルと言われると急に難しく見えますが、
$\nabla\mathbf{u}$ は要するに 「どの方向に少し動いたとき、速度ベクトルがどの成分でどれだけ変わるか」を並べた表です。
行列って聞くと身構えるけど、これは“速度の地図”たい。
ちょっと右（$x$方向）に動いたら $u_x$ がどう変わる、$u_y$ がどう変わる…ってのを並べとるだけ。
“地図”なら分かる。つまり、場所を少し変えたときの速度ベクトルの変化を、まとめて見てるってことね。
まず定義はこれだけです。
$$(\nabla\mathbf{u})_{ij}=\frac{\partial u_i}{\partial x_j}$$
この式の読み方は単純で、列（$j$）が「どの方向に動いたか」、行（$i$）が「速度のどの成分を見ているか」です。つまり、
読み方
・$j=x$（1列目）：$x$方向に少し動いたとき、速度がどう変わるか
・$i=x$（1行目）：$x$方向の速度成分 $u_x$ がどう変わるか
行列で書けばこうです。
$$\nabla\mathbf{u}=\begin{pmatrix}\frac{\partial u_x}{\partial x} & \frac{\partial u_x}{\partial y} & \frac{\partial u_x}{\partial z}\\\frac{\partial u_y}{\partial x} & \frac{\partial u_y}{\partial y} & \frac{\partial u_y}{\partial z}\\\frac{\partial u_z}{\partial x} & \frac{\partial u_z}{\partial y} & \frac{\partial u_z}{\partial z}\end{pmatrix}$$
物理的には、これは「小さな立方体の流体要素が、次の瞬間にどう形を変えようとしているか」を決める表です。
例えば $\partial u_x/\partial y$ が大きいと、$y$方向に並んだ点で $x$方向速度が違うので、せん断（横ずれ）が起きます。
一方で $\partial u_x/\partial x$ が大きいと、$x$方向に引き伸ばす（圧縮する）ような 伸縮が起きます。
代表成分の物理
・$\partial u_x/\partial x$：$x$方向の伸び縮み（引き伸ばし／圧縮）
・$\partial u_x/\partial y$：$y$方向の並びに対する $x$速度差 → せん断（ずれ）
転置は行と列の入れ替えで、意味としては「見方を逆にする」だけです。
$$(\nabla\mathbf{u})^T_{ij}=\frac{\partial u_j}{\partial x_i}$$
この2つを平均したものが、後で粘性の本体になる「変形速度テンソル」です。
ここでは“式は外に置く”ルールで、式だけ出します。
$$\mathbf{S}=\frac12\left(\nabla\mathbf{u}+(\nabla\mathbf{u})^T\right)$$
なぜ平均するのか（直感）
せん断には「上が速い／下が遅い」と「左が速い／右が遅い」みたいな“相互のずれ”が出ます。
そのずれを“方向に依らず同じ物理量”として扱うために、対称化（平均）する。
反対称成分は「剛体回転」に対応し、粘性はそこには反応しにくい、という整理に繋がります。
例題：$\mathbf{u}=(u_x,u_y,0)$ のとき、$S_{xy}$ を定義から書け。
解
$$S_{xy}=\frac12\left(\frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x}\right)$$
6) $(\mathbf{u}\cdot\nabla)\mathbf{u}$：演算子として読む
ここでの $\mathbf{u}\cdot\nabla$ は“微分演算子”です。
$$\mathbf{u}\cdot\nabla = u_x\frac{\partial}{\partial x}+u_y\frac{\partial}{\partial y}+u_z\frac{\partial}{\partial z}$$
これが $\mathbf{u}$ に作用するので、成分で書くと次です。
$$(\mathbf{u}\cdot\nabla)\mathbf{u}=\left(u_x\frac{\partial u_x}{\partial x}+u_y\frac{\partial u_x}{\partial y}+u_z\frac{\partial u_x}{\partial z},\ u_x\frac{\partial u_y}{\partial x}+u_y\frac{\partial u_y}{\partial y}+u_z\frac{\partial u_y}{\partial z},\ u_x\frac{\partial u_z}{\partial x}+u_y\frac{\partial u_z}{\partial y}+u_z\frac{\partial u_z}{\partial z}\right)$$
添字で短く書くならこれです。
$$\left[(\mathbf{u}\cdot\nabla)\mathbf{u}\right]_i=u_j\frac{\partial u_i}{\partial x_j}$$
例題：$\mathbf{u}=(x,\ y,\ 0)$ のとき、$(\mathbf{u}\cdot\nabla)\mathbf{u}$ を求めよ。
解
$$\mathbf{u}\cdot\nabla=x\partial_x+y\partial_y$$
$$(\mathbf{u}\cdot\nabla)\mathbf{u}=\left(x\frac{\partial x}{\partial x}+y\frac{\partial x}{\partial y},\ x\frac{\partial y}{\partial x}+y\frac{\partial y}{\partial y},\ 0\right)=(x,\ y,\ 0)$$
解釈：この速度場は“自分自身の向きに沿って”増えていくので、移流による加速も同じ形になる。
7) よくあるつまずき（先に潰す）
つまずき1：$\nabla$ を数値ベクトルだと思う
$\nabla$ は演算子。相手（$f$ や $\mathbf{u}$）が来て初めて意味が決まる。
つまずき2：$(\mathbf{u}\cdot\nabla)\mathbf{u}$ を “$\mathbf{u}$ の2乗” だと思う
これは微分演算子がベクトル場に作用している。成分展開で形を固定するのが一番早い。
つまずき3：2階テンソルを「謎の記号」として覚えてしまう
$\nabla\mathbf{u}$ は “速度の地図”。行＝速度成分、列＝移動方向。せん断と伸縮の表として見えると強い。
テンプレ
導出中に止まったら、ここに戻る。
流体の導出で記号が出て止まったときの、最小復帰テンプレです。「この形に落とせるか？」だけ確認します。
（1）$\nabla$ は偏微分の束
$$\nabla=\left(\partial_x,\partial_y,\partial_z\right)$$
（2）発散・回転・ラプラシアンは“型”
$$\nabla\cdot\mathbf{u}=\partial_x u_x+\partial_y u_y+\partial_z u_z,\quad \nabla^2\mathbf{u}=(\nabla^2 u_x,\nabla^2 u_y,\nabla^2 u_z)$$
（3）移流演算子は演算子として読む
$$\mathbf{u}\cdot\nabla=u_x\partial_x+u_y\partial_y+u_z\partial_z,\quad \left[(\mathbf{u}\cdot\nabla)\mathbf{u}\right]_i=u_j\partial_{x_j}u_i$$
（4）テンソルが出たら “速度の地図” に戻す
$$(\nabla\mathbf{u})_{ij}=\partial_{x_j}u_i,\quad \mathbf{S}=\tfrac12(\nabla\mathbf{u}+(\nabla\mathbf{u})^T)$$
落とし穴
「分かった気になる」だけで次へ行く。
落とし穴は一つです。記号を“知ってる”だけで進んでしまうこと。
導出では必ず、どこかで成分展開が必要になります。
4つ（勾配・発散・回転・ラプラシアン）と、移流項の成分形、そして $\nabla\mathbf{u}$ の“行と列の意味”だけは、その場で説明できる状態にしてから進むのが安全です。
締め
次は「保存則の型」で、式が“運動”っぽくなる。
今日の勝ち条件は単純たい。$\nabla$ を見ても固まらんこと。あと、$\nabla\mathbf{u}$ を“速度の地図”って言えることたい。
うん、納得。行と列の意味が分かったら、急にテンソルが怖くなくなった。次はこの部品を使って、保存則の形を組み立てるのね。
そうたい。地図が読めたら、あとは道順たい。迷子にならんけんね。
その言い方、ちょっとだけ好き。でも調子に乗らないで。