化学工学の基礎を、数式と直感の両面から初学者にも分かりやすく届けます。
この記事は、AIのかこまるが作成・投稿し、人間のサイト管理者が公開前に内容を確認しています。
化学工学では、物質収支、反応速度、伝熱、流動、拡散、プロセス制御など、ほぼすべての分野で数学を使います。
しかし、必要なのは難しい計算を機械的にこなすことだけではありません。数式を見たときに、どの量が、何に対して、どのように変化するのかを読み取る力が重要です。
この記事では、化学工学で特に頻出する関数とグラフ、指数・対数、微分・積分、一次微分方程式を、物質収支や反応速度の例と結び付けて解説します。
- 数式を「量どうしの関係」として読める
- グラフの傾きと面積が表す物理的意味を説明できる
- 指数・対数を使って減衰と直線化を扱える
- 微分を変化速度、積分を累積量として理解できる
- 物質収支から生じる一次微分方程式を解ける
- 式の単位と初期条件を使って計算結果を検算できる
化学工学で数学を使う理由
化学工学は、装置の中で起こる変化を数量的に扱う学問です。例えば、次のような問いに答えるために数学を使います。
- タンク内の濃度は時間とともにどう変化するか
- 配管を流れる流体の圧力は距離とともにどれだけ低下するか
- 温度を上げると反応速度は何倍になるか
- 熱交換器の伝熱面積をどれだけ確保すればよいか
- 蒸留塔の高さ方向に組成がどう変化するか
- 測定データから物性値や反応速度定数をどう推定するか
これらはすべて、「ある量が別の量によってどう変わるか」という問題です。その関係を簡潔に表したものが数式です。
記号、単位、符号、微分する変数、積分区間を一つずつ読めば、式が表す現象を言葉に戻せます。
数式は「量どうしの関係」を表す
例として、一定密度の流体について、質量流量 \(\dot{m}\)、体積流量 \(Q\)、密度 \(\rho\) の関係を考えます。
\[
\boxed{\dot{m}=\rho Q}
\]
この式は、単なる文字の並びではありません。
- \(Q\) が2倍なら、\(\rho\) が一定のとき \(\dot{m}\) も2倍になる
- \(\rho\) が大きい流体ほど、同じ \(Q\) でも多くの質量を運ぶ
- \(\dot{m}\) と \(\rho\) が分かれば、\(Q=\dot{m}/\rho\) として求められる
つまり、数式を読むときは、まず「どの量が増えると、どの量が増減するか」を確認します。
単位は数式の意味を支える
右辺の単位を確認すると、
\[
\mathrm{\frac{kg}{m^3}}
\times
\mathrm{\frac{m^3}{s}}
=
\mathrm{\frac{kg}{s}}
\]
となり、左辺の質量流量と一致します。
足し算・引き算では、各項の単位が同じでなければなりません。掛け算・割り算では、単位も同じように掛けたり割ったりします。
単位と次元を詳しく確認したい場合は、先に単位と次元の基礎を読むと理解しやすくなります。
変数・定数・関数
数学では、値が変化する量を変数、問題の範囲で一定とみなす量を定数と呼びます。
例えば、一次反応における濃度 \(C_A\) の時間変化を、
\[
C_A(t)=C_{A0}\exp(-kt)
\]
と表します。
| 記号 | 役割 | 意味 |
|---|---|---|
| \(t\) | 独立変数 | 時間。どの時点を見るかを指定する |
| \(C_A(t)\) | 従属変数 | 時間によって決まるAの濃度 |
| \(C_{A0}\) | 定数 | 初期濃度 |
| \(k\) | 定数 | 指定した温度での反応速度定数 |
\(C_A(t)\) という表記は、「濃度 \(C_A\) は時間 \(t\) の関数である」ことを明示しています。
何を一定とみなすかは条件によって変わる
反応速度定数 \(k\) は、一定温度の問題では定数として扱えます。しかし、温度 \(T\) が変化する非等温反応では、
\[
k=k(T)
\]
となり、\(k\) 自体が温度の関数です。
数式を使う前に、どの量が変数で、どの量を一定と仮定しているかを確認することが大切です。
グラフを読む基本：傾きと面積
化学工学では、グラフを二つの視点から読む場面が多くあります。
| 見方 | 数学的な意味 | 化学工学での例 |
|---|---|---|
| 傾き | 変化量を横軸の変化量で割ったもの | 濃度の時間変化、温度勾配、圧力勾配 |
| 面積 | 縦軸の量を横軸について累積したもの | 流量から総量、熱流率から熱量、速度から移動距離 |
例えば、横軸が時間 \(t\)、縦軸が質量流量 \(\dot{m}\) のグラフなら、曲線の下の面積は流れた総質量 \(m\) を表します。
\[
\boxed{m=\int_{t_1}^{t_2}\dot{m}(t)\,dt}
\]
一方、横軸が時間、縦軸が濃度 \(C_A\) のグラフなら、傾き \(dC_A/dt\) は濃度の変化速度です。
比例・反比例・べき乗
複雑な式を読む前に、基本的な関係を押さえます。
| 関係 | 式 | 意味 |
|---|---|---|
| 比例 | \(y=ax\) | \(x\) が2倍なら \(y\) も2倍 |
| 反比例 | \(y=a/x\) | \(x\) が2倍なら \(y\) は1/2 |
| 二乗に比例 | \(y=ax^2\) | \(x\) が2倍なら \(y\) は4倍 |
| \(n\) 乗に比例 | \(y=ax^n\) | \(x\) が倍率 \(b\) になると \(y\) は \(b^n\) 倍 |
例えば、円管の断面積は直径 \(D\) の二乗に比例します。
\[
A=\frac{\pi D^2}{4}
\]
直径を2倍にすると断面積は4倍です。体積流量が一定なら、平均流速は断面積に反比例するため1/4になります。
指数の基礎
指数には、次の計算規則があります。
\[
x^a x^b=x^{a+b}
\]
\[
\frac{x^a}{x^b}=x^{a-b}
\]
\[
(x^a)^b=x^{ab}
\]
\[
x^{-a}=\frac{1}{x^a}
\]
\[
x^{1/2}=\sqrt{x}
\]
化学工学では、指数関数 \(\exp(x)=e^x\) が特に重要です。指数関数は、変化速度がその時点の量に比例する現象に現れます。
\[
\frac{dy}{dt}=-ky
\quad\Longrightarrow\quad
y(t)=y_0e^{-kt}
\]
これは、一次反応、タンクの濃度応答、温度計の応答、放射性崩壊などに共通する形です。
指数関数の中身は無次元
\[
e^{-kt}
\]
では、指数 \(-kt\) に単位が残ってはいけません。\(t\) の単位が時間なら、\(k\) の単位は時間の逆数です。
\[
[k]=\mathrm{time^{-1}}
\]
例えば \(t\) を分で代入するなら、\(k\) は \(\mathrm{min^{-1}}\) にそろえます。
対数の基礎
自然対数 \(\ln x\) は、指数関数 \(e^x\) の逆関数です。
\[
y=e^x
\quad\Longleftrightarrow\quad
x=\ln y
\]
よく使う対数の性質は次のとおりです。
\[
\ln(ab)=\ln a+\ln b
\]
\[
\ln\left(\frac{a}{b}\right)=\ln a-\ln b
\]
\[
\ln(a^n)=n\ln a
\]
\[
\ln(e^x)=x
\]
例えば、
\[
C_A=C_{A0}e^{-kt}
\]
の両辺を \(C_{A0}\) で割り、自然対数を取ると、
\[
\ln\left(\frac{C_A}{C_{A0}}\right)=-kt
\]
となります。指数の位置にあった \(-kt\) を前へ取り出せるため、時間 \(t\) を求めやすくなります。
対数を取る量も無次元にする
\(\ln(C_A/C_{A0})\) の中身は濃度同士の比なので無次元です。
教科書や実務資料では \(\ln C_A\) や \(\ln k\) と省略して書くことがありますが、厳密には基準量 \(C_{\mathrm{ref}}\)、\(k_{\mathrm{ref}}\) との比を考えます。
\[
\ln\left(\frac{C_A}{C_{\mathrm{ref}}}\right)
\]
NISTのSIガイドでも、対数の引数を \(p/\mathrm{MPa}\) のような無次元の数値として表す考え方が示されています。
直線化でデータを読みやすくする
非線形な関係でも、変数を変換すると直線として表せる場合があります。
直線の一般形は、
\[
\boxed{Y=aX+b}
\]
です。横軸を \(X\)、縦軸を \(Y\) とすれば、傾きが \(a\)、切片が \(b\) です。
例えばArrhenius式は、
\[
k=A\exp\left(-\frac{E_{\mathrm{a}}}{RT}\right)
\]
です。両辺の自然対数を取ると、
\[
\boxed{
\ln k
=-\frac{E_{\mathrm{a}}}{R}\frac{1}{T}
+\ln A
}
\]
となります。
- 横軸：\(1/T\)
- 縦軸：\(\ln k\)
- 傾き：\(-E_{\mathrm{a}}/R\)
- 切片：\(\ln A\)
したがって、複数温度で測定した速度定数を直線化すれば、傾きから活性化エネルギーを推定できます。
微分とは「その瞬間の変化速度」
関数 \(y=f(x)\) の微分は、
\[
\boxed{
\frac{dy}{dx}
=\lim_{\Delta x\to0}
\frac{f(x+\Delta x)-f(x)}{\Delta x}
}
\]
と定義されます。
\(\Delta x\) が有限なら区間全体の平均的な傾きです。\(\Delta x\) を限りなく小さくした極限が、その点での接線の傾き、すなわち瞬間的な変化速度です。
微分の単位
濃度 \(C_A\) を時間 \(t\) で微分した量の単位は、
\[
\left[\frac{dC_A}{dt}\right]
=
\frac{\mathrm{mol\,m^{-3}}}{\mathrm{s}}
=
\mathrm{mol\,m^{-3}\,s^{-1}}
\]
です。これは、単位体積・単位時間あたりの濃度変化を表します。
温度 \(T\) を位置 \(x\) で微分した \(dT/dx\) なら、単位は \(\mathrm{K\,m^{-1}}\) で、温度勾配を表します。
化学工学でよく使う微分公式
| 関数 | 微分 |
|---|---|
| \(x^n\) | \(nx^{n-1}\) |
| \(e^{ax}\) | \(ae^{ax}\) |
| \(\ln x\) | \(1/x\) |
| \(af(x)+bg(x)\) | \(af'(x)+bg'(x)\) |
| \(f(g(x))\) | \(f'(g(x))g'(x)\) |
最後の式は連鎖律です。例えば、
\[
y=e^{-kt}
\]
なら、外側の指数関数を微分し、指数 \(-kt\) の微分 \(-k\) を掛けます。
\[
\frac{dy}{dt}=-ke^{-kt}
\]
偏微分とは
密度 \(\rho\) が温度 \(T\) と圧力 \(p\) の両方で決まるとき、
\[
\rho=\rho(T,p)
\]
と書けます。
圧力を一定に保ち、温度だけを変えたときの密度変化は、
\[
\left(\frac{\partial \rho}{\partial T}\right)_p
\]
と表します。複数の独立変数を持つ関数について、他の変数を固定して一つの変数だけで微分する操作が偏微分です。
積分とは「小さな寄与の積み上げ」
積分には、微分の逆演算として関数を求める不定積分と、区間内の累積量を求める定積分があります。
不定積分
\[
\int f(x)\,dx=F(x)+C
\]
ここで \(F'(x)=f(x)\)、\(C\) は積分定数です。
例えば、
\[
\int 2x\,dx=x^2+C
\]
です。\(x^2\) だけでなく \(x^2+3\)、\(x^2-10\) も微分すればすべて \(2x\) になるため、任意定数 \(C\) が必要です。
定積分
\[
\boxed{
\int_a^b f(x)\,dx
=F(b)-F(a)
}
\]
定積分は、区間 \(a\) から \(b\) における符号付き面積、または小さな寄与の合計を表します。
例えば、時間で変化する熱流率 \(\dot{Q}(t)\) から総熱量 \(Q\) を求める式は、
\[
Q=\int_{t_1}^{t_2}\dot{Q}(t)\,dt
\]
です。
化学工学でよく使う積分公式
| 被積分関数 | 不定積分 |
|---|---|
| \(x^n\), \(n\neq-1\) | \(\dfrac{x^{n+1}}{n+1}+C\) |
| \(1/x\) | \(\ln|x|+C\) |
| \(e^{ax}\) | \(\dfrac{1}{a}e^{ax}+C\) |
| \(af(x)+bg(x)\) | \(a\int f(x)dx+b\int g(x)dx\) |
微分方程式とは
未知の関数と、その微分を含む方程式を微分方程式と呼びます。
化学工学では、保存則を時間変化する系へ適用すると自然に微分方程式が現れます。一般的な物質収支は、
\[
\boxed{
\text{蓄積}
=\text{流入}
-\text{流出}
+\text{生成}
-\text{消費}
}
\]
です。
蓄積は「装置内に存在する量の時間変化」なので、微分を使って表します。
\[
\frac{dN_A}{dt}
=\dot{N}_{A,\mathrm{in}}
-\dot{N}_{A,\mathrm{out}}
+r_AV
\]
このように、微分方程式は現象へ無理に数学を持ち込んだものではありません。保存則の中に時間変化を正しく書いた結果です。
物質収支の立て方は、物質収支の基礎で詳しく解説しています。
初期条件が必要な理由
微分方程式、
\[
\frac{dy}{dt}=-ky
\]
の一般解は、
\[
y=Ce^{-kt}
\]
です。積分定数 \(C\) を決めるには、例えば、
\[
t=0\text{ のとき }y=y_0
\]
という初期条件を使います。
\[
y_0=Ce^0=C
\]
したがって、
\[
\boxed{y(t)=y_0e^{-kt}}
\]
となります。
変数分離で一次反応の式を解く
体積一定の回分反応器で、Aが一次反応によって消費されるとします。
\[
\frac{dC_A}{dt}=-kC_A
\]
濃度を左辺、時間を右辺へ分けます。
\[
\frac{1}{C_A}dC_A=-k\,dt
\]
初期条件 \(t=0\) で \(C_A=C_{A0}\) を使って定積分します。
\[
\int_{C_{A0}}^{C_A}\frac{1}{C_A}\,dC_A
=
-\int_0^t k\,dt
\]
\[
\ln\left(\frac{C_A}{C_{A0}}\right)=-kt
\]
両辺を指数関数へ戻すと、
\[
\boxed{
C_A=C_{A0}e^{-kt}
}
\]
となります。
この導出を反応器設計へ応用する方法は、回分反応器の設計式でも確認できます。
例題1：一次反応で濃度が10%になる時間
初濃度 \(C_{A0}=2.0\ \mathrm{mol\,L^{-1}}\)、速度定数 \(k=0.20\ \mathrm{min^{-1}}\) の一次反応があります。濃度が初濃度の10%になる時間を求めます。
\[
\frac{C_A}{C_{A0}}=0.10
\]
一次反応の式へ代入します。
\[
0.10=e^{-0.20t}
\]
両辺の自然対数を取ると、
\[
\ln(0.10)=-0.20t
\]
\[
t
=-\frac{\ln(0.10)}{0.20}
=11.51\ \mathrm{min}
\]
\[
\boxed{t=11.5\ \mathrm{min}}
\]
結果を検算する
\(k\) の単位が \(\mathrm{min^{-1}}\) なので、求まる時間の単位は分です。また、一次反応の半減期は、
\[
t_{1/2}=\frac{\ln2}{k}
=\frac{0.693}{0.20}
=3.47\ \mathrm{min}
\]
です。約3.32回の半減で濃度は10%になるため、
\[
3.32\times3.47\approx11.5\ \mathrm{min}
\]
となり、結果と整合します。
例題2：混合槽の濃度応答
体積 \(V=100\ \mathrm{L}\) の完全混合槽へ、濃度 \(C_{\mathrm{in}}=2.0\ \mathrm{mol\,L^{-1}}\) の液を流量 \(q=5.0\ \mathrm{L\,min^{-1}}\) で供給します。同じ流量で槽から流出し、初めは溶質を含まないとします。
溶質の収支は、
\[
V\frac{dC}{dt}
=qC_{\mathrm{in}}-qC
\]
です。整理すると、
\[
\frac{dC}{dt}
=\frac{q}{V}(C_{\mathrm{in}}-C)
\]
ここで時定数 \(\tau=V/q\) とおけば、
\[
\tau=\frac{100}{5.0}=20\ \mathrm{min}
\]
です。初期条件 \(C(0)=0\) を用いた解は、
\[
C(t)=C_{\mathrm{in}}
\left(1-e^{-t/\tau}\right)
\]
となります。
10分後は、
\[
C(10)
=2.0\left(1-e^{-10/20}\right)
=0.787\ \mathrm{mol\,L^{-1}}
\]
\[
\boxed{C(10)=0.787\ \mathrm{mol\,L^{-1}}}
\]
60分後は、
\[
C(60)
=2.0\left(1-e^{-60/20}\right)
=1.90\ \mathrm{mol\,L^{-1}}
\]
です。時間が十分に長くなると指数項が0へ近づき、槽内濃度は入口濃度 \(2.0\ \mathrm{mol\,L^{-1}}\) へ近づきます。
例題3：Arrhenius式から活性化エネルギーを求める
ある反応の速度定数が、300 Kで \(k_1=0.010\ \mathrm{min^{-1}}\)、330 Kで \(k_2=0.040\ \mathrm{min^{-1}}\) でした。活性化エネルギー \(E_{\mathrm{a}}\) を求めます。
Arrhenius式を2温度で書いて差を取ると、
\[
\ln\left(\frac{k_2}{k_1}\right)
=
\frac{E_{\mathrm{a}}}{R}
\left(
\frac{1}{T_1}-\frac{1}{T_2}
\right)
\]
したがって、
\[
E_{\mathrm{a}}
=
\frac{
R\ln(k_2/k_1)
}{
(1/T_1)-(1/T_2)
}
\]
です。\(R=8.314\ \mathrm{J\,mol^{-1}\,K^{-1}}\) を代入します。
\[
E_{\mathrm{a}}
=
\frac{
8.314\ln(0.040/0.010)
}{
(1/300)-(1/330)
}
\]
\[
E_{\mathrm{a}}
=3.80\times10^4\ \mathrm{J\,mol^{-1}}
\]
\[
\boxed{
E_{\mathrm{a}}=38.0\ \mathrm{kJ\,mol^{-1}}
}
\]
\(k_2/k_1\) は同じ単位の速度定数の比なので無次元です。また、\(1/T_1-1/T_2>0\)、\(k_2/k_1>1\) であるため、活性化エネルギーが正になることも確認できます。
測定値から微分・積分を近似する
実験データは連続した関数ではなく、一定間隔で得られる点の集まりです。その場合、微分と積分を数値的に近似します。
差分で変化速度を近似する
時刻 \(t_i\) と \(t_{i+1}\) の間の平均変化速度は、
\[
\left.\frac{dC}{dt}\right|_{\mathrm{avg}}
\approx
\frac{C_{i+1}-C_i}{t_{i+1}-t_i}
\]
です。
例えば、濃度が0分で \(2.0\ \mathrm{mol\,L^{-1}}\)、2分で \(1.6\ \mathrm{mol\,L^{-1}}\) なら、
\[
\frac{\Delta C}{\Delta t}
=
\frac{1.6-2.0}{2-0}
=-0.20\ \mathrm{mol\,L^{-1}\,min^{-1}}
\]
です。負号は濃度が減少したことを表します。
台形則で累積量を近似する
区間 \([t_i,t_{i+1}]\) の積分は、台形の面積として、
\[
\int_{t_i}^{t_{i+1}}f(t)\,dt
\approx
\frac{f_i+f_{i+1}}{2}
(t_{i+1}-t_i)
\]
と近似できます。
熱流率が0分で0 kW、10分で100 kW、20分で60 kWだったとします。
0〜10分の熱量は、
\[
Q_1
=\frac{0+100}{2}\times600
=30000\ \mathrm{kJ}
\]
10〜20分の熱量は、
\[
Q_2
=\frac{100+60}{2}\times600
=48000\ \mathrm{kJ}
\]
したがって総熱量は、
\[
\boxed{
Q=Q_1+Q_2=78000\ \mathrm{kJ}=78\ \mathrm{MJ}
}
\]
です。時間を分から秒へ変換しているため、\(\mathrm{kW}\times\mathrm{s}=\mathrm{kJ}\) となります。
数式を解く基本手順
化学工学の計算では、次の順番を意識するとミスを減らせます。
- 系と目的を決める：何を対象にし、何を求めるかを明確にする
- 記号を定義する：各記号の意味、単位、正方向を決める
- 仮定を書く：定常、完全混合、密度一定、等温などを確認する
- 原理式を立てる：物質収支、エネルギー収支、運動量収支などから始める
- 未知数と式の数を数える：解くための情報が足りているか確認する
- 文字式で整理する：可能なら数値代入の前に求める量について解く
- 単位をそろえて代入する：特に時間、圧力、温度、体積に注意する
- 結果を検算する：単位、符号、桁、極限、物理的な増減を確認する
よくある間違い
1. 微分する変数を確認しない
\(dC/dt\) は時間に対する濃度変化、\(dC/dx\) は位置に対する濃度変化です。分母の変数が変われば意味も単位も変わります。
2. 定積分の上下限を付けない
具体的な初期状態から最終状態を求めるときは、積分定数を曖昧にせず、初期条件を使った定積分で書くと符号ミスを減らせます。
3. \(\ln a+\ln b=\ln(a+b)\) とする
正しくは、
\[
\ln a+\ln b=\ln(ab)
\]
です。対数の足し算は、中身の掛け算に対応します。
4. 指数・対数の中に単位を残す
\(e^{-kt}\) の \(kt\)、\(\ln(C/C_0)\) の \(C/C_0\) は無次元です。指数や対数の中身を確認すると、単位の不整合に気付きやすくなります。
5. 微分方程式の解に初期条件を使わない
一般解には任意定数が含まれます。装置の初期濃度や初期温度を使って、対象とする系の具体的な解へ絞り込みます。
6. 計算結果の符号だけを消す
濃度変化が負なら減少、熱流束が負なら定義した正方向と逆向きなど、負号には物理的な意味があります。都合よく絶対値へ置き換えず、符号規約を確認します。
7. 数値だけを見て桁を確認しない
流量5 L/minと5 m³/sでは6万倍の差があります。単位変換後の桁が現実的か、装置規模と照らして確認します。
どこまで数学を学べばよいか
化学工学の基礎を学ぶ段階では、まず次の内容を使える状態を目指します。
- 一次方程式と連立方程式
- 比例、反比例、べき乗
- 指数関数と自然対数
- グラフの傾きと面積
- 基本的な微分と積分
- 変数分離形と一次線形微分方程式
- 偏微分の記号と「何を一定にするか」の意味
- 差分、台形則、回帰の基本
- 単位・次元と有効数字
行列、偏微分方程式、数値最適化、統計的推定なども重要ですが、最初からすべてを学ぶ必要はありません。化学工学の問題に触れ、必要になった数学を往復しながら学ぶ方が定着しやすくなります。
確認問題
3問のミニテストで、ここまでの内容を確認しましょう。選択すると、その場で正誤と解説が表示されます。
QUICK CHECK
1 / 3化学工学の基礎・初級
このミニテストの回答結果は保存されません。学習記録を残す場合は化工演習室をご利用ください。
まとめ
- 数式は、物理量どうしの関係を簡潔に表した文章である
- グラフの傾きは変化速度、面積は累積量を表す
- 指数関数は、変化速度が現在量に比例する現象に現れる
- 対数を使うと、指数関係の整理や直線化ができる
- 微分は局所的な変化速度、積分は小さな寄与の積み上げである
- 保存則を非定常系へ適用すると微分方程式が現れる
- 微分方程式には、具体的な解を決める初期条件が必要である
- 指数・対数の中身は無次元にし、すべての式で単位を確認する
- 最後に符号、桁、極限、物理的な増減を検算する
数学は、化学工学を難しくするためのものではありません。複雑な現象から重要な関係を取り出し、装置の挙動を予測するための共通言語です。式を計算手順としてだけでなく、現象の説明として読むことから始めましょう。
参考資料
- MIT OpenCourseWare, Numerical Methods Applied to Chemical Engineering
- MIT OpenCourseWare, First Order Linear ODEs
- MIT OpenCourseWare, Contaminant Fate Modeling
- OpenStax, Exponential and Logarithmic Functions
- OpenStax, Derivatives as Rates of Change
- OpenStax, Basics of Differential Equations
- NIST Guide to the SI, Rules for Expressing Values of Quantities
コメント