Modeling of UTSG in the Pressurized Water Reactor Using Accurate Formulae of Thermodynamic Properties
Abstract: In nuclear power plants, steam generators are important parts. This paper introduces a U-tube steam generator (UTSG) model implemented completely using MATLAB environment. The UTSG consists of four regions: upward and downward primary regions, upward and downward metal tube regions, and secondary regions, which contain heat transfer region, steam separation region, and subcooled water region. Governing equations are derived by applying energy and mass conservation equations in all regions. Accurate functions that describe the relationships between thermodynamic properties of the saturated steam are introduced instead of interpolation method that is widely used. Steady state and one transient case are presented as well.

1. Introduction

Steam generators (SGs) of the vertical recirculating type (UTSG) are the dominant types that are used in the pressurized water reactors (PWRs). SGs are heat exchangers containing thousands of tubes to transfer heat from the primary coolant which represents the hot side to the secondary coolant which represents the cold side to produce steam, that produce electricity by powering turbine generators. A large number of nuclear power plants (NPPs) have from 2 to 6 SGs; however some have up to 12 SGs [1]. In nuclear power plants overall the world, there are around 1300 SGs being in service.

Water/Steam thermodynamic properties are required in equation derivation of SG models. These properties are available in the form of plots or tables. Extracting data directly from these plots or tables is not suitable for dynamic modeling. Therefore, most of the researchers expressed them in linear mathematical forms, which can be effective for simulation purposes. In UTSG mathematical models, the relations of the specific volumes/densities of the vapor and liquid, and their specific enthalpy of the vapor and liquid and saturated temperature that depend on the saturation pressure (Psat) are required. Typically, such relationships can be found for pressure intervals in thermodynamic properties tables [2]. Many of these models develop approximating formulae using an interpolation method. In addition, the relationships of the derivatives of these parameters are also needed [3]. That is the reason that most of the models use linear approximating formulae.

Recently, many researchers reported the modeling of the UTSG using different platforms. For example, the derivation of the UTSG model was done using MATLAB as a part of a complete power plant simulator [4]. They used steam tables from Ref. [5] to get water/steam properties using interpolation functions programmed in C++ and FORTRAN languages. Reference [6] used an external software package to calculate the physical parameters of water and vapor. Then, a dynamic model for steady-state and transient responses was programmed on a personal computer using MATLAB/SIMULINK. While in Ref. [7], a dynamic model of AP1000 steam generator introduced using the following assumption: linear approximate profiles for volume weighted specific volume and specific enthalpy are used in secondary side tube bundle region. Other researchers reported the using of MATLAB to model Steam Generators as in [8] [9] [10] [11] [12] to do the same for calculating water/steam properties.

In the present work, a mathematical model for the UTSG is introduced using a home-made computer model written fully using MATLAB programming environment. It is well known that if there are many measurement data, MATRIXx environment provides a quick solution to the least square’s problems, so matrix calculations can be performed easily using commercial software packages such as MATHCAD and MATLAB [1]. MATLAB is widely used in modeling of steam generators because of its interactive programmable matrix calculator with graphics capabilities [13].

The Steam Generator is divided into four regions according to Ali’s second intermediate model (model C) for the UTSG, which is utilized in this work using a MATLAB Program [14]. The regions are: upward and downward primary regions which represent primary water pumped from the reactor cooling pump going up then down inside thousands of inverted U tubes, and during this movement heat transferred through upward and downward metal tube regions. Finally, the heat is transferred to the secondary regions where the cold feed water uses this heat to generate steam; then steam is separated during its way to turbine. Governing equations are derived by applying energy, mass, and momentum conservation equations on all regions.

The most popular design for steam generators is UTSG [1]. Figure 1 represents a schematic diagram of an inverted U-shaped vertical steam generator, in which all parts of the steam generator are contained within a cylindrical vessel and the pre-cooling water coming from the reactor core enters the inlet at a pressure of about 15 MPa facing upward at the inverted U. This flowing water transfers the heat to the secondary water in the steam generator through the pipe walls without mixing. The primary cooling water enters the steam generator at approximately 583 - 600 K and exits at 530 - 560 K. Also, the steam mass ratio varies from 17% - 33%. The higher this percentage, the higher the efficiency of the steam generator and the better for the turbine because this wet steam cannot be sent to the high -pressure turbine, which increases the need for high techniques to remove the liquid from the steam as much as possible [15].

Figure 1. U-tube steam generators for nuclear power plant [1].

2. The UTSG Model Equations

2.1. Assumptions

The first step in the UTSG modeling is to make appropriate assumptions. These assumptions are introduced to smoothing calculations and minimizing the complexity. The basic assumptions on which our model is constructed are:

1) One dimensional flow for primary and secondary cooling water.

2) Primary flow (Wpi) is constant.

3) Steam exit flow (Wso) is proportion to saturated pressure (Psat).

4) Thermal conductivity of the inverted metal tube is constant.

5) Heat transfer coefficients are assumed to be constant during transients.

6) Heat transfer between the tube bundle area and the downcomer is neglected.

7) Flow of exit Steam (Wso) = flow of feedwater (Wfi).

2.2. UTSG Parameters and Design State Values

The current mathematical model of UTSG is based on the geometrical parameters and the numerical constants for the four regions listed in Table 1 and Table 2, respectively [4]. In addition, Table 3 shows symbols and units of the variables used in the thermodynamic equations in the current model.

Table 1. Parameters of UTSG.

Table 2. Design steady state values of UTSG.

Table 3. Symbols and units of the variables used in the current model.

2.3. State Equations for UTSG Model

Nodalization scheme of the UTSG model is shown in Figure 2. The model governing the equations is derived by applying energy, mass, and momentum conservation equations.

Primary coolant equations derivation is done by applying the energy conservation equation to both PRIMARY_UP and PRIMARY_DOWN, which represent the moving of primary coolant up and down inside the inverted U-tube respectively (Figure 2).

$\frac{\text{d}}{\text{d}t}\left({M}_{p}\ast {C}_{p}\ast {T}_{1}\right)={W}_{pi}\ast {C}_{p}\left({T}_{p}-{T}_{1}\right)-{U}_{pm}\ast {A}_{pm}\ast \left({T}_{1}-{T}_{3}\right)$ (1)

${M}_{p}\ast {C}_{p}\ast \frac{\text{d}{T}_{1}}{\text{d}t}={W}_{pi}\ast {C}_{p}\ast {T}_{p}-{W}_{pi}\ast {C}_{p}\ast {T}_{1}-{U}_{pm}\ast {A}_{pm}\ast {T}_{1}+{U}_{pm}\ast {A}_{pm}\ast {T}_{3}$ (2)

${U}_{pm}=\frac{1}{\frac{1}{{h}_{i}}+\frac{{R}_{in}}{k}\ast \mathrm{log}\left(\frac{{R}_{avg}}{{R}_{in}}\right)}$ (3)

so

$\frac{\text{d}{T}_{1}}{\text{d}t}=\left(\frac{-{W}_{pi}}{{M}_{p}}-\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{p}\ast {C}_{p}}\right){T}_{1}+\left(\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{p}\ast {C}_{p}}\right){T}_{3}+\left(\frac{{W}_{pi}}{{M}_{p}}\right){T}_{pi}$ (4)

In the same way, we have

$\frac{\text{d}{T}_{2}}{\text{d}t}=\left(\frac{-{W}_{pi}}{{M}_{p}}-\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{p}\ast {C}_{p}}\right){T}_{2}+\left(\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{p}\ast {C}_{p}}\right){T}_{4}+\left(\frac{{W}_{pi}}{{M}_{p}}\right){T}_{1}$ (5)

Figure 2. Schematic diagram of the UTSG.

Also, U-Tube metal region equations derivation is done by applying the energy conservation equation on both METAL_UP and METAL_DOWN, which represent the inverted U-metal tube up flow and downflow respectively as shown in Figure 2.

$\frac{\text{d}}{\text{d}t}{M}_{3}\ast {C}_{m}\ast {T}_{3}={U}_{pm}\ast {A}_{pm}\ast \left({T}_{1}-{T}_{3}\right)-{U}_{ms}\ast {A}_{ms}\ast \left({T}_{3}-{T}_{s}\right)$ (6)

${U}_{ms}=\frac{1}{\frac{1}{{h}_{d}}+\frac{{R}_{out}}{k}\ast \mathrm{log}\left(\frac{{R}_{out}}{{R}_{avg}}\right)}$ (7)

$\frac{\text{d}{T}_{3}}{\text{d}t}=\left(\frac{-{U}_{pm}\ast {A}_{pm}-{U}_{ms}\ast {A}_{ms}}{{M}_{3}\ast {C}_{m}}\right){T}_{3}+\left(\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{3}\ast {C}_{m}}\right){T}_{1}+\left(\frac{{U}_{ms}\ast {A}_{ms}}{{M}_{3}\ast {C}_{m}}\right){T}_{s}$ (8)

In the same way, we have

$\frac{\text{d}{T}_{4}}{\text{d}t}=\left(\frac{-{U}_{pm}\ast {A}_{pm}-{U}_{ms}\ast {A}_{ms}}{{M}_{4}\ast {C}_{m}}\right){T}_{4}+\left(\frac{{U}_{pm}\ast {A}_{pm}}{{M}_{4}\ast {C}_{m}}\right){T}_{2}+\left(\frac{{U}_{ms}\ast {A}_{ms}}{{M}_{4}\ast {C}_{m}}\right){T}_{s}$ (9)

Secondary coolant region equations derivation is done applying energy and mass conservation equations on the three secondary coolant regions: Effective heat transfer region (EHT), upper tank region (ACCU_T), and Inlet annular tank region (IA_T) as shown in Figure 2.

EHT region of secondary cooling water is the area of heat accumulation in the metal tubes. Because of the heat build-up, the secondary cooling water is pushed up, and the energy conservation equation can be written as follows:

$\frac{\text{d}}{\text{d}t}\left({h}_{s}\ast {M}_{s}\right)={\stackrel{˙}{Q}}_{ms1}+{\stackrel{˙}{Q}}_{ms2}+{W}_{1}\ast {C}_{ps}\ast {T}_{5}-{W}_{2}\ast {h}_{e}$ (10)

$\begin{array}{c}{V}_{s}\frac{\text{d}}{\text{d}t}\left({\rho }_{s}\ast {h}_{s}\right)={U}_{ms}\ast {A}_{ms}\ast \left({T}_{3}-{T}_{s}\right)+{U}_{ms}\ast {A}_{ms}\ast \left({T}_{4}-{T}_{s}\right)\\ \text{\hspace{0.17em}}+{W}_{1}\ast {C}_{ps}\ast {T}_{5}-{W}_{2}\ast {h}_{e}\end{array}$ (11)

Now the following equations are deduced for the average thermodynamic properties of the secondary coolant in the effective heat exchange mass.

$L={L}_{sb}+{L}_{b}$ (12)

${T}_{s}=\frac{{T}_{5}+{T}_{sat}}{2}\ast \frac{{L}_{sb}}{L}+{T}_{sat}\ast \frac{{L}_{b}}{L}$ (13)

${\rho }_{s}=\frac{{\rho }_{d}+{\rho }_{f}}{2}\ast \frac{{L}_{sb}}{L}+\frac{{L}_{b}}{\left({\upsilon }_{f}+\frac{{x}_{e}}{2}\ast {\upsilon }_{fg}\right)\ast L}$ (14)

${h}_{s}={C}_{ps}\ast \frac{{T}_{5}+{T}_{sat}}{2}\ast \frac{{L}_{sb}}{L}+\left({h}_{f}+\frac{{x}_{e}}{2}\ast {h}_{fg}\right)\ast \frac{{L}_{b}}{L}$ (15)

${h}_{e}={h}_{f}+{x}_{e}\ast {h}_{fg}$ (16)

The heat transfer in the ACCU_T zone is not carried out immediately when the secondary coolant meets the metal tubes. It is in this region specifically that the secondary cooling water is in the state of saturation, as a mixture of water and steam. And dryers are used to separate the steam from the water to allow the steam to flow into the turbine while the rest of the water is returned to the inlet ring tank to be mixed with the feed water. Now by applying the continuity equations only to the mixture of water vapor and mass of vapor respectively, and by applying both the continuity equations and the energy balance to the mass of water only:

$\frac{\text{d}\left({V}_{r}\ast {\rho }_{r}\right)}{\text{d}t}={W}_{2}-{W}_{3}$ (17)

$\frac{\text{d}\left({L}_{dw}\ast {A}_{dw}\ast {\rho }_{f}\right)}{\text{d}t}={W}_{fi}+\left(1-{x}_{e}\right)\ast {W}_{3}-{W}_{1}$ (18)

$\begin{array}{c}\frac{\text{d}\left({M}_{dw}\ast {C}_{ps}\ast {T}_{6}\right)}{\text{d}t}={W}_{fi}\ast {M}_{dw}\ast {C}_{ps}\ast {T}_{fi}+\left(1-{x}_{e}\right)\ast {W}_{3}\ast {C}_{ps}\ast {T}_{sat}\\ \text{\hspace{0.17em}}-{W}_{1}\ast {C}_{ps}\ast {T}_{6}\end{array}$ (19)

$\frac{\text{d}\left(\left({V}_{dr}-{A}_{dw}\ast {L}_{dw}\right)\ast {\rho }_{g}\right)}{\text{d}t}={x}_{e}\ast {W}_{3}-{W}_{so}$ (20)

Subcooled water in IA_T region of UTSG is a mixture of both recycled water from the previous zone and the main feed water. By applying the energy balance:

$\frac{\text{d}\left({M}_{d}\ast {C}_{ps}\ast {T}_{5}\right)}{\text{d}t}={W}_{1}\ast {C}_{ps}\ast {T}_{6}-{W}_{1}\ast {C}_{ps}\ast {T}_{5}$ (21)

We can rewrite Equation (21) as:

$\frac{\text{d}{T}_{5}}{\text{d}t}=\frac{{W}_{1}}{{M}_{d}}\ast {T}_{6}-\frac{{W}_{1}}{{M}_{d}}\ast {T}_{5}$ (22)

From Equations (18) and (19):

$\frac{\text{d}{T}_{6}}{\text{d}t}=\left(\frac{{W}_{3}\left({x}_{e}-1\right)\ast {W}_{1-{W}_{fi}}}{{L}_{dw}\ast {A}_{dw}\ast {\rho }_{f}}\right){T}_{6}+\frac{{W}_{fi}}{{L}_{dw}\ast {A}_{dw}\ast {\rho }_{f}}\ast {T}_{fi}+\frac{\left(1-{x}_{e}\right)\ast {W}_{3}}{{L}_{dw}\ast {A}_{dw}\ast {\rho }_{f}}\ast {T}_{sat}$ (23)

From Equation (20):

$\begin{array}{l}{W}_{3}=\frac{1}{{x}_{e}}\ast \left\{\left(\left({V}_{dr}-{A}_{dw}\ast {L}_{dw}\right)\ast \frac{\text{d}{\rho }_{g}}{\text{d}P}+\frac{{\rho }_{g}}{{\rho }_{f}}\ast {A}_{dw}\ast {L}_{dw}\ast \frac{\text{d}{\rho }_{f}}{\text{d}P}\right)/\left(1+\frac{\left(1-{x}_{e}\right)}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\ast \frac{\text{d}{P}_{sat}}{\text{d}t}-\frac{{\rho }_{v}}{{\rho }_{L}}\ast \frac{\left({W}_{fi}-{W}_{1}\right)}{{G}_{2}}+\frac{{W}_{so}}{{G}_{2}}{A}_{g}\right\}\end{array}$ (24)

From Equations (17) and (24):

$\begin{array}{c}{W}_{2}=\left\{\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right\}\ast \frac{\text{d}{P}_{sat}}{\text{d}t}-\left\{\frac{{V}_{r}\ast {v}_{fg}}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right\}\ast \frac{\text{d}{x}_{e}}{\text{d}t}\\ \text{\hspace{0.17em}}-\frac{1}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\ast \frac{{W}_{fi}-{W}_{1}}{\left(1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\right)}+\frac{1}{{x}_{e}}\ast \frac{{W}_{so}}{\left(1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\right)}\end{array}$ (25)

From Equation (18):

$\frac{\text{d}{L}_{dw}}{\text{d}t}=\left(\frac{-1}{{\rho }_{f}}\ast \frac{\text{d}{\rho }_{f}}{\text{d}P}\ast \frac{\text{d}{P}_{sat}}{\text{d}t}\right)\ast {L}_{dw}+\frac{1}{{\rho }_{f}\ast {A}_{dw}}\ast \left({W}_{fi}+\left(1-{x}_{e}\right)\ast {W}_{3}-{W}_{1}\right)$ (26)

From Equations (10), (12), (24), (25), and (26), the state equations of both Psat and xe can be written as:

$\begin{array}{l}\frac{\text{d}{x}_{e}}{\text{d}t}=\frac{\left({V}_{s}\ast {G}_{5}+\left(\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)\right)\ast {U}_{ms}\ast {A}_{ms}}{{G}_{1}}\ast \left({T}_{3}+{T}_{4}\right)\\ \text{\hspace{0.17em}}+\left({W}_{1}\ast {C}_{ps}+{V}_{s}\ast {\rho }_{s}\ast \frac{{C}_{ps}\ast {L}_{sb}}{2\ast L}\ast \frac{{W}_{1}}{{M}_{d}}-{U}_{ms}\ast {A}_{ms}\ast \frac{{L}_{sb}}{L}\right)\ast {T}_{6}\\ \text{\hspace{0.17em}}+\left(\frac{\left({W}_{1}\ast {C}_{ps}+{V}_{s}\ast {\rho }_{s}\ast \frac{{C}_{ps}\ast {L}_{sb}}{2\ast L}\ast \frac{{W}_{1}}{{M}_{d}}-{U}_{ms}\ast {A}_{ms}\ast \frac{{L}_{sb}}{L}\right)\ast {G}_{7}}{{G}_{1}}\right)\ast {T}_{5}+{G}_{8}\end{array}$ (27)

and

$\begin{array}{l}\frac{\text{d}{P}_{sat}}{\text{d}t}=\left({V}_{s}\ast \left(\frac{{L}_{b}}{L}\ast \frac{\frac{{v}_{fg}}{2}}{{\left({v}_{f}+\frac{{x}_{e}}{2}\ast {v}_{fg}\right)}^{2}}\right)\right)\ast \frac{{U}_{ms}\ast {A}_{ms}}{{G}_{1}}\ast \left({T}_{3}+{T}_{4}\right)\\ \text{\hspace{0.17em}}-\left(\frac{\left({V}_{s}\ast \frac{{L}_{b}}{2L}\ast {\upsilon }_{fg}/{\left({\upsilon }_{f}+\frac{{x}_{e}}{2}\ast {\upsilon }_{fg}\right)}^{2}+{G}_{2}\right)\ast {V}_{s}\ast {\rho }_{s}\ast \left({C}_{ps}\ast {L}_{sb}\right)\ast {W}_{1}}{2\ast L\ast {M}_{d}\ast {G}_{1}}\right)\ast {T}_{6}\\ \text{\hspace{0.17em}}+{G}_{3}\ast {T}_{5}+{G}_{4}\end{array}$ (28)

where

$\begin{array}{l}{G}_{1}={V}_{s}\ast \left\{\frac{{L}_{sb}}{2\ast L}\ast \frac{\text{d}{\rho }_{f}}{\text{d}P}-\frac{{L}_{b}}{L}\ast \frac{\left(\frac{\text{d}{v}_{f}}{\text{d}P}+\frac{{x}_{e}}{2}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+\frac{{x}_{e}}{2}\ast {v}_{fg}\right)}^{2}}\right\}\\ +\left\{\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right\}\ast {V}_{s}\\ \ast \left({\rho }_{s}\ast \left(\frac{{L}_{b}}{2\ast L}\ast {h}_{fg}\right)-{h}_{s}\ast \left(\frac{{L}_{b}}{L}\ast \frac{\frac{{v}_{fg}}{2}}{{\left({v}_{f}+\frac{{x}_{e}}{2}*{v}_{fg}\right)}^{2}}\right)\right)\end{array}$

$\begin{array}{l}-{h}_{e}\ast \left(\frac{{V}_{r}\ast {v}_{fg}}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)+\left\{{V}_{s}\ast \left(\frac{{L}_{b}}{L}\ast \frac{\frac{{v}_{fg}}{2}}{{\left({v}_{f}+\frac{{x}_{e}}{2}\ast {v}_{fg}\right)}^{2}}\right)+\left(\frac{{V}_{r}\ast {v}_{fg}}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)\right\}\\ \ast \left\{{V}_{s}\ast \left({\rho }_{s}\ast {G}_{4}+{h}_{s}\ast {G}_{5}\right)+{h}_{e}\ast \left(\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)\right\}\end{array}$ (29)

${G}_{2}=\left(\frac{{h}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\ast \frac{{W}_{fi}-{W}_{1}}{1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}}-2\ast {U}_{ms}\ast {A}_{ms}\ast \left(\frac{{L}_{sb}}{2\ast L}+\frac{{L}_{b}}{L}\right)\right)\ast {T}_{sat}$ (30)

$\begin{array}{c}{G}_{3}=\frac{\left({V}_{s}\ast \left(\frac{{L}_{b}}{L}\ast \frac{\frac{{v}_{fg}}{2}}{{\left({v}_{f}+\frac{{x}_{e}}{2}\ast {v}_{fg}\right)}^{2}}\right)+\frac{{V}_{r}\ast {v}_{fg}}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)}{{G}_{1}}\\ \text{\hspace{0.17em}}\ast \left({W}_{1}\ast {C}_{ps}+{V}_{s}\ast {\rho }_{s}\ast \left(\frac{{C}_{ps}\ast {L}_{sb}}{2\ast L}\right)\ast \frac{{W}_{1}}{{M}_{d}}-{U}_{ms}\ast {A}_{ms}\ast \frac{{L}_{sb}}{L}\right)\end{array}$ (31)

${G}_{4}=\frac{{C}_{ps}\ast {L}_{sb}}{2\ast L}\ast \frac{\text{d}{T}_{sat}}{\text{d}P}+\frac{{L}_{b}}{L}\ast \left(\frac{\text{d}{h}_{f}}{\text{d}P}+\frac{{x}_{e}}{2}\ast \frac{\text{d}{h}_{fg}}{\text{d}P}\right)$ (32)

${G}_{5}=\frac{{L}_{sb}}{2\ast L}\ast \frac{\text{d}{\rho }_{f}}{\text{d}P}-\frac{{L}_{b}}{L}\ast \frac{\left(\frac{\text{d}{v}_{f}}{\text{d}P}+\frac{{x}_{e}}{2}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+\frac{{x}_{e}}{2}\ast {v}_{fg}\right)}^{2}}$ (33)

${G}_{6}=\left(\left({V}_{dr}-{A}_{dw}\ast {L}_{dw}\right)\ast \frac{\text{d}{\rho }_{g}}{\text{d}P}+\frac{{\rho }_{g}}{{\rho }_{f}}\ast {A}_{dw}\ast {L}_{dw}\ast \frac{\text{d}{\rho }_{f}}{\text{d}P}\right)/\left(1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\right)$ (34)

${G}_{7}={V}_{s}\ast {G}_{5}+\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}$ (35)

${G}_{8}=-\frac{\begin{array}{l}\left({V}_{s}\ast \left({\rho }_{s}\ast \left(\frac{{C}_{ps}\ast {L}_{sb}}{2\ast L}\ast \frac{\text{d}{T}_{sat}}{\text{d}P}+\frac{{L}_{b}}{L}\ast \left(\frac{\text{d}{h}_{f}}{\text{d}P}+\frac{{x}_{e}}{2}\ast \frac{\text{d}{h}_{fg}}{\text{d}P}\right)\right)+{h}_{s}\ast {G}_{5}\right)\begin{array}{c}\stackrel{}{}\\ \\ \end{array}\\ +{h}_{e}\ast \left(\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)\right)\ast \left({W}_{1}+\frac{1}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\ast \frac{{W}_{fi}-{W}_{1}}{1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}}\right)-{G}_{2}\ast {G}_{7}\end{array}}{{G}_{1}}$

$\text{\hspace{0.17em}}+\frac{\left({V}_{s}\ast \left({\rho }_{s}\ast {G}_{4}+{h}_{s}\ast {G}_{5}\right)+{h}_{e}\ast \left(\frac{1}{{x}_{e}}\ast {G}_{6}-\frac{{V}_{r}\ast \left(\frac{\text{d}{v}_{f}}{\text{d}P}+{x}_{e}\ast \frac{\text{d}{v}_{fg}}{\text{d}P}\right)}{{\left({v}_{f}+{x}_{e}\ast {v}_{fg}\right)}^{2}}\right)\right)-{h}_{e}\ast {G}_{7}}{{G}_{1}\ast {x}_{e}\ast \left(1+\frac{1-{x}_{e}}{{x}_{e}}\ast \frac{{\rho }_{g}}{{\rho }_{f}}\right)}\ast {W}_{so}$ (36)

2.4. Accurate Formulae of Steam/Water Properties

From Equations (15) to (36), we need to calculate the terms: ${\rho }_{V},{\rho }_{L},{h}_{V},{h}_{L}$ and ${T}_{sat}$, which are functions of Psat. Also we need their derivatives: $\frac{\text{d}{\rho }_{V}}{\text{d}{P}_{sat}}$, $\frac{\text{d}{\rho }_{L}}{\text{d}{P}_{sat}}$, $\frac{\text{d}{h}_{V}}{\text{d}{P}_{sat}}$, $\frac{\text{d}{h}_{L}}{\text{d}{P}_{sat}}$ and $\frac{\text{d}{T}_{sat}}{\text{d}{P}_{sat}}$.

We used MATLAB Curve Fitting Toolbox to provide an accurate functions or equations of steam/water properties taken from [16]. Also, the derivatives of functions are calculated easily. The MATLAB Curve Fitting Toolbox uses Linear and Nonlinear Regression to fit the data. It introduces extensive fitting capabilities. These functions are constant and can be used easily by itself or by their derivatives or integrals. These accurate functions introduced more accuracy for calculating of the thermodynamic properties of steam /water. Also, these functions give very fast computation, because steam properties are often called millions of times within the run of the model.

The formulae which are fitted are:

1) Vapor Saturated Density in terms of Psat, ρV(Psat) as shown in Figure 3 and Table 4.

2) Liquid Saturated Density in terms of PsatL(Psat) as shown in Figure 4 and Table 5.

3) Vapor Saturated Enthalpy in terms ofPsat,hV(Psat) as shown in Figure 5 and Table 6.

Figure 3. Fitting of saturated vapor density in terms of Psat.

Table 4. (a) General form of fitted curve of saturated vapor density in terms of Psat, (b) Coefficients of the equation with confidence Interval. (c) Goodness of fit.

Figure 4. Fitting of saturated liquid density in terms of Psat.

Table 5. (a) General form of fitted curve of Saturated Liquid Density in terms of Psat, (b) Coefficients of the equation with confidence Interval. (c) Goodness of fit.

4) Liquid Saturated Enthalpy in terms of Psat,hV(Psat) as shown in Figure 6 and Table 7.

5) Saturated Temperature in terms of Psat,,Tsat(Psat,) as shown in Figure 7 and Table 8.

It is clear that the results of the statistical data in Tables 4-8 show good compatibility according to four important statistical indicators; the first one is the (SSE), which represents the sum of squares due to the error and is considered a measure of the total deviation, the second one is R-Square, which is a measure of the variance of the data, the third is the Adjusted R-Square, which is closer to

Figure 5. Fitting of saturated vapor enthalpy in terms of Psat.

Table 6. (a) General form of fitted curve of saturated vapor enthalpy in terms of Psat, (b) Coefficients of the equation with confidence interval. (c) Goodness of fit.

Figure 6. Fitting of saturated liquid enthalpy in terms of Psat.

Table 7. (a) General form of fitted curve of Saturated Liquid Enthalpy in terms of Psat, (b) Coefficients of the equation with confidence Interval. (c) Goodness of fit.

Figure 7. Fitting of saturated temperature in terms of Psat.

Table 8. (a) General form of fitted curve of Saturated Temperature in terms of Psat, (b) Coefficients of the equation with confidence Interval. (c) Goodness of fit.

one, which indicates better fit to the data, and finally, the RMSE, which expresses the standard deviation.

3. State Variables of UTSG

The above Equations ((4), (5), (8), (9), (22), (23), (26)-(28)) define the following main UTSG variables:

· T1: temperature of the upward direction primary coolant in inverted U-Tube.

· T2: temperature of the downward direction primary coolant in inverted U-Tube.

· T3: temperature of the metal tube containing upward direction primary water.

· T4: temperature of the metal tube containing downward direction primary water.

· T5: recirculated water temperature before mixing.

· T6: annular tank temperature.

· Ldw: water level of recirculated water plus feedwater (mixing region).

· x: steam quality at the exit.

· Psat: secondary coolant saturation pressure.

Now we can arrange these equations in matrix form to be suitable for MATLAB as the following:

$\frac{\text{d}X}{\text{d}t}=A\ast X+B$ (42)

and state vector is:

$X=\left[\begin{array}{c}{T}_{1}\\ {T}_{2}\\ {T}_{3}\\ {T}_{4}\\ {T}_{5}\\ {T}_{6}\\ {L}_{dw}\\ x\\ {P}_{sat}\end{array}\right]$ (43)

and state matrix is:

$A=\left[\begin{array}{ccccccccc}{A}_{1,1}& 0& {A}_{1,3}& 0& 0& 0& 0& 0& 0\\ {A}_{2,1}& {A}_{2,2}& 0& {A}_{2,4}& 0& 0& 0& 0& 0\\ {A}_{3,1}& 0& {A}_{3,3}& 0& 0& 0& 0& 0& 0\\ 0& {A}_{4,2}& 0& {A}_{4,4}& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& {A}_{5,5}& 0& 0& 0& 0\\ 0& 0& 0& 0& {A}_{6,5}& {A}_{6,6}& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& {A}_{7,7}& 0& 0\\ 0& 0& {A}_{8,3}& {A}_{8,4}& {A}_{8,5}& {A}_{8,6}& 0& 0& 0\\ 0& 0& {A}_{9,3}& {A}_{9,4}& {A}_{9,5}& {A}_{9,6}& 0& 0& 0\end{array}\right]$ (44)

$B=\left[\begin{array}{c}{B}_{1}\\ 0\\ {B}_{3}\\ {B}_{4}\\ {B}_{5}\\ 0\\ {B}_{7}\\ {B}_{8}\\ {B}_{9}\end{array}\right]$ (45)

To make the model works independently, the model must be maintained in a closed-loop situation. This is done by adding the following two control conditions:

1) The flow rate of the exit steam is equal to the flow rate of thefeedwater (Wso = Wfi).

2) The flow rate of the exit steam is proportional to the saturated pressure of the steam (Wso/Psat = constant).

Now, by solving Equation (42), and taking into consideration the two control conditions, the UTSG model is presented. It was totally implemented in MATLAB environment. In the next section, some results of running the model in both steady state and transient state will be presented.

4. Results

Simulations are carried out by running the MATLAB model for 3 hours (1800 sec). Simulations are programmed on a personal computer (Intel Core™ i5-11600T Processor (12M Cache, up to 4.10 GHz)). The first step is using the initial design data from Table 2 to simulate the steady-state condition. Figures 8-11 show the steady state simulation of the thermodynamic properties T1, T2, T3, T4, Psat, and x, respectively.

Simulation of the transient performance of the UTSG is also done. As shown in Figure 12, the transient case increases the steam flow by 5% starting from initial time ti = 500 seconds tell final time tf = 1500 seconds. Figure 13 through Figure 16 show the transient state simulation of the response of thermodynamic properties T1, T2, T3, T4, Psat, and x. respectively.

5. Conclusion

An accurate mathematical modeling of the U-Tube steam generator, which is introduced in this work, offers a good start to assist a complete simulation for

Figure 8. The Tp-up (Tp-dn) inside the U-tube in steady state.

Figure 9. The Tm-up (Tm-dn) inside the U-tube in the steady state.

Figure 10. The Psat, in steady state.

Figure 11. The Xe, in steady state.

Figure 12. The steam flow rate.

Figure 13. The Tp-up (Tp-dn) inside the U-tube in Steam Flow Increase 5%.

Figure 14. The Tm-up (Tm-dn) inside the U-tube in Steam Flow Increase 5%.

Figure 15. The Psat, in steam flow increase 5%.

Figure 16. The Xe, in steam flow increase 5%.

the pressurized water reactor. The main advantage of using such highly precision expressions for deducing the thermodynamic relationships is the reduction in the simulation time as these functions provide a very fast calculation and control the safety usage of the PWR. This is because the properties of the vapor are called millions of times within the framework of the model while running and used the interpolation technique with the steam tables. The evidence of the high precision is shown through the statistical quantities computed in the fitting process with R-square being less than 1 in all cases, and can be easily used on their own or through their derivatives or integrals. These computed functions increase the precision for the calculation of the thermodynamic properties of steam/water. Moreover, this model was tested with both static persistence and 5% transient changes and gave typical and reasonable results that reflect the dynamic characteristics of the steam generator in real time. It shows a good adaptability in a wide range of industrial control.

Cite this paper: Alramady, A. , Al-Sharif, S. and Nafee, S. (2021) Modeling of UTSG in the Pressurized Water Reactor Using Accurate Formulae of Thermodynamic Properties. Journal of Applied Mathematics and Physics, 9, 947-967. doi: 10.4236/jamp.2021.95065.
References

[1]   Riznic, J. (2017) Introduction to Steam Generators from Heron of Alexandria to Nuclear Power Plants. In: Steam Generators for Nuclear Power Plants, Woodhead Publishing, Cambridge, 3-33.
https://doi.org/10.1016/B978-0-08-100894-2.00001-7

[2]   Nikolov, A., Pasheva, V. and Gerova, G. (2021) A Simulator of Vertical U-Tube Steam Generator in SCILAB Environment. AIP Conference Proceedings, 2333, Article ID: 090031.
https://doi.org/10.1063/5.0041750

[3]   Gerova, G. and Pasheva, V. (2017) A Simplified Mathematical Model of Vertical U-Tube Steam Generator in MATLAB Environment. AIP Conference Proceedings, 1910, Article ID: 020016.
https://doi.org/10.1063/1.5013953

[4]   Guimarães, L.N.F., Oliveira, N.S. and Borges, E.M. (2008) Derivation of a Nine Variable Model of a U-Tube Steam Generator Coupled with a Three-Element Controller. Applied Mathematical Modelling, 32, 1027-1043.
https://doi.org/10.1016/j.apm.2007.02.022

[5]   Haar, L., Gallagher, J.S. and Kell, G.S. (1984) NBS/NRC Steam Tables Thermodynamic and Transport Properties and Computer Programs for Vapor and Liquid States of Water in SI Units.

[6]   Dong, Z., et al. (2009) Dynamic Model for Control System Design and Simulation of a Low Temperature Nuclear Reactor. Nuclear Engineering and Design, 239, 2141-2151.
https://doi.org/10.1016/j.nucengdes.2009.05.006

[7]   Wan, J., Wu, S., Nuerlan, A. and Zhao, F. (2017) Dynamic Modeling of AP1000 Steam Generator for Control System Design and Simulation. Annals of Nuclear Energy, 109, 648-657.
https://doi.org/10.1016/j.anucene.2017.05.016

[8]   Ahmed, A.S.E., Elhosseini, M.A. and Arafat Ali, H. (2018) Modelling and Practical Studying of Heat Recovery Steam Generator (HRSG) Drum Dynamics and Approach Point Effect on Control Valves. Ain Shams Engineering Journal, 9, 3187-3196.
https://doi.org/10.1016/j.asej.2018.06.004

[9]   Li, X., Xu, E., Ma, L., Song, S. and Xu, L. (2019) Modeling and Dynamic Simulation of a Steam Generation System for a Parabolic Trough Solar Power Plant. Renewable Energy, 132, 998-1017.
https://doi.org/10.1016/j.renene.2018.06.094

[10]   Dong, Z. and Pan, Y. (2018) A Lumped-Parameter Dynamical Model of a Nuclear Heating Reactor Cogeneration Plant. Energy, 145, 638-656.
https://doi.org/10.1016/j.energy.2017.12.153

[11]   Sohn, J.J. and Seong, P.H. (2010) A Steam Generator Model Identification and Robust H∞ Controller Design with ν-Gap Metric for a Feedwater Control System. Annals of Nuclear Energy, 37, 180-195.
https://doi.org/10.1016/j.anucene.2009.11.005

[12]   Zhang, Y.S. and Ma, Y.Y. (2010) A Simplified Lumped Parameter Model for U-Tube Steam Generator. 2010 International Conference on Electrical and Control Engineering, Wuhan, 25-27 June 2010, 253-256.
https://doi.org/10.1109/iCECE.2010.68

[13]   Ordys, A.W., Pike, A.W., Johnson, M.A., Katebi, R.M. and Grimble, M.J. (1994) Process Models. In: Modelling and Simulation of Power Generation Plants, Springer, London, 117-215.
https://doi.org/10.1007/978-1-4471-2114-5_4

[14]   Ali, M.R.A. (1976) Lumped Parameter, State Variable Dynamic Models for U-Tube Recirculation Type Nuclear Steam Generator. PhD Dissertation, The University of Tennessee, Knoxville.

[15]   Kolev, N.I. (2015) Multiphase Flow Dynamics 5. Springer International Publishing, Berlin.
https://doi.org/10.1007/978-3-319-15156-4

[16]   Release on the IAPWS Formulation 2011 for the Thermal Conductivity of Ordinary Water Substance, September 2011.
http://www.iapws.org/relguide/ThCond.html

Top