240 likes | 348 Views
Natural convection flow in a cask storage ware house For Toyo Engineering. PHOENICS Version : 3.5.1 DETAILS : BFC grid Three-Dimensional steady flow with heat transfer Buoyancy-influenced flow Surface to surface radiation included NX*NY*NZ= 288*74*42=895,104 NOTES :
E N D
Natural convection flow in a cask storage warehouseFor Toyo Engineering • PHOENICS Version : 3.5.1 • DETAILS : • BFC grid • Three-Dimensional steady flow with heat transfer • Buoyancy-influenced flow • Surface to surface radiation included • NX*NY*NZ=288*74*42=895,104 • NOTES : • The purpose of this model is to establish • the temperature distribution and the air flow rate .
Geometry x-y plane nx×nz=288×42 x-y plane nx×ny=288×74
Q A2 A1 A1 Q A2 outlet outlet outlet outlet External temperature 23.2℃ External temperature 23.2℃ Radiation model in the cask cask Air out Heat source : 22.6[kw] Air in
Velocity FlowRate[m3/s] per one cask Experimentaldata 0.28 0.3 Result of PHOENICS
Natural convection flow in a nuclear reactorFor Denchuken • PHOENICS Version : 3.6 • DETAILS : • Cylindrical-polar grid • Three-Dimensional steady or transient flow with heat transfer • Buoyancy-influenced flow • Surface to surface radiation included • NX*NY*NZ=40*45*80=144,000 • NOTES : • When an accident happens in the nuclear system, eg the cooling system may fail. Hence, heat will be released solely • by natural convection of air flow. The purpose of this model is to establish the distribution of temperature and the air flow rate under these circumstances.
Air outlet (connected to stack) Air inlet (connected to stack) PRACS IHX(Heat sink Q=-28.55[MW]) Argon gas zone Liquid Na convection zone Pump(11m3/min) Steady::on Transient::off Air convection zone Radial shield Reflector Core (Heat source Q=30[MW]) Structure of porous media zone Air Flow Liquid Na Flow
Boundary conditions Heat source and Heat sink: Steady:-28.55[MW] Transient:0[WM] Steady:30[MW] Transient: change profile
Pump power : Steady:11[m3/s] Transient: change profile.
Body force (Buoyancy): Liquid Na convection zone Air convection zone
Resistance force in the porous medium zone: where K: Pressure resistance coefficient. (This is different in each zone) ρ:density U:velocity
Q A2 A1 A1 Q A2 Radiation : 安全容器 コレクタ 原子炉容器 Q1 Q2 radiation
i Ai[m2] ki[-] n 1 1.0 0.3 2 2 1.0 0.25 2 POTOP 3 - 0.657 2 ΔPf9 4 5.28 1.3 1 5 3.48 0.76 1 ΔPf8 Δh Δhout 6 1.26 2.6 1 7 - 0.3174 2 ΔPf1 ΔPf2 Δhin ΔPf3 POUT PIN 8 1.0 1.5 2 ΔPf7 9 1.0 7.0 ΔPf6 2 ΔPf4 ΔPf5 Pressure value in the air convection zone boundary : The value of pressure in the air convection zone boundary is renewed while it is calculated from the next experience equation. i=6 to 10 i=1 to 5 PITOP
Maximum temperature [℃] Experimentaldata 550 540.6699 Result of PHOENICS Temperature distribution (steady)
Pressure drop [MPa] in Na convection zone Experimentaldata 2.5 2.51 Result of PHOENICS Pressure distribution (steady)
Air flow rate in the air convection zone. Maximum temperature in the liquid Na convection zone.
Transient simulation in a vapor turbine • PHOENICS Version : 3.5.0 • DETAILS : • Cartesian grid • Three-Dimensional transient flow • NX*NY*NZ=222×218×98=4,742,808
Air :Fixed temperature 下流側流入境界 Inlet of vapor Inlet of vapor Outlet Nozzle Box : Fixed temperature
温度出力点 ①(0°) 鉛直断面 ②(30°) ③(60°) ④(90°) 水平断面 ⑤(120°) ⑥(150°) ⑦(180°) SLICE POINT1 POINT2