Module 13: Critical Flow Phenomenon Joseph S. Miller, PE
Total Page:16
File Type:pdf, Size:1020Kb
Fundamentals of Nuclear Engineering Module 13: Critical Flow Phenomenon Joseph S. Miller, P.E. 1 2 Objectives: Previous Lectures described single and two-phase fluid flow in various systems. This lecture: 1. Describe Critical Flow – What is it 2. Describe Single Phase Critical Flow 3. Describe Two-Phase Critical Flow 4. Describe Situations Where Critical Flow is Important 5. Describe origins and use of Some Critical Flow Correlations 6. Describe Some Testing that has been Performed for break flow and system performance 3 Critical Flow Phenomenon 4 1. What is Critical Flow? • Envision 2 volumes at different pressures suddenly connected • Critical flow (“choked flow”) involves situation where fluid moves from higher pressure volume at speed limited only by speed of sound for fluid – such as LOCA • Various models exist to describe this limiting flow rate: • One-phase vapor, one phase liquid, subcooled flashing liquid, saturated flashing liquid, and two-phase conditions 5 2. Critical Single Phase Flow Three Examples will be given 1.Steam Flow 2.Ideal Gas 3.Incompressible Liquid 6 Critical Single Phase Flow - Steam • In single phase flow: sonic velocity a and critical mass flow are directly related: g v(P,T)2 a2 = c dv(P,T) dP S 2 g G = a2 ρ(P,T )2 = c crit dv(P,T ) dP S • Derivative term is total derivative of specific volume evaluated at constant entropy • Tabulated values of critical steam flow can be found in steam tables 7 Example Critical Steam Flow Calculation • Assume 2 in2 steam relief valve opens at 1000 psi • What is steam mass discharge rate? • Assume saturated system with Tsat = 544.61°F • f = 50.3 • Wcrit = f P A =(50.3lb-m/hr)(1000 psi)(2 in2) = 100,600 lb-m/hr = 27.94 lb-m/sec. 8 Mass flow rate of a gas at choked conditions for Ideal Gas All gases flow from upstream higher stagnation pressure sources to downstream lower pressure sources. There are several situations in which choked flow occurs, such as: change of cross section (as in a convergent- divergent nozzle or flow through an orifice plate. When the gas velocity is choked, the equation for the mass flow rate in SI metric units is: m = or for an ideal gas m = where: m= mass flow rate, kg/s C = discharge coefficient, dimensionless (usually about 0.72) A = discharge hole cross-sectional area, m² k= cp/cv of the gas cp = specific heat of the gas at constant pressure cv = specific heat of the gas at constant volume ρ = real gas density at P and T, kg/m³ P= absolute upstream stagnation pressure, Pa M= the gas molecular mass, kg/kmole (also known as the molecular weight) R= Universal gas law constant = 8314.5 (N·m) / (kmole·K) T= absolute gas temperature, K Z= the gas compressibility factor at P and T, dimensionless 9 Critical Single Phase Flow – Incompressible Liquid (Lahey and Moody, Ref. 1) Where throat pressure, pr, is equal to the downstream receiver pressure, pR, and Gl, is the so-called Bernoulli mass flux. 10 Critical Single Phase Flow – Ideal Gas(Lahey and Moody, Ref. 1) (cont) 11 Ideal Gas Flow Rate 12 3. Two-phase Critical Flow Using Moody Equilibrium Flow Model (Refs. 1 & 2) 13 Two-phase Critical Flow Using Moody Equilibrium Flow Model (Refs. 1 & 2) 14 Two-phase Critical Flow Using Moody Equilibrium Flow Model (Refs. 1 & 2) (cont) That is, for a known stagnation state, local static pressure, p, and slip ratio, S, the flow rate unit area is uniquely determined. From numerous experimental and theoretical studies, the slip 1/2 ratio, S, has been determined to range from S = 1.0 (homogeneous) to S = (ρf / ρg) . For condition imposed by Moody, the slip ratio is evaluated as 1/3 S = S(p) = (ρf / ρg) . (9.53) Equation (9.53) shows that at a maximum G, the slip ration is a function of pressure only. Employing Eq. (9.53) in Eq. (9.52) and using saturated steam-water properties, Fig. 9-10a can be obtained, which gives Gc in terms of p0 and h0 . 15 Two-phase Critical Flow Using Moody Equilibrium Flow Model (Refs. 1 & 2) (cont) 16 Example Calculation From page 8 for sat steam at 1000 psia: Wcrit = f P A =(50.3lb-m/hr)(1000 psi)(2 in2) = 100,600 lb-m/hr = 27.94 lb-m/sec. What is this value when using Moody table on page 16? 2 Gc(P = 1000, h = 1193) = 2000 lbm/sec-ft 2 2 2 Wcrit/A = (27.94 lbm/sec)/ (2 in / 144 in /ft ) 2 Wcrit/A = 2012 lbm/sec-ft Checks 17 Two-phase Critical Flow Using Moody Homogeneous Equilibrium Model (Ref. 1 & 2) (cont) 18 Two-phase Critical Flow Comparing Moody, HEM and Bernoulli (Refs. 1 & 2) (cont) 19 Henry-Fauske/ Homogeneous Equilibrium Critical Flow Transition Model 20 Brunell Critical Flow for Flashing Water • As with anything dealing with two phase flow: • Assumed Slip Ratio affects two-phase critical flow correlations. • For example: Burnell developed critical flow correlation for flow of flashing water in short pipes and orifices 1/2 Gcrit = (2gcρl(P – (1-C)Psat)) • Other critical flow correlations exist for different assumed regimes 21 4. Where is Critical Flow Important • Used in LOCA break flow calculations • Main Steam Line Breaks • Maximum flow in a orifice • Safety Valve Opening LOCA break flow calculations • Small Breaks • Large Breaks • Intermediate Breaks • Include Break side and discharge side, especially to determine flow into a containment or sub compartment Break Flow from Discharge Pipe for double-ended Rupture Reference: NEDO 24548 ASSUMPTIONS: •The initial Velocity of the fluid in the pipe is zero. When considering both sides of the break, the effects of initial velocities would tend to cancel out. •Constant reservoir pressure at 1000 psia at saturated liquid conditions. •Initially fluid conditions inside the pipe on both sides of the break are similar. • Wall thickness of the pipe is small compared to the diameter. •Subcompartment pressure ≈ 0 psig. • Quasi-steady mass flux is calculated using the Moody steady slip flow model with subcooling. Nomenclature and Definitions ABR - Break Area, ft AL - Minimum cross-sectional area between the vessel and the break. This is the sum of the area of parallel flow paths. c1 - Sonic speed in the fluid D - Pipe inside diameter at the break location, ft FI - Inventory flow multiplier FI = 0.75 for saturated steam FI = 0.50 for liquid 2 gc - Proportionality constant (= 32.17 lbm-ft/lbf-sec ) 2 G1 - Mass Flux, lbm/ft -sec 2 Gc - Maximum mass flux, lbm/ft -sec ho - Reservoir or vessel enthalpy, BTU/lbm hp - Initial enthalpy of the fluid in the pipe, BTU/lbm LI - Inventory length. The distance between the break and the nearest area decrease of AL . - Mass Flowrate, lbm/sec - Mass flowrate during the inventory period, lbm/sec Po - Reservoir or vessel pressure, 1000 psia Psat - Saturation pressure for liquid with an enthalpy of hp t - Time, sec tI - Length of the inventory period, sec v - specific volume of the fluid initially in the pipe, ft 3/ lbm 3 VI - Volume of the pipe between the break and AL , ft X - Separation distance of the break, ft Sonic Speed in Water versus Water Pressure Moody Equilibrium Critical Flow with Slip Extended into Subcooled Region Short Term Release Rate Calculations Short Term Release Rate Calculations Btu h := 542.6 From ASME Steam Tables f lbm ho := hf FI := 0.50 for liquid from NEDO-24548 pg.11 D := 2.5 ft Inside diameter of pipe at break down location (Assumed) DL := 0.5 ft Diameter of restriction (Assumed) LI := 26 ft The distance between break and nearest restriction Po := 1000 psi Initial pressure Short Term Release Rate Calculations (cont) lbm G := 7500 Maximum mass flux from Moody Plots c 2 ft G1 := Gc mass flux from Moody Plots 2 ft c := 1.6⋅10 sonic speed in the fluid (from Sonic Speed in Water vs. Water Pressure Figure 1 s 3. ft v:= 0.02159 specific volume of fluid initially in the pipe from ASME Steam Tables lbm 1 2 ABR := π⋅D 2 4 ABR = 4.909 ft Break Area 1 2 2 A := π⋅D A = 0.196 ft Mimimum cross sectional area between vessel and break L 4 L L 3 VI := ABR⋅LI VI = 127.627 ft Volume of the pipe between the break and AL Short Term Release Rate Calculations (cont) AL = 0.04 FI = 0.5 ABR AL If > FI ABR (2⋅LI) tI := tI = 0.325 sec c1 AL If < FI ABR VI t := I. t = 0.321 sec ABR⋅G1⋅FI⋅v I. Short Term Release Rate Calculations (cont) lb 4 MI := G1⋅ABR⋅FI MI = 1.841 × 10 sec Flow during tI 3 lb M := A ⋅G M = 1.473 × 10 Steady state flow through ss L 1 ss sec restriction after tI 4 2 .10 4 1.5 .10 4 M 1 .10 5000 Flow (lb/sec) Discharge Pipe In 0 0 0.2 0.4 0.6 0.8 tI Time after Break (sec) 31 5. Describe Origins and Use of Critical Flow Correlations 32 5. Describe Origins and Use of Critical Flow Correlations (cont.) 33 5. Describe Origins and Use of Critical Flow Correlations (cont.) 5. Describe Origins and Use of Critical Flow Correlations (cont.) 35 5. Describe Origins and Use of Critical Flow Correlations (cont.) Use of Critical Flow Correlations in Nuclear Applications • Most common use of HEM, Moody, and Henry Fauske • Used in LOCA break flow calculations • Main Steam Line Breaks • Maximum flow in a orifice • Safety Valve Opening Use of Critical Flow Correlations in Nuclear Applications (cont.) 10CFR50 Appendix K--ECCS Evaluation Models I.