A LATIN computational strategy for multiphysics problems: application to poroelasticity (original) (raw)
Related papers
A multi-time-scale strategy for multiphysics problems: application to poroelasticity
2003
Usually, multiphysics phenomena and coupled-field problems lead to computationally intensive structural analysis. Strategies to keep these problems computationally affordable are of special interest. For coupled fluidstructure problems, for instance, partitioned procedures and staggered algorithms are often preferred to direct analysis. In a previous paper, a new strategy derived from the LArge Time INcrement (LATIN) method was described. This strategy was applied to the consolidation of saturated porous soils, which is a highly coupled fluid-solid problem. The feasibility of the method and the comparison of its performance with that of a standard partitioning scheme (the so-called ISPP method) was presented. Here, we go one step further and use the LATIN method to take into account the different time scales which usually arise from the different physics. We propose a multi-timescale strategy which improves the existing method.
A computational strategy suitable for multiphysics problems
2005
Multiphysics phenomena and coupled-field problems usually lead to computationally intensive structural analyses. Strategies to keep these problems computationally affordable are of special interest. For coupled fluid-structure problems for instance, partitioned procedures and staggered algorithms are often preferred to direct analysis (also called the monolithic approach), from a computational efficiency point of view. Recently, a mixed domain decomposition method has been designed for parallel computing environments, and a multi-level approach embedding a homogenization procedure makes it suitable for highlyheterogeneous problems. From the generalization of the concept of geometric interfaces between substructures to an interface between different physics, the Large Time INcrement method (LATIN) allows building an approach suited for solving coupled multiphysics problems. The proposed application concerns the consolidation of porous saturated soil, i.e. a coupled fluid-solid problem in the domain. The feasability of the method and its performance comparison with a standard partitioning scheme (the so-called ISPP) has been presented in a previous paper. As an improvement, the further step is to take into account different time scales arising from multiphysics problem. Thus, the present paper proposes a time multiscale strategy.
A computational strategy for poroelastic problems with a time interface between coupled physics
International Journal for Numerical Methods in Engineering, 2008
This paper deals with a computational strategy suitable for the simulation of multiphysics problems and based on the LArge Time INcrement (LATIN) method. One of the main issues in the design of advanced tools for the simulation of such problems is to take into account the different time and space scales that usually arise with the different physics. Here, we focus on using different time discretizations for each physics by introducing an interface with its own discretization. The proposed application concerns the simulation of a 2-physics problem: the fluid-structure interaction in porous media.
In this paper, we solve a problem describing the mechanical de- formations of a porous medium in the presence of a monophasic linear flow or a two phase nonlinear flow with the purpose of modelizing subsidence of hydrocarbon reservoirs. An essential characteristics of this problem is that the mechanical deformation and the flow have different time scales. In petroleum industry, one uses different very efficient simulators for theflow problem and the mechanical deformations, which enables to handle complex models. Therefore it is necessary to be able to combine as efficiently as possiblethe exploitation of these simulators. We propose two alternative splitting approaches. The first one is the staggered algorithm used by engineers, which amounts to a Gauss- Seidel method in the one phase linear case. The second approach is based upon the preconditioned conjugate gradient method. We use a numerical multi-scale method in both of these algorithms. We compare these two approaches and we s...
On a decoupled algorithm for poroelasticity
In this paper we present a reformulation of the poroelasticity system and a decoupled algorithm for the numerical solution. This reformulation enables us to construct a highly efficient multigrid method, confirmed by a realistic experiment.
Numerical convergence study of iterative coupling for coupled flow and geomechanics
In this paper we consider algorithms for modeling complex processes in porous media that include fluid and structure interactions. Numerous field applications would benefit from a better understanding and integration of porous flow and solid deformation. Important applications in environmental and petroleum engineering include carbon sequestration, surface subsidence, pore collapse, cavity generation, hydraulic fracturing, thermal fracturing, wellbore collapse, sand production, fault activation, and waste disposal, while similar issues arise in biosciences and chemical sciences as well.
A Fully-Conservative Finite Volume Formulation for Coupled Poro- Elastic Problems
2018
Solid mechanics is a research field that deals with the mechanical behavior of a wide variety of materials undergoing external loads. Among the various types of solids, porous materials, for instance, can be found in applications such as soil and rock mechanics, biomechanics, ceramics, etc. These applications are studied in the field of poromechanics, which is a specific branch of the solid mechanics that considers all types of porous materials. An important characteristic of such materials is that they contain a network of interconnected pore channels saturated with a fluid. In most situations the mechanical behavior of the porous matrix and the fluid flow through the pore channels are two tightly coupled phenomena interfering with each other. When the fluid moves from one region to another in the porous matrix it changes the pressure field inside the pore channels, which is perceived by the porous matrix as a force imbalance. As a consequence, the porous matrix tends to deform in ...
Advances in Geo-Energy Research, 2021
Research on the scientific and engineering problems of porous media has drawn increasing attention in recent years. Digital core analysis technology has been rapidly developed in many fields, such as hydrocarbon exploration and development, hydrology, medicine, materials and subsurface geofluids. In summary, science and engineering research in porous media is a complex problem involving multiple fields. In order to encourage communication and collaboration in porous media research using digital core technology in different industries, the 5 th International Conference on Digital Core Analysis & the Workshop on Multiscale Numerical and Experimental Approaches for Multiphysics Problems in Porous Media was held in Qingdao from April 18 to 20, 2021. The workshop was jointly organized by the China InterPore Chapter, the Research Center of Multiphase Flow in Porous Media at the China University of Petroleum (East China) and the University of Aberdeen with financial support from the National Sciences Foundation of China and the British Council. Due to the current pandemic, a hybrid meeting was held (participants in China met in Qingdao, while other participants joined the meeting online), attracting more than 150 participants from around the world, and the latest multi-scale simulation and experimental methods to study multi-field coupling problems in complex porous media were presented.
A fully coupled multiphase flow and geomechanics solver for highly heterogeneous porous media
Journal of Computational and Applied Mathematics, 2013
This paper introduces a fully coupled multiphase flow and geomechanics solver that can be applied to modeling highly heterogeneous porous media. Multiphase flow in deformable porous media is a multiphysics problem that considers the flow physics and rock physics simultaneously. To model this problem, the multiphase flow equations and geomechanical equilibrium equation must be tightly coupled. Conventional finite element modeling of coupled flow and geomechanics does not conserve mass locally since it uses continuous basis functions. Mixed finite element discretization that satisfies local mass conservation of the flow equation can be a good solution for this problem. In addition, the stabilized finite element method for discretizing the saturation equation minimizes numerical diffusion and provides better resolution of saturation solution.
We develop a computational model to study the interaction of a fluid with a poroelastic material. The coupling of Stokes and Biot equations represents a prototype problem for these phenomena, which feature multiple facets. On one hand, it shares common traits with fluid-structure interaction. On the other hand it resembles the Stokes-Darcy coupling. For these reasons, the numerical simulation of the Stokes-Biot coupled system is a challenging task. The need of large memory storage and the difficulty to characterize appropriate solvers and related preconditioners for the equations at hand are typical shortcomings of classical discretization methods applied to this problem, such as the finite element method for spatial discretization and finite differences for time stepping. The application of loosely coupled time advancing schemes mitigates these issues, because it allows to solve each equation of the system independently with respect to the others, at each time step. In this work, we develop and thoroughly analyze a loosely coupled scheme for Stokes-Biot equations. The scheme is based on Nitsche's method for enforcing interface conditions. Once the interface operators corresponding to the interface conditions have been defined, time lagging allows us to build up a loosely coupled scheme with good stability properties. The stability of the scheme is guaranteed provided that appropriate stabilization operators are introduced into the variational formulation of each subproblem. The error of the resulting method is also analyzed, showing that splitting the equations pollutes the optimal approximation properties of the underlying discretization schemes. In order to restore good approximation properties, while maintaining the computational efficiency of the loosely coupled approach, we consider the application of the loosely coupled scheme as a preconditioner for the monolithic approach. Both theoretical insight and numerical results confirm that this is a promising way to develop efficient solvers for the problem at hand.