The development of a comprehensive annular flow modeling package for two-phase three-field transient safety analysis codes.

摘要

The goal of the current study was to assemble a physically-based and self-consistent annular flow modeling package that is amenable to implementation in three-field analysis environments and accurately captures the variation in entrainment and interfacial drag within co-current and counter-current regimes over the pressure range of interest (atmospheric to 2000-psia). The models that were incorporated into the newly proposed modeling package were either based on those developed in previous studies or developed uniquely within the current study. It is important to note that the current study used COBRA-TF to provide the baseline modeling package as a means for comparison and as a vehicle for assessing the newly proposed modeling packages; however, the proposed packages are amenable to implementation into any other three-field analysis tool.;The proposed modeling packages for co-current and counter-current annular flow are outlined in Chapters 4 and 5, respectively. The co-current modeling package: (1) applies an interfacial shear model that explicitly accounts for the presence of interfacial waves, (2) idealizes the structure of the interface in a manner that is consistent with both the interfacial shear model and other visual observations, (3) includes three mechanistic-based entrainment rate models (roll wave stripping, Kelvin-Helmholtz lifting, and liquid bridge breakup) that calculate a theoretical entrainment rate for a single wave based on the physical structures and controlling phenomena as they are currently understood for each mechanism, and (4) provides a functional relationship between the actual and theoretical entrainment rates based on comparisons to experimental data to account for any deficiencies that exist in the theoretical model. This methodology improves the physical basis of the modeling package while simultaneously leveraging the available experimental data to ensure the modeling package is able to accurately reflect the experimental data.;Meanwhile, the three-field Counter-Current Flow Limitation (CCFL) model developed in the current study is based on an empirical model that has been shown to suitably correlate specific sets of data over a wide range of flow path dimensions and geometries. The resulting correlation provides a quantitative description of the experimentally determined flooding curve. The proposed model compares the flow conditions predicted by the code to the results of the userspecified CCFL correlation to determine if the standard set of momentum equations should be replaced with a newly developed set of CCFL momentum equations. The proposed model also provides appropriate entrainment rate models (pool and excess film) and necessary criterion to exit the model in a stable manner. In general, this approach provides flexibility to the code user and again leverages the available experimental data to improve the predictive capability of the code since a universal model has yet to be determined for this phenomenon. While not entirely mechanistic, this approach ensures the proper amount of liquid flow can penetrate these regions, which is preeminent to achieving accurate predictions of coolant and temperature distributions for Loss-of-Coolant Accident (LOCA) scenarios. Overall the development of this model is a unique aspect of the current study because of the explicit treatment of the entrained field, which previously suggested models did not consider because they were aimed at two-field analysis environments.;The results of the current study indicate that the inclusion of these newly proposed modeling packages for both co-current and counter-current annular flow has provided increased accuracy in the predictions of phenomena that are of interest to reactor safety analyses. In particular, the mean relative error in entrained fraction was reduced from 20.2% (underprediction) to 4.5% (overprediction) and the mean relative error in axial pressure gradient was reduced from 108.2% to 7.6% (both overprediction) for co-current upward annular flow situations following the implementation of these packages into COBRA-TF and the code-to-data agreement of several different parameters within the counter-current flow regime was improved significantly. It was also shown that the proposed co-current annular modeling package: (1) provided reasonable estimates of a variety of more fundamental annular flow parameters such as wave spacing, velocity, and intermittency, and (2) was able to capture the general behavior within the developing flow region. Both these results provide confidence that the proposed modeling package reasonably reflects the underlying physics of the annular regime. Moreover, the current study is one of the few works that has examined the predictive capabilities of transient analysis codes within the developing, or non-equilibrium, annular flow region. (Abstract shortened by UMI.)