1次元移流方程式でのCIP解法とは? わかりやすく解説

Weblio 辞書 > 辞書・百科事典 > ウィキペディア小見出し辞書 > 1次元移流方程式でのCIP解法の意味・解説 

1次元移流方程式でのCIP解法

出典: フリー百科事典『ウィキペディア(Wikipedia)』 (2021/04/02 08:51 UTC 版)

CIP法」の記事における「1次元移流方程式でのCIP解法」の解説

CIP法1次元移流方程式高精度に解く解法である。1次元移流方程式は次式で与えられる。 ∂ f ∂ t + c ∂ f ∂ x = 0 {\displaystyle {\frac {\partial f}{\partial t}}+c{\frac {\partial f}{\partial x}}=0} ここで、cは移流速度である。 CIP法では、格子点の値 g ( = ∂ f ∂ x ) {\displaystyle g(={\frac {\partial f}{\partial x}})} についても同時に移流計算を行うことが特徴である。空間微分値 g {\displaystyle g} に対す移流方程式は、上の移流方程式空間に関して微分することで得られ、以下のようになる。 ∂ g ∂ t + c ∂ g ∂ x = 0 {\displaystyle {\frac {\partial g}{\partial t}}+c{\frac {\partial g}{\partial x}}=0} 時刻 n {\displaystyle n} における値 f {\displaystyle f} とその微分値 g {\displaystyle g} が格子点上の点 i {\displaystyle i} 、 i u p {\displaystyle iup} (点 i u p {\displaystyle iup} は点 i {\displaystyle i} の上流点、つまり移流速度 c > 0 {\displaystyle c>0} なら i u p = i − 1 {\displaystyle iup=i-1} である)において既知とすると、この2点間のプロファイル(つまり形状)は以下のように3次多項式表される。ここで、上付き添字時刻を、下付き添字格子点番号をあらわす。 F i n ( x ) = a i ( x − x i ) 3 + b i ( x − x i ) 2 + g i n ( x − x i ) + f i n {\displaystyle F_{i}^{n}(x)=a_{i}(x-x_{i})^{3}+b_{i}(x-x_{i})^{2}+g_{i}^{n}(x-x_{i})+f_{i}^{n}} このようにプロファイル補間関数3次関数表現することがCubicInterpolated Pseudo-Particle Schemeたる所以である。ここで、係数 a i {\displaystyle a_{i}} 、 b i {\displaystyle b_{i}} は、 a i = g i n + g i u p n D 2 + 2 ( f i nf i u p n ) D 3 {\displaystyle a_{i}={\frac {g_{i}^{n}+g_{iup}^{n}}{D^{2}}}+{\frac {2(f_{i}^{n}-f_{iup}^{n})}{D^{3}}}} b i = 3 ( f i u p nf i n ) D 22 g i n + g i u p n D {\displaystyle b_{i}={\frac {3(f_{iup}^{n}-f_{i}^{n})}{D^{2}}}-{\frac {2g_{i}^{n}+g_{iup}^{n}}{D}}} のようになる。ただし、移流速度 c > 0 {\displaystyle c>0} のとき D = − Δ x {\displaystyle D=-\Delta x} 、 i u p = i − 1 {\displaystyle iup=i-1} であり、移流速度 c < 0 {\displaystyle c<0} のとき D = Δ x {\displaystyle D=\Delta x} 、 i u p = i + 1 {\displaystyle iup=i+1} である。 適合条件式により、 F i n ( 0 ) = f i n {\displaystyle F_{i}^{n}(0)=f_{i}^{n}} 、 ∂ F i n ( 0 ) ∂ x = g i n {\displaystyle {\frac {\partial F_{i}^{n}(0)}{\partial x}}=g_{i}^{n}} 、 F i n ( D ) = f i u p n {\displaystyle F_{i}^{n}(D)=f_{iup}^{n}} 、 ∂ F i n ( D ) ∂ x = g i u p n {\displaystyle {\frac {\partial F_{i}^{n}(D)}{\partial x}}=g_{iup}^{n}} が成り立つので、上式において係数 a i {\displaystyle a_{i}} 、 b i {\displaystyle b_{i}} が求まるこのように格子点上の点において微分値 g {\displaystyle g} も与えられるので、格子間のプロファイル3次多項式補間することができ、精度の高い計算が可能となる。 対象とする問題移流方程式であるので、次の時刻 n + 1 {\displaystyle n+1} での値 f i n + 1 {\displaystyle f_{i}^{n+1}} と微分g i n + 1 {\displaystyle g_{i}^{n+1}} は、この2点間のプロファイルを c Δ t {\displaystyle c\Delta t} だけ移動することで得られる。つまり、 ξ = − c Δ t {\displaystyle \xi =-c\Delta t} として次式のようになるf i n + 1 = F i n ( ξ ) = a i ξ 3 + b i ξ 2 + g i n ξ + f i n {\displaystyle f_{i}^{n+1}=F_{i}^{n}(\xi )=a_{i}\xi ^{3}+b_{i}\xi ^{2}+g_{i}^{n}\xi +f_{i}^{n}} g i n + 1 = ∂ F i n ( ξ ) ∂ x = 3 a i ξ 2 + 2 b i ξ + g i n {\displaystyle g_{i}^{n+1}={\frac {\partial F_{i}^{n}(\xi )}{\partial x}}=3a_{i}\xi ^{2}+2b_{i}\xi +g_{i}^{n}}

※この「1次元移流方程式でのCIP解法」の解説は、「CIP法」の解説の一部です。
「1次元移流方程式でのCIP解法」を含む「CIP法」の記事については、「CIP法」の概要を参照ください。

ウィキペディア小見出し辞書の「1次元移流方程式でのCIP解法」の項目はプログラムで機械的に意味や本文を生成しているため、不適切な項目が含まれていることもあります。ご了承くださいませ。 お問い合わせ



英和和英テキスト翻訳>> Weblio翻訳
英語⇒日本語日本語⇒英語
  

辞書ショートカット

すべての辞書の索引

「1次元移流方程式でのCIP解法」の関連用語

1
16% |||||

1次元移流方程式でのCIP解法のお隣キーワード
検索ランキング

   

英語⇒日本語
日本語⇒英語
   



1次元移流方程式でのCIP解法のページの著作権
Weblio 辞書 情報提供元は 参加元一覧 にて確認できます。

   
ウィキペディアウィキペディア
Text is available under GNU Free Documentation License (GFDL).
Weblio辞書に掲載されている「ウィキペディア小見出し辞書」の記事は、WikipediaのCIP法 (改訂履歴)の記事を複製、再配布したものにあたり、GNU Free Documentation Licenseというライセンスの下で提供されています。

©2025 GRAS Group, Inc.RSS