Modification and application of water film model in COCOSYS for PWR's passive containment cooling

Abstract In this paper the physical model describing water film behaviors in German containment code system COCOSYS has been modified by taking into consideration the film breakup and subsequent phenomena as well as the effect of film interfacial shear stress created by countercurrent air flow. The modified model has extended its capability to predict particular water film behaviors including breakup at a critical film thickness based on minimum total energy criterion, the formation of rivulets according to total energy equilibrium as well as subsequent performance of rivulets according to several assumptions and observations from experiments. Furthermore, the modification considers also the change of velocity distribution on the cross section of film/rivulets due to shear stress. Based on the geometry of AP1000 and Generic Containment, simulations predicting containment pressure variation during accidents with operation of passive containment cooling system have been carried out. With the new model, considerably larger peak pressures are observed by comparing with those predicted with original water film model within a certain range of water film flow rate. Sensitivity analyses also point out that contact angle between water rivulets and steel substrate plays a significant role in the film cooling.

[1]  Mohamed S. El-Genk,et al.  On the breakup of a thin liquid film subject to interfacial shear , 2004, Journal of Fluid Mechanics.

[2]  R. P. Roy,et al.  A study of thin water film flow down an inclined plate without and with countercurrent air flow , 1989 .

[3]  D. G. Kröger,et al.  Countercurrent gas–liquid flow in inclined and vertical ducts — I: Flow patterns, pressure drop characteristics and flooding , 2000 .

[4]  Thomas J. Hanratty,et al.  The interfacial drag and the height of the wall layer in annular flows , 1976 .

[5]  Hans-Josef Allelein,et al.  Simulation of AP1000's passive containment cooling with the German Containment Code System COCOSYS , 2013 .

[6]  Jingqi Yuan,et al.  Portraying the Countercurrent Flow on Packings by Three-Dimensional Computational Fluid Dynamics Simulations , 2008 .

[7]  A. Doniec,et al.  Flow of a laminar liquid film down a vertical surface , 1988 .

[8]  C. Spengler,et al.  COCOSYS: Status of development and validation of the German containment code system , 2008 .

[9]  Nicola Forgione,et al.  Statistical Characteristics of a Water Film Falling Down a Flat Plate at Different Inclinations and Temperatures , 2002 .

[10]  S. Bankoff,et al.  INITIAL BREAKDOWN OF A HEATED LIQUID FILM IN COCURRENT TWO-COMPONENTANNULAR FLOW: II. RIVULET AND DRYPATCH MODELS‡ , 1980 .

[11]  Goon-Cherl Park,et al.  An experimental study on evaporative heat transfer coefficient and applications for passive cooling of AP600 steel containment , 2001 .

[12]  Spiros V. Paras,et al.  Counter-current gas-liquid flow in a vertical narrow channel—Liquid film characteristics and flooding phenomena , 2006 .

[13]  S. Maharudrayya,et al.  Pressure losses in laminar flow through serpentine channels in fuel cell stacks , 2004 .

[14]  Mohamed S. El-Genk,et al.  Minimum thickness of a flowing down liquid film on a vertical surface , 2001 .

[15]  Franz Mayinger,et al.  Experimental and analytical study of countercurrent flow limitation in vertical Gas/liquid flows , 1992 .

[16]  D. Bharathan,et al.  Air-water countercurrent annular flow , 1979 .

[17]  Michele Ciofalo,et al.  Natural convection cooling of a hot vertical wall wet by a falling liquid film , 2009 .

[18]  L. Hopf,et al.  Turbulenz bei einem Flusse , 1910 .

[19]  Andrzej Doniec,et al.  Laminar flow of a liquid rivulet down a vertical solid surface , 1991 .

[20]  黄鸿鼎,et al.  FLOODING IN VERTICAL TUBES , 1989 .

[21]  Xu Cheng,et al.  Experimental study of water film falling and spreading on a large vertical plate , 2012 .

[22]  A. R. Balakrishnan,et al.  Effect of tube diameter on flooding , 2001 .

[23]  J. Mikielewicz,et al.  Minimum thickness of a liquid film flowing vertically down a solid surface , 1976 .

[24]  G. Brunner,et al.  Film Thickness, Flow Regimes, and Flooding in Countercurrent Annular Flow of a Falling Film at High Pressures , 2001 .

[25]  T. L. Schulz,et al.  Westinghouse AP1000 advanced passive plant , 2006 .

[26]  L. B. Fore,et al.  Interfacial friction in gas–liquid annular flow: analogies to full and transition roughness , 1999 .

[27]  James A. Drallmeier,et al.  Developing Laminar Gravity-Driven Thin Liquid Film Flow Down an Inclined Plane , 2010 .

[28]  Y. Sudo.  Mechanism and Effects of Predominant Parameters Regarding Limitation of Falling Water in Vertical Countercurrent Two-Phase Flow , 1996 .

[29]  Luis E. Herranz,et al.  Experimental analysis of heat transfer within the AP600 containment under postulated accident conditions , 1998 .

[30]  G. Wallis.  One Dimensional Two-Phase Flow , 1969 .

[31]  W. Murgatroyd,et al.  CRITERIA FOR THE BREAK-UP OF THIN LIQUID LAYERS FLOWING ISOTHERMALLY OVER SOLID SURFACES. Nuclear Engineering Laboratory Memorandum Q 5 , 1961 .

[32]  Jesse Phillips,et al.  Application of the MELCOR code to design basis PWR large dry containment analysis. , 2009 .

[33]  P. Stephan,et al.  Measurement of water falling film thickness to flat plate using confocal chromatic sensoring technique , 2009 .

[34]  S. Bankoff.  Minimum thickness of a draining liquid film , 1971 .

[35]  Thomas J. Hanratty,et al.  Interfacial drag and film height for vertical annular flow , 1985 .