The Gas Assisted Gravity Drainage (GAGD) process has become one of the most important processes to enhance oil recovery in both secondary and tertiary recovery stages and through immiscible and miscible modes. Its advantages came from the ability to provide gravity-stable oil displacement for improving oil recovery, when compared with conventional gas injection methods such as Continuous Gas Injection (CGI) and Water – Alternative Gas (WAG). Vertical injectors for CO2 gas were placed at the top of the reservoir to form a gas cap which drives the oil towards the horizontal oil producing wells which are located above the oil-water-contact. The GAGD process was developed and tested in vertical wells to increase oil recovery in reservoirs with bottom water drive and strong water coning tendencies. Many physical and simulation models of GAGD performance were studied at ambient and reservoir conditions to investigate the effects of this method to enhance the recovery of oil and to examine the most effective parameters that control the GAGD process. A prototype 2D simulation model based on the scaled physical model was built for CO2-assisted gravity drainage in different statement scenarios. The effects of gas injection rate, gas injection pressure and oil production rate on the performance of immiscible CO2-assisted gravity drainage-enhanced oil recovery were investigated. The results revealed that the ultimate oil recovery increases considerably with increasing oil production rates. Increasing gas injection rate improves the performance of the process while high pressure gas injection leads to less effective gravity mediated recovery.
This paper is concerned with the numerical solutions of the vorticity transport equation (VTE) in two-dimensional space with homogenous Dirichlet boundary conditions. Namely, for this problem, the Crank-Nicolson finite difference equation is derived. In addition, the consistency and stability of the Crank-Nicolson method are studied. Moreover, a numerical experiment is considered to study the convergence of the Crank-Nicolson scheme and to visualize the discrete graphs for the vorticity and stream functions. The analytical result shows that the proposed scheme is consistent, whereas the numerical results show that the solutions are stable with small space-steps and at any time levels.
In this study the simple pullout concrete cylinder specimen reinforced by a single steel bar was analyzed for bond-slip behavior. Three-dimension nonlinear finite element model using ANSYS program was employed to study the behavior of bond between concrete and plain steel reinforcement. The ANSYS model includes eight-noded isoperimetric brick element (SOLID65) to model the concrete cylinder while the steel reinforcing bar was modeled as a truss member (LINK8). Interface element (CONTAC52) was used in this analysis to model the bond between concrete and steel bar. Material nonlinearity due to cracking and/or crushing of concrete, and yielding of the steel reinforcing bar were taken into consideration during the analysis. The accuracy of this
... Show MoreIn this study the simple pullout concrete cylinder specimen reinforced by a single steel bar was analyzed for bond-slip behavior. Three-dimension nonlinear finite element model using ANSYS program was employed to study the behavior of bond between concrete and plain steel reinforcement. The ANSYS model includes eight-noded isoperimetric brick element (SOLID65) to model the concrete cylinder while the steel reinforcing bar was modeled as a truss member (LINK8). Interface element (CONTAC52) was used in this analysis to model the bond between concrete and steel bar. Material nonlinearity due to cracking and/or crushing of concrete, and yielding of the steel reinforcing bar were taken into consideration during the analysis. The accuracy of t
... Show MoreCurrent numerical research was devoted to investigating the effect of castellated steel beams without and with strengthening. The composite concrete asymmetrical double hot rolled steel channels bolted back to back to obtain a built-up I-shape form are used in this study. The top half part of the steel is smaller than the bottom half part, and the two parts were connected by bolting and welding. The ABAQUS/2019 program employed the same length and conditions of loading for four models: The first model is the reference without castellated and strengthening; the second model was castellated without strengthened; the third model was castellated and strengthened with reactive powder concrete encased in the
... Show MoreIn this work, we use the explicit and the implicit finite-difference methods to solve the nonlocal problem that consists of the diffusion equations together with nonlocal conditions. The nonlocal conditions for these partial differential equations are approximated by using the composite trapezoidal rule, the composite Simpson's 1/3 and 3/8 rules. Also, some numerical examples are presented to show the efficiency of these methods.
This paper reports experimental and computational fluid dynamics (CFD) modelling studies to investigate the effect of the swirl intensity on the heat transfer characteristics of conventional and swirl impingement air jets at a constant nozzle-to-plate distance ( L = 2 D). The experiments were performed using classical twisted tape inserts in a nozzle jet with three twist ratios ( y = 2.93, 3.91, and 4.89) and Reynolds numbers that varied from 4000 to 16000. The results indicate that the radial uniformity of Nusselt number (Nu) of swirl impingement air jets (SIJ) depended on the values of the swirl intensity and the air Reynolds number. The results also revealed that the SIJ that was fitted with an insert of y = 4.89, which correspo
... Show MoreNumerical simulations are carried out to assess the quality of the circular and square apodize apertures in observing extrasolar planets. The logarithmic scale of the normalized point spread function of these apertures showed sharp decline in the radial frequency components reaching to 10-36 and 10-34 respectively and demonstrating promising results. This decline is associated with an increase in the full width of the point spread function. A trade off must be done between this full width and the radial frequency components to overcome the problem of imaging extrasolar planets.
This paper is dealing with non-polynomial spline functions "generalized spline" to find the approximate solution of linear Volterra integro-differential equations of the second kind and extension of this work to solve system of linear Volterra integro-differential equations. The performance of generalized spline functions are illustrated in test examples
In this study, experimental and numerical applied of heat distribution due to pulsed Nd: YAG laser surface melting. Experimental side was consists of laser parameters are, pulse duration1.3
The aim of this paper is to evaluate the rate of contamination in soils by using accurate numerical method as a suitable tool to evaluate the concentration of heavy metals in soil. In particular, 2D –interpolation methods are applied in the models of the spread the metals in different direction.The paper illustrates the importance of the numerical method in different applications, especially nvironment contamination. Basically, there are many roles for approximating functions. Thus, the approximating of function namely the analytical expression may be expressed; the most common type being is polynomials, which are the easy implemented and simplest methods of approximation. In this paper the divided difference formula is used and extended
... Show More