薄層乾燥 Page モデルシミュレーター — 食品・農産物の乾燥曲線

Q: Page モデルと Lewis（Newton）モデルの違いは何ですか？

Lewis モデルは MR = exp(-k·t) と1次反応で書ける最も単純な薄層乾燥モデルで、Page モデル MR = exp(-k·t^n) で n=1 に固定したものに相当します。実際の食品・農産物では乾燥初期に表面水分が急速に蒸発し、その後の内部拡散律速期にカーブの曲率が変化するため n≠1 になることが多く、Page モデルではこの曲率を指数 n で吸収できます。米やとうもろこしでは n≈1.0–1.3、葉物野菜や果実スライスでは n≈0.7–1.0 が報告されており、Lewis よりはるかに実測値とよく一致します。

Q: k と n はどう決めればよいですか？

実測の MR–時間データを ln(-ln MR) = ln k + n·ln t と二重対数で線形化し、最小二乗で傾き n と切片 ln k を求めるのが標準手順です。直接非線形最小二乗（Levenberg–Marquardt 等）で MR = exp(-k·t^n) に当てはめてもよく、こちらの方が小数点の精度が高くなります。同じ材料でも温度・風速・厚みごとに別個の (k, n) を求め、温度依存性は k = k0·exp(-Ea/RT) と Arrhenius 則で外挿します。Ea は食品で 20–50 kJ/mol 程度になることが多いです。

Q: 実効拡散係数 Deff の物理的な意味は何ですか？

Deff は Fick の第2法則を球または無限平板に適用したときに、内部の水分移動全体（液状水拡散、毛細管現象、表面拡散、蒸気拡散）をひとつの等価拡散係数にまとめたものです。本ツールでは球近似で Deff = k·r²/π² として Page 定数から逆算しています。食品では Deff は 10⁻¹¹ から 10⁻⁸ m²/s のレンジに収まり、米や穀類で 10⁻¹⁰ オーダー、果実スライスで 10⁻⁹ オーダーが典型値。Deff が温度とともに Arrhenius 則で増加するかを確認することで、データの妥当性をチェックできます。

Q: 褐変リスクスコアはどう読めばよいですか？

本ツールの褐変スコアは「空気温度が 50°C を超えた量 ×0.1」と「初期含水率が 100%db 未満の場合に +1」を加算した簡易プロキシで、Maillard 反応・酵素的褐変・カラメル化を厳密に予測するものではありません。スコア 1 以下は概ね安全（果実低温乾燥や凍結乾燥代替）、1–2 は要注意（コーヒー豆焙煎前乾燥・りんごスライス）、2 以上は色変化を許容する用途（穀類乾燥・乾物）と読みます。色管理が必須な場合は実機で L*a*b* 値を測定し、温度上限を別途設定してください。

薄層乾燥 Page モデルシミュレーター — 食品・農産物の乾燥曲線

米・小麦・りんご・コーヒー豆など食品・農産物の薄層乾燥を、半経験的な Page モデル MR = exp(-k·t^n) で予測します。初期含水率・平衡含水率・温度・風速を変えて、半乾燥時間・実効拡散係数・必要エネルギー・褐変リスクをリアルタイムに把握できます。

パラメータ設定

初期含水率 M₀

%db

乾燥前の含水率（乾量基準）。200%db = 水/乾物 2:1

平衡含水率 Mₑ

%db

空気の相対湿度と温度で決まる到達下限

Page 定数 k

1/min

乾燥動力学の時間定数。温度・風速で増加

Page 指数 n

n>1 で初期急乾燥、n<1 で遅い乾燥

空気温度 T

°C

空気流速 v

m/s

サンプル質量 m

計算結果

—

半乾燥時間 t(MR=0.5) (min)

—

5%乾燥時間 t(MR=0.05) (min)

—

実効拡散係数 Deff (m²/s)

—

蒸発水分 (g)

—

必要エネルギー (kJ)

—

褐変リスクスコア

—

試料断面アニメーション — 含水率カラーマップ

球状試料の中心→表面の含水率分布と空気流をリアルタイムに表示。色：青＝湿潤、赤＝乾燥近傍。経過時間と MR がライブ更新されます。

乾燥曲線 — MR vs 時間（Page モデル）

Page 定数 k の温度依存性（Arrhenius 風）

理論・主要公式

$$MR = \frac{M - M_e}{M_0 - M_e} = \exp(-k\,t^{n})$$

Page モデル。MR：モイスチャ比（無次元）、M：時刻 t の含水率、M₀：初期含水率、Mₑ：平衡含水率、k：動力学定数 [1/min]、n：Page 指数（食品種類・温度依存）。

$$t_{1/2} = \left(\frac{\ln 2}{k}\right)^{1/n}, \qquad t_{5\%} = \left(\frac{\ln 20}{k}\right)^{1/n}$$

半乾燥時間（MR=0.5）と 5% 残存時間（MR=0.05）。Page モデルから解析的に得られる。

$$D_{\mathrm{eff}} \approx \frac{k\,r^{2}}{\pi^{2}}, \qquad E_{\mathrm{evap}} = m_w\cdot\lambda_v$$

球近似の実効拡散係数 Deff（r：等価球半径、密度 1300 kg/m³ 仮定）と蒸発所要エネルギー（mw：蒸発水分質量、λv≈2.26 kJ/g の潜熱）。

乾燥動力学 Page モデルによる薄層乾燥

🙋

「Page モデル」って初めて聞きました。乾燥って単純に「時間が経てば水が抜ける」だけじゃないんですか？

🎓

うん、まさにそこが Page モデルが生まれた背景だよ。最初は Lewis さん（1921年）が「乾燥は1次反応」と仮定して MR = exp(-k·t) という単純な式を出した。でも米とかコーヒー豆を実測すると、初期は急に乾いて後半は粘る、という S 字に近いカーブが出てしまう。そこで Page さん（1949年）が指数を t^n に一般化して、たった2パラメータ (k, n) で実測カーブのほとんどを再現できるようにしたんだ。半経験式なんだけど、これがめちゃくちゃ強力で、今でも食品工学の論文で一番引用されているモデルだよ。

🙋

n の値で曲線の形が変わるのが面白いです。n を上げると初期にガッと落ちて、下げるとだらーっと長引きますね。これって物理的にどういう意味があるんですか？

🎓

いい観察だね。n が 1 を超える材料は「初期に表面水（自由水）が大量にあって蒸発しやすい」ものが多い。例えば葉物野菜の表面、果実のスライス断面なんかは n が 1.0〜1.3 になりやすい。逆に n が 1 未満の材料は「最初から内部拡散律速」で、稲穀粒のように硬い殻に覆われていたり、コーヒー豆のように内部組織が緻密だったりすると n が 0.6〜0.9 になる。だから n を測ること自体が、その材料の水分移動メカニズムを推定するヒントになるんだ。

🙋

空気温度を上げるとグラフの右下にある k の Arrhenius プロットが伸びてきます。実際の現場ではどこまで温度を上げていいんですか？

🎓

それが「乾燥は早ければ早いほど良い」とは限らない理由なんだ。確かに k = k₀·exp(-Ea/RT) で温度を 10°C 上げると k は 1.5〜2 倍になる。でも食品の場合、60°C を超えると Maillard 反応（褐変）が加速し、70°C を超えると熱に弱いビタミン C やフレーバー成分が壊れる。だから「乾燥時間が許容範囲に収まる最低温度」を狙うのが基本。右上の褐変リスクスコアはそのトレードオフを荒く見るための指標で、スコア 2 以上は色重視の用途では NG だよ。

🙋

実効拡散係数 Deff が「球近似で k から逆算」になっていますが、これって精度大丈夫なんですか？平板や円柱で乾かす場合はどうすれば？

🎓

球近似は穀粒・豆類・浆果系で 30〜50% の精度で当てはまる目安だね。形状が違うときは式の係数が変わって、無限平板なら Deff = k·L²/π²（L は半厚）、無限円柱なら Deff = k·r²/5.78 を使う。本ツールは「乾燥時間とエネルギーをサクッと見積もる」位置づけだから球で固定しているけど、論文用に精度を上げたいときは実際の形状で Crank の解（Fick の第2法則の長級数解）を使うか、CFD/熱物質移動シミュレーションで物体内部の水分プロファイルを直接解く必要があるよ。

🙋

必要エネルギーが意外と大きい（100g で 441 kJ）のに驚きました。これが食品工場の電気代に響くんですよね？

🎓

そう、乾燥は食品プロセスで最もエネルギーを食う工程の一つで、製造コストの 20〜30% を占めることもある。だから熱風乾燥にヒートポンプを組み合わせて排気から熱回収したり、減圧乾燥で潜熱を下げたり、過熱蒸気乾燥で空気酸化を避けつつ熱伝達を上げたり、と省エネ技術が日進月歩で発展している。本ツールが出すのはあくまで「水を蒸発させる理論最小エネルギー」で、実際の機械では熱効率 30〜60% 程度に落ちるから、実エネルギーは 2〜3 倍くらいに見積もっておくと現実的だね。

よくある質問

Lewis モデルは MR = exp(-k·t) と1次反応で書ける最も単純な薄層乾燥モデルで、Page モデル MR = exp(-k·t^n) で n=1 に固定したものに相当します。実際の食品・農産物では乾燥初期に表面水分が急速に蒸発し、その後の内部拡散律速期にカーブの曲率が変化するため n≠1 になることが多く、Page モデルではこの曲率を指数 n で吸収できます。米やとうもろこしでは n≈1.0–1.3、葉物野菜や果実スライスでは n≈0.7–1.0 が報告されており、Lewis よりはるかに実測値とよく一致します。

実測の MR–時間データを ln(-ln MR) = ln k + n·ln t と二重対数で線形化し、最小二乗で傾き n と切片 ln k を求めるのが標準手順です。直接非線形最小二乗（Levenberg–Marquardt 等）で MR = exp(-k·t^n) に当てはめてもよく、こちらの方が小数点の精度が高くなります。同じ材料でも温度・風速・厚みごとに別個の (k, n) を求め、温度依存性は k = k0·exp(-Ea/RT) と Arrhenius 則で外挿します。Ea は食品で 20–50 kJ/mol 程度になることが多いです。

Deff は Fick の第2法則を球または無限平板に適用したときに、内部の水分移動全体（液状水拡散、毛細管現象、表面拡散、蒸気拡散）をひとつの等価拡散係数にまとめたものです。本ツールでは球近似で Deff = k·r²/π² として Page 定数から逆算しています。食品では Deff は 10⁻¹¹ から 10⁻⁸ m²/s のレンジに収まり、米や穀類で 10⁻¹⁰ オーダー、果実スライスで 10⁻⁹ オーダーが典型値。Deff が温度とともに Arrhenius 則で増加するかを確認することで、データの妥当性をチェックできます。

本ツールの褐変スコアは「空気温度が 50°C を超えた量 ×0.1」と「初期含水率が 100%db 未満の場合に +1」を加算した簡易プロキシで、Maillard 反応・酵素的褐変・カラメル化を厳密に予測するものではありません。スコア 1 以下は概ね安全（果実低温乾燥や凍結乾燥代替）、1–2 は要注意（コーヒー豆焙煎前乾燥・りんごスライス）、2 以上は色変化を許容する用途（穀類乾燥・乾物）と読みます。色管理が必須な場合は実機で L*a*b* 値を測定し、温度上限を別途設定してください。

実世界での応用

米・穀類の収穫後乾燥：収穫直後の籾は含水率 20〜25%wb（≈25〜33%db）あり、長期保管のためには 14%wb 以下まで乾燥する必要があります。Page モデルで k と n を求め、循環式乾燥機や横流式乾燥機の風量・温度・滞留時間を設計します。割れ米（胴割れ）を防ぐため、急乾燥を避けて 30〜40°C で時間をかける「省エネ低温乾燥」が近年は主流。AACC 法・AOAC 法に乾燥試験の標準手順が定められています。

コーヒー豆・カカオ豆の発酵後乾燥：パーチメントコーヒーは収穫後の発酵を経て含水率 60%wb 程度から 11〜12%wb まで天日または機械乾燥します。Page モデルで n が 0.7〜0.9 と1未満になることが多く、内部の緻密組織が拡散律速になっていることを示唆します。乾燥が速すぎると割れや風味劣化（オフフレーバー）を招くため、CFD と組み合わせた均一乾燥プロファイル設計が重要。

果実・野菜のスライス乾燥（ドライフルーツ・乾燥野菜）：りんご・マンゴー・トマトなどのスライス乾燥は、初期含水率 400〜800%db から目標 15〜25%db まで落とします。空気温度 50〜70°C、風速 1〜3 m/s が一般的で、褐変防止のため抗酸化処理（亜硫酸塩浸漬）や予備乾燥（湿度制御）を組み合わせます。Page モデルの (k, n) からエネルギー原単位を推定し、ヒートポンプ乾燥との比較に使えます。

医薬品錠剤・顆粒の流動層乾燥：製薬では造粒後の顆粒を流動層乾燥機で含水率数% まで仕上げます。Page モデルの k を温度・気流速度の関数として相関化し、バリデーション試験のスケールアップに活用。GMP の観点から、各バッチで MR の対数プロットが直線になるか（モデル適合度 R²>0.99）が品質管理指標になります。

よくある誤解と注意点

第一の落とし穴は、「Page モデルは万能ではなく、薄層（thin-layer）の仮定が成立する条件でのみ使える」という点です。本モデルは試料内部の温度勾配と外気との抵抗を無視し、内部水分のみを律速段階と見なします。試料厚みが数 mm を超える厚層（穀粒の堆積層、果実の塊乾燥）では、層内部の空気温湿度変化が無視できなくなり、Page モデルの k がパラメータ化できなくなります。厚層乾燥では Thompson モデルや層モデル（多孔質体の物質移動 PDE を解く）に切り替える必要があります。

第二に、「平衡含水率 Mₑ を一定値と仮定する」ことの危険性。Mₑ は空気の相対湿度・温度で決まる吸湿等温線上の値で、本来は時間とともに変化する空気状態の関数です。バッチ乾燥機では乾燥が進むにつれ排気の湿度が下がり、Mₑ も低下します。Mₑ を初期値で固定して MR を計算すると、後半の乾燥速度を過小評価し、必要時間を 10〜20% 短く見積もる誤りが起きます。連続式乾燥機では Mₑ をほぼ一定と見なせますが、バッチでは GAB モデルや BET モデルで湿度依存性を組み込むべきです。

第三に、「実効拡散係数 Deff の単位を勘違いしやすい」こと。本ツールでは k [1/min] から Deff [m²/s] を換算する際、内部で 60 で割って単位を整合させています（k [1/s] = k [1/min] / 60）。論文で Deff を引用するときは、出典の単位系（cm²/s なのか m²/s なのか、min ベースなのか s ベースなのか）を必ず確認してください。同じ材料・条件で報告された Deff が 1〜2 桁ずれている場合、ほぼ間違いなく単位ミスです。また、本ツールの球近似 r²/π² は形状係数で、平板なら L²/π²·(1/4)、円柱なら r²/5.78 と係数が変わるため、形状を変えたら式自体を入れ替える必要があります。

薄層乾燥 Page モデルシミュレーター — 食品・農産物の乾燥曲線

乾燥動力学 Page モデルによる薄層乾燥

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点