4.4 Article

On the accuracy of a nonlinear finite volume method for the solution of diffusion problems using different interpolations strategies

Journal

Publisher

WILEY
DOI: 10.1002/fld.3850

Keywords

diffusion problems; heterogeneous and anisotropic media; non linear finite volume; cell-centered; unstructured meshes; interpolation strategies

Funding

  1. Fundacao de Amparo a Ciencia e Tecnologia do Estado de Pernambuco (FACEPE)
  2. Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq)
  3. CENPES-PETROBRAS (SIGER - Rede Tematica de Simulacao e Gerenciamento de Reservatorios)
  4. Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior (CAPES)
  5. Agencia Nacional do Petroleo (ANP)

Ask authors/readers for more resources

In this paper, we consider a nonlinear finite volume method to solve the steady-state diffusion equation in nonhomogeneous and non-isotropic media. The method is nonlinear even if the original problem is linear. In its original form, the scheme is monotone, because the coefficient matrix is monotone under certain assumptions and, as a consequence, whenever the analytic operator demands, it preserves the positivity of numerical solutions. On the other hand, the scheme is unable to reproduce piecewise linear solutions exactly. In order to recover this interesting feature, we use two different interpolation strategies. In this case, even though we are unable to prove monotonicity, we show some numerical evidences that the combined method has an improved behavior, producing second order accurate solutions, even for nonhomogeneous and strongly anisotropic media. Copyright (c) 2013 John Wiley & Sons, Ltd.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.4
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Engineering, Multidisciplinary

A very high-order flux reconstruction approach coupled to the MPFA-QL finite volume method for the numerical simulation of oil-water flows in 2D petroleum reservoirs

G. Galindez-Ramirez, F. R. L. Contreras, D. K. E. Carvalho, P. R. M. Lyra

Summary: In this paper, a new high-order numerical methodology based on unstructured quadrilateral meshes is proposed for modeling oil-water displacements in highly heterogeneous and anisotropic petroleum reservoirs. The method achieves comparable accuracy to lower-order counterparts while reducing computational cost, and its accuracy, efficiency, and robustness are demonstrated through representative examples.

APPLIED MATHEMATICAL MODELLING (2022)

Article Engineering, Multidisciplinary

An algebraic multiscale solver for the simulation of two-phase flow in heterogeneous and anisotropic porous media using general unstructured grids (AMS-U)

Artur Castiel Reis de Souza, Darlan Karlo Elisiario de Carvalho, Jose Cicero Araujo dos Santos, Ramiro Brito Willmersdorf, Paulo Roberto Maciel Lyra, Michael G. Edwards

Summary: This paper introduces a new multiscale finite volume framework for simulating multi-phase flow in heterogeneous and anisotropic porous media. The framework allows the use of geophysical grid defined properties on high-definition grids, addressing the issues of basis function leakage and mass conservation. The accuracy of the framework is validated through comparisons with direct simulations on fine-scale, and the results demonstrate its ability to produce well-resolved solutions for complex geological formations found in petroleum reservoir problems.

APPLIED MATHEMATICAL MODELLING (2022)

Article Engineering, Civil

Assessment by finite element modeling of pipelines with corrosion defects based on River-Bottom Profile model

Rodolfo M. S. Cabral, Adriano D. M. Ferreira, Julio T. Pimentel, Marco A. F. da Silva Cabral, Paulo R. M. Lyra, Silvana M. B. Afonso, Ramiro B. Willmersdorf

Summary: This paper presents a new methodology for predicting pipelines failure pressure using axisymmetric modeling and the River Bottom Profile (RBP) of corrosion defects. This practical tool allows pipeline engineers in the oil and gas industry to quickly and accurately assess pipeline integrity.

ENGINEERING STRUCTURES (2022)

Letter Computer Science, Interdisciplinary Applications

An algebraic dynamic multilevel and multiscale method with non-uniform mesh resolution and adaptive algebraic multiscale solver operator for the simulation of two-phase flows in highly heterogeneous petroleum reservoirs

Jose Cicero Araujo dos Santos, Paulo Roberto Maciel Lyra, Joao Paulo Rodrigues de Andrade, Artur Castiel Reis de Souza, Ricardo Jorge Morais de Lira Filho, Darlan Karlo Elisiario de Carvalho

Summary: This paper proposes an adaptive flow-based dual volume agglomeration strategy for correcting non-physical terms in the coarse transmissibility matrix of classical multiscale finite volume (MsFV) methods. It also presents a framework to handle non-uniform levels at each coarse control volume, in order to reduce the size of coarse scale matrices. The proposed methodologies are applied to approximate pressure solutions in an Implicit Pressure Explicit Saturation (IMPES) strategy, and their accuracy and efficiency are demonstrated through testing with challenging benchmark problems.

JOURNAL OF COMPUTATIONAL PHYSICS (2022)

Article Computer Science, Interdisciplinary Applications

Numerical simulation of oil-water displacements in naturally fractured reservoirs using a hybrid grid MPFA-streamline (SLS-HyG) formulation

J. C. Teixeira, D. K. E. Carvalho, T. M. Cavalcante, K. C. L. Angelim, P. R. M. Lyra

Summary: This paper presents a numerical formulation for simulating two-phase flow in naturally fractured oil reservoirs using unstructured quadrilateral meshes and Hybrid-Grid MPFA-Streamline method. It can handle irregular polygonal grids and decouples transport equations using streamline-based method to solve advective saturation problems.

COMPUTERS AND GEOTECHNICS (2022)

Article Mathematics, Applied

A Multipoint Flux Approximation with a Diamond Stencil and a Non-Linear Defect Correction Strategy for the Numerical Solution of Steady State Diffusion Problems in Heterogeneous and Anisotropic Media Satisfying the Discrete Maximum Principle

T. M. Cavalcante, R. J. M. Lira Filho, A. C. R. Souza, D. K. E. Carvalho, P. R. M. Lyra

Summary: In this paper, we solve the steady state diffusion equation in 3D domains using the MPFA-DNL method, which guarantees the Discrete Maximum Principle by introducing a non-linear defect correction strategy. The method is locally conservative and capable of handling arbitrary anisotropic diffusion tensors and unstructured meshes while maintaining second order convergence rates.

JOURNAL OF SCIENTIFIC COMPUTING (2022)

Review Water Resources

A high-resolution multidimensional finite volume scheme coupled to a nonlinear two-point flux approximation method for the numerical simulation of groundwater contaminant transport using unstructured 2D meshes

Fernando R. L. Contreras, Uewerton A. O. Vaz, Gustavo L. S. S. Pacheco, Alessandro R. E. Antunes, Paulo R. M. Lyra, Darlan K. E. Carvalho

Summary: This paper proposes a novel full finite volume method to solve the advection-dispersion transport equation, combining various numerical methods to ensure accuracy and robustness of the numerical solution. Numerical experiments show that the method can provide accurate solutions when simulating groundwater processes with complex physical and geological properties.

ADVANCES IN WATER RESOURCES (2023)

Article Energy & Fuels

A very high order Flux Reconstruction (FR) method for the numerical simulation of 1-D compositional fluid flow model in petroleum reservoirs

M. E. S. Galindo, I. V. Lacerda, G. Galindez-Ramirez, P. R. M. Lyra, D. K. E. Carvalho

Summary: Compositional reservoir simulation is a crucial tool for modeling fluid flow in complex petroleum reservoirs, especially for volatile reservoir fluids or those involving enhanced oil recovery. Simple black-oil models are inadequate in these cases. The compositional model involves solving a large system of partial differential equations that describe mass conservation, Darcy's law, and fugacity constraints. However, the complexity and computational demands of the compositional problems are high due to the large number of equations and constraints. To improve accuracy and reduce computational costs, higher-order methods can be used to approximate the advective flux terms in the reservoir's multicomponent transport.

GEOENERGY SCIENCE AND ENGINEERING (2023)

Article Computer Science, Interdisciplinary Applications

Coupling finite elements of class C1 on composite curved meshes for second order elliptic problems

Ashish Bhole, Herve Guillard, Boniface Nkonga, Francesca Rapetti

Summary: Finite elements of class C-1 are used for computing magnetohydrodynamics instabilities in tokamak plasmas, and isoparametric approximations are employed to align the mesh with the magnetic field line. This numerical framework helps in understanding the operation of existing devices and predicting optimal strategies for the international ITER tokamak. However, a mesh-aligned isoparametric representation encounters issues near critical points of the magnetic field, which can be addressed by combining aligned and unaligned meshes.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS (2024)

Article Computer Science, Interdisciplinary Applications

Efficient hyperbolic-parabolic models on multi-dimensional unbounded domains using an extended DG approach

Federico Vismara, Tommaso Benacchio

Summary: This paper introduces a method for solving hyperbolic-parabolic problems on multidimensional semi-infinite domains. By dividing the computational domain into bounded and unbounded subdomains and coupling them using numerical fluxes at the interface, accurate numerical solutions are obtained. In addition, computational cost can be reduced by tuning the parameters of the basis functions.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS (2024)

Article Computer Science, Interdisciplinary Applications

Extension of moving particle simulation by introducing rotational degrees of freedom for dilute fiber suspensions

Keigo Enomoto, Takato Ishida, Yuya Doi, Takashi Uneyama, Yuichi Masubuchi

Summary: We have developed a novel Moving Particle Simulation (MPS) method to accurately reproduce the motion of fibers in sheared liquids. By introducing the micropolar fluid model, we address the issue of fibers being aligned with the flow direction in conventional MPS simulations. Our method is capable of accurately reproducing the fiber motion predicted by Jeffery's theory.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS (2024)