粒子法による流れの数値解析

Size: px
Start display at page:

Download "粒子法による流れの数値解析"

Transcription

1 Numercal Analyss of Flow usng Partcle Method Sech KOSHIZUKA Los Alamos PAF Partcle-and-Force MAC Marker-and- Cell MAC PIC Partcle-n-Cell E-mal: kosh@utnl.jp PIC Los Alamos Harlow 1 PAF SPH Smoothed Partcle Hydrodynamcs 2 SPH MPS Movng Partcle Sem-mplct , 7 8 9

2 MPS r e 2 MPS 2.1 ρ t = 0 Du Dt = 1 ρ P + ν 2 u + f 1 2 MPS MPS 2.2 MPS w r e w (r) = r 1 0 r < r e 3 0 r e r r 3 r e 1 1 r = 0 w MPS r e MPS 3 j n = w ( rj r ) 4 j 1 n φ φ = d φ j φ ( ) ( n 0 rj r r j r 2 w rj r ) 5 j

3 232 j r j r 2 r j r r e j r e j j 2 2 j w d SPH MPS SPH MPS MPS 3 MPS φ = 2d [ (φ ) ( λ n 0 j φ w rj r ) ] 6 j 3 j λ λ = 3 r j r 2 ( w rj r ) j w ( rj r ) 7 j MPS 12

4 233 MPS MPS MPS SMAC Smplfed MAC 2 u r u = u n + t [ ν 2 u n + f n] 8 r = r n + t u 9 n 0 n n 0 n 0 n 0 = n + n 10 u u = t n+1 P 11 ρ n n 0 t + u = P n+1 = ρ n n 0 t 2 n P n+1 11 u n+1 = u + u r n+1 = r + t u n 0 13 n n 0 PIC 14 MPS 3 r = 0 2

5 234 2 n 0 4 PAF SPH SPH SPH 2.4 MPS n <βn 0 16 β P n+1 = 0 β β β 4 MPS MPS 15 DEM Dstnct Element Method 3 MPS MPS 4

6 MPS 5

7 MPS 6 1 MPS ALE Arbtrary Lagrangan-Euleran MAFL Meshless Advecton usng Flowdrectonal Local-grd MPS MPS MAFL ALE MPS MAFL 96 C 110 C 7 1 x y2 MPS MAFL MPS MPS MAFL ALE 9, MPS 8 X 1% 17

8 越塚誠一 図6 円筒水槽の変形と液面振動の相互作用 図 7 サブクール核沸騰における単一気泡の成長と離脱 図8 蒸気爆発における溶融液滴の細粒化過程 237

9 238 1 MPS x y MPS 4 MPS MPS MPS MPS 9 13

10 239 1 Harlow, F. H. : PIC and Its Progeny, Comput. Phys. Comm Monaghan, J. J. : An Introducton to SPH, Comput. Phys. Comm Koshzuka,S.,Tamako,H.&Oka,Y.: APartcle Method for Incompressble Vscous Flow wth Flud Fragmentaton, Comput. Flud Dyn. J Koshzuka, S. & Oka, Y. : Movng-Partcle Sem- Implct Method for Fragmentaton of Incompressble Flud, Nucl. Sc. Eng :, Koshzuka, S., Nobe, A. & Oka, Y. : Numercal Analyss of Breakng Waves usng the Movng Partcle Sem-mplct Method, Int. J. Numer. Meth. Fluds ,,, :, Yoon, H. Y., Koshzuka, S. & Oka, Y. : Drect Calculaton of Bubble Growth, Departure, and Rse n Nucleate Pool Bolng, Int. J. Multphase Flow Koshzuka, S., Ikeda, H. & Oka, Y. : Numercal Analyss of Fragmentaton Mechansms n Vapor Explosons, Nucl. Eng. Des Chkazawa, Y., Koshzuka, S. & Oka, Y. : A Partcle Method for Elastc and Vsco-plastc Structures and Flud-structure Interactons, Comput. Mech Chkazawa, Y., Koshzuka, S. & Oka, Y. : Numercal Analyss of Three-dmensonal Sloshng n an Elastc Cylndrcal Tank usng Movng Partcle Sem-mplct Method, Comput. Flud Dyn. J Ghonem, A. F. & Sherman, F. S. : Grd-free Smulaton of Dffuson usng Random Walk Methods, J. Comput. Phys Koshzuka, S., Ohta, K. & Oka, Y. : Development of a 3-D Calculaton Scheme usng Movng Partcle Sem-mplct Method for Thermal Hydraulcs, Proc. 6th Int. Conf. Nucl. Eng. ICONE-6, San Dego, May 10-15, 1998, ICONE ,, :, Koshzuka, S., Chkazawa, Y. & Oka, Y. : Partcle Method for Flud and Sold Dynamcs, Proc. Frst MIT Conf. on Computatonal Flud and Sold Mechancs, Boston, June 12-15, 2001, Koshzuka, S. & Oka, Y. : Applcaton of Movng Partcle Sem-mplct Method to Nuclear Reactor Safety, Comput. Flud Dynamcs J Lu, J., Koshzuka, S. & Oka, Y. : Investgaton on Energetcs of Ex-vessel Vapor Exploson Based on Spontaneous Nucleaton Fragmentaton, J. Nucl. Sc. Technol Nomura, K., Koshzuka, S., Oka, Y. & Obata, H. : Numercal Analyss of Droplet Breakup Behavor usng Partcle Method, J. Nucl. Sc. Technol

FRP SPH(Smoothed Partcle Hydrodynamcs) (3) MPS(Movng Partcle sem-mplct) (4) MPS (5) (6) (7, 8) (9) (10) (11) Tama and Koshzuka LSMPS(Least Squares Mov

FRP SPH(Smoothed Partcle Hydrodynamcs) (3) MPS(Movng Partcle sem-mplct) (4) MPS (5) (6) (7, 8) (9) (10) (11) Tama and Koshzuka LSMPS(Least Squares Mov Transactons of JSCES, Paper No.20160015 Ansotropc Hgh Vscosty Flud Analyss Usng a Partcle Method for Evaluatng CFRTP Press Moldng Process 1 1 1 2 2 Ryosaku SHINO, Tasuku TAMAI, Sech KOSHIZUKA, Akra MAKI

More information

IPSJ SIG Technical Report Vol.2017-HPC-158 No /3/9 OpenACC MPS 1,a) 1 Moving Particle Semi-implicit (MPS) MPS MPS OpenACC GPU 2 4 GPU NVIDIA K2

IPSJ SIG Technical Report Vol.2017-HPC-158 No /3/9 OpenACC MPS 1,a) 1 Moving Particle Semi-implicit (MPS) MPS MPS OpenACC GPU 2 4 GPU NVIDIA K2 OpenACC MPS 1,a) 1 Movng Partcle Sem-mplct (MPS) MPS MPS OpenACC GPU 2 4 GPU NVIDIA K20c GTX1080 P100(PCIe) P100(NVlnk) 5 OpenACC 3.5 3 Fortran 29.0 74.5 GPU 1. MPS [1] 1 MPS MPS CUDA GPU [2] [3] [4] OpenACC

More information

混相流 Vol.31 No.4

混相流 Vol.31 No.4 論 文 Movng Partcle Sem-mplct 法による ノズル内キャビテーションと液体噴流のシミュレーション * Numercal Smulaton of Cavtaton n a Nozzle and Lqud Jet Usng Movng Partcle Sem-mplct Method 関 根 章 裕 ** SEKINE Akhro 吉 村 一 樹 *** YOSHIMURA Kazuk

More information

untitled

untitled - 37 - - 3 - (a) (b) 1) 15-1 1) LIQCAOka 199Oka 1999 ),3) ) -1-39 - 1) a) b) i) 1) 1 FEM Zhang ) 1 1) - 35 - FEM 9 1 3 ii) () 1 Dr=9% Dr=35% Tatsuoka 19Fukushima and Tatsuoka19 5),) Dr=35% Dr=35% Dr=3%1kPa

More information

u t x p u t q () p u C q (3) ls x x I I u 0 (4) t x p 0 q f (5) x xs q (6) * x xs 計 算 を,CIP 法 に 代 わりにTHINC/WLIC 法 を 用 いて 実 施 する. 固 相 にLagrange 粒 子 を 配

u t x p u t q () p u C q (3) ls x x I I u 0 (4) t x p 0 q f (5) x xs q (6) * x xs 計 算 を,CIP 法 に 代 わりにTHINC/WLIC 法 を 用 いて 実 施 する. 固 相 にLagrange 粒 子 を 配 3 次 元 固 気 液 多 相 乱 流 数 値 モデルDOLPHIN-3D への 高 精 度 界 面 捕 獲 法 THINC/WLIC 法 の 導 入 川 崎 浩 司 1 松 野 哲 弥 1 正 会 員 株 式 会 社 ハイドロソフト 技 術 研 究 所 ( 530-616 大 阪 市 北 区 中 之 島 3-3-3 中 之 島 ダイビル6F) E-mal: kawasak@hydro-soft.co.jp

More information

平成18年度

平成18年度 平成 18 年度 修士論文 粒子法によるキャビティ流れの数値シミュレーション 高知工科大学大学院工学研究科基盤工学専攻知能機械システム工学コース知能流体力学研究室 矢野敦大 目次 第 1 章 緒言 -1-1. 1 はじめに -1-1. 差分法と粒子法の比較 -- 1. 3 研究目的 -3- 第 章 基礎方程式 -4-. 1 オイラーの方法とラグランジュの方法 -4-. 基礎方程式 -5-. 3 無次元化

More information

7 OpenFOAM 6) OpenFOAM (Fujitsu PRIMERGY BX9, TFLOPS) Fluent 8) ( ) 9, 1) 11 13) OpenFOAM - realizable k-ε 1) Launder-Gibson 15) OpenFOAM 1.6 CFD ( )

7 OpenFOAM 6) OpenFOAM (Fujitsu PRIMERGY BX9, TFLOPS) Fluent 8) ( ) 9, 1) 11 13) OpenFOAM - realizable k-ε 1) Launder-Gibson 15) OpenFOAM 1.6 CFD ( ) 71 特集 オープンソースの大きな流れ Nonlinear Sloshing Analysis in a Three-dimensional Rectangular Pool Ken UZAWA, The Center for Computational Sciences and E-systems, Japan Atomic Energy Agency 1 1.1 ( ) (RIST) (ORNL/RSICC)

More information

シリコン結晶化過程の分子動力学

シリコン結晶化過程の分子動力学 1-65 13 16 96154 1 4 1.1 5 1.1.1 5 1.1. 5 1.1.3 5 1.1.4 CVD 7 1. 8 1..1 8 1.. SPE 10 1.3 11 1.1 13. 14..1 Tesoff 14.. Lennad-Jones 17.3 18.4 19.5 0.6 Langevn.7 3.7.1 SPE 3.7. 4.7.3 6 3 SPE 8 3.1 9 3. [111]

More information

IEEE HDD RAID MPI MPU/CPU GPGPU GPU cm I m cm /g I I n/ cm 2 s X n/ cm s cm g/cm

IEEE HDD RAID MPI MPU/CPU GPGPU GPU cm I m cm /g I I n/ cm 2 s X n/ cm s cm g/cm Neutron Visual Sensing Techniques Making Good Use of Computer Science J-PARC CT CT-PET TB IEEE HDD RAID MPI MPU/CPU GPGPU GPU cm I m cm /g I I n/ cm 2 s X n/ cm s cm g/cm cm cm barn cm thn/ cm s n/ cm

More information

Web Two-phase Flow Analyses Using Interface Volume Tracking Tomoaki Kunugi Kyoto University 1) 2) 3)

Web Two-phase Flow Analyses Using Interface Volume Tracking Tomoaki Kunugi Kyoto University   1) 2) 3) Web 11 3 2003 8 Two-phase Flow Analyses Using Interface Volume Tracking Tomoaki Kunugi Kyoto University E-mail: kunugi@nucleng.kyoto-u.ac.jp 1) 2) 3) Lagrangian 4) MAC(Marker and Cell) 5) (VOF:Volume of

More information

program_japanese

program_japanese Eng Eng Eng - 40 - Eng Eng - 41 - Eng Eng - 42 - Eng - 43 - Eng - 44 - Eng Eng - 45 - Eng - 46 - - 47 - Eng - 48 - - 49 - - 50 - - 51 - - 52 - - 53 - - 54 - - 55 - - 56 - - 57 - - 58 - - 59 - - 60 - -

More information

Explicit MPS Algorithm for Large-scale Free Surface Flow Analysis

Explicit MPS Algorithm for Large-scale Free Surface Flow Analysis 先駆的科学計算に関するフォーラム 014 大規模分散メモリ並列 MS 法の開発と 市街地浸水津波解析 開催日 014 年 8 月 5 日 ( 火 )-8 月 6 日 ( 水 ) 会場福岡市博多区リファレンス駅東ビル 5 階会議室 V- 室谷浩平 ( 東京大学 ) 目次 MS 陽解法 並列計算 津波解析 MS 陽解法とMS 半陰解法の比較 目次 MS 陽解法 並列計算 津波解析 MS 陽解法とMS 半陰解法の比較

More information

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical calculation method of the gradient as a differential

More information

数値計算:有限要素法

数値計算:有限要素法 ( ) 1 / 61 1 2 3 4 ( ) 2 / 61 ( ) 3 / 61 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61 u(x) 0 L x ( ) 6 / 61 P(0) P(L) f d dx ( EA du dx

More information

IPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of

IPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of GPU 1,a) 1,b) 1,c) 1,d) GPU 1 GPU Structure Of Array Array Of Structure 1. MPS(Moving Particle Semi-Implicit) [1] SPH(Smoothed Particle Hydrodynamics) [] DEM(Distinct Element Method)[] [] 1 Tokyo Institute

More information

2005 2006.2.22-1 - 1 Fig. 1 2005 2006.2.22-2 - Element-Free Galerkin Method (EFGM) Meshless Local Petrov-Galerkin Method (MLPGM) 2005 2006.2.22-3 - 2 MLS u h (x) 1 p T (x) = [1, x, y]. (1) φ(x) 0.5 φ(x)

More information

非線形長波モデルと流体粒子法による津波シミュレータの開発 I_ m ρ v p h g a b a 2h b r ab a b Fang W r ab h 5 Wendland 1995 q= r ab /h a d W r ab h

非線形長波モデルと流体粒子法による津波シミュレータの開発 I_ m ρ v p h g a b a 2h b r ab a b Fang W r ab h 5 Wendland 1995 q= r ab /h a d W r ab h 土木学会論文集 B2( 海岸工学 ) Vol. 70, No. 2, 2014, I_016-I_020 非線形長波モデルと流体粒子法による津波シミュレータの開発 Development of a Tsunami Simulator Integrating the Smoothed-Particle Hydrodynamics Method and the Nonlinear Shallow Water

More information

Fig. Division of unbounded domain into closed interior domain and its eterior domain. Zienkiewicz [5, 6] Burnett [7, 8] [3] The conjugated Ast

Fig. Division of unbounded domain into closed interior domain and its eterior domain. Zienkiewicz [5, 6] Burnett [7, 8] [3] The conjugated Ast 7 6 pp. 635 643 635 43..Rz; 43.4.Rj * 3 3 Unbounded problems, Finite element method, Infinite element, Hybrid variational principle, Fourier series. Boundary Element Method: BEM BEM Finite Element Method:

More information

Undulator.dvi

Undulator.dvi X X 1 1 2 Free Electron Laser: FEL 2.1 2 2 3 SACLA 4 SACLA [1]-[6] [7] 1: S N λ [9] XFEL OHO 13 X [8] 2 2.1 2(a) (c) z y y (a) S N 90 λ u 4 [10, 11] Halbach (b) 2: (a) (b) (c) (c) 1 2 [11] B y = n=1 B

More information

untitled

untitled @=@tur (substantal dervatve) Euler Naver-Stkes d @ dt @t ur t u 2 (streamlne).(a) x ( ) A (prperty ) ' B 'ur' ' t (path lne) ru =0 d dt ' = @ @t ' r(u') (cnservatve frm) (nncnservatve frm) (gradent frm)

More information

タイトル

タイトル Flud Flow Smulton wth Cellulr Automt 00N2100008J 2002 225 1. cellulr utomton n t+ 1 t t = f( r, L, + r (1 t t+ 1 f t r t +1 Prllel Vrtul Mchne Messge-Pssng Interfce 1 2. 2. 1 t t 0 t = 1 = 1 = t = 2 2

More information

QW-3414

QW-3414 MA1312-C P 1 2 3 A E L D E D A A E D A D D D D D E A C A C E D A A A C A C A C E E E D D D A C A C A A A A C A C A C E E C C E D D C C C E C E C C E C C C E D A C A C A C E L B B

More information

土壌環境行政の最新動向(環境省 水・大気環境局土壌環境課)

土壌環境行政の最新動向(環境省 水・大気環境局土壌環境課) 201022 1 18801970 19101970 19201960 1970-2 1975 1980 1986 1991 1994 3 1999 20022009 4 5 () () () () ( ( ) () 6 7 Ex Ex Ex 8 25 9 10 11 16619 123 12 13 14 5 18() 15 187 1811 16 17 3,000 2241 18 19 ( 50

More information

syuryoku

syuryoku 248 24622 24 P.5 EX P.212 2 P271 5. P.534 P.690 P.690 P.690 P.690 P.691 P.691 P.691 P.702 P.702 P.702 P.702 1S 30% 3 1S 3% 1S 30% 3 1S 3% P.702 P.702 P.702 P.702 45 60 P.702 P.702 P.704 H17.12.22 H22.4.1

More information

,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), Mellor and Yamada 1974) 4 2 3, 2 4,

,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), Mellor and Yamada 1974) 4 2 3, 2 4, Mellor and Yamada1974) The Turbulence Closure Model of Mellor and Yamada 1974) Kitamori Taichi 2004/01/30 ,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), 4 1 4 Mellor and Yamada 1974) 4 2 3, 2

More information

H. Kuninaka and H. Hayakawa Phys. Rev. Lett. 93, 543 (4) M. Y. Louge & E. Adams Phys. Rev. E 65, 33 () Y.Tanaka et. al. Europhys. Lett. 63, 46 (3) K. Okumura et. al. Europhys. Lett. 6, 37 (3) R +X h R

More information

-37-

-37- ) ) ) ) ) ) (a) (b) -36- -37- (a) ) ) (b) ) -3 3 LIQCA ) ) a) Oka 999 3 LIQCA -38- u-p formulationoka 994 Newmark 3 FEM -39- b) (i) 3 FEM Zhang 9 3 (ii) () Tatsuoka 98Dr=9% 89% RL=.6DA=7.5% 998 999 Dr=5%

More information

Keywords: corotational method, Rigid-Bodies-Spring model, accuracy, geometrical nonlinearity

Keywords: corotational method, Rigid-Bodies-Spring model, accuracy, geometrical nonlinearity Keywords: corotational method, Rigid-Bodies-Spring model, accuracy, geometrical nonlinearity vo2=voi+(sina1+sina2)l/2+{f1(sin2a1+sin2a2) +Fz1(sina1cosai+sina2cosa2)}l/(2EA) woe=w01-l+(cosai+cosa2)l/2+{fy1(sinaicosai

More information

163 KdV KP Lax pair L, B L L L 1/2 W 1 LW = ( / x W t 1, t 2, t 3, ψ t n ψ/ t n = B nψ (KdV B n = L n/2 KP B n = L n KdV KP Lax W Lax τ KP L ψ τ τ Cha

163 KdV KP Lax pair L, B L L L 1/2 W 1 LW = ( / x W t 1, t 2, t 3, ψ t n ψ/ t n = B nψ (KdV B n = L n/2 KP B n = L n KdV KP Lax W Lax τ KP L ψ τ τ Cha 63 KdV KP Lax pair L, B L L L / W LW / x W t, t, t 3, ψ t n / B nψ KdV B n L n/ KP B n L n KdV KP Lax W Lax τ KP L ψ τ τ Chapter 7 An Introduction to the Sato Theory Masayui OIKAWA, Faculty of Engneering,

More information

τ p ω πτ p ω π τ p (t) = 2 2 t 2 exp(i t)exp 8 2 S(,t) = s( ) (t )d d 2 E x dz 2 = 2 E x z E x = E 0 e z, = + j = 1 2 0 tan = 0, v = c r 10 11 Horn Circulator Net Work Analyzer t H = E t E = H E t B =

More information

1.indd

1.indd 14 15 6 april 6 16 17 18 april 18 april 18 19 28 april 20 21 28 april 22 23 21 may 24 25 10 june 21 may 26 27 10 june 28 29 12 june 30 31 12 25 june 32 33 25 34 35 4 july 37 36 38 39 25 july 25 july 40

More information

mr0805.indd

mr0805.indd 62 MAY 2008 MAY 2008 63 64 MAY 2008 MAY 2008 65 66 MAY 2008 MAY 2008 67 68 MAY 2008 MAY 2008 69 70 MAY 2008 MAY 2008 71 72 MAY 2008 MAY 2008 73 74 MAY 2008 MAY 2008 75 76 MAY 2008 MAY 2008 77 78 MAY 2008

More information

野岩鉄道の旅

野岩鉄道の旅 29th 5:13 5:34 5:56 6:00 6:12 6:20 6:21 6:25 6:29 6:31 6:34 6:38 6:40 6:45 6:52 6:56 7:01 7:07 7:11 7:32 7:34 7:50 7:58 8:03 8:17 8:36 8:44 5:50 5:54 6:15 6:38 6:39 6:51 6:59 6:59 7:03 7:08 7:08 7:11 7:15

More information

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x

Trapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )

More information

gr09.dvi

gr09.dvi .1, θ, ϕ d = A, t dt + B, t dtd + C, t d + D, t dθ +in θdϕ.1.1 t { = f1,t t = f,t { D, t = B, t =.1. t A, tdt e φ,t dt, C, td e λ,t d.1.3,t, t d = e φ,t dt + e λ,t d + dθ +in θdϕ.1.4 { = f1,t t = f,t {

More information

1 capillary-gravity wave 1) 2) 3, 4, 5) (RTI) RMI t 6, 7, 8, 9) RMI RTI RMI RTI, RMI 10, 11) 12, 13, 14, 15) RMI 11) RTI RTI y = η(x, t) η t

1 capillary-gravity wave 1) 2) 3, 4, 5) (RTI) RMI t 6, 7, 8, 9) RMI RTI RMI RTI, RMI 10, 11) 12, 13, 14, 15) RMI 11) RTI RTI y = η(x, t) η t 1 capillary-gravity wave 1) ) 3, 4, 5) (RTI) RMI t 6, 7, 8, 9) RMI RTI RMI RTI, RMI 10, 11) 1, 13, 14, 15) RMI 11) RTI RTI 3 3 4 y = η(x, t) η t + ϕ 1 η x x = ϕ 1 y, η t + ϕ η x x = ϕ y. (1) 1 ϕ i i (i

More information

ohpr.dvi

ohpr.dvi 2003/12/04 TASK PAF A. Fukuyama et al., Comp. Phys. Rep. 4(1986) 137 A. Fukuyama et al., Nucl. Fusion 26(1986) 151 TASK/WM MHD ψ θ ϕ ψ θ e 1 = ψ, e 2 = θ, e 3 = ϕ ϕ E = E 1 e 1 + E 2 e 2 + E 3 e 3 J :

More information

untitled

untitled Stacking sequence optimization of composite wing using fractal branch and bound method Orbiting Plane : HOPE-X (JAXA) Fractal Branch and Bound Method (FBBM) Fractal structure of design space 5 y V a 9º

More information

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2 2005 9/8-11 2 2.2 ( 2-5) γ ( ) γ cos θ 2πr πρhr 2 g h = 2γ cos θ ρgr (2.1) γ = ρgrh (2.2) 2 cos θ θ cos θ = 1 (2.2) γ = 1 ρgrh (2.) 2 2. p p ρgh p ( ) p p = p ρgh (2.) h p p = 2γ r 1 1 (Berry,1975) 2-6

More information

3 風工学シンポ の場合や 7),),9), 計算領域外側の境界条件が不明確 9),1) であったりする 本研究では, 並進移動ならびに回転移動する物体に固定された非慣性系座標による LES の基礎方程式を誘導し, 橋梁断面の空力特性検討を想定した物体の風直交方向ならびに捩れの空力不安定振動を行った

3 風工学シンポ の場合や 7),),9), 計算領域外側の境界条件が不明確 9),1) であったりする 本研究では, 並進移動ならびに回転移動する物体に固定された非慣性系座標による LES の基礎方程式を誘導し, 橋梁断面の空力特性検討を想定した物体の風直交方向ならびに捩れの空力不安定振動を行った 3 風工学シンポ第 3 回風工学シンポジウム (1) 非慣性系座標の LES による一様流中角柱の空力不安定振動の計算 LARGE EDDY SIMULAION OF AEROELASIC INSABILIY OF RECANGULAR CYLINDER IN UNIFORM WIND USING NON-INERIAL COORDINAE SYSEM 野田博 1) 中山昭彦 ) Hrosh NODA

More information

Vol.1( ) No JASCOME Trefftz ( ) SIMULATION OF SLOSHING PHENOMENON BY INDIRECT TREFFTZ METHOD (EXTENSION OF SIMULATION SCHEME) 1), 2),

Vol.1( ) No JASCOME Trefftz ( ) SIMULATION OF SLOSHING PHENOMENON BY INDIRECT TREFFTZ METHOD (EXTENSION OF SIMULATION SCHEME) 1), 2), Vol.1( 2001 7 ) No.01-070611 JASCOME Trefftz ( ) SIMULATION OF SLOSHING PHENOMENON BY INDIRECT TREFFTZ METHOD (EXTENSION OF SIMULATION SCHEME) 1), 2), 3) 4), Yoichi IKEDA, Jun ichi KATSURAGAWA, Eisuke

More information

wakate2005-tobari-ver2.ppt

wakate2005-tobari-ver2.ppt 8 005 3 17 Outlne 1. Introducton. 3. 4. Introducton Image of MUSES-C on engne JAXA b = 4ε 9 q m V d 3 a large-scale) MHD MPD chokng VASIMR Super-Alfvénc flow plasma detachment HITOP(HIgh densty TOhoku

More information

noted by Kazuma MATSUDA (force) (mechanics) (statics) (dynamics) O F OA O (point of application) (line of action) (vector) F F F

noted by Kazuma MATSUDA (force) (mechanics) (statics) (dynamics) O F OA O (point of application) (line of action) (vector) F F F noted by Kazuma MATSUDA 2016 5 29 1. (force) (mechancs) (statcs) (dynamcs) 1.1 1 1 O F OA O (pont of applcaton) (lne of acton) (vector) F F F F 1.2 ( ) 1 (kg-wt) (Internatonal System of Unt, SI) kg SI

More information

ウェーブレットによる経済分析

ウェーブレットによる経済分析 E-malmasakazu.nada@boj.or.jp E-malkouchrou.kamada@boj.or.jp wavelet J.J. Morlet D. D. Gabor uncertanty prncple Conway and Frame Schlecher wave let DWT: dscrete wavelet transformcwt: contnuous wavelet

More information

> σ, σ j, j σ j, σ j j σ σ j σ j (t) = σ (t ) σ j (t) = σ () j(t ) n j σ, σ j R lm σ = σ j, j V (8) t σ R σ d R lm σ = σ d V (9) t Fg.. Communcaton ln

> σ, σ j, j σ j, σ j j σ σ j σ j (t) = σ (t ) σ j (t) = σ () j(t ) n j σ, σ j R lm σ = σ j, j V (8) t σ R σ d R lm σ = σ d V (9) t Fg.. Communcaton ln IIC-- Dstrbuted Cooperatve Atttude Control for Multple Rgd Bodes wth Communcaton Delay Yoshhro achbana, oru Namerkawa (Keo Unversty) Abstract hs paper descrbes dstrbuted cooperatve atttude consensus and

More information

Evidence for jet structure in hadron product by e+e-

Evidence for jet structure in hadron product by e+e- G. Hanson et al. Phys. Rev. Lett. 5 (1975) 1609 Physcs Colloquum July 7th, 008 Evdence for Jet Structure n Hadron Producton by e + e - Annhlaton Contents: 1. Introducton. Exerment at SLAC. Analyss 4. Results

More information

橡jttc2.PDF

橡jttc2.PDF 1 ( ) 1 GA GA GA MOGA (Multple-Objectve Genetc Algorthm) GA GA GA MOGA GA GA MOGA GA GA 3.1MOGA ( ) x x j f = f, f, 1 2 L, f q x x j x j f ( x ) f ( x ) f ( x ) f ( x ) L f ( x ) f ( x ) ( ) ( ) 1 1 j

More information

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting

More information

From Evans Application Notes

From Evans Application Notes 3 From Evans Application Notes http://www.eaglabs.com From Evans Application Notes http://www.eaglabs.com XPS AES ISS SSIMS ATR-IR 1-10keV µ 1 V() r = kx 2 = 2π µν x mm 1 2 µ= m + m 1 2 1 k ν = OSC 2

More information

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P 9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)

More information

A

A A04-164 2008 2 13 1 4 1.1.......................................... 4 1.2..................................... 4 1.3..................................... 4 1.4..................................... 5 2

More information

光学

光学 Control of Refractive Indices of Optical Polymers: Theoretical Prediction and Molecular Design Method Shinji ANDO The fundamental theory and a method to predict the refractive indices and dispersions of

More information

arxiv: v1(astro-ph.co)

arxiv: v1(astro-ph.co) arxiv:1311.0281v1(astro-ph.co) R µν 1 2 Rg µν + Λg µν = 8πG c 4 T µν Λ f(r) R f(r) Galileon φ(t) Massive Gravity etc... Action S = d 4 x g (L GG + L m ) L GG = K(φ,X) G 3 (φ,x)φ + G 4 (φ,x)r + G 4X (φ)

More information

JFE.dvi

JFE.dvi ,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho

More information

untitled

untitled 1 17 () BAC9ABC6ACB3 1 tan 6 = 3, cos 6 = AB=1 BC=2, AC= 3 2 A BC D 2 BDBD=BA 1 2 ABD BADBDA ABC6 BAD = (18 6 ) / 2 = 6 θ = 18 BAD = 12 () AD AD=BADCAD9 ABD ACD A 1 1 1 1 dsinαsinα = d 3 sin β 3 sin β

More information

o ξς ξ, ς u, w u ( ξ, ς ), w ( ξ, ς ) { α } u ( ξ, ς ) = { ξ ξς ς}{ α } T α α3 α 4 { }{ } T ( ξ ς = ξ ξς ς α α α α w, ) 5 6 7 8 = { }{ α α α α } T = { b}{ α α α α } T u 3 4 u 3 a ab b α α α α 3 4 u 3 4

More information

[Ver. 0.2] 1 2 3 4 5 6 7 1 1.1 1.2 1.3 1.4 1.5 1 1.1 1 1.2 1. (elasticity) 2. (plasticity) 3. (strength) 4. 5. (toughness) 6. 1 1.2 1. (elasticity) } 1 1.2 2. (plasticity), 1 1.2 3. (strength) a < b F

More information

Microsoft PowerPoint - OpenSourceNov07m.pptx

Microsoft PowerPoint - OpenSourceNov07m.pptx PPU PU による圧縮性 非圧縮性流体計算の驚異的な高速化 東京工業大学学術国際情報センター 青木尊之 (eneral-purpose raphcs Processng Unt) PU を画像処理以外の目的で使う PU の魅力 高性能 : ハイエンド PU はピーク TFlops 超 低価格 : ハイエンドでも数万円 (4 万円 ) 手軽さ : 事務所の PC にも装着できる 単一 PU 高消費電力

More information

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' =

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' = y x = α + β + ε =,, ε V( ε) = E( ε ) = σ α $ $ β w ( 0) σ = w σ σ y α x ε = + β + w w w w ε / w ( w y x α β ) = α$ $ W = yw βwxw $β = W ( W) ( W)( W) w x x w x x y y = = x W y W x y x y xw = y W = w w

More information

IPSJ SIG Technical Report Vol.2013-MUS-100 No /9/1 1,a) [1 9] 1 3 MIDI MIDI (Interonset interval) 1 NTT NTT Communication Science Labo

IPSJ SIG Technical Report Vol.2013-MUS-100 No /9/1 1,a) [1 9] 1 3 MIDI MIDI (Interonset interval) 1 NTT NTT Communication Science Labo 1,a) 1 1 1. [1 9] 1 3 MIDI MIDI (Interonset nterval) 1 N N Communcaton Scence Laboratores, Atsug, Kanagawa 43 0198, Japan he Insttute of Statstcal Mathematcs, achkawa, okyo 190 856 a) ohsh.yasunor@lab.ntt.co.jp

More information

工学的な設計のための流れと熱の数値シミュレーション

工学的な設計のための流れと熱の数値シミュレーション 247 Introduction of Computational Simulation Methods of Flow and Heat Transfer for Engineering Design Minoru SHIRAZAKI Masako IWATA Ryutaro HIMENO PC CAD CAD 248 Voxel CAD Navier-Stokes v 1 + ( v ) v =

More information

201711grade1ouyou.pdf

201711grade1ouyou.pdf 2017 11 26 1 2 52 3 12 13 22 23 32 33 42 3 5 3 4 90 5 6 A 1 2 Web Web 3 4 1 2... 5 6 7 7 44 8 9 1 2 3 1 p p >2 2 A 1 2 0.6 0.4 0.52... (a) 0.6 0.4...... B 1 2 0.8-0.2 0.52..... (b) 0.6 0.52.... 1 A B 2

More information

C:/大宮司先生関連/大宮司先生原稿作業用1/圧縮性流れの解法3.dvi

C:/大宮司先生関連/大宮司先生原稿作業用1/圧縮性流れの解法3.dvi 8 aldwn-lomax hen 8. +ru = t (8.a) u+r(uu+) =r +f t (8.b) e +rhu = r ( u);r q+f u t (8.c) u e = (+u 2 =2) H = h+u 2 =2=(e+)= q f (8.) (comressble Naver-Stokes equatons) T h = += = RT d = c v dt dh = c

More information

CD納品用.indd

CD納品用.indd 1. 1.1 1.2 1.3 1.4 2. 3. 3.1 3.2 3.3 4. 4.1 4.2 6. Journal of the Institute of Nuclear Safety System Journal of the Institute of Nuclear Safety System Journal of the Institute of Nuclear Safety System

More information

, 3, STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO

, 3, STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO , 3, 2012 9 STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO 1 950-2181 2 8050 2 950-2181 2 8050 Numerical computation of river flows

More information

/ Christopher Essex Radiation and the Violation of Bilinearity in the Thermodynamics of Irreversible Processes, Planet.Space Sci.32 (1984) 1035 Radiat

/ Christopher Essex Radiation and the Violation of Bilinearity in the Thermodynamics of Irreversible Processes, Planet.Space Sci.32 (1984) 1035 Radiat / Christopher Essex Radiation and the Violation of Bilinearity in the Thermodynamics of Irreversible Processes, Planet.Space Sci.32 (1984) 1035 Radiation and the Continuing Failure of the Bilinear Formalism,

More information

Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue

Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Date 2017 URL http://hdl.handle.net/2433/229150 Right

More information

Bifurcation Structure of a Novel Car-following Model with Relative Velocity Effect 1 Akiyasu Tomoeda 2 Tomoyuki Miyaji 3 2011 AUTO Hopf Kota Ikeda 1 [1, 2] [3, 4, 5]. [4] 1 1 / / JST CREST 2 / 3 / 54 (,

More information

橡超弦理論はブラックホールの謎を解けるか?

橡超弦理論はブラックホールの謎を解けるか? 1999 3 (Can String Theory Solve the Puzzles of Black Holes?) 305-0801 1-1 makoto.natsuume@kek.jp D-brane 1 Schwarzschild 60 80 2 [1] 1 1 1 2 2 [2] 25 2.2 2 2.1 [7,8] Schwarzschild 2GM/c 2 Schwarzschild

More information

第5章 偏微分方程式の境界値問題

第5章 偏微分方程式の境界値問題 October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ

More information

K E N Z U 01 7 16 HP M. 1 1 4 1.1 3.......................... 4 1.................................... 4 1..1..................................... 4 1...................................... 5................................

More information

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising ,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising Model 1 Ising 1 Ising Model N Ising (σ i = ±1) (Free

More information

10_4.dvi

10_4.dvi Vol.44, No.1, 1/7 28 Synchronzed Control for Blateral Teleoperaton wth Dfferent Confguratons and Communcaton Delays Hsanosuke Kawada,KoueYoshda and Toru Namerkawa Ths paper addresses the problem of the

More information

A 99% MS-Free Presentation

A 99% MS-Free Presentation A 99% MS-Free Presentation 2 Galactic Dynamics (Binney & Tremaine 1987, 2008) Dynamics of Galaxies (Bertin 2000) Dynamical Evolution of Globular Clusters (Spitzer 1987) The Gravitational Million-Body Problem

More information

MS#sugaku(ver.2).dvi

MS#sugaku(ver.2).dvi 1 1,,,.,,,,,.,.,,,.,, Computer-Aided Design).,,, Boltzmann,, [1]., Anderson, []., Anderson, Schrödinger [3],[4]., nm,,.,,,.,, Schrödinger.,, [5],[6].,,.,,, 1 .,, -. -, -., -,,,. Wigner-Boltzmann Schrödinger,

More information