粒子径分布・ロジン-ラムラー分布解析とは
🧑🎓
粒子径分布って、粉体の大きさのバラツキを表すものだと思うんですけど、d50とかd90って具体的に何を表しているんですか？
🎓
ざっくり言うと、全体の何%がその大きさ以下かを示す値だよ。例えばd50は「中央径」で、粒子を大きさ順に並べた時、真ん中（体積で50%）の粒子の大きさだ。シミュレーターのグラフで、累積分布曲線が50%の高さと交わるところの粒径がd50だね。上の「分布タイプ」を切り替えて、d50の値がどう変わるか見てみよう。
🧑🎓
なるほど！ロジン-ラムラー分布の「n」ってスライダーがありますけど、これは何を変えてるんですか？「分布の急峻さ」って書いてありますが。
🎓
その通り、nが大きいほど粒子の大きさが揃っている（分布が狭い）ことを意味するんだ。例えば、噴霧乾燥で作った微粒子はnが大きく、粗い粉砕で作った砂はnが小さいことが多い。右のグラフでnのスライダーを動かすと、累積分布カーブの傾きが急になったり緩やかになったりするのがわかるよ。d50は同じでも、d10とd90の値が大きく変わるはずだ。
🧑🎓
「ザウター平均径」や「比表面積」もリアルタイムで計算されてますね。これって何に使われるんですか？
🎓
実務ではすごく重要だよ。ザウター平均径d32は、燃焼や溶解の反応速度に直結する。例えば燃料噴霧の液滴が小さい（d32が小さい）ほど燃えやすいんだ。比表面積は、触媒や塗料の性能評価に必須だ。「代表径 d'」のスライダーを動かして、d32と比表面積がどう反応するか確かめてみて。粒子が全体的に小さくなると、表面積は劇的に増えることが実感できるはずだ。
物理モデルと主要な数式
ロジン-ラムラー分布 (Rosin-Rammler)：粉砕や噴霧など、機械的処理で生成される粒子群の粒径分布をよく記述する経験式です。累積通過率Q(d)は以下の式で表されます。
$$Q(d) = 1 - \exp\left[-\left(\frac{d}{d'}\right)^n\right]$$
$d$: 粒子径, $Q(d)$: 粒径d以下の粒子の体積割合（累積通過率）, $d'$: 特性径（Q(d)=63.2%となる粒径）, $n$: 分布係数（大きいほど分布が狭く均一）
対数正規分布 (Log-Normal)：自然に生成される微粒子（エアロゾルなど）や結晶成長で見られる分布です。粒径の対数が正規分布に従うと仮定します。
$$Q(d) = \Phi\left[\frac{\ln(d/d_{50})}{\sigma}\right]$$
$d_{50}$: 中央径（メジアン径）, $\sigma$: 幾何標準偏差（分布の幅）, $\Phi$: 標準正規分布の累積分布関数
ザウター平均径 (Sauter Mean Diameter, d32)：粒子群の体積と表面積の比を保つ等価球の直径で、反応工学で特に重要です。
$$d_{32}= \frac{\sum d_i^3}{\sum d_i^2}$$
比表面積：単位質量あたりの総表面積です。粒子を球と仮定し、密度を$\rho$とすると、$S_v = \frac{6}{\rho d_{32}}$ の関係があります。
実世界での応用
粉砕・製粉工程：鉱石の粉砕やセメントの製造では、製品の粒度分布が反応性や強度を決定します。ロジン-ラムラー分布を用いてミルの性能評価や最適な粉砕条件の設定を行います。
噴霧乾燥・塗装：食品や医薬品の噴霧乾燥では、液滴の粒径分布(d32)が乾燥速度と製品の粒子形状に影響します。塗装では、塗料粒子の分布が仕上がり表面の平滑性を左右します。
燃焼工学：ディーゼルエンジンやガスタービンの燃料噴霧では、ザウター平均径(d32)が燃料と空気の混合速度、ひいては燃焼効率と排ガス特性を支配する重要なパラメータとなります。
粉体プロセス：製薬での錠剤圧縮、セラミックスの成形では、原料粉末の粒度分布が充填性や焼結後の密度に直接影響するため、d10, d50, d90を用いた厳密な管理が行われます。
よくある誤解と注意点
まず、「d50が平均径だ」という思い込みは捨てましょう。d50は中央径（メジアン）であり、算術平均径とは異なります。例えば、非常に細かい粉と粗い粒子が混ざったバイモーダル分布では、d50はその「混ざり具合」の真ん中を示すだけで、粒子全体の「平均的な大きさ」の直感とはズレることがあります。実務では、d50だけでなくd10とd90を必ずセットで確認し、分布の幅を把握することが鉄則です。
次に、分布モデルの安易な選択。ロジン-ラムラー分布は粉砕物に、対数正規分布は自然発生エアロゾルに、という大まかな指針はありますが、測定データに無理やりフィッティングするのは危険です。例えば、粉砕工程でも初期の粗砕品はロジン-ラムラーに従わないことが多々あります。ツールでパラメータをいじる前に、まずは実測データのプロット形状を観察し、「なぜこの分布が適しているのか」を工学的に考える癖をつけましょう。
最後に、比表面積計算の前提条件。ツールで算出される比表面積は、すべての粒子が球で、かつ表面が滑らかという理想仮定に基づいています。しかし、実際の触媒粒子は多孔質ですし、フレーク状の顔料は表面積が桁違いに大きい。このツールの値は「完全に球状だった場合の理論値」として扱い、実測値との乖差から粒子形状の複雑さを推察する材料にするのが賢い使い方です。
関連する工学分野
このツールの計算ロジックは、燃焼工学と深く結びついています。ディーゼルエンジンやガスタービンの燃焼室では、燃料の噴霧液滴径分布が燃焼効率と排ガス特性を支配します。ここで活躍するのがザウター平均径d32です。d32が小さいほど蒸発速度が速く、燃焼は完全かつ迅速に進みます。逆に、d32が大きすぎると未燃燃料（スス）の原因に。シミュレーションツールで分布パラメータとd32の関係を体感することは、燃焼設計の第一歩です。
また、製剤工学（ファーマシューティカル・エンジニアリング）でも核心的な役割を果たします。吸入剤の粒子径分布は、肺のどの部位まで薬剤が到達するかを決定します。例えば、d50が5μm以下の粒子は肺胞まで到達できますが、それより大きいと気道に沈着してしまいます。分布係数nが大きい（分布が狭い）製剤は、薬効と副作用のバランスを精密に制御できるため、高度な品質管理が求められます。
さらに粉体塗装や3Dプリンティング（粉末床溶融結合）では、粒子の流動性と充填密度が品質を左右します。これらは粒径分布に強く依存します。単一径（nが極端に大きい）だと充填率が悪く、適度に幅のある（最密充填に近い）分布を設計する必要があります。ツールでd10, d50, d90を調整しながら、最適な分布形状を探る作業は、まさに材料設計そのものと言えるでしょう。
発展的な学習のために
まず次の一歩は、「モーメント平均径」の概念を理解することです。d50（中央径）やd32（ザウター平均径）は、実は一般式「モーメント平均径 $$D[p, q] = \left( \frac{\sum d_i^p}{\sum d_i^q} \right)^{1/(p-q)}$$」の特殊なケースです。d32はp=3, q=2の場合。d10（体積基準の10%径）とは別物の「数平均径」はp=1, q=0の場合です。この式をマスターすれば、あらゆる平均径の定義が統一的に理解でき、文献で異なる平均径が出てきても混乱しなくなります。
数学的背景を深めたいなら、確率密度関数（PDF）と累積分布関数（CDF）の関係をグラフと数式で行き来できるようになりましょう。このツールが描いている累積分布曲線Q(d)を微分すると、確率密度関数q(d)が得られます。ロジン-ラムラー分布のPDFは $$q(d) = \frac{n}{d'} \left( \frac{d}{d'} \right)^{n-1} \exp\left[-\left(\frac{d}{d'}\right)^n\right]$$ となります。PDFのグラフのピーク（最頻値）とd50の位置がどうずれるかを考えると、分布の非対称性に対する感覚が磨かれます。
実務に直結する次のトピックとしては、粒度測定法の原理と限界を学ぶことを強くお勧めします。レーザー回折法、動的光散乱法、画像解析法…それぞれが出力する「径」の定義は異なります（例えば回折法は体積等価径）。異なる装置で測ったデータを単純比較するのは危険です。このツールでさまざまな分布を生成し、「もしこの分布の粉体を異なる原理の装置で測ったら、報告されるd50はどう変わるか？」と思考実験してみると、測定技術の本質的な理解に繋がります。