CIP法の二つの顔：Constrained Interpolation Profile と Continuous Interior Penalty

Introduction

CIP という略語には，少なくとも二つの重要な意味があります．一つは，移流問題において値と導関数を用いてプロファイルを復元する Constrained Interpolation Profile です．もう一つは，有限要素法において要素境界での導関数ジャンプを抑える Continuous Interior Penalty です．本連載では，この二つを「プロファイルをどのように作り，運び，正則化するか」という観点から整理します．論文投稿前の研究構想メモとして書いたものです．今後，定式化・安定性評価・数値実験の結果に応じて内容を更新します．

1．なぜ今 CIP を見直すのか

移流方程式の数値計算では，解の形をなるべく壊さずに運ぶことが重要です．鋭い層，段差，波の山，あるいは局所的なピークを含む問題では，単純な低次補間を使うと数値拡散により解が丸くなります．一方で，高次補間を安易に使うと，非物理的な振動が生じることがあります．

この問題に対して，古典的な Constrained Interpolation Profile 型の CIP 法は，値だけでなく導関数も用いて，セル内のプロファイルを高精度に復元するという考え方を採ります．一方，有限要素法の世界で CIP と呼ばれる Continuous Interior Penalty 法は，連続有限要素空間を用いながら，要素間の導関数ジャンプにペナルティを課して安定化する方法です．

名前は同じ CIP ですが，出発点はかなり違います．前者は「プロファイルを復元して運ぶ」方法，後者は「有限要素解の不安定な成分を抑える」方法です．しかし，どちらも導関数，プロファイル，正則化という共通語を持っています．ここに，短いが面白い研究テーマがあります．

2．二つの CIP

次の二つは，通常は別々の文脈で語られます．しかし，1次元移流方程式を共通の舞台にすると，両者の距離は意外に近くなります．Constrained Interpolation Profile は，導関数情報を使ってプロファイルを作ります．Continuous Interior Penalty は，導関数ジャンプを抑えることで，プロファイルとして使える情報を整えます．この見方を採ると，後者は単なる安定化項ではなく，プロファイル生成・正則化のための機構として見えてきます．

Constrained Interpolation Profile

値と導関数を同時に使い，セル内の補間プロファイルを作る方法です．典型的には，節点値と節点導関数から三次 Hermite 型の多項式を構成し，特性の足で関数値を評価します．移流問題に対して，低拡散で高解像度な計算を狙う方法です．

Continuous Interior Penalty

連続有限要素空間を用いながら，要素境界での導関数ジャンプを安定化項として加える方法です．移流優勢問題や移流拡散問題において，Galerkin 法だけでは不足する安定性を補う役割を持ちます．

3．共通モデル：1次元移流方程式

まずは，空間1次元の線形移流方程式から始めます． \[ \partial_t u + b\,\partial_x u = 0, \quad x\in (0,1),\quad t>0, \]

ここで，速度 \(b\) は簡単のため正の定数とします．この方程式の解は，特性曲線に沿って運ばれます．したがって，時間刻み \(\Delta t\) に対して，形式的には

\[ u(x,t^{n+1}) = u(x-b\Delta t,t^n) \]

と書けます．問題は，特性の足 \(x-b\Delta t\) が一般にメッシュの節点とは一致しないことです．したがって，数値解を更新するためには，どこかで補間または再構成が必要になります．

この「特性の足でどう評価するか」という一点に，CIP 型の考え方が現れます．値だけで評価するのか．導関数も使うのか．導関数を点値として持つのか．弱形式で持つのか．あるいは，有限要素の安定化によって導関数プロファイルを整えてから使うのか．本連載の中心は，まさにこの分岐です．

4．古典的 CIP の見方：値と勾配を運ぶ

古典的な Constrained Interpolation Profile 型の方法では，各節点で

\[ u_i^n \simeq u(x_i,t^n), \quad g_i^n \simeq \partial_x u(x_i,t^n) \]

のように，値 \(u_i^n\) と導関数 \(g_i^n\) を保持します．そして，区間 \([x_i,x_{i+1}]\) 上で，これらの値を満たす三次多項式を作ります．これは Hermite 型の再構成と見ることができます．

この方法の魅力は，セル内の形状をある程度保ったまま，特性の足で値を評価できる点にあります．値だけを線形補間するよりも多くの情報を使うため，低拡散な計算が期待できます．

注意． ただし，点値としての導関数自由度を使うことは，有限要素解析の立場からは少し扱いにくい面があります．弱解を対象にしたい場合や，低正則な解を扱いたい場合には，点値導関数を前提にすること自体が強すぎる可能性があります．

5．有限要素法から見た問い

有限要素法の立場から見ると，古典的 CIP には次の問いが自然に出てきます．

  1. 導関数プロファイルを，点値ではなく弱形式で保持できないか．
  2. semi-Lagrangian 型の更新を，射影や質量行列を含む有限要素の言葉で書けないか．
  3. 保存形の移流方程式に対して，質量保存と高解像度性を両立できないか．
  4. Continuous Interior Penalty を，導関数プロファイルの正則化として解釈できないか．

この問いのうち，前半は weak-moment CIP finite elements につながります．後半は CIP-regularised Hermite reconstruction in conservative Lagrange–Galerkin FEM につながります．つまり，この連載は，二つの短報候補のための前段階でもあります．

6．弱モーメントという方向

第2回以降で詳しく扱いますが，重要な方針は，導関数情報を点値として持つのではなく，弱い意味で持つことです．たとえば，セル \(I_j\) 上で補助変数 \(q_h\) を導入し，

\[ \int_{I_j} q_h \varphi_h\,dx \simeq \int_{I_j} \partial_x u_h\,\varphi_h\,dx \]

のような関係で，導関数プロファイルを弱形式として管理する考え方です．これはまだ最終的な定式化ではありませんが，点値導関数に依存しない CIP 的再構成へ進むための自然な出発点です．

この方向では，CIP を単なる補間公式ではなく，有限要素空間，射影，弱モーメント，再構成からなる一つの構造として捉え直します．ここが，通常の「CIP 法の紹介」とは異なる点です．

7．Continuous Interior Penalty をどう見るか

Continuous Interior Penalty 型の安定化では，要素境界における導関数ジャンプを罰します．1次元のイメージでは，内部節点 \(x_i\) に対して，

\[ \sum_i \gamma h_i^2 \,[\partial_x u_h]_{x_i}\,[\partial_x v_h]_{x_i} \]

のような項が現れます．通常，これは移流優勢問題に対する安定化項として説明されます．しかし，本連載では，もう一つの見方を採ります．

見方の転換． Continuous Interior Penalty は，要素間の高周波な導関数ジャンプを抑えることで，有限要素解からより滑らかな導関数プロファイルを取り出すための正則化と見ることができます．

もしこの見方が有効なら，Continuous Interior Penalty によって整えられた導関数情報を，Hermite 型再構成に使うことができます．すると，Constrained Interpolation Profile と Continuous Interior Penalty の間に，「プロファイル」という共通の橋が架かります．少し大げさに言えば，CIP と CIP がここで出会います．略語が同じなのは偶然かもしれませんが，研究テーマとしては都合がよい偶然です．

8．本連載の予定

本連載では，次の順番で話を進めます．

  1. 第1回： CIP 法の二つの顔を整理する．
  2. 第2回： 1次元移流方程式と semi-Lagrangian 更新を導入する．
  3. 第3回： 古典的 CIP の Hermite 再構成を確認する．
  4. 第4回： 点値導関数ではなく，弱モーメントでプロファイルを持つ方法を考える．
  5. 第5回： Continuous Interior Penalty を導関数プロファイルの正則化として見る．
  6. 第6回： 保存形 Lagrange–Galerkin 有限要素法との接続を整理する．
  7. 第7回： 数値実験の最小コードを作り，古典的 CIP，弱モーメント型，CIP 正則化型を比較する．

9．この段階での研究メモ

現時点では，この構想を大きな主力論文としてではなく，1次元から始める短報型の研究として育てるのが現実的です．まずは，導関数プロファイルをどのように定義し，どのように特性写像上で使うかを明確にする必要があります．その後で，安定性，保存性，誤差分解を順に確認します．

特に重要なのは，古典的 CIP の直感的な強さを尊重しつつ，有限要素法の弱形式に載せ直すことです．この変換がうまくいけば，CIP 法を「工学的に有効な補間テクニック」から，「解析可能な有限要素再構成手法」へ近づけられます．

次回は，1次元移流方程式を題材として，特性曲線，semi-Lagrangian 更新，および Lagrange–Galerkin 型の弱形式を整理します．ここからようやく，数式が少しだけ本気を出します．ただし，まだ暴れません．

参考文献・関連資料

  1. H. Takewaki, A. Nishiguchi and T. Yabe, “Cubic interpolated pseudo-particle method (CIP) for solving hyperbolic-type equations,” Journal of Computational Physics, 61, 261--268, 1985. ScienceDirect
  2. T. Yabe, F. Xiao and T. Utsumi, “The constrained interpolation profile method for multiphase analysis,” Journal of Computational Physics, 169, 556--593, 2001. ScienceDirect
  3. E. Burman and A. Ern, “Continuous interior penalty hp-finite element methods for advection and advection-diffusion equations,” Mathematics of Computation, 76, 1119--1140, 2007. UCL Discovery
  4. H. Ishizaka, “Ongoing Papers, Next Topics,” Publications/Talks page. ISHIZAKA Hiroki Website

ページの先頭へ