presentation_ups by xiaoyounan

VIEWS: 1 PAGES: 37

									Thermal Hydraulic Studies for
     PFBR using PHOENICS


         U. PARTHA SARATHY

  Indira Gandhi Centre for Atomic Research
                 Kalpakkam
              May 3-5 th 2004
         PROTOTYPE FAST BREEDER REACTOR
                     (PFBR)
       Power - 500 MWe, 1250 MWth
       Fuel – Mixture of UO2 (79 %) and PuO2 (21 %)
       Coolant – Sodium (liquid metal) in Pry and Secy Circuits
                         – Water in Tertiary Circuit
       High Temperatures
       High Velocities
       Problems –
          High temperatures leading to creep, fatigue damage
          Flow induced vibrations
          Thermal striping
            Gas entrainment

2         U. Partha Sarathy, IGCAR   05/05/2004
                                           PFBR Primary Circuit

    Inner
    Vessel                              PUMP
                                                                          IHX

                                                                  Hot Pool




                                                       Nuclea
                                                       r heat
                                                                   CORE



                              Cold Pool              Grid Plate




3            U. Partha Sarathy, IGCAR   05/05/2004
                         Schematic PFBR Flow Sheet




    •Primary Circuit •Secondary Circuit      •Steam/Water circuit
4    U. Partha Sarathy, IGCAR   05/05/2004
HYDRAULIC ANALYSIS OF
    GRID PLATE- e Page
                 HYDRAULIC ANALYSIS OF GRID PLATE


                                            • Consists of 1758 sleeves
                                            • Receives flow from four pipes
                                            • Distributes flow to various
                                                subassemblies


                                            Objectives

                                               Flow and pressure distribution

                                               Pressure drop in GP

                                               Velocity over sleeves


6   U. Partha Sarathy, IGCAR   05/05/2004
         HYDRAULIC ANALYSIS OF GRID PLATE
                                                 Modelling
                                                    2-D model in cylindrical co-ordinates
                                                     (r- θ)
                    SECTION A - A                   Sleeves modeled through porosity in
                                       
                                                     radial and circumferential directions
                                                     (Porous body formulation)
                                                    Inlet as Velocity BC
                                                    Outlets as mass sinks
    
                                                    Pressure drop due to sleeves
                                                     modeled through Zukauskas
                                                     correlation
                                                    Addition of resistance terms in the
                                                     momentum equation using ‘ground’
                        PLAN



                                                     subroutine.
        Schematic of Grid Plate
                                                    K-E Turbulence model
7        U. Partha Sarathy, IGCAR   05/05/2004
                         Results of Grid Plate Analysis
                                               Results
                                                Predicted   ΔP is 4.6 m
                                                 of sodium
                                                Similar to that
                                                 extrapolated from 1:3
                                                 scale air
                                                 experiments.
                                                Pressure  contours
                                                 are concentric –
                                                 uniform flow through
                                                 fuel SA
    Flow Distribution in Grid Plate
                                                Maximum     cross flow
                                                 velocity is 8.5 m/s

8      U. Partha Sarathy, IGCAR   05/05/2004
           Thermal Analysis of
Hot and Cold Pools- Title Page
     Thermal Analysis of Hot and Cold Pools


                                                    Objectives
                                                       Inner Vessel
                                                        temperature distribution
                                                       Stratification In sodium
                                                        pools
                                                       Hot pool free surface
                                                        velocity & temperature
                                             CORE




10   U. Partha Sarathy, IGCAR   05/05/2004
          CFD Model and Boundary Conditions
                                             Modelling
                                                2-D model in cylindrical co-
                                                 ordinates (r-z)
                                                Core is modeled as a block
                                                Porous body approximation for
                                                 immersed components – IHX,
                                                 Pump
                                                Mass sink at IHX & PUMP
                                                 inlets
                                                Velocity BC at IHX and Core
                                                 outlets
                                                Conjugate thermal hydraulic
                                                 analysis of hot & cold pools
                                                 including IV
                                                K-E Turbulence model

11   U. Partha Sarathy, IGCAR   05/05/2004
      Flow Distribution in Hot and Cold Pools




        Good mixing in hot and cold pools

12   U. Partha Sarathy, IGCAR   05/05/2004
                                                   Results
                                                    Tmax   in IV is 534 OC
                                                    ΔT   across thickness is 64 K
                                                    Max hot pool free surface
     Temperature Distribution in Inner Vessel
                                                     temperature is 572 OC




         Hot Pool Free Surface Temperature
                    Distribution
13         U. Partha Sarathy, IGCAR   05/05/2004
Flow Distribution in SG Inlet Plenum-
                         Title Page
                                                                          TOP PLATE
                                                      = 180°
                                                                            = 0°
                                                                    R




                                                                                      460
                                                                          TUBE BUNDLE
                                                                          SHROUD


                                                                          OUTER SHELL




                                                                                      1200
                                                                             INLET


                                                                 OD 356     NOZZLE
                                                                 ID 520




                                                                                      800
                                                                  EXIT



      Schematic of PFBR SG                    3/5 scale model of SG Inlet Plenum

     Objective:To identify flow distribution devices and reduce
       maximum radial velocity over tubes from FIV considerations.
15    U. Partha Sarathy, IGCAR   05/05/2004
                                          TOP PLATE
                  = 180°

                                R
                                             = 0°
                                                             Modelling




                                                      460
                                                              3/5   scale model
                                          TUBE BUNDLE
                                          SHROUD
                                                              3-D cylindrical
                                          OUTER SHELL

                                                               coordinates




                                                      1200
                                                              180 O    symmetric model
                                             INLET


                             OD 356
                                                              K-E     turbulence model
                                             NOZZLE
                             ID 520
                                                      800
                                                              Inlet   as velocity BC


                              EXIT


     3/5 scale model of SG Inlet Plenum




16       U. Partha Sarathy, IGCAR     05/05/2004
                                    INNER
                                    SHELL
                                                                                                                             =0


                                                Flow distribution in Inlet window region at
                                                            1430 mm from inlet
                                                                                       Theta = 15             Theta = 45            Theta = 60
                                                                                       Theta = 75             Theta = 90

                                                                         1700

                                  INLET
                                                                         1600




                                                 Height from Inlet, mm
                                                                         1500



                                                                         1400



                      7.88 E + 00                                        1300
             Min. : 1.87 E - 01
             Max. : 5.19 E + 00                                          1200
                                                                                -1   -0.5     0     0.5        1       1.5   2        2.5        3
      Flow distribution in SG Inlet                                                                       Velocity, m/s

     Plenum – Basic Configuration                Fig. 5 Radial Velocity Profile Along the Window with Basic Configuration
                                                                           Radial Velocity Profile along the Window
                                                                                   with Basic Configuration
17     U. Partha Sarathy, IGCAR    05/05/2004
                             4.0

                             3.5

                             3.0
     Axial Velocity, m/s

                                                           Basic Configuration
                             2.5

                             2.0

                             1.5

                             1.0

                             0.5

                             0.0
                                      0                30          60         90       120   150   180
                                                                        Angle, degrees

                                   Axial Velocity in the Annulus mm – Basic Configuration
                               Axial Velocity in the Annulus at 575at a height of 575 mm from
                                                  inlet with various Porous Plates



18                         U. Partha Sarathy, IGCAR   05/05/2004
           TOP PLATE
                                                                                                                                                                          R = 260
               = 180°                                      = 0°                                                                                                        R = 245
                                        R                                                                                                                                          R = 219
           60 % porosity




                                                      230
                                                      230
           20 % porosity




                                                                                                                           Po
                                                                                                                                         R = 193




                                                                                                                            d4 ity =
                                                                                                                             ros
                                                      460




                                                                                                                                                                                                          5
                                                                                                                               (Ø 60%
           TUBE BUNDLE




                                                                                                                                 25
                                                                                                                                         R = 178




                                                                                                                                    )
           SHROUD                                                                                           d5 (Ø24)




                                                                                                                                                   65%
                                                                                                           Porosity = 58%




                                                                                                                                                                                                    55
          Porous plate                                                                                                                                         d2 (Ø25)




                                                                                                                                            sity =
                                                                                                                                        Poro 25)
                                                                        1200




                                                                                                                                                                                         45
                                                                                                                                                         Porosity = 63%




                                                                                                                                        d3 (Ø




                                                                                                                                                                                                         58
                                                                                                                                                                     5%
                                                                                                                                                                  y=9




                                                                                                                                                                                              48
                                                                                                      d6 (Ø22.5)
                                                                                                                                                               sit




                                                                                                                                                                           15°
           OUTER SHELL                                                                                                                                    Poro




                                                                                                                              °
                                                                                                    Porosity = 55%




                                                                                                                          (Typ)
                                                                                                                          30
                                                      450




                                                                                                                                                                                                              2
                                                                                                                                                                                                4
                                                                INLET             (  = 180°)   Radial pitch = 26 mm                                                         (  = 0°)        34         34
                                                                                                Plate thickness = 10 mm
                                                                                                Material = Carbon steel
                                  OD 356
                                                          NOZZLE
                                  ID 520
                                                                               Porous plate used as a Flow distribution devices

                                                                        800




                                      EXIT


     3/5 scale model of SG Inlet plenum with
            Flow distribution devices


             Porous body formulation                                                        for
         porous plate and porous shell

19         U. Partha Sarathy, IGCAR          05/05/2004
                                                Basic Configuration                      74.5 % to 24.4 %
                                                75 % to 55 %                             80 % to 55 %
                                              4.5

                                               4

                                              3.5

                        Axial Velocity, m/s    3

                                              2.5

                                               2

                                              1.5

                                                1

                                              0.5

                                               0
                                                    0           30    60    90     120      150   180
                                                                     Angle, degrees
                    Fig. 5 Axial Velocity in the Annulus at 575 m m
                        from inlet w ith different 575 mm from Inlet with
                     Axial Velocity in the Annulus atPorous Plates
                                                                 Different Porous Plates

20   U. Partha Sarathy, IGCAR                      05/05/2004
                                                    Theta = 15 deg             Theta = 75 deg          Theta = 135 deg




                              1600
     Height from Inlet, mm




                              1500


                              1400


                              1300


                              1200
                                       -2.5         -2      -1.5       -1   -0.5     0     0.5     1   1.5      2        2.5   3
                                                                                   Velocity, m/s
                             Fig. 9 Radial Velocity Profile Along the Window with Porous Plate (80 % to 55 %)




21                             U. Partha Sarathy, IGCAR   05/05/2004
                                                            TOP PLATE

                                                                = 180°                      = 0°
                                                                                  R

                                        Inner shell         60 % porosity




                                                                                       230
                                                                                       230
                                                            20 % porosity




                                                                                       460
                                                            TUBE BUNDLE
                                                            SHROUD




                                                                                                         1200
                                                            OUTER SHELL




                                                                                       450
                                                                                                 INLET


                                                                            OD 356
                                                                                         NOZZLE
                                                                            ID 520




                                                                                                         800
                                        Porus plate


                                                                                EXIT


                                        Inlet         3/5 scale model of SG Inlet plenum with Flow
                                                                  distribution devices




                    6.42 E + 00
             Min. : 4.61 E - 02
             Max. : 4.44 E + 00
                                                                                                                (= 0)
     Flow distribution in SG Inlet plenum
     with Flow distribution devices       Flow distribution in Inlet window
                                            region at 1430 mm from inlet
22      U. Partha Sarathy, IGCAR   05/05/2004
                               Basic configuration          With Porous Plate      With Porous Plate and Shell

                               1700



       Height from Inlet, mm   1500



                               1300



                               1100



                                900



                                700
                                      -1    -0.5        0       0.5      1      1.5      2      2.5      3
                                                                   Velocity, m/s

                                       Velocity Profile Along the Window at 135 deg


23   U. Partha Sarathy, IGCAR              05/05/2004
     RESULTS
        Combination of graded porous plate and porous shell render as
         uniform flow both axially and circumferentially.
        The distributions of porosity in the plate and shell have been
         identified.
        Maximum radial velocity is 0.75 m/s (average is 0.45 m/s)
         whereas the same is 3 m/s in basic configuration




24   U. Partha Sarathy, IGCAR   05/05/2004
Inter-Wrapper flow Studies-Title
                          Page
                      Inter-Wrapper flow Studies -
                              Steady State
     Inter Wrapper flow                           Objectives
                                                  •Effect of IWF on SA
                                                  clad hotspot
                                                  •Flow distribution in
                                                  IWS
                                                  •To develop a model
                                                  for studying various
                                                  design basis events
                                                  which will give detailed
                                                  temperature
         Sub-Assembly                             distribution in hot and
         Steel hexagonal                          cold pools
         Wrapper

26        U. Partha Sarathy, IGCAR   05/05/2004
                     Sodium Flow in Primary Circuit
                                             DHX




                                                   CORE




27   U. Partha Sarathy, IGCAR   05/05/2004
                                                Modeling
                                                 2-D cylindrical
                                                  coordinates (r-z)
                                                 Inlets   as velocity BC
                                                 Outlets   as mass sink
                                                 Porous  body formulation
                                                  for core and other
                                                  immersed structures
                                                 Coupling  with 1-D model
                                                  for neutronics, heat
                                                  transfer calculations in
                                                  core, IHX, DHX etc.
     CFD model for IWS and Hot and
             Cold Pools
28      U. Partha Sarathy, IGCAR   05/05/2004
                                                                                      Sub Assembly Top




                                                                      1290
                                                   SA OUTLET REGION




                                                                                                               SA SODIUM




                                                                                                                              Heat transfer to IWS
                                                  TOP AXIAL BLANKET




                                                                      300




                                                                                                        CLAD
                                                                      1000




                                                                                           FUEL 2
                                                                             FUEL 1
                                                        ACTIVE CORE




                                                                      300
                                               BOTTOM AXIAL BLANKET




                                                                      1010
                                                    SA INLET REGION


                                                                                Sub Assembly Bottom


                                                                                         SSA        -   Storage subassembly
                                                                                         IWS        -   Inter-wrapper space
                                                                                                    -   calculated temperature


                                               Schematic of the SA Computational
     Schematic of Fuel SA                                    Model
29     U. Partha Sarathy, IGCAR   05/05/2004
                                                     CORE FLOW

                                                   SA OUTLET TEMP.
                                                 HEAT TRANSFERRED
                                                       TO IWS
                                                  IHX PRIMARY FLOW
                                                 PRIMARY PUMP FLOW
                                                 IHX and DHX PRIMARY
                            1 D CODE                                        PHOENICS
                                                OUTLET TEMPERATURES
                    PRIMARY HYDRAULICS
                                                                         TWO DIMENSIONAL
                                CORE                                      MODEL OF HOT

                                IHX                                      AND COLD POOLS

                                                                          INCLUDING IWS,
                          DECAY HEAT              IWS TEMPERATURE
                                                                            IV and MV
                      REMOVAL SYSTEM
                                                IHX PRIMARY INLET TEMP

                                               DHX PRIMARY INLET TEMP

                                               PRIMARY PUMP INLET TEMP




                  Exchange of Results between 1-D and 2-D
                  PHOENICS Models for Boundary Conditions


30   U. Partha Sarathy, IGCAR     05/05/2004
                                                    START



                                                    iter = 0



                                             Guess the temperatures
                                               of Pump inlet, IWS
                                               Hot and Cold pools



                                                 1-D code
                                                   ID
                                                  steady state


                                                         BC


                                              Heat transfer to IWS      IWS temperature
                                             Flow and temperatures    Inlet temperatures of
                                               of SA, IHX and DHX     IHX, Pump and DHX

                                                                                   BC

                                                                          iter = iter +1

                                                                                  NO

                                                       2D
                                                     2-D                   Check for
                                                  PHOENICS               convergence


                                                                                  YES


                                                                             STOP




                   Flow Chart for Coupled 1D Code – PHOENICS code
                                     Calculations
31   U. Partha Sarathy, IGCAR   05/05/2004
                                             m/s




                m/s

       Flow Distribution in Hot                    Temperature Contours in
           and Cold pools                          Hot and Cold pools

32   U. Partha Sarathy, IGCAR   05/05/2004
                                             425


                                             415
         555
                                               405




                                                         395 OC




                                                                  m/s




                                                   m/s


              Temperature and Velocity Distribution in Inter-
                            Wrapper Space
33   U. Partha Sarathy, IGCAR   05/05/2004
                                                 Results
                                                 •SSA outlet
                                                 temperature increases
                                                 by about 2 K
                                                 •Total heat transferred
                                                 to IWS is 370 kW
     Temperature Distribution in IV
                                                 •Axial temperature
                                                 gradient of hot/cold
                                                 interface is 150 K/m




      Temperature Distribution in MV
34       U. Partha Sarathy, IGCAR   05/05/2004
            Inter-Wrapper flow Studies -
        Transient Analysis (under progress)
      Station            blackout incident

      All   pumps trip

      Primary               circuit flow coasts down

      Secondary                    circuits not available

      Reactor               trips only at 2.5 s

      Temperature                      inside SA goes up

      Good          amount of heat is taken away by the IWF



35       U. Partha Sarathy, IGCAR   05/05/2004
            Results

Transient Evolution of Temperatures in
         Hot and Cold Pools

								
To top