Abstract: Relaxation process of magnetohydrodynamics (MHD) inside a sphere is investigated by a newly developed spherical grid system, Yin-Yang-Zhong grid. An MHD fluid with low dissipation rates is confined by a perfectly conducting, stress-free, and thermally insulating spherical boundary. The Reynolds number Re and the magnetic Reynolds number Rm is the same: Re=Rm=8600. Starting from a simple and symmetric state in which a ring-shaped magnetic flux without flow, a dynamical relaxation process of the magnetic energy is numerically integrated. The relaxed state has a characteristic structure of the magnetic field and the flow field with four vortices.

Abstract: Nanoparticles are particles that are between 1 and 100 nanometers in size. They present possible dangers to the environment due to the high surface to volume ratio, which can make the particles very reactive or catalytic. Furthermore, rapid increase in the implementation of nanotechnologies has released large amount of the nanowaste into the environment. In the last two decades, transport of nanoparticles in the subsurface and the potential hazard they impose to the environment have attracted the attention of researchers. In this work, we use numerical simulation to investigate the problem regarding the transport phenomena of nanoparticles in anisotropic porous media. We consider the case in which the permeability in the principal direction components will vary with respect to time. The interesting thing in this case is the fact that the anisotropy could disappear with time. We investigate the effect of the degenerating anisotropy on various fields such as pressure, porosity, concentration and velocities.

Abstract: Fully-Implicit (FI) Methods are often employed in the numerical simulation of large-scale subsurface flows in porous media. At each implicit time step, a Newton-like method is used to solve the FI discrete nonlinear algebraic system. The linear solution process for the Newton updates is the computational workhorse of FI simulations. Empirical observations suggest that the computed Newton updates during FI simulations of multiphase flow are often sparse. Moreover, the level of sparsity observed can vary dramatically from iteration to the next, and across time steps. In several large scale applications, it was reported that the level of sparsity in the Newton update can be as large as 99\%. This work develops a localization algorithm that conservatively predetermines the sparsity pattern of the Newton update. Subsequently, only the flagged nonzero components of the system need be solved. The localization algorithm is developed for general FI models of two phase flow. Large scale simulation results of benchmark reservoir models show a 10 to 100 fold reduction in computational cost for homogeneous problems, and a 4 to 10 fold reduction for strongly heterogeneous problems.

Abstract: Unconventional reservoirs are typically comprised of a multicontinuum stimulated formation, with complex fracture networks that have a wide range of length scales and geometries. A timely topic in the simulation of unconventional petroleum resources is in coupling the geomechanics of the fractured media to multiphase fluid flow and transport. We propose a XFEM-EDFM method which couples geomechanics with multiphase flow in fractured tight gas reservoirs. A proppant model is developed to simulate propped hydraulic fractures. The method is verified by analytical solutions. A simulation example with the configuration of two multiple-fractured horizontal wells is investigated. The influence of stress-dependent fracture permeability on cumulative production is analyzed.