Un método de elementos finitos incondicionalmente estable en norma uniforme para resolver las ecuaciones de Euler 2D
Visualitza/Obre
Estadístiques de LA Referencia / Recolecta
Inclou dades d'ús des de 2022
Cita com:
hdl:2099/7355
Tipus de documentArticle
Data publicació1993
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
Abstract
Este trabajo presenta dos algoritmos de tipo transporte e interpolación con Elementos
Finitos para la resolución numérica de las ecuaciones de Euler para flujos incompresibles bidimensionales en todo el espacio IR'. En la primera versión de nuestro algoritmo, la vorticidad es discretizada mediante elementos finitos triangulares de primer grado. En la segunda, mediante elementos finitos triangulares de segundo grado. La velocidad se obtiene en ambos casos calculando exactamente el producto de convolución del núcleo de Biot-Savart con una aproximación lineal a trozos sobre cada triángulo de la vorticidad discreta. Se prueba que el primer algoritmo es incondicionalmente estable y convergente con una precisión de primer orden, en norma uniforme. Sin embargo, en la práctica la precisión alcanzada resulta escasa, debido a la difusión numérica introducida en el paso de interpolación. En el caso del segundo
algoritmo, los ensayos numéricos muestran un notable incremento de la precisión, incluso para tiempos largos. Sin embargo, en este caso el algoritmo deja de ser estable en norma uniforme. We introduce two Finite Element transport-interpolation algorithms to solve the twodimensional Euler equations in the wole R2. In the first of these algorithms, the vorticity is discretized with triangular finite elements of degree one, and of degree two in the second one. The velocity is computed by convolution of the Biot-Savart kernel with a piecewise affine interpolate of the vorticity. We prove that the first algorithm is unconditionally stable in
uniform norm, with first order accuracy. However, in practice its precision is rather low, due to the numerical diffusion introduced in the interpolation step. The second algorithm is shown numerically to produce a remarkable increase of precision, even for long integration times. However, in this case the algorithm is no longer stable in uniform norm.
ISSN1886-158X
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
Article01.pdf | 1,217Mb | Visualitza/Obre |