境界条件の設定
境界条件0フォルダも以下のように領域ごとに作成する必要があります。
あとでもとに戻せるように0.origの中で設定ファイルを作成しましょう。
主要なファイルだけ記載しておきます。
0/coldwater/U
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | dimensions [ 0 1 -1 0 0 0 0 ]; internalField uniform (0.008 0 0); boundaryField { #includeEtc "caseDicts/setConstraintTypes"; hotwater_to_solid { type fixedValue; value uniform (0 0 0); } hotwaterXmax { type pressureInletOutletVelocity; value uniform (0.008 0 0); } hotwaterXmin { type fixedValue; value uniform (0.008 0 0); } /* hotwaterYmin { type noSlip; } */ Zmin { type wedge; } Zmax { type wedge; } } |
0/coldwater/p_rgh
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | dimensions [ 1 -1 -2 0 0 0 0 ]; internalField uniform 101325; boundaryField { #includeEtc "caseDicts/setConstraintTypes"; hotwater_to_solid { type fixedFluxPressure; value uniform 101325; } hotwaterXmin { type fixedFluxPressure; value uniform 101325; } hotwaterXmax { type fixedValue; value uniform 101325; } /* hotwaterYmin { type fixedValue; value uniform 101325; } */ Zmin { type wedge; } Zmax { type wedge; } } |
0/coldwater/T
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | dimensions [ 0 0 0 1 0 0 0 ]; internalField uniform 293; boundaryField { #includeEtc "caseDicts/setConstraintTypes"; hotwater_to_solid { type compressible::turbulentTemperatureCoupledBaffleMixed; value uniform 293; Tnbr T; kappaMethod fluidThermo; kappa none; } hotwaterXmax { type zeroGradient; } hotwaterXmin { type fixedValue; value uniform 313.15; } /* hotwaterYmin { type zeroGradient; } */ Zmin { type wedge; } Zmax { type wedge; } } |
0/coldwater/U
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | dimensions [ 0 1 -1 0 0 0 0 ]; internalField uniform (0.008 0 0); boundaryField { #includeEtc "caseDicts/setConstraintTypes"; coldwater_to_solid { type fixedValue; value uniform (0 0 0); } coldwaterXmax { type pressureInletOutletVelocity; value uniform (0.008 0 0); } coldwaterXmin { type fixedValue; value uniform (0.008 0 0); } coldwaterYmax { type noSlip; } Zmin { type wedge; } Zmax { type wedge; } } |
0/coldwater/p_rgh
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | dimensions [ 1 -1 -2 0 0 0 0 ]; internalField uniform 101325; boundaryField { #includeEtc "caseDicts/setConstraintTypes"; coldwater_to_solid { type fixedFluxPressure; value uniform 101325; } coldwaterXmax { type fixedValue; value uniform 101325; } coldwaterXmin { type fixedFluxPressure; value uniform 101325; } coldwaterYmax { type fixedFluxPressure; value uniform 101325; } Zmin { type wedge; } Zmax { type wedge; } } |
0/coldwater/T
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | dimensions [ 0 0 0 1 0 0 0 ]; internalField uniform 293; boundaryField { #includeEtc "caseDicts/setConstraintTypes"; coldwater_to_solid { type compressible::turbulentTemperatureCoupledBaffleMixed; value uniform 293; Tnbr T; kappaMethod fluidThermo; kappa none; } coldwaterXmax { type zeroGradient; } coldwaterXmin { type fixedValue; value uniform 293.15; } coldwaterYmax { type zeroGradient; } Zmin { type wedge; } Zmax { type wedge; } } |
0/solid/T
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 | dimensions [ 0 0 0 1 0 0 0 ]; internalField uniform 293; boundaryField { #includeEtc "caseDicts/setConstraintTypes"; solid_to_hotwater { type compressible::turbulentTemperatureCoupledBaffleMixed; Tnbr T; kappaMethod solidThermo; kappa none; value uniform 293; } solid_to_coldwater { type compressible::turbulentTemperatureCoupledBaffleMixed; Tnbr T; kappaMethod solidThermo; kappa none; value uniform 293; } solidXmin { type zeroGradient; } solidXmax { type zeroGradient; } Zmin { type wedge; } Zmax { type wedge; } } |
軸対称の面は「type wedge;」を指定しています。
物性値の変更
constantの中の物性値情報を以下のように書き換えておきます。
hotwater(coldwater)/thermophysicalProperties
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 18; } equationOfState { rho 1000; } thermodynamics { Cp 4181; Hf 0; } transport { mu 959e-6; Pr 6.62; } } |
solid/thermophysicalProperties
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 | thermoType { type heSolidThermo; mixture pureMixture; transport constIso; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 27; //[g/mol] } thermodynamics { Cp 903; Hf 0; } transport { kappa 237; } equationOfState { rho 2702; } } |
あとはふく射については考慮しないので以下のようにオフにしておきます。
hotwater(coldwater solid)/radiationProperties
1 2 3 | radiation off; radiationModel none; |
重力はyの負の方向にかけておきます。
constant/g
1 2 | dimensions [0 1 -2 0 0 0 0]; value (0 -9.81 0); |
さらに領域の指定をします。
constant/regionProperties
1 2 3 4 5 | regions ( fluid (hotwater coldwater) solid (solid) ); |
Follow @t_kun_kamakiri