Two?parameter dynamical scaling analysis of single?phase natural circulation in a simple rectangular loop based on dilation transformation

Jia?Ning Xu · Xiang?Bin Li · Yu?Sheng Liu · Shuai Yang · De?Chen Zhang · Qiao Wu

Abstract Scaling analysis is widely used to design scaled-down experimental facilities through which the prototype phenomena can be effectively evaluated. As a new method, dynamic system scaling (DSS) must be verified as a rational and applicable method.A DSS method based on dilation transformation was evaluated using single-phase natural circulation in a simple rectangular loop. The scaled-down cases were constructed based on two parameters—length ratio and dilation number—and the corresponding transient processes were simulated using the Relap5 computational code. The results show that this DSS method can simulate the dynamic flow characteristics of scaled-down cases. The transient deviation of the temperature difference and mass flow rate of the scaled cases decrease with increases in the length ratio and dilation number. The distortion of the transient temperature difference is smaller than that of the mass flow; however, the overall deviation is within a reasonable range.

Keywords Dynamical system scaling analysis · Single-phase natural circulation · Transient scaling deviation · Dilation transformation

1 Introduction

During the development and application of third-generation nuclear power technology, passive safety is an important factor to consider when designing safety facilities for nuclear power plants. In the AP1000 advanced passive Pressurized Water Reactor (PWR) system, designed by Westinghouse,the use of passive safety technology increases the safety margin of a nuclear power plant [1]. CAP1400 (China Advanced PWR 1400), a large-scale advanced nuclear power plant project in China, adopts passive safety technology similar to AP1000 [2]. This technology improves system safety over that of second-generation nuclear power plants, which are being widely built, while reducing cost and simplifying the structure. The HPR1000 (Hua-long Pressurized Reactor)nuclear power plant project independently developed by China also applies a reliable passive safety design concept[3].

Passive safety systems, such as passive core cooling and passive containment cooling, are widely used in thirdgeneration nuclear power plants and require rigorous safety performance assessments to verify the reliability and stability of the facilities before they are used. From an economic point of view, using full-scale prototype models to verify the reliability of passive safety features is unrealistic. Therefore,a validation experiment using a scaled-down experimental facility is a reasonable and widely used method. Many such scaled-down test facilities are already in use. Examples include the ACME (Advanced Core-cooling Mechanism Experiment) facility [2], APEX (Advanced Plant Experiment) facility [4], ATLAS (Advanced Thermal—hydraulic Test Loop for Accident Simulation) test loop [5], and ROSA(Rig of Safety Assessment) large-scale test facility [6].

In recent years, many scaling methods have been reported[7, 8]. Hsu et al. [9] analyzed the applicability of the linear scaling method based on the mass—energy balance equations and flow pressure drop relationship in the simulation of small-break LOCA accidents. Nahavandi et al. [10] derived the corresponding series of scaling criteria using three scaling methods: time-reducing scaling laws, time-preserving volumetric scaling laws, and time-preserving idealized model/prototype scaling laws. They discussed the rationality of these criteria for simulating and predicting the application of nuclear reactor systems. The power-volume method is suitable for establishing full-height test facilities [11]. Ishii et al. [12] deduced the corresponding scaling criteria for single-phase and two-phase natural circulation flows and indicated that it can be applied to experimentally simulate the natural circulation of pressurized water reactors in a LOFT(Loss Of Fluid Test) facility. Ishii developed a three-level scaling method to guide the design of an overall test facility[13]. This method obtains scaling criteria through both topdown and bottom-up scaling analysis. Zuber et al. [14, 15]developed Hierarchical Two-Tiered Scaling (H2TS) method suitable for more complex multiphase thermal experimental systems. By nondimensionalizing the conservation equations at multiple levels, a similar standard consisting of geometric and fluid parameter scaling factors was obtained. In addition,Zuber proposed the fractional scaling analysis method and applied it to analyze the loss of coolant accidents [16, 17].

Scale analysis based on the above scaling criteria often focuses on the static parameter deviation of the scaled model, whereas the actual situation it simulates is mainly a dynamic process. This type of simulation can lead to scaling errors if the actual parameters change over time throughout the dynamics. To address this defect, recent natural circulation simulation studies have sought improvement by focusing on transient processes. Kuran et al. [18] simulated the start-up transient of a boiling water reactor following a natural cycle at low pressure and flow using a PUMA (Purdue University Multi-Dimensional Integral Test Assembly facility). Lakshmanan et al. [19] confirmed that the Natural Circulation Test Facility (CNTF) model developed by Relap5 could predict transient cases of natural circulation under low pressure and power. Dixit simulated two-phase natural circulation instability in an integrated modular water reactor[20]. Basu used rectangular natural circulation to study the dynamic response to a change in the input power of a singlephase natural circulation system [21]. Reyes developed a new dynamic system scaling (DSS) method to simulate the dynamic process and evaluate the dynamic distortion of the scaled model [22]. In this scale-up process, five types of similarity criteria are generated, for which different arrays of similar criteria can be obtained. The dynamic deviation generated during the scaling process can be effectively evaluated by comparing the original dynamic process with the scaling process. Li et al. [23] used the DSS method to evaluate the scaling distortion of gravity-driven drainage systems and explored the dynamic changes of key parameters in the gravity drainage process under different scaled-down cases.In addition, a similarity criterion array based on the DSS identity transformation law for single-phase natural circulation is derived. Using a simple loop model, dynamical scaling analysis theory is enriched by comparing the difference between the simulation results of the DSS method and that of the traditional H2TS method [24].

To further understand the DSS method and expand its application in the field of natural circulation, a scaling criterion based on dilation transformation was deduced using the two-parameter affine transformation concept. Based on single-phase natural circulation in a simple rectangular loop,the corresponding two-parameter dynamical scaling analysis similarity arrays under different combinations of length ratios and dilation numbers are obtained. The Relap5 code is used to establish a natural-circulation loop model. Different scaling cases are developed according to the calculated similar array, and the dynamical scaling characteristics are evaluated under different length ratios and dilation numbers.

2 Two?parameter scaling analysis based on dilation transformation

2.1 Physical model

For the convenience of analysis, a simple rectangular loop model was built. The simplification of the model involves several assumptions.

(1) Flow is one-dimensional along the loop axis, and the fluid properties of the cross section remain uniform.

(2) Boussinesq approximation is applied to loops.

(3) Heat loss is negligible throughout the rest of the piping,except for the piping sections used to simulate heating and cooling.

(4) Pipes have the same diameter and cross-sectional area throughout the circuit.

A simple rectangular loop is illustrated in Fig. 1. The circuit is mainly composed of heating, cooling, ascending and descending pipes, and horizontal pipe sections. The heating section is used to simulate the heating of the water in the core of the nuclear reactor, whereas the cooling section is used to simulate the heat exchange and cooling process between the primary coolant circuit and the secondary side of the steam generator. The density of the heated fluid decreases, and the density increases after cooling,thereby realizing circulation in a rectangular loop. To create a dynamic process for analysis, in the heating section,the thermal component provides the circuit with an initial power of 12,000 W. In this case, the natural circulation in the loop reaches a steady state. The power is then increased by 5% at 12,600 W, after which the cycle gradually reaches the second steady state. The dynamic process to be analyzed lies between these two states. The main parameters of the simple rectangular natural-circulation loop are listed in Table 1.

Fig. 1 Sketch of the rectangular loop

Table 1 Prototype model parameters

2.2 Similarity criterion

Using the above simple rectangular loop model, a similarity criterion array based on the dilation transformation can be derived.

The momentum conservation equation for single-phase natural circulation is given by:

Here,Lis the total length of the circuit,Ais the crosssectional area of the pipe,Wis the mass flow rate,βfis the thermal expansion coefficient,ρis the fluid density, ΔTis the temperature difference at the outlet of the hot and cold sections,His the height difference between the centers of the cold and heat sources,fis the friction loss coefficient,dis the pipe diameter, andkis the flow loss coefficient.

Reyes provides a detailed description of the theoretical analysis of the DSS method [22]. Here, we use the DSS method to process the momentum balance equation above.

Define

In Eq. (11),tM,Wrepresents the model reference time(more details can be found in [25] and [26]), the subscripts M and P represent the corresponding parameters of the model and prototype, and R represents the ratio of the parameters of the model and prototype,

βTrepresents the dimensionless temperature difference between the hot and cold sections of the rectangular loop,andωTrepresents the corresponding change caused by the temperature change.

The constant scaling factor is given as

Clearly, the main parameter ratios of the scaled models are dilation numberλand length ratioLR. Forλ=1 , a scaling criterion based on the identity transformation of the DSS method can be obtained.

By selecting the range of dilation numbers and length ratios, different scaling ratios can be obtained by changing the combination of the two variables. The two-parameter scaling curves based on the above criteria array are shown in Fig. 2.

By fixing the dilation numberλ=0.5 , Fig. 2a shows the similarity criterion array formed by the parameters of the natural-circulation loop as the length scaling changes. For a fixed length ratioLR=0.5 , Fig. 2b shows the similarity criterion array formed by the loop parameters as the dilation number changes. The two graphs from this figure reflect the scaling law of the natural-circulation loop parameters from the two parameters, dilation number and length ratio. In general, the pipe diameter, heating power, and heat flux density of the scaled model increase with dilation number or length ratio. Considering the two cases, “fixed length ratio, dilation number gradually increases” and “fixed dilation number, length ratio gradually increases”, in the latter case, the rate of increase of the parameters of the loop is significantly greater than in the former case. This phenomenon is particularly evident when the scaling ratio is large. However, the same conclusions hold for other dilation and length scales.

Fig. 2 (Color online) Scaling number curve with fixed dilation number and length ratio

In addition, combining the two images reveals that the ratio of the heating power of the scaled-down model to that of the prototype case is smaller in the criterion array derived from the dilation transformation based on the DSS analysis(less than 0.15). This is conducive to the construction of a scale-model experimental bench, which would benefit from the smaller power. This is one of the advantages of DSS.

2.3 Numerical model

The simple rectangular-loop natural circulation model shown in Fig. 3 was established using the Relap5 computational code. Relap5 is a transient analysis program for lightwater reactors, developed by Idaho National Laboratory for the US Nuclear Regulatory Commission, and is widely used in thermal hydraulic analysis and calculation [27]. It should be noted that Relap5 uses a one-dimensional computational method for a simple natural loop.

Fig. 3 System nodalization

The heating pipe section was simulated using Pipe160,whereas the cooling section was completed by heat exchange between the Pipe110 and Pipe120 components. The heat of the fluid in the main circulation loop was removed by two time-dependent control bodies: TDV200 and TDV220 components. The distance between Pipe110 and Pipe160 represents the height difference between the heating and cooling sections. This provides the driving force for natural circulation in the circuit. In addition, the Pipe180 and TDV182 components are set on the loop to simulate the pressurizer to ensure that the pressure in the loop is maintained at the set value during the cycle. Considering that only the natural circulation process is simulated, no pump components are set up in this simple circuit. Before the calculation, the node sensitivity of the simple rectangular loop model is analyzed, which ensures that the calculation is independent of the number of nodes.

3 Results and discussion

According to the array of dilation transformation and scaling criteria, the two-parameter scaling criteria, based on the length ratio and dilation number, are calculated, summarized in Table 2.

To ensure the rationality of the scaled-down case setting,after length ratioLRis selected, dilation numberλshould be chosen to ensure that the ratio of each scaled-down parameter is not greater than 1.

First, we analyzed the static errors of the flow velocity and temperature difference between the prototype and model cases.

As shown in Figs. 4 and 5, compared with the prototype cases of natural circulation under steady-state cases, all cases based on the scaling of the two parameters can achieve good agreement regarding the flow parameters. The calculationerrors of the steady-state flow velocity and temperature difference between the hot and cold sections can be controlled within a level of no greater than 7%. This shows that the setting of the length ratio—dilation number two-parameter scaling case is reasonable, based on the above criteria.

Table 2 Two-parameter scaling number group based on DSS method

Fig. 4 (Color online) Error in flow velocity

Fig. 5 (Color online) Error in temperature difference

Figures 6 and 7 show the transient mass flow curves of the prototype case and each scaled-down model case.Based on the two-parameter analysis criteria, a comparison of transient mass flow rate is carried out from the two parameters. In Fig. 6, the length ratio is given, and the transient change in the mass flow rate of the model operating cases corresponding to different dilation numbers is compared. Figure 7 shows a comparison of the transient trend of the mass flow rate under different length ratios for a given dilation number. Notably, based on the DSS method, the processing time will be distorted for different scaled models, which is reflected in the time difference required for the fluid to circulate in a rectangular loop.To consider the influence of this process distortion on the analysis, the horizontal ordinate in the figure is defined ast?, which has the following relationship,

Fig. 6 (Color online) Transient mass flow rate comparison (length scaling)

Fig. 7 (Color online) Transient mass flow rate comparison (dilation number scaling)

wheret0denotes the time required for the fluid to complete a cycle. The ordinateβWrepresents the non-dimensional mass flow rate, andβWis defined as (W?W1,0)∕(W2,0?W1,0) .W1,0andW2,0are the mass flow rates at the initial and final stages, respectively.

The same method is used to obtain the transient change curve of the normalized temperature difference between the hot and cold sections for each case, shown in Figs. 8 and 9.βTis defined as (ΔT?ΔT1,0)∕(ΔT2,0?ΔT1,0) , where ΔT1,0and ΔT2,0are the temperature differences at the initial and final stages. Overall, the temperature difference shows a significantly different profile from that of the mass flow rate. Briefly, the transient temperature difference trends significantly in the early stages of the power step; that is,the normalized temperature difference jumps greatly to between 1.1 and 1.4 with the increase in power and then rapidly decreases to approximately 0.7 to 1.0. This oscillation process is completed in the first cycle of the loop after the power jump. After this time, the transient normalized temperature difference rises slowly, eventually reaching a steady state. Figures 8 and 9 show the evolution of the change in the normalized temperature difference as a function of the two parameters. The main variable in Fig. 8 is the dilation number, whereas it is the length ratio in Fig. 9.The conclusion is similar to that for the mass flow rate,in that increasing the dilation number and length scale produces simulation results of the scaled-down cases that match more closely with those of the prototype case. When the dilation number is small, the effect of the length scale on the results is less significant.

By comparison, the normalized temperature difference for the scaled-down cases containing a small length ratio and large dilation number (Fig. 8a—c) can be very close to that of the prototype case. However, for the normalized mass flow rate (Fig. 6a—c), even if a larger dilation number is selected,the time to reach the second steady state in the scaled-down cases is still significantly less than in the prototype case.This may be due to a change in the Reynolds number. A small length ratio corresponds to a small Reynolds number,which may cause errors in the calculation.

Fig. 8 (Color online) Transient temperature comparison (length scaling)

Fig. 9 (Color online) Transient temperature comparison (dilation number scaling)

Fig. 10 (Color online) Relative error of the mass flow rate (length scaling)

Fig. 11 (Color online) Relative error of the mass flow rate (dilation number scaling)

Figures 10 and 11 show the comparison of the transient deviation of the normalized mass flow rate. The dynamic deviation of the normalized mass flow rate is defined as(W?M?W?P)∕W?M×100% . Figure 10 shows the results with a fixed length ratio, comparing the scaled-down cases corresponding to different dilation numbers. During the first cycle, the normalized mass flow rate of all scaled-down cases oscillate significantly. For the same length ratio, a reasonably larger dilation number can reduce the error.AtLR=0.25 , if the dilation number is greater than 1.5,the maximum error during the oscillation can be less than 20%. AtLR=0.5 , the dilation number must be greater than 1.0 to achieve the same effect. This rule can be seen in Fig. 10a—c. For all cases, when natural circulation reaches 4.5 cycles, the transient error of the mass flow rate decreases to less than 10%, gradually reaching 0. Figure 11 shows the comparison of the scaled-down facility with different length ratios and a fixed dilation number. As shown in Fig. 11a, when the dilation number is selected to be small, the length ratio has little effect on the transient deviation of the mass flow rate. This is consistent with the conclusion drawn in Fig. 7a. For a fixed dilation number,selecting a larger length ratio can also reduce the transient error of the normalized mass flow rate, which can be seen in Fig. 11b—d. Of course, the length ratio should be selected in accordance with the cases of the scaled-down models.

The normalized deviations in the temperature difference between the hot and cold sections are similar to those of the mass flow rate (Figs. 12 and 13). Correspondingly,under a fixed length ratio, a larger dilation number yields a smaller error than that of the mass flow rate under the same parameters.

In addition, a comparison of the error between the normalized mass flow rate and the normalized temperature difference between the hot and cold sections shows an evident phenomenon: if a small dilation number is selected, the difference between the dynamic error curves corresponding to different length ratios is very small (Figs. 11a and 13a).The same conclusion is obtained for the normalized mass flow rate and normalized temperature difference between the hot and cold sections. When the dilation number is gradually increased, this phenomenon disappears (Figs. 11b—d and 13b—d). These results are consistent with the previous analysis of Fig. 11.

In the H2TS scaling analysis method, a series of criteria arrays was obtained by the dimensionless transformation of relevant control equations. In single-phase flow,this mainly includes the Richardson number, heat source number, and Stanton number. According to the scaling method, the theoretical ratios are all 1, and the specific scaling criteria are derived from these criteria. Therefore,in the H2TS method, the Richardson number, heat source number, and Stanton number are decisive criterion numbers. However, in the DSS method, the control equation of the system is treated using the entire method, by which the similarity criteria are obtained. This differs from the H2TS method.

Fig. 12 (Color online) Relative error of the temperature difference(length scaling)

Fig. 13 (Color online) Relative error of the temperature difference(dilation number scaling)

Fig. 14 (Color online) Steadystate error of the Richardson number

To facilitate comparison with the traditional scaling analysis method, the relevant criterion numbers are defined as follows:

Equations (41)—(43) represent the Richardson, heat source, and Stanton numbers, respectively. Notably, for the Richardson and Stanton numbers, their theoretical values are not 1 but will change with the length ratio and dilation number. The final scaling criteria were not determined based on their proportions. Therefore, in the DSS method,the Richardson number and Stanton number are considered non-decisive criterion numbers. The theoretical ratios of the criterion numbers for the different scaled-down cases are listed in Table 3. Here, the theoretical ratio of the heat source number is 1 for all scaled-down cases, whereas the Richardson number and Stanton number follow their respective ratios.

Fig. 15 (Color online) Steadystate error of the Stanton number

Fig. 16 (Color online) Steadystate error of the heat source number

Figures 14, 15 and 16 show the steady-state errors of the criterion numbers for each scaled-down case when the loop reached the second steady state. The error in the Richardson number is shown in Fig. 14. In all the cases, the static deviation is within ± 15%. The deviation in the Stanton number is smaller, as shown in Fig. 15. All scaled-down cases are controlled within ± 7%. It is easy to explain from Eqs. (41) and (43) that the error sources of the Richardson number and Stanton number are the deviation of the flow velocity and temperature difference between the hot and cold sections. An error analysis of multiple groups of independent variables reveals that the error of the Richardson number should be slightly larger. As shown in Fig. 16, the variation in the heat source number is very small compared to the first two criterion numbers, within ± 0.5%. This is because the number of heat sources, as a decisive parameter, is obtained during the derivation of the scaling criterion of the DSS method, and its theoretical ratio is always 1. The Richardson and Stanton numbers are composed of other criterion numbers. The error in the criterion number varies with the error of the independent variable. In addition, for each length scale that has been determined, there is a corresponding dilation number, making the steady-state heat source number of the scaled-down cases the same as that in the prototype case.

4 Conclusion

Based on the dilation transformation in the DSS method,the corresponding scaling criteria are deduced, and the twoparameter dilation scaling criteria array is obtained by considering the two parameters of the length scaling ratio and dilation number as the independent variables. Two-parameter scaled-down cases were constructed, and the natural circulation transient process was calculated using the Relap5 code. The static and dynamic deviations associated with the natural circulation flow parameters were evaluated, and the following conclusions were drawn.

(1) The scaling criterion obtained using the dilation transformation method can be used to effectively simulate single-phase natural circulation in a rectangular loop.The DSS method considers the distortion of the process time in the dynamic process simulation of the scaleddown model and can evaluate the dynamic process. In addition, the error in the static parameters of natural circulation is small.

(2) All scaled-down models based on two parameters,length scaling and dilation number, can simulate the variation trends of the mass flow rate and temperature difference under the prototype case. When the length ratio is fixed, the results obtained with the scaled-down cases are closer to those of the prototype case with an increase in the dilation number. When the dilation number is fixed, the scaled-down model case conform to the prototype case with an increase in the length ratio.Considering the factors of the two parameters, selecting a moderate length ratio and slightly larger dilation number can make the scaled-down model fit the prototype more closely. In the case of the same combination of dilation number and length ratio, the deviation of the temperature difference is smaller than that of the mass flow rate.

(3) Overall, the static deviations of the Richardson number,Stanton number, and heat source number can be maintained within a reasonable range. Because the Richardson and Stanton numbers are non-decisive criterion numbers, the corresponding errors are significantly greater than that obtained when the heat source number is used as the decisive criterion number.The above conclusions were obtained by analyzing the natural circulation of a simple rectangular loop using the dilation criterion of the DSS dynamical scaling method. To facilitate modeling, calculation, and analysis, ideal simplifications were made when the model was built using Relap5.Non-ideal cases such as pipe heat loss, wall thickness error,and two-parameter flow of fluid in the pipes need to be considered during the construction of the scaled model test bench. Additionally, there are other scaling criteria for the DSS method. It is especially noted that theβstrain transformation andωstrain transformation are essentially singlephase dilation transformations. Further comparative analyses can be performed for an in-depth evaluation of this approach.The generalization of the conclusions from the simple natural-circulation loop to the reactor system is also worthy of further analysis. Furthermore, two-phase natural circulation occurs frequently. In addition to single-phase circulation, the applicability of the DSS method to two-phase natural circulation should be further considered, such that this method can be better applied in actual scaling experiments.

Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Jia-Ning Xu and Xiang-Bin Li. The first draft of the manuscript was written by Jia-Ning Xu, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.