Christopher Lemon, H. Cao, M. D. E. Szyndel, Eduard Khramchenkov
{"title":"An Adaptive Coloring Scheme for Graphics Processing Unit Preconditioners","authors":"Christopher Lemon, H. Cao, M. D. E. Szyndel, Eduard Khramchenkov","doi":"10.2118/212248-ms","DOIUrl":"https://doi.org/10.2118/212248-ms","url":null,"abstract":"\u0000 A single modern graphics processing unit (GPU) typically has the memory bandwidth equivalent to many central processing unit (CPU) nodes. This makes GPU hardware appealing for linear solvers that tend to require high memory bandwidth and fast inter-core communication. Reservoir simulators are designed to handle a wide range of simulation models, and to obtain peak performance the linear solver must be well suited to the resulting linear systems. This fact can lead to disappointing performance when shifting the linear solver from CPU to GPU. To fully utilize the capabilities of the latest GPU devices, we must transition from coarse-grained to fine-grained parallel preconditioners. To enable such high levels of parallelism in the linear solver a common approach is to employ a multicolor reordering of the linear system. Depending on the specific properties of the simulation model, this process can cause a significant weakening of the parallel preconditioner, resulting in much slower convergence. In some situations, this slow convergence can cause an order of magnitude increase in the linear iteration count, and result in the GPU linear solver performing worse than the CPU version. In this paper we analyze the impact on performance of employing different coloring schemes for different simulation models and we identify how the coloring can be automatically adapted for the properties of each simulation model. In this way, the performance improvements expected on the GPU can be realized for a wider range of simulations.","PeriodicalId":225811,"journal":{"name":"Day 1 Tue, March 28, 2023","volume":"22 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"129899407","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Hybrid-Dimensional Darcy's Law: A Non-Conforming Reinterpreted Discrete Fracture Model (RDFM) for the Compressible Miscible Displacement and Multicomponent Gas Flow in Fractured Media","authors":"Ziyao Xu, Yang Yang","doi":"10.2118/212164-ms","DOIUrl":"https://doi.org/10.2118/212164-ms","url":null,"abstract":"\u0000 In this work, we develop a non-conforming reinterpreted discrete fracture model for the compressible miscible displacement and multicomponent gas flow in porous media containing high-permeability fractures and/or low-permeability barriers based on the hybrid-dimensional Darcy's law established in our previous work.\u0000 The key idea of the model is to describe the permeability of codimension-one fractures and barriers by the Dirac-delta functions. When there are only fractures, delta functions are added to the permeability tensor on the right-hand side of the Darcy's law. In contrast, when there are only barriers, delta functions are added to the inverse of the permeability tensor, which represents the resistance to fluids, on the left-hand side of the Darcy's law. When both appear, delta functions are contained on both sides by the principle of superposition. Thereby, we establish partial differential equations (PDEs) to model fluid flow in fractured porous media, which exempts any requirements on meshes.\u0000 We adopt the discontinuous Galerkin (DG) method to discretize the model in space and the second order implicit pressure explicit concentration (SIMPEC) method to march in time. The resulting non-conforming discrete fracture model is local mass conservative, flexible for complex geometry and easy to implement. The good performance of the method is demonstrated by several numerical examples.","PeriodicalId":225811,"journal":{"name":"Day 1 Tue, March 28, 2023","volume":"252 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"127541721","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}