2.23次元非圧縮性熱流体計算プログラム CSMOS の定評価.pdfVIP

2.2 3 次元非圧縮性熱流体計算プログラム COSMOS の測定評価 （株）豊田中央研究所　　堀之内 成明 2.2.1 はじめに 自動車の設計?開発における熱流体 CAE (Computer-Aided Engineering) は，日常的なものとして 様々な製品?部品を対象に用いられてきているが，実際には，非定常な乱流に加えて，熱，騒音，運動 （移動境界），反応などが連成するような現象も多く，これらに対して信頼性の高い結果を得るために は，計算格子幅や，時間ステップ間隔を十分に小さくする必要が出てくるため，その結果として必然 的に計算規模が大きくなる． このような大規模熱流体シミュレーションは，従来ではベクトル演算方式の計算機で行われること が多く，プログラムもこの上で高速計算が可能となるようにチューニング，ならびにアルゴリズム開 発がなされてきた．本稿では，このようなプログラムを，ベクトル計算機での性能も維持しつつ，「京」 をはじめとするマルチコアのスカラー型計算機においても高い実効性能を満たすことを目的として， 高速化に向けた改良を試みる． 2.2.2 プログラム概要 ここでは，（株）豊田中央研究所で開発した非定常非圧縮性乱流計算プログラム COSMOS を対象に チューニングを行う．このプログラムは，空力や空力騒音といった問題に適用するため，自動車まわ りの高 Reynolds 数非定常流れを精度よく計算することを目的に開発された [1]． 2.2.2.1 プログラムの特徴 本プログラムは，以下の特徴を有する． ? Large Eddy Simulation (LES) による非定常乱流の計算に関して，以下の独自技術を開発し，実 問題に適用し得る計算精度と実用性を両立させている． – Grid-scale (GS) 成分計算の離散化スキームとして，質量，運動量のみならず，運動エネル ギーまでをも高い精度で保存させ，数値粘性を用いなくても安定な計算を可能とする「改 良型コロケーション格子法 [2]」 – Subgrid-scale (SGS) 成分の乱流モデルとして，様々な問題に対してモデル係数のチューニ ングが不要で，かつダンピング関数を必要としない「混合時間スケール SGS モデル [3]」 ? 計算格子には，計算精度に優れているというメリットから「境界適合型構造格子」を用いる．こ れは，碁盤目状の格子を，対象となる物体の表面に沿わせ，且つ表面近くに格子点を集中させた ものである．また，三角形（2 次元の場合）やテトラ（3 次元の場合）を用いる非構造格子に比 べて，プログラミングに際して間接参照（リストアクセス）の必要性が軽減でき，計算効率の 観点からも有利である．但し，車両などの複雑形状に対しては，格子生成に多くの工数やノウ ハウを必要とするというデメリットもあり，そのような計算領域全体を一つの格子で覆うこと 自体が現実的には困難であったり，計算格子の歪みが大きくなり質が低下するといった課題も ある．これを解決するために，局所的な形状に適合させた高品質な部分的構造格子（以下，部 分格子）をパッチワーク的に重ね合わせる「重合格子法」を採用している．車両形状に対する 重合格子の例を図 1 に示す． ? 非圧縮性流れの基礎方程式は，連続の式と Navier-Stokes 方程式（運動方程式）であるが，SMAC 法に準じて，これらから圧力に関する Poisson 方程式を導いて解く．時間方向の離散化には， Crank-Nicolson 法に基づく陰解法を採用しており，これから速度場に関する連立一次方程式が 導かれる．一方，圧力の Pois