Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
渡辺 勢也*; 河原 淳*; 青木 尊之*; 杉原 健太; 高瀬 慎介*; 森口 周二*; 橋本 博公*
Engineering Applications of Computational Fluid Mechanics, 17(1), p.2211143_1 - 2211143_23, 2023/00
被引用回数:1 パーセンタイル:52.66(Engineering, Multidisciplinary)津波氾濫や大雨の斜面災害では、多くの浮遊物や流木、丸太などが流れに含まれる。剛体衝突による構造物の被害は、水圧による被害よりも深刻である。浮遊物体を含む自由表面流れを研究するためには、大規模計算が可能で高性能な自由表面流のシミュレーションコードの開発が必要となる。本研究ではキュムラント格子ボルツマン法と粒子ベースの剛体シミュレーションを組み合わせた単相自由表面モデルを提案する。剛体間の接触相互作用は個別要素法で計算される。解析精度の向上と計算の高速化のために8分木ベースの局所細分化格子法を導入し、自由表面や固体表面近傍には高解像度の格子を割り当てた。提案モデルの精度検証のために、八戸工業大学と神戸大学の15m水槽と70m水槽で2種類の津波実験を行った。シミュレーションの結果、漂流速度、捕捉木片数、積層角について実験と良い一致を示した。
小野寺 直幸; 井戸村 泰宏; 長谷川 雄太; 下川辺 隆史*; 青木 尊之*
計算工学講演会論文集(CD-ROM), 27, 4 Pages, 2022/06
我々の研究グループでは、風況デジタルツインの実現に向けて、風況解析コードCityLBMを開発している。CityLBMは、計算領域周辺の境界条件をメソスケール気象データに同化させるナッジング法を導入することで、現実の風況を反映した解析が可能である。しかしながら、従来のナッジング法では、ナッジング係数が一定のため、大気状態が変化するような長時間解析の乱流強度を再現できない問題点が挙げられる。そこで、本研究では、パーティクルフィルタを用いた動的なナッジング・パラメータの最適化手法を提案する。CityLBMの検証として、米国オクラホマシティの風況実験に対する解析を実施した。シミュレーションと観測のそれぞれで得られる乱流強度の誤差を低減するようにナッジング係数を更新した結果、シミュレーションにおいて終日の大気境界層を再現できることを確認した。
長谷川 雄太; 青木 尊之*; 小林 宏充*; 井戸村 泰宏; 小野寺 直幸
Parallel Computing, 108, p.102851_1 - 102851_12, 2021/12
被引用回数:3 パーセンタイル:41.82(Computer Science, Theory & Methods)GPUスーパコンピュータに対して格子ボルツマン法(LBM: lattice Botltzmann method)およびforest-of-octreesに基づくブロック構造型の局所細分化格子(LMR: local mesh refinement)を用いた空力解析コードを実装し、その性能を評価した。性能評価の結果、従来の空間充填曲線(SFC; space-filling curve)に基づく領域分割アルゴリズムでは、本空力解析において袖領域通信のコストが過大となることがわかった。領域分割の改善手法として本稿では挿し木法を提案し、領域分割の局所性とトポロジーを改善し、従来のSFCに基づく手法に比べて通信コストを1/31/4に削減した。強スケーリング測定では、最大で1.82倍の高速化を示し、128GPUで2207MLUPS(mega-lattice update per second)の性能を達成した。弱スケーリング測定では、8
128GPUで93.4%の並列化効率を示し、最大規模の128GPU計算では44.73億格子点を用いて9620MLUPSの性能を達成した。
小野寺 直幸; 井戸村 泰宏; 長谷川 雄太; 中山 浩成; 下川辺 隆史*; 青木 尊之*
Boundary-Layer Meteorology, 179(2), p.187 - 208, 2021/05
被引用回数:13 パーセンタイル:75.07(Meteorology & Atmospheric Sciences)汚染物質の拡散解析手法CityLBMは、GPUスーパーコンピュータ上において、適合細分化格子(AMR)法を適用する事で、数kmの解析領域の実時間解析が可能である。本論文では、CityLBMの検証としてオクラホマ市で実施された野外拡散実験(JU2003)に対する解析を実施した。計算条件として、Weather Research and Forecasting(WRF)モデルを用いた風況条件および、建物と植生を考慮した地表面データをCityLBMに与えることで、JU2003の実験条件を再現した。さらにアンサンブル計算の実施により、乱流の不確実性を軽減した。汚染物質の時間平均濃度および最大値を実験測定値と比較した結果、アンサンブル計算により解析精度を向上すると共に、2m解像度・4km四方の解析では、24個の計測値に対して70%の高い割合でFactor2を満たす事を確認した。
長谷川 雄太; 青木 尊之*; 小林 宏充*; 井戸村 泰宏; 小野寺 直幸
計算工学講演会論文集(CD-ROM), 26, 6 Pages, 2021/05
Forest-of-octreesに基づく局所格子細分化法(LMR)を導入した格子ボルツマン法(LBM)に基づく空力解析コードに対し、挿し木法による領域分割の改善手法を提案した。従来の空間充填曲線に基づく領域分割法は、適合格子細分化法(AMR)やLMRで広く用いられているものの、GPUスパコン向けに実装された本空力解析コードにおいては袖領域通信が増大し計算のボトルネックとなるうることが確認された。本研究で提案する挿し木法は、粗い等間隔格子状の領域分割と細かい空間充填曲線に基づく分割のハイブリッドによる手法である。挿し木法により、領域分割の局所性と幾何形状が改善しており、通信量が従来の空間充填曲線に基づく手法に比べて3分の1に削減された。8GPU並列による性能検証では、コード全体で1.23倍の高速化が確認された。また、強スケーリングにおいてさらに性能の改善が見られ、128GPUの強スケーリングにおいては、従来手法に比べて1.82倍の高速化を示し、2207MLUPS (mega-lattice update per second)の計算性能を達成した。
小野寺 直幸; 井戸村 泰宏; 長谷川 雄太; 下川辺 隆史*; 青木 尊之*
計算工学講演会論文集(CD-ROM), 26, 3 Pages, 2021/05
本研究では、二相流体解析コードJUPITER-AMRに対して、圧力ポアソン方程式に対する混合精度前処理手法を開発した。マルチグリッド前処理手法として、3段のVサイクルの幾何学的MG法およびキャッシュを再利用したSOR(CR-SOR)法を適用した。原子力工学問題での性能測定として、バンドル体系に対する多相流体解析を実施した。計算速度として、単精度演算を適用する事で、倍精度演算の75%へと削減すると共に、強スケーリング性能においては、32台から96台のGPUを利用する事で1.88倍を実現した。
小野寺 直幸; 井戸村 泰宏; 長谷川 雄太; 山下 晋; 下川辺 隆史*; 青木 尊之*
Proceedings of International Conference on High Performance Computing in Asia-Pacific Region (HPC Asia 2021) (Internet), p.120 - 128, 2021/01
被引用回数:0 パーセンタイル:0.01(Computer Science, Hardware & Architecture)本研究では、二相流体解析コードJUPITERに対して、マルチグリッド前処理付き共役勾配(MG-CG)法を開発した。MG法は、3段のVサイクルMG法に基づいて構築し、各段に対して、RB-SOR法およびGPUのキャッシュを再利用したCR-SORを開発・適用した。性能測定として、バンドル体系に対する気液二相流体解析を行った。RB-SOR法およびCR-SOR法を適用したMG-CG法では、MG法を適用しないPCG法と比較して、収束までの反復回数を15%と9%以下に削減するとともに、3.1倍, 5.9倍の計算速度が達成された。以上の結果から、本研究で開発したMG-CG法は、GPUを用いたスーパーコンピュータ上にて、効率的に大規模な二相流体解析が可能であることが示された。
小野寺 直幸; 井戸村 泰宏; 朝比 祐一; 長谷川 雄太; 下川辺 隆史*; 青木 尊之*
第34回数値流体力学シンポジウム講演論文集(インターネット), 2 Pages, 2020/12
本研究では、二相流体解析コードJUPITERにおいて、圧力ポアソン方程式に対するマルチグリッド前提共役勾配(MG-CG)ソルバーを開発した。プログラムの開発言語として、C++およびCUDAを用いることで、様々なコンピュータプラットフォームに対応した。CG解法の主な計算カーネルは、GPU, CPU、およびARM上において、ルーフライン性能の0.40.75と妥当な性能を達成した。一方で、SpMVカーネルでは、ARM上において、大幅な性能劣化が確認された。その原因を調査したところ、SpMVカーネル内にて関数呼び出しを行うことで、コンパイラの最適化が働かないことが確認された。
小野寺 直幸; 井戸村 泰宏; Ali, Y.*; 山下 晋; 下川辺 隆史*; 青木 尊之*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.210 - 215, 2020/10
本研究では、ブロック型局所細分化(AMR)法に基づくPoisson解法のGPU高速化を実施した。ブロック型AMR法はGPUに適したデータ構造であり、複雑な構造物で構成された原子炉等の解析に必須な解析手法である。これに、最新の前処理手法であるマルチグリッド(MG)法を共役勾配(CG)法へと組み合わせることで、計算の高速化を実現した。MG-CG法を構成する計算カーネルをGPUスーパーコンピュータであるTSUBAME3.0上にて測定した結果、ベクトル-ベクトル和、行列-ベクトル積、およびドット積の帯域幅は、ピークパフォーマンスの約60%となり、良好なパフォーマンスを実現した。更に、MG法の前処理手法として、3段のVサイクル法および各段に対してRed-Black SOR法を適用した手法を用いて、格子点の大規模問題の解析を実施した結果、元の前処理付きCG法と比較して、反復回数を30%未満に削減すると共に、2.5倍の計算の高速化を達成した。
小野寺 直幸; 井戸村 泰宏; Ali, Y.*; 下川辺 隆史*; 青木 尊之*
計算工学講演会論文集(CD-ROM), 25, 4 Pages, 2020/06
原子力機構では3次元多相流体解析手法としてJUPITERを開発している。本研究では、JUPITERの圧力Poisson方程式解法として、適合細分化格子(AMR)を用いたマルチグリッド前提条件付き共役勾配法(P-CG)を開発した。計算の高速化として、全ての計算カーネルはCUDAを用いて実装すると共に、GPUスーパーコンピュータ上にて高い性能を発揮する様に最適化した。開発したマルチグリッド圧力Poisson解法は、オリジナルのP-CG法と比較して約1/7の反復回数で収束することが確認された。また、TSUBAME3.0上で8から216GPUまでの強スケーリング性能測定により、更なる3倍の高速化が達成された。
青木 尊之*; 長谷川 雄太
自動車技術, 74(4), p.18 - 23, 2020/04
LESに基づくCFD計算を用いて自転車競技の空力解析を行った。単独での走行および24人の集団走行では、算出された抗力は風洞実験と良く一致した。競合する2集団の走行について、集団内の選手の配置を複数検討した。72人の選手の集団走行として、GPUスーパコンピュータで22.3億格子を用いた大規模空力解析を実施した。
稲垣 厚至*; Wangsaputra, Y.*; 神田 学*; Ycel, M.*; 小野寺 直幸; 青木 尊之*
SOLA (Scientific Online Letters on the Atmosphere) (Internet), 16, p.120 - 124, 2020/00
被引用回数:1 パーセンタイル:4.12(Meteorology & Atmospheric Sciences)都市境界層を対象とした風況解析により、内層および外層のスケーリングと乱流強度分布の類似性を検討した。計算条件として、現実的な建物形状の上に発達する中立条件を仮定すると共に、計算領域19.2km4.8km
高さ1kmに対して2m格子を設定した。乱流強度分布は計算領域内で局所的に定義できる。内層と外層に対してスケーリング則を当てはめることで、表面形状に関係なく、内層と外層内の乱流強度予測のばらつきを減らすことが可能となる。スケーリングされたプロファイル間のばらつきは、各層のスケーリングパラメーターの不一致に起因するが、長さまたは速度の比率からなる無次元パラメーターを導入することで、それらの類似性を示した。
下川辺 隆史*; 遠藤 敏夫*; 小野寺 直幸; 青木 尊之*
Proceedings of 2017 IEEE International Conference on Cluster Computing (IEEE Cluster 2017) (Internet), p.525 - 529, 2017/09
ステンシルに基づくCFDコードは、規則的なメモリアクセスを持つため、GPUで高い性能を得ることができる。しかしながら、GPUはCPUと比較して、メモリ容量が小さいため、CPUと同様の大きさの問題を解くことができない。そこで、本研究では、CPUのホストメモリとCPUのデバイスメモリの局所性を向上させることが可能な、テンポラルブロッキング法を用いることで、GPUのメモリ容量を超える大きさの計算を可能とした。本研究で開発したフレームワークでは、複雑なコーディングは必要とせずに、テンポラルブロッキング法を含む並列計算用のコードを生成できる。フレームワークを用いて開発した気流解析コードでは、TSUBAME2.5において、GPUのメモリ容量の2倍の計算規模においても、通常のメモリ容量の計算の80%程度の実効性能を達成した。
稲垣 厚至*; 神田 学*; Ahmad, N. H.*; 八木 綾子*; 小野寺 直幸; 青木 尊之*
Boundary-Layer Meteorology, 164(2), p.161 - 181, 2017/08
被引用回数:30 パーセンタイル:73.03(Meteorology & Atmospheric Sciences)本研究では、東京都市部に対して、中立安定状態における大気境界層に対する数値解析を行なった。GPUを用いた並列計算を実施することで、19.2km4.8km
1kmの領域に対して2m解像度のラージエディ・シミュレーションが可能となった。大規模計算結果より、境界層上部の乱流統計量や境界層全域におよび特徴的なストリーク構造等の再現が可能であることが確認された。
谷 啓二*; 西尾 敏; 飛田 健次; 筒井 広明*; 三又 秀行*; 飯尾 俊二*; 青木 尊之*
電気学会論文誌,A, 129(9), p.569 - 574, 2009/09
低アスペクト比トカマク炉のリップル損失をモンテカルロコードで解析を行った。リップル損失はアスペクト比(A)に強く依存し、Aが2.5以上の領域ではAの8.8乗に比例することがわかった。このことは低アスペクト比トカマク炉においてはトロイダルコイルの個数を少なくできる可能性を有しており保守を容易にできることを意味する。
高瀬 和之; 吉田 啓之; 玉井 秀定; 小瀬 裕男*; 青木 尊之*; Xu, Z.*
日本機械学会2006年度年次大会講演論文集, Vol.2, p.39 - 40, 2006/09
次世代型軽水炉の燃料集合体を対象にして、開発中であるシミュレーションを主体とした炉心熱設計手法を使って複雑な二相流挙動の予測評価を行った。使用したコードは著者らが開発している界面追跡法を改良したTPFITと商用コードとして世界的に有名であるFLUENTである。TPFITは気液二相流現象の高精度予測が可能であり、FLUENTは非構造格子による複雑流路形状の解析が可能である。本研究によって、燃料棒表面を流れる液膜の挙動,燃料チャンネルを流れる気泡の挙動,スペーサ領域で飛散する液滴の挙動などの詳細を定量的に把握することが可能になり、熱設計精度向上に関しての見通しが得られた。
内海 隆行*; 青木 尊之*; Koga, J. K.; 山極 満
Communications in Computational Physics, 1(2), p.261 - 275, 2006/04
CIP(Constrained Interpolation Profile)基底関数法は、単純な多項式基底関数により、物理量を格子点での値と微分値を用いて近似する。また、関数に対する非線形作用については、微分代数が適用される。さらに、基底関数のスカラー積を導入し、微分方程式の両辺の差(剰余)が基底に直交することを要請することにより、線形及び非線形偏微分方程式は関数値とその空間微分値に対する常微分方程式に簡約化される。この方法は、結合非線形シュレディンガー方程式に対して安定で拡散性が少なく正確な結果を与えることがわかった。
高瀬 和之; 小瀬 裕男*; 吉田 啓之; 叶野 琢磨; 青木 尊之*
日本機械学会関東支部第12期総合講演会講演論文集, p.229 - 230, 2006/03
革新的水冷却炉の炉心では、燃料棒は1mm程度の間隔で三角格子状に稠密に配置され、その間隙を冷却材である水が流れる。水は燃料棒の加熱によって沸騰し、蒸気と水の気液二相流となる。したがって、炉心熱設計を行う場合、相変化を含む気液二相流の挙動を正確に把握する必要がある。本報では、稠密体系に配置された燃料棒まわりに形成される大量の気泡が合体,分裂を繰り返しながら下流へと移行する挙動を数値的に調べた。その結果、微小気泡は時間とともに合体を繰り返しながら大きな気泡に成長するが、稠密流路の体系では気泡は燃料棒間隔の最も狭い領域よりもそれに隣り合う流路面積の広い領域を流れやすい傾向にあることがわかり、シミュレーションによる炉心熱設計手法開発についての見通しを高くできた。
谷 啓二; 飛田 健治; 飯尾 俊二*; 筒井 広明*; 西尾 敏; 青木 尊之*
電気学会論文誌,A, 125(11), p.938 - 942, 2005/11
低アスペクト比トカマク炉VECTORにおけるアルファのリップル損失を軌道追跡モンテカルロコードを用いて検討した。リップル損失は、アスペクト比は小さくになるにしたがって急激に減少する。このため、VERCTORではアルファ粒子は非常によく閉じ込められる。低アスペクト炉における良好なアルファ粒子閉じ込めにより、プラズマ外側端近傍の真空容器に冷却機構を設け、トロイダル磁場コイル内径を30%程度大きくすることで、トロイダル磁場コイル数を6程度まで削減できることがわかった。
小川 慧*; 青木 尊之*; 今井 陽介*; 高瀬 和之
日本機械学会第18回計算力学講演会講演論文集, p.795 - 796, 2005/11
原子炉や化学プラントなどにおいて気液二相流の挙動を詳細にシミュレーションすることは重要である。本研究ではCIP法を用いた気泡流のシミュレーションを行い、計算結果の数値データを可視化することによって高精度な予測を行うことができることを示した。特に可視化手法の一つでありフォトリアリスティックな画像が得られるレイトレーシング法は、気液界面が時間的に変化する気泡流挙動の詳細に理解するうえでたいへん有効であることがわかった。また、可視化に伴うデータ処理には膨大な時間が必要であるが、グリッドコンピューティングシステムの利用によって描画時間を大幅に短縮できる可能性を示すことができた。