Source author record

Hamdi A. Tchelepi

Hamdi A. Tchelepi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2020arXiv

Algebraically stabilized Lagrange multiplier method for frictional contact mechanics with hydraulically active fractures

Accurate numerical simulation of coupled fracture/fault deformation and fluid flow is crucial to the performance and safety assessment of many subsurface systems. In this work, we consider the discretization and enforcement of contact conditions at such surfaces. The bulk rock deformation is simulated using low-order continuous finite elements, while frictional contact conditions are imposed by means of a Lagrange multiplier method. We employ a cell-centered finite-volume scheme to solve the fracture fluid mass balance equation. From a modeling perspective, a convenient choice is to use a single grid for both mechanical and flow processes, with piecewise-constant interpolation of Lagrange multipliers, i.e., contact tractions and fluid pressure. Unfortunately, this combination of displacement and multiplier variables is not uniformly inf-sup stable, and therefore requires a stabilization technique. Starting from a macroelement analysis, we develop two algebraic stabilization approaches and compare them in terms of robustness and convergence rate. The proposed approaches are validated against challenging analytical two- and three-dimensional benchmarks to demonstrate accuracy and robustness. These benchmarks include both pure contact mechanics problems and well as problems with tightly-coupled fracture flow.

preprint2020arXiv

Dissipation-Based Continuation Method for Multiphase Flow in Heterogeneous Porous Media

In reservoir simulation, solution of the coupled systems of nonlinear algebraic equations that are associated with fully-implicit (backward Euler) discretization is challenging. Having a robust and efficient nonlinear solver is necessary in order for reservoir simulation to serve as the primary tool for managing the recovery processes of large-scale reservoirs. Here, we develop a continuation method based on the use of a dissipation operator. We focus on nonlinear two-phase flow and transport in heterogeneous formations in the presence of viscous, gravitational, and capillary forces. The homotopy is constructed by adding numerical dissipation to the coupled discrete conservation equations. A continuation parameter is introduced to control the amount of dissipation. Numerical evidence of multi-dimensional models and detailed analysis of single-cell problems are used to explain how the dissipation operator improves the nonlinear convergence of the coupled system of equations. An adaptive strategy to determine the dissipation coefficient is proposed. The dissipation level is computed locally for each cell interface. We demonstrate the efficiency of the dissipation-based continuation (DBC) nonlinear solver using several examples, including 1D scalar transport and 2D heterogeneous problems with fully-coupled flow and transport. The DBC solver has better convergence properties compared with the standard damped-Newton solvers used in reservoir simulation.

preprint2020arXiv

MeshfreeFlowNet: A Physics-Constrained Deep Continuous Space-Time Super-Resolution Framework

We propose MeshfreeFlowNet, a novel deep learning-based super-resolution framework to generate continuous (grid-free) spatio-temporal solutions from the low-resolution inputs. While being computationally efficient, MeshfreeFlowNet accurately recovers the fine-scale quantities of interest. MeshfreeFlowNet allows for: (i) the output to be sampled at all spatio-temporal resolutions, (ii) a set of Partial Differential Equation (PDE) constraints to be imposed, and (iii) training on fixed-size inputs on arbitrarily sized spatio-temporal domains owing to its fully convolutional encoder. We empirically study the performance of MeshfreeFlowNet on the task of super-resolution of turbulent flows in the Rayleigh-Benard convection problem. Across a diverse set of evaluation metrics, we show that MeshfreeFlowNet significantly outperforms existing baselines. Furthermore, we provide a large scale implementation of MeshfreeFlowNet and show that it efficiently scales across large clusters, achieving 96.80% scaling efficiency on up to 128 GPUs and a training time of less than 4 minutes.

preprint2020arXiv

Nonlinear Acceleration of Sequential Fully Implicit (SFI) Method for Coupled Flow and Transport in Porous Media

The sequential fully implicit (SFI) method was introduced along with the development of the multiscale finite volume (MSFV) framework, and has received considerable attention in recent years. Each time step for SFI consists of an outer loop to solve the coupled system, in which there is one inner Newton loop to implicitly solve the pressure equation and another loop to implicitly solve the transport equations. Limited research has been conducted that deals with the outer coupling level to investigate the convergence performance. In this paper we extend the basic SFI method with several nonlinear acceleration techniques for improving the outer-loop convergence. Specifically, we consider numerical relaxation, quasi-Newton (QN) and Anderson acceleration (AA) methods. The acceleration techniques are adapted and studied for the first time within the context of SFI for coupled flow and transport in porous media. We reveal that the iterative form of SFI is equivalent to a nonlinear block Gauss-Seidel (BGS) process. The effectiveness of the acceleration techniques is demonstrated using several challenging examples. The results show that the basic SFI method is quite inefficient, suffering from slow convergence or even convergence failure. In order to better understand the behaviors of SFI, we carry out detailed analysis on the coupling mechanisms between the sub-problems. Compared with the basic SFI method, superior convergence performance is achieved by the acceleration techniques, which can resolve the convergence difficulties associated with various types of coupling effects. We show across a wide range of flow conditions that the acceleration techniques can stabilize the iterative process, and largely reduce the outer iteration count.

preprint2016arXiv

Stochastic Galerkin Framework with Locally Reduced Bases for Nonlinear Two-Phase Transport in Heterogeneous Formations

The generalized polynomial chaos method is applied to the Buckley-Leverett equation. We consider a spatially homogeneous domain modeled as a random field. The problem is projected onto stochastic basis functions which yields an extended system of partial differential equations. Analysis and numerical methods leading to reduced computational cost are presented for the extended system of equations. The accurate representation of the evolution of a discontinuous stochastic solution over time requires a large number of stochastic basis functions. Adaptivity of the stochastic basis to reduce computational cost is challenging in the stochastic Galerkin setting since the change of basis affects the system matrix itself. To achieve adaptivity without adding overhead by rewriting the entire system of equations for every grid cell, we devise a basis reduction method that distinguishes between locally significant and insignificant modes without changing the actual system matrices. Results are presented for problems in one and two spatial dimensions, with varying number of stochastic dimensions. We show how to obtain stochastic velocity fields from realistic permeability fields and demonstrate the performance of the stochastic Galerkin method with local basis reduction. The system of conservation laws is discretized with a finite volume method and we demonstrate numerical convergence to the reference solution obtained through Monte Carlo sampling.