## Search

Now showing items 1-10 of 66

JavaScript is disabled for your browser. Some features of this site may not work without it.

Author

Sun, Shuyu (66)

Salama, Amgad (28)El-Amin, Mohamed (24)Kou, Jisheng (7)Allen, Rebecca (4)View MoreDepartmentEnvironmental Science and Engineering Program (66)Physical Sciences and Engineering (PSE) Division (66)Computational Transport Phenomena Lab (65)Earth Science and Engineering Program (36)Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division (8)View MoreJournalJournal of Computational and Applied Mathematics (5)Journal of Computational Physics (5)Journal of Petroleum Science and Engineering (3)Computers & Fluids (2)Computers & Mathematics with Applications (2)View MorePublisherElsevier BV (25)Society of Petroleum Engineers (SPE) (14)Springer Nature (6)Wiley-Blackwell (4)Informa UK Limited (3)View MoreSubjectTwo-phase flow (7)Porous media (6)Finite volume method (4)Porous medium (4)Finite element method (3)View MoreTypeArticle (49)Conference Paper (17)Year (Issue Date)2017 (1)2015 (6)2014 (17)2013 (13)2012 (17)View MoreItem AvailabilityMetadata Only (63)Open Access (3)

Now showing items 1-10 of 66

- List view
- Grid view
- Sort Options:
- Relevance
- Title Asc
- Title Desc
- Issue Date Asc
- Issue Date Desc
- Submit Date Asc
- Submit Date Desc
- Results Per Page:
- 5
- 10
- 20
- 40
- 60
- 80
- 100

Computing and Comparing Effective Properties for Flow and Transport in Computer-Generated Porous Media

Allen, Rebecca; Sun, Shuyu (Geofluids, Hindawi Publishing Corporation, 2017-02-13) [Article]

We compute effective properties (i.e., permeability, hydraulic tortuosity, and diffusive tortuosity) of three different digital porous media samples, including in-line array of uniform shapes, staggered-array of squares, and randomly distributed squares. The permeability and hydraulic tortuosity are computed by solving a set of rescaled Stokes equations obtained by homogenization, and the diffusive tortuosity is computed by solving a homogenization problem given for the effective diffusion coefficient that is inversely related to diffusive tortuosity. We find that hydraulic and diffusive tortuosity can be quantitatively different by up to a factor of ten in the same pore geometry, which indicates that these tortuosity terms cannot be used interchangeably. We also find that when a pore geometry is characterized by an anisotropic permeability, the diffusive tortuosity (and correspondingly the effective diffusion coefficient) can also be anisotropic. This finding has important implications for reservoir-scale modeling of flow and transport, as it is more realistic to account for the anisotropy of both the permeability and the effective diffusion coefficient.

Efficient numerical methods for simulating surface tension of multi-component mixtures with the gradient theory of fluid interfaces

Kou, Jisheng; Sun, Shuyu; Wang, Xiuhua (Computer Methods in Applied Mechanics and Engineering, Elsevier BV, 2015-08) [Article]

Surface tension significantly impacts subsurface flow and transport, and it is the main cause of capillary effect, a major immiscible two-phase flow mechanism for systems with a strong wettability preference. In this paper, we consider the numerical simulation of the surface tension of multi-component mixtures with the gradient theory of fluid interfaces. Major numerical challenges include that the system of the Euler-Lagrange equations is solved on the infinite interval and the coefficient matrix is not positive definite. We construct a linear transformation to reduce the Euler-Lagrange equations, and naturally introduce a path function, which is proven to be a monotonic function of the spatial coordinate variable. By using the linear transformation and the path function, we overcome the above difficulties and develop the efficient methods for calculating the interface and its interior compositions. Moreover, the computation of the surface tension is also simplified. The proposed methods do not need to solve the differential equation system, and they are easy to be implemented in practical applications. Numerical examples are tested to verify the efficiency of the proposed methods. © 2014 Elsevier B.V.

Parallel simulation of wormhole propagation with the Darcy-Brinkman-Forchheimer framework

Wu, Yuanqing; Salama, Amgad; Sun, Shuyu (Computers and Geotechnics, Elsevier BV, 2015-07-09) [Article]

The acid treatment of carbonate reservoirs is a widely practiced oil and gas well stimulation technique. The injected acid dissolves the material near the wellbore and creates flow channels that establish a good connectivity between the reservoir and the well. Such flow channels are called wormholes. Different from the traditional simulation technology relying on Darcy framework, the new Darcy-Brinkman-Forchheimer (DBF) framework is introduced to simulate the wormhole forming procedure. The DBF framework considers both large and small porosity conditions and should output better simulation results than the Darcy framework. To process the huge quantity of cells in the simulation grid and shorten the long simulation time of the traditional serial code, a parallel code with FORTRAN 90 and MPI was developed. The experimenting field approach to set coefficients in the model equations was also introduced. Moreover, a procedure to fill in the coefficient matrix in the linear system in the solver was described. After this, 2D dissolution experiments were carried out. In the experiments, different configurations of wormholes and a series of properties simulated by both frameworks were compared. We conclude that the numerical results of the DBF framework are more like wormholes and more stable than the Darcy framework, which is a demonstration of the advantages of the DBF framework. Finally, the scalability of the parallel code was evaluated, and we conclude that superlinear scalability can be achieved. © 2015 Elsevier Ltd.

Numerical and dimensional analysis of nanoparticles transport with two-phase flow in porous media

El-Amin, Mohamed; Salama, Amgad; Sun, Shuyu (Journal of Petroleum Science and Engineering, Elsevier BV, 2015-04) [Article]

In this paper, a mathematical model and numerical simulation are developed to describe the imbibition of nanoparticles-water suspension into two-phase flow in a porous medium. The flow system may be changed from oil-wet to water-wet due to nanoparticles (which are also water-wet) deposition on surface of the pores. So, the model is extended to include the negative capillary pressure and mixed-wet relative permeability correlations to fit with the mixed-wet system. Moreover, buoyancy and capillary forces as well as Brownian diffusion and mechanical dispersion are considered in the mathematical model. An example of countercurrent imbibition in a core of small scale is considered. A dimensional analysis of the governing equations is introduced to examine contributions of each term of the model. Several important dimensionless numbers appear in the dimensionless equations, such as Darcy number Da, capillary number Ca, and Bond number Bo. Throughout this investigation, we monitor the changing of the fluids and solid properties due to addition of the nanoparticles using numerical experiments.

Unconditionally stable methods for simulating multi-component two-phase interface models with Peng-Robinson equation of state and various boundary conditions

Kou, Jisheng; Sun, Shuyu (Journal of Computational and Applied Mathematics, Elsevier BV, 2015-03) [Article]

In this paper, we consider multi-component dynamic two-phase interface models, which are formulated by the Cahn-Hilliard system with Peng-Robinson equation of state and various boundary conditions. These models can be derived from the minimum problems of Helmholtz free energy or grand potential in the realistic thermodynamic systems. The resulted Cahn-Hilliard systems with various boundary conditions are fully coupled and strongly nonlinear. A linear transformation is introduced to decouple the relations between different components, and as a result, the models are simplified. From this, we further propose a semi-implicit unconditionally stable time discretization scheme, which allows us to solve the Cahn-Hilliard system by a decoupled way, and thus, our method can significantly reduce the computational cost and memory requirements. The mixed finite element methods are employed for the spatial discretization, and the approximate errors are also analyzed for both space and time. Numerical examples are tested to demonstrate the efficiency of our proposed methods. © 2015 Elsevier B.V.

Speeding up the flash calculations in two-phase compositional flow simulations - The application of sparse grids

Wu, Yuanqing; Kowitz, Christoph; Sun, Shuyu; Salama, Amgad (Journal of Computational Physics, Elsevier BV, 2015-03) [Article]

Flash calculations have become a performance bottleneck in the simulation of compositional flow in subsurface reservoirs. We apply a sparse grid surrogate model to substitute the flash calculation and thus try to remove the bottleneck from the reservoir simulation. So instead of doing a flash calculation in each time step of the simulation, we just generate a sparse grid approximation of all possible results of the flash calculation before the reservoir simulation. Then we evaluate the constructed surrogate model to approximate the values of the flash calculation results from this surrogate during the simulations. The execution of the true flash calculation has been shifted from the online phase during the simulation to the offline phase before the simulation. Sparse grids are known to require only few unknowns in order to obtain good approximation qualities. In conjunction with local adaptivity, sparse grids ensure that the accuracy of the surrogate is acceptable while keeping the memory usage small by only storing a minimal amount of values for the surrogate. The accuracy of the sparse grid surrogate during the reservoir simulation is compared to the accuracy of using a surrogate based on regular Cartesian grids and the original flash calculation. The surrogate model improves the speed of the flash calculations and the simulation of the whole reservoir. In an experiment, it is shown that the speed of the online flash calculations is increased by about 2000 times and as a result the speed of the reservoir simulations has been enhanced by 21 times in the best conditions.

Investigation of thermal energy transport from an anisotropic central heating element to the adjacent channels: A multipoint flux approximation

Salama, Amgad; Sun, Shuyu; El-Amin, Mohamed (Annals of Nuclear Energy, Elsevier BV, 2015-02) [Article]

The problem of heat transfer from a central heating element pressed between two clad plates to cooling channels adjacent and outboard of the plates is investigated numerically. The aim of this work is to highlight the role of thermal conductivity anisotropy of the heating element and/or the encompassing plates on thermal energy transport to the fluid passing through the two channels. When the medium is anisotropic with respect to thermal conductivity; energy transport to the neighboring channels is no longer symmetric. This asymmetry in energy fluxes influence heat transfer to the coolant resulting in different patterns of temperature fields. In particular, it is found that the temperature fields are skewed towards the principal direction of anisotropy. In addition, the heat flux distributions along the edges of the heating element are also different as a manifestation of thermal conductivity anisotropy. Furthermore, the peak temperature at the channel walls change location and magnitude depending on the principal direction of anisotropy. Based on scaling arguments, it is found that, the ratio of width to the height of the heating system is a key parameter which can suggest when one may ignore the effect of the cross-diagonal terms of the full conductivity tensor. To account for anisotropy in thermal conductivity, the method of multipoint flux approximation (MPFA) is employed. Using this technique, it is possible to find a finite difference stencil which can handle full thermal conductivity tensor and in the same time enjoys the simplicity of finite difference approximation. Although the finite difference stencil based on MPFA is quite complex, in this work we apply the recently introduced experimenting field approach which construct the global problem automatically.

Quasistatic analysis on configuration of two-phase flow in Y-shaped tubes

Zhong, Hua; Wang, Xiaoping; Salama, Amgad; Sun, Shuyu (Computers & Mathematics with Applications, Elsevier BV, 2014-12) [Article]

We investigate the two-phase flow in a horizontally placed Y-shaped tube with different Young's angle and width in each branch. By using a quasistatic approach, we can determine the specific contact position and the equilibrium contact angle of fluid in each branch based on the minimization problem of the free energy of the system. The wettability condition and the width of the two branches play important roles in the distribution of fluid in each branch. We also consider the effect of gravity. Some fluid in the upper branch will be pulled down due to the competition of the surface energy and the gravitational energy. The result provides some insights on the theory of two-phase flow in porous media. In particular, it highlights that the inhomogeneous wettability distribution affects the direction of the fluid penetrating a given porous medium domain. It also sheds light on the current debate whether relative permeability may be considered as a full tensor rather than a scalar.

Numerical investigation of natural convection in two enclosures separated by anisotropic solid wall

Salama, Amgad; El-Amin, Mohamed; Sun, Shuyu (International Journal of Numerical Methods for Heat & Fluid Flow, Emerald, 2014-10-28) [Article]

Purpose: The problem of natural convection in two cavities separated by an anisotropic central solid wall is considered numerically. When the thermal conductivity of the central wall is anisotropic, heat flux and temperature gradient vectors are no longer coincidence. This apparently has interesting influences on the heat and fluid flow patterns in this system. The paper aims to discuss these issues.Design/methodology/approach: In this work, several flow patterns have been investigated covering a wide range of Rayleigh number up to 108. Several thermal conductivity anisotropy scenarios of the central wall have been investigated including 0, 30, 60, 120 and 1501 principal anisotropy directions. The governing equations have been solved using control volume approach.Findings: Probably the most intriguing is that, for some anisotropy scenarios it is found that the temperature at the same elevation at the side of the central wall which is closer to the colder wall is higher than that at the side closer to the hot wall. Apparently this defies intuition which suggests the reverse to have happened. However, this behavior may be explained in light of the effect of anisotropy. Furthermore, the patterns of streamlines and temperature fields in the two enclosures also changes as a consequence of the change of the central wall temperatures for the different anisotropy scenarios.Originality/value: This work discusses a very interesting topic related to heat energy exchange among two compartments when the separating wall is anisotropic. In some anisotropy scenarios, this leads to more uniform distribution of Nusselt number than the case when the wall is isotropic. Interesting patterns of natural convection is investigated.

EMSE: Synergizing EM and seismic data attributes for enhanced forecasts of reservoirs

Katterbauer, Klemens; Hoteit, Ibrahim; Sun, Shuyu (Journal of Petroleum Science and Engineering, Elsevier BV, 2014-10) [Article]

New developments of electromagnetic and seismic techniques have recently revolutionized the oil and gas industry. Time-lapse seismic data is providing engineers with tools to more accurately track the dynamics of multi-phase reservoir fluid flows. With the challenges faced in distinguishing between hydrocarbons and water via seismic methods, the industry has been looking at electromagnetic techniques in order to exploit the strong contrast in conductivity between hydrocarbons and water. Incorporating this information into reservoir simulation is expected to considerably enhance the forecasting of the reservoir, hence optimizing production and reducing costs. Conventional approaches typically invert the seismic and electromagnetic data in order to transform them into production parameters, before incorporating them as constraints in the history matching process and reservoir simulations. This makes automatization difficult and computationally expensive due to the necessity of manual processing, besides the potential artifacts. Here we introduce a new approach to incorporate seismic and electromagnetic data attributes directly into the history matching process. To avoid solving inverse problems and exploit information in the dynamics of the flow, we exploit petrophysical transformations to simultaneously incorporate time lapse seismic and electromagnetic data attributes using different ensemble Kalman-based history matching techniques. Our simulation results show enhanced predictability of the critical reservoir parameters and reduce uncertainties in model simulations, outperforming with only production data or the inclusion of either seismic or electromagnetic data. A statistical test is performed to confirm the significance of the results. © 2014 Elsevier B.V. All rights reserved.

The export option will allow you to export the current search results of the entered query to a file. Different formats are available for download. To export the items, click on the button corresponding with the preferred download format.

By default, clicking on the export buttons will result in a download of the allowed maximum amount of items. For anonymous users the allowed maximum amount is 50 search results.

To select a subset of the search results, click "Selective Export" button and make a selection of the items you want to export. The amount of items that can be exported at once is similarly restricted as the full export.

After making a selection, click one of the export format buttons. The amount of items that will be exported is indicated in the bubble next to export format.