Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
長谷川 雄太; 小野寺 直幸; 井戸村 泰宏
第25回計算工学講演会論文集(CD-ROM), 4 Pages, 2020/06
都市部における風況と汚染物質拡散のアンサンブル計算を行うため、局所細分化格子ボルツマン法を用いたGPUペースのCFDコードを開発した。本コードはPascalまたはVolta世代のGPUアーキテクチャ向けに最適化されており、数km四方の計算領域・数m解像度の格子において実時間で風況解析を行うことができる。開発したコードを用いて、オクラホマシティで行われたフィールド実験JU(Joint Urban)2003の再現計算を行った。計算では、風況は実験値とよく一致し、また、トレーサ濃度のアンサンブルの平均値および最大値がFactor2(計算値と実験値の比が1/2倍
2倍)の条件を満たすことを確認した。
=0 long-range magnetic order in centennialite CaCu
(OD)
Cl
0.6DO; A Spin-
perfect kagome antiferromagnet with 
-
-
飯田 一樹*; 吉田 紘行*; 中尾 朗子*; Jeschke, H. O.*; Iqbal, Y.*; 中島 健次; 河村 聖子; 宗像 孝司*; 稲村 泰弘; 村井 直樹; et al.
Physical Review B, 101(22), p.220408_1 - 220408_6, 2020/06
鉱物センテニアライトCaCu(OD)Cl 0.6DOの結晶構造と磁気構造をシンクロトロンX線回折と中性子回折測定に密度汎関数理論(DFT)と疑似フェルミオン汎関数繰り込み群(PFFRG)の計算を組み合わせることで調べた。CaCu(OD)Cl 0.6DOではCu
イオンは反強磁性と幾何学的に完全なカゴメネットワークを形成する。CuとCaイオン間のサイト間無秩序は見つからなかった。CaCu(OD)Cl 0.6DOは
=7.2K以下で磁気秩序を示し、負のベクトルスピンキラリティーを持つ=0の磁気構造が現れる。0.3Kでの秩序モーメントは0.58(2)
Bに抑えられている。我々のDFT計算では、この系が量子臨界点の近くにあり、--カゴメ反強磁性体の優れた実現であることを示唆している。
青木 尊之*; 長谷川 雄太
自動車技術, 74(4), p.18 - 23, 2020/04
LESに基づくCFD計算を用いて自転車競技の空力解析を行った。単独での走行および24人の集団走行では、算出された抗力は風洞実験と良く一致した。競合する2集団の走行について、集団内の選手の配置を複数検討した。72人の選手の集団走行として、GPUスーパコンピュータで22.3億格子を用いた大規模空力解析を実施した。
稲垣 厚至*; Wangsaputra, Y.*; 神田 学*; Y
cel, M.*; 小野寺 直幸; 青木 尊之*
SOLA (Scientific Online Letters on the Atmosphere) (Internet), 16, p.120 - 124, 2020/00
被引用回数：0 パーセンタイル：100(Meteorology & Atmospheric Sciences)都市境界層を対象とした風況解析により、内層および外層のスケーリングと乱流強度分布の類似性を検討した。計算条件として、現実的な建物形状の上に発達する中立条件を仮定すると共に、計算領域19.2km
4.8km高さ1kmに対して2m格子を設定した。乱流強度分布は計算領域内で局所的に定義できる。内層と外層に対してスケーリング則を当てはめることで、表面形状に関係なく、内層と外層内の乱流強度予測のばらつきを減らすことが可能となる。スケーリングされたプロファイル間のばらつきは、各層のスケーリングパラメーターの不一致に起因するが、長さまたは速度の比率からなる無次元パラメーターを導入することで、それらの類似性を示した。
mesons as a probe of Casimir effect for chiral symmetry breaking石川 力*; 中山 勝政*; 末永 大輝*; 鈴木 渓
Physical Review D, 100(3), p.034016_1 - 034016_14, 2019/08
被引用回数：0 パーセンタイル：100(Astronomy & Astrophysics)本論文では、ゼロ温度および有限温度におけるカイラル対称性の自発的破れに対する有限体積効果を検証するために、中間子がプローブとして役立つことを示した。まず、2+1フレーバーの構成子クォークを含む線形シグマ模型を用いることで、
平均場に対するCasimir効果を解析した。この解析では、反周期境界条件でカイラル対称性が回復し、周期境界条件で対称性の破れが増幅するという結果が得られた。さらに、有限温度・体積平面における平均場の相図を示した。中間子に対しては、カイラルパートナー構造に基づく有効模型を構築した。ここで、中間子質量の体積依存性は、平均場によって与えられる。
中間子は
平均場を含むため、中間子と比べて体積変化に対する応答が鈍いことが判明した。高温・周期境界条件においては、中間子の質量シフトに異常が見られることを発見し、この振る舞いは将来の格子QCDシミュレーションによる検証で役立つことが期待される。さらに、コンパクト化空間次元の数の依存性も調べた。
小野寺 直幸; 井戸村 泰宏; 河村 拓馬; 上澤 伸一郎; 山下 晋; 吉田 啓之
Proceedings of 27th International Conference on Nuclear Engineering (ICONE-27) (Internet), 6 Pages, 2019/05
福島第一原子力発電所の廃炉の方法の一つとして、乾式法が挙げられる。日本原子力研究開発機構(JAEA)では、空冷性能のCFD評価手法としてJUPITERコードを開発している。しかしながら、JUPITERコードにおいて、複雑な原子炉内構造物を捉えた解析を実施するためには、非常の多くの計算資源と計算時間が必要となる。このような問題に対して、本研究ではGPUスーパーコンピュータに適した格子ボルツマン法に基づくCityLBMコードを開発している。CityLBMにてDry methodを模擬したJAEAの実験に対して検証計算を行なった結果、JUPITERコードと同様の結果が得られることが示された。また、同じ解像度および同数の並列数にて計算速度の比較を行った結果、GPUを用いたCityLBM法は、CPUを用いたJUPITERの1/6の計算時間にて解析が行えることが示された。以上の結果より、CityLBMは熱流動解析コードの有効な手法の一つであることが示された。
小野寺 直幸; 井戸村 泰宏; Ali, Y.*; 下川辺 隆史*
Proceedings of 9th Workshop on Latest Advances in Scalable Algorithms for Large-Scale Systems (ScalA 2018) (Internet), p.9 - 16, 2018/11
被引用回数：2 パーセンタイル：12.32計算の高速化に向けて適合細分化格子(AMR)法を適用した格子ボルツマン法(LBM)に対して、通信削減マルチタイムステップ法(CRMT)を提案した。本手法はテンポラルブロッキング法に基づく定式化を行うことで、GPU計算で大きなボトルネックとなる通信回数の削減が可能となる。東京工業大学のTSUBAMEおよび東京大学のReedbushスーパーコンピュータにて性能測定を実施した結果、通信コストが64%に削減され、200GPUまでの弱および強スケーリング結果が改善された。以上の高速化により、2km四方の計算領域に対して1m解像度の風速5msの実時間解析が可能であることが示された。
3(111) symmetrical tilt grain boundary in 
-iron海老原 健一; 鈴土 知明
Modelling and Simulation in Materials Science and Engineering, 26(6), p.065005_1 - 065005_10, 2018/09
被引用回数：1 パーセンタイル：86.78(Materials Science, Multidisciplinary)照射誘起粒界リン偏析の見積もりは原子炉圧力容器鋼の脆化を評価する上で重要な要素であるが、粒界へのリン移動の物理的過程は依然として不明である。このことから、分子動力学を用いて3(111)対称傾角粒界へのリン移動を分子動力学シミュレーションによって評価した。結果として、粒界近傍1nmの範囲で、自己格子間原子が粒界に押し出されることで、鉄-リン混合ダンベルのリン原子と八面体格子間リン原子が置換型原子になることを見出した。空孔-リン複合体も解離し、空孔はリン原子を引きずることなく粒界に吸収された。この結果は、従来考えられている偏析プロセスとは異なることから、それについて新しい視点が必要であると示唆している。
小野寺 直幸; 井戸村 泰宏
Proceedings of 26th International Conference on Nuclear Engineering (ICONE-26) (Internet), 7 Pages, 2018/07
放射性物質の拡散計算は、核セキュリティの観点から非常に重要である。近年、低消費電力および演算性能の高さよりGPUによる科学技術計算が注目されている。本研究では、GPUを用いた格子ボルツマン法による計算コード手法を構築することで、高精度かつ高速な風況解析手法を構築した。東京工業大学のスーパーコンピュータTSUBAMEによる弱スケーリングの性能測定を行った結果、4から144台のGPUに対して良いスケーリングが得られるとともに、CPUとの比較ではノードあたり30倍程度の性能を達成した。産業技術総合研究所の風洞実験に対する検証計算では、AMR法を適用することで全てを細かい格子で計算するのと比較して10%程度の格子点数で解析が可能となるとともに、実験による観測と良く一致している結果が得られた。
高野 和也; 丸山 修平; 羽様 平; 宇佐美 晋
Proceedings of Reactor Physics Paving the Way Towards More Efficient Systems (PHYSOR 2018) (USB Flash Drive), p.1725 - 1735, 2018/04
2010年に実施した、もんじゅ炉心確認試験における炉心反応度の照射依存性について評価した。ゼロ出力で実施した炉心確認試験において、
Puの崩壊に伴う反応度低下以外に、照射量増加に伴う正の反応度増加を確認した。照射依存の反応度増加は炉心確認試験開始から約1ヶ月(10
fissions/cm
)でほぼ飽和する。照射依存の反応度増加は、自己照射損傷に伴い蓄積した格子欠陥が炉心起動中の核分裂片照射により回復したことに起因すると仮定すると、運転前にMOX燃料に蓄積した自己照射損傷に伴う格子欠陥の約47%が回復したことに相当する。
内堀 昭寛; 渡部 晃*; 高田 孝; 大島 宏之
日本機械学会論文集(インターネット), 84(859), p.17-00394_1 - 17-00394_6, 2018/03
Na冷却高速炉の蒸気発生器において伝熱管破損時に形成される隣接伝熱管周りのウェステージ環境を評価するため、Na側で生じる圧縮性多成分多相流及びNa-水化学反応を対象とした機構論的数値解析コードSERAPHIMを開発している。従来のSERAPHIMコードは差分法を用いているが、本研究では、伝熱管の存在する複雑形状領域に対して解析精度を向上することを目的に非構造格子に対応した解析手法を開発し、SERAPHIMコードに組み込んだ。組み込み後SERAPHIMコードにより不足膨張噴流実験の解析を実施した結果、解析結果における圧力分布が実験結果と一致することを確認した。また、Na中へ水蒸気が噴出する現象を対象とした試解析も実施し、妥当と考えられる解析結果を得た。
酒井 卓郎; 飯倉 寛; 山田 尚人*; 佐藤 隆博*; 石井 保行*; 内田 正哉*
QST-M-8; QST Takasaki Annual Report 2016, P. 140, 2018/03
In this work, we report on the application of a UV/EB curable resin, in proton beam writing (PBW) for the fabrication of fine neutron optics devices. The resin is a liquid at room temperature; therefore, it mixes easily with functional materials, such as neutron absorbers. However, 20-30 m thick neutron absorber must be required as the grating materials. So PBW is a promising technique to fabricate the neutron devices. The fabrication process is similar to lithography. The neutron absorber is gadolinium oxide (GdO). Pattern exposures are performed using 3 MeV proton beams approximately 1 m in diameter. The irradiated areas are quickly cured by polymerization. The fabricated gratings were observed using an optical microscope. The microscopic images show that the patterns of the fabricated structures agree well with the original one. In addition, a lateral direction of the sample was observed using the scanning electron microscope (SEM). The thickness of the sample is approximately 35 m. In conclusion, the fabrication of grating structures is successful and the results also proved that the UV/EB curable resin is very promising material for use in PBW micromachining.
小野寺 直幸; 井戸村 泰宏
Lecture Notes in Computer Science 10776, p.128 - 145, 2018/00
被引用回数：6 パーセンタイル：4.68本研究では、局所細分化格子を適用した格子ボルツマン法を開発した。計算コードは、東京工業大学のGPUベースのスーパーコンピュータTSUBAME3.0を用いて開発を行い、最新のPascalアーキテクチャに対して最適化を行なった。1から36ノードを用いた弱スケーリングの性能測定では、NVIDIA TESLA P100を用いたGPU計算がBroadwellによるCPU計算の10倍以上の高速化が達成された。
松川 瞬*; 板倉 賢一*; 早野 明; 鈴木 幸司*
Journal of MMIJ, 133(11), p.256 - 263, 2017/11
LIDAR(Laser Imaging Detection and Ranging)は、岩盤表面を点群の形式で取得することができる。先行研究では、必要なパラメータを手動で設定して、点群から岩盤の不連続面を取得するアルゴリズムが開発されてきた。DiAna(Discontinuity Analysis)アルゴリズムは、点群を格子状に分割して半自動的に岩盤の面を推定するアルゴリズムである。しかし、適切な格子サイズを決定するには熟達した技術が必要である。そこで本研究では、場所によって適切な格子サイズを自動決定するVBS(Variable-Box Segmentation)アルゴリズムを開発した。VBSアルゴリズムは、小さな格子を結合して適切なサイズの格子を作り、面を推定する。VBSアルゴリズムの性能は、DiAnaアルゴリズムと比較して評価した。その際、手動で岩盤表面を推定して作成した参照面との類似度を用いた。比較結果では、VBSアルゴリズムはDiAnaアルゴリズムよりも参照面に類似した面を推定した。よって、VBSアルゴリズムは場所により自動的に適切な格子サイズを決定し、適切に面を推定した。
下川辺 隆史*; 遠藤 敏夫*; 小野寺 直幸; 青木 尊之*
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
被引用回数：5 パーセンタイル：68.61(Meteorology & Atmospheric Sciences)本研究では、東京都市部に対して、中立安定状態における大気境界層に対する数値解析を行なった。GPUを用いた並列計算を実施することで、19.2km4.8km1kmの領域に対して2m解像度のラージエディ・シミュレーションが可能となった。大規模計算結果より、境界層上部の乱流統計量や境界層全域におよび特徴的なストリーク構造等の再現が可能であることが確認された。
CoSbO
伊藤 沙也*; 栗田 伸之*; 田中 秀数*; 河村 聖子; 中島 健次; 伊藤 晋一*; 桑原 慶太郎*; 加倉井 和久*
Nature Communications (Internet), 8, p.235_1 - 235_6, 2017/08
被引用回数：31 パーセンタイル：5.81(Multidisciplinary Sciences)A spin-1/2 triangular-lattice Heisenberg antiferromagnet (TLHAF) is a prototypical frustrated quantum magnet, which exhibits remarkable quantum many-body effects that arise from the synergy between spin frustration and quantum fluctuation. The ground-state properties of a spin-1/2 TLHAF are theoretically well understood. However, magnetic excitations are less well understood and the theoretical consensus is limited. The experimental study of the magnetic excitations in spin-1/2 TLHAFs has also been limited. Here we show the whole picture of magnetic excitations in the spin-1/2 TLHAF BaCoSbO investigated by inelastic neutron scattering. Significantly different from the theory, the excitation spectra have a three-stage energy structure. The lowest first stage is composed of dispersion branches of single-magnon excitations. The second and third stages are dispersive continua accompanied by columnar continuum extending above 10 meV, which is six times larger than the exchange interaction J=1.67 meV. Our results indicate the necessity of a new theoretical framework.
土田 紀之*; Nagahisa, N.*; Harjo, S.
Materials Science & Engineering A, 700, p.631 - 636, 2017/07
被引用回数：1 パーセンタイル：86.86(Nanoscience & Nanotechnology)This study investigated room-temperature creep tests using a low-carbon TRIP steel to clarify the TRIP effect on tensile deformation behavior under constant load. Nominal strain and strain rate increased with an increase in applied stress and those were almost stagnated at holding times of about 
s. The volume fractions of deformation-induced martensite (
) at a given true strain obtained from constant load creep tests were larger than those from tensile tests. From the in situ neutron diffraction experiments during the constant load creep tests, the phase strain of the austenite (
) phase in the creep tests was found to be larger than that in the tensile tests at the same applied stress. This means that the phase strain or true stress of the phase in the TRIP steel was associated with the difference in the volume fraction of between the creep and the tensile tests.
關 義親; 篠原 武尚; Parker, J. D.*; 矢代 航*; 百生 敦*; 加藤 宏祐*; 加藤 秀実*; Sadeghilaridjani, M.*; 大竹 淑恵*; 鬼柳 善明*
Journal of the Physical Society of Japan, 86(4), p.044001_1 - 044001_5, 2017/03
被引用回数：12 パーセンタイル：19.41(Physics, Multidisciplinary)パルス中性子源において効率的な位相イメージングを行うために、複数波長で機能する多色Talbot-Lau干渉計を設計・開発した。J-PARC物質・生命科学実験施設のエネルギー分解型中性子イメージング装置「螺鈿」において、中心波長0.25, 0.50, 0.75nmそれぞれに由来する異なるビジビリティーのモアレ縞を同一セットアップで確認することで、その動作原理の実証を行った。さらに、入射波長の波長分解能を18%から50%まで変化させながら、モアレ縞のビジビリティーの減衰を測定し、高波長分解能が実現できるパルスビームの利点を示した。さらに、中心波長0.5nmを用いて、アルミニウム、鉛、銅のロッドをテストサンプルとしたイメージングも行い、吸収像に加えて、境界がエンハンスされた微分位相像、ビジビリティー像も取得することに成功した。この干渉計では、アナライザー用吸収格子の製作方法として、初めて金属ガラスのインプリント法を採用した。中心波長でのモアレ縞のコントラストは68%と良好であり、従来のガドリニウムの斜め蒸着法を用いた格子と比較すると、形状の精密なコントロールが可能である。
-Iron Using First-Principles-Based Kinetic Monte Carlo海老原 健一; 鈴土 知明; 山口 正剛
Materials Transactions, 58(1), p.26 - 32, 2017/01
被引用回数：3 パーセンタイル：66.31(Materials Science, Multidisciplinary)原子炉圧力容器鋼における照射による粒界リン偏析を数値シミュレーションによって評価することを目指し、原子レベルの計算に基づく拡散レート方程式の開発を行っている。本研究では、第一原理計算で評価した格子間原子対の移動モデルおよび移動障壁エネルギーを用いたキネティックモンテカルロシミュレーションによって、鉄原子とリン原子の混合格子間原子対の拡散係数を評価した。評価した混合格子間原子対の拡散係数は、八面体格子間サイトを移動するリン原子の拡散係数とほぼ同じであり、空孔によるリン輸送の拡散係数よりかなり大きい値であった。さらに、八面体格子間サイトのリン原子を取り入れるように修正した拡散レート方程式に評価した拡散係数を組み入れ、照射誘起粒界リン偏析の数値シミュレーションをしたところ、八面体格子間サイトのリン原子に対する粒界における従来の境界条件が不適切であることが明らかとなった。今後、適切な境界条件を取り入れる必要がある。
