コンピュータシミュレーションによる葉巻の燃焼過程の技術パラメータ入力と燃焼時間予測の方法

タバコ科学における数値シミュレーション研究の中で、葉巻の燃焼過程のシミュレーションは、極めて挑戦的なマルチフィジックス連成問題と常にみなされてきた。比較的均一で構造が単純な紙巻きタバコとは異なり、葉巻は充填煙葉（Filler）、結合煙葉（Binder）、包覆煙葉（Wrapper）から構成される複雑な多孔質媒体システムである。この構造上の不均質性は、酸素拡散経路、熱伝導効率、燃焼フロント形状の高度な複雑化をもたらす。研究者にとって、これらの複雑な物理現象をいかに精確な数学モデルへ変換し、計算によって得られた燃焼時間が実験室での実測値と高度に一致することを確保するかが、シミュレーション精度を評価する最も中心的な基準である。

多孔質媒体フレームワークにおける物理モデリングの基礎

本章では、葉巻を多孔質媒体連続体として扱うモデリングの基礎を解説し、運動量方程式、浸透率の層別設定、有効熱伝導率の補正を扱います。

葉巻の燃焼を正確に記述するには、まず葉巻を「固体ブロック」とみなす単純化された思考を捨て、多孔質媒体連続体モデル（Porous Media Continuum Model）を採用しなければならない。このフレームワークでは、葉巻は固体粒子（煙草繊維、灰分）と孔隙（空気または熱分解ガスで満たされている）からなる複合体系として扱われる。

運動量方程式の処理においては、ナビエ・ストークス（Navier-Stokes）方程式を単純に使用することはできず、流体が煙葉層を通過する際の圧力降下を記述するためにダルシー・フォルヒハイマー（Darcy-Forchheimer）項を導入しなければならない。式は以下のとおりである：

$$\nabla p = -\frac{\mu}{K} \mathbf{u} - \frac{1}{2} \rho C_F |\mathbf{u}| \mathbf{u}$$ ここで、$K$ は浸透率（Permeability）、$C_F$ は慣性抵抗係数である。実際のモデリングでは、浸透率の設定は燃焼形態（円錐形燃焼かトンネル式燃焼か）を決定する鍵となる。葉巻の部位ごとの特性に応じて、$K$ を層状に設定する必要がある：充填層の浸透率は通常高く（約 $10^{-10}$ から $10^{-9} \text{ m}^2$）、包覆層は葉繊維の緻密さにより、その浸透率はしばしば一桁低い。熱伝導過程も、有効熱伝導率（Effective Thermal Conductivity、$k_{eff}$）の補正を考慮する必要がある。煙葉層内部には大量の気相対流と放射が存在するため、単一の固体熱伝導率ではエネルギー伝達速度を著しく過小評価してしまう。我々は通常、体積加重法で推定する： $$k_{eff} = \varepsilon k_g + (1-\varepsilon) k_s$$

ここで $\varepsilon$ は空隙率、$k_g$ と $k_s$ はそれぞれ気相と固相の熱伝導率である。燃焼フロント付近では、温度の激しい変動により、この等価処理は局所的な空隙率の時間発展とリアルタイムに連成されなければならない。

10⁻¹⁰–10⁻⁹ m²

充填層の浸透率の典型的な値の範囲

12%

初期モデルにおける葉巻全体の統一含水率設定

18.5% → 3.2%

水分蒸発修正の導入前後における予測誤差の変化

2.26 MJ/kg

水の気化潜熱の参考値

パラメータ入力の不確実性：水分と密度の連成の罠

本章では、含水率の空間的不均一性と熱シンク効果がモデルの歪みの重要な罠となることを明らかにします。

長年の実践的なモデリングの中で、モデルの失敗を最も引き起こしやすいのは複雑な動力学方程式ではなく、基礎的な物理パラメータ入力の「過度な単純化」であることに気づいた。

2023年にある高級キューバ葉巻の燃焼特性を対象とした研究で、我々は厄介な難題に直面した：シミュレーションで得られた燃焼時間は実験による実測値より約18.5%速く、燃焼フロントの温度分布は不自然な「スパイク」状を示し、実際に観測された緩やかな燃焼波面と全く一致しなかった。

深い遡及分析の結果、問題の根源は煙草の含水率（Moisture Content、MC）の空間分布の不均一性と、それがもたらす熱シンク効果を見落としていたことにあると判明した。初期モデルでは、計算を簡単にするために、葉巻全体の含水率を一定の12%に統一して設定していた。しかし、実際の生産における葉巻では、包覆層が環境湿度の影響をより直接受けるため、含水率は外側から内側への勾配分布を示すことが多い。さらに重要なのは、水分の蒸発が巨大な吸熱過程であることだ（水の気化潜熱は約 $2.26 \text{ MJ/kg}$）。

含水率が高い場合、熱分解反応が放出するエネルギーはまず水分蒸発に使われ、燃焼フロントに明らかな「温度緩衝帯」を形成する。改良後のモデルでは、アレニウス（Arrhenius）速度式に基づく過渡的水分蒸発モデルを導入し、確率場法（Stochastic Field Method）を採用して含水率の空間的不均一性をシミュレーションした。水分蒸発と熱分解反応を分離し、局所的な蒸発熱項を導入することで、シミュレーションの予測誤差は18.5%から3.2%以内へと急速に低下した。この経験から学んだ教訓は、煙草のようなバイオマス材料を扱う際、水分は決して単なる質量分率の変数ではなく、動的な熱力学の調節器であるということだ。

燃焼動力学と移動境界の数学的処理

本章では、熱分解と酸化が交互に進行する燃焼動力学、および移動境界問題の数学的処理と適応メッシュ技術を解説します。

葉巻燃焼の核心は、熱分解（Pyrolysis）と酸化（Oxidation）の交互進行にある。完全な燃焼サイクルは通常、以下の論理チェーンに従う：固体煙葉が加熱される $\rightarrow$ 水分が蒸発する $\rightarrow$ 有機物が熱分解して揮発分を生成する $\rightarrow$ 揮発分が気相に入り酸化反応を起こす $\rightarrow$ 残留チャー（Char）が酸素拡散制御下で表面酸化を起こす。

動力学パラメータの入力時には、多段反応モデルを採用しなければならない。単一の活性化エネルギーでは、複雑な煙草の化学組成を記述できない。我々は通常、多成分熱分解モデルを採用し、煙草をセルロース類、リグニン類、残留灰分に分け、それぞれに異なるアレニウスパラメータ（$A, E_a$）を付与する。

燃焼時間を予測する最も中心的なアルゴリズムは、本質的に「移動境界問題」（Moving Boundary Problem）を解くことであり、数学的にはステファン問題（Stefan Problem）に類似している。燃焼フロントの位置 $\mathbf{s}(t)$ は時間の経過とともに移動し、その移動速度 $\dot{r}$ は局所的な質量輸送とエネルギー収支に依存する：

$$\rho_s (1-\varepsilon) \frac{d\mathbf{s}}{dt} = \dot{m}_{pyro} + \dot{m}_{ox}$$ ここで $\dot{m}_{pyro}$ は熱分解速度、$\dot{m}_{ox}$ はチャー酸化速度である。実際のCFD求解プロセスでは、この極めて薄く激しい反応層を捕捉するために、適応メッシュ細分化技術（Adaptive Mesh Refinement、AMR）を採用しなければならない。メッシュ解像度が燃焼フロントの勾配をカバーするのに不十分な場合、数値拡散が反応速度を人為的に「ならして」しまい、予測される燃焼時間が長くなる。このような高非線形問題を扱う際、時間ステップの選定もまた芸術的である——クーラン数（Courant Number）制約を満たす必要があると同時に、熱拡散の時間スケールも考慮しなければ、計算は燃焼爆発の瞬間に発散しやすい。 ## 予測アルゴリズムの統合と工学的実装最終的に正確な燃焼時間を得るために、我々は微視的反応から巨視的時程に至る積分モデルを構築した。まず、過渡CFDシミュレーションを通じて、各微要素位置における異なる時刻の温度場 $T(x,t)$、酸素濃度場 $C_{O_2}(x,t)$、局所含水率場 $M(x,t)$ を取得する。次に、これらの場変数に基づき、局所燃焼速度 $\dot{r}(x,t)$ を用いて、葉巻の幾何学的長さ $L$ にわたって時空間積分を実行する。総燃焼時間 $t_b$ を予測する式は、以下のように簡略化して表せる： $$t_b = \int_0^L \frac{1}{\bar{\dot{r}}(x)} dx$$

ここで $\bar{\dot{r}}(x)$ は軸方向に沿った平均線速度である。しかし工学的実践では、葉巻の燃焼過程に「点火・消火・再点火」のパルス的特性が存在するため（特に吸引プロセスをシミュレーションする場合）、この積分プロセスは時間次元で累積されなければならない。各吸引サイクル（Puff）内の燃焼進行と、2回の吸引の間の陰燃（Smoldering）段階を精密に時系列で接続する必要がある。

「陰燃段階」の予測を解決する際、もう一つの重要なポイントは境界条件の設定である。吸引停止後、葉巻は低酸素環境下での緩慢な酸化状態にある。このとき環境境界条件を単純に常温空気と設定すると、シミュレーションされた陰燃温度が高くなりすぎ、予測される燃焼時間が短くなる。環境対流熱伝達係数 $h$ を考慮した動的境界を導入し、葉巻表面の放射放熱も考慮する必要がある。

まとめと展望

本章では、過渡CFDの場変数から燃焼時間の積分予測に至る工学的実装経路を紹介します。

上記の複雑なパラメータ入力とアルゴリズム構築を通じて、我々は葉巻の燃焼時間の精確な予測を実現できるだけでなく、さらに重要なことに、シミュレーションによって燃焼過程の物理的詳細を観察できる——例えば、包覆葉の浸透率が低すぎるために生じる「燃焼不均一」現象や、水分勾配による「燃焼停滞」現象などである。

現在のシミュレーション技術はより高次元の方向へ発展している。マイクロCT（Micro-CT）スキャン技術の普及に伴い、煙葉の微視的孔隙構造データを極めて精確に取得し、それを直接CFDモデルに導入して、「連続体近似」から「実在の幾何学構造」への飛躍を実現できる。同時に、機械学習アルゴリズムと従来のCFDを連成し、既知の実験データを用いてニューラルネットワークを訓練して複雑な反応速度パラメータを予測することが、煙草の燃焼動力学パラメータの取得難と不確実性の大きさを解決する次の技術的突破口になるかもしれない。