Articulo de referencia

Tipos de malla

Una malla es una representación de un dominio geométrico mayor mediante celdas discretas más pequeñas. Las mallas se utilizan comúnmente para calcular soluciones de ecuaciones d...

Una malla es una representación de un dominio geométrico mayor mediante celdas discretas más pequeñas. Las mallas se utilizan comúnmente para calcular soluciones de ecuaciones diferenciales parciales , generar gráficos por computadora y analizar datos geográficos y cartográficos. Una malla divide el espacio en elementos (o celdas o zonas ) sobre los cuales se pueden resolver las ecuaciones, lo que a su vez aproxima la solución en todo el dominio. Los límites de los elementos pueden estar restringidos a coincidir con los límites internos o externos del modelo. Los elementos de mayor calidad (mejor forma) poseen mejores propiedades numéricas, y lo que constituye un elemento "mejor" depende de las ecuaciones generales que rigen el modelo y de la solución particular para la instancia del modelo.

Formas celulares comunes

Bidimensional

Formas celulares bidimensionales básicas

Hay dos tipos de formas celulares bidimensionales que se utilizan comúnmente. Estas son el triángulo y el cuadrilátero .

Los elementos con baja capacidad de cálculo tendrán ángulos internos agudos , bordes cortos o ambas cosas.

Triángulo

Esta forma de celda consta de 3 lados y es uno de los tipos de malla más simples. Una malla de superficie triangular siempre es rápida y fácil de crear. Es más común en cuadrículas no estructuradas .

Cuadrilátero

Esta forma celular es la básica de cuatro lados, como se muestra en la figura. Es la más común en las redes estructuradas.

Los elementos cuadriláteros generalmente no pueden ser ni volverse cóncavos.

Tridimensional

Formas celulares tridimensionales básicas

Los elementos tridimensionales básicos son el tetraedro , la pirámide cuadrilátera , el prisma triangular y el hexaedro . Todos ellos tienen caras triangulares y cuadriláteras.

Los modelos bidimensionales extruidos pueden representarse completamente mediante prismas y hexaedros como triángulos y cuadriláteros extruidos.

En general, las caras cuadriláteras en 3 dimensiones pueden no ser perfectamente planas. Una cara cuadrilátera no plana puede considerarse un volumen tetraédrico delgado compartido por dos elementos vecinos.

Tetraedro

Un tetraedro tiene 4 vértices, 6 aristas y está delimitado por 4 caras triangulares. En la mayoría de los casos, se puede generar automáticamente una malla volumétrica tetraédrica.

Pirámide

Una pirámide de base cuadrilátera tiene 5 vértices, 8 aristas y está delimitada por 4 caras triangulares y 1 cara cuadrilátera. Estos elementos se utilizan eficazmente como elementos de transición entre elementos de caras cuadradas y triangulares, así como en mallas y rejillas híbridas.

Prisma triangular

Un prisma triangular tiene 6 vértices, 9 aristas y está delimitado por 2 caras triangulares y 3 cuadriláteras. La ventaja de este tipo de capa es que resuelve la capa límite de manera eficiente.

Hexaedro

Un cuboide , un cubo topológico , tiene 8 vértices, 12 aristas y 6 caras cuadriláteras, lo que lo convierte en un tipo de hexaedro . En el contexto de las mallas, un cuboide se suele llamar hexaedro , hex o ladrillo . [ 1 ] Para la misma cantidad de celdas, la precisión de las soluciones en mallas hexaédricas es la más alta.

Las zonas de pirámide y prisma triangular pueden considerarse computacionalmente como hexaedros degenerados, donde algunas aristas se han reducido a cero. También pueden representarse otras formas degeneradas de un hexaedro.

Celdas avanzadas (poliedro)

Un elemento poliédrico (dual) tiene cualquier número de vértices, aristas y caras. Generalmente requiere más operaciones de cálculo por celda debido al número de vecinos (normalmente 10). [ 2 ] Sin embargo, esto se compensa con la precisión del cálculo.

Clasificación de cuadrículas

Red estructurada
Red no estructurada

Redes estructuradas

Las mallas estructuradas se caracterizan por una conectividad regular. Los elementos posibles son cuadriláteros en 2D y hexaedros en 3D. Este modelo es altamente eficiente en cuanto al espacio, ya que las relaciones de vecindad se definen mediante la disposición del almacenamiento. Otras ventajas de las mallas estructuradas sobre las no estructuradas son una mejor convergencia y una mayor resolución. [ 3 ] [ 4 ] [ 5 ]

Redes no estructuradas

Una malla no estructurada se caracteriza por una conectividad irregular. No se puede representar fácilmente como una matriz bidimensional o tridimensional en la memoria de la computadora . Esto permite cualquier elemento posible que un solucionador pueda utilizar. En comparación con las mallas estructuradas, para las cuales las relaciones de vecindad son implícitas, este modelo puede ser muy ineficiente en cuanto a espacio, ya que requiere el almacenamiento explícito de dichas relaciones. Los requisitos de almacenamiento de una malla estructurada y de una malla no estructurada se encuentran dentro de un factor constante. Estas mallas suelen emplear triángulos en 2D y tetraedros en 3D. [ 6 ]

redes híbridas

Una malla híbrida contiene una mezcla de partes estructuradas y no estructuradas. Integra las mallas estructuradas y no estructuradas de manera eficiente. Las partes regulares de la geometría pueden tener mallas estructuradas, mientras que las partes complejas pueden tener mallas no estructuradas. Estas mallas pueden ser no conformes, lo que significa que las líneas de la malla no necesitan coincidir en los límites de los bloques. [ 7 ]

Calidad de la malla

Se considera que una malla tiene mayor calidad si permite calcular una solución más precisa con mayor rapidez. La precisión y la velocidad están en tensión. Disminuir el tamaño de la malla siempre aumenta la precisión, pero también incrementa el coste computacional.

La precisión depende tanto del error de discretización como del error de solución. En cuanto al error de discretización, una malla dada es una aproximación discreta del espacio y, por lo tanto, solo puede proporcionar una solución aproximada, incluso cuando las ecuaciones se resuelven de forma exacta. (En el trazado de rayos en gráficos por computadora , el número de rayos emitidos es otra fuente de error de discretización). En cuanto al error de solución, para las EDP se requieren muchas iteraciones sobre toda la malla. El cálculo se interrumpe prematuramente, antes de que las ecuaciones se resuelvan de forma exacta. La elección del tipo de elemento de malla afecta tanto al error de discretización como al error de solución.

La precisión depende tanto del número total de elementos como de la forma de cada elemento individual. La velocidad de cada iteración aumenta (linealmente) con el número de elementos, y el número de iteraciones necesarias depende del valor de la solución local y del gradiente en relación con la forma y el tamaño de los elementos locales.

Precisión de la solución

Una malla gruesa puede proporcionar una solución precisa si esta es constante, por lo que la precisión depende del caso particular. Se puede refinar selectivamente la malla en las zonas donde los gradientes de la solución son elevados, aumentando así la fidelidad en esas áreas. La precisión, incluyendo los valores interpolados dentro de un elemento, depende del tipo y la forma de dicho elemento.

Tasa de convergencia

Cada iteración reduce el error entre la solución calculada y la solución real. Una mayor velocidad de convergencia implica un menor error con menos iteraciones.

Una malla de baja calidad puede omitir características importantes, como la capa límite del flujo de fluidos. El error de discretización será grande y la tasa de convergencia se verá afectada; es posible que la solución no converja en absoluto.

Independencia de la red eléctrica

Una solución se considera independiente de la malla si el error de discretización y de solución es suficientemente pequeño tras un número suficiente de iteraciones. Esto es fundamental para obtener resultados comparativos. Un estudio de convergencia de malla consiste en refinar los elementos y comparar las soluciones refinadas con las soluciones iniciales. Si un mayor refinamiento (u otros cambios) no modifica significativamente la solución, la malla es una «malla independiente».

Decidir el tipo de malla

Asimetría basada en el volumen equilátero

Si la precisión es primordial, la malla hexaédrica es la más recomendable. La densidad de la malla debe ser suficientemente alta para capturar todas las características del flujo, pero no tan alta como para capturar detalles innecesarios, sobrecargando así la CPU y consumiendo más tiempo. Cuando hay una pared, la malla adyacente es lo suficientemente fina para resolver el flujo de la capa límite, y generalmente se prefieren las celdas cuadriláteras, hexagonales y prismáticas a las triangulares, tetraédricas y piramidales. Las celdas cuadriláteras y hexagonales se pueden estirar cuando el flujo está completamente desarrollado y es unidimensional.

Representa la asimetría de un cuadrilátero.

En función de la asimetría, la suavidad y la relación de aspecto, se puede determinar la idoneidad de la malla. [ 8 ]

Oblicuidad

La asimetría de una cuadrícula es un indicador preciso de su calidad y adecuación. Una asimetría elevada compromete la precisión de las regiones interpoladas. Existen tres métodos para determinar la asimetría de una cuadrícula.

Basado en el volumen equilátero

Este método solo es aplicable a triángulos y tetraedros, y es el método predeterminado.

 Oblicuidad = tamaño celular óptimo - tamaño celular tamaño celular óptimo{\displaystyle {\text{ Asimetría }}={\frac {\text{ tamaño óptimo de celda - tamaño de celda }}{\text{tamaño óptimo de celda}}}}
Cambio de salto suave y amplio

Basado en la desviación del ángulo equilátero normalizado

Este método se aplica a todas las formas celulares y faciales, y casi siempre se utiliza para prismas y pirámides.

 Asimetría (para un cuadrilátero) =máximo[θmáximo9090,90θmin90]{\displaystyle {\text{ Asimetría (para un cuadrilátero) }}=\max {\left[{\frac {\theta _{\text{máx}}-90}{90}},{\frac {90-\theta _{\text{mín}}}{90}}\right]}}

Sesgo equiangular

Otra medida común de calidad se basa en la asimetría equiangular.

 Sesgo equiángulo =máximo[θmáximoθmi180θmi,θmiθminθmi]{\displaystyle {\text{ Equiangle Skew }}=\max {\left[{\frac {\theta _{\text{max}}-\theta _{e}}{180-\theta _{e}}},{\frac {\theta _{e}-\theta _{\text{min}}}{\theta _{e}}}\right]}}

dónde:

  • θmáximo{\displaystyle \theta _{\text{max}}\,}es el ángulo más grande en una cara o célula,
  • θmin{\displaystyle \theta _{\text{min}}\,}es el ángulo más pequeño en una cara o célula,
  • θmi{\displaystyle \theta _{e}\,}es el ángulo para cara o celda equiangular, es decir, 60 para un triángulo y 90 para un cuadrado.

Un coeficiente de asimetría de 0 es el mejor posible, y un coeficiente de asimetría de 1 casi nunca es preferible. Para celdas hexagonales y cuadriláteras, el coeficiente de asimetría no debe superar 0,85 para obtener una solución bastante precisa.

Representa los cambios en la relación de aspecto.

Para las celdas triangulares, la asimetría no debe exceder 0,85 y para las celdas cuadriláteras, la asimetría no debe exceder 0,9.

Suavidad

El cambio de tamaño también debe ser gradual. No debe haber saltos bruscos en el tamaño de la celda, ya que esto podría provocar resultados erróneos en los nodos cercanos.

Relación de aspecto

Es la relación entre el lado más largo y el más corto de una celda. Idealmente, debería ser igual a 1 para garantizar los mejores resultados. Para flujos multidimensionales , debería ser cercana a uno. Además, las variaciones locales en el tamaño de la celda deben ser mínimas; es decir, los tamaños de celdas adyacentes no deberían variar en más del 20 %. Una relación de aspecto elevada puede generar un error de interpolación de magnitud inaceptable.

Generación y mejora de mallas

Véase también generación de mallas y principios de generación de cuadrículas . En dos dimensiones, el volteo y el suavizado son herramientas poderosas para adaptar una malla deficiente a una buena malla. El volteo consiste en combinar dos triángulos para formar un cuadrilátero y luego dividir el cuadrilátero en la dirección opuesta para producir dos nuevos triángulos. El volteo se utiliza para mejorar medidas de calidad de un triángulo, como la asimetría. El suavizado de mallas mejora las formas de los elementos y la calidad general de la malla ajustando la ubicación de los vértices de la malla. En el suavizado de mallas, se conservan características esenciales como el patrón no nulo del sistema lineal , ya que la topología de la malla permanece invariante. El suavizado laplaciano es la técnica de suavizado más utilizada.

Véase también

Referencias

  1. "Elementos del hexaedro" (PDF) . Archivado del original (PDF) el 24 de febrero de 2015. Consultado el 13 de abril de 2015 .
  2. "Copia archivada" (PDF) . Archivado del original (PDF) el 06-12-2013 . Recuperado el 10-01-2018 .{{cite web}}: CS1 mantenimiento: copia archivada como título ( enlace )
  3. "Calidad y control: dos razones por las que las redes estructuradas no van a desaparecer" .
  4. Castillo, JE (1991), "Aspectos matemáticos de la generación de cuadrículas", Sociedad de Matemáticas Industriales y Aplicadas, Filadelfia
  5. George, PL (1991), Generación automática de mallas
  6. Mavriplis, DJ (1996), "Generación de mallas y adaptabilidad para geometrías y flujos complejos", Manual de mecánica de fluidos computacional
  7. Bern, Marshall; Plassmann, Paul (2000), "Generación de mallas", Manual de geometría computacional. Elsevier Science
  8. "Meshing, Lección 7" . Andre Bakker . Consultado el 10 de noviembre de 2012 .
  • Logotipo de Wikimedia CommonsContenido multimedia relacionado con los tipos de malla en Wikimedia Commons.