## 平成24年度日本シミュレーション学会学会賞 （敬称略）

### 論文賞（1件）

Maria Savchenko*, Luis A. Diago*, Vladimir Savchenko**, Olga Egorova* and
Ichiro Hagiwara*（*Tokyo Institute of Technology, **Hosei University）

“Mesh Segmentation Using the Platonic Solids”

（日本シミュレーション学会論文誌第3巻第1号）

### 研究賞（2件）

Toshiyuki Hirano（The University of Tokyo）

“Massively Parallel Computation Strategies for Canonical Molecular Orbital Calculations of Proteins”

（第30回日本シミュレーション学会大会(JSST2011) 44）

Teruou Takayama（Yamagata University）

“Time-Domain FEM Analysis of Shielding Current Density in High-Temperature Superconducting Film with Crack: Application to Permanent Magnet Method”

（第30回日本シミュレーション学会大会(JSST2011) 20）

### 奨励賞（3件）

Seiki Saito（Nagoya University）

“Progress of Binary-Collision-Approximation-Based Simulation for Cumulative Structural Changes of Target Materials”

（第30回日本シミュレーション学会大会(JSST2011) 11）

Muhammad Tariq（Waseda University）

“Analyzing Network Reliability in the IEEE 802.15.4 Multi-hop Wireless Sensor Networks”

（第30回日本シミュレーション学会大会(JSST2011) 46）

Ichiro Sakurai（Tohoku University）

“MPTP Injection Yields Oscillating Pattern of Neuronal Activity through the Loop Circuit Formation in the Basal Ganglia”

（第30回日本シミュレーション学会大会(JSST2011) 58）

※所属は論文掲載時

### Outstanding Presentation Award(11件)

1.Speaker：Taku Komura (The University of Edinburgh)

Coauthor：Jake Lever

Session：Regular Session

T i t l e：Controllable Fire using Textured Forces

**Abstract：**Fluid dynamics can produce realistic looking fire effects, which are heavily used in animation and films. However, the parameters of the various underlying physicals equations are not intuitive enough to be controlled easily. As a result, animators face problems when editing the fine details of the fire, especially the turbulence and growth at the fire surface. In this paper, we propose a new method to enable animators to interactively edit such fine details using textured forces. Texture patterns created procedurally or manually sketched by artist are converted into forces that are added at the interface of the burning fuel, where turbulence is known to occur due to fuel igniting. Using our system, the animator can easily control the details of the fire while continuing to produce realistic looking results. Our system can be integrated into existing GPU fluid solvers to run in real-time. As a result, it is applicable for interactive applications such as 3D computer games.

2 .Speaker：Tetsuji Matsuo (Kyoto University)

Coauthor：Takeshi Mifune(Kyoto University)

Session：OS1

T i t l e：3D and 4D Space-Time Grids for Electromagnetic Field Analysis Using Finite Integration Method

**Abstract：**A finite integration method on 3D and 4D space-time grids is studied for
electromagnetic wave computation, where a non-uniform time-step distribution is naturally introduced. A dual grid based on the Hodge duality and the Lorentz metric is introduced to provide a simple constitutive equation for electromagnetic variables. An improved space-time grid shape is proposed to suppress unphysical wave-reflection.

3.Speaker：Daisuke Matsuoka (Japan Agency for Marine-Earth Science and Technology)

Coauthor：Fumiaki Araki (Japan Agency for Marine-Earth Science and Technology)

Shinichiro Kida (Japan Agency for Marine-Earth Science and Technology)

Hideharu Sasaki (Japan Agency for Marine-Earth Science and Technology)

Bunmei Taguchi (Japan Agency for Marine-Earth Science and Technology)

Session：OS2

T i t l e：Visualization for Ocean General Circulation Model via Multi-dimensional Transfer Function and Multivariate Analysis

**Abstract：**With advances in supercomputing technology, ultra-high resolution largescale simulation study has been the mainstream in the field of ocean science. The dataset from such simulation is not only large in size, but also contains a lot of useful information. However, as the geophysical fluid simulation, it is difficult to understand such information hiding in the dataset due to the complexity of ocean dynamics and the number of physical variables. In order to obtain scientific knowledge from large-scale simulation data, it is important to effectively extract and to efficiently express the characteristic feature. The aim of this study is how to solve the above problem. In this study, a multi-dimensional transfer function to emphasis the feature such as ocean currents and vortices from OFES (Ocean general circulation model for the Earth Simulator) dataset is proposed. Furthermore, a multivariate analysis to extract such features is developed. As the results, ocean currents such as the Kuroshio, Oyashio and KE region can efficiently extracted and can effectively be visualized from high-resolution simulation dataset.

4.Speaker：Hiroshi Tamura (Chuo University)

Coauthor：Kazuyuki Miyakita(Niigata University)

Keisuke Nakano(Niigata University)

Masakazu Sengoku(Niigata University)

Shoji Shinoda(Waseda University)

Session：OS3

T i t l e：Generalization of an edge coloring of graphs and channel assignment problem
on multi-hop wireless networks

**Abstract：**Abstract: Multi-hop wireless networks are autonomous systems of mobile nodes connected by wireless links and consist of nodes and wireless links. Multihop wireless networks are useful in various applications when cellular infrastructure is neither available nor realistic or when ubiquitous communication services are required without the presence or use of a fixed infrastructure. For multi-hop wireless networks, we assign channels to communication between two terminals. This is a kind of edge coloring problem in graph theory. In the previous paper, we defined a new edge coloring problem called middle edge coloring and discussed the property of the coloring. In this paper, we generalize the coloring and consider the meaning of the generalization.

5.Speaker：Kazuyasu Sugiyama (The University of Tokyo,RIKEN Computational Science Research Program)

Session：OS5

T i t l e：A Novel Massively Parallel Computing Method for Fluid-Structure Coupling
Analysis

**Abstract：**An efficient and scalable numerical method for massively parallel computing of fluid-structure interaction systems has been developed for biomedical applications. The basic equation set is monolithically formulated in a full Eulerian framework, and solved on a fixed mesh. The discretization is based on a finite-difference volume-of-fluid method, corresponding to a standard numerical approach for simulating incompressible multiphase flows. Further, a novel algorithm based on an artificial compressibility method is proposed. The parameters therein are adaptively determined to guarantee the computed velocity field to be nearly incompressible. The applicability and effectiveness of the proposed method are demonstrated.

6.Speaker：Shinya Maeyama (Japan Atomic Energy Agency)

Coauthor：Tomohiko Watanabe(National Institute for Fusion Science)

Yasuhiro Idomura(Japan Atomic Energy Agency)

Motoki Nakata(Japan Atomic Energy Agency)

Masanori Nunami(National Institute for Fusion Science)

Akihiro Ishizawa(National Institute for Fusion Science)

Session：OS6

T i t l e：Parallelization of a local gyrokinetic Vlasov simulation code for peta-scale
computing

**Abstract：**One of the important phenomena in magnetically-confined fusion plasma is plasma turbulence, which causes particle and heat transport and degrades
plasma confinement. To address multi-scale turbulence including temporal and spatial scales of electrons and ions, we extend our gyrokinetic Vlasov simulation code GKV to run efficiently on peta-scale supercomputers. A key numerical technique is the parallel two-dimensional Fast Fourier Transform (FFT) required for evaluating the nonliear E×B advection term, where masking of the cost of inter-node transpose communications is essential to improve strong scaling. To mask communication costs, overlaps of computations and transpose communications is implemented with the help of the hybrid parallelization of message passing interface and open multiprocessing. Integrated overlaps including FFTs and calculations of the nonlinear E×B term show better scaling than simple overlaps of transpose and FFTs. The masking of communication costs significantly improves strong scaling of the GKV code, and makes substantial speed-up toward multi-scale turbulence simulations.

7.Speaker：Tomoaki Matsumoto (Hosei University)

Session：OS7

T i t l e：Protostellar collapse of magneto-turbulent cloud cores: formation of protoplanetary disks and outflows

**Abstract：**We investigate formation of protoplanetary disks and outflows inside
collapsing turbulent molecular cloud cores by resistive MHD simulations. By using a self-gravitational adaptive mesh refinement (AMR) code, SFUMATO, the collapse of the cloud core is followed in a wide dynamic range; we resolve both a cloud cores scale and a protoplanetary disk scale. A protostar is modeled with a Lagrangian sink particle. The resistivity yields the ohmic dissipation, which is solved by an implicit scheme. We followed collapse of the turbulent molecular cloud cores and accretion onto the protostars (sink particles) up to ∼ 1000 yr after the protostar formation. The simulations show that the strong turbulence and weak magnetic fields promote a relatively rapid growth of a protoplanetary disk around the sink particle. The strong magnetic field models exhibit the cavities in infalling envelopes around the protostars. All the models exhibit protostellar
outflows.

8.Speaker：Yoichiro Yagi (Okayama University of Science)

Coauthor：Yoshinobu Naoshima(Okayama University of Science)

Session：OS8

T i t l e：Large-Scale All-Electron Quantum Chemical Calculation Toward a Sweet-Tasting Protein, Brazzein, and Its Mutants

**Abstract：**Many naturally occurring sweet-tasting proteins have been isolated and are largely expected to be available as a replacement for natural sugars and artificial sweeteners in the near future. We have carried out a large-scale all-electron quantum chemical calculation on des-pGlu brazzein, a sweettasting
protein, and four different mutants by employing a density functional method program, ProteinDF, in order to examine a relationship between the sweetness of protein sweeteners and their electronic properties such as electrostatic potential and frontier molecular orbitals, HOMO and LUMO. ProteinDF all-electron canonical molecular orbital calculations toward des-pGlu brazzein and two different sweet-tatsting mutants, Asp29Lys and Glu41Lys, indicated that a number of amino acidresidues including neutral amino acids, Tyr8 and Tyr51, have positive charge and especially for two mutants the positive charge widely spreads over each mutant. On the other hand, similar all-electron calculations for Arg33Ala and Arg43Ala mutants having a taste like water revealed that the number of amino acid residues with positive charge in these two mutants is less than that in des-pGlu brazzein and Asp29Lys and Glu41Lys mutants. In addition, we found that some positively charged amino acid residues occur in the HOMO and LUMO of des-pGlu brazzein and two sweet-tasting mutants.

9.Speaker：Satoshi Tanaka (Chuo University)

Coauthor：Kazuo Kashiyama (Chuo University)

Akira Kageyama(Kobe University)

Nobuaki Ohno

Session：OS9

T i t l e：Development of Mesh Modification System for Finite Element Simulations
Using CAVE Environments

**Abstract：**This paper presents an interactive mesh modification system for three
dimensional finite element simulaations using CAVE environments. The present system is developed by the VR programming languages, Open GL and CAVE library.Users can check the details of three dimensional mesh structures and can modify the shape of mesh idealization in VR space interactively by using the controller.The mesh refinement method is introduced as a new mesh modification method in this paper. Furthermorehe, the node movement condition which is applicable to the arbitrary domain is also developed for the modificatioan method by node reloacation method. The present system is applied to the three dimensional finite element analysis of over-dentures using dental implants. The results show that the mesh quality is improved by the application of the present
system.

10.Speaker：Boonruang Wangsilabatra (Pathumwan Institute of Technology)

Coauthor：Punyapat Chofabunthi (Pathumwan Institute of Technology)

Session：OS10

T i t l e：Single-Phase to Three-Phase Sinusoidal Waveform Converter for Driving a 3-Phase Induction Motor Based on SPWM Inverter

11.Speaker：Vasily Lubashevskiy (The University of Tokyo)

Coauthor：Taro Kanno (The University of Tokyo)

Kazuo Furuta (The University of Tokyo)

Session：OS11

T i t l e：Emergency Reaction of a Resource Redistribution System to Shock Events Depending on Its Topology

**Abstract：**This work presents the mechanism of emergency reaction of a resource
redistribution system to a large scale disaster which functioning can be regarded as an emergent process. A certain concept of the full recovering of a region after a shock event is constructed and three main stages of this recovering are singled out.The main attention of the work is focused on the initial stage called the
emergency reaction stage. Its main goal is to recover the affected region up to the minimal survival conditions as fast as possible. A simple algorithm implementing this stage is considered. To be specific and to demonstrate its efficiency two limit cases are studied in detail. Namely, the first one is determined by a uniform distribution of relatively small resource stocks over the area under consideration. The second case matches a few large stocks located at main cities distant from one another. The developed mathematical model justifies the results meeting also the general argumentation, which enables us to expect the model to be efficient in
studying complex situations characterized multiscale distributions of resource stocks.

### Student Presentation Award(7件)

1．.Speaker：Fitri Yakub (UTM KL)

Coauthor：Abdul Qadir Adam(IIUM)

Session：Regular Session

T i t l e：ANFIS Based Controller for Variable Speed Induction Generator

**Abstract：**This paper introduces work of Adaptive Neural Fuzzy Inference Systems
(ANFIS) based on voltage control of Induction Generator for wind turbine. Wind power is developing rapidly on both European and global levels over the past 15 years, increased from around 2.5 GW in 1992 to just over 200 GW at the end of last year. Owing to ongoing improvements in turbine efficiency and higher fuel prices, wind power is becoming economically competitive with conventional power production. The optimization of wind turbines has enabled wind energy to become increasingly competitive with conventional energy sources. These developments of wind energy raise a number of challenges to be dealt with now and in the future. The dispersion of wind energy in the grid raises questions about the compatibility of the wind turbine power production with the grid. In this research work, the terminal voltage control of induction generator driven by a variable speed wind turbine is implemented. Wind powered induction generators have an input of unstable which is turbines are equipped with line-connected to the induction generators. At constant frequency, the induction generator operates in a small range of speeds; therefore, it operates with a small range of slips with respect to synchronous speed. Whereas, when variation of the wind turbine speed is wide, induction is no more capable to supply constant voltage and frequency to the grid. In proposed grid connected induction generator, ANFIS based control system is developed that keeps the output ac voltage of induction generator at constant magnitude and frequency even the speed of the wind turbine is varied.

2．.Speaker：Yuki Sato (Hokkaido University)

Coauthor：Hajime Igarashi (Hokkaido University)

Yuta Watanabe (Hokkaido University)

Session：OS1

T i t l e：Finite Element Analysis of Contactless Power Transmission System
Using Meander Line Antennas

**Abstract：**This paper presents finite element analysis and optimization of contactless power transmission systems composed of a pair of meander line antennas. In the electromagnetic field analysis, the circuit equations are coupled with the Maxwell equations. The shape optimization of the meander line antennas is performed using micro genetic algorithm (μGA). The optimized antennas are shown to have the resonant frequency identical with the driving frequency. Moreover, the higher power transmission rate increases with the number of folding of the antennas.

3.Speaker：Seiki Saito (Nagoya University)

Coauthor：Atsushi Ito(National Institute for Fusion Science)

Arimichi Takayama(National Institute for Fusion Science)

Hiroaki Nakamura(Nagoya University and National Institute for Fusion
Science)

Session：OS2

T i t l e：Finite Temperature Effect on Structural Changes of Single Crystalline Graphite under Plasma Irradiation

**Abstract：**Under bombardment of hydrogen plasmas whose thermal energies are
higher than several eV, the structures of single crystalline graphite materials dynamically changes. The channeling effect gradually disappears because of the amorphousization of the target material. In this paper, the effect of the finite temperature of the single crystalline graphite on the channeling effect under hydrogen atom injection is investigated by binarycollision-approximation-based simulation. The result shows that the thermal vibration of target atoms reduces the channeling effects. It is also seen that the thermal vibration is not effective after single crystalline structure changes to amorphous structure by hydrogen injection.

4.Speaker：Katsuya Mouri (National Defense Academy)

Coauthor：Keita Ochiyama (National Defense Academy)

Naobumi Michishita (National Defense Academy)

Yoshihide Yamada (National Defense Academy)

Session：OS3

T i t l e：Bandwidth Characteristics of Folded Structures with Very Small Normal-Mode Helical Antennas

**Abstract：**Recently, normal-mode helical antennas (NMHA) are receiving interests
because of high antenna efficiencies in vely small sizes. However, NMHA have serious defect of very narrow bandwidths. The typical structure that is effective for bandwidth expansion was the folded dipole structure. In previous researches of folded helical antennas, bandwidth expansion were not suceeded. Therefore, methods of expanding bandwidths in folded NMHAs are investigated in this paper. Here, electromagnetic simulations play very important roles to uderstand physical phenomina of antennas and to obtai accurate resilts. First of all, self-resonant structures and bandwidth characteristics of the original NMHA are explained. Next the design methods of the folded NMHA are shown. In order to achive folded antenna mode, to maintain self-resonant conditions at NMHs are shown very
important. For this purpose, the antenna structural tuning chart is clarified.
When antenna structures are well selected, bandwidth expansion of almost
1.7 times are achieved. Also, imput impedance improvements are accomplished.

5.Speaker：Shunsuke Kohno (Kyoto University)

Session：OS5

T i t l e：Molecular Dynamics Simulation of nanobubbles in the bulk and on the
substrate

**Abstract：**We performed molecular dynamics simulations of shrinking nanobubbles in the bulk liquid or on the substrate in binary (liquid and gas) or ternary (liquid, gas and wall) systems. In the simulations, after we controlled the system temperature low and produced nanobubbles, controlling temperature is suddenly raised. After the temperature is raised, nanobubbles begin to shrink immediately. Shrinking process of nanobubbles in the bulk liquid is well described by Rayleigh-Plesset equation. Nanobubbles on the substrate have much longer lifetime than in the bulk liquid. One cause of this long lifetime is the inhomogeneous pressure in the system. We employed not only the wall consists of fixed particles but also the wall consists of thermally vibrating particles in order
to induce the transient heat flux from wall atoms to liquid and gas atoms. Anisotropic motion of gas atoms in the bubble proposed in preceding study is not observed in performed simulations.

6.Speaker：Takashi Masuhara (Kobe University)

Coauthor：Satofumi Souma(Kobe University)

Matsuto Ogawa(Kobe University)

Session：OS8

T i t l e：A First Principles Analysis of Electronic Structure of InGaZnO4
with Modified Becke-Johnson Exchange Potential

**Abstract：**We have investigated the electronic structure of a bulk crystalline
InGaZnO4 by using first principles calculations with the modified Becke-
Johnson (MBJ) potential. We have found that the underestimation of the bandgap and the effective mass of electron, which are inherent to the conventional first principles calculations based on local (spin) density approximation (L(S)DA), is
significantly improved in the present treatment. The obtained values of both the bandgap and the effective mass are found to be comparable to the experimental data. The significant role of the oxygen partial density states (DOS) in the total
DOS in the material is also verified in our calculation.

7.Speaker：Sipaseuth Suvilath (King Mongkut’s Institute of Technology Ladkrabang)

Coauthor：Taworn Benjanarasuth(King Mongkut’s Institute of Technology Ladkrabang)

Noriyuki Komine(Tokai University)

Khamphong Khongsomboun(NUOL)

Session：OS10

T i t l e：Experiments on IMC-Based PID Controller Design for a Two-link SCARA Robot

8.Speaker：Dingding Chao (The University of Tokyo)

Coauthor：Taro Kanno Chao (The University of Tokyo)

Kazuo Furuta Chao (The University of Tokyo)

Session：OS11

T i t l e：Virtual Environment for Participatory Simulation in Tourism