Consistent Fluid Property Evaluation in RELAP5-3D
Previous investigations have developed and tested a methodology to improve the consistency between used to determine the fluid property derivatives used in the numerical solution scheme and the solution scheme itself. Within a bounding rectangle of the fluid property table for a fluid state point, linear interpolants for both fluid property values and their derivatives are calculated. A weighted average between these interpolants the analytical approximations of the fluid properties and derivatives is formed with a weighting factor between 0.0 and 1.0. This average has shown potential to reduce mass error relating to the fluid properties in RELAP5-3D. This linear interpolation and weighted averaging methodology has been implemented in RELAP5-3D with two additional options that adjust the weighting factor in the average according to the location within the bounding rectangle. The first calculates the weighting factor based only on the position between the left and right internal energy endpoints. The second also factors the position between the upper and lower pressure grid points into the calculation of the weighting factor. The new coding passes verification testing and expanded testing of the various options. Graphs are presented from three kinds of test cases, a two volume basic model, Edward’s O’Brien blowdown, and typical PWR. Some trends of the effect of crossing fluid property table pressure and internal energy grid line are identified and discussed for the simple and Edwards test cases. Preliminary results show that it the consistent fluid methodology can reduce mass error in simple cases, but in more complex input models, there is no clear best.