An advanced fully-implicit solver for heterogeneous porous media based on foam-extend

IF 3.4 2区物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computer Physics Communications Pub Date : 2025-09-09 DOI:10.1016/j.cpc.2025.109842

Roberto Lange, Gabriel M. Magalhães, Franciane F. Rocha, Hélio Ribeiro Neto

{"title":"An advanced fully-implicit solver for heterogeneous porous media based on foam-extend","authors":"Roberto Lange, Gabriel M. Magalhães, Franciane F. Rocha, Hélio Ribeiro Neto","doi":"10.1016/j.cpc.2025.109842","DOIUrl":null,"url":null,"abstract":"<div><div>Multiphase flow in porous media is present in many engineering applications, including hydrogeology, oil recovery, and CO<sub>2</sub> sequestration. Accurate predictions of fluid behavior in these systems can improve process efficiency while mitigating environmental and health risks. Commercial simulators and open source software, such as the <span>porousMultiphaseFoam</span> repository based on the OpenFOAM framework, have been developed to model this type of problem. However, simulating heterogeneous porous media with heterogeneous porosity and permeability distributions poses significant numerical challenges. We introduce <span>coupledMatrixFoam</span>, an OpenFOAM-based solver designed for enhanced numerical stability and robustness. <span>coupledMatrixFoam</span> integrates the Eulerian multi-fluid formulation for phase fractions with Darcy's law for porous media flow, applying a fully implicit, block-coupled solution for pressure and phase fractions. The solver is based on foam-extend 5.0, leveraging the latest <span>fvBlockMatrix</span> developments to improve computational efficiency. This approach enables a significant increase in time step sizes, particularly in cases involving capillary pressure effects and other complex physical interactions. This work details the formulation, implementation and validation of <span>coupledMatrixFoam</span>, including comparisons with <span>porousMultiphaseFoam</span> that uses a segregated approach, to assess performance improvements. Additionally, a scalability analysis is conducted, demonstrating the solver's ability for high-performance computing (HPC) applications, which are essential for large-scale, real-world simulations.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> coupledMatrixFoam</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/3d3xdh4x89.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link:</em> <span><span>https://gitlab.com/wikki.brasil/porousmedia</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> GPLv3</div><div><em>Programming language:</em> C++</div><div><em>Supplementary material:</em> Available in the repository <span><span>porousMedia</span><svg><path></path></svg></span>.</div><div><em>Nature of problem:</em> This software solves multiphase flow in porous media.</div><div><em>Solution method:</em> Fully implicit solver based on the Euler-Euler multifluid formulation combined with Darcy's law developed in the OpenFOAM framework, that is based on the Finite Volume Method (FVM). Implicit coupling of phase fraction and pressure equations allowing significantly larger time steps. Complex physical phenomena are accounted for, including capillary pressure effects, gravitational forces, compressibility, and heterogeneous media properties. The porous medium is modeled as a stationary phase, and nonlinear terms are linearized using Taylor series expansions. Parallel computation is supported through OpenFOAM's standard domain decomposition approach.</div><div><em>Additional comments including restrictions and unusual features:</em> Compatibility with foam-extend (OpenFOAM) version 5.0.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"317 ","pages":"Article 109842"},"PeriodicalIF":3.4000,"publicationDate":"2025-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465525003443","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

Abstract

Multiphase flow in porous media is present in many engineering applications, including hydrogeology, oil recovery, and CO₂ sequestration. Accurate predictions of fluid behavior in these systems can improve process efficiency while mitigating environmental and health risks. Commercial simulators and open source software, such as the porousMultiphaseFoam repository based on the OpenFOAM framework, have been developed to model this type of problem. However, simulating heterogeneous porous media with heterogeneous porosity and permeability distributions poses significant numerical challenges. We introduce coupledMatrixFoam, an OpenFOAM-based solver designed for enhanced numerical stability and robustness. coupledMatrixFoam integrates the Eulerian multi-fluid formulation for phase fractions with Darcy's law for porous media flow, applying a fully implicit, block-coupled solution for pressure and phase fractions. The solver is based on foam-extend 5.0, leveraging the latest fvBlockMatrix developments to improve computational efficiency. This approach enables a significant increase in time step sizes, particularly in cases involving capillary pressure effects and other complex physical interactions. This work details the formulation, implementation and validation of coupledMatrixFoam, including comparisons with porousMultiphaseFoam that uses a segregated approach, to assess performance improvements. Additionally, a scalability analysis is conducted, demonstrating the solver's ability for high-performance computing (HPC) applications, which are essential for large-scale, real-world simulations.

Program summary

Program Title: coupledMatrixFoam

CPC Library link to program files: https://doi.org/10.17632/3d3xdh4x89.1

Developer's repository link: https://gitlab.com/wikki.brasil/porousmedia

Licensing provisions: GPLv3

Programming language: C++

Supplementary material: Available in the repository porousMedia.

Nature of problem: This software solves multiphase flow in porous media.

Solution method: Fully implicit solver based on the Euler-Euler multifluid formulation combined with Darcy's law developed in the OpenFOAM framework, that is based on the Finite Volume Method (FVM). Implicit coupling of phase fraction and pressure equations allowing significantly larger time steps. Complex physical phenomena are accounted for, including capillary pressure effects, gravitational forces, compressibility, and heterogeneous media properties. The porous medium is modeled as a stationary phase, and nonlinear terms are linearized using Taylor series expansions. Parallel computation is supported through OpenFOAM's standard domain decomposition approach.

Additional comments including restrictions and unusual features: Compatibility with foam-extend (OpenFOAM) version 5.0.

查看原文本刊更多论文

基于泡沫扩展的非均质多孔介质全隐式求解器

多孔介质中的多相流存在于许多工程应用中，包括水文地质、石油开采和二氧化碳封存。准确预测这些系统中的流体行为可以提高工艺效率，同时降低环境和健康风险。商业模拟器和开源软件，如基于OpenFOAM框架的porousmmultiphasefoam存储库，已经被开发出来对这类问题进行建模。然而，模拟具有非均质孔隙度和渗透率分布的非均质多孔介质存在显著的数值挑战。我们介绍了coupledMatrixFoam，一个基于openfoam的求解器，旨在增强数值稳定性和鲁棒性。coupledMatrixFoam将相分数的欧拉多流体公式与多孔介质流动的达西定律相结合，为压力和相分数应用了完全隐式的块耦合解决方案。求解器基于foam-extend 5.0，利用最新的fvBlockMatrix开发来提高计算效率。这种方法可以显著增加时间步长，特别是在涉及毛细管压力效应和其他复杂物理相互作用的情况下。这项工作详细介绍了coupledMatrixFoam的配方、实施和验证，包括与使用分离方法的porousMultiphaseFoam进行比较，以评估性能改进。此外，还进行了可扩展性分析，展示了求解器在高性能计算（HPC）应用程序中的能力，这对于大规模的真实世界模拟至关重要。程序摘要程序标题：coupledMatrixFoamCPC库链接到程序文件：https://doi.org/10.17632/3d3xdh4x89.1Developer's存储库链接：https://gitlab.com/wikki.brasil/porousmediaLicensing条款：gplv3编程语言：c++补充材料：在存储库中提供porousMedia。问题性质：该软件解决多孔介质中的多相流问题。求解方法：基于OpenFOAM框架中开发的基于有限体积法（FVM）的欧拉-欧拉多流体公式结合达西定律的全隐式求解器。相分数和压力方程的隐式耦合允许显著较大的时间步长。复杂的物理现象，包括毛细管压力效应，重力，压缩性和非均质介质性质。将多孔介质建模为固定相，非线性项采用泰勒级数展开进行线性化。OpenFOAM的标准领域分解方法支持并行计算。其他注释包括限制和不寻常的功能：与foam-extend (OpenFOAM) 5.0版本的兼容性。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

求助全文

约1分钟内获得全文求助全文

来源期刊

Computer Physics Communications 物理-计算机：跨学科应用

CiteScore

12.10

自引率

3.20%

发文量

287

审稿时长

5.3 months

期刊介绍： The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper. Computer Programs in Physics (CPiP) These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged. Computational Physics Papers (CP) These are research papers in, but are not limited to, the following themes across computational physics and related disciplines. mathematical and numerical methods and algorithms; computational models including those associated with the design, control and analysis of experiments; and algebraic computation. Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.