GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation

This paper designs a highly parallel Nested Factorization (NF) to solve large linear equations generated in reservoir numerical simulation problems. The NF method is a traditional linear solution preprocessing method for reservoir numerical simulation problems, and has regained attention in recent y...

Full description

Bibliographic Details
Main Authors: Xin Shi, Yuan Di
Format: Article
Language:English
Published: Universidad Nacional de Colombia 2019-07-01
Series:Earth Sciences Research Journal
Subjects:
cpr
Online Access:https://revistas.unal.edu.co/index.php/esrj/article/view/81669
id doaj-20b798d66a5e44369b51a5c0ed3905ea
record_format Article
spelling doaj-20b798d66a5e44369b51a5c0ed3905ea2020-11-25T01:13:46ZengUniversidad Nacional de ColombiaEarth Sciences Research Journal1794-61902019-07-0123324925710.15446/esrj.v23n3.8166951764GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical SimulationXin Shi0Yuan Di1Peking University - College of EngineeringPeking University - College of EngineeringThis paper designs a highly parallel Nested Factorization (NF) to solve large linear equations generated in reservoir numerical simulation problems. The NF method is a traditional linear solution preprocessing method for reservoir numerical simulation problems, and has regained attention in recent years due to its potential to extend to parallel architectures such as GPUs (Graphics Processor Units). The parallel algorithm of this paper is based on the MPNF (Massively Parallel Nested Factorization) framework proposed by Appleya (Appleyard, Appleyard, Wakefield, & Desitter, 2011). The MPNF algorithm designed in this paper focuses on its efficient implementation on the GPU parallel architecture. Its features include: using a custom matrix structure to achieve merge access, improving access bottlenecks and improving the efficiency of the SpMV algorithm. It is also applicable to the two-stage preprocessing method CPR. (Constrain Pressure Residual) pressure solution and global preprocessing stage; the MPNF method is extended to the solution of 2.5-dimensional unstructured grid problem. The parallel algorithm in this paper has been integrated into the reservoir numerical simulator. For the SPE10 (million grid, highly heterogeneous) standard example, the GPU-based parallel NF algorithm is in the structured grid model and the equivalent 2.5-dimensional non- On the structured grid model, compared with the serial version of the NF method, the acceleration ratios of 19.8 and 17.0 times were obtained respectively; compared with the mainstream serial solution method, the efficiency was also improved by 2 to 3 times.https://revistas.unal.edu.co/index.php/esrj/article/view/81669nested decompositiongpu parallellinear solutioncpr
collection DOAJ
language English
format Article
sources DOAJ
author Xin Shi
Yuan Di
spellingShingle Xin Shi
Yuan Di
GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
Earth Sciences Research Journal
nested decomposition
gpu parallel
linear solution
cpr
author_facet Xin Shi
Yuan Di
author_sort Xin Shi
title GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
title_short GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
title_full GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
title_fullStr GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
title_full_unstemmed GPU Parallelization Nested Decomposition Method for Solving Large Linear Systems in Reservoir Numerical Simulation
title_sort gpu parallelization nested decomposition method for solving large linear systems in reservoir numerical simulation
publisher Universidad Nacional de Colombia
series Earth Sciences Research Journal
issn 1794-6190
publishDate 2019-07-01
description This paper designs a highly parallel Nested Factorization (NF) to solve large linear equations generated in reservoir numerical simulation problems. The NF method is a traditional linear solution preprocessing method for reservoir numerical simulation problems, and has regained attention in recent years due to its potential to extend to parallel architectures such as GPUs (Graphics Processor Units). The parallel algorithm of this paper is based on the MPNF (Massively Parallel Nested Factorization) framework proposed by Appleya (Appleyard, Appleyard, Wakefield, & Desitter, 2011). The MPNF algorithm designed in this paper focuses on its efficient implementation on the GPU parallel architecture. Its features include: using a custom matrix structure to achieve merge access, improving access bottlenecks and improving the efficiency of the SpMV algorithm. It is also applicable to the two-stage preprocessing method CPR. (Constrain Pressure Residual) pressure solution and global preprocessing stage; the MPNF method is extended to the solution of 2.5-dimensional unstructured grid problem. The parallel algorithm in this paper has been integrated into the reservoir numerical simulator. For the SPE10 (million grid, highly heterogeneous) standard example, the GPU-based parallel NF algorithm is in the structured grid model and the equivalent 2.5-dimensional non- On the structured grid model, compared with the serial version of the NF method, the acceleration ratios of 19.8 and 17.0 times were obtained respectively; compared with the mainstream serial solution method, the efficiency was also improved by 2 to 3 times.
topic nested decomposition
gpu parallel
linear solution
cpr
url https://revistas.unal.edu.co/index.php/esrj/article/view/81669
work_keys_str_mv AT xinshi gpuparallelizationnesteddecompositionmethodforsolvinglargelinearsystemsinreservoirnumericalsimulation
AT yuandi gpuparallelizationnesteddecompositionmethodforsolvinglargelinearsystemsinreservoirnumericalsimulation
_version_ 1725160146827476992