Granular Matter / Yield Stress / Jamming

砂は、固体にも
流体にもなる。

砂時計の砂は流れ、握れば形を保ちます。
水でも金属でもないこの中間の振る舞いは、粒状体(granular matter)と呼ばれ、現代物理でもまだ完全には理解されていない題材です。

  • 粒は熱で動かないので、力をかけないと流れない
  • ある「降伏応力」を超えた瞬間、固体から流体に切り替わる
  • 砂山の角度・ジャミング・ブラジルナッツ効果は、すべて粒の接触ネットワークが鍵になる

Particle

粒の集まり

砂の 1 粒は剛体で、摩擦と非弾性衝突を持ちます。普通の分子のように熱で動き回ったりしません。

Threshold

降伏応力

小さな力では固体のままで、ある臨界を超えた瞬間に急に流れ出す。これが流体との大きな違いです。

Network

接触のネットワーク

粒同士の接触がつながっていれば固体、崩れれば流体。砂の正体は、このネットワークの強さです。

Core Idea

砂が流れるかどうかは、粒の接触ネットワークが壊れるかどうかで決まる。

普通の流体は応力に比例して流れますが、砂は「ある力までは動かず、その先で突然崩れる」塑性体的なふるまいをします。

Mechanism

なぜ砂は流体とも固体とも違うのか

砂のふるまいは、水のような流体や金属のような固体だけでは説明できません。粒のサイズ・摩擦・非弾性衝突・熱で動かないという 4 つが重なって、粒状体ならではの振る舞いを生みます。

Step 1

粒は剛体で大きい

1 粒のサイズが分子よりずっと大きいため、変形よりも転がる・ぶつかるが基本になります。

Step 2

摩擦と非弾性

接触点で摩擦が働き、衝突するとエネルギーが熱や音に変わって失われていきます。

Step 3

熱で動かない

分子と違い、室温で勝手に動き出すほどの熱は持っていないため、外から力を入れないと流れません。

Step 4

力を超えると突然流れる

せん断応力がある臨界を超えた瞬間、接触ネットワークが崩れて流体的にふるまいます。

Newtonian Fluid

τ = η · du/dy

水や空気のような普通の流体では、せん断応力は速度勾配に比例します。どれだけ弱い力でも、それに合わせて流れます。

Granular Matter

τ < τy → 動かない
τ ≥ τy → 流れる

砂はビンガム流体や摩擦モデルのように、降伏応力 τy を超えるまで動かず、超えてから流れ始めます。

Models

砂をどう数理化するか

砂を扱う方法は、目的によって粒度が変わります。1 粒ずつ追うのか、流れの場として平均化するのか、あるいは流れやすさを 1 つの関数にまとめるのか、レイヤーで整理できます。

Model 1

粒子シミュレーション (DEM)

1 粒ずつ位置と速度を追い、衝突・摩擦を直接扱う方法です。一番物理に近いですが、粒数が多いと計算量が爆発します。

  • 個々の挙動が見える
  • 摩擦・回転・形状を入れやすい
  • 粒数が増えると重い

Model 2

連続体モデル

砂を密度・速度の場として扱い、Navier-Stokes 方程式を拡張して流れを記述します。マクロな流れの予測に向きます。

  • 大規模な流れに適用可能
  • 降伏応力を粘性の式に組み込む
  • 粒同士の局所的な乱れは消える

Model 3

μ(I) レオロジー

有効摩擦係数 μ を、慣性数 I(粒径・流速・圧力で作られる無次元数)の関数として与える、現代の標準的なモデルです。

  • 多くの実験を統一的に説明
  • μ = μ(I) という単純な形
  • 準静的な極限と高速流の両方を含む

最初の一歩

セルオートマトン的な落下モデル

「重力で落ち、行き先がふさがれていれば斜めに崩れる」という 2 行ルールでも、砂時計や山の挙動の雰囲気を再現できます。

中間レベル

2D の DEM

粒を円として、ばね・ダンパ・摩擦を入れた DEM は、ブラジルナッツ効果やジャミングを再現する定番です。

本格レベル

μ(I) を組み込んだ連続体

ダム崩壊や砂丘の流れなど、大規模な砂の運動を扱うときの研究レベルのモデルです。

Phenomena

砂が起こす代表的な現象

砂の物理を面白くしているのは、いくつかの典型現象です。すべて「接触ネットワークが壊れるかどうか」に戻って理解できます。

Key Word

安息角とジャミング

砂山には自然に立ち上がる限界角度(およそ 30〜35°)があり、これを安息角と呼びます。出口でアーチができて流れが止まるのがジャミングで、どちらも接触ネットワークが粒の運動を支える現象です。

1

安息角

砂山が崩れずに保てる最大の角度。重力と摩擦の釣り合いで決まります。

2

ジャミング

出口でアーチが組まれ、それ以上流れなくなる詰まり現象です。

3

砂の対流

容器を振動させると、内部で対流が起き、大粒が表面に出てきます(ブラジルナッツ効果)。

4

砂丘の形成

風で粒が跳ねる saltation と、表面の不安定性が組み合わさって、波模様の砂丘が育ちます。

Interactive Demo 1

セルオートマトンで砂を落とす

下の実演は、資料にあった「セルオートマトン型」砂シミュレーションです。マウスドラッグで砂を落とすと、真下が空いていれば落ち、ふさがっていれば斜めに崩れる、という 2 つのルールだけで山が育ちます。

砂の総量 0
山の高さ 0%
状態 流れている

マウスドラッグで砂を加えると、重力で落ち、斜めに崩れて山が育ちます。「上から降らせる量」を上げると、砂時計のように継続的に降り続けます。

明るい砂色のセルが粒、暗い箇所が空気です。仕切りを置くと、その上に山ができたり、横から崩れる様子を観察できます。

Interactive Demo 2

ブラジルナッツ効果を見る

容器に大小の粒を混ぜて振ると、大きな粒がだんだん上に出てきます。この対流的な並べ替えが「ブラジルナッツ効果」で、ミューズリーの中で大きなナッツが上に集まる現象としても知られます。

大粒の平均高さ 0%
容器の振動 中程度
経過フレーム 0

赤い大粒は、最初は容器の下の方に置かれています。振動を上げると、小粒が空隙へ落ち込み、大粒は徐々に上へ押し上げられます。

黄色が小粒、赤が大粒です。容器を上下に揺らすと、振動の慣性と空隙落下の組み合わせで、大粒が表面に浮かんできます。

Regimes

流れの強さで変わる 3 つの相

砂は単に「固体/流体」の 2 択ではなく、流れの慣性数 I によって、ゆるやかに性格が変わります。これが μ(I) レオロジーの背景です。

I が小さい

準静的

流れがほぼ止まっている領域で、安息角や山の自然崩壊で支配的になる相です。

I が中くらい

密な流れ

ホッパーやチューブの中を砂が連続的に流れる状態。μ が I に応じて少しずつ増える領域です。

I が大きい

気体的な相

粒同士が接触ネットワークを失い、衝突主体で運動する高速流。砂塵嵐や雪崩の先端などが該当します。

One Sentence

砂のレオロジーは、慣性数 I が「どの相」にいるかで決まる

静止した山も、ゆっくり崩れる砂時計も、激しい砂塵嵐も、同じ μ(I) の枠組みで連続的に説明できます。

Next Step

砂から、他のテーマへつなげる

粒状体の話題は、シミュレーション・統計・パターン形成・カオスといった他のページとも自然につながります。気になる方向に進んでください。