著者
村田 健史 上田 裕子 上岡 功治 臼井 英之 岡田 雅樹 大村 善治 松本 紘
出版者
一般社団法人電子情報通信学会
雑誌
電子情報通信学会論文誌. B, 通信 (ISSN:13444697)
巻号頁・発行日
vol.85, no.8, pp.1411-1425, 2002-08-01
被引用文献数
7

現在,プラズマ粒子シミュレーションは,プラズマ観測・応用の様々な分野において用いられている.本研究では,既存のコードをもとに,オブジェクト指向開発技法によりプラズマ粒子シミュレーションコードを設計した.設計したモデルは,プラズマ粒子コードの基本的な機能を有するスケルトンモデルである.オブジェクト指向分析設計方法論によるプラズマ粒子コードの設計については,これまで,詳しい報告がなかった.したがって本論文では,本設計をもとに独自のプラズマ粒子コードを設計できるよう,オブジェクトクラスの洗い出しからデザインまでを詳細に説明した.更に,このデザインに基づき,オブジェクト指向プラズマ粒子シミュレーションコードを,C++言語により実装した.実装コードの基本動作及びコードの再利用性を確認し,更に,既存の手続き型言語によるコードとの,処理速度の比較を行った.その結果,実装したオブジェクト指向コードは,Fortranによる既存コードと比較し,最大3倍程度の計算時間を必要とするが,高い再利用性を有することがわかった.
著者
秋山 隼太 小路 真史 三宅 洋平 大村 善治 中島 浩
出版者
一般社団法人情報処理学会
雑誌
研究報告ハイパフォーマンスコンピューティング(HPC) (ISSN:09196072)
巻号頁・発行日
vol.2010, no.8, pp.1-11, 2010-02-15
参考文献数
6

本論文では,粒子・流体ハイブリッドプラズマシミュレーションの,負荷分散技法 OhHelp を用いた並列化について述べる.すでに OhHelp を適用して良好な結果が得られている全粒子シミュレーションに比べ,ハイブリッドシミュレーションは電磁場の計算負荷が相対的に大きいため,計算と通信のバランスを大幅に見直した実装を行った.特に Cyclic Leapfrog 法による電磁場計算に関する,通信回数削減と計算量増加のトレードオフポイントを見出すことが可能な設計とした.また実用的なシミュレーションに不可欠な,スナップショットやダンプファイルの出力方式も検討し,並列 I/O を用いて並列性能と利便性の両立を図る設計・実装を行った.性能評価の結果,256 プロセスでの実行で 241-456 倍の台数効果が得られること,電磁場計算では計算量増加を抑えることが効果的であること,およびスナップショットと Weak/Strong Scalability との関係が明らかになった.This paper describes a parallel implementation of particle-fluid hybrid plasma simulation with our load balancing method OhHelp. In hybrid simulation, the cost to simulate the progress of electromagnetic field is more significant than that in full-particle simulation whose OhHelp'ed parallelization has already been proved efficient. Thus in this work we revisited the issue of the cost balance between computation and communication, especially for Cyclic Leapfrog method and the trade-off between reducing the number of communications and increasing computational amount. We also designed and implemented parallel-I/O for snapshot and dump, being essential for practical use of our simulator, to reconcile parallel performance and convenience of users. Our evaluation exhibited that the speedup with 256 process is 241- to 456-fold and that suppressing computational cost is the first priority in Cyclic Leapfrog. We also obtained valuable insights about the relationship between weak/strong scalability and snapshot frequency.
著者
臼井 英之 杉崎 由典 冨田 清司 大村 善治 三宅 洋平 青木 正樹
出版者
情報処理学会
雑誌
情報処理学会論文誌コンピューティングシステム(ACS) (ISSN:18827829)
巻号頁・発行日
vol.1, no.2, pp.250-260, 2008-08-21

プラズマ粒子シミュレーションで電磁界成分を更新する際には電流値が必要であり,そのために,個々の粒子の運動量を各空間格子点に集める必要がある.しかし,粒子が空間的にランダムに分布しているため,電流計算の並列演算による高速化は容易ではなく工夫を要する.本論文では,粒子の位置情報を利用して各スレッドに粒子を明示的に割り当てるスレッド並列化アルゴリズムを新しく提案し,OpenMPを用いた実装によりその有効性を検証した.動作検証により,提案手法のCPU台数効果はシミュレーション内の空間格子数の影響を受け,粒子数密度の影響はないことが分かった.特に,各スレッドに割り当てられた空間格子配列がキャッシュに収まりきる程度に細分化される場合,並列台数効果を得やすいことを明らかにした.特に並列台数10前後の場合,その台数効果はスーパリニアとなり,自動並列化コンパイラを用いた電流ルーチン実装に比べて高速になることを明らかにした.また,本提案手法は,各スレッドで全粒子を走査する冗長的な並列化方法であるため,従来アルゴリズムで用いられていた作業領域用配列は不要となり,シミュレーションに必要なメモリ容量を大幅に節約できることを示した.In Particle-In-Cell (PIC) plasma simulations, we calculate the current density to advance the electromagnetic fields. One of the ways to obtain the current density is to gather the velocity moment of each particle to the adjacent grid points. The current calculation is not basically parallelized because the particle positions, which are random in the simulation space, are independent of the array number of current density. In the present paper, we propose a new parallelization method which explicitly associates particles to threads by using OpenMP and evaluate the performance of the proposed method. We clarified that the scalability performance is affected by the number of spatial grid points and is independent of the number of particle per grid. In the proposed method, each thread is in charge of a part of the array of current density divided with the number of thread. When the memory size of the array allocated to each thread becomes small and close to the data cache size of CPU, we found that the scalability performance shows super-linear characteristics and the execution needs less time than the case of using the automatic parallelization compiler. In addition, each thread redundantly scans the particle array to obtain the information of the particle positions for assigning the corresponding particles in charge. Because of this redundant parallelization, we do not have to use work arrays and can save the memory consumed for simulations.
著者
園部 太郎 佐藤 孝宏 奥村 与志弘 広田 勲 津田 冴子 小石 和成 大村 善治
出版者
京都大学東南アジア研究所
雑誌
Kyoto Working Papers on Area Studies: G-COE Series
巻号頁・発行日
vol.36, pp.1-15, 2009-03

A possibility of sustainable biomass power generation from rice fields in Thailand has been investigated from view points of potential power generation, their chemical properties, current status of power generation using rice production residue and the field survey at the rice husk power plant and rice mill of Chainat in Thailand. It is found that rice husk has been mostly used in this region for power generation mostly because of its easiness of collection and transportation as compared with rice straw. With comparison of chemical properties of rice husk and rice straw, it is also suggested that rice husk can be converted to heat relatively easier compared with rice straw. In order to supply the electricity from biomass power plant steadily in Thailand, it can be suggested to increase the fuel variation by utilization of other biomass resources such as rice straw.
著者
松本 紘 BOUGERET Jea ANDERSON Rog 小嶋 浩嗣 GURNETT Dona 村田 健史 笠原 禎也 八木谷 聡 臼井 英之 大村 善治 岡田 敏美 筒井 稔 橋本 弘蔵 長野 勇 木村 磐根 BOUGRET Jean-Louis ANDERSON Roger r. GURNETT D.A. BOUGERET J.L ANDERSON R.R
出版者
京都大学
雑誌
国際学術研究
巻号頁・発行日
1994

平成7年度には、GEOTAIL衛星は、地球から30Re付近の近地球軌道にあり、WIND衛星も主に、昼間側の太陽風の定常観測状態にあった。一方、同年度8月には、ロシアの衛星INTERBALLが、3月には、米国の衛星POLARが打ち上げられ、ISTP衛星による磁気圏の総合観測体制がほぼ整ったといえる。これらの衛星のうち、INTERBALL、POLAR衛星は、打ち上げ後、まもないということで、具体的な共同観測については、来年度に行われる予定であり、本年度は、主に、WIND衛星との共同観測を昨年度までのAKEBONO、Freja,ULYSSES衛星との共同観測に加えて重点的に行った。以下に、交付申請書の調査研究実施計画の項目に従って研究成果を列挙する。1.まず、惑星間衝撃波の観測でGEOTAILとWIND衛星で同時に観測を行った例において、WIND衛星で観測された磁場やプラズマの変化とそのGEOTAILでのある時間遅れでの観測、そしてそれに対応するプラズマ波動の強度の変化について解析を行った。その中には、衝撃波の到来とともにGEOTAILがバウショックを何度もよぎる現象がみられるものがあり、惑星間衝撃波の影響によりバウショックの位置が変化している様子を観測することができた。2.磁気圏昼間側のショック領域全面で発生しているといわれている2fpエミッションの観測をWIND、GEOTAIL両衛星を用いて行い、その発生時間や周波数変化の時間差から、その発生領域がやはりショック全面にあることが確認された。現在その位置的な偏りについても、より多くのデータを集めて解析を行っている。3.GEOTAILによって磁気圏内部で観測された「振幅変調をうけた電子プラズマ波」と同様な波形がWIND衛星によって太陽風中でも観測されていることがわかった。GEOTAILでの観測では、その波動の伝搬方向は外部磁場に対して平行、垂直の両者があることがわかっていたが、現在までのところWINDの方では平行伝搬のみがみつかっている。4.POLAR衛星の打ち上げに伴う共同観測体制を整えるための情報交換をアイオワ大学と行っている。5.POLAR衛星の打ち上げが遅れて本年度の3月になったため、具体的な共同観測は来年に執り行われることになる。6.本研究課題に関連して投稿された論文リストは、本報告書の研究発表欄に列挙する。以上が、交付申請書に書かれていた計画に対応する報告であるが、上述の他に、以下の項目についても共同研究を行った。1.極域で観測されるイオンサイクロトロン波とイオンコニックス分布との相関をAKEBONO衛星とFreja衛星の共同観測で明らかにした。2.極域で観測されるAKRの観測をGEOTAIL、WIND衛星で共同して行い、その観測が衛星の位置によってどのように変化してみられるかの評価をを行い、AKRの伝搬特性についての解析をおこなっている。3.太陽バースト伝搬をGEOTAIL、WIND衛星で同時に観測し、その強度を比較することにより、両者の受信機の較正を行った。