EBE Método EBE de los elementos finitos para sistemas de gran tamaño y su aplicación en la Física Digital de Rocas

• Benavides, Francisco José ^[1]
  1. [1] Universidad Nacional de Costa Rica

     Universidad Nacional de Costa Rica

     Heredia, Costa Rica

     
• Localización: Ingeniería: Revista de la Universidad de Costa Rica, ISSN-e 2215-2652, ISSN 1409-2441, Vol. 29, Nº. 2, 2019, págs. 44-58
• Idioma: español
• Títulos paralelos:
  • EBE finite element methdology for large scale systems applied to digital rock physics
• Enlaces
  • Texto completo
• Resumen
  • español

    El método de los elementos finitos (FEM) es una técnica numérica ampliamente utilizada en física e ingeniería, que aproxima soluciones de ecuaciones diferenciales en derivadas parciales en dominios arbitrarios. En general, el paso final de esta técnica consiste en la resolución de un sistema lineal de la forma  en el que la matriz  es dispersa y cuyo ancho de banda (número de filas no nulas en la matriz) depende del soporte de las funciones de forma definidas en los elementos.  Cuando este número de elementos es muy grande, inclusive las representaciones en estructuras de datos eficientes para matrices dispersas son capaces de consumir en su totalidad la memoria de un computador poderoso. En este artículo se describe el método EBE (a.k.a. Element By Element), un truco computacional diseñado para implementar el método FEM sin utilizar la representación explicita de esta matriz. Se describe también un contexto en el que es necesaria esta implementación, en el área de la petrofísica digital de rocas, para estimar las propiedades elásticas de una muestra de roca a partir de imágenes tomográficas.

  • English

    The finite element method (FEM) is a numerical technique that estimates solutions of partial differential equations on arbitrary domains. It has been widely used to solve problems in physics and engineering. In general, the final step of this technique consists of a linear system of equations  in which the matrix  is sparse and its bandwidth depends on the finite element shape functions support. When the number of elements is large, even the efficient data structure sparse matrix representations can consume the entire computer memory. In this article, we describe a technique to solve these large-scale problems without explicitly representing this matrix. This computational trick is known as EBE (Element By Element). We also describe an application in which such kind of implementation is necessary, in the field of digital rock physics, to estimate the elastic coefficient of rock samples using micro-tomographic images.