Recent Submissions

  • Numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow

    Zhu, Guangpu; Kou, Jisheng; Sun, Shuyu; Yao, Jun; Li, Aifen (arXiv, 2018-04-17) [Preprint]
    In this paper, we consider the numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow. The nonlinearly coupled model consists of two Cahn-Hilliard type equations and incompressible Navier-Stokes equations. Using the Invariant Energy Quadratization (IEQ) approach, the governing system is transformed into an equivalent form, which allows the nonlinear potentials to be treated efficiently and semi-explicitly. we construct a first and a second-order time marching schemes, which are extremely efficient and easy-to-implement, for the transformed governing system. At each time step, the schemes involve solving a sequence of linear elliptic equations, and computations of phase variables, velocity and pressure are totally decoupled. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient and unconditionally energy stable. Using our schemes, we investigate the effect of surfactants on droplet deformation and collision under a shear flow. The increase of surfactant concentration can enhance droplet deformation and inhibit droplet coalescence.
  • Global mass conservation method for dual-continuum gas reservoir simulation

    Wang, Yi; Sun, Shuyu; Gong, Liang; Yu, Bo (Journal of Natural Gas Science and Engineering, Elsevier BV, 2018-03-17) [Article]
    In this paper, we find that the numerical simulation of gas flow in dual-continuum porous media may generate unphysical or non-robust results using regular finite difference method. The reason is the unphysical mass loss caused by the gas compressibility and the non-diagonal dominance of the discretized equations caused by the non-linear well term. The well term contains the product of density and pressure. For oil flow, density is independent of pressure so that the well term is linear. For gas flow, density is related to pressure by the gas law so that the well term is non-linear. To avoid these two problems, numerical methods are proposed using the mass balance relation and the local linearization of the non-linear source term to ensure the global mass conservation and the diagonal dominance of discretized equations in the computation. The proposed numerical methods are successfully applied to dual-continuum gas reservoir simulation. Mass conservation is satisfied while the computation becomes robust. Numerical results show that the location of the production well relative to the large-permeability region is very sensitive to the production efficiency. It decreases apparently when the production well is moved from the large-permeability region to the small-permeability region, even though the well is very close to the interface of the two regions. The production well is suggested to be placed inside the large-permeability region regardless of the specific position.
  • A fully discrete energy stable scheme for a phase filed moving contact line model with variable densities and viscosities

    Zhu, Guangpu; Chen, Huangxin; Sun, Shuyu; Yao, Jun (arXiv, 2018-01-26) [Preprint]
    In this paper, a fully discrete scheme which considers temporal and spatial discretizations is presented for the coupled Cahn-Hilliard equation in conserved form with the dynamic contact line condition and the Navier-Stokes equation with the generalized Navier boundary condition. Variable densities and viscosities are incorporated in this model. A rigorous proof of energy stability is provided for the fully discrete scheme based on a semi-implicit temporal discretization and a finite difference method on the staggered grids for the spatial discretization. A splitting method based on the pressure stabilization is implemented to solve the Navier-Stokes equation, while the stabilization approach is also used for the Cahn-Hilliard equation. Numerical results in both 2-D and 3-D demonstrate the accuracy, efficiency and decaying property of discrete energy of the proposed scheme.
  • Physics-preserving averaging scheme based on Grunwald-Letnikov formula for gas flow in fractured media

    Amir, Sahar; Sun, Shuyu (Journal of Petroleum Science and Engineering, Elsevier BV, 2018-01-02) [Article]
    The heterogeneous natures of rock fabrics, due to the existence of multi-scale fractures and geological formations, led to the deviations from unity in the flux-equations fractional-exponent magnitudes. In this paper, the resulting non-Newtonian non-Darcy fractional-derivatives flux equations are solved using physics-preserving averaging schemes that incorporates both, original and shifted, Grunwald-Letnikov (GL) approximation formulas preserving the physics, by reducing the shifting effects, while maintaining the stability of the system, by keeping one shifted expansion. The proposed way of using the GL expansions also generate symmetrical coefficient matrices that significantly reduces the discretization complexities appearing with all shifted cases from literature, and help considerably in 2D and 3D systems. Systems equations derivations and discretization details are discussed. Then, the physics-preserving averaging scheme is explained and illustrated. Finally, results are presented and reviewed. Edge-based original GL expansions are unstable as also illustrated in literatures. Shifted GL expansions are stable but add a lot of additional weights to both discretization sides affecting the physical accuracy. In comparison, the physics-preserving averaging scheme balances the physical accuracy and stability requirements leading to a more physically conservative scheme that is more stable than the original GL approximation but might be slightly less stable than the shifted GL approximations. It is a locally conservative Single-Continuum averaging scheme that applies a finite-volume viewpoint.
  • Thermodynamically consistent modeling and simulation of multi-component two-phase flow with partial miscibility

    Kou, Jisheng; Sun, Shuyu (Computer Methods in Applied Mechanics and Engineering, Elsevier BV, 2017-12-09) [Article]
    A general diffuse interface model with a realistic equation of state (e.g. Peng-Robinson equation of state) is proposed to describe the multi-component two-phase fluid flow based on the principles of the NVT-based framework which is an attractive alternative recently over the NPT-based framework to model the realistic fluids. The proposed model uses the Helmholtz free energy rather than Gibbs free energy in the NPT-based framework. Different from the classical routines, we combine the first law of thermodynamics and related thermodynamical relations to derive the entropy balance equation, and then we derive a transport equation of the Helmholtz free energy density. Furthermore, by using the second law of thermodynamics, we derive a set of unified equations for both interfaces and bulk phases that can describe the partial miscibility of multiple fluids. A relation between the pressure gradient and chemical potential gradients is established, and this relation leads to a new formulation of the momentum balance equation, which demonstrates that chemical potential gradients become the primary driving force of fluid motion. Moreover, we prove that the proposed model satisfies the total (free) energy dissipation with time. For numerical simulation of the proposed model, the key difficulties result from the strong nonlinearity of Helmholtz free energy density and tight coupling relations between molar densities and velocity. To resolve these problems, we propose a novel convex-concave splitting of Helmholtz free energy density and deal well with the coupling relations between molar densities and velocity through very careful physical observations with a mathematical rigor. We prove that the proposed numerical scheme can preserve the discrete (free) energy dissipation. Numerical tests are carried out to verify the effectiveness of the proposed method.
  • Thermodynamically consistent simulation of nonisothermal diffuse-interface two-phase flow with Peng-Robinson equation of state

    Kou, Jisheng; Sun, Shuyu (arXiv, 2017-12-06) [Preprint]
    In this paper, we consider a diffuse-interface gas-liquid two-phase flow model with inhomogeneous temperatures, in which we employ the Peng-Robinson equation of state and the temperature-dependent influence parameter instead of the van der Waals equation of state and the constant influence parameter used in the existing models. As a result, our model can characterize accurately the physical behaviors of numerous realistic gas-liquid fluids, especially hydrocarbons. Furthermore, we prove a relation associating the pressure gradient with the gradients of temperature and chemical potential, and thereby derive a new formulation of the momentum balance equation, which shows that gradients of the chemical potential and temperature become the primary driving force of the fluid motion. It is rigorously proved that the new formulations of the model obey the first and second laws of thermodynamics. To design efficient numerical methods, we prove that Helmholtz free energy density is a concave function with respect to the temperature under certain physical conditions. Based on the proposed modeling formulations and the convex-concave splitting of Helmholtz free energy density, we propose a novel thermodynamically stable numerical scheme. We rigorously prove that the proposed method satisfies the first and second laws of thermodynamics. Finally, numerical tests are carried out to verify the effectiveness of the proposed simulation method.
  • Linearly decoupled energy-stable numerical methods for multi-component two-phase compressible flow

    Kou, Jisheng; Sun, Shuyu; Wang, Xiuhua (arXiv, 2017-12-06) [Preprint]
    In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure a consistency relationship with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods have the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.
  • Direct numerical simulation of noninvasive channel healing in electrical field

    Wang, Yi; Sun, Shuyu (Advances in Mechanical Engineering, SAGE Publications, 2017-11-25) [Article]
    Noninvasive channel healing is a new idea to repair the broken pipe wall, using external electric fields to drive iron particles to the destination. The repair can be done in the normal operation of the pipe flow without any shutdown of the pipeline so that this method can be a potentially efficient and safe technology of pipe healing. However, the real application needs full knowledge of healing details. Numerical simulation is an effective method. Thus, in this research, we first established a numerical model for noninvasive channel healing technology to represent fluid–particle interaction. The iron particles can be attached to a cracking area by external electrostatic forces or can also be detached by mechanical forces from the fluid. When enough particles are permanently attached on the cracking area, the pipe wall can be healed. The numerical criterion of the permanent attachment is discussed. A fully three-dimensional finite difference framework of direct numerical simulation is established and applied to different cases to simulate the full process of channel healing. The impact of Reynolds number and particle concentration on the healing process is discussed. This numerical investigation provides valuable reference and tools for further simulation of real pipe healing in engineering.
  • Discrete-fracture-model of multi–scale time-splitting two–phase flow including nanoparticles transport in fractured porous media

    El-Amin, Mohamed; Kou, Jisheng; Sun, Shuyu (Journal of Computational and Applied Mathematics, Elsevier BV, 2017-11-23) [Article]
    In this article, we consider a two-phase immiscible incompressible flow including nanoparticles transport in fractured heterogeneous porous media. The system of the governing equations consists of water saturation, Darcy’s law, nanoparticles concentration in water, deposited nanoparticles concentration on the pore-wall, and entrapped nanoparticles concentration in the pore-throat, as well as, porosity and permeability variation due to the nanoparticles deposition/entrapment on/in the pores. The discrete-fracture model (DFM) is used to describe the flow and transport in fractured porous media. Moreover, multiscale time-splitting strategy has been employed to manage different time-step sizes for different physics, such as saturation, concentration, etc. Numerical examples are provided to demonstrate the efficiency of the proposed multi-scale time splitting approach.
  • Molecular Dynamics Simulation Study of Carbon Dioxide, Methane, and Their Mixture in the Presence of Brine

    Yang, Yafan; Nair, Arun Kumar Narayanan; Sun, Shuyu (The Journal of Physical Chemistry B, American Chemical Society (ACS), 2017-10-03) [Article]
    We perform molecular dynamics simulation study of CO2, methane, and their mixture in the presence of brine over a broad range of temperature (311–473 K), pressure (up to about 100 MPa), and NaCl concentration (up to about 14 wt %). The general decrease in the interfacial tension (IFT) values of the CH4–brine system with pressure and temperature is similar to that obtained for the corresponding CH4–water system. The IFT of methane and brine is a linearly increasing function of salt concentration, and the resulting slopes are dependent on the pressure. A similar behavior as methane is observed for such systems containing CO2 and CO2–CH4 mixture. The IFT of CO2 and brine increases linearly with increasing salt content; however, the resulting slopes are independent of pressure. The simulations show that the presence of CO2 decreases the IFT values of the CH4–water and CH4–brine systems, but the degree of reduction depends on the amount of CO2 in each sample, which is consistent with experimental evidence. These IFT values show a linear correlation with the amount of CO2, and the resulting slopes are dependent on the temperature and pressure. Furthermore, our results for the mole fractions of the different species in the CO2–CH4–water system at 323 K and 9 MPa are in agreement with those of experiments. The mole fractions of methane and CO2 in the water-rich phase decrease with increasing salt concentration, whereas that of H2O in the methane- or CO2-rich phases remains almost unaffected in all of the studied cases. Our results could be useful because of the importance of carbon dioxide sequestration and shale gas production.
  • A stable algorithm for calculating phase equilibria with capillarity at specified moles, volume and temperature using a dynamic model

    Kou, Jisheng; Sun, Shuyu (Fluid Phase Equilibria, Elsevier BV, 2017-09-30) [Article]
    Capillary pressure can significantly affect the phase properties and flow of liquid-gas fluids in porous media, and thus, the phase equilibrium calculation incorporating capillary pressure is crucial to simulate such problems accurately. Recently, the phase equilibrium calculation at specified moles, volume and temperature (NVT-flash) becomes an attractive issue. In this paper, capillarity is incorporated into the phase equilibrium calculation at specified moles, volume and temperature. A dynamical model for such problem is developed for the first time by using the laws of thermodynamics and Onsager's reciprocal principle. This model consists of the evolutionary equations for moles and volume, and it can characterize the evolutionary process from a non-equilibrium state to an equilibrium state in the presence of capillarity effect at specified moles, volume and temperature. The phase equilibrium equations are naturally derived. To simulate the proposed dynamical model efficiently, we adopt the convex-concave splitting of the total Helmholtz energy, and propose a thermodynamically stable numerical algorithm, which is proved to preserve the second law of thermodynamics at the discrete level. Using the thermodynamical relations, we derive a phase stability condition with capillarity effect at specified moles, volume and temperature. Moreover, we propose a stable numerical algorithm for the phase stability testing, which can provide the feasible initial conditions. The performance of the proposed methods in predicting phase properties under capillarity effect is demonstrated on various cases of pure substance and mixture systems.
  • Acceleration of Gas Flow Simulations in Dual-Continuum Porous Media Based on the Mass-Conservation POD Method

    Wang, Yi; Sun, Shuyu; Yu, Bo (Energies, MDPI AG, 2017-09-12) [Article]
    Reduced-order modeling approaches for gas flow in dual-porosity dual-permeability porous media are studied based on the proper orthogonal decomposition (POD) method combined with Galerkin projection. The typical modeling approach for non-porous-medium liquid flow problems is not appropriate for this compressible gas flow in a dual-continuum porous media. The reason is that non-zero mass transfer for the dual-continuum system can be generated artificially via the typical POD projection, violating the mass-conservation nature and causing the failure of the POD modeling. A new POD modeling approach is proposed considering the mass conservation of the whole matrix fracture system. Computation can be accelerated as much as 720 times with high precision (reconstruction errors as slow as 7.69 × 10−4%~3.87% for the matrix and 8.27 × 10−4%~2.84% for the fracture).
  • Energy Stability Analysis of Some Fully Discrete Numerical Schemes for Incompressible Navier–Stokes Equations on Staggered Grids

    Chen, Huangxin; Sun, Shuyu; Zhang, Tao (Journal of Scientific Computing, Springer Nature, 2017-09-01) [Article]
    In this paper we consider the energy stability estimates for some fully discrete schemes which both consider time and spatial discretizations for the incompressible Navier–Stokes equations. We focus on three kinds of fully discrete schemes, i.e., the linear implicit scheme for time discretization with the finite difference method (FDM) on staggered grids for spatial discretization, pressure-correction schemes for time discretization with the FDM on staggered grids for the solutions of the decoupled velocity and pressure equations, and pressure-stabilization schemes for time discretization with the FDM on staggered grids for the solutions of the decoupled velocity and pressure equations. The energy stability estimates are obtained for the above each fully discrete scheme. The upwind scheme is used in the discretization of the convection term which plays an important role in the design of unconditionally stable discrete schemes. Numerical results are given to verify the theoretical analysis.
  • Modeling and Analysis of Magnetic Nanoparticles Injection in Water-Oil Two-Phase Flow in Porous Media under Magnetic Field Effect

    El-Amin, Mohamed; Saad, Adel; Salama, Amgad; Sun, Shuyu (Geofluids, Hindawi Limited, 2017-08-28) [Article]
    In this paper, the magnetic nanoparticles are injected into a water-oil, two-phase system under the influence of an external permanent magnetic field. We lay down the mathematical model and provide a set of numerical exercises of hypothetical cases to show how an external magnetic field can influence the transport of nanoparticles in the proposed two-phase system in porous media. We treat the water-nanoparticles suspension as a miscible mixture, whereas it is immiscible with the oil phase. The magnetization properties, the density, and the viscosity of the ferrofluids are obtained based on mixture theory relationships. In the mathematical model, the phase pressure contains additional term to account for the extra pressures due to fluid magnetization effect and the magnetostrictive effect. As a proof of concept, the proposed model is applied on a countercurrent imbibition flow system in which both the displacing and the displaced fluids move in opposite directions. Physical variables, including waternanoparticles suspension saturation, nanoparticles concentration, and pore wall/throat concentrations of deposited nanoparticles, are investigated under the influence of the magnetic field. Two different locations of the magnet are studied numerically, and variations in permeability and porosity are considered.
  • Complexation Behavior of Polyelectrolytes and Polyampholytes

    Nair, Arun Kumar Narayanan; Jimenez, Arturo Martinez; Sun, Shuyu (The Journal of Physical Chemistry B, American Chemical Society (ACS), 2017-07-25) [Article]
    We perform grand canonical Monte Carlo simulations to study the pH titrations of isolated polyampholytes and polyelectrolyte-polyampholyte complexes in dilute solutions. Our simulations indicate that the electrostatic interactions promote the coexistence of opposite charges along the polyampholyte chain during titration. The repulsion between excess charges typically dominates the electrostatic interaction and leads to polymer stretching. Salt ions can screen the repulsion between excess charges as well as the fluctuation-induced attraction between opposite charges, and therefore make the variation between titration curves of polyampholytes and the ideal (no electrostatic interactions) curves less significant. We observe that this screening of charge repulsion decreases the chain size. The presence of pearl-necklace configuration of polyampholytes is diminished by the addition of salt. Similar simulations for the polyelectrolyte-polyampholyte system show that the resulting complexes are generally stable in the low pH region. In comparison to ideal case, electrostatic interactions strongly influence the acid-base properties of polyampholyte chains in the adsorbed state by reducing the presence of the coexistence domain of both positive and negative charges in the titration curves. We attribute the complex formation between polyelectrolyte and polyampholyte chains in the high pH region to, e.g., the high salt content. The pH variation leads to abrupt transition between adsorbed and desorbed states. Independent of charge sequence, a polyampholyte chain in a complex is usually located at one of the ends of the polyelectrolyte chain.
  • Reduced Fracture Finite Element Model Analysis of an Efficient Two-Scale Hybrid Embedded Fracture Model

    Amir, Sahar; Chen, Huangxin; Sun, Shuyu (Procedia Computer Science, Elsevier BV, 2017-06-09) [Article]
    A Hybrid Embedded Fracture (HEF) model was developed to reduce various computational costs while maintaining physical accuracy (Amir and Sun, 2016). HEF splits the computations into fine scale and coarse scale. Fine scale solves analytically for the matrix-fracture flux exchange parameter. Coarse scale solves for the properties of the entire system. In literature, fractures were assumed to be either vertical or horizontal for simplification (Warren and Root, 1963). Matrix-fracture flux exchange parameter was given few equations built on that assumption (Kazemi, 1968; Lemonnier and Bourbiaux, 2010). However, such simplified cases do not apply directly for actual random fracture shapes, directions, orientations …etc. This paper shows that the HEF fine scale analytic solution (Amir and Sun, 2016) generates the flux exchange parameter found in literature for vertical and horizontal fracture cases. For other fracture cases, the flux exchange parameter changes according to the angle, slop, direction, … etc. This conclusion rises from the analysis of both: the Discrete Fracture Network (DFN) and the HEF schemes. The behavior of both schemes is analyzed with exactly similar fracture conditions and the results are shown and discussed. Then, a generalization is illustrated for any slightly compressible single-phase fluid within fractured porous media and its results are discussed.
  • Dual-mixed finite elements for the three-field Stokes model as a finite volume method on staggered grids

    Kou, Jisheng; Sun, Shuyu (Procedia Computer Science, Elsevier BV, 2017-06-09) [Article]
    In this paper, a new three-field weak formulation for Stokes problems is developed, and from this, a dual-mixed finite element method is proposed on a rectangular mesh. In the proposed mixed methods, the components of stress tensor are approximated by piecewise constant functions or Q1 functions, while the velocity and pressure are discretized by the lowest-order Raviart-Thomas element and the piecewise constant functions, respectively. Using quadrature rules, we demonstrate that this scheme can be reduced into a finite volume method on staggered grid, which is extensively used in computational fluid mechanics and engineering.
  • Molecular Simulation Study of Montmorillonite in Contact with Variably Wet Supercritical Carbon Dioxide

    Kadoura, Ahmad Salim; Nair, Arun Kumar Narayanan; Sun, Shuyu (The Journal of Physical Chemistry C, American Chemical Society (ACS), 2017-03-07) [Article]
    We perform grand canonical Monte Carlo simulations to study the detailed molecular mechanism of intercalation behavior of CO2 in Na-, Ca-, and Mg- montmorillonite exposed to variably hydrated supercritical CO2 at 323.15 K and 90 bar, The simulations indicate that the intercalation of CO2 strongly depends on the relative humidity (RH). The intercalation of CO2 in the dehydrated interlayer is inhibited, followed by the swelling of the interlayer region due to uptake of water and CO2 as the RH increases. In all of the hydrated clay samples, the amount of the intercalated CO2 generally decreases as a function of increasing RH, which is attributed mainly to the weakening of the interaction between CO2 and clay. At low RH values, Ca- and Mg- montmorillonite are relatively more efficient in capturing CO2. The amount of CO2 trapped in all clay samples shows similar values above RH of similar to 60%. Molecular dynamics simulations show that the diffusion coefficient of each species generally increases with increasing RH due to the associated expansion of the interlayer distance of the clay. For all the hydrated samples, the diffusion coefficients of CO2 and water in the interlayers are mostly comparable due to the fact that CO2 molecules are well solvated. The diffusion of CO2 in each hydration state is mostly independent of the type of cation in accordance with the fact that CO2 molecules hardly migrate into the first hydration shell of the interlayer cations.
  • POD-Galerkin Model for Incompressible Single-Phase Flow in Porous Media

    Wang, Yi; Yu, Bo; Sun, Shuyu (Open Physics, Walter de Gruyter GmbH, 2017-01-25) [Article]
    Fast prediction modeling via proper orthogonal decomposition method combined with Galerkin projection is applied to incompressible single-phase fluid flow in porous media. Cases for different configurations of porous media, boundary conditions and problem scales are designed to examine the fidelity and robustness of the model. High precision (relative deviation 1.0 x 10(-4)% similar to 2.3 x 10(-1)%) and large acceleration (speed-up 880 similar to 98454 times) of POD model are found in these cases. Moreover, the computational time of POD model is quite insensitive to the complexity of problems. These results indicate POD model is especially suitable for large-scale complex problems in engineering.
  • An Efficient Two-Scale Hybrid Embedded Fracture Model for Shale Gas Simulation

    Amir, Sahar; Sun, Shuyu (Journal of Petroleum Science and Engineering, Elsevier BV, 2016-12-27) [Article]
    Natural and hydraulic fractures existence and state differs on a reservoir-by-reservoir or even on a well-by-well basis leading to the necessity of exploring the flow regimes variations with respect to the diverse fracture-network shapes forged. Conventional Dual-Porosity Dual-Permeability (DPDP) schemes are not adequate to model such complex fracture-network systems. To overcome this difficulty, in this paper, an iterative Hybrid Embedded multiscale (two-scale) Fracture model (HEF) is applied on a derived fit-for-purpose shale gas model. The HEF model involves splitting the fracture computations into two scales: 1) fine-scale solves for the flux exchange parameter within each grid cell; 2) coarse-scale solves for the pressure applied to the domain grid cells using the flux exchange parameter computed at each grid cell from the fine-scale. After that, the D dimensions matrix pressure and the (D-1) lower dimensional fracture pressure are solved as a system to apply the matrix-fracture coupling. HEF model combines the DPDP overlapping continua concept, the DFN lower dimensional fractures concept, the HFN hierarchical fracture concept, and the CCFD model simplicity. As for the fit-for-purpose shale gas model, various fit-for-purpose shale gas models can be derived using any set of selected properties plugged in one of the most popularly used proposed literature models as shown in the appendix. Also, this paper shows that shale extreme low permeability cause flow behavior to be dominated by the structure and magnitude of high permeability fractures.

View more