TokyoUnivITC

Size: px

Start display at page:

Download "TokyoUnivITC"

こうだいのあき
5 years ago
Views:

1 1

2 2

3 3

4 + r T r ( T )=S T solve(fvm::ddt(t) + fvm::div(phi,t) - fvm::laplacian(dt,t) == fvoptions(t)); 4

5 5

6 標準ソルバの解析例 MRFSimpleFoam 回転攪拌槽の流れ interdymfoam 撹拌槽内の流れ pisofoam LESによるバックステップ流れ図出典 OpenFOAMチュートリアルドキュメント作成プロジェクト東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 6

7 移動格子ソルバの解析例ピストン押し込み流れスクリューの回転流れ場翼型の6自由度剛体運動図出典 OpenFOAMチュートリアルドキュメント作成プロジェクト東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 7

8 8

9 9

10 10

11 11

12 OpenFOAM (2014/2) の新機能図出典 The OpenFOAM Foundation 境界適合六面体メッシャーsnappyHexMeshの改善新自動六面体メッシャーfoamyHexMesh 並列計算時の線型ソルバGAMGの高速化 VOF法多相流ソルバにおけるMulti-dimensional limiter for explicit solution(mules)の半陰解版 Arbitary Mesh Interface(AMI)の改善拡張密粒子流れ用のDiscrete Particle Modelling MultiPhase Particle-in-Cell法移動格子での球面線型補間を用いたモーフィング東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 12

13 13

14 14

[cfmesh,salome,gmsh 商用メッシャー等] 解析ソルバ [icofoam等] 様々な結果後処理必要あれば格子変換領域統合(並列計算時)

15 OpenFOAMの代表的な解析手順前処理(格子生成など) 解析格子生成初期設定 [blockmesh [setfields等] snappyhexmeshなど] 後処理(可視化など) 可視化 [ParaView Visit 商用可視化ツール等] 領域分割(並列計算時) [decomposepar] または格子生成 [cfmesh,salome,gmsh 商用メッシャー等] 解析ソルバ [icofoam等] 様々な結果後処理必要あれば格子変換領域統合(並列計算時) [gmshtofoam等] [reconstructpar] [sample等] 東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 15

16 16

17 : ( ) parafoam OpenFOAM 17

18 U =1m/s d =0.1m =0.01m 2 /s 1m/s Re = du = m icofoam : ( ) parafoam OpenFOAM 0.1m 18

19 Oakleaf- FX lecture.tar.gz tar zxf lecture.tar.gz cd ~/lecture/cavity tree. 0 U p blockmesh.pjm constant polymesh blockmeshdict boundary transportproperties icofoam.pjm system controldict fvschemes fvsolution 19

20 (*.PJM) 20

21 blockmesh icofoam 21

22 22

23 (OpenFOAM!! 3 z ) キャビティケースの格子分割 20分割 (Δy=5mm) 1分割 (Δz=10mm) 20分割(Δx=5mm) 注)2次元なので z方向は1分割(幅は任意) 図出典: 大嶋拓也 (新潟大学) キャビティ流れの解析 parafoamの実習第一回OpenFOAM講習会東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 23

24 converttometers 0.1; vertices ( ); (0 0 0) 0 (1 0 0) (1 1 0) (0 1 0) ( ) ( ) ( ) ( ) 7 24

25 blocks ( hex ( ) ( ) simplegrading (1 1 1) ( 1: ) ); edges ( ); 25

26 boundary ( movingwall { type wall; faces ( ( ) ); } fixedwalls ( { } type wall; faces ( ( ) ( ) ( ) ); ); frontandback { type empty; faces ( ( ) ( ) ); } 26

27 more constant/polymesh/blockmeshdict more ( more) SPC b h q 27

28 28

29 blockmesh.pjm #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" # ( tutorial) #PJM - g gt00 #( ) #PJM - L "node=1" # #PJM - j # # Program execution # module load OpenFOAM/2.3.0 #OpenFOAM/2.3.0 module source $WM_PROJECT_DIR/etc/bashrc #OpenFOAM blockmesh > log.blockmesh #blockmesh 29

30 pjsub blockmesh.pjm [INFO] PJM 0000 pjsub Job JOB_ID submitted. JOB_ID pjstat JOB_ID JOB_NAME STATUS PROJECT RSCGROUP JOB_ID blockmesh. QUEUED gt00 lecture JOB_ID JOB_NAME STATUS PROJECT RSCGROUP JOB_ID blockmesh. RUNNING gt00 lecture No unfinished job found. 30

31 tree 0 U p blockmesh.pjm blockmesh.pjm.o constant polymesh blockmeshdict boundary faces neighbour owner points transportproperties icofoam.pjm log.blockmesh system controldict fvschemes fvsolution (JOB_ID ) blockmesh more blockmesh.pjm.o* *() JOB_ID 31

32 more log.blockmesh Mesh Information boundingbox: (0 0 0) ( ) npoints: 882 ncells: 400 nfaces: 1640 ( ) ninternalfaces: 760 ( ) Patches ( ) patch 0 (start: 760 size: 20) name: movingwall patch 1 (start: 780 size: 60) name: fixedwalls patch 2 (start: 840 size: 800) name: frontandback End 32

33 33

34 34

35 ~/lecture/ oakleaf- fx.cc.u- tokyo.ac.jp ~ $ cd ~ cd $ scp - r tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/lecture./ tyyxxx= scp rsync passphrase $ cd lecture/cavity 35

36 36

37 ~/lecture/cavity $ touch pv.foam touch 37

38 38

39 39

40 more 0/U 40

41 dimensions [ ]; internalfield uniform (0 0 0); 41

42 boundaryfield { movingwall { type fixedvalue; value uniform (1 0 0); } fixedwalls { type fixedvalue; value uniform (0 0 0); } frontandback { type empty; } } 42

43 more 0/p dimensions [ ]; internalfield uniform 0; 43

44 boundaryfield { movingwall { } type zerogradient; 0 fixedwalls { } type zerogradient; } frontandback { type empty; 2 } 44

45 45

46 system/controldict application icofoam; startfrom starttime; ( latesttime ) starttime 0; stopat endtime; ( nextwrite ) endtime 0.5; deltat 0.005; ( [s]) writecontrol timestep; writeinterval 20; (20time step=0.1s ) writeformat ascii; (ascii, binary) binary FX10 Big endian binary intel Little endianasciibinary writeformat ascii foamformatconvert ascii Little endian writeprecision 6; ( ascii ) writecompression off; (off, on) timeformat general; timeprecision 6; runtimemodifiable true; 46

47 47

48 icofoam.pjm ( ) #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" # ( tutorial) #PJM - g gt00 #( ) #PJM - L "node=1" # #PJM - j # # Program execution # module load OpenFOAM/2.3.0 #OpenFOAM/2.3.0 module source $WM_PROJECT_DIR/etc/bashrc #OpenFOAM icofoam> log # icofoam log pjsub icofoam.pjm 48

49 tree U p phi uniform: time icofoam.pjm.o (JOB_ID ) log 49

50 more log Build : f5222ca19ce6 Exec Date Time Host PID Case : icofoam : Jan : 00:00:00 : "a01t70081" : : /home/******/lecture/cavity nprocs : 1 sigfpe : Enabling floating point exception trapping (FOAM_SIGFPE). filemodificationchecking : Monitoring run- time modified files using timestampmaster allowsystemoperations : Disallowing user- supplied system call operations 50

51 Starting time loop ( ) Time = ( ) Courant Number mean: 0 max: 0 smoothsolver: Solving for Ux, Initial residual = 1, Final residual = e- 06, No Iterations 19 Ux( x ) smoothsolver: Initial residual: Final residual: No Iterations 1 smoothsolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 Uy ( y ) ICPCG: Solving for p, Initial residual = 1, Final residual = e- 07, No Iterations 35 p( ) DICPCG: 51

52 time step continuity errors : sum local = e- 09, global = e- 21, cumulative = e- 21 sum local : global : ( ) cumulative : global ExecutionTime = 0.05 s ClockTime = 0 s ExecutionTime:, ClockTime : I/O Time = 0.01 Courant Number mean: max: Time = 2 ( ) 52

53 53

54 54

55 ~/lecture/cavity oakleaf- fx.cc.u- tokyo.ac.jp ~/lecture/cavity/ $ rsync - auv tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/lecture/cavity/./ /( ) a=archive() u=update(), v=verbose() receiving file list... done./ log sent 592 bytes received bytes bytes/sec total size is speedup is 1.23 icofoam (rsync ) 55

56 56

57 57

58 58

59 59

60 60

61 61

dambreakfineケース Koshizukaら[1]によるダム崩壊(dam break)の実験をvof(volume of Fluid) 法の二相流ソルバinterFoamを用いて二次元層流モデルで解析したもの解析図出典 SM ただで始める流体解析第11回オープンCAE勉強会@岐阜実験[1] OpenFOAM 実験[1] OpenFOAM [1] Koshizuka, S., H.

62 dambreakfineケース Koshizukaら[1]によるダム崩壊(dam break)の実験をvof(volume of Fluid) 法の二相流ソルバinterFoamを用いて二次元層流モデルで解析したもの解析図出典 SM ただで始める流体解析実験[1] OpenFOAM 実験[1] OpenFOAM [1] Koshizuka, S., H. Tamako and Y. Oka : A Particle Method for Incompressible Viscous Flow with Fluid Fragmentation, CFD Journal,Vol.4, No.1, pp.29-46, 1995 東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 62

U p g F r U =0 @ U @t @ @t apple + r ( UU) r = rp + g + F + r ( U)+r [(1 ) U U r (= U l U g ) (= appler ) r]=0 [2] Weller, H.G., 2008.

63 U p g F r U apple + r ( UU) r = rp + g + F + r ( U)+r [(1 ) U U r (= U l U g ) (= appler ) r]=0 [2] Weller, H.G., A new approach to VOF-based interface capturing methods for incompressible and compressible flows. Technical Report No. TR/HGW/04 [3] A. Albadawi, D.B. et al., Influence of surface tension implementation in Volume of Fluid and coupled Volume of Fluid with Level Set methods for bubble growth and detachment, International Journal of Multiphase Flow, Volume 53, pp l g 63

64 0.584 m 水柱 m m m m m m 64

65 cd ~/lecture/dambreakfine/ tree. 0 U alpha.water.org p_rgh Allrun- p.pjm Allrun.PJM constant g polymesh blockmeshdict boundary transportproperties turbulenceproperties system controldict decomposepardict fvschemes fvsolution setfieldsdict 65

66 66

67 alpha.water setfields system/setfieldsdict 1 0 defaultfieldvalues// ( volscalarfieldvalue alpha.water 0 ); regions // ( boxtocell // { } ); box (0 0-1) ( ); // fieldvalues // ( volscalarfieldvalue alpha.water 1 ); 67

68 atmosphere ( ) m rightwall leftwall 水柱 m 0/U boundaryfield // { : lowerwall // (leftwall,rightwall ) { type fixedvalue; value uniform (0 0 0); 0.584} m atmosphere // { type pressureinletoutletvelocity; value uniform (0 0 0); } } m m m lowerwall 68

69 ( ) m rightwall leftwall 水柱 m m atmosphere m m lowerwall 0/p_rgh boundaryfield // { : lowerwall // (leftwall,rightwall ) { type fixedfluxpressure; value uniform 0; } 0.584atmosphere m // { type totalpressure; p0 uniform 0; U U; phi phi; rho rho; psi none; gamma 1; value uniform 0; } } 69

70 水柱 atmosphere ( ) m rightwall leftwall 0/alpha.water boundaryfield // { : lowerwall // (leftwall,rightwall ) { 0.584} m type zerogradient; m m m m lowerwall } atmosphere // { type inletoutlet; inletvalue uniform 0; value uniform 0; } 70

71 constant/g dimensions [ ]; // [m/s 2 ] value ( ); // (y ) constant/turbulenceproperties simulationtype laminar; // ( ) 71

72 constant/transportproperties phases (water air); // water // ( ) { transportmodel Newtonian; // nu nu [ ] 1e- 06; // [m 2 /s] rho rho [ ] 1000; // [kg/m 3 ] //transportmodel Newtonian2 CrossPowerLawCoeffs { } BirdCarreauCoeffs { } } air // ( ) { transportmodel Newtonian; nu nu [ ] 1.48e- 05; rho rho [ ] 1; } sigma sigma [ ] 0.07; // [N/m](=[kg/s 2 ]) 72

73 system/controldict Application interfoam; // starttime 0; // [s] endtime 1; // [s] deltat 0.001; // [s] ( ) writeinterval 0.05; // [s] adjusttimestep yes; // // // maxco 1; // maxalphaco 1; // (0< <1) maxdeltat 1; // 73

74 blockmesh cp 0/ alpha.water.org 0/alpha.water setfields interfoam 74

75 4. [ParaView] 2. [pjsub, blockmesh, setfields, interfoam] 75

76 Allrun.PJM ( ) #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" # ( tutorial) #PJM - g gt00 #( ) #PJM - L "node=1" # #PJM - j # # Program execution # module load OpenFOAM/2.3.0 #OpenFOAM/2.3.0 module source $WM_PROJECT_DIR/etc/bashrc #OpenFOAM blockmesh > log.blockmesh # cp 0/alpha.water.org 0/alpha.water # setfields > log.setfields # interfoam > log # pjsub Allrun.PJM 76

77 more log.blockmesh Mesh Information boundingbox: (0 0 0) ( ) npoints: ncells: 7700 nfaces: ninternalfaces: Patches ( ) ( ) ( ) patch 0 (start: size: 86) name: leftwall patch 1 (start: size: 86) name: rightwall patch 2 (start: size: 110) name: lowerwall patch 3 (start: size: 90) name: atmosphere patch 4 (start: size: 15400) name: defaultfaces End 77

78 more log Create time Create mesh for time = 0 PIMPLE: Operating solver in PISO mode Reading field p_rgh Reading field U Reading/calculating face flux field phi Reading transportproperties Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian (0%) more SPC b h q 78

79 Courant Number mean: max: Interface Courant Number mean: max: deltat = ( 1 ) Time = 0.15 (0.05[s] ) smoothsolver: Solving for alpha.water, Initial residual = , Final residual = e - 09, No Iterations 3 alpha.water : Initial residual: Final residual: No Iterations: Phase- 1 volume fraction = Min(alpha1) = e- 09 Max(alpha1) = 1 : DICPCG: Solving for p_rgh, Initial residual = e- 05, Final residual = e- 08, No Itrations 44 p_rgh Poisson time step continuity errors : sum local = e- 08, global = e- 09, cumulative = sum local : global : ( ) cumulative : global ExecutionTime = 0.77 s ClockTime = 1 s ExecutionTime:, ClockTime : I/O 79

80 ~/lecture/dambreakfine/ oakleaf- fx.cc.u- tokyo.ac.jp ~/lecture/cavity $ ~/lecture/dambreakfine $ rsync - auv tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/lecture/dambreakfine/./ $ touch pv.foam ParaViewpv.foam 80

81 81

82 82

83 演習課題課題1: 水相率を補間して表示させる課題2: コマ送りして 0.1秒毎の実験可視化画像[1]と比較する課題3: 並列計算の結果と比較するため全時間ステップの可視化画像を保存する方法 Fileメニュー/Save Animation/Save Animationボタン/File name: serial Files of type: PNG(静止画のほうが比較が容易であるため)/OK serial.連番.pngのpngファイルが作成されるので画像ビューワ等で表示する [1] Koshizuka, S., H. Tamako and Y. Oka : A Particle Method for Incompressible Viscous Flow with Fluid Fragmentation, CFD Journal,Vol.4, No.1, pp.29-46, 1995 東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 83

84 84

OpenFOAMの並列計算手法 1. 格子生成 2. 領域分割 (decomposepar) [1]櫻井片桐ら OpenFOAMへの疎行列計算ライブラリ Xabclibの適用と評価オープンCAEシンポジウム2014 3.

85 OpenFOAMの並列計算手法 1. 格子生成 2. 領域分割 (decomposepar) [1]櫻井片桐ら OpenFOAMへの疎行列計算ライブラリ Xabclibの適用と評価オープンCAEシンポジウム MPI並列でソルバを実行 (MPI+OpenMPのハイブリット並列は標準では未実装ただし櫻井片桐らによる線型ソルバのスレッド並列化(1MPI+OpenMP)の研究例有り[1] 現在 MPI+OpenMP実装を開発中) 4. 領域毎の解析結果を再構築 (reconstructpar) プロセッサ間MPI通信 (フラット MPI通信) プロセッサ間通信は通常ボトルネックとなる processor 2 processor 3 processor 0 processor 1 dambreakfineチュートリアル東京大学情報基盤センター第45回 FX10 スーパーコンピュータシステムお試しアカウント付き並列プログラミング講習会 OpenFOAM初級入門 85

86 4. [ParaView] 2. [pjsub,blockmesh, setfields, decomposepar, interfoam, recconstructpar] 86

87 system/decomposepardict numberofsubdomains 4; // method simple; // simplecoeffs// { } n delta 0.001; ( ); // hierarchicalcoeffs // { } n ( ); order xyz; // delta 0.001; simple hierarchical metis scotch metis scotch manual 87

3 4 2 4 metis scotch 1 2 3 4 1 3 1 4 scotchcoeffs

88 metis scotch scotchcoeffs metiscoeffs { processorweights ( ); // }

89 Allrun- p.pjm ( MPI ) #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" # ( tutorial) #PJM - g gt00 #( ) #PJM - L "node=1" # ( 12 ) #PJM - - mpi "proc=4" #(MPI 192 ) #PJM - j # # Program execution # module load OpenFOAM/2.3.0 source $WM_PROJECT_DIR/etc/bashrc blockmesh > log.blockmesh # cp 0/alpha.water.org 0/alpha.water # setfields > log.setfields # decomposepar > log.decomposepar # mpiexec - stdout log- np4 - np 4 interfoam - parallel #MPI reconstructpar > log.reconstructpar # 89

90 pjsub Allrun- p.pjm more log.decomposepar processor 2 processor 3 ( ) Processor 0 0 Number of cells = ( ) Number of faces shared with processor 1 = 44 Number of faces shared with processor 2 = 46 Number of processor patches = 2 Number of processor faces = 90 (=44+46) processor 0 Number of boundary faces = 3958 processor 1 Processor 1 1( processor1 ) Number of cells = 1915 ( 0 ) ( ) 90

91 ( ) Number of processor faces = 178 ( ) Max number of cells = 1935 ( % above average 1925) Max number of processor patches = 2 (0% above average 2) Max number of faces between processors = 90 (1.1236% above average 89) Time = 0 0 Processor 0: field transfer Processor 1: field transfer Processor 2: field transfer Processor 3: field transfer End. processor 2 processor 0 processor 3 processor 1 91

92 processor 2 processor 3 constant/ polymesh/ 0/ ( ) alpha.water,u,p_rgh processor0/ constant polymesh/ processor 0 processor 1 0/ alpha.water,u,p_rgh processor1/ processor0 processor2/ processor3/ 92

93 more log- np4 nprocs : 4 Slaves : 3 ( - 1) (.PID( ID) "a01t " PID "a01t " "a01t " ) Pstream initialized with: () floattransfer : 0 nprocssimplesum : 0 commstype : nonblocking polling iterations : 0 93

94 processor0/ 0.05/ alpha*,u,p,p_rgh,phi : 0.4/ processor1/ 0.05/ alpha*,u,p,p_rgh,phi processor2/ 0.05/ alpha*,u,p,p_rgh,phi processor3/ 0.05/ alpha*,u,p,p_rgh,phi () 0.05/ alpha*,u,p,p_rgh,phi 0.1/ alpha*,u,p,p_rgh,phi : 0.4/ alpha*,u,p,p_rgh,phi 94

95 ~/lecture/dambreakfine/ oakleaf- fx.cc.u- tokyo.ac.jp ~/lecture/dambreakfine $ rsync - auv tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/lecture/dambreakfine/./ 95

96 96

97 1: File /Save Animation/Save Animation /File name: parallel Files of type: PNG( )/OK parallel..png PNG 2: ParaView (Decomposed Case) Case Type/Decomposed Case Apply 3: Decomposed Case1 (Reconstructed Case) 4: (Strong scaling) (log) (log- np4) ExecutionTime t1,t4 =t1/t4 =(t1/t4)/4 97

98 98

99 [1] K Iwamoto, K., 2002."Database of Fully Developed Channel Flow,"THTLAB Internal Report, No. ILR-0201.THTLAB, Dept. of Mech. Eng., The Univ. of Tokyo. Lx Ly Lz = 5π 2 2π Reτ = uτ δ/ µ = 110 [-] uτ [m/s] δ [m] ( =Ly / 2 ) µ [m 2 /s 2 ] x x z 99

100 SP SP = TS / TP TS TP E P EP = SP / P 100 [%] T8 < T12 100

101 101

102 cd ~/lecture/channelretau110 tree Oakleaf-FX batchscript blockmesh.sh blockmesh decomposepar.sh decomposepar solve.sh (pimplefoam) benchmark.conf (benchmark.sh ) README bin benchmark.sh fippplot.py fipp (Python) fipptable.sh fipp plot.py (Python) table.sh template channelretau

103 cd Oakleaf- FX emacs - nw benchmark.conf vi,gedit narray( n ) n (nx=120,ny=65,nz=48) narray n # narray=(\ # "00001" \ # "00002" \ # "00004" \ # "00008" \ ) 103

104 mpiarray(mpi mpi ) MPI mpi mpiarray mpi 0 5 mpiarray=( "00192" "00128" "00064" "00032" "00016" ) simulationtypesarray( ) simulationtypesarray (simulationtypes )- LESModel_(LESModel )- delta_(delta )" LES Smagorinsky # simulationtypesarray=(\ "laminar- LESModel_laminar- delta_cuberootvol" \ # "LESModel- LESModel_Smagorinsky- delta_vandriest" \ ) 104

105 solversarray() solversarray PCG- preconditioner_(preconditioner )" (PCG ) GAMG- smoother_(smoother ) (GAMG ) MPI GAMGPCG # solversarray=(\ # "PCG- preconditioner_fdic" \ # "PCG- preconditioner_dic" \ # "PCG- preconditioner_diagonal" \ # "GAMG- smoother_fdic" \ # "GAMG- smoother_dic" \ # "GAMG- smoother_dicgaussseidel" \ # "GAMG- smoother_gaussseidel" \ # "GAMG- smoother_nonblockinggaussseidel" \ # "GAMG- smoother_symgaussseidel" \ ) 105

106 ../bin/benchmark.sh >& log.benchmark.sh & tail - f log.benchmark.sh ^C() pjstat 106

107 107

108 ~/lecture/channelretau110/oakleaffx10/batchscript/solve.sh!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" #PJM - g "gt00" #PJM -j #PJM - S # Program execution # module load OpenFOAM/2.3.0 source $WM_PROJECT_DIR/etc/bashrc getnumberofprocessors getapplication bash. $WM_PROJECT_DIR/bin/tools/RunFunctions env log=log.${pjm_jobid} fipp fipp - d $log.fippdir - C - Ihwm,call - H - L shared \ mpiexec - - stdout $log - np $(getnumberofprocessors) $(getapplication) - parallel 108

109 ~/lecture/channelretau110/bin/fipptable.sh fipppx - A - Ihwm,call - d > ( 0.74M, 16MPI, laminar, PCG- FDIC ) Performance monitor event:statistics ( ) ************************************************************************************* Application - performance monitors ************************************************************************************* Elapsed(s) MFLOPS MFLOPS/PEAK(%) MIPS MIPS/PEAK(%) Application Mem throughput Mem throughput Elapsed(s) _chip(mb/s) /PEAK(%) SIMD(%) Application ( ) # ##### 26% <4> Foam::FDICPreconditioner::precondition(Foam::Field<double> &, const Foam::Field<double> &, unsigned char) const [57] 109

110 ../bin/table.sh table.csv../bin/plot.py plot.pdf../bin/fipptable.sh fipptable.csv../bin/fippplot.py ~/lecture/dambreakfine $ cd../channelretau110/oakleaf- ~/lecture/channelretau110/oakleaf- FX $ rsync - auv tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/lecture/channelretau110/ Oakleaf- FX/./ - - exclude=n_* n_ (*.pdf)(*.csv) 110

111 = (T11step-T1step)/10 111

112 112

113 113

114 DNS basic combustion compressible discretemethods electromagnetics financial heattransfer incompressible mesh multiphase lagrangian resources stressanalysis 114

115 cd cp - r /usr/local/openfoam/2.3.0/openfoam /tutorials/. cd tutorials ls Allclean Allrun Alltest DNS basic combustion compressible discretemethods electromagnetics financial heattransfer incompressible lagrangian mesh multiphase resources stressanalysis cd incompressible/icofoam/cavity ls 0 constant system cp ~/lecture/foamruntutorials.pjm./ 115

116 foamruntutorials.pjm #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" #PJM - g gt00 #PJM - L "node=1" #PJM - - mpi "proc=1" 1MPI( ) #PJM - j # Program execution # module load OpenFOAM/2.3.0 source $WM_PROJECT_DIR/etc/bashrc foamruntutorials (Allrun blockmesh ) pjsub foamruntutorials.pjm pjstat tail - f log.icofoam 116

117 cd ~/tutorials cd incompressible/simplefoam/motorbike ls 0.org Allclean Allrun constant system more Allrun runparallel snappyhexmesh 6 -overwrite runparallel patchsummary 6 runparallel potentialfoam 6 -nofunctionobjects -writep runparallel $(getapplication) 6 runparallel $application 6 runparallel( bash ) 6 117

118 cp ~/lecture/foamruntutorials.pjm./ emacs - nw foamruntutorials.pjm #!/bin/bash # pjsub option # #PJM - L "rscgrp=lecture" #PJM - g gt00 #PJM - L "node=1" #PJM - - mpi proc= 1 6" 6MPI #PJM - j # Program execution # module load OpenFOAM/2.3.0 source $WM_PROJECT_DIR/etc/bashrc foamruntutorials 118

119 emacs - nw system/controldict endtime 5 100; lecure 15 pjsub foamruntutorials.pjm pjstat (RUNNING) tail - f foamruntutorials.pjm.o* foamruntutorials.pjm.ojobid (foamruntutorials mpiexec stdout ) 119

120 @ ~/ $ cd $ mkdir tutorials $ cd tutorials $ rsync - auv tyyxxx@oakleaf- fx.cc.u- tokyo.ac.jp:~/tutorials/./ - - exclude=processor* (processor*) $ cd incompressible/simplefoam/motorbike $ touch pv.foam ParaView 120

121 foamruntutorials ~/.bashrc module load OpenFOAM/2.3.0 #OpenFOAM/2.3.0 module source $WM_PROJECT_DIR/etc/bashrc #OpenFOAM foamcleantutorials 121

122 tokyo.ac.jp/system/fx10/ office.jp/pages/ws_openfoam_ office.jp/pages/ws_openfoam2_

Microsoft PowerPoint OpenFOAMの使い方(柴田).ppt [互換モード]

Microsoft PowerPoint OpenFOAMの使い方(柴田).ppt [互換モード] OpenFOAM 勉強会 for beginner 2 期第 1 回 OpenFOAM の使い方柴田貴裕目標既に OpenFOAM はインストール済みさまざまなチュートリアルに取り組めるようにチュートリアルの実行の方法を scalartransportfoam を例に用いて見ていくインストールは OpenCFD の HP の方法に従えば比較的容易にできる ) OpenCFD 社の HP http://www.openfoam.com/