Author:  Peng, Qiujin 
Title:  Numerical solutions of a diffusive interface model with PengRobinson equation of state 
Degree:  Ph.D. 
Year:  2015 
Subject:  Fluid dynamics  Mathematical models. Complex fluids  Analysis. Complex fluids  Mathematics. Hong Kong Polytechnic University  Dissertations 
Department:  Dept. of Applied Mathematics 
Pages:  xi, 118 pages : color illustrations 
Language:  English 
InnoPac Record:  http://library.polyu.edu.hk/record=b2837218 
URI:  http://theses.lib.polyu.edu.hk/handle/200/8403 
Abstract:  This work is concerned with mathematical modeling and numerical simulations of the steady state and the movements of complex fluids involved in oil exploitation practice. Capillary pressure caused by surface tension at the interface between every two adjacent different phases of the mixture is viewed as the leading force in oil recovery from fractured oil reservoirs. Therefore, the interface between contiguous phases has become a critical mathematical modeling aspect. The diffuse interface theory, or the phase field model, or the gradient theory, has been widely applied to model or understand the interface between different phases of oil mixture. Based on the assumption that the density of every substance is continuous over the whole fluid region, the total Helmholtz free energy often contains the homogeneous part F0 (n) and the gradient contribution part Fv (n). The derivative of the total homogeneous free energy, f0 (n), or of the gradient part of free energy, fv(n), varies from substance to substance. Based on the original total free energy, the equilibrium state and the kinetic processes could be determined according to thermodynamic principles. As for the fluid system related to the oil recovery process, we apply the homogeneous free energy density and the parameters of the gradient part of the free energy density provided by the widely used PengRobinson equation of state (EOS). The fourthorder parabolic equation is derived and solved numerically by a convexsplitting scheme, the CrankNicolson scheme and a second order linearization scheme to describe the evolution processes of onecomponent, twophase substances. The theoretical analyses of these numerical schemes have been obtained to demonstrate their mass conservation, energy stability, unique solvability and convergence.The EulerLagrange (EL) equation derived from this expression of total Helmholtz free energy to determine the equilibrium state of the fluid systems has also been studied. In this study, it is solved numerically by the original Newton iteration and a convexsplitting based Newton iterative method. Its theoretical analysis remains a part of our future work. Numerical experiments have been carried out for both the fourthorder equation approach and the EulerLagrange equation approach. Our computational results match well with laboratory experimental data and are in good agreement with the wellknown YoungLaplace equation. 
Files  Size  Format 

b28372189.pdf  860.2Kb 


As a bona fide Library user, I declare that:  


By downloading any item(s) listed above, you acknowledge that you have read and understood the copyright undertaking as stated above, and agree to be bound by all of its terms. 