22.312 ENGINEERING OF NUCLEAR REACTORS - MIT · PDF file22.312 ENGINEERING OF NUCLEAR REACTORS...

5
22.312 ENGINEERING OF NUCLEAR REACTORS Monday, December 14 th , 2009, 9:00am-12:00 pm J. Buongiorno, 2009 OPEN BOOK FINAL EXAM 3 HOURS Problem 1 (15%) – Sizing the shell of a spherical containment i) The principal stresses for a thin spherical shell are: σ r = -(p i +p o )/2 (1) σ θ = σ ϕ = (p i -p o )R c /(2 t c ) where p i =1.9 MPa and p o =0.1 MPa, R c =12.5 m and t c is the shell thickness. Hook’s law yields: ε θ =u/R c =1/E[σ θ - ν (σ ϕ +σ r )] (2) where E = 184 GPa and ν = 0.33. Substituting Eq. (1) in Eq. (2), setting u=1 cm and solving for t c , one gets: t c =R c (1-ν) (p i -p o )/[2E u/R c -ν(p i +p o )]3.7 cm Since R c /t c >10, the thin shell asumption of accurate. ii) The primary membrane general stress intensity for this case is: P m =(σ θ - σ r )102 MPa σ θ and σ r were calculated from Eq. (1) (thin shell assumption still applies), for t c =8 cm. The ASME limit is S m =110 MPa. Therefore the margin is S m /P m 1.075, or 7.5%. Problem 2 (25%) – Reduction of containment pressure after LOCA Conservation of energy for the containment: ss decay CV Q Q t E & & = (3) where 2 . 0 0 066 . 0 = t Q Q decay & & , 0 Q & =1000 MW, and ss Q & =20 MW. Integrating Eq. (3): 1

Transcript of 22.312 ENGINEERING OF NUCLEAR REACTORS - MIT · PDF file22.312 ENGINEERING OF NUCLEAR REACTORS...

22.312 ENGINEERING OF NUCLEAR REACTORS

Monday, December 14th, 2009, 9:00am-12:00 pm

J. Buongiorno, 2009

OPEN BOOK FINAL EXAM 3 HOURS

Problem 1 (15%) – Sizing the shell of a spherical containment

i) The principal stresses for a thin spherical shell are:

σr = -(pi+po)/2 (1) σθ = σϕ = (pi-po)Rc/(2 tc)

where pi=1.9 MPa and po=0.1 MPa, Rc=12.5 m and tc is the shell thickness. Hook’s law yields: εθ=u/Rc=1/E[σθ - ν (σϕ+σr)] (2) where E = 184 GPa and ν = 0.33. Substituting Eq. (1) in Eq. (2), setting u=1 cm and solving for tc, one gets: tc=Rc(1-ν) (pi-po)/[2E u/Rc-ν(pi+po)]≈3.7 cm Since Rc/tc>10, the thin shell asumption of accurate. ii) The primary membrane general stress intensity for this case is:

Pm=(σθ - σr)≈102 MPa

σθ and σr were calculated from Eq. (1) (thin shell assumption still applies), for tc=8 cm. The

ASME limit is Sm=110 MPa. Therefore the margin is Sm/Pm≈1.075, or 7.5%.

Problem 2 (25%) – Reduction of containment pressure after LOCA

Conservation of energy for the containment:

ssdecayCV QQt

E && −=∂

∂ (3)

where 2.0

0 066.0 −= tQQdecay&& , 0Q& =1000 MW, and ssQ& =20 MW. Integrating Eq. (3):

1

J. Buongiorno, 2009

− = & 0.066E2 E1 Q 0.80 t

0.8 2 −Q& sst2 (4)

Expanding the left-hand side, one gets:

0.066M acv,a (T2 −T1) +M w{[u f (T2 )(1− x2 ) + ug (T2 )x2 ]−[u f (T1)(1− x + u 01) g (T1)x1]} = Q&0 t .8 −Q& t

0.8 2 ss 2

(5) where Ma, cv,a, T1, Mw and x1 are all known from the problem statement. The following equation holds for the control volume: Vtot = M w[v f (T2 )(1− x2 ) + vg (T2 )x2 ] (6) The containment pressure at t2, P2=0.5 MPa, is the sum of the partial pressures of water and air:

M R TP P (T ) 22 = 2 + a a

sat (7) Vtot −M w(1− x2 )v f (T2 )

Therefore, Eqs. (5), (6) and (7) are 3 equations in the only unknowns t2, T2 and x2. Actually solving the equations, one finds t2≈14300 s, T2≈140.4 °C and x2≈0.035. Problem 3 (45%) – Superheated Boiling Water Reactor

i) T-s diagram:

ii) Taking the whole system as a control volume, the conservation of energy yields: 0 = Q& +m& FW (hFW − hsup ) ⇒ m& = Q&FW /(hsup − hFW ) (8)

T 6.0 MPa

s

E F

A B

C=D • • •

2

where Q& =1000 MW and hFW and hsup are the specific enthalpy of the feedwater and superheated steam, respectively. The difference hsup-hFW can be expressed as follows: hsup hFW cp,g (Tsup Tsat ) hfg cp, f (Tsat TFW )

J. Buongiorno, 2009

−++−=− ≈2936 kJ/kg where TFW= 230°C and Tsup=510°C. Therefore, Eq. (8) yields FWm& ≈340.6 kg/s. iii) The acceleration pressure drop is

]11[,,

2++ −=Δ

inmoutmacc GP

ρρ (9)

where AmG&

= ≈1800 kg/m2s, m& =2270 kg/s, A=1.26 m2 and

fg

m xxρααρ

ρ

)1()1(

122

−−

+≡+ (10)

Since at the inlet there is only the liquid phase, it is finm ρρ =+, , while at the outlet x=0.15 and the

void fraction can be found from the fundamental relation of two-phase flow:

Sx

x

f

g ⋅−

⋅+=

11

1

ρρα ≈0.69

where S=2, per the problem statement. Then it is +outm,ρ ≈240.3 kg/m3 from Eq. (10), and finally

Eq. (9) yields ΔPacc≈ 9,235 Pa

iv) Since the heat flux is axially constant, dryout would occur at the outlet (Point B). The critical quality at the outlet is found to be xcr≈0.344 from the CISE-4 correlation with Lb=3 m, and the coefficients a=0.5987 and b=2.2255, calculated for P=6 MPa, Pc=22.1 MPa, G=1800 kg/m2s>G*=1211 kg/m2s, De=0.02 m. Then the critical power of the A→B channels is ])([ ,, fgcrAsatfpABcr hxTTcmQ +−= && ≈1311 MW,

where TA= 268°C. So, the AB

ABcr

QQ

CPR &

&,= ≈2.12, with ])([ , fgBAsatfpAB hxTTcmQ +−= && ≈618 MW

being the operating power of the A→B channels, where xB=0.15. Problem 4 (15%) – Thermodynamic analysis of a new power cycle To be thermodynamically feasible, the cycle must not violate the 1st and 2nd law of thermodynamics.

3

Taking the whole power cycle as the control volume, the conservation of energy (1st law) becomes: 0 = Q& −W& +m& sea (hin − hout ) (11) where steady-state was assumed and Q& =1000 MW, W& =400 MW, m& sea =15000 kg/s, (hin − hout ) = csea (Tin −Tout ) , csea =4000 J/kg°C and Tin =288 K (15°C) and Tout =298 K (25°C). Using these numbers, Eq. 11 is identically satisfied. Therefore, the cycle does not violate the 1st law. With the same choice of control volume, the 2nd law becomes:

J. Buongiorno, 2009

genoutinsear

SssmTQ &&&

+−+= )(0 ⇒ r

inoutseagen TQssmS&

&& −−= )( (12)

where Tr=723 K (450°C) and in

outseainout T

Tcss ln=− . Then Eq. 12 yields genS& = 665 kW/K > 0,

therefore the cycle does not violate the 2nd law either.

4

MIT OpenCourseWarehttp://ocw.mit.edu

22.312 Engineering of Nuclear ReactorsFall 2015

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.