Métricas de uso común para cuantificar la irregularidad de una malla triangular


9

Digamos que tiene una malla triangular en un plano plano. Esto se ha diseñado para resolver eventualmente algún problema en la mecánica, por ejemplo.

Una malla de triángulos equiláteros es la mejor en la medida en que las distancias entre los vértices y entre los centroides son las mismas en todas partes. Esto hace que las interpolaciones y el cálculo de gradientes sean una tarea fácil y precisa. Sin embargo, debido a restricciones y circunstancias, no siempre es posible trabajar en una malla de todos los triángulos equiláteros.

Entonces, las preguntas se refieren a una malla de elementos triangulares de forma arbitraria.

Sobre elementos de malla individuales . ¿Qué métricas se usan comúnmente para cuantificar la diferencia de un triángulo genérico de alguna forma equilátera ideal subyacente?

Sobre toda la malla . ¿Qué métricas se utilizan para cuantificar la irregularidad de una malla de triángulos arbitrarios en general? Estas métricas deberían indicar qué tan codificada está la malla.

Gracias por pensarlo.

Nota Todas las contribuciones de la comunidad de elementos finitos han sido muy apreciadas. Para esta pregunta, tenga en cuenta que el interés es cuantificar las diferencias puramente en la geometría (triángulos arbitrarios vs equiláteros). El efecto posterior sobre los errores de interpolación y condicionamiento están fuera del alcance. Por supuesto, estos pueden ser perspicaces y relevantes, complican el manejo matemático.


44
¿Has revisado esta pregunta ? Y de esa publicación: "¿Qué es un buen elemento finito?".
nicoguaro

3
Creo que la relación de áreas / radios entre el círculo y el círculo podría funcionar. La relación de valores propios de los ángulos jacobianos, mínimo y máximo, también.
nicoguaro

44
Uno de los documentos más famosos de Shewchuck cubre este tema en profundidad: ¿Qué es un buen elemento finito lineal?
Paul

@nicoguaro Gracias. No estoy especialmente interesado en FEM, sino en cuantificar la diferencia en la forma de los elementos. ¿Podría por favor explicar las proporciones de radios, por ejemplo? ¿Eso es independiente del tamaño? En otras palabras, se agradecerá si puede enumerar sus opciones en una respuesta para que cualquier otra persona desarrolle.
XavierStuvw

También puede mirar el ángulo mínimo en cualquiera de los elementos de malla. La idea es que esto quiere ser lo más grande posible
KyleW

Respuestas:


4

Como han dicho @Nicoguaro y @Paul en los comentarios a la publicación de preguntas, hay muchas maneras de hacer este tipo de cosas, y no estoy seguro de si existe un único "mejor" enfoque.


De un estudio de revisión de Jonathan Richard Shewchuck en Berkley, una respuesta es:

ingrese la descripción de la imagen aquí

Consulte el documento original (versión 31/12/2002) para obtener simbología, terminología, características especiales y posiblemente más (por ejemplo, tetraedros). El capítulo 6 trata sobre medidas de calidad. El documento vinculado es la versión extendida, y en la página web del JRS también hay una versión abreviada.


Personalmente, soy un fanático de la métrica de "longitud de volumen". Es un buen indicador escalar robusto de calidad simple (isotrópica) y es económico de calcular. En dos dimensiones:

una=4 433UNAmirmetros2

UNAmirmetrosuna=1una<0 0

Para evaluar la calidad de una triangulación no estructurada, es típico observar histogramas de tales métricas de calidad de elementos. Hay muchas implementaciones de este tipo de cosas por ahí, pero una recta de avance MATLABcódigo base mío es aquí .

Además de las puntuaciones de longitud de volumen, los histogramas de los ángulos de los elementos y el grado de vértice también se calculan de forma predeterminada.


¿Por qué eres fanático de esta métrica? ¿Fue bueno predecir la precisión de las simulaciones que hiciste con las mallas?
BrunoLevy

@BrunoLevy: Bueno, como una simple opción "predeterminada" para los símplex: se generaliza de forma robusta a dimensiones más altas, es barata de calcular, está numéricamente bien acondicionada, proporciona un indicador de "enredo" re. orientación, y es un simple indicador de "solo geometría", según la pregunta. ¿Es un buen indicador de la calidad de la simulación? Bueno, eso depende de lo que estés haciendo. Si estás interesado en las mallas isotrópicas, diría que sí. Configuraciones anisotrópicas altamente dependientes de la dirección, luego no, no directamente, aunque en tales casos todavía se puede usar después de una transformación de coordenadas adecuada.
Darren Engwirda

1
También es sencillo, por lo que puede ejecutarlo con una formulación arbitraria implícita de Lagrangian-Eulerian. Con un pequeño esfuerzo, puede generalizarlo a mallas anisotrópicas.
like

@likask: Sí, un buen punto: puede ser una buena función de costo para el suavizado y la optimización de mallas.
Darren Engwirda

He agregado un extracto del trabajo de Shewcuck que amplía el alcance de la respuesta de Darren. Esto resume varios comentarios también. Gracias a todos los colaboradores de esta publicación.
XavierStuvw

4

No creo que exista una respuesta a esta pregunta en general , porque todo depende del uso previsto para la malla. Por ejemplo, si está haciendo dinámica de fluidos computacional, es posible que desee tener una malla que sea extremadamente anisotrópica cerca de la capa límite. Ahora, si está haciendo electromagnetismo computacional, la mejor malla probablemente será completamente diferente.

Hay en la literatura muchas definiciones diferentes para un criterio de "calidad de malla". La mayoría de ellos favorecerá las mallas con triángulos que sean lo más equiláteros posible. También se puede mencionar la idea de maximizar el ángulo más pequeño (que se realiza mediante la triangulación de Delaunay para un conjunto fijo de puntos). Está justificado por el análisis de Jonathan Shewchuk mencionado en uno de los comentarios, que relaciona este ángulo con el número de condición de la matriz de rigidez para la ecuación de Laplace discretizada con elementos P1, pero nuevamente, dependiendo del uso previsto, la buena malla de alguien puede ser alguien Malla de los demás.

No creo que tenga sentido "cuantificar las diferencias puramente en la geometría (triángulos arbitrarios vs equiláteros)": antes de medir si los triángulos son equiláteros y decidir qué "desviación wrt equilateralidad" es la mejor, es necesario averiguar si "triángulos equiláteros" es lo que queremos, ¡y no siempre es así! Todo proviene de la "interpolación y condicionamiento" que usted menciona. Sí, como dijiste "complica el manejo matemático" pero sin él, no es posible hacer la diferencia entre los criterios objetivos para una aplicación dada y los criterios que no tienen ningún sentido.

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.