This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: ENU 4134 Project #2, Fall 2011 D. Schubring 1 Development & Use of Simple SingleChannel Analysis Code DUE: Part A (November 23, graded in order received, early turnin advised), Presentations (November 30th and December 2nd), Part B (December 7) ALLOWED COLLABORATION: You may collaborate and divide labor within your group however you choose. One set of deliverables is required per group. No collaboration outside of your group is permitted. 1 Assignment Develop a singlechannel analysis code using the numerical method derived/discussed in class, with the following features, in the programming/scripting language of your choice (MATLAB ac cepted, TK/EES not accepted): User selectable values: L , L e , D , Pitch , T m,in , m , P nominal , q max , type of reactor (BWR vs. PWR, for CHF limit applied), D ci , D fo , subcooled boiling flag. These user selectable values should be in a seperate file (an input file) as the main program. Finite volume size: 400 finite volumes along the axial length. Power shape: cosine, with extrapolated length. Evaluation of properties: use the digitized version of T&Ks Table E1 (on the website) and lineraly interpolate. Note : for density, linearly interpolate the values of specific volume and invert, do not invert these values and interpolate that result. The computed will be different and so will your final answers. Convection to coolant: use the Weisman correlation for singlephase calculations, with the Schrock and Grossman correlation (based on Weismans htc lo ) for twophase. Subcooled boiling model: compute x e from an energy balance and use Equation 1222 in T&K to estimate x . Take Z D as the first location where T co > T sat ( i.e. , neglect the distinctions among bubble departure, boiling incipience, etc.). Use x , not x e , in htc 2 , CHF, and dp/dz correlations. When the subcooled boiling model flag is off, x = x e . CHF: use the correlation of Bowring, adapted by multiplying through by the same as for Weisman. Compute the local CHFR (termed DNBR in a PWR) whenever any vaporization has occurred in the coolant. Alert the user when this is below 1.3 at any point for a PWR and 1.9 at any point for a BWR. Cladding: assume a constant conductivity k c = 15 W m 1 K 1 . Gap conductance: use Equation 8107a, with f = c = 1 and eff = R ci R fo . This will require iteration; select an appropriate convergence criterion on T fo . Assume a heliumfilled gap. Fuel conductivity: integrate Equation 816a and use a conductivity integral approach. This will require iteration; select an appropriate convergence criterion on T max . Pressure drop: in each finite volume, estimate this using the Cheng and Todreas equation for singlephase friction (if singlephase) or the HEM model if twophase. dx/dz comes from the ENU 4134 Project #2, Fall 2011 D. Schubring 2 energy balance. In twophase, take Cheng and Todreas for f lo and assume f lo = f TP . This is....
View
Full
Document
This note was uploaded on 10/22/2011 for the course ENU 4134 taught by Professor Schubring during the Fall '11 term at University of Florida.
 Fall '11
 Schubring

Click to edit the document details