ダムケラー数シミュレーター 戻る
反応工学

ダムケラー数シミュレーター — 反応と滞留時間

化学反応器の性能を支配する無次元数「ダムケラー数 Da」を体験するツールです。反応速度定数・滞留時間・入口濃度を変えると、CSTR と PFR の転化率と出口濃度がリアルタイムで分かり、反応と流れのどちらが律速になっているかを直感的につかめます。

パラメータ設定
反応器タイプ
完全混合の CSTR か、押し出し流れの PFR か
反応次数
反応速度の濃度依存性
反応速度定数 k
1次は 1/s、2次は m³/(mol·s) の単位
滞留時間 τ
s
反応器体積 ÷ 体積流量。大きいほど長く反応
入口濃度 C₀
mol/m³
2次反応では Da と転化率に直接効く
計算結果
ダムケラー数 Da
転化率 X (%)
出口濃度 C (mol/m³)
入口反応速度 r₀ (mol/m³·s)
他方式の転化率(比較)(%)
律速の判定
反応器内の流れと濃度 — アニメーション

左から原料(青)が入り、反応して生成物(橙)に変わりながら右へ流れます。色の濃さが濃度を表し、CSTR は均一、PFR は入口から出口へグラデーションになります。

転化率 vs ダムケラー数(CSTR・PFR 比較)
反応器内の濃度プロファイル
理論・主要公式

$$Da=k\tau\ (\text{1次}),\qquad Da=kC_0\tau\ (\text{2次})$$

ダムケラー数 Da は「反応速度」を「対流による物質輸送速度」で割った無次元数。k:反応速度定数、τ:滞留時間、C₀:入口濃度。

$$X_{CSTR}=\frac{Da}{1+Da},\qquad X_{PFR}=1-e^{-Da}\quad(\text{1次})$$

1次反応の転化率 X。同じ Da でも PFR は CSTR より高い転化率を達成する(PFR は高濃度域を通るため)。

$$Da(1-X)^2=X\ \Rightarrow\ X=\frac{1+2Da-\sqrt{1+4Da}}{2Da}\quad(\text{2次CSTR})$$

2次反応 CSTR の転化率は物質収支の2次方程式を解いて求める。出口濃度は C = C₀(1−X)。

ダムケラー数とは

🙋
「ダムケラー数」って初めて聞きました。レイノルズ数みたいな無次元数の仲間ですか?
🎓
そう、化学反応器の世界でいちばん大事な無次元数だよ。ざっくり言うと「反応がどれだけ速いか」を「流れがどれだけ速く流体を運び去るか」で割ったもの。Da = kτ で、k が反応の速さ、τ が滞留時間——流体が反応器の中に居座る時間だ。Da が大きいほど、流体が出ていく前にしっかり反応できる、というわけ。
🙋
なるほど。じゃあ Da が小さいとどうなるんですか?左のスライダーで τ をうんと小さくしてみました。
🎓
いまの操作で転化率がガクッと下がっただろう? Da ≪ 1 は「滞留時間律速」——流体が反応する暇もなく反応器を通り抜けてしまう状態だ。原料がほとんど未反応のまま出口に出てくる。直す手は2つ。反応器を大きくして τ を稼ぐか、流量を絞って τ を伸ばすか。逆に Da ≫ 1 なら「反応律速」で、ほぼ完全に反応しきる。だから Da を見れば「反応器が小さすぎないか」がひと目で分かるんだ。
🙋
反応器タイプを CSTR から PFR に切り替えると、同じ条件なのに転化率が上がりました。これはどうしてですか?
🎓
いい発見だ。CSTR は「完全混合」、つまりタンクの中身が一瞬で混ざって全部が出口濃度になる。反応速度は濃度が高いほど速いから、CSTR は常に「いちばん薄くて遅い」状態で運転していることになる。一方 PFR は管の中を押し出し流れで進むので、入口は濃くて速く、出口に近づくほど薄くなる。高濃度=高速の区間を通る分、同じ Da でも PFR のほうが転化率は高い。グラフの2本の曲線がそれを示しているよ。
🙋
2次反応に切り替えると、入口濃度 C₀ を上げただけで転化率が変わりました。1次のときは変わらなかったのに。
🎓
そこが反応次数の効いてくるところだ。1次反応は Da = kτ で、C₀ は式に出てこない——濃度が2倍になっても反応速度も2倍になり、転化率という「割合」は変わらない。でも2次反応は Da = kC₀τ。速度が濃度の2乗で効くので、C₀ を上げると Da が増えて転化率も上がる。だから2次反応では「濃く仕込む」こと自体が立派な設計手段になるんだ。
🙋
じゃあ Da をとにかく大きくすれば良いツールが作れそうですね?
🎓
気持ちは分かるけど、そう単純でもない。Da が 10 を超えると転化率は 99% 台でほぼ頭打ち——それ以上 τ を伸ばしても改善はわずかなのに、反応器はどんどん大きく・高コストになる。しかも滞留時間が長すぎると、目的物がさらに反応してしまう逐次反応や副反応が進んで、選択率が落ちることもある。実務では Da を 1〜10 くらいの「拮抗領域」に置く設計が多い。グラフの曲線が急に立ち上がる、いちばんおいしいところだね。

よくある質問

ダムケラー数 Da は「反応の速さ」を「流れ(対流)による物質輸送の速さ」で割った無次元数です。1次反応では Da = kτ、2次反応では Da = kC₀τ で計算します(k は反応速度定数、τ は滞留時間、C₀ は入口濃度)。Da ≪ 1 なら流体は反応する前に反応器を出てしまい、Da ≫ 1 なら反応器内でほぼ完全に反応します。Da は反応器を大きくすべきか流量を下げるべきかを判断する出発点になります。
同じ Da なら PFR(管型反応器)の転化率は常に CSTR(連続槽型反応器)より高くなります。CSTR は完全混合のため、反応器全体が低い出口濃度で運転され、反応速度が濃度に依存する反応では速度が最も遅い状態で進みます。一方 PFR は入口から出口へ濃度が連続的に下がり、高濃度(=高速)の領域を通るため、同じ滞留時間でより多く反応できます。1次反応では CSTR は X=Da/(1+Da)、PFR は X=1−e^(−Da) です。
転化率は Da で決まるので、Da = kτ(1次)を上げれば改善します。滞留時間 τ を上げる手段は反応器を大きくする・流量を下げることで、設備投資や生産量に直結します。反応速度定数 k を上げる手段は温度を上げる・触媒を使うことで、運転コストや選択率に影響します。Da ≪ 1(滞留時間律速)なら τ を増やすほうが効きやすく、すでに Da が大きい領域では τ を増やしても転化率はほとんど飽和します。
Da ≪ 0.1 だと転化率がごくわずかで、未反応原料が大量に出口へ流れます。反応器を大型化するか流量を下げる必要があります。逆に Da ≫ 10 ではほぼ完全転化となり、それ以上 Da を上げても転化率の改善はほぼ見込めません。むしろ滞留時間が長すぎると副反応・逐次反応が進んで目的生成物の選択率が下がる、装置が無駄に大きくなる、といった弊害が出ます。実務では Da を 1〜10 程度の拮抗領域に置く設計が多く見られます。

実世界での応用

化学プラントの反応器設計:連続プロセスで反応器の体積を決める最初の計算がダムケラー数です。目標転化率(例えば 90%)から必要な Da を逆算し、反応速度定数 k と要求生産量(=体積流量)から反応器体積 V = τ·Q を求めます。CSTR にするか PFR にするか、PFR を多段にするかといった構成の選択も、Da と転化率の関係曲線を見ながら判断します。

反応器タイプの選定:同じ Da で PFR が CSTR より高い転化率を出すという性質は、装置選定の基本原理です。高い転化率が必要なら PFR、強い発熱で温度制御が重要なら混合のよい CSTR、というトレードオフを Da ベースで定量比較します。実際には CSTR を直列に並べて PFR に近い性能を狙う「CSTR カスケード」もよく使われます。

燃焼・燃焼器の設計:ガスタービンや工業炉の燃焼器では、化学反応の時間スケールと流体が燃焼室に留まる時間スケールの比としてダムケラー数が使われます。Da が小さい(流れが速すぎる)と火炎が吹き消える「ブローオフ」が起き、Da が大きいと安定燃焼します。燃焼器の保炎設計はまさにダムケラー数の管理です。

環境・廃水処理:活性汚泥法の曝気槽や触媒燃焼による排ガス処理装置も、滞留時間と反応速度のバランスで分解率が決まります。処理しきれない場合、Da が小さすぎる=槽が小さい/流量が多すぎるサインで、増設や流量調整を Da で評価します。バイオリアクターの設計指標としても重要です。

よくある誤解と注意点

まず多いのが、「ダムケラー数の定義は一つだけだと思い込む」ことです。本ツールが扱うのは「反応速度 ÷ 対流輸送速度」の Da(しばしば Da_I と呼ばれる)ですが、ほかにも「反応速度 ÷ 拡散輸送速度」を表す Da_II が存在します。物質移動が律速になる不均一系触媒反応や、多相反応では Da_II が主役になります。論文や教科書でダムケラー数を見たら、それがどの輸送過程との比なのかを必ず確認してください。定義を取り違えると、まったく違う結論に至ります。

次に、「滞留時間 τ を平均値だけで扱う」こと。本ツールの PFR は理想的な押し出し流れ、CSTR は理想的な完全混合を仮定しています。しかし実際の反応器には流れの偏り・短絡・デッドゾーンがあり、流体粒子ごとに反応器に留まる時間がばらつきます。この「滞留時間分布(RTD)」を無視して平均 τ だけで設計すると、実機の転化率が計算値を下回ることがあります。非理想流れの評価にはトレーサー試験による RTD 測定が欠かせません。

最後に、「Da が大きいほど良い設計だ」という思い込み。Da を上げれば転化率は上がりますが、それは目的物の収率や選択率が上がることを意味しません。目的物がさらに反応する逐次反応(A→B→C で B が欲しい場合など)では、滞留時間を長くしすぎると B が C に変わってしまい、最適 Da は「ほどほど」の値になります。さらに発熱反応では Da を上げすぎると除熱が追いつかず、反応暴走(ランナウェイ)の危険もあります。Da は転化率の指標であって、収率・安全性・コストまで保証する万能指標ではない、という点を忘れないでください。

使い方ガイド

  1. 反応速度定数k(m³/mol·s)と滞留時間τ(s)を入力欄に設定します。CSTR・PFR両反応器について、ダムケラー数Da = kC₀τで自動計算されます
  2. 初期濃度C₀(mol/m³)を指定すると、入口反応速度r₀ = kC₀が決定され、各反応器における転化率X(%)と出口濃度C(mol/m³)がリアルタイムで算出されます
  3. Da値の大小により反応律速か拡散律速かを判定し、CSTR・PFR間の転化率差異を比較グラフで確認します。一次反応を仮定し、CSTR: X = Da/(1+Da)、PFR: X = 1-exp(-Da)の式を適用します

具体的な計算例

硝酸化反応(尿素からの)でk = 0.015 m³/mol·s、τ = 120 s、C₀ = 500 mol/m³の条件でシミュレーションを実行した場合、Da = 0.015 × 500 × 120 = 900となります。CSTR転化率は99.89%、出口濃度0.55 mol/m³に対し、PFR転化率は99.99%、出口濃度0.05 mol/m³となり、PFRの方が高い転化率を示します。r₀ = 0.015 × 500 = 7.5 mol/m³·sです

実務での注意点

  1. 触媒充填層リアクターではDa > 10で反応支配、0.1 < Da < 10で中間領域、Da < 0.1で拡散支配となるため、触媒粒径と流速の最適化が必須です
  2. 水処理プラント(生物反応槽)ではCSTRが採用され、HRT(滞留時間)を4~8時間に設定することでBOD除去率90%以上を達成します
  3. 有機合成での逐次反応では、中間生成物の蓄積を避けるためPFRが有利ですが、発熱反応の場合は温度制御が困難なため多段CSTR直列が標準です