Scaling, Uncertainty, and 3D Coupled Code Calculations in Nuclear Technology
View this Special IssueReview Article  Open Access
Hideaki Ikeda, Takanori Fukahori, Yuichiro Kubo, Hideo Soneda, Shinya Mizokami, "BWR Stability Issues in Japan", Science and Technology of Nuclear Installations, vol. 2008, Article ID 358616, 11 pages, 2008. https://doi.org/10.1155/2008/358616
BWR Stability Issues in Japan
Abstract
The present paper reviews activates relevant to the boiling water reactor (BWR) stability phenomenon, which has a coupled neutronic and thermalhydraulic nature, from the viewpoint of model and code developments and their applications to the BWR stability solution methodology in Japan.
1. introduction
The core power oscillation phenomenon inherently exists in BWR cores [1], as generally called by the BWR stability or instability. The BWR instability is possible even at the normal plant operation conditions, and significant core power oscillations may threaten core fuel integrity due to the fuel cladding dryout occurrence and/or due to the strong PCMI (pelletcladding mechanical interaction). Therefore, an accurate prediction for the onset of BWR instability is indispensable for the safety of BWR core design and operation. Hence, numerous efforts have been paid to understand the complicated BWR instability mechanism and to develop the advanced analysis models.
The stability problem has become an important concern on safety of BWR operations, in particular, after the instability incident at LaSalle2. It should be emphasized that the applied analysis code predicted a stable core condition while instability actually occurred. Therefore, GE and US BWROG (BWR Owners’ Group) have improved the stability analysis models which can be adequately applicable to the actual core design and operation, and have developed the longterm stability solution methodologies with several modifications in the plant installation.
Also in Japan, similar activities have been proceeded by the BWR plant/fuel vendors and utilities to exclude any instability concern. Main goals in Japanese activities are as follows: (1) to analytically investigate the complicated BWR instability mechanism, the power oscillation onset/growth, and formation of the limit cycle oscillation, by using the threedimensional timedomain code; (2) to empirically define the stability performance of the employed fuel design, and to assess the accuracy of calculation results by stability analysis codes using the experimental data, and (3) to establish the stability solution methodology, in which the selected control rod insertion (SRI) system is installed to automatically exclude the operated core from possibly unstable core condition.
The present paper describes the BWR stability issues in Japan. Researches related to the phenomena identification, models, and codes applicable to the design analysis and stability solution methodologies are described. Authors suppose that understanding the basis of the BWR stability issues can be useful for future improvements in the BWR stability solution methodology based on the advanced analysis models and codes. In the last section of the present paper, an outline of the ongoing research on the advanced BWR stability solution methodology is to be introduced, which employs the bestestimate analysis code and the statistical approach in the safety evaluation methodology.
2. BWR Instabilities
The BWR instability can be subcategorized into the three phenomena: (1) channel instability (density wave oscillation); (2) core instability (global core power oscillation); and (3) regional instability (powers in two halves of a core oscillate with an outofphase mode).
2.1. Channel Instability
The channel instability is equivalent to the coolant density wave oscillation in a boiling channel, where the channel pressure drop is kept constant by any constraint [2, 3]. As shown in Figure 1, the coolant void sweeps in the boiling region, which significantly affects the 2phase pressure drop, consequently leads to the coolant mass flow oscillation at the channel inlet. Hence, the channel instability can be invoked in a channel, where the 2phase pressure drop is relatively larger than the singlephase pressure drop, for such conditions as (1) higher channel power and lower flow rate, (2) lower inlet coolant subcooling, (3) downskewed axial power shape, (4) numbers of fuel rods and of fuel spacers which tend to generate the larger pressure drop in the 2phase boiling region. In general, however, excitation of the channel instability can be suppressed by many other stable channels via the neutronic coupling effect among fuel bundles in an actual core.
2.2. Core Instability
The coupled neutronic and thermalhydraulic power oscillation can be categorized into the global instability and into the regional instability. In the first mode, the global core power oscillates inphase, while in the regional oscillating mode, the power in a half core oscillates in an outofphase mode with respect to the other half. The core power oscillation is mainly driven by the negative coolant void feedback with the finite time del ay due to the fuel heat conduction [2]. This power oscillation can be actually excited by synchronizing with the mentioned density wave oscillation, as schematically described in Figure 2. a range from 0.3 to 0.6?Hz [4, 5], which are correlated with the wave propagation velocity through the core fuel channel.
The core power oscillation becomes unstable under the lower flow and higher power core operation condition, corresponding to the density wave oscillation behavior. Large negative void feedback and faster fuel heat conduction make the core state unstable. In addition, the past investigation using frequencydomain stability analysis codes revealed interesting sensitivity with respect to the core power distribution, as shown in Figure 3 [6]. As for the radial power shape, fuel bundles with high power peaking factors tend to reduce the channel stability in the entire core, resulting in the core instability. The sensitivity regarding the axial power shape has more complicated nature as described below. The downskewed shape leads to the longer boiling length, which makes the frequency of the density wave oscillation greater than the time constant in the fuel heat conduction. This mismatch tends to result in the stable core power oscillation. On the other hand, the flat and/or the middleskewed shapes make the greater influence of neutronics in the high void region of the core, inducing the core instability due to increase in the negative void feedback.
(a) Sensitivity of radial core power shape
(b) Sensitivity of axial core power shape
2.3. Regional Instability
The basic phenomenon dominating the regional instability is similar to that for the core instability, and the coupled neutronic and thermalhydraulic oscillation can be individually excited in two halves of a core with an outofphase mode. Previous researchers proposed that the regional instability is equivalent to the oscillation of the higher harmonics (1st azimuthal mode) of the neutron flux distribution, while the core instability is to the oscillation of the fundamental mode (see Figure 4) [7]. Hashimoto derived the socalled ‘modal point neutron kinetics equations in order to analytically represent the phenomenon, in stead of the ordinary point kinetics equations [8]: where is the order of the higher harmonic mode (); , , and are the coreaveraged neutron flux, del ayed neutron precursor, and delayed neutron fraction, respectively. The other variables and notations are defined in the original paper [8]. Physically, represents the subcriticality of the th harmonic mode, which is mathematically corresponding to the eigenvalue separation, and is a negative value in the above definition. Hashimoto [8] and Takeuchi et al. [9] pointed out that a smaller absolute value of the subcriticality makes the feedback gain of the regional oscillation larger, which is correlated to the first term of the righthand side of (1), inducing the regional instability.
(a) Fundamental mode
(b) Higher harmonics (1st Azimuthal) mode
As mentioned above, powers in two halves of a core oscillate with an outofphase mode, therefore, significant oscillations cannot be observed in the coreaveraged power and inlet coolant flow responses. This results in that the hydraulic flow response via the recirculation loop is less sensitive to the regional stability.
3.BWR Stability Analysis Codes, Verifications, and Applications
Several stability analysis codes have been developed so as to investigate the BWR instability phenomena in detail, and to apply on the BWR core design in Japan. The analysis codes can be mainly classified into the two categories, the frequencydomain code and the timedomain code. Features of the frequencydomain stability analysis code and of the timedomain code are summarized in Table 1, respectively.

3.1. ReducedOrder FrequencyDomain Codes
In general, the frequencydomain code employs the reducedorder model like the point neutron kinetics, to mathematically simplify the phenomenological representation, and to attain the faster computation time. In addition, the decay ratio, representing the stability degree of an oscillation, is determinable by the unique methodology based on the system transfer functions. These features are favorable in the design analysis. All the equations representing the physical phenomena are linearized for small perturbations to yield the system transfer functions via the Laplace transformation, which characterize the channel, core, and regional stabilities. The primary physical equations employed in a representative frequencydomain code are the followings: (1)mass, energy, and momentum equations for 2phase mixture boiling coolant flow;(2)radial onedimensional fuel heat conduction equations; and(3)point neutron kinetics equations. The thermalhydraulic behavior in a core is modeled with the parallel channel geometry, and the fuel heat conduction is accounted in each hydraulic calculation node. As for the regional stability analysis, the point neutron kinetics equation is replaced by the modal point kinetics equation as mentioned in the previous section.
Figure 5 shows a sample of verification result for the frequencydomain code, which is currently applied on the BWR core design analysis. The code is able to derive good correlations over the wide stability range for the core stability analysis as well as for the regional stability analysis, while the code models are conservative a priori.
(a)
(b)
3.2. ThreeDimensional TimeDomain Codes
As described above, the frequencydomain codes generally employ the simpler fundamental equation set in order to avoid mathematical difficulties in derivation of the system transfer functions representing the coupled neutronic and thermalhydraulic phenomena in a BWR. The timedomain code, on the other hand, adopts the more sophisticated physical models, like the spatial neutron kinetics model. In fact, their implementation on a code is simple and straightforward, while it consumes larger computational time than the mentioned reducedorder model. Since 1990s, however, the significant advance in computation technologies has facilitated development of timedomain codes that employ the complicated threedimensional and multigroup neutron diffusion kinetics model [1, 10–14]. Several timedomain codes developed by Japanese organizations are listed in Table 2. The most significant advantage in these codes is that the detailed spatial kinetics behavior in a core can be explicitly simulated, namely, both the core stability and the regional stability can be evaluated using a single threedimensional timedomain code without any modification. However, users have to pay attention to the applied numerical time step size, which is sensitive to the simulated oscillation and decay ratio [15, 16].

Furthermore, a simulator has been implemented on the recent timedomain codes in order to accomplish the more realistic dynamic simulation reflecting the actual core state including the fuel history data thus being seamlessly consistent to the static core design [17, 18]. Figure 6 shows a sample of verification for the threedimensional timedomain stability analysis code, SIMULATEKinetics, using the Ringhals1 stability test data [19, 20]. It can be confirmed that the code is basically targeting on the bestestimate stability analysis on the contrary to conservative approach applied in the frequencydomain code. The Ringhals1 cycle14 PT9 stability test, where a regional instability was observed, was accurately simulated as shown in Figure 7 [19]. In addition, the results of numerical simulation demonstrated that the observed regional instability is equivalent to an oscillation of the higher harmonics mode (1st azimuthal, defined by (3), and that modal reactivities ( and defined by (4) are dominant in the regional event as shown in Figure 8.
(a)
(b)
(a) Simulated modal responses
(b) Simulated modal reactivity responses
A feature of the threedimensional timedomain code is that it is applicable to the analytical investigation of the limit cycle oscillation which is driven by the complicated nonlinear effects [21–23]. Figure 9 schematically describes the del iberated mechanism in the formation of limit cycle oscillation. The nonlinear behavior in the 2phase boiling flow tend to increase the coreaveraged void fraction and the negative void feedback, which suppresses the growing neutron flux oscillation due to nonlinearity in the neutron kinetics, resulting in the limit cycles. Any power shift observed in the measured core power responses and/or in the numerically simulated power responses (see Figure 10) is due to the above nonlinearities.
As for another scientific interest on the regional limit cycle oscillation, the bifurcation behavior observed via the spectrum analysis of the measured core power responses [24], Farawila theoretically proposed that the nonlinear interaction in the modal reactivities defined by (4) plays an important role in this phenomenon [25]. In addition, Ikeda et al. have numerically demonstrated that the nonlinearity excites the different higher harmonics of the coreaveraged and regional power responses, respectively, as shown in Figure 11 [23], which was obtained by applying a spectrum analysis to the simulated fundamental and higher modal responses (refer to Figure 8).
(a) Harmonics excitation in measured modal responses [24]
(b) Harmonics excitations in simulated modal responses
4. Current BWR Stability Solution Methodology
Since the instability incident at LaSalle2 [26], GE, and USBWROG has developed several longterm stability solution methodologies [27, 28]. Also in Japan, a similar stability solution methodology was established, where the adequate stability margin must be ensured in the core design process, and the selected control rod insertion (SRI) system is equipped to exclude the BWR core from the unstable operation region (stability exclusion region) as shown in Figure 12 [29].The SRI system is activated to suppress the core power when the core coolant recirculation pumps are tripped and the core goes into the preliminary determined stability exclusion region. The stability exclusion is to be determined by using stability design codes certified via the regulatory assessments, with the conservative stability criteria (decay ratio is less than 0.8). Consequently, this methodology is targeting on that the BWR instabilities are not possible in the operated core in Japan.
5. Research on Advanced BWR Stability Solution Methodology
The current stability solution methodology is effectively contributing to safety of BWR plant operations in Japan. However, considering the recent occurrences of BWR instabilities [30, 31], authors suppose that any improvement may be indispensable for the future stability solution methodology, which is able to correspond to the recent modifications in the existent BWR plants as the extended core thermal poweruprate [32] with the advanced fuel designs [33–35]. An approach to resolve this concern is that sufficient stability margin is to be introduced, namely, the plant operable region is limited by the wider stability exclusion region, which can be determined by using the current conservative stability analysis code, as shown in Figure 13. This approach, however, possibly leads to the economical loss by consuming longer time for the plant startup operation. Reactors with the larger stability exclusion region are generally allowed to adopt few continuous withdrawals of control rods under the lower power condition. This is because the continuouswithdrawal operation induces significant increase in core power at the fixed core flow condition, possibly removing the core into the prohibited stability exclusion region. Therefore, a lot of control rod operations, which must be conducted slowly and intermittently to maintain the fuel mechanical integrity, are required under the higher power condition to attain the target control rod pattern at the rated power operation. Consequently, the overall plant startup time tends to become longer in the BWR plant with the larger stability exclusion region.
In order to reasonably enhance the operable region even under the poweruprated core, a joint research group organized by several Japanese industrial and academic organizations has started a development of an advanced stability solution methodology based on the bestestimate code system [36]. Basis of the present research is to apply the original regulatory criterion with respect to the BWR instability [37], that is, “exceeding specified acceptable fuel design limits (SAFDLs) are not possible”, not prevention from the instability occurrence. From the viewpoint of the applicable SAFDLs on the BWR instability, the PCMI and the material fatigue via the power oscillation possibly make no significant affect on the fuel integrity, because temperature responses of the fuel pellet and cladding are negligibly small as shown in Figure 14. Therefore, occurrence of the core coolant boiling transition (BT) can be a primary cause for the fuel failure under the BWR instabilities. So as to accurately and mechanistically predict the BT onset even under the BWR instabilities, the research group is applying an advanced code system based on the bestestimate plant simulator, TRACBF1/ENTRÉE [13], and the 2fluid/3field subchannel code, NASCA [38]. As schematically described in Figure 15, TRACBF1/ENTRÉE provides the pinbypinbased power responses in each fuel bundle; the subchannel thermalhydraulic behavior and BT onset on the local rods are evaluated by NASCA with the boundary conditions supplied by the TRACBF1/ENTRÉE.
The research group is also investigating the possibilities to introduce the statistical safety evaluation methodology [39] so as to establish the reasonable conservatism in the stability exclusion region determined by using the above bestestimate code system. The research, in particular, currently pays a lot of efforts to establish the phenomena identification ranking table (PIRT) applicable to BWR instabilities including the subchannel thermal hydraulics, based on the existent stability PIRTs [40–42]. This is the basis of the uncertainty evaluation for the bestestimate BWR stability analysis.
6. Conclusions
Many efforts have been paid to research on BWR stability issues in Japan, as introduced in the present paper. The industrial organizations have developed and improved the BWR stability analysis using computational tools specific for the reducedorder frequencydomain and threedimensional timedomain codes. The first category is currently applied to the BWR stability design analysis, while the latter one has been exploited to understand the complicated phenomena related to BWR stability. The current stability solution methodology based on the SRI system with the stability exclusion region is successfully preventing the occurrence of BWR instabilities in Japan. However, authors suppose that the future application of the extended core power uprate requires further improvements to the current solution methodology in order to reasonably minimize the stability exclusion region. A Japanese research group is currently proposing to apply the bestestimate analysis code with the statistical safety evaluation methodology. This will allow better evaluation of the stability exclusion region, and will be consequently applied to the BWR plants with the extended core power uprate.
References
 F. D'Auria, “State of the Art Report on Boiling Water Reactor Stability,” 1997, OCDE/GD(97)13, AEN/NEA. View at: Google Scholar
 R. T. Lahey Jr. and F. J. Moody, The ThermalHydraulics of a Boiling Water Nuclear Reactor, American Nuclear Society, Chicago, Ill, USA, 2nd edition, 1993. View at: Google Scholar
 J. MarchLeuba, “Density wave instabilities in boiling water reactors,” Tech. Rep. NUREG/CR 6003 ORNL/TM12130, Oak Ridge National Laboratory, 1992. View at: Google Scholar
 L. A. Carmichael, “Transient and stability tests at peach bottom atomic power station unit 2 at end of cycle 2,” Electr Power Res Inst Rep EPRI NP 564, 1978. View at: Google Scholar
 B. G. Bergdahl and R. Oguma, “BWR stability investigation in ringhals 1 measurement data from October 26, 1989,” in Proceedings of International Workshop on Boiling Water Reactor Stability, vol. 142, New York, NY, USA, 1990. View at: Google Scholar
 “Standard for assessment of coupled neutronic and thermalhydraulic stability in BWR power plants,” to appear in Atomic Energy Society of Japan. View at: Google Scholar
 J. MarchLeuba and E. D. Blakeman, “Mechanism for outofphase power instabilities in boiling water reactors,” Nuclear Science and Engineering, vol. 107, no. 2, pp. 173–179, 1991. View at: Google Scholar
 K. Hashimoto, “Linear modal analysis of outofphase instability in boiling water reactor cores,” Annals of Nuclear Energy, vol. 20, no. 12, pp. 789–797, 1993. View at: Publisher Site  Google Scholar
 Y. Takeuchi, Y. Takigawa, and H. Uematsu, “Study on boiling water reactor regional stability from the viewpoint of higher harmonics,” Nuclear Technology, vol. 106, no. 3, pp. 300–314, 1994. View at: Google Scholar
 O. Yokomizo, “Analysis of BWR core stability with spacedependent neutron feedback,” Transactional Analysis, vol. 44, p. 560, 1983. View at: Google Scholar
 Y. Takigawa, Y. Takeuchi, S. Ebata, C. K. Chan, C. Tricoli, and S. Tsunoyama, “Caorso limit cycle oscillation analysis with threedimensional transiet code TOSDYN2,” Nuclear Technology, vol. 79, no. 2, pp. 210–227, 1987. View at: Google Scholar
 Y. Takeuchi, Y. Takigawa, H. Uematsu, S. Ebata, J. C. Shaug, and B. S. Shiralkar, “TRACG transient analysis code—threedimensional kinetics model implementation and applicability for spacedependent analysis,” Nuclear Technology, vol. 105, no. 2, pp. 162–182, 1994. View at: Google Scholar
 A. Hotta, M. Zhang, and H. Ninokata, “BWR regional instability analysis by TRAC/BF1ENTRÉEII: application to ringhals unit1 stability test,” Nuclear Technology, vol. 135, no. 1, pp. 17–38, 2001. View at: Google Scholar
 “Present status and subjects on the evaluation of neutronic and thermalhydraulic stability in BWR,” 2001, Atomic Energy Society of Japan. View at: Google Scholar
 J. G. M. Andersen, “TRACG time domain analysis of thermal hydraulic stability sensitivity to numerical method and comparison to data,” 1989, Stability Symposium. View at: Google Scholar
 A. Hotta, M. Honma, H. Ninokata, and Y. Matsui, “BWR regional instability analysis by TRAC/BF1ENTRÉEI: application to densitywave oscillation,” Nuclear Technology, vol. 135, no. 1, pp. 1–16, 2001. View at: Google Scholar
 G. M. Grandi and K. S. Smith, “BWR atability analysis with SIMULATE3K,” in Proceedings of Physics of Nuclear Reactors, Seoul, Korea, October 2002. View at: Google Scholar
 M. Tamitani, T. Iwamoto, and B. R. Moore, “Development of kinetics model for BWR core simulator AETNA,” Journal of Nuclear Science and Technology, vol. 40, no. 4, pp. 201–212, 2003. View at: Publisher Site  Google Scholar
 H. Ikeda and A. Hotta, “Development of a timedomain threedimensional core analysis code and its application to spacedependent BWR stability analyses,” in Proceedings of Joint International Conference on Mathematical Methods and Supercomputing for Nuclear Applications, vol. 1624, Saratoga Springs, NY, USA, October 1997. View at: Google Scholar
 T. Lefvert, “OECD/NEA Nuclear Science Committee BWR Stability Benchmark Preliminary Specifications,” 1993, NEA/NSC/DOC(93)13. View at: Google Scholar
 J. MarchLeuba, D. G. Cacuci, and R. B. Perez, “Nonlinear dynamics and stability of boiling water reactors: part 1—qualitative analysis,” Nuclear Science and Engineering, vol. 93, no. 2, pp. 111–123, 1986. View at: Google Scholar
 J. MarchLeuba, D. G. Cacuci, and R. B. Perez, “Nonlinear dynamics and stability of boiling water reactors: part 2—quantitative analysis,” Nuclear Science and Engineering, vol. 93, no. 2, pp. 124–136, 1986. View at: Google Scholar
 H. Ikeda, T. Ama, K. Hashimoto, and T. Takeda, “Nonlinear behavior under regional neutron flux oscillations in BWR cores,” Journal of Nuclear Science and Technology, vol. 38, no. 5, pp. 312–323, 2001. View at: Google Scholar
 G. Th. Analytis and J. K.H. Karlsson, “On the existence of multiple harmonics in neutronic spectra of regional oscillations in the ringhals1 boiling water reactor,” Nuclear Science and Engineering, vol. 131, no. 2, pp. 286–292, 1999. View at: Google Scholar
 Y. M. Farawila, “Application of modal neutron kinetics to boiling water reactor oscillation problems,” Nuclear Science and Engineering, vol. 129, no. 3, pp. 261–272, 1998. View at: Google Scholar
 USNRC, “AEOD concerns regarding the power oscillation event at LaSalle 2 (BWR5),” AEOD Special Report S803, 1988. View at: Google Scholar
 General Electric Company, “BWR Owners' Group Long Term Stability Solutions Licensing Methodology,” 1991, NEDO31960. View at: Google Scholar
 General Electric Company, “BWR Owners' Group Long Term Stability Solutions Licensing Methodology,” 1995, NEDO31960A. View at: Google Scholar
 “Present status and subjects on the evaluation of neutronic and thermalhydraulic stability in BWR,” 2001, Atomic Energy Society of Japan. View at: Google Scholar
 M. Kruners, “Analysis of Instability Event in Oskarshamn3,” February 1998, SIMULATE3K', SKI/RAPORT 98:43. View at: Google Scholar
 M. May, “Nine Mile Point 2 Instability Event Briefing,” August 2003, BWR Owners' Group Mtg., Presentation to USNRC. View at: Google Scholar
 General Electric Company, “Generic Guidelines and Evaluations for GE BWR Thermal Power Optimization,” 2002, NEDO32938 Rev.1. View at: Google Scholar
 G. A. Watford, “General Electric Fuel Bundle Designs,” 2000, NEDO31152 Rev.7. View at: Google Scholar
 N. L. Garner and P. Urban, “ATRIUM 10: ten years of operational experience,” Nuclear Engineering International, vol. 47, no. 578, pp. 30–34, 2002. View at: Google Scholar
 S. Helmersson, “Advances in BWR fuel development,” Nuclear Europe Worldscan, vol. 21, 2001. View at: Google Scholar
 A. Hotta, “Research on hyper ABWR plants of extended core power density,” 2006, Innovation and Viable Nuclear Energy Technology (IVNET) Development Project, Ministry of Economy, Trade and Industory. View at: Google Scholar
 “Regulatory Guide for Safety Design in Light Water Reactor Power Facilities,” 1990, Nuclear Safety Committee of Japan. View at: Google Scholar
 H. Shirai, H. Ninokata, A. Hotta, and T. Hara, “Analytical study on detailed void distributions inside BWR fuel bundle under turbine trip event considering timedependent pin power distributions,” Journal of Nuclear Science and Technology, vol. 39, no. 11, pp. 1152–1161, 2002. View at: Publisher Site  Google Scholar
 B. Boyack, “Quantifying Reactor Safety Margins—Application of Code Scaling, Applicability and Uncertainty Methodology to a LargeBlake LossofCoolant Accident,” 1989, NUREG/CR5249. View at: Google Scholar
 SCIENTECH, “BWR PIRT and Assessment Matrices for BWR LOCA and nonLOCA Events,” 1999, SCIENRC39399. View at: Google Scholar
 USNRC, “PIRTs for Power Oscillations without Scram in BWRs containing High Fuel Burnup,” 2001, NUREG/CR6743. View at: Google Scholar
 General Electric Company, “Licensing Application Framework for BWR Stability Limit Analysis,” 2003, NEDO33095. View at: Google Scholar
Copyright
Copyright © 2008 Hideaki Ikeda et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.