Skip Navigation



IMA Journal of Applied Mathematics Advance Access published online on January 7, 2008
IMA Journal of Applied Mathematics, doi:10.1093/imamat/hxm066
This Article
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Holmås, H.
Right arrow Articles by Schulkes, R.
Right arrow Search for Related Content
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Analysis of a 1D incompressible two-fluid model including artificial diffusion

H. Holmås{dagger}

Scandpower Petroleum Technology, Gåsevikvn. 2-4, 2027 Kjeller, Norway

T. Sira

Institute for Energy Technology, PO Box 40, 2027 Kjeller, Norway

M. Nordsveen

Scandpower Petroleum Technology, Gåsevikvn. 2-4, 2027 Kjeller, Norway

H. P. Langtangen

Simula Research Laboratory, PO Box 134, 1325 Lysaker, Norway and Department of Informatics, University of Oslo, PO Box 1080 Blindern, 0316 Oslo, Norway

R. Schulkes

Hydro Oil and Energy Research Centre, PO Box 2561, 3907 Porsgrunn, Norway

{dagger} Email: havard.holmas{at}sptgroup.com

Received on June 25, 2007; Accepted on October 16, 2007

This article examines a 1D incompressible two-fluid model including artificial tensor diffusion. The aim is to obtain a formulation that provides convergent numerical solutions for all flow conditions within the stratified and the stratified wavy flow regime. With appropriate simplifications, the two-fluid model reduces to one momentum balance, one mass conservation and two algebraic equations. It has previously been established that a formulation that is well posed in possessing exclusively real characteristics can be obtained by including an axial diffusion term in the momentum balance. In this article, however, we demonstrate that this is not sufficient to obtain a system suitable for numerical simulations. Although the unbounded growth rates of the standard two-fluid model are eliminated, linear stability theory predicts that infinitesimal wavelengths still experience finite growth. This entails that grid refinement always will result in new unstable wavelengths being resolved. On the other hand, if artificial axial diffusion is added to both the mass and the momentum equations as suggested here, a cut-off wavelength is established below which all wavelengths are stable. Thus, a numerically converging model is formed, which retains the long-wavelength properties of the standard two-fluid model. The conclusions of the mathematical analysis are substantiated by numerical simulations of 1D gravity waves.

Keywords: multiphase pipe flow; two-fluid model; diffusion; simulation; stability; well posedness.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?




Disclaimer:
Please note that abstracts for content published before 1996 were created through digital scanning and may therefore not exactly replicate the text of the original print issues. All efforts have been made to ensure accuracy, but the Publisher will not be held responsible for any remaining inaccuracies. If you require any further clarification, please contact our Customer Services Department.