Subgrid scale stabilized finite elements for low speed flows

Abstract

La descripción del flujo de fluidos involucra la solución de las ecuaciones de Navier-Stokes compresible, un problema muy complejo cuya estructura matemática no es del todo comprendida. Por lo tanto, mediante análisis asintótico, se pueden derivar modelos simplificados bajo ciertas hipótesis sobre el problema hechas en términos de parámetros adimensionales que miden la importancia relativa de los diferentes procesos físicos. Los flujos a baja velocidad se pueden describir por diferentes modelos que incluyen las ecuaciones de Navier Stokes incompresible cuya matemática es mucho mas conocida. Sin embargo, algunos flujos importantes no se pueden considerar incompresibles debido a la presencia de efectos térmicos. En esta clase de problemas se pueden derivar otra clase de ecuaciones simplificadas: las ecuaciones de Boussinesq y las ecuaciones de bajo numero de Mach.<br/><br/>La complejidad de estos problemas matemáticos hace que su solución numérica sea muy difícil. En estos problemas el método de los elementos finitos es inestable, lo que en la práctica implica soluciones numéricas que presentan oscilaciones nodo a nodo de naturaleza no física. En las ecuaciones de Navier Stokes incompresible, dos fuentes bien conocidas de inestabilidad son la condición de incompresibilidad y la presencia del término convectivo. Muchas técnicas de estabilización utilizadas hoy en día se basan en la separación de escalas, descomponiendo la incógnita en una parte gruesa inducida por la discretización del domino y una parte fina de subescala. Modelar la subescala y su influencia conduce a un problema modificado para la escala gruesa que resulta estable.<br/><br/>Aunque las técnicas de estabilización son ampliamente utilizadas hoy en día, importantes problemas permanecen abiertos. Contribuyendo a su comprensión, en este trabajo se analizan varios aspectos del modelado de las subescalas. Para problemas escalares de segundo orden, se encuentra la dependencia de la subescala con el tamaño de la malla en el caso general de mallas anisótropas. Estas ideas son extendidas a sistemas de ecuaciones para considerar el problema de Oseen. También se analiza el modelado de las subescalas en problemas transitorios, obteniendo un mejor esquema de integración temporal para el problema de escala gruesa. Para considerar flujos a baja velocidad, se presenta la extensión de estas técnicas a problemas no lineales acoplados, lo que esta íntimamente relacionado con el problema del modelado de la turbulencia, que es un tema en si mismo.<br/><br/>Los flujos acoplados térmicamente, aparte del interés intrínseco que merecen, son importantes desde un punto de vista ingenieril. Una solución precisa del problema de flujo es necesaria para definir las cargas térmicas sobre las estructuras, que en muchos casos responden fuertemente, haciendo el problema acoplado. Esta clase de problemas, que motivaron este trabajo, incluyen la respuesta estructural en el caso de un incendio.

A general description of a fluid flow involves the solution of the compressible Navier-Stokes equations, a very complex problem whose mathematical structure is not well understood. Therefore, simplified models can be derived by asymptotic analysis under some assumptions on the problem, made in terms of dimensionless parameters that measure the relative importance of different physical processes. Low speed flows can be described by several models including the incompressible Navier Stokes equations whose mathematical structure is much better understood. However many important flows cannot be considered as incompressible, even at low speed, due to the presence of thermal effects. In such kind of problems another class of simplified equations can be derived: the Boussinesq equations and the Low Mach number equations.<br/>The complexity of these mathematical problems makes their numerical solution very difficult. For these problems the standard finite element method is unstable, what in practice means that node to node oscillations of non physical nature may appear in the numerical solution. In the incompressible Navier Stokes equations, two well known sources of numerical instabilities are the incompressibility constraint and the presence of the convective terms. Many stabilization techniques used nowadays are based on scale separation, splitting the unknown into a coarse part induced by the discretization of the domain and a fine subgrid part. The modelling of the subgrid scale and its influence leads to a modified coarse scale problem that now can be shown to be stable. <br/>Although stabilization techniques are nowadays widely used, important problems remain open. Contributing to their understanding, several aspects of the subgrid scale modelling are analyzed in this work. For second order scalar problems, the dependence of the subgrid scale on the mesh size, in the general anisotropic case, is clarified. These ideas are extended to systems of equations to consider the Oseen problem. The modelling of the subgrid scales in transient problems is also analyzed, leading to an improved time discretization scheme for the coarse scale problem. To consider low speed flow models, the extension of these techniques to nonlinear and coupled problems is presented, something that is intimately related to the problem of turbulence modelling, which a entire subject on its own right. <br/>Thermally coupled flow problems, despite the intrinsic interest they deserve, are important from an engineering point of view. An accurate solution of a flow problem is needed to define thermal loads on structures which, in many cases have a strong response, making the problem coupled. This kind of problems, that motivated this work, include the problem of a structural response in the case of fires.