混合対流
理論と物理
混合対流とは
先生、混合対流って強制対流と自然対流が同時に起きている状態ですよね？
そのとおり。外部からの流れ（ファンやポンプ）と浮力による流れが同程度の強さで共存する場合を混合対流（mixed convection）と呼ぶ。たとえば垂直管内の上昇流で管壁を加熱する場合、ポンプによる強制流れと浮力による上昇流が重畳する。
どちらが支配的かはどうやって判断するんですか？
Richardson数 $Ri$ で判断する。
$Ri \ll 1$ なら強制対流支配、$Ri \gg 1$ なら自然対流支配、$Ri \sim O(1)$ なら混合対流だ。実務的には $0.1 < Ri < 10$ の範囲が混合対流領域と考えてよい。
助勢流と対向流
浮力と強制流れの方向関係は重要ですか？
非常に重要だ。加熱垂直管内の上昇流では浮力が流れを助ける方向に働く（助勢流 / aiding flow）。逆に下降流では浮力と流れが逆向きになる（対向流 / opposing flow）。
助勢流ではNu数が純粋な強制対流より増加する。対向流では流れの減速・逆流・再層流化が起こりうるため、Nu数が複雑に変化する。特に対向流の再層流化現象（laminarization）はCFDで正確に予測するのが難しいことで知られているよ。
水平管の場合はどうですか？
水平管では浮力が二次流れ（縦渦）を生む。管断面の上部に高温流体が集まり、下部に低温流体が溜まる。この非対称な温度分布を正確に予測するには3D計算が必須で、2D軸対称近似は使えない。
Richardson数が「どっちが強い？」を決める
混合対流で鍵になるRichardson数（Ri = Gr/Re²）は、「浮力が強制流れの何倍強いか」を示す比率だ。Ri ≪ 1なら強制対流が支配的でシンプルに解けるが、Ri ≈ 1付近では浮力と慣性力が拮抗し、流れが複雑に絡み合う。電子機器の筐体内部がまさにこの厄介なゾーンに入りやすい——ファンで気流を作りつつ、基板の温度差で浮力流も起きる。現場では「まず強制対流だけで解いてみて、壁面温度が高くなりすぎたらBoussinesq項を追加」という段階的アプローチが多い。Ri ≫ 1の領域は自然対流が独壇場で、ファンを止めた緊急冷却時の計算がこれにあたる。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
Boussinesq近似の適用範囲
混合対流のCFDでは密度の扱いが重要だと聞きました。
浮力項を正しく扱うために、密度の温度依存性をどうモデル化するかが鍵だ。最も単純なのがBoussinesq近似で、密度を
と線形近似し、浮力項にのみこの変動を反映する。運動方程式の他の項では密度を一定とする。Ansys FluentではOperating ConditionsでGravityを設定し、MaterialのDensityをBoussinesqに設定する。
Boussinesq近似はいつ使えないんですか？
温度差が大きくて $\beta \Delta T > 0.1$〜$0.2$ 程度になると精度が悪化する。空気の場合、$\Delta T > 30$度C程度で注意が必要だ。この場合は理想気体やpolynomial密度を使って密度の非線形温度依存性を直接扱うべきだ。FluentのIncompressible Ideal Gas設定が便利だよ。
乱流モデルの注意点
混合対流特有の乱流モデル選択の注意点はありますか？
浮力による乱流生成/減衰効果が重要になる。k-ε系モデルでは浮力生成項 $G_b = -g_i \frac{\mu_t}{\rho Pr_t} \frac{\partial \rho}{\partial x_i}$ が追加される。FluentではViscous ModelのOptionsで「Full Buoyancy Effects」をONにすることを強く推奨する。
また、対向流での再層流化現象を予測するにはTransition SSTモデルが有効だ。標準的な乱流モデルでは層流化を予測できず、Nu数を過大評価する。
メッシュの要件は強制対流と異なりますか？
壁面垂直方向のメッシュ要件は同等（$y^+ \approx 1$ 推奨）だが、浮力による二次流れを解像するために管断面方向のメッシュも十分に細かくする必要がある。水平管の場合、断面方向に少なくとも40〜60分割は必要だ。粗すぎると二次流れが解像できずNu数の非対称性を過小評価する。
Boussinesq近似——「密度は一定、でも浮力は計算する」という矛盾
Boussinesq近似は一見するとおかしな仮定だ。「密度ρは一定とする、ただし浮力項だけは温度依存の密度変化を使う」——これは矛盾に見える。でも実はこの近似は物理的に正しい。温度差ΔTが小さい場合（ΔT ＜ 20〜30℃程度が目安）、密度変化は浮力力として流体を動かすには十分だが、慣性力や連続式への影響は無視できる程度に小さい。この分離によって非線形性が大幅に緩和され、計算が安定する。一方、高温炉内や太陽光を受ける建物外壁のような大きなΔT環境では、密度の全変化を考慮する「全密度モデル（ideal gas等）」が必要になる。モデル選択の目安はβΔT ＜ 0.1を確認することだ。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
建築環境での混合対流
建築のHVAC設計でも混合対流は重要ですか？
非常に重要だ。オフィスの空調では天井からの給気（強制対流）と室内の発熱体（PC、人体、照明）による浮力が共存する。天井吹出し口からの冷気が十分な運動量を持たないと、暖気層（warm layer）が室上部に停滞して温度成層が発生する。
それをCFDで予測するわけですね。
そのとおり。Fluentの場合、定常RANSではRealizable k-ε + Enhanced Wall Treatmentが室内環境CFDの標準選択肢だ。OpenFOAMではbuoyantSimpleFoamにkOmegaSSTを組合せる。人体の発熱は80〜120W程度のheat sourceとしてモデル化し、PMV-PPDによる快適性評価まで含めるのが実務的なワークフローだよ。
電子機器筐体の熱設計
電子機器の筐体内部も混合対流ですか？
ファン付きの筐体では強制対流が主体だけど、ファンなし（自然空冷）やファン故障時は浮力駆動になる。設計上は両方のケースを評価する必要があり、これはまさに混合対流問題だ。
実務では3Dの筐体モデルを作成し、各部品をvolume heat sourceとしてモデル化する。基板はorthotropic（異方性）熱伝導体として面内/面外の熱伝導率を分けて入力する。STAR-CCM+やFluentのCHT機能で固体-流体を同時に解くのが標準的だよ。
筐体内の自然対流って収束しにくいイメージがあります。
鋭い指摘だ。密閉筐体内の自然対流は定常解が存在しない場合がある（Ra数が高いと非定常振動流になる）。定常計算で残差が振動する場合は、非定常計算に切り替えて時間平均を取るべきだ。FluentならTransient設定でadaptive time steppingを使うのが実用的だよ。
検証手法
混合対流CFDの検証にはどんなベンチマークがありますか？
垂直管内の混合対流ではJacksonらの実験データが標準参照だ。矩形キャビティの混合対流ではHaidari et al.のデータがある。建築分野ではIEA Annex 20のベンチマークケース（室内環境CFD）が広く使われている。
データセンター冷却設計——混合対流CFDで消費電力30%削減
現代のハイパースケールデータセンターでは、強制冷却（床下給気）と自然対流（ホットアイル/コールドアイル）が複雑に混合する典型的な混合対流環境だ。GoogleとFacebookが公開した事例では、CFDによる気流最適化でPUE（電力使用効率）を1.5から1.15に改善し、冷却電力を約30%削減した実績がある。解析上の難点はサーバーラックの多孔質体モデル化——実際のティル開口率とCFD設定の整合が肝で、開口率±10%の誤差がホットスポット温度予測に±8℃の差を生む。現場計測との比較検証が不可欠だ。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
Ansys Fluentでの浮力設定
Fluentで混合対流を解くときの設定手順を教えてください。
重要なポイントを順に。(1) Operating Conditions: Operating Densityを適切に設定。Boussinesq近似なら参照温度での密度。(2) Gravity: 重力ベクトルを正しく設定（デフォルトは0）。(3) Material: DensityをBoussinesqまたはIncompressible Ideal Gasに設定。(4) Viscous Model: Full Buoyancy EffectsをON。
Operating Densityって何ですか？
Boussinesq近似では圧力の計算から静水圧成分を除去するためにOperating Densityを使う。参照温度での密度に設定するのが標準。これを設定しないと圧力場に大きな静水圧成分が残り、数値精度が悪化する。Incompressible Ideal Gasを使う場合はOperating Densityは不要だ。
STAR-CCM+での浮力設定
STAR-CCM+の場合は？
Physics ModelsでGravity ModelをONにし、重力ベクトルを指定する。密度はBoussinesq ModelまたはIdeal Gas（非圧縮性近似付き）を選択。乱流モデルでBuoyancy ProductionをONにするのを忘れないこと。Reference ValuesでReference Temperatureを適切に設定するのも重要だ。
OpenFOAMでの浮力設定
OpenFOAMの場合はどうですか？
buoyantSimpleFoam（定常）またはbuoyantPimpleFoam（非定常）を使う。constant/gファイルで重力加速度を設定。密度はconstant/thermophysicalPropertiesのequationOfStateでBoussinesq（Boussinesq近似）またはperfectGas（理想気体）を指定する。Boussinesqの場合、betaパラメータ（体膨張係数）の入力が必要だ。
圧力のsolverはSIMPLEでいいですか？
浮力がある場合はSIMPLE系でも解けるが、収束が遅いことがある。Pseudo transient（FluentのPseudo Transient設定）やPIMPLE法（OpenFOAMのbuoyantPimpleFoamの定常モード活用）が安定する場合が多い。圧力のunder-relaxationは0.3程度から始めて、収束を見ながら調整しよう。
浮力駆動流の収束判定で注意することは？
残差だけでなく、モニタリングポイントの温度と速度が定常に達していることを確認する。浮力流は残差が10のマイナス3乗程度で振動しても、積分量は収束している場合がある。逆に残差が下がっていても積分量が変動しているなら収束していない。
混合対流ソルバーの「浮力フラグ」を忘れずに
商用ソルバーで混合対流を解く際、意外と見落とされるのが「浮力項の有効化」だ。Fluent、STAR-CCM+、CFX のいずれも、エネルギー方程式とNavier-Stokes方程式を連立させるだけでは不十分で、温度差による密度変化（Boussinesq近似または全密度法）を明示的にオンにする必要がある。これを忘れると、熱は伝わっても流体が動かないという「熱だけ解けている状態」になる。また浮力方向（重力ベクトル）の設定ミスも多く、g = (0, -9.81, 0) のつもりが (0, 9.81, 0) になっていると浮力が逆向きになり、暖かい空気が下に溜まる非物理的な結果が出る。ソルバー選定前のチェックリストに「浮力モデルの設定手順の複雑さ」を加えておくことを勧める。
選定で最も重要な3つの問い
- 「何を解くか」：混合対流に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
再層流化（Laminarization）
対向流の再層流化って具体的にどういう現象ですか？
加熱された垂直管内の上昇流（助勢流）で、浮力による加速が壁面近傍の速度プロファイルを一様化させ、乱流エネルギーの生成が抑制されて流れが局所的に層流化する現象だ。Nu数が急激に低下するため、過剰な壁面温度上昇の原因になりうる。
再層流化のパラメータとして $Bo^ = Gr^ / (Re^{3.425} Pr^{0.8})$ が使われ、$Bo^* > 6 \times 10^{-7}$ 程度で層流化が始まるとされている。
RANSで再層流化を予測できますか？
通常のk-εやk-ωでは予測できない。Transition SSTモデルを使えばある程度予測可能だが、再層流化の開始位置が実験と10〜30%ずれることがある。正確な予測にはDNS（Direct Numerical Simulation）やWall-Resolved LESが必要だ。
LESによる混合対流解析
LESで混合対流を解くメリットは何ですか？
混合対流では浮力プルームの非定常的な揺動や乱流との相互作用が重要で、RANSではこれを時間平均化してしまう。LESでは浮力による大規模渦構造や温度変動を直接解像できるので、非定常熱荷重の評価やThermal fatigue寿命予測に有用だ。
計算コストの目安は？
垂直管内の混合対流LESでは、$Re_{\tau}$ に基づく壁面解像のセル数が圧倒的に多くなる。$Re = 10000$ 程度の管流で数千万セル、$Re = 50000$ なら数億セルが必要。HPCクラスターでも1ケース数日〜数週間かかる。WMLESで壁面モデルを使えばセル数を1/10程度に削減できるよ。
原子力分野での応用
原子力で混合対流が重要なのはなぜですか？
超臨界流体の混合対流——CO₂冷却炉の浮力逆転現象
超臨界CO₂（sCO₂）は次世代原子炉や太陽熱発電の冷媒として注目されているが、臨界点付近（31℃、7.38MPa）では密度・比熱・粘性が急変し、通常の混合対流理論が破綻する。特に垂直加熱管の上昇流では「浮力補助から浮力抑制への遷移」が発生し、局所熱伝達係数が急低下する「伝熱劣化」現象が起きる。この現象はBuoyancy Parameter Bo* = Gr/Re³⁻⁴ で整理できるが、CFDでは圧力依存の物性変化をReal-Gas EOSとLow-Mach数ソルバで扱う必要があり、商用ツールの標準設定では再現できないケースがある。
トラブルシューティング
計算が発散する
混合対流の計算を始めたら、数十反復で発散してしまいます。
浮力流は初期条件との乖離が大きいと速度場が急激に変動して発散しやすい。以下の手順を試そう。(1) まずgravityをOFFにして等温流れ場を収束させる。(2) gravityをONにして温度差を段階的に増加させる（FluentならPatch機能で温度を設定）。(3) Pseudo Transientまたは非定常計算で安定化する。
under-relaxation factorはどう調整すべきですか？
密度のunder-relaxationを追加して0.8程度に設定する（Fluentでは通常1.0がデフォルトでリストに表示されない。Solution > Methods > Under-Relaxation FactorsでDensityを追加する）。Body Forceのunder-relaxationも0.8程度に下げるとよい。
温度場が非物理的な対称性を持つ
水平管の混合対流なのに、温度分布が上下対称になってしまいます。
浮力効果がモデルに反映されていない可能性がある。確認ポイント: (1) gravityベクトルの向きが正しいか。(2) DensityがConstantではなくBoussinesqまたはIdeal Gasになっているか。(3) Operating DensityまたはReference Densityが適切に設定されているか。(4) Full Buoyancy EffectsがONか。
全部確認したけど対称です。
非常に低いRi数（$Ri < 0.01$）なら浮力効果がほぼ無視できるので対称は正常。$Ri > 0.1$ で対称なら、初期条件が完全対称になっていて、浮力による摂動が数値的に発達しない可能性がある。わずかな非対称初期摂動（温度場に0.1度程度のランダムノイズ）を加えるか、非定常計算に切り替えるとよい。
OpenFOAMのbuoyantSimpleFoamで残差が振動する
buoyantSimpleFoamで残差が周期的に振動して収束しません。
浮力流では定常解が存在しない場合がある。まずRa数を確認しよう。$Ra > 10^8$ 程度（密閉キャビティの場合）では流れが本質的に非定常になる。buoyantPimpleFoamに切り替えて非定常計算し、時間平均を取るのが正しいアプローチだ。
非定常計算の時間刻みの目安は？
浮力流の特性速度 $u_{buoy} = \sqrt{g \beta \Delta T L}$ とセルサイズ $\Delta x$ からCFL数が1以下になるように設定する。adaptive time steppingが使えるならMaxCo 0.5程度で自動調整させるのがよい。OpenFOAMのcontrolDictでadjustTimeStepをyesに設定し、maxCoを指定すればよい。
混合対流の収束が遅い——その本当の理由
混合対流の定常解析で「残差が10⁻³から先に全く下がらない」という経験をしたエンジニアは多いはずだ。これは単なる設定ミスではなく、物理的な理由がある。Ri ≈ 1付近の混合対流は流れが本質的に「不安定」な状態にあり、複数の安定解が共存できる。つまり定常解が存在しない、あるいは複数存在するという状況になっている。この場合、定常ソルバーを使い続けても収束しない——解そのものが振動している可能性が高い。対策として非定常計算（URANS）に切り替え、十分な時間積分後に時間平均を取る手法が有効だ。また緩和係数を0.1〜0.3程度に下げると振動が抑制されることもある。「定常で解けない混合対流は非定常に切り替えろ」が経験則だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——混合対流の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告