Mac!!!ガジェットとはちょっとした小物・気のきいた小道具のこと。Macのコネタを紹介。スマートフォンなどなど。ishogaki0はじめに縁あって流体解析なんぞやっていて,オープンソースのOpenFOAMをよく使うので,それについてメモ書き。ややこしいコードなので間違ってるかもしれませんが,そのときはご指摘ください。ブログの本筋からは大分外れててすいません。}pEqn.flux()関数の呼び出しがあります。Globalによるソースコードの読解で調べてみました。GDBとGlobalGDBGDBはGCCのデバッガです。OpenFOAMでGDBを使えるようにするには,コンパイルオプションWM_COMPILE_OPTIONをDebugにしてAllwmakeします。詳細は下記なんかがわかりやすいです。Globalはソースコードにタグをつけてくれるソフトです。関数のソースがどこにあるかなどの情報が見やすくなります。EclipseのF3を押しても同様のことができますが,候補が一覧で表示されるので,個人的にはEclipseよりも使いやすいです。除外の仕方は以下を参考にしました。http://d.hatena.ne.jp/ohtorii/20110219/1298092604pisoFOAMのソースコードを見るはじめにGlobalでflux()関数を見るGlobalでpisoFoam.Cを見ると,flux()へのリンクができています。リンクをたどると以下のような候補が表示されます。limitedSurfaceInterpolationScheme<Type>::fluxごちゃごちゃしてますが,pEqnがfvMatrixのオブジェクトなので,fvMatrix.Cのflux()を見てみます。}途中boundaryFieldの量を各パッチで計算したりしています。Nozakiの資料にあるように,flux()は圧力勾配から来る流束の修正を行っているように見えます。}faceFluxCorrectionPtr_の示す値が追加されています。これが何をやっているかを調べてみます。faceFluxCorrectionPtr_についてfaceFluxCorrecitonPtr_のゲッターとしてfaceFluxCorrectionPtr()関数があります。}pisoFoamのfaceFluxCorrectionPtr()関数の詳細について次で調べてみます。GDBでデバッグしてみるここからGDBによるデバッグでfaceFluxCorrectionPtr()関数について調べてみます。デバッグのためのセットアップGDBの結果を見るのにはEmacsが便利だと思います。以下の設定を.emacs.d/init.elや.emacsにします。Emacs+GDBは以下参照。gdbでGDBを起動します。起動時にはpisoFOAMが立ち上がるようにします。bufferを閉じておきます。しないとデバッグが途中で停止します。flowのチュートリアルをもとにしてデバッグします。fvSchemesはそのままです。非直交補正のについて調べるので,laplacianスキーム,snGradスキームがcorrectedになっているか確認して下さい。//fvSolutionは以下のようにしました。0;}メッシュ数が多いとデバッグで見ていくときに何かと大変なので,セル数5個のメッシュを作成しました。以下がそのメッシュを生成するためのblockMeshdictです。blockMeshでメッシュを生成します。初期条件,境界条件は適当に設定します。//GDBでブレークポイントを設定する実際にブレイクポイントを設定して,pisoFoamのfaceFluxCorrectionPtr()関数の挙動を見てみます。fvMatrix.Hはsrc/finiteVolume/fvMatrices/fvMatrix/fvMatrix.Hにありますが,計算時に参照されるのはsrc/finiteVolume/lnInclude/fvMatrix.Hです。ここにブレイクポイントを置きます。runで計算を走らせます。コールスタックを抜き出したのが以下。pisoFoam.C:102#1をクリックするとその呼出が表示されます。gaussLapalacianSchemes.Cの該当箇所を見ると,gaussLapacianSchemes.CC++declareFvmLaplacianScalarGamma(scalar);12declareFvmLaplacianScalarGamma(scalar);となっています。これはマクロでその定義はすぐ上にあります。}64行目以降が非直交補正に関係有るところです。fvSchemesでfluxRequiredを指定してあるので,68行目でfaceFluxCorrectionPtr()関数に対して,snGradの非直交補正による補正量が代入されています。74行目でsource項に非直交補正による補正を代入しています。このようにGDBを使うことで,OpenFOAMの深い関数呼び出しの関係を調べられます。途中マクロ展開で関数呼び出しを実装することがOpenFOAMの場合には多くあります。このため,関数にあたりをつけて,適切にブレイクポイントを設定する必要があります。次に71行目のcorrection()関数についてGlobalで見てみます。gammaMagSf*this->tsnGradScheme_().correction(vf)GlobalでtsnGradScheme_().correction()を調べるcorrection()関数の候補はたくさん表示されます。ここではcorrectedSnGradを使っているので,correctedSnGrad.Cのcorrection()を見てみます。実装が以下で,correction()関数の計算には102行目にあるfullGradCorrection()関数が関わっています。}fullGradCorrection()はその上に定義されています。}mesh.nonOrthCorrectionVectors()とlinearで補間されたgradの値の内積を計算しています。)にある,非直交補正の式と対応します。以上の内容はNozakiさんの資料にも説明がありますので,そちらを参考にするとわかりやすいです。まとめflux()関数を使った流束の補正とfaceFluxCorrectionPtr()関数から始まる非直交補正の詳細について調べてみました。GDBとGlobalでいろいろたどれるので,使ってみてください。WordPress.com.。