後向きステップ流れ
理論と物理
概要
先生、後向きステップ流れってどんな問題なんですか？
流路の途中に段差（ステップ）があり、そこで流れが剥離して再付着する現象だ。剥離・再付着流れの最も基本的なベンチマーク問題として、CFDコードの検証に古くから使われている。Armaly et al. (1983) の実験データが有名だね。
段差の後ろで渦ができるイメージですか？
そうだ。ステップ下流に再循環領域（剥離泡）ができる。その長さ $x_r$ をステップ高さ $h$ で割った再付着長さ $x_r/h$ がReynolds数の関数になる。これが最も重要な検証指標だ。
支配方程式
支配方程式はNavier-Stokesですよね？
非圧縮性Navier-Stokes方程式と連続の式だ。
Reynolds数はステップ高さ $h$ と入口の平均流速 $U$ で定義する。
再付着長さとReの関係はどうなるんですか？
層流域（$Re < 400$ 程度）では再付着長さはReにほぼ比例する。$x_r/h \approx 0.06 \times Re$ という近似が知られている。Armalyらの実験では拡大比 $ER = (H+h)/H = 1.94$ で、Re=100 のとき $x_r/h \approx 5$、Re=400 で $x_r/h \approx 14$ だ。
Reが大きくなると再付着点が遠くなるんですね。
ただし $Re > 400$ 前後で3次元効果が顕著になり、2D計算だけでは実験と合わなくなる。これが有名な「Armaly問題の2D-3D遷移問題」だ。
流れの構造
再循環領域以外にも構造があるんですか？
Reが大きくなると、ステップ対面の上壁側にも二次的な剥離泡が現れる。さらにステップ角部にも小さな渦が形成される。流れ場の全体構造はReに強く依存するんだ。
| Re範囲 | 流れの特徴 |
|---|---|
| Re < 200 | 下壁の主再循環のみ |
| 200 < Re < 400 | 上壁にも二次剥離泡出現 |
| Re > 400 | 3次元不安定性、スパン方向変動 |
| Re > 1000 | 乱流遷移、非定常渦放出 |
2次元計算で十分なのはRe=400くらいまでなんですね。勉強になります。
再付着点の位置が「乱流モデル選びの踏み絵」になる理由
後向きステップ流れで「剥離後に流れがステップ下流の壁に再付着するまでの距離（再付着長さ）」は、乱流モデルの性能を比較する定番指標です。実験値はステップ高さの約6〜8倍の位置に再付着します。ところが標準k-εモデルで解くと再付着が遅れて9〜11倍あたりになりがちで、SSTモデルを使うと精度が上がる。「同じ問題を解いているのに答えが違う」この差が、乱流モデル選択の実務的な根拠です。後向きステップはシンプルな形状にもかかわらず「剥離・再循環・再付着」が全部入っているので、新しい乱流モデルを開発したとき最初に試す試験台になっています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値手法
後向きステップ流れの数値解法で気をつけることは何ですか？
この問題では圧力と速度の連成が重要だ。非圧縮流れの定番アルゴリズムであるSIMPLE系やProjection法を使う。
圧力-速度連成
SIMPLE法ってどんな仕組みですか？
Semi-Implicit Method for Pressure-Linked Equations の略だ。Patankar & Spalding (1972) が提案した。手順は以下の通り。
1. 仮の速度場を運動量方程式から求める
2. 圧力補正方程式（Poisson型）を解く
3. 速度と圧力を補正する
4. 収束するまで繰り返す
派生手法としてSIMPLEC、PISO（非定常向き）がある。定常計算ならSIMPLE/SIMPLEC、非定常ならPISOが標準だ。
空間離散化
対流項のスキームは何がいいですか？
後向きステップでは再循環があるので、上流差分（1次精度）だと数値拡散で再付着長さが長くなりすぎる。最低でも2次精度が必要だ。
| スキーム | 精度 | 安定性 | 再付着長さへの影響 |
|---|---|---|---|
| 1次風上 | 1次 | 高 | 過大評価（数値拡散大） |
| 2次風上 | 2次 | 中 | 適切 |
| QUICK | 3次 | やや低 | 適切 |
| 中心差分 | 2次 | 低 | 振動の恐れ |
メッシュ設計
メッシュはどこを細かくすべきですか？
壁面第一層の $y^+$ は、層流なら不要だが、乱流計算なら $y^+ < 1$（壁関数なしの場合）が望ましい。構造格子なら拡大比1.1〜1.2で壁面から離れるように分布させる。
出口を遠く取るのが大事なんですね。30hですか、結構長い。
短すぎると出口境界条件の影響で再付着長さが変わってしまう。検証時は出口位置の感度も確認すべきだ。
ステップ直後メッシュの「密度格差」問題
後向きステップの数値計算で見逃されやすいのが「ステップコーナー直後のメッシュ遷移」問題です。剥離点（ステップ端）は速度勾配が最大になる場所なので密なメッシュが必要ですが、そこから再付着点に向かって急激に粗くすると数値拡散が増えて再循環領域が過大・過小評価されます。経験則では「ステップ高さの1/10以下の要素サイズをステップ端周辺に使い、下流方向に膨張率1.1以下で広げる」が安全策です。また、2Dモデルで再付着長さを実験と合わせたとしても、3D効果（コーナー流れ）を無視しているためにその設定が3Dモデルで通用しないことも多い。「2D検証→3D本番」では設定を引き継がないほうが無難です。
風上差分（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. 形状定義
- 拡大比 $ER = (H+h)/H$ を決める。Armaly問題なら $ER = 1.94$（$h = 0.0049$ m、$H = 0.0052$ m）
- 入口助走区間: $10h$ 以上（完全発達流を入口条件にするなら不要）
- 出口区間: $30h$ 以上
2. 境界条件
- 入口: 完全発達放物線速度分布 $u(y) = U_{max}[1-(2y/H)^2]$ または一様流速
- 出口: 自然流出条件（$\partial u/\partial x = 0$, $p = 0$）
- 壁面: no-slip条件
3. メッシュ収束確認
| メッシュ水準 | セル数（2D） | $x_r/h$（Re=200） |
|---|---|---|
| 粗 | 5,000 | 10.8 |
| 中 | 20,000 | 10.2 |
| 細 | 80,000 | 10.05 |
| 参照値 | - | 10.0 |
粗いメッシュだと再付着長さが8%もずれるんですか。
再付着長さは壁面せん断応力がゼロになる位置で決まるから、メッシュ感度が高い。少なくとも3水準のメッシュでRichardson外挿するのが望ましい。
検証のポイント
計算結果の妥当性はどうやって確認しますか？
壁面摩擦係数の符号が変わる場所を見ればいいんですね。
その通り。$C_f > 0$ なら順流、$C_f < 0$ なら逆流。この符号変化点が再付着点だ。
燃焼器の「ステップ」——意図的な再循環が炎を安定させる
後向きステップ流れは「避けるべき流れ現象」として語られることが多いですが、ガスタービンや工業炉の燃焼器設計では逆に「意図的にステップで再循環域を作る」ことがあります。再循環域には高温ガスが滞留するため、新たに入ってくる燃料・空気の着火源になる——つまり「炎の種火」として機能します。この再循環域がないと冷たい流れが炎を吹き消してしまいます（フレームアウト）。航空エンジンの燃焼器設計でCFDが必須になった背景の一つは、このステップ後の再循環形状を最適化するためです。「問題現象」が「有用な現象」になるトポロジーの好例です。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツールでの実装
主要なCFDソフトでの設定方法を教えてください。
後向きステップはほぼ全てのCFDソフトにチュートリアルがある。ツール別のポイントを整理しよう。
Ansys Fluent
OpenFOAM
simpleFoam（定常）、icoFoam/pimpleFoam（非定常）blockMesh で構造格子生成。grading で壁面近傍を細分化divSchemes に Gauss linearUpwind grad(U) を設定fixedValue で放物線分布。codedFixedValue が便利OpenFOAMだと blockMesh でステップ形状を作るんですね。
ブロックを2つ組み合わせるだけで作れる。上流側の高さが $H$、下流側が $H+h$ になるようにブロックを定義する。
STAR-CCM+
ツール間比較
| 項目 | Fluent | OpenFOAM | STAR-CCM+ |
|---|---|---|---|
| 設定の容易さ | GUI操作 | テキスト辞書 | GUI+Java |
| メッシュ柔軟性 | 非構造/構造 | blockMesh/snappyHex | ポリヘドラル/Directed |
| 再付着長さ精度 | 同等 | 同等 | 同等 |
| 並列性能 | 良好 | 良好 | 良好 |
どのソフトでも精度は同等なんですね。メッシュと設定次第ということですか。
その通り。ソルバーの違いよりもメッシュ品質と対流スキームの選択が結果を左右する。
ベンダー比較の「後向きステップ検証」が信頼の証明書になる理由
CFDソルバーのベンダーがパンフレットやウェブサイトに後向きステップの検証結果を掲載しているのをよく見かけます。理由は「この1問題が乱流計算能力の総合評価として業界で広く認知されているから」です。剥離・再循環・再付着が全部含まれ、実験データも豊富（Driver & Seegmiller 1985、Armaly et al. 1983など）で比較しやすい。ソルバー選定のときに「後向きステップのベンチマーク結果を見せてください」と頼むのは、実はかなり有効な評価手法です。どのモデルで何%誤差か、をきちんと開示しているベンダーは技術力と誠実さの両方で信頼できる傾向があります。
選定で最も重要な3つの問い
- 「何を解くか」：後向きステップ流れに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端トピック
後向きステップ流れの研究って、まだ進んでいるんですか？
基本問題だからこそ、新手法の検証に使われ続けている。
乱流モデルの検証
高Re（Re > 5000）ではRANSモデルの性能差が顕著に出る。Driver & Seegmiller (1985) の実験（$Re_h = 37,400$）が乱流検証の標準だ。
| 乱流モデル | 再付着長さ予測 | 実験値との乖離 |
|---|---|---|
| $k$-$\varepsilon$ 標準 | 過小（$x_r/h \approx 5.5$） | -15% |
| $k$-$\omega$ SST | 良好（$x_r/h \approx 6.2$） | -4% |
| RSM（Reynolds応力） | 良好（$x_r/h \approx 6.4$） | -1% |
| LES（Smagorinsky） | 良好（$x_r/h \approx 6.5$） | 0% |
| 実験値 | $x_r/h \approx 6.5$ | - |
標準 $k$-$\varepsilon$ はかなりずれるんですね。
逆圧力勾配が強い剥離流れでは $k$-$\varepsilon$ 標準モデルは渦粘性を過大評価しがちだ。SST以上のモデルを使うべきだよ。
LES/DNSによる研究
Le, Moin & Kim (1997) がDNSで $Re_h = 5,100$ を計算した。再循環領域内の乱流統計量（レイノルズ応力テンソルの各成分）が得られ、RANSモデル改良の基礎データになっている。
最近ではWall-Modeled LES（WMLES）で工学的Re数（$Re_h > 10^5$）の計算も可能になりつつある。
3次元効果とスパン方向不安定性
先ほどRe=400あたりで3次元効果が出るという話がありましたが、詳しく教えてください。
Barkley et al. (2002) の線形安定性解析によると、$Re \approx 748$ で3次元モード（スパン方向波長 $\lambda \approx 6.9h$）が最初に不安定になる。これは定常分岐で、Hopf分岐ではない。つまり3次元だが定常のパターンが出現するんだ。
非定常じゃなくて定常の3次元パターンなんですか。面白い。
この分岐構造の理解は、2D RANS計算の適用限界を見極めるうえで非常に重要だ。
Coffee Break よもやま話
LESで見える「再循環域の非定常ゆらぎ」が設計を変える
RANSで後向きステップを解くと再付着長さは一定値に収束しますが、LESや実験では再付着点が上下流にランダムに揺れ続けています。この非定常ゆらぎは単なるノイズではなく、物体への動的な圧力変動荷重を生み出します。建物の壁面に流れが再付着する都市風工学の問題では、この変動荷重が疲労破壊の原因になるため、RANSでの定常解析では設計が甘くなるリスクがある。「ステップ後ろに何かを取り付けたい」という設計案件では、LESで時系列データを取ることが強く推奨される理由はここにあります。
Le, Moin & Kim (1997) がDNSで $Re_h = 5,100$ を計算した。再循環領域内の乱流統計量（レイノルズ応力テンソルの各成分）が得られ、RANSモデル改良の基礎データになっている。
最近ではWall-Modeled LES（WMLES）で工学的Re数（$Re_h > 10^5$）の計算も可能になりつつある。
先ほどRe=400あたりで3次元効果が出るという話がありましたが、詳しく教えてください。
Barkley et al. (2002) の線形安定性解析によると、$Re \approx 748$ で3次元モード（スパン方向波長 $\lambda \approx 6.9h$）が最初に不安定になる。これは定常分岐で、Hopf分岐ではない。つまり3次元だが定常のパターンが出現するんだ。
非定常じゃなくて定常の3次元パターンなんですか。面白い。
この分岐構造の理解は、2D RANS計算の適用限界を見極めるうえで非常に重要だ。
LESで見える「再循環域の非定常ゆらぎ」が設計を変える
RANSで後向きステップを解くと再付着長さは一定値に収束しますが、LESや実験では再付着点が上下流にランダムに揺れ続けています。この非定常ゆらぎは単なるノイズではなく、物体への動的な圧力変動荷重を生み出します。建物の壁面に流れが再付着する都市風工学の問題では、この変動荷重が疲労破壊の原因になるため、RANSでの定常解析では設計が甘くなるリスクがある。「ステップ後ろに何かを取り付けたい」という設計案件では、LESで時系列データを取ることが強く推奨される理由はここにあります。
トラブルシューティング
トラブルシューティング
後向きステップの解析でありがちなトラブルを教えてください。
典型的な問題と対策を整理しよう。
再付着長さが参照値と合わない
原因と対策:
- メッシュ不足: 壁面近傍の解像度不足。壁直交方向のセル数を増やす
- 1次精度スキーム: 数値拡散で再循環が潰される。2次精度以上に変更
- 出口が近すぎる: 下流 $30h$ 以上確保。感度チェックで $40h$ と比較
- 入口条件の不一致: 完全発達流を仮定しているか、一様流か。文献と条件を揃える
入口条件って結果にそんなに影響しますか？
大きく影響する。一様流入口だと助走区間が必要で、入口距離が短いと発達しきらない速度分布でステップに到達してしまう。放物線プロファイルを直接与えるか、十分な助走区間（$20h$ 以上）を設けよう。
計算が収束しない
定常計算が発散する場合:
- Reが高すぎて実際は非定常→非定常計算に切り替える
- 緩和係数が大きすぎる→圧力 0.2〜0.3、速度 0.5〜0.7 に下げる
- メッシュのアスペクト比が極端→壁面近傍のグレーディングを見直す
非定常計算のCFL条件:
2D計算と実験の不一致
Re=600くらいで実験と全然合わないんですが...
それは3次元効果だ。Re > 400 では3D計算が必要。2D計算は再付着長さを過大評価する傾向がある。スパン方向に周期境界条件を設けた3D計算に切り替えよう。スパン幅は $8h$ 以上が推奨だ。
チェックリスト
| 確認項目 | 対処法 |
|---|---|
| 再付着長さのメッシュ収束 | 3水準以上でRichardson外挿 |
| 出口距離の十分性 | $30h$ 以上、感度確認 |
| 入口速度分布の妥当性 | 放物線 or 十分な助走区間 |
| 対流スキームの精度 | 2次以上 |
| 適切なRe範囲 | 2D: Re < 400、3D: それ以上 |
まずはメッシュ収束と出口距離を確認すれば、大半の問題は解決できそうですね。
その通り。基本に忠実にやれば、後向きステップは非常に教育的で信頼できるベンチマークだ。
「再付着点が実験より遠い」は必ずしも誤りではない
後向きステップで「再付着長さが実験のXH（ステップ高さX倍）より長く出た」場合、すぐ「計算が間違い」と思いがちですが、まず確認すべきことがあります。実験のRe数・入口境界層状態（層流か乱流か、その厚さ）が正確に再現されているかです。入口が完全発達乱流プロファイルでなく一様流だと、せん断層の発達が遅れて再付着点が実験より遠くなることがあります。つまり「計算は正しく、実験再現の境界条件が不正確」というケース。入口から一定の助走距離を設けてから乱流を発達させる「ドライバーセクション」の設定は、後向きステップ解析のこだわりポイントの一つです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——後向きステップ流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告