液膜モデル
理論と物理
概要
先生、液膜モデルって何ですか？
壁面上に形成される薄い液体の膜の流動・蒸発・飛散を計算するモデルだ。自動車のウィンドシールド上の雨水、航空機翼の着氷、エンジン内壁の燃料液膜、塗装コーティングなど、壁面を流れる液膜の挙動を予測する。
VOF法で壁面の液膜を解くのとは違うんですか？
液膜の厚さが数十μm〜数mmと非常に薄いので、VOF法で直接解像するにはメッシュを極端に細かくする必要があり、現実的でない。液膜モデルは壁面上の2Dシェル方程式で液膜を記述するため、3Dメッシュに依存せず効率的に計算できる。
支配方程式
液膜の方程式を教えてください。
液膜の質量保存（膜厚の変化）を記述する式は次の通りだ。
$h$ は液膜厚さ、$\bar{\mathbf{u}}_f$ は膜厚方向に平均した液膜速度、$\nabla_s$ は壁面に沿った勾配演算子だ。右辺のソース項はそれぞれ液滴衝突（impingement）、蒸発、飛散（splash）による質量変化を表す。
液膜速度はどうやって求めるんですか？
薄膜近似（lubrication theory）を使う。液膜内部の速度分布は壁面でno-slip、液膜表面でせん断力（気流からのせん断応力 $\tau_g$）のバランスから放物線分布になる。膜厚方向に平均すると、
第1項は圧力勾配と重力の壁面接線成分による駆動、第2項は気流せん断による駆動だ。液膜のエネルギー方程式も同様に薄膜近似で解き、蒸発率を計算する。
液滴-壁面相互作用
液滴が壁面に衝突したときの挙動はどうモデル化するんですか？
衝突レジームはWeber数と壁面温度で決まる。
| レジーム | 条件 | 挙動 |
|---|---|---|
| Stick | $We < We_{cr,low}$ | 壁面に付着 |
| Rebound | 高温壁面 | 弾性的に反射 |
| Spread | 中程度の$We$ | 広がって液膜形成 |
| Splash | $We > We_{cr,high}$ | 飛散して二次液滴生成 |
Stanton-Rutland モデルや Bai-Gosman モデルが代表的で、Fluent や STAR-CCM+ に実装されている。
液膜の薄さが生む複雑さ——µmスケールの支配方程式
壁面液膜（Wall Film）は厚さが1〜1000 µmという極めて薄い液体層で、航空機アイシング・エンジン壁面冷却・胃の粘液層まで多様な場所に現れます。薄膜近似（Thin Film Approximation）では厚さ方向の速度分布を放物線と仮定することで、3次元ナビエ-ストークス方程式を2次元の薄膜方程式に還元できます。液膜表面に生じるMarangoni対流（温度・濃度勾配による表面張力差駆動流れ）は塗装工程での塗膜均一性や熱交換器の液膜不均一化に直結し、実用上の重要性が高い現象です。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
液膜モデルの数値的な解法を教えてください。
液膜は壁面のサーフェスメッシュ上で解かれる。3DのCFDメッシュとは独立に、壁面境界フェイスの2D接続情報を使って膜厚・速度・温度の輸送方程式を解く。
気相CFD ↔ 液膜モデルの連成は以下の情報交換で行う。
| 気相 → 液膜 | 液膜 → 気相 |
|---|---|
| 壁面せん断応力 $\tau_g$ | 蒸発による質量ソース |
| 壁面近傍の温度・濃度 | 蒸発による熱ソース |
| DPM液滴の壁面衝突 | 液膜からの飛散液滴 |
| 壁面圧力分布 | 液膜表面のラフネス効果 |
DPMの液滴が壁面に衝突して液膜になって、またちぎれて液滴になるんですね。
ツール別の実装
| ツール | 液膜モデル名 | 主な機能 |
|---|---|---|
| Ansys Fluent | Eulerian Wall Film | 液膜流動、蒸発、splash、DPM連成 |
| STAR-CCM+ | Thin Film Model | 液膜流動、熱伝達、蒸発、飛散 |
| OpenFOAM | regionFaModel | 有限面積法、基本的な液膜流動 |
| Ansys CFX | Wall Film（限定的） | 基本的な液膜追跡 |
FluentとSTAR-CCM+が充実しているんですね。
液膜モデルは自動車・航空宇宙産業のニーズが大きいため、この2つのツールが最も成熟した実装を持っている。OpenFOAMの regionFaModel は有限面積法（Finite Area Method）ベースで、研究目的のカスタマイズに適している。
薄膜数値解法——壁面曲率と重力の統一処理
壁面液膜のCFD実装では、複雑形状の壁面上での液膜流れに対してShell要素アプローチが有効です。壁面の法線方向を積分した積算方程式（Integral Method）により、液膜厚さhと平均速度の輸送方程式が導出されます。ANSYS Fluentのwall film modelは重力・圧力勾配・せん断応力・蒸発・凝縮をすべてソース項として統一的に扱い、エンジン壁面油膜の挙動予測に広く使われています。ただし厚い液膜（h > 1 mm程度）や乱流液膜では薄膜近似が破綻し、3次元VOFへの切り替えが必要になることがあります。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
液膜モデルを使った解析の手順を教えてください。
自動車のウォーターマネジメント（フロントガラスの雨水挙動）を例にとろう。
1. 外部流れ解析: 車両周りの空気流れを定常で計算
2. 雨滴投入: DPMで雨滴を車両前方から投入
3. 液膜形成: 壁面衝突モデルで液膜を形成
4. 液膜流動: 気流せん断と重力で液膜が流動
5. 飛散: 高せん断領域で液膜がちぎれて二次液滴が発生
6. 後処理: 液膜厚さ分布、視界妨害領域の評価
メッシュの注意点
液膜モデルで特に気をつけるべきメッシュのポイントはありますか？
液膜は壁面メッシュ上で解かれるので、壁面の面メッシュ品質が重要だ。
| 項目 | 推奨 | 理由 |
|---|---|---|
| 壁面メッシュサイズ | 1〜5 mm | 液膜の流動パターンを解像 |
| 壁面メッシュ品質 | スキューネス < 0.7 | 2D輸送方程式の精度 |
| 壁面の曲率対応 | 十分なメッシュ密度 | 液膜の流れ方向の正確性 |
| 3Dメッシュ壁面第1層 | 通常の壁関数推奨 | せん断応力の正確な伝達 |
液膜が薄いところと厚いところでメッシュを変える必要はありますか？
液膜の厚さ自体はメッシュに依存しない（壁面上のスカラー変数として計算される）ので、液膜厚さでメッシュを変える必要はない。ただし液膜が集中する領域（窪み、エッジ等）は壁面メッシュを細かくしたほうが良い。
検証手法
液膜モデルの結果を検証するにはどうすればいいですか？
代表的な検証実験を挙げよう。
| 実験 | 条件 | 計測量 |
|---|---|---|
| 傾斜平板上の液膜流 | 重力駆動、Nusselt解と比較 | 膜厚、流量 |
| 液膜に対する気流せん断 | 風洞実験 | 液膜速度、飛散開始風速 |
| 単一液滴壁面衝突 | 高速カメラ撮影 | スプレッド径、splash液滴径 |
| 車両走行時の水膜 | 実車風洞 | 液膜流動パターン |
自動車エンジンの壁面油膜——燃費と排気規制を繋ぐCFD
直噴エンジンのシリンダー壁に衝突した燃料液滴が油膜を形成し、この油膜の蒸発遅延が未燃炭化水素（HC）排出量を増加させます。Euro 7規制（2025年以降）の粒子数規制値は厳しく、壁面油膜蒸発のCFD精度が実際の認証試験通過に直結しています。BMWやトヨタが公開したベンチマーク事例では、壁温が20℃から80℃に変わると油膜蒸発速度が3倍になることがCFDで確認されており、暖機過程の精密シミュレーションがコールドスタートHC低減策の鍵となっています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
液膜モデルに対応しているツールを比較してください。
| ツール | モデル名 | DPM連成 | 蒸発 | 用途実績 |
|---|---|---|---|---|
| Ansys Fluent | Eulerian Wall Film | 完全対応 | 対応 | 自動車、航空機着氷 |
| STAR-CCM+ | Thin Film Model | 完全対応 | 対応 | 自動車、エンジン |
| OpenFOAM | regionFaModel | 基本対応 | 限定的 | 学術研究 |
| CONVERGE | Film Model | 対応 | 対応 | エンジン内壁液膜 |
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 自動車ウォーターマネジメント | Fluent, STAR-CCM+ | DPM + Film + 外部流れの統合 |
| 航空機着氷（Icing） | Fluent (FENSAP-ICE) | 専用の着氷モジュール |
| エンジンシリンダー壁面液膜 | CONVERGE, STAR-CCM+ | 動メッシュ + 液膜 |
| コーティング・塗装 | Fluent | 液膜の膜厚分布予測 |
| 学術・モデル開発 | OpenFOAM | Finite Area Methodのカスタマイズ |
FENSAP-ICEって何ですか？
Numerica（現Ansys傘下）が開発した航空機着氷シミュレーション専用モジュールだ。液膜の流動・凍結・氷形状成長を一貫して計算できる。FAA（米国連邦航空局）の認証プロセスでも使われている。Fluent 2020以降で統合されている。
着氷シミュレーションでは液膜が凍るところまで計算するんですね。
その通り。過冷却水滴が翼面に衝突して液膜を形成し、気流と熱交換しながら凍結する。Messinger モデルやShallcross モデルで凍結率を計算する。液膜モデルと相変化モデルの連成が鍵だ。
FENSAP-ICE——航空機アイシングCFDの業界標準
航空機着氷（アイシング）解析では液膜（水膜・氷膜）のCFDが飛行安全認証の根幹をなします。BoeingやAirbusが採用しているFENSAP-ICE（現ANSYS Fluent Ice）は、外流・液滴衝突・液膜流れ・着氷の4段連成解析を専用に設計したツールです。FAA（米国連邦航空局）の規格ではCFDによるアイシング予測の精度要件が明記されており、FENSAPのvalidationパッケージがそのまま型式証明の提出資料に使われています。新型航空機の開発では着氷試験飛行（実機icing test）の前に必ずCFDで設計検証することが業界慣行となっています。
選定で最も重要な3つの問い
- 「何を解くか」：液膜モデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
液膜モデルの最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
液膜-気流連成のLES
壁面液膜と気流の相互作用をLESで高精度に解く研究が進んでいる。気流のせん断応力の時間変動が液膜の波立ちやリガメント形成に与える影響を直接捕捉する。
VOFとの切り替え
液膜が厚くなったらVOFに切り替えたりできますか？
実はこれが活発な研究テーマだ。薄い液膜は液膜モデル（Thin Film）で効率的に計算し、液膜が厚くなったり大きな液塊が形成されたらVOF法に自動遷移する手法が研究されている。STAR-CCM+のFluid Film + VOFハイブリッドモデルがこの方向性だ。
液膜の不安定性と飛散
液膜がちぎれて液滴が飛散するメカニズムのモデリングが重要課題だ。Kelvin-Helmholtz不安定性やRayleigh-Taylor不安定性に基づく飛散モデルの研究が進んでいる。
噴霧の一次分裂と似た物理ですね。
まさにそう。液膜のエッジや高せん断部からのリガメント形成と分裂は、液柱の一次分裂と同じ物理だ。Film strippingモデルの改良が航空機着氷や自動車防水設計で重要になっている。
接触角・濡れ性の影響
液膜の移動接触線（液膜先端のエッジ）の挙動は接触角に依存する。前進・後退接触角のヒステリシスモデルや、表面粗さの影響を取り入れた濡れ性モデルの研究が進んでいる。自動車のソイリング（泥汚れ）予測に不可欠だ。
液膜のMarangoni対流——表面張力勾配が生む薄膜流れ
温度や濃度の分布が液膜表面に存在すると、表面張力の差（Marangoni効果）によって表面に沿った流れが生じます。この「Marangoni対流」は食器洗いで洗剤を落とすとコーヒーが素早く広がる現象と同じ原理です。半導体ウェハの液膜コーティング（スピンコーティング）では、Marangoni対流が膜厚均一性を±1%以下に制御するために積極的に利用されています。CFDでMarangoni対流を正確に再現するには表面張力の温度・濃度微分を精密に入力する必要があり、このデータが文献にない新規流体では信頼性が著しく下がります。
トラブルシューティング
トラブルシューティング
液膜モデルでよくあるトラブルを教えてください。
順番に見ていこう。
1. 液膜が壁面に留まらず消える
症状: 液滴が壁面に衝突しても液膜が形成されない。
対策:
- Wall Film Modelが該当壁面で有効になっていることを確認
- DPMの壁面境界条件が「wall-film」に設定されていることを確認（「reflect」や「escape」ではなく）
- 液滴のWeber数が低すぎてstickレジームに入っていないか確認
2. 液膜厚さが非物理的に大きい
液膜が分厚くなりすぎるんですが…
対策:
- 液膜の蒸発モデルが有効になっているか確認（蒸発なしだと溜まり続ける）
- 液膜の排出境界条件（Film Edge）が適切に設定されているか確認
- せん断応力が液膜に正しく伝達されているか確認（壁関数の解像度）
3. 液膜が不安定に振動する
対策:
- タイムステップを小さくする
- 壁面メッシュ品質を改善する（スキューネスの高い面を修正）
- 液膜厚さの初期値を小さな正の値（$10^{-6}$ m）で初期化
4. DPM-Film連成で液滴が壁面を通過する
対策:
- DPMのMaximum Number of Stepsを増やす
- DPMの積分タイムステップが気相タイムステップより十分小さいか確認
- 壁面メッシュの法線方向が正しいか確認
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Wall Film zoneの定義を忘れやすい。Named Selectionで明示的に設定 |
| STAR-CCM+ | Thin Film Physics Modelを壁面regionに追加する必要がある |
| OpenFOAM | regionFaModelの壁面パッチ指定とメッシュ接続に注意 |
| CONVERGE | Film modelのactivation conditionとmesh embeddingの整合性確認 |
液膜が分離する——フィルムドライアウトの予測失敗
壁面液膜CFDで最も実害の大きい予測ミスが「フィルムドライアウト」の見逃しです。液膜が薄くなってゼロになる点（ドライパッチ）では局所的に熱伝達率が急落し、金属表面温度が数百℃急上昇することがあります。CFDでは液膜厚さの最小値制限（min film thickness）パラメータが設定されていることが多く、これが大きすぎるとドライアウトが起きないように計算され続ける「偽の安全」が生じます。航空機アイシング解析でフィルムドライアウトを見逃すと、翼面に氷が局所的に堆積して揚力が低下するという危険な結果につながります。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——液膜モデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告