A Method of Limiting Intermediate Values of Volume Fraction in Iterative Two-Fluid Computations

Abstract

Multidimensional computational analysis of fluid flow is usually done by segmented iterative methods, as the equations sets generated are too large to permit simultaneous solution. Frequently the need arises to compute values for variables which must remain bounded for physical reasons. In two-phase computation, for example, the volume fraction is restricted to values between 0 and 1, but iterative procedures often return intermediate values which violate these bounds. It is fairly straightforward to prevent negative values, however no satisfactory method of imposing the upper limit has been published. A method of smoothly applying the limit in reversible fashion is outlined in this note.