Change search

Numerical Solvers for Transient Two-Phase Flow
Norwegian University of Science and Technology, Faculty of Engineering Science and Technology, Department of Energy and Process Engineering.
2013 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

Certain numerical methods have been well developed for solving one-dimensional two-phase flow (e.g. gas and liquid) problems in the literatures during the last two decades. Based on the existing methods, the present work compares the computational efficiency, accuracy, and robustness of various numerical schemes by predicting the numerical solutions of fluid properties for a specific case to find the proper numerical method.

One of the numerical schemes introduced in this work is a practical, semi-implicit upwind method used for fluid flow simulations in different flow patterns,stratified flow and slug flow. This method implements the iterative and non-iterative schemes using a two-fluid model that consists of sets of non-hyperbolic equations. A numerical error term is applied in the pressure equation to maintain the volume balance of the two-phase flow model. If the temperature varies, the discretised energy equations use similar error terms as in the pressure equation. In some cases, the small values of the numerical errors are negligible and do not influence the numerical results. These errors are, however, important factors to consider when maintaining the stability and robustness of the above numerical schemes for strong non-linear cases. The computational efficiency ofthe non-iterative scheme, where the inner iterations are deactivated, is better than the iterative scheme.

Different grid arrangements are compared with respect to computational accuracy and efficiency. A staggered structured grid implements the same semi-implicit upwind method as in the non-iterative scheme; the non-staggered grid arrangement uses an existing flux-splitting scheme (Evje and Flåtten, 2003) as a reference.

All the above schemes produce numerical solutions with a single precision that normally satisfy the requirements of computational accuracy of industrial two-phase pipe flows. However, if one pursues a higher-order accuracy scheme, e.g. a Roe-averaged algorithm, the governing equations should be strictly a hyperbolic system of partial differential equations, which is achieved by introducing the nonviscous force terms in the two-fluid model (LeVeque, 2002).By properly incorporating the non-conservative terms in the formulation of the numerical fluxes, the capability of the Roe-averaged algorithm is demonstrated by capturing shock waves.

Results from the present research include the following. A one-dimensional scheme that solves a system of discretised equations with the staggered semi-implicit upwind method is presented and validated for its computational efficiencyand robustness. This scheme can be widely used in the industry with sufficient accuracy. The other first-order semi-implicit numerical schemes producestable numerical results, especially in the dynamic cases of two-phase flow, except when the gas phase nearly disappears or appears in pipes. The Roe-averaged algorithm is recommended due to the high-resolution numerical results obtained, but at the costs of computational time and effort.

Place, publisher, year, edition, pages
Trondheim NTNU: Skipnes Kommunikasjon as , 2013. , 148 p.
Series
Doctoral Theses at NTNU, ISSN 1503-8181 ; 2013:202
Keyword [en]
Numerical schemes, two-phase flow, one dimension
National Category
Mathematical modeling and numerical methods
Identifiers
ISBN: 978-82-471-4519-7ISBN: 978-82-471-4519-7 (electronic version)OAI: oai:DiVA.org:ntnu-23608DiVA: diva2:677279
Public defence
2013-12-03, EL2, Trondheim, 13:47 (English)
Supervisors
Available from: 2013-12-17 Created: 2013-12-09 Last updated: 2014-01-08Bibliographically approved
List of papers
1. Iterative versus Non-iterative Numerical Schemes forOne-Dimensional Compressible Transient Flow PartI: Isothermal Flow
Open this publication in new window or tab >>Iterative versus Non-iterative Numerical Schemes forOne-Dimensional Compressible Transient Flow PartI: Isothermal Flow
Identifiers
urn:nbn:no:ntnu:diva-23642 (URN)
Available from: 2013-12-17 Created: 2013-12-17 Last updated: 2013-12-17Bibliographically approved
2. Iterative versus Non-iterative Numerical Schemes forOne-Dimensional Compressible Transient Flow PartII: Thermal Flow
Open this publication in new window or tab >>Iterative versus Non-iterative Numerical Schemes forOne-Dimensional Compressible Transient Flow PartII: Thermal Flow
Identifiers
urn:nbn:no:ntnu:diva-23643 (URN)
Available from: 2013-12-17 Created: 2013-12-17 Last updated: 2013-12-17Bibliographically approved
3. Staggered and Non-staggered Numerical Schemes for Two-Phase Transient Flow
Open this publication in new window or tab >>Staggered and Non-staggered Numerical Schemes for Two-Phase Transient Flow
Identifiers
urn:nbn:no:ntnu:diva-23644 (URN)
Available from: 2013-12-17 Created: 2013-12-17 Last updated: 2013-12-17Bibliographically approved

Open Access in DiVA

File information
File name FULLTEXT03.pdfFile size 1399 kBChecksum SHA-512
ea7765daba4518a2002ccd4e3595247b3ca9a69bab3ab7ddf33e10dec764fcdc1203e5005a9f66482a9b039f3848015fd4bba4194b91f535c384654034b4a873
Type fulltextMimetype application/pdf
By organisation
Department of Energy and Process Engineering
On the subject
Mathematical modeling and numerical methods