Pseudocritical Rapid Energy Dissipation Analysis of Base-Load Electrical Demand Reduction on Nuclear Steam Supply System ()
1. Introduction
A considerable reduction of base-load electricity demand generated by nuclear power plant (NPP) for industrial and domestic applications, classified as common cause failure (CCF), may cause the temperature of the reactor coolant to increase and corresponding pressure increase in Pressurizer and Steam Generators. The short fall in terms of electricity demand triggers an integrated plant control system (IPCS) of nuclear instrumentation and control (I & C) systems into action, to compensate for the demand reduction of the base-load, since the amount of heat released from the nuclear chain reaction must be proportional to the amount of heat taking away for generation of electrical energy.
The Pressurizer in NPP maintains the pressure of reactor coolant system (RCS), preserves the threshold parameters through the steady state operations, and regulates the pressure during the transient process of the Reactor. The Pressurizer pressure control system (PPCS) forms part of the IPCS of the nuclear steam supply system (NSSS) of the Pressurized Water Reactor (PWR) [1]. It offers three main functions: 1) to protect the reactor from trip, 2) to protect the reactor from changes in reactivity and 3) to protect the activation of pressure relief valve [2]. The introduction of Pseudocritical rapid energy dissipation (PRED) pressure control, seeks to maintain the safety of the Pressurizer and the RCS.
The Pressurizer operates on three main functions that take place concurrently; 1) the dynamics of temperature variation, 2) the rise and fall of pressure values and 3) the variation of water/steam levels. The rod control system regulates the in-core rod position, the power and the power distribution. The input of the pressurizer pressure is the deviation from a set-point, while the output drives proportional spray controller, to actuate heaters and power relief valve within fixed operating points. The other system controllers provide ex-core integrated plant control by synchronization of fluid flow and control signal activities.
Most of the reactors in operation today are about 40 - 50 years older. The assessment of existing Instrumentation and control (I & C) systems, and its modernization needs in terms of obsolete systems and the impact of failure rate, the NSSS-I & C appeared to be of high priority for modernization [3], the NSSS has reach its operational limits and condition (OLC). The nuclear I & C system serves as nerve center of NPP [4], to coordinate activities of thousands of components and equipment [5], that allow the plant operators to monitor the safe operation of NPP. The advances made by the evolution of digital technology have improved operations of automation and control systems, and yet the nuclear industry remains conservative with the use of analogue systems, hence the need for the proposed modern PRED I & C system.
The level control of pressurizer represent balance of injection of water inflow into the reactor coolant system and water letdown into the chemical and volume control system (CVCS) [6]. The two-phase simulation model of compressibility involving the thermodynamic equilibrium revealed sub-critical evaporation with large variation of liquid compressibility factor and density [7].
A two-phase turbulence model where interfacial turbulence accounts for pseudoturbulence in liquid bubble-induced mixing [8]. The hyperbolic transition model of two-phase fast depressurization is highlighted on six equations of two-phase model, for accurate tabulation of equation of state (EoS) for thermodynamic equilibrium recovery, for depressurization of water-steam system [9]. Which referred to PWR with water as coolant with highly non-linear behavior of two-phase, model provided thermodynamic separation between water and steam independent of each other.
Modern instrumentation and control system was used to examine the effect heat transfer of fluid properties, MATLAB/Simulink library tools was used to: 1) design a model of pseudocritical rapid energy dissipation I & C system, where thermophysical properties of the two-phase fluid characteristics with introduction of three thermodynamic property sensors. Analyze pseudocritical saturation line, the effect of specific internal energy, pressure vector and other parameterization of fluid content of Pressurizer. Evaluated the effect of heat transfer coefficient (Prandtl number) and momentum of the transport of the kinematic viscosity, as pseudocritical rapid energy dissipation from pressurizer.
2. Theories
The pressurizer relief line connect the pressurizer safety and relief valve to the pressurizer relief tank with the input coming from refueling water and storage tank.
The energy conservation in the pressurizer sometimes referred to in the study as pipe or open-channel can be expressed as [10].
(1)
where, M is the fluid mass inside the pressurizer,
and
are the mass flow rate into the pipe through port A and B, while uI is the specific internal energy of the fluid inside the pressurizer, ϕA is the energy flow rate into the pressurizer through port A. While ϕB is the energy flow rate into the pressurizer through port B, and QH is the heat flow rate into the pressurizer through the pressurizer wall represented by port subscript H is the thermal conductive insulator.
The energy conservation in the pressurizer is given by the energy balance equation as expressed in “Equation (2)” [11]:
(2)
where: PA is the pressure at port A, S is the cross-sectional area of the pressurizer, VA is the specific volume of the fluid at port A, Fvisc,A is the viscous friction force is port A, IA is the fluid inertia at port A which is further expressed as in “Equation (3)”,
(3)
where L is the length of the pressurizer in the half pipe adjacent to port B.
Similarly, for port B, the energy conservation in the pressurizer can be expressed as,
(4)
The heat flow rate inside the walls of the pressurizer and the internal fluid volume was modelled as
(5)
where, hcoeff is the average heat transfer coefficient, Ssurf is the surface area, TH is the pressurizer wall temperature and TI is the temperature of the fluid in the pressurizer.
The heat transfer coefficient rely on the phase of the fluid. With subcooled liquid and superheated vapour conditions, [12] the coefficient was modelled as:
(6)
The asterisk represent a value of the phase under consideration (liquid or vapour). While Nu is the average Nusselt number in the pipe, kI is the average thermal conductivity in the pressurizer, and Dh is the hydraulic diameter.
(7)
The subscript M refers to specific value of two-phase mixture and SL is the saturated liquid.
The Nusselt number have effect on convective and conductive heat transfer coefficient, for laminar flow assumed as constant, while the Reynolds number acts on the inertia and viscous forces applied when the value is smaller than the laminar flow [13].
The Gnielinski correlation [14], expressed in “Equation (8)” was used to compute the Nusselt number in the liquid and vapour phases:
(8)
where, f is the friction factor of the pressurizer, Re is the Reynolds number and PrI is the Prandtl number.
The friction factor was evaluated as:
(9)
where, εr is the roughness of the pressurizer inner walls. The Reynolds number was expressed as:
(10)
where the subscript Avg represent an average value between the ports, S is the cross-sectional area, vI is the specific volume and vI is the kinematic viscosity.
For the two-phase mixture, the Nusselt number for turbulent flow correlation expressed by Cavallini and Zecchin [15] was utilized as in “Equation (11)”.
(11)
where subscript SL represent the saturated liquid and SV saturated vapor, xI is the vapour quality and v is the specific volume.
The Reynolds number of the saturated liquid is expressed as
(12)
The conservation of mass balance in the pressurizer was modelled according to Equation (13)
(13)
where, ρ is the fluid density,
is the pressure inside the pressurizer, V is the volume of fluid,
is the mass flow rate into the pipe through port A,
is the mass flow rate into the pipe through port B and εM is the correction term density partial derivation of transition phase boundaries.
The density partial derivative using cubic polynomial function of subcooled blended with superheated vapour in two-phase mixture domain was accounted for according to “Equation (14)”.
(14)
where the εM is the cubic polynomial function introduces the numerical error has always a single inflection point. M is the fluid mass in the pressurizer and is expressed as where vI as the specific volume of the fluid in the pressurizer, τ is the phase change time constant.
(15)
In the laminar regime, the Reynolds number limit value in port A was specified by:
(16)
and the half pipe adjacent to port B by:
(17)
where, fshape is the pipe shape factor, Leff is the effective pipe length and Dh is the hydraulic diameter of the pipein the turbulent regime, the Reynolds number limit for port A was specified according to “Equation (18)”.
(18)
while in the half pipe adjacent to port B, it was specified as
(19)
In Equations “(18)” and “(19)”, fA is the Darcy friction factor for turbulent flow in the half pipe adjacent to port A and fB is the Darcy friction factor for turbulent flow in the half pipe adjacent to port B.
The Darcy friction factor [16] for the turbulent to port A was modelled by utilizing Haaland equation for ports A and B as expressed in Equations “(20)” and “(21)” respectively.
(20)
The half pipe adjacent to port B
(21)
The Darcy friction factor formulates the friction losses in the pressurizer similar to the losses in pipe flow and other open-channel flow. The Haaland equation explicitly approximates experimental results of laminar and turbulent flow in pipe, open channels or pressurizer.
3. Methods
The pressure control system (Figure 1) sends a signal to the actuator valve, either to open or to shut the valve. The monitoring system compare the output value and apply proportional-integral-derivative (PID) values based on the errors from the feedback values to the summer, which adds the set-point figures to the error values.
The designed model of PRED is based on two-phase fluid domain which consists of elements, sources and sensors linked with fluid within the loop network. The schematic in Figure 2 incorporated three new sensors: 1) pressure & Internal energy sensor, 2) thermodynamic property sensor and 3) two-phase fluid property sensor. In addition, the controlled pressure source also known as pseudocritical rapid dissipation valve (PRDV), a physical signal (P) that control the pressure source and vent out excess energy buildup in the pressurizer in order to restore the transient state to the desire steady state operation set-point values.
The PRED control system designed model consist of two-phase (2P) fluid domain, linked with the fluid within the loop network, fulfilled the requirements for the configuration solver to compiled and initialized the model design, updated the library link blocks, evaluated the block parameters, construct the system equations for the physical network and initialized the design ready for simulation.
The PRED functions involve distinct interactions of the multivariable multiple-input multiple-output (MIMO) system is quite complex. There are 1) solid and liquid, where the liquid chemical content (particles) interact with the solid materials such as the pressurizer walls and the heating elements. 2) Solid and gaseous, here the solid materials of the pressurizer and the piping network (the carriers for steam) interact with the vapour. 3) Liquid and vapour, where the volume of the liquid has one phase and interact with the volume of vapour as the second phase. 4) The Combine internal components within the pressurizer where the chain effect of these interactions and other factors have significant influence on the pressurizer.
Figure 2. Pseudocritical rapid energy dissipation.
4. Results
The thermodynamic analysis of fluid properties has been widely studied in 1D graphics. However, such representation does not provide detail illustration compared to the 2-D contour and the 3-D graphics. The thermodynamic properties sensor in PRED provided the plot of the magnitude and phase of the temperature, specific enthalpy, specific volume and specific entropy against time at reference zero of the fluid in one-dimensional (1D) display (Figure 3).
The temperature values simulated for 10 seconds appeared to be constant, however, the range of temperature measured in degree Celsius ranges from 280 degree Celsius to 330 degree Celsius for specified time.
The spectrum analyzer plotted the noise power spectral density (Figure 4) that illustrated the peak as finder and distortion measurement of specified number of harmonics with the corresponding power measured in dBm/Hz. The
Figure 3. Plot of Thermodynamic properties.
Figure 4. Plot of distortion measurement.
device under test (DUT) was the pressure and internal energy sensor where selected peak finder provided corresponding peak values and frequencies was determined. The Gaussian reference plot provided complementary cumulative distribution function (CCDF), where the average maximum power in decibels was measured.
The Spectrum Analyzer displayed the results of the fundamental frequencies of the total harmonics distortion (THD), where the frequency variation was considered as the sum of the internal energy of the water, the steam produced, the pressure and the volume of fluid within the pressurizer. The distortion measured by six harmonics at different frequencies with different power levels measured in decibels equals −400.38 dBm.
4.1. Effect of Pseudocritical Saturation
The flow dynamics of the two-phase have effect on the fluid inertia, the viscous friction losses and the convective heat transfer within the pressurizer. The volume of the fluid is assumed to be constant, where the pressure and temperature of the fluid affect the thermal conservation of the pressurizer.
The proposed PRED function coordinates activities of the heaters to compensate for the heat loss of the pressurizer to achieve equilibrium between steam and water. The subcooled liquid has normalized internal energy defined as:
(22)
where, ū is the normalized internal energy of the fluid, u is the specific internal energy of the fluid, ūmin is the minimum specific internal energy and
is the specific internal energy of the liquid at saturation as shown in Equations “(22)” to “(24)”.
The changes in the internal energy of the pressurizer depends on specific volume of the liquid, the temperature and the pressure, due to the energy balance of the chemical reaction (Figure 5). The Pressure vector (0.001 and 100) coordinates the pressure values that provided the grid surface points. The pressure of the liquid is within the minimum values, why the pressure above the pseudocritical saturation approaches the maximum value.
At the critical pressure of 172 bar plus 22.0640 bar with atmospheric pressure of 0.101325 bar, where the critical point with saturated liquid internal energy of x = 2071 kJ/kg and pressure values of y = 24.77 bar (typical value 196.77 bar).
The specific volume (m3/kg) of the liquid (Table 1) represented the pressure vector and internal energy at specific grid point, where the boundary of saturation above or at the critical pressure is equal to vapour specific volume.
At the critical pressure the properties of the water gradually changes from liquid with small compressibility and high density to gaseous with low density and high compressibility. The physical properties of the liquid such as specific enthalpy, specific heat and density also changes with the changing temperatures, resulting in subcooled liquid and superheated steam as shown Equation (6) and Equations “(22)” and “(23)”.
Figure 5. Effect of Pseudocritical saturation on specific volume.
Table 1. Specific volume as variable.
The superheated vapour of the normalized internal energy was expressed as,
(23)
where, ūmax is the maximum specific internal energy and
is the specific internal energy of saturated vapour.
The specific entropy around the saturation boundary as per unit mass of liquid and vapour at saturation, where the liquid phase and vapour phase changes to mixture state in Figure 6.
Where the pressure vector at the upper part of the pressurizer registered 62.8 bar, the base remain at the minimum range with values of 0.001417 bar. While the specific internal energy ranges between 117 and 3933 kJ/kg/K and the specific entropy (kJ/kg/K) registered variable figures at different levels (Table 2).
The normalized internal energy of the two-phase mixture was evaluated as defined below:
(24)
Figure 6. Liquid and vapour mixer state.
Table 2. Internal energy of specific entropy.
For the initial fluid energy set for specific enthalpy with the default value of 1500 kJ/kg for the inertial specific internal energy. Equations “(22) to (24)” as the expression of normalized internal energy at all pressures where the minimum valid specific internal energy is minus one, the liquid saturation boundary is zero, the vapour saturation boundary is plus one, and the maximum valid specific internal energy is plus two.
Where the parameterization of the fluid properties in terms of pressure and normalized internal energy (ū), a linear transformation of specific internal energy (u), as the heat capacity and temperature of the subcooled liquid.
The internal energy axis of temperature in degree Celsius of the different regions of subcooled liquid and superheated vapor in two-phase mixture of saturated specific internal energy and pressure (Figure 7).
The Pseudocritical line or the Pseudocritical point is the point when the pressure above the critical pressure and the temperature reaches the maximum value of specific heat at particular pressure.
The interpolation of x, y and z values of the mesh grid where x = 2071 kJ/kg, y = 24.77 kJ/kg and Thermal conductivity level of 0.4047 mW/m/K represent the
Figure 7. Effect of temperature on internal energy and pressure vector.
pseudocritical point. Where the pressure values at point B and point C is equal to 19.63 bar above the 172 bar, with different specific internal energy values and different Prandtl numbers (Figure 8). Whiles the pressure vector values at points E, F and G with y values equals 0.002257 bars, with variable specific internal energy values.
The contour lines showing specific entropy (kJ/kg/K) that illustrate the dependence of the critical pressure cp and the critical volume cv. Where gases was considered to have constant specific heat over a specified range of temperatures.
The Point A serves as the pseudocritical rapid dissipation point (PRDP) where the specific internal energy value of x = 2071 kJ/kg and the pressure value for y = 24.77 bar plus 172 bar. Where the fluid becomes supercritical, the liquid and the vapour ceases to exist separate phases. The boundary between the merger of liquid and vapour represent the pseudocritical line. Where the specific internal energy values at points B, C and D correspond to 533.6; 2594 and 3810 kJ/kg respectively as shown in Table 3.
The liquid Prandtl number provided a table of matrix 25 × 100 that corresponds to pressure vector and normalized internal energy, where the saturated boundary above the critical pressure is clipped. These numerical values of specific internal energy and pressure vector correspond to the point x = 2071 kJ/kg and y = 24.77 bar plus 172 bar as the critical point (Figure 9).
The effect of Prandtl number on laminar flow regime of Newtonian fluid arises due to viscous stress and correlation of local strain overtime, as a heat transfer between the movement of fluid and that of the solid bodies of the pressurizer walls and the heating coil element.
Table 3. Thermal conductivity selected grid points.
Figure 8. Effect of thermal conductivity on internal energy.
Figure 9. Effect of Prandtl number on specific internal energy.
The Prandtl numbers along the pseudocritical saturation line remain infinity, while the level values at points g = 0.8909 and h = 5.036 dimensionless, with corresponding values of specific internal energy and pressure vector as shown in Table 4 and Table 5.
The resultant manipulations provided graphical illustrations of two dimension (2-D) contours where the pressure above the critical point provided pseudocritical line, the specific internal energy over the supercritical portion of the saturated internal energy critical value, pressure vector and other parameterization.
Effect of Pseudocritical Saturation in 3-D
The dynamics of the system is subject to temperature changes and the mass flow rate of the liquid into the pressurizer and the outlet of vapour from the pressurizer. The internal energy of the gas according to Waal’s is not a function of the specific volume. The total change in internal energy depends on the change in volume, temperature and pressure. The energy production in the fluid is due to the rate of change of chemical reaction, where the total energy balance depends on the first law of thermodynamic (energy transport).
Under steady state operation, the specific enthalpy seeks to balance the energy equation in terms of incoming and outgoing variable parameters. These changes include specific volume of the fluid and the pressure of the fluid or vapour, where low enthalpy corresponds to high mass liquid flow rate and high volumetric flow rate. There are other regimes such as supersonic and hypersonic velocity associated with gas flow from pressurizer or steam generator. The thermodynamic propertied of the liquid and vapour differ from each other.
The specific entropy as thermodynamic property, where the default substance (water) used to analyze the entropy per unit mass. The entropy determines the amount of heat transfer used to perform the task of steam generation within the Pressurizer. Here the specific entropy of the two phase as dependent of pressure, temperature and the level of the fluid in the Pressurizer.
The kinematic viscosity at point a remains as low as 0.1276 mm2/s, that at point e is extremely as high at the value of 26,940 mm2/s (Table 6). The kinematic
Table 4. Effect of Prandtl number above infinity.
Table 5. Prandtl number below threshold.
Table 6. Kinematic viscosity on Pseudocritical line.
viscosity of the fluid in the pressurizer is the inherent shear stress to flow, where the kinematic of liquid decreases with temperature and increases with temperature for gasses.
The thermal conductivity of specific liquid depends on the nature of that liquid, where thermal conductivity of water differs from other solid and liquids (Figure 10). The thermal conductivity of water changes due to temperature changes that correspond to the changes of atomic vibration frequency of the liquid.
In addition, the analysis of the thermal conductivity of liquid in NSSS pressurizer, where the thermal conductivity (W/m/K) values is a fraction of the thermal properties with the values from points a to e ranging between 0.4381 to 0.7361 as shown in Table 7.
The thermal conductivity of solid and liquids have correlation with the Prandlt number and that of Nusselt number. Depending on the laminar flow, the transitional flow and the turbulent flow, where the effect of pressure is minimal. The increase in temperature of the heater decreases the conductivity of the liquid due to the viscosity of the water, where the correlation of liquid and solid are function of temperature.
The thermal transport properties data provide specific heat, thermal conductivity and the Prandtl number that shows the large peak close to the critical point (Figure 11). The two-phase fluid properties offer visual representation of the function of pressure and specific internal energy. The Prandtl number (dimensionless) provided the correlation between viscosity of the fluid and that of the thermal conductivity.
The Prandtl number defines the momentum of the transport of the kinematic viscosity of the liquid as the ratio of the heat transport. At the critical point, the specific internal energy value is 2059 kJ/kg, while the critical pressure is 22.05 + 172 bar and the corresponding Prandlt number is 13.38.
The thermal transport properties of the critical pressure range near the critical point, where the peak value is large or sharp with the width as a fraction, which is either above or below the critical pressure, a region clip at peak values using selective option. A small pressure range around the critical point can be limited by clipping. The action of pseudocritical rapidly energy dissipate from pressurizer indicating the Prandtl number as thermal expansion coefficient, which truncate the effect of excess energy at the critical pressure point and prevent the water of the reactor coolant system from boiling.
Figure 10. Effect of transport thermal conductivity.
Figure 11. Pseudocritical rapid energy dissipation from Pressurizer.
5. Conclusions
The “Pseudocritical rapid energy dissipation (PRED) from Pressurizer in Nuclear steam supply system”, the study affirmed the significance of heat transfer coefficient between conductive and convective heat transfer with the Nusselt number as the ratio of convective heat transfer and conductive heat transfer. The Reynolds number provided the thermal diffusion as the ratio of the inertia forces to the viscous forces. While the Prandtl number offered the ratio of momentum diffusivity and the diffusivity boundaries layers in 2-D contour and 3-D illustrations.
The additional data provided from the lookup table for the parameterization of Density, Isothermal bulk modulus, Isobaric thermal expansion coefficient, Specific internal energy, Specific heat at constant pressure, Kinematic viscosity and Thermal conductivity plotted against Temperature, Pressure and specific internal energy with heat transfer coefficient (Prandtl) and pressure vector as correlation variables.
The thermodynamic properties of two-phase fluid characteristics change into subcooled liquid and superheated steam. At the boundaries when the fluid becomes supercritical as the liquid and the vapour ceases to exist in separate phases, where the merger of the liquid and vapour is the point of pseudocritical saturation, the liquid becomes subcooled and the vapour tends to be superheated.
The action of the controlled pressure source acting as the rapid dissipation control valve provided the Pseudocritical Rapid Energy Dissipation from the Pressurizer, the corresponding effect restores the system pressure back to the threshold values. Where the proportional heaters are energized and provide block of heat that changes the water into steam to reestablish the set-point values, for the safety of the Pressurizer as it maintains the reactor coolant system and keeps the water from boiling even at high temperatures.
Recommendation for future works in 2-D and 3-D Pressurizer level control, where the Boron Recycle System (BRS) and the Boron Thermal Regeneration System (BTRS) forms part of the nuclear steam supply system, parameters that determine the chemical concentration of boron and offers a balance in terms of chemical and volume control system (CVCS).
Acknowledgements
The data provided by the lookup table imbedded in Matlab/Simulink/Simscape (Acknowledgement 1994-2019 The MathWorks, Inc. ARPACK, Copyright © 2001, Rice University. Developed by D.C. Sorensen, R.B. Lehoucq, C. Yang, and K. Maschhoff). All rights reversed.