Interior point methods for an algebraic system involving complementarity equations for geomechanical fractures
Fuente:
arXiv
Enregistré dans:
| Auteur principal: | |
|---|---|
| Format: | Preprint |
| Publié: |
2024
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866912181505753088 |
|---|---|
| author | Hoang, Trung Hau |
| author_facet | Hoang, Trung Hau |
| contents | Many applications like subseismic fault modeling, fractured reservoir modeling and interpretation/validation of fault connectivity involve the solution to an elliptic boundary value problem in a background medium perturbed by the presence of cracks that take the form of one or many pieces of surface (with boundary). When the background medium can be considered as homogeneous, boundary integral equations appear as a method of choice for the numerical solution to fractures problems. With such an approach, the problem is reformulated as a fully non-local equation posed at the surface of cracks. Discretization of boundary integral resulting in the so-called Boundary Element Method (BEM) leads to densely populated matrices due to the full non-locality of the operators under consideration. After the discretization process, geologists are faced with a system of equations that turns out difficult to solve numerically. Many empirical algorithms have been proposed by geologists to solve this system of equations. Unfortunately, none of them is guaranteed to converge in theory (in particular when faults (fractures) intersect each other forming a geometrically highly irregular structure). In practice, none of them appears to be either robust or efficient. We investigate another approach, referred to as interior point methods, for which convergence can be ensured (even if faults are too close). Interior point methods have proved their efficiency in a wide variety of domains, most notably for linear programming. Here, even though we do not have any optimization problem, we can adapt ideas from interior point methods for the numerical resolution of the system considered. The numerical results obtained demonstrate computational efficiency and accuracy, highlighting the robustness and effectiveness of the implemented methods. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2501_04710 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Interior point methods for an algebraic system involving complementarity equations for geomechanical fractures Hoang, Trung Hau Optimization and Control Numerical Analysis G.1 Many applications like subseismic fault modeling, fractured reservoir modeling and interpretation/validation of fault connectivity involve the solution to an elliptic boundary value problem in a background medium perturbed by the presence of cracks that take the form of one or many pieces of surface (with boundary). When the background medium can be considered as homogeneous, boundary integral equations appear as a method of choice for the numerical solution to fractures problems. With such an approach, the problem is reformulated as a fully non-local equation posed at the surface of cracks. Discretization of boundary integral resulting in the so-called Boundary Element Method (BEM) leads to densely populated matrices due to the full non-locality of the operators under consideration. After the discretization process, geologists are faced with a system of equations that turns out difficult to solve numerically. Many empirical algorithms have been proposed by geologists to solve this system of equations. Unfortunately, none of them is guaranteed to converge in theory (in particular when faults (fractures) intersect each other forming a geometrically highly irregular structure). In practice, none of them appears to be either robust or efficient. We investigate another approach, referred to as interior point methods, for which convergence can be ensured (even if faults are too close). Interior point methods have proved their efficiency in a wide variety of domains, most notably for linear programming. Here, even though we do not have any optimization problem, we can adapt ideas from interior point methods for the numerical resolution of the system considered. The numerical results obtained demonstrate computational efficiency and accuracy, highlighting the robustness and effectiveness of the implemented methods. |
| title | Interior point methods for an algebraic system involving complementarity equations for geomechanical fractures |
| topic | Optimization and Control Numerical Analysis G.1 |
| url | https://arxiv.org/abs/2501.04710 |