×

Comparison of Roe-type methods for solving the two-fluid model with and without pressure relaxation. (English) Zbl 1194.76161

Summary: Two strategies for the numerical resolution of a two-fluid model have been investigated. Both methods employ a Roe-type scheme. The first method (Roe4) solves the four-equation, one-pressure, isentropic two-fluid model directly. The second strategy (Roe5) is to add an evolution equation for the volume fraction. In the present case, that results in a five-equation two-pressure model, where it is necessary to employ pressure relaxation to calculate the typical two-phase problems that have been tested: The water faucet and two benchmark shock-tube problems known from the literature. The numerical calculations showed that the Roe4 and Roe5 schemes converge to the same solution when instantaneous pressure relaxation is employed in the Roe5 scheme. This is true both with and without the use of high-resolution flux-limiter functions. However, the Roe5 scheme was found to be significantly more diffusive than the Roe4 scheme. The diffusion is a strong function of the chosen time-step length, the grid size, whether a limiter function is employed or not, and also the liquid speed of sound. As the pressure-relaxation parameter in the Roe5 scheme was increased, the solution gradually approached that obtained using instantaneous pressure relaxation. Furthermore, the results indicate that the approach of two pressures and instantaneous pressure relaxation does not provide an easy way to overcome the problem of complex eigenvalues in the one-pressure two-fluid model.

MSC:

76M12 Finite volume methods applied to problems in fluid mechanics
76N15 Gas dynamics (general theory)

Software:

CATHARE
Full Text: DOI

References:

[1] Abgrall, R., How to prevent pressure oscillations in multicomponent flow calculations: a quasi conservative approach, J Comput Phys, 125, 1, 150-160 (1996) · Zbl 0847.76060
[2] Abgrall, R.; Saurel, R., Discrete equations for physical and numerical compressible multiphase mixtures, J Comput Phys, 186, 2, 361-396 (2003) · Zbl 1072.76594
[3] Baer, M. R.; Nunziato, J. W., A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials, Int J Multiphase Flow, 12, 6, 861-889 (1986) · Zbl 0609.76114
[4] Bestion, D., The physical closure laws in the CATHARE code, Nucl Eng Des, 124, 3, 229-245 (1990)
[5] Chung, M. S.; Chang, K. S.; Lee, S. J., Numerical solution of hyperbolic two-fluid two-phase flow model with non-reflecting boundary conditions, Int J Eng Sci, 40, 7, 789-803 (2002) · Zbl 1211.76078
[6] Coquel, F., A numerical method using upwind schemes for the resolution of two-phase flows, J Comput Phys, 136, 2, 272-288 (1997) · Zbl 0893.76052
[7] Cortes, J.; Debussche, A.; Toumi, I., A density perturbation method to study the eigenstructure of two-phase flow equation systems, J Comput Phys, 147, 2, 463-484 (1998) · Zbl 0917.76047
[8] Drew, D. A.; Passman, S. L., Theory of multicomponent fluids, Applied mathematical sciences, vol. 135 (1999), Springer-Verlag: Springer-Verlag New York
[9] Evje, S.; Flåtten, T., Hybrid flux-splitting schemes for a common two-fluid model, J Comput Phys, 192, 1, 175-210 (2003) · Zbl 1032.76696
[10] Evje, S.; Flåtten, T., Hybrid central-upwind schemes for numerical resolution of two-phase flows, ESAIM - Math Model Num, 29, 2, 253-273 (2005) · Zbl 1130.76057
[11] Flåtten T. Hybrid flux splitting schemes for numerical resolution of two-phase flows. Dr. ing. thesis, Norwegian University of Science and Technology, Department of Energy and Process Engineering, Trondheim, 2003. ISBN 82-471-5670-9.; Flåtten T. Hybrid flux splitting schemes for numerical resolution of two-phase flows. Dr. ing. thesis, Norwegian University of Science and Technology, Department of Energy and Process Engineering, Trondheim, 2003. ISBN 82-471-5670-9.
[12] Harten, A., High resolution schemes for hyperbolic conservation laws, J Comput Phys, 49, 3, 357-393 (1983) · Zbl 0565.65050
[13] Karni, S., Compressible two-phase flows by central and upwind schemes, ESAIM - Math Model Num, 38, 3, 477-493 (2004) · Zbl 1079.76045
[14] Lallemand, M. H.; Chinnayya, A.; Le Metayer, O., Pressure relaxation procedures for multiphase compressible flows, Int J Numer Meth Fl, 49, 1, 1-56 (2005) · Zbl 1073.76056
[15] LeVeque, R. J., Finite volume methods for hyperbolic problems (2002), Cambridge University Press: Cambridge University Press Cambridge (UK) · Zbl 1010.65040
[16] Munkejord ST. Analysis of the two-fluid model and the drift-flux model for numerical calculation of two-phase flow. Doctoral thesis, Norwegian University of Science and Technology, Department of Energy and Process Engineering, Trondheim, 2005. ISBN 82-471-7338-7.; Munkejord ST. Analysis of the two-fluid model and the drift-flux model for numerical calculation of two-phase flow. Doctoral thesis, Norwegian University of Science and Technology, Department of Energy and Process Engineering, Trondheim, 2005. ISBN 82-471-7338-7.
[17] Munkejord, S. T., Partially-reflecting boundary conditions for transient two-phase flow, Commun Numer Meth En, 22, 7, 781-795 (2006) · Zbl 1115.76053
[18] Munkejord, S. T.; Papin, M., The effect of interfacial pressure in the discrete-equation multiphase model, Comput Fluids, 36, 4, 742-757 (2007) · Zbl 1177.76448
[19] Murrone, A.; Guillard, H., A five equation reduced model for compressible two phase flow problems, J Comput Phys, 202, 2, 664-698 (2005) · Zbl 1061.76083
[20] Niu, Y. Y., Advection upwinding splitting method to solve a compressible two-fluid model, Int J Numer Meth Fl, 36, 3, 351-371 (2001) · Zbl 1044.76041
[21] Paillère, H.; Corre, C.; García Cascales, J. R., On the extension of the AUSM+ scheme to compressible two-fluid models, Comput Fluids, 32, 6, 891-916 (2003) · Zbl 1040.76044
[22] Ramshaw, J. D.; Trapp, J. A., Characteristics, stability, and short-wavelength phenomena in two-phase flow equation systems, Nucl Sci Eng, 66, 1, 93-102 (1978)
[23] Ransom, V. H., Faucet flow, (Hewitt, G. F.; Delhaye, J. M.; Zuber, N., Numerical benchmark tests. Numerical benchmark tests, Multiphase science and technology, vol. 3 (1987), Hemisphere/Springer: Hemisphere/Springer Washington (USA)), 465-467
[24] Ransom, V. H.; Hicks, D. L., Hyperbolic two-pressure models for two-phase flow, J Comput Phys, 53, 1, 124-151 (1984) · Zbl 0537.76070
[25] Roe, P. L., Approximate Riemann solvers, parameter vectors, and difference schemes, J Comput Phys, 43, 2, 357-372 (1981) · Zbl 0474.65066
[26] Sainsaulieu, L., Finite-volume approximation of two phase-fluid flows based on an approximate Roe-type Riemann solver, J Comput Phys, 121, 1, 1-28 (1995) · Zbl 0834.76070
[27] Saurel, R.; Abgrall, R., A multiphase Godunov method for compressible multifluid and multiphase flow, J Comput Phys, 150, 2, 425-467 (1999) · Zbl 0937.76053
[28] Saurel, R.; LeMetayer, O., A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation, J Fluid Mech, 431, 239-271 (2001) · Zbl 1039.76069
[29] Soo, S. L., Multiphase fluid dynamics (1990), Science Press: Science Press Beijing · Zbl 0173.52901
[30] Tiselj, I.; Petelin, S., Modelling of two-phase flow with second-order accurate scheme, J Comput Phys, 136, 2, 503-521 (1997) · Zbl 0918.76050
[31] Toro, E. F., Riemann solvers and numerical methods for fluid dynamics (1999), Springer-Verlag: Springer-Verlag Berlin · Zbl 0923.76004
[32] Toumi, I., An upwind numerical method for two-fluid two-phase flow models, Nucl Sci Eng, 123, 2, 147-168 (1996)
[33] Toumi, I.; Kumbaro, A., An approximate linearized Riemann solver for a two-fluid model, J Comput Phys, 124, 2, 286-300 (1996) · Zbl 0847.76056
[34] Trapp, J. A.; Riemke, R. A., A nearly-implicit hydrodynamic numerical scheme for two-phase flows, J Comput Phys, 66, 1, 62-82 (1986) · Zbl 0622.76110
[35] van Leer, B., Towards the ultimate conservative difference scheme IV. New approach to numerical convection, J Comput Phys, 23, 3, 276-299 (1977) · Zbl 0339.76056
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.