Mustapha Gourma

Use of OpenMP to parallelise a one-dimensional multiphase flow code

The project of parallelising a Fortran code designed to simulate one-dimensional multiphase flow simulations has been successfully completed. Parallelisation of such a code is of crucial importance due to its use for simulating flow in long pipes. The choice was made to use OpenMP rather than MPI parallelisation, and reasons for the choice are given in detail, together with explanations about the tasks involved. Speed-ups were observed to be better for complex mathematical models rather than simple ones.

Y-dependent perturbation analysis of wave growth for 2D volume of fluid model

The wave growth analysis is one of various methodologies used to validate 1D multiphase flow codes [1], and it has also been used as a first stage validation of two-dimensional volume of fluid (VOF) CFD simulations of pipeline multiphase flow [2]. A mathematical perturbation analysis to first order was carried out and shown to be in good agreement with simulation results [2]. In this paper a more in-depth perturbation analysis is carried out, taking into account the y-dependency of perturbation parameters and re-evaluating the surface tension term.

Validation of 2D CFD for Two-Phase Transient Flow in a Channel and Comparison With 1D Model

Severe limitations of the use of three-dimensional computational fluid dynamics codes (CFD) arise when trying to simulate multiphase flow in long pipes due to time constraints. 1D codes for two-phase flow, based on two-fluid models, are fast but are known to be accurate only when the velocities are within the Kelvin-Helmholtz inviscid limit [1]. An alternative is to carry out a two-dimensional CFD simulation of a channel based on the Volume of Fluid (VOF) model. 2D CFD has a wider applicability range compared to 1D, it does not have the issue of ill-posedness and it also has better turbulence models built in. Again compared to 1D the 2D VOF model has a better interface description and wall treatment. In this paper a novel method is introduced that allows swift simulations of pipeline two-phase flow in the stratified and slug flow regime, by approximating the pipe as a channel and with a methodology that solves the problem of the interfacial velocity differences, inherent in the volume of fluid model. An initial validation using the wave growth problem has already been carried out [2]. Here a set consisting of 92 experimental cases in the slug flow regime has been simulated with 2D CFD, and the simulation results showed a good agreement with experimental results. Discussions in the paper include also the question of the range of applicability for 2D CFD, and the advantages and disadvantages compared to 3D CFD and also to 1D code based on the two-fluid model. Shear stresses are then extracted from the 2D CFD simulations and used to recalibrate the friction factors [3] used in the 1D code.Copyright