OpenFOAM による熱移動と流れの数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/069101 このサンプルページの内容は, 初版 1 刷発行時のものです.
OpenFOAM R OpenCFD OpenFOAM OpenFOAM R OpenCFD R OpenCFD R OPENFOAM R is a registered trade mark of OpenCFD Limited, producer and distributor of the OpenFOAM software. This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software and owner of the OPENFOAM R and OpenCFD R trade marks. The authers and the editor are not associated to OpenCFD.
i OpenFOAM 1. OpenFOAM 2. OpenFOAM 3. OpenFOAM 4. OpenFOAM 34 OpenFOAM 1 OpenFOAM OpenFOAM Linux OpenFOAM Linux B Linux OpenFOAM OpenFOAM [19] Versteeg [17] 2 4 5
ii 2 OpenFOAM OpenFOAM 3 D 4 OpenFOAM 5 OpenFOAM A OpenFOAM DEXCS2015 for OpenFOAM B Linux Linux C OpenFOAM () ParaView D FreeCAD OpenFOAM ParaViewFreeCAD OpenFOAM 2.4.0 3.0.x (x 0 1) Linux B Linux Web Linux 2 OpenFOAM C ParaView ParaView 3 D 4 () () OpenFOAM Linux ParaView Linux ParaView
iii CAE IHI OCAEL 2,3 CAE OCAEL OCSEˆ2
iv 1 1 1.1 OpenFOAM 1 1.2 OpenFOAM 2 1.3 OpenFOAM 2 1.4 3 1.5 OpenFOAM 4 1.6 OpenFOAM 4 1.7 5 2 OpenFOAM 8 2.1 OpenFOAM 8 2.2 12 2.3 14 2.4 21 2.5 24 3 33 3.1 OpenFOAM 33 3.2 blockmesh 36 3.3 snappyhexmesh 42 3.4 55 3.5 56 3.6 57 3.7 57 4 OpenFOAM 59 4.1 59 4.2 59 4.3 65 4.4 66 4.5 70 4.6 80 4.7 88
v 4.8 93 4.9 97 4.10 101 4.11 102 4.12 104 4.13 108 4.14 109 5 OpenFOAM 128 5.1 128 5.2 132 5.3 133 5.4 135 5.5 138 5.6 140 5.7 150 A DEXCS2015 for OpenFOAM 166 A.1 166 A.2 166 A.3 VirtualBox 166 A.4 166 A.5 DEXCS2015 for OpenFOAM 170 A.6 DEXCS 176 A.7 OpenFOAM 177 A.8 177 B Linux 179 B.1 179 B.2 Linux 179 B.3 181 B.4 182 B.5 182 B.6 184 B.7 185 B.8 185
vi B.9 185 B.10 187 B.11 187 B.12 Linux 196 C ParaView 197 C.1 197 C.2 197 C.3 ParaView 197 C.4 199 C.5 200 C.6 200 C.7 200 C.8 201 C.9 202 C.10 202 C.11 203 C.12 203 C.13 204 C.14 204 C.15 204 C.16 205 C.17 206 C.18 206 C.19 207 C.20 207 C.21 207 C.22 207 D FreeCAD 208 D.1 208 D.2 FreeCAD 208 D.3 213 D.4 216 218 220
1 1 1.1 OpenFOAM OpenFOAM (Open source Field Operation And Manipulation ) GNU General Public License (GPL) (Computational Fluid Dynamics: CFD) C++ C++ solve ( fvm::ddt(t) + fvm::div(phi, T) - fvm::laplacian(dt, T) ); (field) OpenFOAM OpenFOAM OpenFOAM
2 1 Imperial College Nabla FOAM 2004 OpenCFD OpenFOAM 2011 8 OpenCFD Silicon Graphics (SGI) OpenFOAM OpenFOAM Foundation 2012 9 OpenCFD ESI Group OpenFOAM OpenCFD 1.2 OpenFOAM / / (CHT) (/) 1.3 OpenFOAM OpenFOAM Linux C++ OpenFOAM OS Linux Linux C++
59 4 OpenFOAM OpenFOAM 2 4.1 $FOAM TUTORIALS tut $ tut 4.2 4.2.1 5.1.1 (p.128)5.1.3 (p.128) simplefoam constant transportproperties constant/transportproperties FoamFile version 2.0; format ascii; class dictionary; location "constant";
60 4 OpenFOAM object transportproperties; transportmodel Newtonian; nu nu [ 0 2-1 0 0 0 0 ] 1e-05; CrossPowerLawCoeffs nu0 nu0 [ 0 2-1 0 0 0 0 ] 1e-06; nuinf nuinf [ 0 2-1 0 0 0 0 ] 1e-06; m m [ 0 0 1 0 0 0 0 ] 1; n n [ 0 0 0 0 0 0 0 ] 1; BirdCarreauCoeffs nu0 nu0 [ 0 2-1 0 0 0 0 ] 1e-06; nuinf nuinf [ 0 2-1 0 0 0 0 ] 1e-06; k k [ 0 0 1 0 0 0 0 ] 0; n n [ 0 0 0 0 0 0 0 ] 1; transportmodel Newtonian powerlaw BirdCarreau CrossPowerLaw HerschelBulkley (Newtonian) (Newtonian ) Newtonian Newtonian nu () nu nu [ 0 2-1 0 0 0 0 ] 1e-05; // [m 2 /s] Boussinesq 5.1.4 (p.129) buoyantboussinesqsimplefoam Boussinesq constant/transportproperties
4.2 61 constant/transportproperties // Thermal expansion coefficient beta beta [0 0 0-1 0 0 0] 3e-03; // [1/K] // Reference temperature TRef TRef [0 0 0 1 0 0 0] 300; // [K] // Laminar Prandtl number Pr Pr [0 0 0 0 0 0 0] 0.9; // // Turbulent Prandtl number Prt Prt [0 0 0 0 0 0 0] 0.7; // ν Pr (/) α α = ν (4.1) Pr 4.2.2 5.1.1 (p.128)5.1.3 (p.128)5.1.4 (p.129)5.1.5 (p.130) buoyantsimplefoam constant/ thermophysicalproperties constant/thermophysicalproperties FoamFile version 2.0; format ascii; class dictionary; location "constant"; object thermophysicalproperties; thermotype type herhothermo; mixture puremixture; specie specie; equationofstate perfectgas; transport const; thermo hconst; energy sensibleenthalpy; mixture specie
62 4 OpenFOAM nmoles 1; molweight 28.9; thermodynamics Cp 1000; Hf 0; transport mu 1.8e-05; Pr 0.7; thermotype type (thermophysical model) hepsithermo herhothermo Rho Psi (compressibility) mixture puremixturemulticomponentmixture reactingmixture puremixture mixture thermotype specie specie specie nmoles 1; // molweight 28.9; //
equationofstate 4.2 63 rhoconst equationofstate rho 1000; // [kg/m 3 ] perfectgas () ρ p T ρ = p RT R (4.2) incompressibleperfectgas () equationofstate pref 101325; // [Pa] p ref ρ = p ref RT (4.3) icopolynomial equationofstate rhocoeffs<8> ( 4.0097-0.016954 3.3057e-05-3.0042e-08 1.0286e-11 0 0 0 ); transport
64 4 OpenFOAM const (/) transport mu 1.8e-05; // [Pa-s] Pr 0.7; // polynomial transport mucoeffs<8> ( 1.5061e-06 6.16e-08-1.819e-11 0 0 0 0 0 ); kappacoeffs<8> ( 0.0025219 8.506e-05-1.312e-08 0 0 0 0 0 ); sutherland Sutherland transport As 1.67212e-06; Ts 170.672; Sutherland μ μ = A s T 1+T s /T (4.4) A s, T s k Eucken (modified Eucken correlation) ( k = μc v 1.32 + 1.77R ) (4.5) C v C v [J/kg-K] thermo
4.3 65 hconst thermodynamics Cp 1000; // [J/kg-K] Hf 0; // [J/kg] hpolynomial thermodynamics Hf 0; // [J/kg] Sf 0; // [J/kg-K] CpCoeffs<8> ( 948.76 0.39171-0.00095999 1.393e-06-6.2029e-10 0 0 0 ); CpCoeffs energy sensibleenthalpy sensibleinternalenergy absoluteenthalpy absoluteinternalenergy 4.3 5.1.6 (p.131) buoyantsimplefoam constant/g constant/g FoamFile version 2.0; format ascii; class uniformdimensionedvectorfield; location "constant"; object g;
128 5 OpenFOAM OpenFOAM 5.1 5.1.1 4.2.1 (p.59)4.2.2 (p.61) (incompressible fluid) (compressible fluid) 0.3 ( 5 ) OpenFOAM (thermophysical properties) 5.1.2 () ρ + (ρu) = 0 (5.1) t u ρ u =0 (5.2) 5.1.3 4.2.1 (p.59)4.2.2 (p.61) ( ) ρu t + (ρuu) = p + [μ u +( u) T ] ( ) 2 3 μ u (5.3) μ
5.1 129 OpenFOAM [ ρu + (ρuu) = p + (μ u)+ μ ( u) T 23 ] t ui (5.4) 2 3 ( 2/3 1/3) [ u + (uu) = p + (ν u)+ ν ( u) T 13 ] t ui (5.5) p ( )ν 3 5.1.4 4.2.1 (p.60)4.2.2 (p.61) E (ρe)+ (ρeu) = (pu)+ (k T) (5.6) t k T E e K E = e + K (5.7) K = 1 2 u u (5.8) (5.6) t (ρe)+ (ρeu)+ (ρk)+ (ρku) = (pu)+ (k T) (5.9) t h h = e + p ρ (5.10) (5.9) [OpenFOAM-3.0.x] ui 2/3
130 5 OpenFOAM t (ρh)+ (ρhu)+ (ρk)+ (ρku) = p + (k T ) (5.11) t t h c p h = h a h 0 = T T 0 c p dt (5.12) h a h 0 T 0 T 0 = 293.15 [K] h h = c p (T T 0 ) (5.13) (5.11) T t α = k/(ρc p ) + (T u) = (α T ) (5.14) OpenFOAM (5.9) (5.11) (k T )= (α e) = (α h) (5.15) α = k/c p Boussinesq (5.14) 5.1.5 4.2.2 (p.61) (equation of state) (ideal gas perfect gas ) ρ = pw RT (5.16) W [kg/mol]r [J/mol-K]
220 A adjustableruntime 99, 100 adjusttimestep 126 Apply 29, 197, 200, 202 206 ascii 100 B backward 82 binary 100 BirdCarreau 60 blockmesh 10, 12, 22, 25, 34, 36, 41 43, 46, 53, 55 blockmeshdict 34, 41, 45, 46 boundary 40, 73 bounded 84 146 buoyantboussinesqpimplefoam 9 buoyantboussinesqsimplefoam 9, 60, 71, 80 buoyantpimplefoam 9, 25, 30, 117 buoyantsimplefoam 9, 61, 65, 71, 78, 79 C calculated 120 castellatedmeshcontrols 47, 50, 53 celllimited 83 cellmdlimited 83 checkmesh 10, 52, 56, 85, 94 Choose Preset 201 class 72 Coloring 26, 28, 31 compressed 100 controldict 16, 99, 100 converttometers 38 corrected 85 CrankNicolson 82 createpatch 25, 55 createturbulencefields 102 CrossPowerLaw 60 cyclic 40, 41, 70, 71 cyclicami 40, 41, 71 D decomposepar 30, 103, 104 deltat 100, 117, 126 diagonal 89, 91 DICGaussSeidel 90 DILU 89 91 dimensions 72 Display 199, 206 distance 54 divschemes 83, 115, 124 dynoneeqeddy 69 E edgegrading 39, 40 edges 40 Edit View Options 206
221 empty 40, 41, 71, 74 endtime 98, 99, 117, 126 equations 94 Euler 82, 124 evince 28 execflowfunctionobjects 107 expansionratio 51 F facelimited 83 facemdlimited 83 faces 41 FDIC 89, 90 featureangle 51 features 49, 54 fields 94 finitevolume 11 firsttime 98 fixed 100 fixedgradient 132 fixedvalue 73, 75, 79, 132 fluent3dmeshtofoam 34 fluentmeshtofoam 10, 34 foamjob 27, 31, 104 foamlog 27 foamtovtk 57 fvschemes 16, 101 fvsolution 16 G GAMG 58, 89, 90, 139 GaussSeidel 90, 91 general 100 General 199 geometry 47, 53 Ghostview 28 gnuplot 104 gradschemes 85 H hepsithermo 62 herhothermo 62 HerschelBulkley 60 homogeneousdynsmagorinsky 69 I icofoam 9 ideasunvtofoam 36 Information 43 inside 54 interfoam 9 internalfield 101, 113, 114 interpolation 135 K kepsilon 67, 102 komega 67 komegasst 67 L laminar 67 150
222 laplacianfoam 9 laplacianschemes 85, 101 Last Frame 28, 31 latesttime 98, 102, 104, 106 layers 51 LESModel 67, 69 LESProperties 16 levels 49, 54 limitedlinear 83, 84 limitedlinearv 84 linearupwind 83 linearupwindv 84 locationinmesh 50 M maxco 100 maxglobalcells 49 mergepatchpairs 42 method 103 mixture 62 mode 54 multicomponentmixture 62 multiphaseinterfoam 9 MUSCL 83 N neighbourpatch 40, 71 Newtonian 60 nextwrite 99 nowritenow 99 nsmoothpatch 50 numberofsubdomains 103 O oneeqeddy 69 Opacity 199 outside 54 P parafoam 22, 41, 43, 57, 74, 106, 197 ParaView 13, 22, 26, 28, 31, 41, 43, 57, 74, 106, 197, 207, 217 patchintegrate 105, 106 PBiCG 89, 91 PCG 89 pimplefoam 9, 67, 82 PIMPLE 9, 93 Pipeline Browser 197, 200 pisofoam 9 PISO 91, 93, 95, 97, 126, 137 Plot Over Line 206 polyline 40 polymesh 16 potentialfoam 9, 101 powerlaw 60 preconditioner 89, 91, 139 prefcell 132 prefpoint 94, 132 prefvalue 94, 132 Properties 26, 28, 29, 31, 199 puremixture 62
223 purgewrite 99, 117 Q QUICK 83, 86 147 R RASModel 67 RASProperties 16 reactingmixture 62 realizableke 67 reconstructpar 104, 106 refinementregions 53 refinementsurfaces 49 reltol 89, 92, 97 Representation 26, 199 renumbermesh 57, 58 RNGkEpsilon 67 runtime 99 runtimemodifiable 99 S Save Animation 32, 207 scalartransportfoam 9 scientific 100 scotch 103 Show Orientation Axes 26 simple 103 simplecoeffs 103 simplefoam 9, 21, 22, 25, 27, 30, 57, 59, 66, 71, 82, 102, 122 124 simplegrading 39, 40 SIMPLE 9, 93 simulationtype 67 Slice 29 slip 79, 133 Smagorinsky 69 smoother 90 smoothsolver 89 91 snap 47, 50 snappyhexmesh 10, 12, 25, 34, 36, 42, 43, 46, 47, 52, 55, 208, 216 Solid Color 26 199, 201 solver 89, 91 specie 62 spline 40 Spreadsheet View 205 startfrom 98 starttime 98 steadystate 82, 115 stopat 99 Styling 199 Surface With Edges 26 surfaceconvert 217 surfacefeatureextract 25, 46, 49 symmetryplane 41, 71, 74 T thermophysicalproperties 16 thermotype 62 timeprecision 100 timestep 99
224 tolerance 89, 97 139 transformpoints 57 transportmodel 60 transportproperties 16, 59, 60 turbulence 150, 156 turbulenceproperties 16 type 40, 74, 123 U UMIST 83 148 uncompressed 100 uncorrected 85 V value 74 vanleer 83 vertices 38 View 205 W wall 41, 71 wallheatflux 106 wallheattransfer 80, 133 wedge 41, 71 writecontrol 99 writeformat 100 writenow 99 Z zerogradient 73 75, 79, 132 13 PIMPLE 9, 93 PISO 9, 93, 95, 97, 126, 137 SIMPLE 9, 93 Pipeline Browser 197, 200 Coloring 26, 28, 31 Display 199, 206 General 199 Information 43 Properties 26, 28, 29, 31, 199 Representation 26, 199 Slice 29 Solid Color 26 Styling 199 Surface With Edges 26 146 39, 40, 51 calculated 120 cyclic 40, 41, 70, 71 cyclicami 40, 41, 71 empty 40, 41, 71, 74 fixedgradient 132 fixedvalue 73, 75, 79, 132 slip 79, 133 symmetryplane 41, 71, 74 wall 41, 71 wallheattransfer 80, 133 wedge 41, 71 zerogradient 73 75, 79, 132 adjusttimestep 126 castellatedmeshcontrols 47, 50, 53 class 72 converttometers 38 deltat 100, 117, 126 dimensions 72 divschemes 83, 115, 124 edgegrading 39, 40 edges 40 equations 94 expansionratio 51 featureangle 51 features 49, 54 fields 94 geometry 47, 53 gradschemes 85 internalfield 101, 113, 114 laplacianschemes 85, 101 layers 51 LESModel 67, 69 levels 49, 54 locationinmesh 50 maxco 100
225 maxglobalcells 49 mergepatchpairs 42 method 103 mode 54 neighbourpatch 40, 71 nsmoothpatch 50 numberofsubdomains 103 preconditioner 89, 91, 139 prefcell 132 prefpoint 94, 132 prefvalue 94, 132 purgewrite 99, 117 RASModel 67 refinementregions 53 refinementsurfaces 49 reltol 89, 92, 97 runtimemodifiable 99 simplegrading 39, 40 simulationtype 67 smoother 90 snap 47, 50 solver 89, 91 specie 62 spline 40 startfrom 98 stopat 99 thermotype 62 timeprecision 100 tolerance 89, 97 transportmodel 60 type 40, 74, 123 value 74 vertices 38 writecontrol 99 writeformat 100 adjustableruntime 99, 100 ascii 100 backward 82 binary 100 bounded 84 celllimited 83 cellmdlimited 83 compressed 100 corrected 85 CrankNicolson 82 diagonal 89, 91 DICGaussSeidel 90 DILU 89 91 distance 54 endtime 98, 99, 117, 126 Euler 82, 124 facelimited 83 facemdlimited 83 FDIC 89, 90 firsttime 98 fixed 100 GAMG 58, 89, 90, 139 GaussSeidel 90, 91 general 100 inside 54 laminar 67 latesttime 98, 102, 104, 106 limitedlinear 83, 84 limitedlinearv 84 linearupwind 83 linearupwindv 84 MUSCL 83 nextwrite 99 nowritenow 99 outside 54 PBiCG 89, 91 PCG 89 polyline 40 QUICK 83, 86 runtime 99 scientific 100 scotch 103 simple 103 simplecoeffs 103 smoothsolver 89 91 starttime 98 steadystate 82, 115 timestep 99 UMIST 83 uncompressed 100 uncorrected 85 vanleer 83 writenow 99 17 ParaView13, 22, 26, 28, 31, 41, 43, 57, 74, 106, 197, 207, 217 39, 40 92, 97, 99, 137, 140, 142 41, 71 foamjob 27, 31, 104 foamlog 27 12 buoyantboussinesqpimplefoam 9 buoyantboussinesqsimplefoam 9, 60, 71, 80 buoyantpimplefoam 9, 25, 30, 117 buoyantsimplefoam 9, 61, 65, 71, 78, 79 icofoam 9 interfoam 9 laplacianfoam 9 multiphaseinterfoam 9 pimplefoam 9, 67, 82
226 pisofoam 9 potentialfoam 9, 101 scalartransportfoam 9 simplefoam 9, 21, 22, 25, 27, 30, 57, 59, 66, 71, 82, 102, 122 124 evince 28 Ghostview 28 gnuplot 104 blockmeshdict 34, 41, 45, 46 boundary 40, 73 controldict 16, 99, 100 faces 41 fvschemes 16, 101 fvsolution 16 LESProperties 16 RASProperties 16 thermophysicalproperties 16 transportproperties 16, 59, 60 turbulenceproperties 16 polymesh 16 Opacity 199 97, 139 12 181, 193, 194 13 Apply 29, 197, 200, 202 206 Choose Preset 201 Last Frame 28, 31 Show Orientation Axes 26 12 Edit View Options 206 Spreadsheet View 205 View 205 Plot Over Line 206 Save Animation 32, 207 Solid Color 199, 201 12 BirdCarreau 60 CrossPowerLaw 60 dynoneeqeddy 69 hepsithermo 62 herhothermo 62 HerschelBulkley 60 homogeneousdynsmagorinsky 69 kepsilon 67, 102 komega 67 komegasst 67 multicomponentmixture 62 Newtonian 60 oneeqeddy 69 powerlaw 60 puremixture 62 reactingmixture 62 realizableke 67 RNGkEpsilon 67 Smagorinsky 69 blockmesh 10, 12, 22, 25, 34, 36, 41 43, 46, 53, 55 checkmesh 10, 52, 56, 85, 94 createpatch 25, 55 createturbulencefields 102 decomposepar 30, 103, 104 execflowfunctionobjects 107 fluent3dmeshtofoam 34 fluentmeshtofoam 10, 34 foamtovtk 57 ideasunvtofoam 36 parafoam 22, 41, 43, 57, 74, 106, 197 patchintegrate 105, 106 reconstructpar 104, 106 renumbermesh 57, 58 snappyhexmesh 10, 12, 25, 34, 36, 42, 43, 46, 47, 52, 55, 208, 216 surfaceconvert 217 surfacefeatureextract 25, 46, 49 transformpoints 57 wallheatflux 106 1, 8, 11, 36 finitevolume 11 mixture 62 16, 21, 66, 67, 69, 71, 102, 131, 132, 150, 156 bounded 146 interpolation 135 laminar 150 QUICK 147 tolerance 139 turbulence 150, 156 UMIST 148
2004 PENGUINITIS http://www.geocities.jp/penguinitis2002/ 1995 1998 1998 1998 2007 2007 2012 OCAEL CAE 160-0007 5-14 2F Tel : 03-3350-0363, Fax : 03-3341-1830 Email : office@opencae.or.jp URL : http://www.opencae.or.jp/ () () OpenFOAM c CAE 2016 2016 6 20 1 1 CAE 1 4 11 102 0071 03 3265 8341FAX 03 3264 8709 http://www.morikita.co.jp/ () PrintedinJapanISBN978 4 627 69101 8