Journal of Applied Mathematics

Volume 2010 (2010), Article ID 476839, 19 pages

http://dx.doi.org/10.1155/2010/476839

## A Characteristic Analysis of One-Dimensional Two-Fluid Model with Interfacial Area Transport Equation

Department of Mechanical and Aerospace Engineering, The Ohio State University, 201 West 19th Avenue, Columbus, OH 43210, USA

Received 6 April 2010; Accepted 8 October 2010

Academic Editor: A. A. Soliman

Copyright © 2010 Xia Wang and Xiaodong Sun. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

An interfacial area transport equation (IATE), proposed to dynamically describe the interfacial structure evolution of two-phase flows, could help improve the predictive capability of the two-fluid model. The present study aims to investigate the well-posedness issue of a one-dimensional two-fluid model with the IATE (named “two-fluid-IATE model” hereafter) using a characteristic analysis. The momentum flux parameters, which take into account the coupling of the volumetric fraction of phase and velocity distributions over the cross-section of a flow passage, are employed. A necessary condition for the system to achieve hyperbolicity under an adiabatic flow condition is identified. A case study is performed for an adiabatic liquid-liquid slug flow, which shows that the hyperbolicity of the two-fluid-IATE model is guaranteed if appropriate correlations of the momentum flux parameters are applied in the two-fluid-IATE model.

#### 1. Introduction

In modeling of two-phase flows, it is realized that the one-dimensional two-fluid model is widely applied in many computer codes [1, 2]. Also, interfacial area transport equations (IATEs), which could capture the interfacial structure evolution through modeling the bubbles/droplets interactions (coalescence and/or disintegration) as well as phase change (boiling, evaporation, or condensation) [3–9], are developed in the framework of the two-fluid model. These IATEs belong to two categories: one-group and two-group models. The one-group IATE models small fluid particles, such as spherical bubbles, and therefore is only applicable to bubbly flows, while the two-group IATE takes into account both small and large fluid particles for applications in flows beyond bubbly flow. For simplicity, the two-fluid model with IATEs is named “two-fluid-IATE model” hereafter.

A successful model proposed for solving physical problems must first satisfy four principles, that is, equipresence, well-posedness, frame invariance, and fulfillment of the second law of thermodynamics. Mathematically, a well-posed problem means that the solution should exist, be unique, and depend continuously on variables. Insights into the field of the well-posedness of the two-fluid-IATE model are of great significance. In contrast to the single-phase Eulerian model, which has real characteristic roots inherently guaranteed by the second law of thermodynamics, the standard form of the partial differential equations in the two-fluid model may not be strictly hyperbolic as an initial value problem [10, 11]. Moreover, the infinite instability growth rate may exist when two-phase flows are modeled with the two-fluid model, which violates the observation that the energy originating from perturbation continuously decays due to the local fluctuations (energy cascade and dissipation), resulting in the naturally steady equilibrium of the flow in different flow regimes. These lead to the possibility of the partial differential equations of the two-fluid model being ill-posed.

In the literature, two fundamental technical approaches were utilized to remedy the ill-posedness of the two-fluid model [12, 13]. The first approach is more realistic, in which additional terms representing physical mechanisms, such as the pressure difference between the phases [14], virtual mass force [15], surface tension force [16, 17], and Reynolds stress [18], are incorporated into the conservation equations. However, difficulties were encountered in both constructing reliable constitutive relations for these terms and resolving the short wavelengths at which these physical mechanisms became dominant. The other approach is to regularize the model mathematically or numerically. For example, the ill-posedness of the two-fluid model can be remedied by either applying a preconditioning mass matrix [19] or adding an artificial numerical damping [12]. The disadvantage of the second approach lies in the possibility of introducing nonphysical dynamic behaviors. Furthermore, Song and Ishii [20] reported that the stable region of a well-posed one-dimensional two-fluid model could be evaluated by introducing momentum flux parameters, which helped preserve the information of the local flow structures, in which the void fraction profile is coupled with the velocity profile over the cross-section of a flow passage.

The current study aims at investigating the hyperbolicity of the one-dimensional two-fluid-IATE model by employing the momentum flux parameters. A characteristic analysis of the partial differential equations in the one-dimensional two-fluid-IATE model for an adiabatic flow was conducted to determine the necessary condition of the stability. The criterion was then developed in terms of the momentum flux parameters and local flow variables, such as the void fraction and phase velocities. Finally, a case study was performed for an adiabatic liquid-liquid two-component slug flow, in which the momentum flux parameters are correlated with the total void fraction of the dispersed phase based on experimental data [9]. The characteristic roots of this two-fluid-IATE model are examined graphically under different flow conditions to reveal the effects of the momentum flux parameters on regularizing the two-fluid-IATE model.

#### 2. Governing Equations of the Two-Fluid-IATE Model

The current paper focuses on discussing the well-posedness of the standard one-dimensional two-fluid-IATE model, which can be obtained directly by area averaging the three-dimensional two-fluid-IATE model. To quantify the coupling effect of the local void fraction and velocity distributions, the momentum flux parameter is defined for field- [20]: where is the velocity of field-. The signs and are introduced to represent area averaging and phase fraction-weighted area averaging, respectively, which are defined for a parameter as and , respectively. Here, is the flow cross-sectional area, and is the void fraction of field-. In the conventional one-dimensional two-fluid model, the momentum flux parameter is assumed to be unity, resulting in an ill-posed one-dimensional two-fluid model.

Dropping the signs and for simplicity, the area-averaged continuity and momentum equations of the two-fluid model and two-group IATEs for an adiabatic incompressible two-phase flow can be written as
Here, the subscripts 1, 2, and 3 denote field-1 (group-1 dispersed fluid particles), field-2 (group-2 dispersed fluid particles), and field-3 (continuous phase), respectively. Here, ,,,,,,, and are the pressure, density of field-, net intergroup mass transfer rate from group-1 to group-2 fluid particles, generalized drag force of field*- *, coefficient accounting for the inter-group void transport at the group boundary, critical bubble diameter at the group boundary, fluid particle Sauter mean diameter, and interfacial area concentration source/sink rate for field-, respectively. In order to highlight the roles that the momentum flux parameters play in the process of regularization, we assume that the density of each phase does not change appreciably, both spatially and temporally.

In a gas-liquid two-phase flow, group-1 and group-2 fluid particles are small (spherical and distorted) and large (cap, slug, and churn-turbulent) bubbles, respectively [6, 7]. Therefore, it is reasonable to assume that * _{}*is identical to

*, since field-1 and field-2 fluid particles are the same fluid with different particle sizes. In addition, there is no mass transfer considered between two different kinds of fluids due to the adiabatic assumption. In other words, no mass transfer takes place either between field-1 and field-3 or between field-2 and field-3. It is also noted that the phase-to-interface pressure difference is neglected, which means a single pressure is shared by all the phases.*

_{}In the literature, the contributions of to the regularization of two-fluid model have attracted attention of many researchers [15, 16]. It was suggested that the inclusion of , including the virtual mass force and drag force, could help render the model into a hyperbolic form. However, as emphasized by Song and Ishii [20], to construct an appropriate correlation of that could be applicable in any circumstances remained a challenge. Noting the argument that usually tends to stabilize the two-fluid model, the analysis of the model hyperbolicity would be conservative if one adopts an algebraic form for since the contribution of to stabilizing the model will then not be credited in the characteristic analysis.

#### 3. Characteristic Analysis

A characteristic analysis for the one-dimensional two-fluid-IATE model was performed to study the hyperbolicity of the model. The analysis provides information on the limiting short wavelength behavior. If a vector is defined, then the governing equations of the system can be rewritten in the form of matrix as where and are the coefficient matrices. The detailed matrices of , , and are given in the appendix. The characteristic analysis of (3.1), which is a set of quasilinear partial differential equations, is analogous to the analysis of linear partial differential equations. The necessary condition for the model to be a well-posed initial value problem is that the set of partial differential equations is hyperbolic, that is, the characteristic equation has real roots for . The complex conjugate values for the eigenvalue lead to the elliptic equations, which are ill-posed mathematically. In this analysis for the system of (2.2), is rearranged and simplified to where the corresponding polynomials are given as Note that is replaced by since four out of the eight roots for are guaranteed to be real roots.

Considering the fact that and are positive and the constrains of being between zero and unity to maintain its physical meaning, we recognize that the coefficient of the highest order of , that is, the 4th order, should be always positive. All roots of are real if the following two requirements are satisfied:(1) has three real roots , and (let );(2), , and .

It is realized that, to solve the eigenvalues, a 3rd-order polynomial equation is required to be solved, from which the expression of its roots are usually very complicated and lengthy. Therefore, this analytical approach is not taken here.

Alternatively, based on available experimental data of a certain two-phase flow configuration, the momentum flux parameters could be constructed as a function of the void fraction (and perhaps other flow variables), reducing the numbers of undetermined variables in . The plots of could therefore be given, from which it becomes straightforward to study graphically the hyperbolicity property of . If there are four intersections between and the real axis, then the characteristic equation will have four real roots for . More details are given in the next section.

#### 4. Results and Discussions: A Case Study

To demonstrate the feasibility of the studies discussed above, is analyzed for a specific two-phase flow as a case study. Since the compressibility of fluids is not considered here for simplicity, it is appropriate to take a liquid-liquid slug flow as an example. One relevant experiment in the literature was performed by Vasavada [9]. In this experiment, two immiscible fluids, namely, water and Therminol 59, flowed upward adiabatically in a vertical pipe with an inner diameter of 25 mm at the room temperature. Therminol 59, an organic engineering fluid, was the dispersed phase with a density of 972 kg/m^{3} at 20°C. Based on the sizes of droplets, the organic droplets were classified into field-1 and field-2, with small droplets being field-1. Water was the continuous phase as field-3 with a density of 998 kg/m^{3} at 20°C. Local flow variables in 17 different flow conditions were measured at three measurement elevations along the flow directions with , 67, and 121, respectively [9]. Table 1 lists the cross-sectional area-averaged data for each flow condition at .

The momentum flux parameter (, and 3) are calculated according to the experimental data [9] when the total void fraction of the organic drops, ), is in the range between 0.1 and 0.6. The results are given in Figure 1, where are plotted against *α*. From the plots, the relations between and *α* could be determined. Figure 1(a) illustrates that is a relatively weak function of *α*. Most of data points fall into the range around 1.5 except two points having a magnitude between 2 and 2.5. After a further investigation of the measured data, it is found that the experimental data from the measurements in the flow conditions associated with these two data points had zero values for group-2 void fraction at several radial locations. This indicates that the flow might be in the bubbly-to-slug transition region instead of slug flows. Considering the difficulties in obtaining accurate measurements in flow transition regions, it might be reasonable to discard those two points and assume that . In Figure 1(b), although it shows a scattered distribution when *α* is small, follows the power law profile of *α* to a good approximation as *α* becomes larger, which corresponds to the slug flow conditions of interest. Other data-fitting methods of constructing correlation of , such as a linear function or higher order polynomials, are also investigated, which generally provide smaller value of , which is a measure of how good the curve fit is with the data and is bounded between zero and unity. A higher value of always provides a better fit of the regression to the data set. Finally, as shown in Figure 1(c), fits a 2nd-order polynomial of *α* quite well.

As a result of the above analysis, three assumptions are made as follows:
Here, *α* stands for the total void fraction of Therminol 59 drops. It should be noted that for a more general analysis, the momentum flux parameters may be better correlated with additional parameters.

In addition, it is reasonable to assume that the velocities of the different phases are the same, that is, , because of the similar densities between water and Therminol 59. (the density difference between the two fluids is about 2.7% at 20°C.) Since the velocity for each phase follows similar turbulent velocity profile along the radial direction of the pipe, this assumption could also be partially verified from comparisons of the area-averaged velocities between group-1 and group-2 drops, as shown in Table 1. In most runs, is approximately equal to , except for three runs, that is, Runs 16, 20, and 27. In these three runs, it is shown that the void fraction of group-2 droplets is significantly smaller than that of group-1 droplets. Also, the average velocity of the dispersed phase is found to be close to that of group-1 drops. These indicate that in Runs 16, 20, and 27, the numbers of group-2 drops measured are quite small, and the flows fall into bubbly-to-cap bubbly transition region. The zero relative motion assumption reduces the two-fluid-IATE model into a nonslip two-fluid model, similar to the homogeneous equilibrium model (HEM).

Furnishing the above assumptions, that is, (4.1) and equal phase velocities, the characteristic function is simplified to a function of the void fraction of group-1 , void fraction of group-2 , characteristic root , and velocity , only. In addition, if one introduces a new parameter, as the dimensionless velocity, and defines it as , then the characteristic equation converts to . The roots of have the same mathematical properties as those of , that is, if has four real roots, so does the original equation .

Figures 2, 3, 4, 5, 6, 7, and 8 provide plots of versus for 17 flow conditions. In these figures, subplot (a) is a plot in a relatively large range of , while subplot (b) shows a plot in a sufficiently narrow range of for the benefit of observing the two roots near which has small variations in the real axis. Figures 2–8 show clearly that has four positive roots, which are consistent with the physical meanings. As illustrated in the plots, the typical shape of is that has one pair of roots that are close to each other and are approaching a less-than-unity value, whereas the other pair intersects at the values of in the order of unity. Attention should be paid to the first pair of roots, which demonstrates an essential role on the stability. For the flow conditions investigated, it is found out that the modified two-fluid-IATE model satisfies the necessary condition of the well-posedness, which indicates that the flows tend to be stable. Therefore, it is mathematically consistent with the HEM.

The advantage of the two-fluid-IATE model with the momentum flux parameters over the conventional two-fluid model can be shown if the relative velocity between the phases is considered. Based on a force balance, the relative velocity between the dispersed phase and the continuous phase, , is obtained as [21]
where , , , and are the bubble diameter, gravitational acceleration, drag force coefficient, and density of the continuous phase (field-3), respectively. It should be noted that the bubbles/drops are assumed spherical in (4.2). The drag force coefficient for distorted bubbles/drops is employed for field-1 as [22]
where is the surface tension. Assuming that the diameter of the slug bubbles is 80% of the pipe inner diameter, the drag force coefficient for slug bubbles/drops is used for field-2 as [22]
As a result, the relative velocity is a function of the total void fraction *α* as shown in Figure 9.

Our study shows that the eigenvalue becomes complex for the one-dimensional two-fluid three-field model with the unity momentum flux parameters if the relative velocities given by (4.2) are taken into account, indicating that the model is ill-posed (no hyperbolic region). In contrast, if the two-fluid-IATE model with the momentum flux parameters is considered, then the stability criteria, which ensures the hyperbolicity of the model, is shown in Figure 10 as a map with the continuous phase velocity as the -axis and group-1 void fraction as the -axis. When the total void fraction is smaller than 0.45, the model is always hyperbolic no matter what the continuous phase velocity is. When the total void fraction is between 0.45 and 0.65, the continuous phase velocity must increase to ensure the model hyperbolicity as the group-1 void fraction increases, indicating that the hyperbolic region becomes smaller. When the total void fraction is larger than 0.65, the model is unstable, which indicates the same mathematical property of the conventional two-fluid model.

Furthermore, the stability criteria can be represented using the slip ratios for field-1 and field-2 dispersed phase with respect to the continuous phase velocity, which are, respectively, defined as The stability criteria for and are plotted in Figure 11. In the current case study for the upward liquid-liquid flow, the dispersed phase is lighter than the continuous phase; therefore, the slip ratios defined by (4.5) are always larger than unity. It is clearly seen that, except the case when the total void fraction is 0.45, the region of hyperbolicity decreases as the group-1 void fraction decreases. When the total void fraction reaches 0.65, the slip ratios slightly larger than 1, which is close to the HEM in order to ensure the model hyperbolicity. When the total void fraction is beyond 0.65, there is no hyperbolicity region if the momentum flux parameters given by (4.1) are adopted. Note that, to ensure model hyperbolicity, the criteria for and shown in Figure 11 should be satisfied simultaneously.

#### 5. Conclusions

This study examines the hyperbolicity property of the one-dimensional incompressible two-fluid-IATE model. In this study, momentum flux parameters are applied to take into account the effect of the void fraction and velocity profiles. These momentum flux parameters have physical meaning based on the local flow distributions and do not introduce any artificial terms. In addition, the proposed two-fluid-IATE model can be applied dynamically over a wide range of flow regimes, up to the churn-turbulent flow, including the flow transition regions. As a case study, a characteristic analysis has been performed to obtain the necessary conditions of the well-posedness of the one-dimensional incompressible two-fluid-IATE model for liquid-liquid slug flows in a circular pipe. Based on the available experimental data, the constitutive correlations of the momentum flux parameters were approximated as functions of the total void fraction of the dispersed phase. It was shown that the inclusion of the momentum flux parameters in the two-fluid-IATE model helps the model maintain hyperbolicity. The model’s hyperbolicity depends heavily on local flow parameters, such as the void fraction and phase velocities. When the relative velocity is considered, stability criteria pose restrictions on the continuous phase velocity to achieve model hyperbolicity.

This analysis only provides the necessary condition of the well-posedness of the two-fluid-IATE model. On the other hand, the question remains whether the instability growth rate for perturbation is unbounded as the amplitude of the perturbation decreases to zero. This requires a study of the dynamic behavior using the linear stability theory. The sufficient condition could be provided by using a wake equation with an accompanying infinitesimal perturbation [11, 12, 23] but has not been conducted for the current liquid-liquid applications. In addition, it remains questionable as how to generalize the results of incompressible model for the applications of gas-liquid two-phase flows, where the gas compressibility needs to be considered. The effect of the compressibility of the gas phase is essential for the validity of two-fluid-IATE applications in general two-phase flows [24]. Despite the aforementioned considerations, the introduction of the momentum flux parameters appears promising to the regularization of two-fluid-IATE

#### Appendix

The detailed matrices of , , and are given as follows: where denotes , denotes , denotes , denotes , denotes , and denotes .

#### Acknowledgment

This work was performed under the auspices of US Nuclear Regulatory Commission (NRC) through a Faculty Development Grant (Award no.: NRC-38-08-945). The authors would like to thank the NRC for the support.

#### References

- Safety Code Development Group, “TRAC-PF1/MOD1, an advanced best-estimate computer program for pressurized water reactor thermal-hydraulic analysis,”
*NUREG/CR-3858*LA-10157-MS, Los Alamos National Laboratory, 1986. View at Google Scholar - RELAP5 Code Development Team, “RELAP5/MOD3 Code manual,”
*NUREG/CR-5535*INEL-95/0174, Idaho National Engineering Laboratory, 1995. View at Google Scholar - G. Kocamustafaogullari and M. Ishii, “Foundation of the interfacial area transport equation and its closure relations,”
*International Journal of Heat and Mass Transfer*, vol. 38, no. 3, pp. 481–493, 1995. View at Google Scholar · View at Scopus - Q. Wu, S. Kim, M. Ishii, and S. G. Beus, “One-group interfacial area transport in vertical bubbly flow,”
*International Journal of Heat and Mass Transfer*, vol. 41, no. 8-9, pp. 1103–1112, 1998. View at Google Scholar · View at Scopus - X. Fu,
*Interfacial area measurement and transport modeling in air-water two-phase flow*, Ph.D. thesis, School of Nuclear Engineering, Purdue University, West Lafayette, Ind, USA, 2001. - M. Ishii and S. Kim, “Development of one-group and two-group interfacial area transport equation,”
*Nuclear Science and Engineering*, vol. 146, no. 3, pp. 257–273, 2004. View at Google Scholar · View at Scopus - X. Sun, S. Kim, M. Ishii, and S. G. Beus, “Modeling of bubble coalescence and disintegration in confined upward two-phase flow,”
*Nuclear Engineering and Design*, vol. 230, no. 1–3, pp. 3–26, 2004. View at Publisher · View at Google Scholar · View at Scopus - M. Ishii and X. Sun, “Interfacial characteristics of two-phase flow,”
*Multiphase Science and Technology*, vol. 18, no. 1, pp. 1–29, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Vasavada,
*Interfacial area transport equation for reduced-gravity two-phase flows*, Ph.D. thesis, School of Nuclear Engineering, Purdue University, West Lafayette, Ind, USA, 2008. - D. Gidaspow, “Modeling of two phase flow,” in
*Proceedings of the 5th International Heat Transfer Conference*, vol. 7, pp. 163–168, Tokyo, Japan, 1974. - A. V. Jones and A. Prosperetti, “On the suitability of first-order differential models for two-phase flow prediction,”
*International Journal of Multiphase Flow*, vol. 11, no. 2, pp. 133–148, 1985. View at Google Scholar · View at Scopus - H. Pokharna, M. Mori, and V. H. Ransom, “Regularization of two-phase flow models: a comparison of numerical and differential approaches,”
*Journal of Computational Physics*, vol. 134, no. 2, pp. 282–295, 1997. View at Publisher · View at Google Scholar · View at Scopus - T. N. Dinh, R. R. Nourgaliev, and T. G. Theofanous, “Understanding the ill-posed two-fluid model,” in
*Proceedings of the 10th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics (NURETH '03)*, Seoul, Republic of Korea, 2003. - W. T. Hancox, R. L. Ferch, W. S. Liu, and R. E. Nieman, “One-dimensional models for transient gas-liquid flows in ducts,”
*International Journal of Multiphase Flow*, vol. 6, no. 1-2, pp. 25–40, 1980. View at Google Scholar · View at Scopus - R. T. Lahey Jr., L. Y. Cheng, D. A. Drew, and J. E. Flaherty, “The effect of virtual mass on the numerical stability of accelerating two-phase flows,”
*International Journal of Multiphase Flow*, vol. 6, no. 4, pp. 281–294, 1980. View at Publisher · View at Google Scholar · View at Scopus - N. Brauner and D. M. Maron, “Stability analysis of stratfied liquid-liquid flow,”
*International Journal of Multiphase Flow*, vol. 18, no. 1, pp. 103–121, 1992. View at Google Scholar · View at Scopus - M.-S. Chung, “Characteristic development of hyperbolic two-dimensional two-fluid model for gas-liquid flows with surface tension,”
*Applied Mathematical Modelling*, vol. 31, no. 3, pp. 578–588, 2007. View at Publisher · View at Google Scholar · View at Scopus - F. H. Harlow and D. Besnard, “Well-posed two-phase flow equations with turbulence transport,”
*Letters in Mathematical Physics*, vol. 10, no. 2-3, pp. 161–166, 1985. View at Publisher · View at Google Scholar · View at Scopus - A. L. Zanotti, C. G. Méndez, N. M. Nigro, and M. Storti, “A preconditioning mass matrix to avoid the III-posed two-fluid model,”
*Journal of Applied Mechanics, Transactions ASME*, vol. 74, no. 4, pp. 732–740, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. H. Song and M. Ishii, “The well-posedness of incompressible one-dimensional two-fluid model,”
*International Journal of Heat and Mass Transfer*, vol. 43, no. 12, pp. 2221–2231, 2000. View at Publisher · View at Google Scholar · View at Scopus - M. Ishii and N. Zuber, “Drag coefficient and relative velocity in bubbly, droplet or particulate flows,”
*AIChE Journal*, vol. 25, no. 5, pp. 843–855, 1979. View at Google Scholar · View at Scopus - M. Ishii and K. Mishima, “Two-fluid model and hydrodynamic constitutive relations,”
*Nuclear Engineering and Design*, vol. 82, no. 2-3, pp. 107–126, 1984. View at Google Scholar · View at Scopus - J. H. Song, “A remedy for the ill-posedness of the one-dimensional two-fluid model,”
*Nuclear Engineering and Design*, vol. 222, no. 1, pp. 40–53, 2003. View at Publisher · View at Google Scholar · View at Scopus - J. Song and M. Ishii, “One-dimensional two-fluid model with momentum flux parameters,”
*Nuclear Engineering and Design*, vol. 205, no. 1-2, pp. 145–158, 2001. View at Publisher · View at Google Scholar · View at Scopus