自由表面流れ
理論と物理
概要
先生、自由表面流れって何ですか？VOF法とは違うんですか？
自由表面流れ特有の物理は何ですか？
表面張力と接触角が小スケールで重要になる。大スケールでは重力と慣性の比であるフルード数 $Fr$ が支配パラメータだ。
表面張力のCSFモデル
表面張力はどうモデル化するんですか？
Brackbill et al.（1992）のCSF（Continuum Surface Force）モデルが最も広く使われている。界面の曲率 $\kappa$ から体積力として計算する。
曲率の計算精度が問題になると聞きました。
CSFの最大の課題が寄生流（parasitic currents / spurious currents）だ。界面の離散化誤差に起因する非物理的な速度場が発生する。毛細管数 $Ca = \mu U / \sigma$ が小さい流れ（表面張力支配）で特に深刻だ。
界面追跡法の比較
| 手法 | 界面表現 | 質量保存 | 界面シャープさ | 計算コスト |
|---|---|---|---|---|
| VOF法 | 体積分率 $\alpha$ | 厳密 | スキーム依存 | 低〜中 |
| Level Set法 | 符号付き距離関数 $\phi$ | 非保存（再初期化必要） | シャープ | 中 |
| CLSVOF | VOF + Level Set | 良好 | シャープ | 中〜高 |
| Phase-Field | 秩序パラメータ | 保存 | 拡散的 | 高 |
| SPH | 粒子 | 保存 | 粒子分解能依存 | 高 |
どの手法を選べばいいんですか？
工業的にはVOF法がコストと精度のバランスで最も広く使われている。界面の曲率精度が重要な場合（マイクロ流体等）はCLSVOFやPhase-Field法が有利だ。大変形・飛散を伴う場合はSPH法も選択肢になる。
自由表面の力学——Bernoulliと毛細管が出会う場所
自由表面（気液界面）は流体力学のなかで最も豊かな物理現象が凝縮した場所です。重力・表面張力・粘性・慣性がすべて競合し、その比を表す無次元数（Bond数、Weber数、Capillary数、Froude数）の組み合わせで支配する物理が決まります。Bond数 Bo = ρgL^2/σ が1を超えると重力が表面張力を圧倒し、逆に小さければ毛細管現象が支配します。マイクロ流体デバイスではBo << 1で自由表面が球面に収束し、大型船舶の波ではBo >> 1で表面張力は無視できます。この20桁以上に及ぶBo数のスケール多様性が、自由表面CFDをあらゆる産業で必要とされる技術にしています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
自由表面流れの数値的なポイントを教えてください。
界面のCourant数管理とメッシュ設計が最重要だ。VOF法では界面Courant数 $Co_\alpha < 0.25$（最大0.5）を守る必要がある。
Open Channel Flow
河川や水路の流れでは、FluentやCFXにOpen Channel Flowの専用設定がある。入口・出口で水位を直接指定でき、水面の初期化が容易だ。
Open Channel Flowの設定ポイントは何ですか？
波の生成と吸収
海洋工学の波浪シミュレーションでは、造波と消波の境界条件が必要だ。
| 手法 | 概要 | ツール対応 |
|---|---|---|
| Stokes波理論 | 入口で波の速度・水位を規定 | Fluent, STAR-CCM+ |
| 5次Stokes波 | 高次の非線形波 | STAR-CCM+ |
| 消波ゾーン | ダンピングで波を吸収 | 全ツール（UDF/Field Function） |
| 数値波動水槽 | 造波板 + 消波ゾーン | OpenFOAM (waves2Foam) |
OpenFOAMのwaves2Foamは有名ですよね。
Jacobsen et al.（DTU, 2012）が開発したライブラリで、波の生成・吸収のための境界条件と弛緩ゾーンが実装されている。海洋構造物の波力解析で広く使われている。
寄生流対策
寄生流を抑えるにはどうすればいいですか？
| 対策 | 効果 | 備考 |
|---|---|---|
| メッシュ細分化 | 高 | 界面付近のセルを十分小さく |
| Height Function法 | 非常に高 | 構造格子向け、曲率精度向上 |
| Sharp Surface Force | 高 | STAR-CCM+で利用可能 |
| CLSVOF | 高 | Level Setで曲率計算を改善 |
| cAlphaの調整 | 中 | OpenFOAMの界面圧縮パラメータ |
Level Set法とVOF法の融合——CLSVOFの登場
Level Set法は界面の幾何学的情報（法線・曲率）を精度よく計算できる一方、質量保存が不完全という弱点を持ちます。VOFは質量保存に優れますが界面形状の記述精度が低い。2000年代にSussmanらが提案したCoupled Level-Set and VOF（CLSVOF）法は両者の長所を組み合わせ、Level Set関数で曲率を計算しながらVOFで質量を保存するハイブリッド手法です。航空機燃料タンクのスロッシング解析への適用では、CLSVOFが単純VOFより25%少ない格子で同等精度を実現したベンチマーク結果があります。
風上差分（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. 形状定義: タンク形状（矩形 or 円筒）、液面高さ
2. メッシュ: 界面付近を細かく（液面高さ方向に20セル以上）
3. 初期条件: setFields（OpenFOAM）or Patch（Fluent）で液面を設定
4. 物性値: 液体と気体の密度・粘度、表面張力
5. 外力: タンクの加振をDynamic Mesh or Moving Reference Frameで設定
6. VOF設定: Explicit + Geo-Reconstruct（Fluent）、MULES（OpenFOAM）
7. タイムステップ: 界面Co < 0.25
AMR（適合格子細分化）の活用
界面付近だけメッシュを細かくしたいんですが。
AMR（Adaptive Mesh Refinement）が非常に有効だ。体積分率の勾配でリファイン基準を設定し、界面追従でメッシュを動的に細分化する。
| ツール | AMR対応 | リファイン基準 |
|---|---|---|
| Fluent | Gradient of Volume Fraction | 自動 |
| STAR-CCM+ | Table-based AMR | カスタマイズ可能 |
| OpenFOAM | dynamicRefineFvMesh | Field-based |
| Basilisk | Octree AMR | ネイティブ |
Basiliskは前にVOF法の記事でも出てきましたね。
Stephane Popinet（Sorbonne大学）が開発した多相流専用OSSソルバーで、octreeベースのAMRが組み込まれている。学術分野では非常に評価が高く、自由表面流れの高精度シミュレーションで多くの実績がある。
接触角の設定
壁面と液面の接触角はどう設定するんですか？
静的接触角はVOF法の壁面境界条件として指定する。動的接触角モデル（前進・後退接触角の差、速度依存）も利用可能だ。マイクロ流体デバイスやインクジェットなど、濡れ性が重要な系では接触角の設定が結果を大きく左右する。
船舶抵抗・耐航性解析——造船所を変えたCFDの自由表面
船舶の造波抵抗（Wave-Making Resistance）は全抵抗の30〜50%を占め、その予測精度が設計の優劣を決めます。2000年代以降、Reynolds-Averaged NS方程式+VOF法のCFDが大型商船設計の必須ツールとなり、造波係数の予測精度が5%以内を達成するようになりました。現代の大型LNGタンカー設計では、船首形状（球状船首バルブ）の最適化をCFDで行うことで燃費を2〜4%改善でき、年間燃料費に換算すると数億円の差に相当します。船級協会が認める解析条件と検証方法が規格化されており、CFDが「試験の補完」から「設計の主役」に変わりつつあります。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
自由表面流れに対応しているツールを比較してください。
| ツール | 界面追跡法 | AMR | Open Channel | 特徴 |
|---|---|---|---|---|
| Ansys Fluent | VOF (HRIC, Geo-Reconstruct, CLSVOF) | 対応 | 対応 | 最も多くのVOFスキーム |
| STAR-CCM+ | VOF (HRIC) | 対応 | 対応 | 海洋工学に強い |
| Ansys CFX | Homogeneous Free Surface | 限定的 | 対応 | 結合型ソルバーの安定性 |
| OpenFOAM (interFoam) | VOF (MULES, isoAdvector) | 対応 | waves2Foam | 完全OSS |
| Flow-3D | TruVOF | 構造格子AMR | 対応 | 鋳造・自由表面特化 |
| SPHFlow / DualSPHysics | SPH法 | N/A（粒子法） | 対応 | 大変形・砕波 |
Flow-3Dって何ですか？
Flow Science社が開発した自由表面流れ特化のCFDソルバーだ。TruVOF法は構造格子上で高精度な界面追跡を行い、充填・凝固・波浪など自由表面を伴う問題で長い実績がある。鋳造シミュレーションでは業界標準の一つだ。
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 船舶造波抵抗 | STAR-CCM+, Fluent | 海洋工学の実績 |
| タンクスロッシング | Fluent, Flow-3D | 動メッシュ + VOF |
| ダムブレイク・洪水 | OpenFOAM, SPH | 大変形流れ |
| 鋳造充填 | Flow-3D, Fluent | 充填 + 凝固の連成 |
| マイクロ流体 | COMSOL, OpenFOAM | 表面張力の精度 |
| 海洋構造物波力 | OpenFOAM (waves2Foam), STAR-CCM+ | 波の生成・吸収 |
SPH法はどういう場面で使うんですか？
津波、砕波、グリーンウォーター（甲板上への海水打ち込み）など、界面が激しく変形・飛散する問題でメッシュベースの手法が苦手とする領域で有効だ。DualSPHysicsはオープンソースのGPU対応SPHコードで、港湾工学の研究で広く使われている。
OpenFOAM vs STAR-CCM+——自由表面解析の実務選択
自由表面流れのCFDツール選択は、コストと機能のトレードオフで決まります。OpenFOAMのinterFoamは無料で世界最大のユーザーベースを持ち、新しいアルゴリズム（isoAdvector等）がいち早く実装されます。STAR-CCM+（Siemens）は一体化されたGUIとCAD連携・自動メッシュ生成が強みで、造船・オフショア分野でのターンアラウンドタイム短縮が評価されています。ANSYS Fluentは汎用性と連成解析（FSI、熱伝達）の統合が強みです。年間ライセンスコストはStar-CCM+が400万円〜、Fluentが500万円〜に対し、OpenFOAMはゼロですが社内教育・サポートコストが無視できません。
選定で最も重要な3つの問い
- 「何を解くか」：自由表面流れに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
自由表面流れの最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
isoAdvector法
Roenby et al.（DTU, 2016）が提案した幾何学的VOF法で、非構造格子上でも界面をiso-surfaceとして正確に再構成できる。OpenFOAM v1806以降に実装されており、従来のMULES法と比較して界面のシャープさと質量保存の両方が向上する。
ハイブリッドEuler-Lagrange法
連続的な液体領域はVOF法で解き、飛散した液滴はDPM（Lagrangian粒子）に自動変換するハイブリッド手法が発展している。Fluent 2020以降のVOF-to-DPM転換機能が代表的だ。
噴霧の一次分裂をVOFで捕まえて、二次分裂はDPMで追うということですね。
まさにその通り。逆にDPMの液滴が壁面に衝突して液膜を形成するDPM-to-VOF転換もある。マルチスケールの多相流を効率的に扱える。
GPU計算
大規模な自由表面流れのCFDをGPUで加速する研究が進んでいる。Fluent 2024以降ではGPUネイティブソルバーがVOF法をサポートし始めている。SPH法はアルゴリズムの特性上GPUとの親和性が高く、DualSPHysicsでは1GPUで数百万粒子のリアルタイム計算が可能だ。
機械学習による波浪予測
AIの活用もあるんですか？
Coffee Break よもやま話
砕波の直接シミュレーション——乱流と自由表面の境界
砕波（Breaking Wave）は自由表面CFDの最難関問題の一つで、波の先端が乱流化してエアーエントレインメント（空気巻き込み）を起こす現象です。RANSモデルでは砕波の乱流エネルギーを過大評価する傾向があり、LESが必要ですが、自由表面+乱流+気泡生成の三者を同時に解くのは計算負荷が桁違いです。東京大学の研究グループが並列VOF-LESで波高1 mクラスの砕波を直接計算し、エアーエントレインメント率が実験値の±15%以内に収まることを示しました。海面CO2交換速度は砕波で2〜3倍に増加するため、地球科学的にも砕波CFDの高精度化は重要課題です。
連続的な液体領域はVOF法で解き、飛散した液滴はDPM（Lagrangian粒子）に自動変換するハイブリッド手法が発展している。Fluent 2020以降のVOF-to-DPM転換機能が代表的だ。
噴霧の一次分裂をVOFで捕まえて、二次分裂はDPMで追うということですね。
まさにその通り。逆にDPMの液滴が壁面に衝突して液膜を形成するDPM-to-VOF転換もある。マルチスケールの多相流を効率的に扱える。
大規模な自由表面流れのCFDをGPUで加速する研究が進んでいる。Fluent 2024以降ではGPUネイティブソルバーがVOF法をサポートし始めている。SPH法はアルゴリズムの特性上GPUとの親和性が高く、DualSPHysicsでは1GPUで数百万粒子のリアルタイム計算が可能だ。
機械学習による波浪予測
AIの活用もあるんですか？
砕波の直接シミュレーション——乱流と自由表面の境界
砕波（Breaking Wave）は自由表面CFDの最難関問題の一つで、波の先端が乱流化してエアーエントレインメント（空気巻き込み）を起こす現象です。RANSモデルでは砕波の乱流エネルギーを過大評価する傾向があり、LESが必要ですが、自由表面+乱流+気泡生成の三者を同時に解くのは計算負荷が桁違いです。東京大学の研究グループが並列VOF-LESで波高1 mクラスの砕波を直接計算し、エアーエントレインメント率が実験値の±15%以内に収まることを示しました。海面CO2交換速度は砕波で2〜3倍に増加するため、地球科学的にも砕波CFDの高精度化は重要課題です。
トラブルシューティング
トラブルシューティング
自由表面流れでよくあるトラブルを教えてください。
順番に見ていこう。
1. 界面がぼやける
症状: 液面が拡散して水と空気の境界がわからなくなる。
対策:
- 界面再構成スキームをHRICからGeo-Reconstruct（Fluent）やisoAdvector（OpenFOAM）に変更
- メッシュを界面付近で細分化（AMR推奨）
- OpenFOAMでは cAlpha を1.0〜1.5に調整
- 界面Courant数を0.25以下に管理
2. 寄生流が発生する
静止しているはずの液滴の周りに渦ができるんですが…
対策:
- メッシュを細かくする（液滴直径に対して最低20セル）
- Height Function法の曲率計算を使う（構造格子）
- FluentではImplicit Body Forceを有効化
- STAR-CCM+ではSharp Interface Methodを検討
3. 液面が初期化できない
対策:
- FluentではPatch操作またはOpen Channel Flow Initializationを使用
- OpenFOAMではsetFieldsユーティリティを使用
- 初期の圧力場を静水圧で初期化する（$p = \rho g h$）
4. Open Channel Flowで水面が振動する
対策:
- 入口と出口の水位を整合させる
- ダンピングゾーンを出口近傍に設定
- タイムステップを十分に小さくする
- Implicit Body Forceを有効にする
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Operating Pressure設定に注意。VOF法では0（ゲージ圧基準）推奨 |
| STAR-CCM+ | VOFのSharpeningファクター設定が界面品質に影響 |
| OpenFOAM | interFoamのalphaEqnSubCyclesとnAlphaSubCyclesの調整 |
| Flow-3D | 構造格子のみ。複雑形状はFAVOR法で表現するため形状解像度に注意 |
波が消える——数値拡散が自由表面を殺す
自由表面CFDで最も多い品質問題は「波が伝わるうちに振幅が減衰する」数値拡散です。一次精度の風上差分（First-Order Upwind）を移流スキームに使うと、波長の数倍進んだだけで波高が半減することがあります。対策は高次精度スキーム（MUSCL、CICSAM等）の使用と、波長に対して十分なメッシュ解像度（波長方向に80格子以上が目安）の確保です。数値波消波器（Numerical Beach）を出口付近に設置しないと反射波が計算域に戻って結果を汚染します。「波が消えた」と感じたらまず時間刻みとメッシュ解像度の系統的なコンバージェンス確認を行うことが鉄則です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——自由表面流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告