H. Watanabe, M. Suzuki, N. Ito

The International Conference on Computational & Experimental Engineering and Sciences, Vol.18, No.2, pp. 49-50, 2011, DOI:10.3970/icces.2011.018.049

Abstract Gas-liquid two-phase flow is a system containing gas and liquid phase. While the gas-liquid two-phase flow is important for applications such as power plants and pump cavitations, it is difficult to study gas-liquid multiphase flow theoretically or numerically since it is multi-scale and multi-physics system involving not only flow but also phase transitions. In order to overcome the difficulties, we study the two-phase flow with full particle simulations. In the full-particle simulation, phase boundaries create and annihilate spontaneously, and therefore, multi-physics phenomena are naturally simulated. We developed a parallel molecular dynamics (MD) simulation code which is scalable up to ten… More >

Hossein M. Shodja^{1, 2, 3}, Alireza Hashemian^{2, 4}

The International Conference on Computational & Experimental Engineering and Sciences, Vol.13, No.3, pp. 57-58, 2009, DOI:10.3970/icces.2009.013.057

Abstract Gradient reproducing kernel particle method (GRKPM) is a meshless technique which incorporates the first gradients of the function into the reproducing equation of RKPM. Therefore, in two-dimensional space GRKPM introduces three types of shape functions rather than one. The robustness of GRKPM's shape functions is established by reconstruction of a third-order polynomial. To enforce the essential boundary conditions (EBCs), GRKPM's shape functions are modified by transformation technique. By utilizing the modified shape functions, the weak form of the nonlinear evolutionary Buckley-Leverett (BL) equation is discretized in space, rendering a system of nonlinear ordinary differential equations (ODEs). Subsequently, Gear's method is… More >

T.Tsukiji¹, Y.Yamamoto²

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.3, pp. 235-242, 2005, DOI:10.3970/cmes.2005.009.235

Abstract The grid free computational model of gas-particle two-phase jet flow using a 3-D vortex method is presented. The calculated results using the present method are compared with the previous experimental and the calculated results using DNS. The interaction between the particle and gas-phase is considered using Lagrangian method. It is found that the present computational model of gas-particle two-phase jet flow using the 3-D vortex method is very useful for the prediction of the physical properties of the two-phase jet flow and for saving the computational time. More >

Armin Iske¹, Martin Käser²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.2, pp. 133-148, 2005, DOI:10.3970/cmes.2005.007.133

Abstract Petroleum reservoir modelling requires effective multiscale methods for the numerical simulation of two-phase flow in porous media. This paper proposes the application of a novel meshfree particle method to the Buckley-Leverett model. The utilized meshfree advection scheme, called AMMoC, is essentially a method of characteristics, which combines an adaptive semi-Lagrangian method with local meshfree interpolation by polyharmonic splines. The method AMMoC is applied to the five-spot problem, a well-established model problem in petroleum reservoir simulation. The numerical results and subsequent numerical comparisons with two leading commercial reservoir simulators, ECLIPSE and FrontSim, show the good performance of our meshfree advection scheme… More >

Yonggang Duan¹, Ting Lu¹, Mingqiang Wei¹, Boming Yu², Zhelun Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.109-110, No.6, pp. 481-504, 2015, DOI:10.3970/cmes.2015.109.481

Abstract Analysis of Buckley-Leverett solution in fractal porous medium does prediction of water saturation profile a favor. On the approximation that porous medium consists of a bundle of tortuous capillaries, a physical conceptual Buckley- Leverett model of transient two-phase flow in fractal porous medium is developed based on the fractal characteristics of pore size distribution. The relationship between water saturation and distance is presented according to Buckley-Leverett solution, and the proposed Buckley-Leverett expression is the function of fractal structural parameters (such as pore fractal dimension, tortuosity fractal dimension, maximum and minimum diameters of capillaries) and fluid properties (such as viscosity, contact… More >

C.T. Lee¹, C.C. Lee²

CMES-Computer Modeling in Engineering & Sciences, Vol.104, No.5, pp. 375-403, 2015, DOI:10.3970/cmes.2015.104.375

Abstract In this article, we demonstrate a numerical 3-D chip, and studied the capillary dynamics inside the microchannel. We applied the level set method on the Navier-Stokes equation which incorporates the surface tension and two-phase flow characteristics. We analyzed the capillary dynamics near the junction of two microchannels. Such a highlighting point is important that it not only can provide the information of interface behavior when fluids are made into a head-on collision, but also emphasize the idea for the design of the chip. In addition, we study the pressure distribution of the fluids at the junction. It is shown that… More >

Ashraf Balabel¹,

CMES-Computer Modeling in Engineering & Sciences, Vol.89, No.4, pp. 283-301, 2012, DOI:10.3970/cmes.2012.089.283

Abstract This paper presents a novel numerical method for solving the twophase flow problems with moving interfaces in either laminar or turbulent flow regimes. The developed numerical method is based on the solution of the Reynolds- Averaged Navier Stokes equations in both phases separately with appropriate boundary conditions located at the interface separating the two fluids. The solution algorithm is performed on a regular and structured two-dimensional computational grid using the control volume approach. The complex shapes as well as the geometrical quantities of the interface are determined via the level set method. The numerical method is firstly validated against the… More >

Yongsheng An¹, Xiaodong Wu¹, Deli Gao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.89, No.2, pp. 123-142, 2012, DOI:10.3970/cmes.2012.089.123

Abstract The accuracy of numerical simulation of a two-phase (oil and water) flow in a fractured horizontal well depends greatly upon the types of grids used in the computation. Cartesian grids have been widely used in recent years, but they have some disadvantages in describing complex structural wells, such as fractured horizontal wells. For example, Cartesian grids are not efficient in describing the main wellbores and the fractures of fractured horizontal wells, and the results can frequently suffer from grid orientation effects, even though a grid-refinement is often introduced to enhance the adaptability of a Cartesian grid. The PEBI (Perpendicular Bisector)… More >

Qin Lou¹, Zhaoli Guo^1,2, Chuguang Zheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.3&4, pp. 175-188, 2011, DOI:10.3970/cmes.2011.076.175

Abstract Due to the mesoscopic and kinetic nature, the lattice Boltzmann equation (LBE) method has become an efficient and powerful tool for modeling and simulating interfacial dynamics of multi-phase flows. In this work we discuss several fundamental properties of two-phase LBE models. Particularly, the effects of force discretization, spurious currents in the vicinity of interfaces, and checkerboard effects with the underlying lattices, are investigated. More >

K. Ibrahim¹, W.A. El-Askary^1,2, A. Balabel¹, I.M. Sakr¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 79-104, 2012, DOI:10.3970/cmes.2012.085.079

Abstract In the present paper, a numerical code has been developed with different turbulence models aiming at simulating turbulent bubbly flows in vertical circular pipes. The mass and momentum conservation equations are used to describe the motion of both phases (water/air). Because of the averaging process additional models are needed for the inter-phase momentum transfer and turbulence quantities for closure. The continuous phase (water) turbulence is represented using different turbulence models namely: two-equation k-ε, extended k-ε and shear-stress transport (SST) k-ω turbulence models which contains additional term to account for the effect of the dispersed phase (air) on the continuous phase… More >