発表のポイント
- 層流から乱流へ切り替わる直前の平板境界層について高解像度の直接数値計算を行い、壁のごく近くに、境界層厚さの0.35〜0.55倍という微小な波長をもつ秩序だった流れ構造（コヒーレント構造）が現れることを示した。この波長は、同じレイノルズ数の遷移過程でこれまで報告されてきた構造より1〜2桁小さい。
- スパン方向速度変動が約45度に傾いたX字型のパターンをもつことを、二点空間相関により統計的に確認した。このスケールと向きは、格子解像度、マッハ数（0.1〜0.3）、上流擾乱の与え方を変えても保たれる。線形安定性解析により、古典的なトルミーン・シュリヒティング波の不安定性では説明できないことも示した。
- 壁面近傍（境界層厚さの11.4%以内）の変動を数値的に抑えると、下流の乱流渦の発達が弱まった。遷移過程がこの薄い層の変動に敏感であることを示す結果である。ただし、この構造が遷移の引き金であると断定できる段階ではなく、生成機構の解明は今後の課題である。
研究の背景
一見無秩序に見える乱れた流れのなかに、秩序だった運動が潜んでいることが実験的に示された（Hama et al. 1957; Kline et al. 1967）。当時これは革命的な発見とされ、乱流に対する見方を大きく変えた。
線形安定性理論に加え、線形化ナビエ・ストークス作用素の非正規性に由来する「非モーダル過渡成長」の枠組みが整った。安定性理論では不安定とならない条件でも遷移が起きる理由が理解されるようになり、バイパス遷移の記述が進んだ。
直接数値計算（DNS）の進展により、ストリーク、斜め波、局所的な乱流塊であるパフ・スラグ、斜めに並ぶ乱流バンドなど、遷移過程に現れる構造が次々と同定された。いずれも境界層厚さと同程度か、それを大きく上回るスケールをもつ構造である。
当研究室では、微小な凹凸を面全体に不規則に分布させた表面（DMR）による遷移遅延効果をDNSと風洞試験で調べてきた。その過程で、遷移が始まる手前の壁のごく近くに、従来報告されていない微細な運動が現れることに気づいた。
この運動を系統的に特徴づけ、既知のどの遷移構造とも合致しないことを、統計解析・線形安定性解析・数値フィルタ試験によって示した。
捉えた構造
入口境界層厚さにもとづくレイノルズ数2,857、マッハ数0.1〜0.3の条件で計算した平板境界層において、遷移が始まる手前の壁のごく近くに、次の特徴をもつ構造が繰り返し現れた。
| 項目 | 本研究で得られた値 | 比較 |
|---|---|---|
| 流れ方向・スパン方向の波長 | 境界層厚さの0.35〜0.55倍 | TS波の中立波長は約55.7倍（100〜160倍の隔たり） |
| 壁からの広がり | 境界層厚さの14.3%以内で明瞭に確認できる（最も強いのは1.43%の高さ、弱いピークは境界層端まで残る） | 既知の遷移構造はいずれも境界層厚さ全体に広がる |
| 向き | スパン方向速度変動が約45度 | 乱流バンドは約25度 |
| 振幅 | スパン方向速度変動で主流速度の約0.05% | 平均速度に対して十分小さく、逆流は生じない |
| 壁単位換算の波長 | 約56〜89 | 乱流境界層のストリーク間隔は約100 |
表1. 本研究で捉えた近壁コヒーレント構造の特徴。長さはいずれも入口境界層厚さで規格化した値。
壁のごく近くでは、負の圧力変動と負の流れ方向速度変動が局所的に重なる。そして、流れ方向にこの低圧域を通り抜けるにつれて、壁垂直方向の速度が下向きから上向きへと反転する。渦による汲み上げに似た、圧力と速度の結合である。
この構造のスケールと向きは、格子解像度を上げても、マッハ数を変えても、上流の擾乱を壁面近傍に与えても自由流側に与えても変わらなかった。また、擾乱を与えない層流計算における数値的な雑音のレベルは、この構造の振幅から4桁（約5,000分の1）以上小さいものであった。これらは、この構造が数値計算に由来する見かけのものではなく、物理的な応答であることを支持する。
約45度に傾いた周期的な流れ構造が現れる。赤が正、青が負の変動を表す。
既知の構造との違い
壁の近くに現れる秩序だった構造はこれまでに多く報告されており、本研究の構造も、見た目だけでは既知のものと混同されやすいところがある。そこで論文では、代表的な既知構造との比較を明示した。
| 既知の構造 | 代表的な波長 （境界層厚さ比） | 本構造との違い |
|---|---|---|
| TS波 | 約55.7 | 二次元的でスパン方向に一様。本構造は三次元的で波長が100〜160分の1 |
| ストリーク | スパン方向が約8.2 流れ方向は10以上 | 流れ方向に細長い。本構造は縦横の波長がほぼ等しい |
| 斜め波 | 10〜100のオーダー | TS波どうしの非線形相互作用で生じる。スケールが桁違いに大きい |
| パフ・スラグ 乱流バンド | 100以上 | すでに乱流化した領域の大規模な間欠パターン |
| 乱流境界層の 45度X字パターン | 約1（境界層厚さと同程度） | 相関の形はよく似るが、十分発達した乱流中の現象で、波長が1〜2桁大きい |
| バックフロー （逆流事象） | 壁単位で20程度 | 粘性底層に限られた稀な偶発事象。本構造は空間的に周期的で、逆流は生じない |
表2. 既知の近壁構造との比較。波長はいずれも入口境界層厚さで規格化した値。
比較にあたっての留保
本構造の波長を壁単位に換算すると約56〜89となり、乱流境界層で古くから知られるストリーク間隔（約100）と同じオーダーになります。ただし本構造は層流に近い境界層に現れ、摩擦速度が乱流時よりはるかに小さいため、壁単位での単純な比較には注意が必要です。本研究の主張は、境界層厚さで規格化したスケールで見たときに、既知のいずれの遷移構造とも合致しない、という点にあります。
また、本構造は流れの側の現象であり、リブレットやDMRのような表面形状そのものではありません。これらの表面技術と直接対応づけられる段階には至っていません。
研究手法
高解像度の直接数値計算（DNS）
三次元圧縮性ナビエ・ストークス方程式を、6次精度コンパクト差分スキームで解いた。レイノルズ数2,857、マッハ数0.1／0.2／0.3。解像度の異なる2種類の計算格子で同じ条件の計算を行い、結果が変わらないことを確認している。風洞実験における壁面粗さや格子乱流に相当する擾乱として、上流域にランダムな体積力を与えた。
二点空間相関による統計的確認
ある瞬間の流れ場について、広い領域の格子点にわたって空間平均を取り、速度・圧力・密度の各成分の二点空間相関を求めた。45度に傾いたX字型のパターンは、スパン方向速度の相関に明瞭に現れる。
線形安定性解析（オア・ゾンマーフェルト方程式）
チェビシェフ選点法により中立安定曲線を求めた。本研究のレイノルズ数におけるTS波の中立波長は境界層厚さの約55.7倍であり、今回の構造は中立曲線の不安定領域の外側に大きく外れる。古典的なモーダル不安定性は起源として除外される。
壁面近傍の数値フィルタ試験
境界層厚さの11.4%以内の変動だけを選択的に抑える数値フィルタを設け、下流の渦構造への影響を調べた。広帯域の低域通過フィルタに加え、本構造の波数のみを狙ったスパン方向スペクトルフィルタも実装している。
非モーダル過渡成長の解析
目標時間を幅広く変えて最適擾乱を計算した。短い時間スケールで最も成長しやすいモードの波長（流れ方向0.571、スパン方向0.457）は、DNSの観測値とおおむね一致する。ただし成長率は1.6倍程度と控えめであり、線形理論だけでは十分に説明できない。
本研究で示していないこと
未解明の点
生成機構。非モーダル過渡成長との波長の一致は示されたが、成長率が控えめであり、線形理論だけでこの構造の形成を説明することはできない。因果関係の特定は今後の課題である。
自律性。上流からの擾乱を与え続けないと構造は減衰する。現時点では自律的に成長する不安定性ではなく、外部擾乱に対する境界層の応答（受容性）として理解される。なお、実際の飛行や走行では大気乱れや表面粗さによる擾乱が常に存在するため、この条件はむしろ現実の流れに近い。
独立検証。本研究は単一のDNSコードによる結果であり、独立な計算コードによる再現や実験による確認は行われていない。ただし、境界層厚さの数パーセントという壁面近傍で主流速度の0.05%程度の変動を捉えることは、現在の計測技術では容易ではない。
DMRとの関係について
本成果は、DMRによる抵抗低減の機構を説明するものではありません。
DNSで最も遷移遅延効果が高かった突起の設置間隔は境界層厚さと同程度であり、本構造の波長の2〜3倍にあたります。実験で効果を示したDMRコーティングの間隔は、さらに小さいものでした。したがって、突起間隔と構造の波長が単純に共鳴しているという説明は成り立ちません。
両者に関係があるのかどうかは未解決の問いであり、DMRの形状と分布を系統的に変えた直接数値計算によって調べていく予定です。
論文情報
2026年9月27日採択。2026年10月16日公開。
研究助成
JSPS科研費（19K14880、23K26030）／ JST CREST（JPMJCR24Q6）／ JST創発的研究支援事業（JPMJFR222R）
数値シミュレーションは、東北大学流体科学研究所 スーパーコンピュータ AFFINITY II を用いて行いました。