ストークスの式とは？沈降速度の導出と計算で陥る罠を徹底解説【2026】

化学工学、粉体プロセス、リチウムイオン電池の電極スラリー調製、さらには土木分野の地盤調査に至るまで、流体中を沈降する微粒子の挙動把握は歩留まりと品質を決定づける中核技術です。その基礎理論として19世紀から活用され続けているのが「ストークスの式」です。しかし、理論のシンプルさに反して、研究室や製造現場では「計算値と実測値が著しく乖離する」「単位換算を誤って致命的な設計ミスに直面する」といったトラブルが後を絶ちません。

微粒子の挙動を正確に予測するためには、単に公式を丸暗記するだけでなく、流体力学的な背景や厳密な前提条件、そして適用可能な限界領域を正しく見極める眼が欠かせません。本稿では、第一線の取材データと実務現場の検証知見をもとに、ストークスの式の完全な導出ステップから、エンジニアが直面する計算の落とし穴、レイノルズ数に基づく適用限界までを体系的に解き明かします。

📌 【この記事の重要ポイントまとめ】

• 要点1：ストークスの式は粘性抵抗力・浮力・重力のつり合いから導かれ、微粒子の終末沈降速度が「粒子径の2乗」に比例することを明確に示す基盤方程式です。
• 要点2：計算ミスの大半は「粘度・密度の単位換算不備」と「レイノルズ数の上限突破（非層流領域への無理な適用）」という初歩的かつ構造的な見落としに起因します。
• 要点3：厳密な適用限界（Re＜0.1〜1）や非球形粒子・濃厚スラリーによる干渉沈降の限界を理解し、補正係数や数値解析モデルへ切り替える判断基準を持つことが実務成功の鍵となります。

微粒子沈降の金字塔「ストークスの式」とは？流体力学の基礎と成立背景

微視的な粒子が静止流体中を沈降していく様子を記述するストークスの式（ストークスの法則）は、1851年にアイルランドの物理学者・数学者であるジョージ・ガブリエル・ストークスによって提唱されました。現代の流体力学の基礎において、極低速かつ粘性が卓越する領域を解析する際の原点として位置づけられています。

流体の運動を包括的に表す支配方程式といえばナビエ・ストークス方程式ですが、非線形項（慣性項）を含むため、厳密解を数学的に導き出すことは極めて困難です。そこでストークスは、「粒子の沈降速度が極めて遅く、流体の慣性力が粘性力に比べて無視できるほど微小である」という極限状態を想定しました。この大胆な物理的仮定がストークス近似（クリーピング流れ・遅い流れの近似）と呼ばれ、ナビエ・ストークス方程式の慣性項をゼロとみなすことで、美しい解析解の導出を可能にしたのです。

この理論の有用性を世界に知らしめた歴史的マイルストーンが、1909年にロバート・ミリカンが行った「ミリカンの油滴実験」です。ミリカンは微小な帯電油滴を電場と重力場で浮遊・落下させ、ストークスの法則を援用して油滴の半径を逆算することで、電気素量（電子1個の電荷）の高精度な測定に成功しました。現代においても、最先端マテリアル開発や大気エアロゾル分析のバックボーンとして脈々と息づいています。

【公式の導出】力のつり合いから終末沈降速度の計算式を導く全手順

流体中を自然落下する微粒子は、落下開始直後は重力加速度によって加速しますが、速度の上昇に伴って流体からの抵抗力が増大し、最終的に加速度がゼロの等速運動へ移行します。この一定速度が終末沈降速度（終端速度）です。単一の球形粒子沈降をモデルとし、微粒子に働く3つの力のつり合いから沈降速度計算式を導出します。

粒子にかかる力は以下の3つに分解されます。

1. 重力（Fg）：下向きに作用する粒子の質量による力。
粒子の直径を $D_p$、粒子密度を $\rho_p$、重力加速度を $g$ とすると、
$F_g = \frac{\pi}{6} D_p^3 \rho_p g$

2. 浮力（Fb）：アルキメデスの原理に基づき、粒子が排除した流体の重さに等しい上向きの力。
流体密度を $\rho_f$ とすると、
$F_b = \frac{\pi}{6} D_p^3 \rho_f g$

3. 粘性抵抗力（Fd）：ストークス近似に基づき、流体の粘性によって粒子運動を妨げる上向きのドラッグフォース。
流体の粘度（流体粘度）を $\mu$、沈降速度を $v$ とすると、半径 $r$（$= D_p / 2$）を用いて $6\pi \mu r v$、すなわち直径表示では、
$F_d = 3\pi \mu D_p v$

終末速度に達した瞬間、下向きの重力と上向きの抗力（浮力＋粘性抵抗力）がつり合います。
$F_g = F_b + F_d$

各大項目を代入して式を整理します。
$\frac{\pi}{6} D_p^3 \rho_p g = \frac{\pi}{6} D_p^3 \rho_f g + 3\pi \mu D_p v$

$3\pi \mu D_p v = \frac{\pi}{6} D_p^3 (\rho_p - \rho_f) g$

両辺を $3\pi \mu D_p$ で除算することで、最終的なストークスの沈降速度式が完成します。
$v_t = \frac{g (\rho_p - \rho_f) D_p^2}{18 \mu}$

この関係式から導き出される最も重要な物理的直感は、「終末沈降速度は粒子径の2乗に比例する」という点です。粒子サイズが2倍になれば沈降速度は4倍になり、10倍になれば沈降速度は100倍に跳ね上がります。逆に、ナノスケールや数マイクロメートルの微粒子になると、沈降速度は極端に低下し、静止水面上にほぼ永久に漂い続けるような挙動を示します。

実務・試験で計算がつまずく決定的な理由｜単位換算と前提条件の落とし穴

製造現場の技術者や工学系受講者がストークスの計算を行う際、結果のオーダーが1000倍や100万倍ずれてしまう事故が頻発します。取材現場でエンジニアから寄せられた声をもとに、計算が破綻する3大要因を浮き彫りにします。

最大の落とし穴は、単位系の不整合（SI単位系への統一漏れ）です。現場の実務データでは、粒径はマイクロメートル（$\mu\text{m}$）、密度は $\text{g/cm}^3$、流体の粘度はセンチポアズ（$\text{cP} = \text{mPa}\cdot\text{s}$）で表記されるケースが一般的です。これらをそのまま公式に代入すると、算出される速度単位は全く意味をなしません。

安全に計算を行うためには、必ずすべての変数をSI基本単位（$\text{m}$, $\text{kg}$, $\text{s}$）に換算してから計算を実行する習慣づけが不可欠です。

• 粒径：$10\,\mu\text{m} \rightarrow 10 \times 10^{-6}\,\text{m}$
• 密度：$2.5\,\text{g/cm}^3 \rightarrow 2500\,\text{kg/m}^3$、水（20℃）：$1.0\,\text{g/cm}^3 \rightarrow 1000\,\text{kg/m}^3$
• 粘度：$1\,\text{cP} \rightarrow 1 \times 10^{-3}\,\text{Pa}\cdot\text{s}\,(\text{kg}/(\text{m}\cdot\text{s}))$
• 重力加速度：$9.8\,\text{m/s}^2$

もうひとつの見落としが「密度差（$\rho_p - \rho_f$）」の符号です。水処理やエマルションの分離工程では、対象微粒子が水より軽いケース（気泡や油滴など）が存在します。この場合、密度差が負の値となり、終末速度はマイナス、すなわち「沈降」ではなく「浮上速度」を表すことになります。現象の向きを前提条件として整理しておかないと、符号の取り違えによる装置設計ミスに直結します。

レイノルズ数の適用限界と流動領域の境界線【データ比較検証】

ストークスの式は無条件で適用できる万能ツールではありません。理論の前提である「慣性力無視・完全な層流」を維持するためには、粒子周りの流れを表す粒子レイノルズ数（$Re_p$）を厳密に照合する必要があります。

粒子レイノルズ数は次の式で定義されます。
$Re_p = \frac{\rho_f v D_p}{\mu}$

学術的にはレイノルズ数の適用限界として $Re_p < 0.1$、工業的な許容範囲を含めても $Re_p \le 1.0$ を超えると、粒子の後流に渦（剥離流）が生じ始め、流体抵抗は粘性摩擦だけでなく形状抗力（圧力抵抗）の影響を強く受けます。その結果、抵抗係数は $C_D = 24/Re$ の比例関係から逸脱し、ストークスの式による沈降予測は急激に過大評価へと転じます。

水中の沈降を例にとると、常温水（粘度約 $1.0 \times 10^{-3}\,\text{Pa}\cdot\text{s}$、密度 $1000\,\text{kg/m}^3$）において密度 $2500\,\text{kg/m}^3$ の砂粒子を落下させた場合、粒径がおよそ $80\,\mu\text{m}$（$0.08\,\text{mm}$）を超えた段階でレイノルズ数は1.0を上回ります。つまり、肉眼で辛うじて粒と認識できるようなサイズ以上の粒子では、もはやストークスの式単体での評価は成立しません。

【実態検証】電池製造から土質力学まで｜現場目線で見えたストークスの式のリアル

ストークスの式が産業界で最もダイナミックに活用されている現場のひとつが、リチウムイオン二次電池の極板製造ラインです。正極活物質や導電助剤（カーボンブラック）、バインダーを溶媒中に分散させたスラリーにおいて、粒子が保管タンク内で予期せぬ沈降を起こすと、塗布重量のバラつきや充放電容量の劣化に直結します。現場では「製造待機時間中に生じる沈降限界」を評価するため、ストークス式をベースに分散剤による見かけ粘度の調整設計が行われています。

また、建設・インフラ分野の基盤である地質調査現場でも、ストークスの式は規格試験の中枢を担っています。日本産業規格に規定される土質力学沈降試験（JIS A 1204「土の粒度試験」における沈降分析・比重計法）では、シルトや粘土といった $75\,\mu\text{m}$ 未満の細粒分を対象に、懸濁液の比重経時変化からストークスの式を用いて粒径測定と粒度分布曲線を算出しています。

しかし、地盤工学コンサルタントや材料工学の試験現場からは、次のようなリアルな課題が指摘されています。「比重計法試験において、水温が1℃変化するだけで水の粘度は約2〜3％変動する。温度管理を怠ると粒径データに無視できない歪みが生じる」「粘土鉱物の平板状粒子を球形と仮定して算出するため、得られた径はあくまで『ストークス相当径』であることを忘れて実地判断を下すと、地盤の圧密沈下予測を誤る」といった現場特有の苦い経験談は少なくありません。

一般に知られていない盲点とネットの誤解｜「万能の沈降式」という幻想を暴く

Web上の技術解説や入門テキストでは単純化されがちなストークスの式ですが、現実の物理系へのアプローチには3つの致命的な死角が存在します。

盲点1：濃厚分散系における「干渉沈降」の無視
ストークスの式は「無限遠まで壁が存在せず、周囲に他の粒子が存在しない単一粒子の沈降」を前提条件としています。しかし、体積分率が数％を超える懸濁液では、ある粒子が沈降する際に周囲の流体を押し上げ、その上昇流が隣接粒子の落下を強烈に押し戻す「干渉沈降（Hindered Settling）」が発生します。体積分率が20〜30％に達するスラリーでは、実際の沈降速度はストークス計算値の数分の一から数十分の一にまで激減します。

盲点2：ナノ領域で立ちはだかる「ブラウン運動」の壁
「粒子が小さければ小さいほどストークス近似の精度は高まる」という認識は半分正しく、半分誤りです。粒径がサブミクロン（約 $0.1\,\mu\text{m}$ 以下）に達すると、流体分子の不規則な熱運動による衝突エネルギー（ブラウン運動）が重力沈降のエネルギーを凌駕します。結果として粒子は沈降せず分散状態を保ち続けるため、沈降速度式による設計自体が無効化されます。

盲点3：容器の「壁面効果（Wall Effect）」による抗力増大
細いメスシリンダーや流路内で沈降試験を行う場合、粒径に対して容器の内径が十分に大きくない（管径比が100倍未満）と、管壁に起因するせん断応力が粒子の落下を抑制します。実務データではフランシスの補正係数などを用いて壁面効果を差し引かなければ、真の材料特性は見えてきません。

【プロの結論】現場導入とモデル選択における判断基準

工学的な課題解決において、ストークスの式を有効活用できるか否かは、「理論の境界線」に対する冷静な自己チェックができるかどうかにかかっています。

ストークスの式を直ちに適用できるケース

• 対象粒径が明確に $0.5\,\mu\text{m} \sim 50\,\mu\text{m}$ の範囲に収まっている
• 懸濁液の固形分濃度が極めて希薄（体積分率1%未満）である
• 流体がニュートン流体であり、局所的なせん断速度による粘度変化がない
• 粒子形状がほぼ等方的な球形、または「球相当径」としての指標化で実用上問題がない

適用を避け、補正モデルや数値解析（CFD）へ切り替えるべきケース

• レイノルズ数 $Re_p > 1.0$ の領域（粗砂、大型の沈降物、気泡など）
• カーボンナノチューブや雲母など、アスペクト比が極端に大きい非球形粒子（形状補正係数や楕円体沈降モデルが必要）
• 固形分濃度が高い濃縮スラリー（Richardson-Zakiの干渉沈降補正式を採用すべき領域）
• 流体がチキソトロピー性や擬塑性を示す非ニュートン流体（塗料、ポリマー溶液、高濃度電極ペーストなど）

数式は単なる計算の道具ではなく、「特定の境界条件のもとで成立する物理モデル」に過ぎません。適用限界を客観的に評価し、前提条件から外れた際には迷わず高次モデルへ移行する柔軟性こそが、実務における工学的トラブルを未然に防ぐ防壁となります。

【ストークスの式】に関するよくある質問（FAQ）

Q1：沈降速度の計算値が実測値より大幅に速いのですが、何が原因ですか？
A1：最も頻度の高い原因は「粒子濃度の高さによる干渉沈降」または「非球形粒子による抵抗増大」です。粒子同士の距離が近い場合、排除された流体の上昇流が抵抗となり、実測速度はストークス計算値の50％以下に落ち込むことがあります。また、粒子が板状や針状の場合も表面積が大きくなり、球体モデルより大幅に遅くなります。

Q2：空気中の微粒子沈降（大気エアロゾル）にもそのまま使えますか？
A2：空気の密度と粘度を代入すれば基本的には適用可能ですが、粒径が $1\,\mu\text{m}$ を下回ると、空気分子の平均自由行程と同等サイズになり、流体が連続体として振る舞わなくなる「すべり（Cunningham補正）」が生じます。サブミクロン粒子の場合は、カニンガム補正係数を導入して抵抗力を低減補正する必要があります。

Q3：レイノルズ数が1を超えているか事前に分からない場合はどう計算すればよいですか？
A3：まず仮定としてストークスの式で沈降速度 $v_t$ を算出し、その結果得られた速度を用いてレイノルズ数 $Re_p = \rho_f v_t D_p / \mu$ を事後検算してください。もし $Re_p \le 1.0$ であればその解を採用し、$Re_p > 1.0$ に達した場合は、アレン域の実験式や抗力係数曲線を用いた逐次近似計算（反復法・ゴールシーク等）へと切り替えるのが定石です。

まとめ：今後の動向と失敗しないための判断基準

ストークスの式は、誕生から170年以上が経過した現在もなお、微粒子プロセスの骨格を支える強力な基礎理論です。リチウムイオン電池の高性能化や微細泥土の環境処理など、最先端の工学課題を解決するアプローチも、すべてはこのシンプルな力のつり合いから始まります。

成否を分けるのは、「公式を単に知っていること」ではなく、「公式が破綻する限界領域を冷徹に把握していること」です。単位系の厳密な統制、レイノルズ数による流れの判別、そして粒子濃度や形状に伴う物理現象の乖離を先回りして検知すること。この工学的リテラシーを備えて臨むことこそが、設計の失敗を回避し、シミュレーションと実務を強固に結びつける最短距離となります。 （出典: ストークス の 式（Yahoo!ニュース））