... with the rate at which data can be brought to the floating-point unit. In reservoir simulation, with the need for sparse linear algebra and unstructured grids, algorithmic intensity is low and memory bandwidth is key. Alternative hardware to the CPU for numerical applications is actively under development...

... as much computations to GPU as possible. Our CUDA implementation covers all reservoir computes, which include property calculation, linearization, linear solver, etc. The well and Field Management still reside on CPU and need minor changes for their interaction with GPU-based reservoir. Importantly...

... subspace, is combined with tra jectory piecewise linearization (TPWL), where solutions with new sets of well controls are represented via linearization around previously simulated (training) solutions. The over-determined system of equations is pro jected into the lowdimensional subspace using a least...

... has to be inevitably rewritten, sometimes almost completely. We demonstrate how to extract complex physics-related computations from the main simulation loop, leaving only an algebraic multilinear interpolation kernel instead. In combination with linear solvers, which usually have made available soon...

... in the subsurface, takes the most of the simulation time. The linearization part of the solution process occupies a significant fraction of that time, especially in compositional models. Moreover, the implementation of the linearization step usually embodies the most substantial, complex, and specific part...

... Our goal in this paper is to develop a very efficient reduced-order modeling technique that can be used as a numerical surrogate model in compositional simulation. The method, referred to as POD-TPWL, is based on trajectory piecewise linearization; i.e., representation of new solutions...

... The key idea of trajectory piecewise linearization is to linearize Eq. 1 around states saved during previous (training) simulations. The goal in this work is to generate training runs using a particular geological model, which we designate as m 0 , and then apply TPWL to simulate models...