Anteproyecto Darwin PDF
Anteproyecto Darwin PDF
Anteproyecto Darwin PDF
PROGRAMA DE MATEMTICA
MAESTRA EN MATEMTICA
Estudiante:
Darwin Lopez Aguas
Director:
Ph.D Miguel Antonio Caro Candezano
Abril 2016
Resumen
este trabajo realizamos un estudio numrico de las ecuaciones de Navier-Stokes para el problema liddriven cavity en el caso bidimensional implementando esquemas upwind de alta resolucin. El sistema
de ecuaciones ser discretizado a travs del mtodo de diferencias finitas aplicando el mtodo de
proyeccin de velocidad para forzar la incompresibilidad del flujo utilizando operadores numricos para
el caso bidimensional. En cuanto al tratamiento discreto de los trminos convectivos del sistema utilizaremos
inicialmente el esquema FOU y posteriormente esquemas upwind de alta resolucin basados en criterios de
limitacin TVD-CBC garantizando estabilidad en las soluciones numricas y poder establecer mayores nmeros
de Reynolds e incrementar la estabiliadad y precisin en los resultados numricos y as comparar los resultados
obtenidos con los de la literatura especializada.
N
Palabras Claves: Navier-stokes; lid-driven cavity; esquemas de alta resolucin; diferencias finitas; TVD-CBC.
Contenido
Resumen
Contenido
iv
Introduccin
Objetivos
vii
1 Preliminares
. . . . . . . . .
16
16
18
19
20
21
2 Resultados Numricos
23
23
iv
REFERENCIAS BIBLIOGRFICAS
CONTENIDO
27
Introduccin
la mecnica de fluidos que utiliza mtodos numricos y algoritmos para resolver y analizar
problemas sobre el flujo de sustancias empleando ordenadores que realizan miles de clculos
para simular la evolucin de fluidos en movimiento en ciertos problemas de inters para
la ingeniera y la ciencia en general. Existen incontables problemas de este tipo que pueden ser
resueltos mediante el clculo numrico pero la necesidad de los esquemas numricos en este trabajo
es indispensable por no contar con soluciones analticas para el sistema de ecuaciones que modela el
problema.
El problema que trataremos numericamente es el lid-driven cavity para el caso bidimensional,
que consiste en un fluido newtoniano incompresible contenido y limitado por un plano cerrado
cuadrado libre de la influencia de fuerzas externas, cuando la tapa o recta superior del cuadrado se
desplaza horizontalmente a velocidad constante conocida, mientras las dos paredes o rectas verticales
y la base o recta inferior del cuadrado permanecen en reposo, analizaremos la evolucin del flujo
confinado en la cavidad transcurrido un tiempo finito definido. Gracias a las condiciones de contorno
y geometra simples del problema, la descripcin matemtica del planteamiento de ste se simplifica
de tal modo que podemos acceder a soluciones numricas de gran ineters a nivel cientfico dado que
hasta el momento no sabemos si exista o no una solucin exacta fsicamente razonable al problema
planteado. El planteamiento trata de un sistema de ecuaciones diferenciales parciales no lineales que
describen la evolucin en el tiempo de un fluido en movimiento, dicho sistema es conocido como las
ecuaciones de Navier-Stokes.
vi
CONTENIDO
Como se puede ver en [1], la no linealidad de las ecuaciones de Navier-Stokes crea la necesidad de un tratamiento numrico especial a los trminos convectivos, la meta siempre es implementar
un esquema numrico que asegure tanto presicin como estabilidad, lo que se logra con los esquemas
upwind de alta resolucin; un esquema upwind es una extensin del mtodo de diferencias finitas
que simula numericamente la direccin de propagacin de la variable transportada en el campo del
flujo, ganando estabilidad y presicin en los resultados a medida que aumenta el orden o resolucin.
Como se puede ver en [2], el fundamento matemtico de los esquemas upwind de alta orden se basan
en los criterios de limitacin TVD-CBC garantizando de manera definitiva la estabilidad y presicin
en los resultados eliminando por completo soluciones oscilatorias sin sentido fsico alguno.
El lid-driven cavity para el caso bidimensional formulado en variables primitivas u,v,p en la
literatura se ha convertido en un problema clsico cuyas investigaciones se han centrado en el anlisis
de la variacin en de los resultados, que dependen del mtodo numrico utilizado, la discretizacin
de los trminos convectivos, el nmero de Reynolds y el tamao de la malla computacional utilizada
para discretizar el problema.
Se puede ver en [1] que utilizando el mtodo de diferencias finitas, mallas de 2048 x 2048,
nmeros de Reynolds de hasta 10000 y diferentes discretizaciones incluyendo upwind de alta orden
a los terminos convectivos, lograron explicar ciertos aspectos sobre los fenomnos que generan la
inestabilidad de las soluciones numricas del problema cuando se trabaja con altos nmeros de
Reynolds.
Se puede ver en [3] tambin el desarrollo de un cdigo computacional rpido y compacto basado
en la optimizacin de la estructura del lenguaje de clculo tcnico matlab utilizando el mtodo de
diferencias finitas para dominios rectangulares del lid-driven cavity para el caso bidimensional.
Este trabajo est basado en el mtodo de diferencias finitas con tratamiento especial a los trminos convectivos cuyos esquemas de discretizacin upwind se fundamentan en criterios de
limitacin TVD-CBC garantizando estabilidad en las soluciones numricas y as comparar algunos
resultados con las literaturas especializadas de [1] y [3].
Objetivos
GENERAL
Estudiar numericamente el lid-driven cavity para el caso bidimensional implementando esquemas
upwind de alta orden
ESPECFICOS
i Aplicar esquemas upwind de alta orden basados en criterios de limitacin TVD-CBC a los trminos
convectivos de las ecuaciones de Navier-Stokes
ii Analizar y comparar los resultados numricos obtenidos en cuanto a nivel de estabilidad y presicin
con los resultados obtenidos en la literatura especializada.
vii
Captulo
Preliminares
Como se puede ver en [4], las ecuaciones de Navier-Stokes modelan el flujo de los fluidos compresibles e incompresibles, turbulentos y laminares. Dependiendo de las propiedades del fluido, se
pueden simplificar considerablemente las ecuaciones, buscando reducir las dificultades inherentes a
su solucin numrica.
Las ecuaciones de Navier-Stokes representan una expresin matemtica que se deduce a partir
de principios fsicos muy familiares:
a. Conservacin de la masa
b. Conservacin del momento(Segunda ley de Newton)
c. Conservacin de la energa (Primera ley de la termodinmica)
1
Captulo 1. Preliminares
(1.1)
2 u 2 u
+
x2
y2
2
v (vv) (uv)
1 p
v 2 v
+
+
=
+
+
t
y
x
y
x2 y2
u (uu) (uv)
1 p
+
+
=
+
t
x
y
x
(1.2)
(1.3)
En el que se consideran y la densidad y viscosidad cinemtica del fluido como constantes con
unidades de medida del sistema internacional.
Una adimensionalizacin de las ecuaciones de Navier-Stokes para flujos incompresibles puede
ser efectuada a partir de las siguientes cantidades adimensionales
p =
p
0 V02
u =
u
V0
v =
v
V0
x =
x
L0
y =
y
L0
t =
tV0
L0
(1.4)
u
v
+
=0
x
y
(1.5)
p
1
u (u u ) (u v )
+
+
=
t
x
y
x
Re
2 u
2 u
+
(x )2 (y )2
2
v (v v ) (u v )
v
p
1
2 v
+
+
= +
+
t
y
x
y
Re (x )2 (y )2
(1.6)
(1.7)
0 V0 L0
(1.8)
2 (u u )
2 (u v ) 2 (v v )
+
=
+
+
+2
+
(x )2 (y )2
Re (x )2 (y )2
t
(x )2
y x
(y )2
Donde
=
v
u
+
=0
x
y
u
x
(1.9)
(1.10)
v
y
u v
+
=0
x y
(1.11)
Captulo 1. Preliminares
2 u 2 u
+
x2
y2
v (vv) (uv)
p
1 2 v 2 v
+
+
=
+
+
t
y
x
y Re x2 y2
2 p 2 p
1 2 2
2 (uu)
2 (uv) 2 (vv)
+
=
+
+2
+
+
x2 y2
Re x2
y2
t
x2
yx
y2
p
1
u (uu) (uv)
+
+
=
+
t
x
y
x Re
1.1.2
(1.12)
(1.13)
(1.14)
Ahora presentamos el desarrollo clave a la obtencin de resultados numricos consistentes del problema, sin entrar en detalle a cerca del proceso de disicretizaciones tanto espaciales como temporales
de las variables, ya que el mtodo de correccin de la presin nos indica el camino a la solucin
numrica independientemente de la discretizacin utilizada.
Empecemos Reescribiendo las ecuaciones (1.12) y (1.13) despejando
2 u 2 u
+
x2
y2
2 v 2 v
+
x2 y2
u
1
=
t
Re
v
1
=
t
Re
u
t
v
t
(uu) (uv)
p
+
x
y
x
(vv) (uv)
p
y
x
y
respectivamente.
(1.15)
(1.16)
Identifiquemos y simbolicemos los trminos viscosos y convectivos de las ecuaciones (1.15) y (1.16)
2 u 2 u
+
x2
y2
(1.17)
(uu) (uv)
+
x
y
(1.18)
2 v 2 v
+
x2 y2
(1.19)
(vv) (uv)
+
y
x
(1.20)
L(u) =
C(u) =
L(v) =
C(v) =
(1.21)
v
1
p
=
L(v) C(v)
t
Re
y
(1.22)
Ahora definamos los operadores y que sean quienes representen las discretizaciones espaciales
de las derivadas de primero y segundo orden respectivamente pero en la discretizacin temporal
utilizamos euler explcito slo por simplicidad en el desarrollo del mtodo de proyeccin ya que como
se mencion anteriormente, el desarrollo del mtodo no depende de la descretizacin utilizada, y
finalmente empleamos superndices n y n + 1 que se refieren a las soluciones en los instantes tn y
tn+1 .
La discretizacin de las ecuaciones (1.21) y (1.22) nos quedan:
un+1 un
1
=
[L(u)]n [C(u)]n pn+1
x
t
Re
(1.23)
1
vn+1 vn
=
[L(v)]n [C(v)]n pn+1
y
t
Re
(1.24)
1
= u + t
[L(u)]n [C(u)]n pn+1
x
Re
n
vn+1 = vn + t
1
[L(v)]n [C(v)]n pn+1
y
Re
(1.25)
(1.26)
(1.27)
Del sistema de ecuaciones discretizadas (1.25)-(1.26) y (1.27) se desconocen los trminos pn+1
y
x
pn+1
y , lo que nos lleva a eliminar dichos trminos al igual que la ecuacin (1.27), y as obtenemos las
siguientes ecuaciones
1
n
n
[L(u)] [C(u)]
F = u + t
Re
1
n
n
n
G = v + t
[L(v)] [C(v)]
Re
(1.28)
(1.29)
Captulo 1. Preliminares
(1.30)
vn+1 G = t(pn+1
y )
(1.31)
Reescribamos y discretizemos ahora la ecuacin de Poisson (1.14) que nos indicar Cunto debe
valer la presin para corregir los valores de F y G y as obtener un+1 y vn+1 garantizando soluciones
numricas consistentes con la ecuacin de continuidad (1.11).
2 p 2 p
1 2 2
(uu) (uv)
(vv) (uv)
+
=
+
+
+
+
+
x2 y2
Re x2
y2
t
x
x
y
y
y
x
(1.32)
Sustituimos ahora los trminos convectivos por C(u) y C(v) y los viscosos por L(u) y L(v) tenemos
2 p 2 p
1 L(u) L(v)
C(u) C(v)
+
=
+
+
+
x2 y2
Re
x
y
t
x
y
(1.33)
Aplicando linealidad de la derivada parcial y factorizando trminos del lado derecho de la ecuacin
(1.33) tenemos:
2 p 2 p
1
+ 2 =
L(u) C(u) +
L(v) C(v)
2
x
y
x Re
y Re
t
El trmino
(1.34)
discretizacin se har de manera detallada siguiendo la convencin adoptada del operador y euler
explcito para la variable temporal
=
t
t
u v
+
x y
(1.35)
(1.36)
s aplicamos la ecuacin 1.27 ntese que el trmino un+1 + vn+1 en la ecuacin 1.36 se anula y la
ecuain nos queda:
un + vn
t
(1.37)
Discreticemos ahora toda la ecuacin (1.34) con la convencin adoptada y sustituyendo el resultado
discretizado de 1.37 nos queda
pn+1
x
pyn+1
1
=
L(u) C(u)
Re
n
1
+
L(v) C(v)
Re
n
+
un + vn
t
(1.38)
La ecuacin (1.38) es una discretizacin de la ecuacin de Poisson que cumple con la ecuacin (1.27)
garantizando la incompresibilidad del flujo en tn+1 , Podemos expresar la ecuacin (1.38) en trminos
de las expresiones F y G , si sustituimos las ecuaciones (1.28) y (1.29) en esta ltima.
pxn+1
pyn+1
=
F un
t
+
G vn
t
+
un + vn
t
(1.39)
Aplicando la linealidad del operador y simplificando trminos del lado derecho de (1.39) finalmente
obtenemos:
1
pxn+1 + pn+1
=
(F + G )
y
t
(1.40)
A partir de la ecuacin (1.40) obtenemos los valores de las presiones paso a paso en el tiempo tn+1
garantizando la incompresiblidad del flujo corrigiendo los valores de F y G , para finalmente clacular
las velocidades u y v en el tiempo tn+1 con las ecuaciones (1.30) y (1.31) que se transforman en:
un+1 = F t(pn+1
x )
(1.41)
vn+1 = G t(pn+1
y )
(1.42)
Captulo 1. Preliminares
Como se puede ver en [4], es importante resaltar el hecho de que el no tomar una malla
escalonada (staggered grid) puede producir soluciones inconsistentes para los campos de presin
debido a que si todas las variables u, v, p se toman en un mismo nodo o punto de la malla se pueden
originar campos oscilatorios que destruyen las soluciones numricas, las cuales deben ser fisicamente
consistentes.
En lo que respecta a la discretizacin temporal debemos tener presente que cada partcula fluida no puede recorrer una distancia mayor que el tamao de la clula computacional en cada paso de
tiempo t, es lo que se denomina condicin CFL.
Para establecer las condiciones de contorno sobre la malla debemos tener cuidado ya que estas
no se pueden imponer directamente, esto por el hecho de que siempre falta alguna variable sobre
el contorno o frontera por lo que se incluyen nuevas incgnitas en el sistema, a las que llamaremos
clulas externas a la malla, teniendo en cuenta la condicin de que el valor interpolado (promediado)
en la frontera sea la condicin de contorno correspondiente como lo ilustra la figura 1.3.
10
Captulo 1. Preliminares
Como se puede ver en [4], para el desarrollo del cdigo utilizamos el mtodo de discretizacin
temporal Euler explcito, aplicando el principio de la sesin 1.1.2. de paso fraccionado para forzar la
incompresibilidad; las ecuaciones que resultan son:
1
F = u + t
[L(u)]n [C(u)]n
Re
1
[L(v)]n [C(v)]n
G = v + t
Re
(1.28)
(1.29)
1
pxn+1 + pn+1
=
(F + G )
y
t
(1.40)
un+1 = F t(pn+1
x )
(1.41)
vn+1 = G t(pn+1
y )
(1.42)
Ahora para la discretizacin espacial igual que [7] teniendo en cuenta la malla escalonada de la
figura (1.3) con nodos xi , i = 1, ..., Nx , yj , j = 1, ..., Ny y uniformemente espaciada con x = xi+1 xi
y y = yj+1 yj . El nmero de puntos Nx y Ny se corresponde con la cantidad de puntos interiores
del dominio, y se utilizan celdas externas para imponer las condiciones de contorno como lo ilustra la
grfica 1.3.
En la distribucin de las variables y sus derivadas en una malla escalonada hay que tener en
cuenta que cada vez que calculemos la primera derivada de una variable, esta cambiar su posicin
en la malla. Sin embargo con la segunda derivada mantendr su posicin. Por ello es importante que
cada vez que sumemos dos trminos, estos tengan la misma posicin en la malla. En caso contrario,
ser necesario interpolar uno de ellos para colocarlo en la posicin adecuada. Para ello utlizamos
interpolacin lineal, por ejemplo, para mover ui,j a Ui,j+ 1 haremos
2
Ui,j+ 1 =
2
Ui,j+1 + ui,j
2
(1.43)
Tanto al principio como al final de cada paso temporal, las variables u, v y p ocuparn las posiciones
que se muestran en la figura 1.4. Esto hace que las variables u, v y p tengan dimensiones distintas.
La tabla 1.2 muestra las dimensiones de cada variable segn consideremos sus puntos interiores o
interiores + exteriores
11
Puntos interiores
(Nx 1) x Ny
(Nx + 1) x (Ny + 2)
Nx x (Ny 1)
(Nx + 2) x (Ny + 1)
Nx x Ny
(Nx + 2) x (Ny + 2)
uni+ 1 ,j
2
2n
u2n
i+1,j ui,j
(1.44)
1
+t
Re
x2
"
Gi,j+ 1
2
vni,j+ 1
2
2n
v2n
i,j+1 vi,j
y2
(un vn )i+1,j+ 1 (un vn )i,j+ 1
2
+
(1.45)
+t
1
Re
n+1
n+1
pn+1
i+1,j 2pi,j + pi1,j
x2
vni+1,j+ 1
2
2vni,j+ 1
2
x2
vni1,j+ 1
2
n+1
n+1
pi,j+1
2pi,j
+ pn+1
i,j1
y2
n+1
ui+
1
,j
2
n+1
vi,j+
1
2
Fi+ 1 ,j
2
Gi,j+ 1
2
1
=
t
vni,j+ 1 +1
2
2vni,j+ 1
2
y2
Fi+ 1 ,j Fi 1 ,j
2
n+1
pn+1
i+1,j pi,j
Gi,j+ 1 Gi,j 1
2
!
, (1.46)
n+1
pn+1
i,j+1 pi,j
!
vni,j+ 1 1
2
(1.47)
(1.48)
12
Captulo 1. Preliminares
Para el desarrollo del primer cdigo computacional, como se puede ver en [4], utilizamos el
esquema FOU para la discretizacin de los trminos convectivos de las ecuaciones (1.44) y (1.45), lo
que significa que adoptaremos un ajuste sobre las velocidades que nos permita acercarnos al caracter
continuo de estas variables ofreciendo mayor estabilidad a las soluciones numricas.
ui 1 + ui+ 1
2
(1.49)
El superndice e se refiere a una estimacin de la velocidad u, as por ejemplo, cuando uei 0 el flujo
va de izquierda a derecha. Podemos entonces definir la velocidad ui del flujo en funcin de uei
ui =
u 1
i
si
uei 0
ui+ 1
si
uei < 0
2
2
(1.50)
La expresin (1.50) se puede expresar en trminos de la funcin signo que permite otorgar el signo
de las velocidades de manera ms prctica
Si = sgn(uei ) =
si
uei 0
si
uei < 0
(1.51)
13
1 + Si
2
ui 1 +
2
1 Si
2
ui+ 1
(1.52)
Para la discretizacin de los trminos viscosos de las ecuaciones (1.44) y (1.45) utilizamos diferencias
finitas adelantadas.
Ahora, como se puede ver en [5], para el clculo de las presiones en la ecuacin (1.46), se
nos presenta un sistema lineal muy particular que resolvemos por mtodos iterativos en la que
el laplaciano de la presin genera una matriz simetrica definida positiva en el sistema, al cual le
aplicamos el mtodo SOR (Mtodo de sobrerelajacin sucesiva) que es una ligera modificacin del
mtodo de Gauss-seidel y permite mejorar la convergencia, debido a las condiciones del sistema lineal.
Para la aplicacin del mtodo SOR tal como vemos en [8], verificamos previamente el valor
optimo para la constante de relajacin el cual est dado analticamente para regiones rectangulares
con condiciones de frontera tipo Dirichlet, condicin que se cumple en el planteamiento del problema,
y se define as:
opt = 2
donde
cos
+ 2 cos
1 + 2
(1.53)
2
(1.54)
x
y
14
Captulo 1. Preliminares
1/2
XX
kRk2 =
R2i,j ,
i
Donde
Ri,j =
n+1
n+1
pi+1,j
2pi,j
+ pn+1
i1,j
x2
(1.55)
n+1
n+1
pn+1
i,j+1 2pi,j + pi,j1
y2
(1.56)
Fi+ 1 ,j
2
Fi 1 ,j
2
Gi,j+ 1
2
Gi,j 1
2
Garantizada la convergencia del mtodo SOR, procedemos entonces a utilizar los resultados del
campo de presiones en las ecuaciones (1.47) y (1.48) para lograr obtener el campo de velocidades
y volver a realizar el ciclo hasta llegar al tiempo final establecido previamente y alcanzar el estado
estable en la solucin.
El diagrama de flujo que ilustra el algoritmo de solucin de las ecuaciones de Navier-Stokes
para el caso bidimensional est en la figura 1.5
15
16
Captulo 1. Preliminares
1.3.1
Con la construccin de los esquemas upwind de alta orden como se hace en [9] se busca mantener
estabilidad en las soluciones as como tambin resolver gradientes elevados de manera simultnea.
Consideremos las posiciones D(Downstream), R(Remote-Upstream) y U(Upstream) en relacin
con el frente f de la clula computacional mostrada en la figura 1.6. Definamos entonces una variable
genrica normalizadad de la siguiente manera.
^=
(R)
(D) (R)
(1.57)
donde (D) y (R) son, respectivamente los valores no normalizados del valor de en los puntos D y
R, por ejemplo, ^(f) = (f) (R)
(D) (R)
17
1.7, los esquemas CDS (Central-differencing Scheme), FOU (First-Order Upwind), QUICK (Quadratic
Upstream Interpolation for convective Kinematics), SOU (Second-Order Upwind) estn representados
en la figura 1.7(a).
Figura 1.7 (a) Esquemas convectivos clsicos de variables no normalizadas, (b) Criterio CBC [10]
Cualquier esquema no lineal, formulado en variables normalizadas (NV) que pase por los puntos
O(0,0), P(1,1), Q(0.5,0.75) de la figura 1.7(a) es de segunda orden de precisin y si tambin pasa
por el punto Q con inclinacin de 0.75 es de tercera orden. Se recomienda que para valores de ^U
menores que 0 o mayores que 1 debe ser utilizado el esquema FOU.
En vista de la importancia de soluciones limitadas al transporte de propiedades fsicas igual que en [9],
propusieron un criterio de limitacin CBC (Convection Boundedness Criterion). En el contexto NV,
un esquema convectivo produce solucin convergente si este est enteramente contenido en la regin
CBC mostrada en la figura 1.7(b), esto es, el esquema debe satisfacer las siguientes condiciones
^f [^U , 1],
^ = 0,
f
^
f = 1,
^ = ^ ,
f
si ^f [0, 1]
si ^U = 0
si ^U = 1
si ^f 6 [0, 1]
(1.58)
A pesar de que el criterio CBC igual que en [9] trata los problemas de estabilidad adecuadamente,
18
Captulo 1. Preliminares
^f [^U , 2U ] y ^f 1,
^f = ^U ,
si ^f [0, 1]
si ^U 6 [0, 1]
(1.59)
Basado en los anteriores criterios de limitacin TVD-CBC es que se crean los diferentes esquemas
upwind de alto orden.
Es importante resaltar como en [9], que esquemas numricos de baja y alta orden como FOU,
SOU, QUICK, Lax-Wendroff no exhiben propiedades de acotamiento a corde a uno de los teoremas
ms importantes del anlisis numrico y la CFD denominado El teorema de Godunov, el cual se
enuncia como en [2]:
Todo esquema numrico lineal monotono para la ecuacin de conveccin es necesariamente de primer orden de precisin.
Los esquemas TOPUS (Third-Order Polynomial Upwind Scheme) definido igual que en [9], EPUS
(Eight-degree Polynomial Upwind Scheme) definido como en [10] y SUPERBEE definido en [11],
cumplen con todas las propiedades que permiten mantener precisin y estabilidad simultneamente
en las soluciones numericas.
1.3.2
El esquema TOPUS
En este trabajo utilizamos el esquema TOPUS, desarrollado por [12], para la discretizacin de los
trminos convectivos, su idea bsica es utilizar una combinacin de esquemas de baja y alta orden a
travs de una funcin de intercambio (limitador). En particular se considera un polinomio de cuarto
grado completamente contenido en la regin CBC y se imponen las condiciones como en [9], esto es,
el esquema debe pasar por los puntos O(0,0), Q(0.5,0.75) (para llegar a segunda orden), P(1,1) y
con pendiente 0.75 en Q (para alcanzar tercera orden) ver figura 1.7(a).
El esquema convectivo TOPUS en funcin del parmetro libre en NV est definido por:
^f =
^4U + (2 + 1)^3U +
^U ,
510
4
^2
U +
+10
4
^U ,
^U [0, 1]
^U 6 [0, 1]
(1.60)
19
Aplicando la definicin de NV, el esquema TOPUS resultante en variables no normalizadas est dado
por
f =
R + (D R ) ^4U + (2 + 1)^3U +
510
4
^2
U +
+10
4
^U ,
U ,
^U [0, 1]
^U 6 [0, 1]
(1.61)
donde:
^U =
1.3.3
U (R)
(D) (R)
El esquema EPUS
Analogamente al esquema TOPUS, el esquema EPUS considera ahora un polinomio de octavo grado
completamente contenido en la regin CBC, imponiendo condiciones igual que en [10], siendo estas
idnticas a las del esquema TOPUS.
El esquema convectivo EPUS en funcin del parmetro libre en NV est entonces definido
por:
^ 8
^ 7
^ 6
^ ,
U
^U [0, 1]
^U 6 [0, 1]
(1.62)
Aplicando la definicin de NV, el esquema EPUS resultante en variables no normalizadas est dado
por:
^ 8
^ 7
^ 6
^ ,
^ 6 [0, 1]
U
donde:
^U =
U (R)
(D) (R)
(1.63)
20
1.3.4
Captulo 1. Preliminares
El esquema SUPERBEE
El esquema SUPERBEE definido en [11] cumple con las restricciones TVD, y la funcin tambin pasa
por el punto Q(0.5,0.75), pero a diferencia de los esquemas TOPUS y EPUS no se considera una
versin de esquema no lineal. En variables no normalizadas se define como sigue:
U ,
2U (R) ,
f =
1
2
1
2
^U 6 [0, 1]
0 < ^U 1
(D) + U ,
3U (R) ,
(D) ,
1
3
1
2
2
3
< ^U
< ^U
3
1
2
2
3
< ^U 1
(1.64)
21
El modelo matemtico adimesionalizado que gobierna la evolucin del flujo dentro de la cavidad
despus de algn tiempo transcurrido, est definido por el sistema:
u v
+
=0
x y
u (uu) (uv)
p
1
+
+
=
+
t
x
y
x Re
(1.65)
2 u 2 u
+
x2
y2
(1.66)
22
Captulo 1. Preliminares
p
1
v (vv) (uv)
+
+
=
+
t
y
x
y Re
2 v 2 v
+
x2 y2
(1.67)
u(x, 1, t) = u0 ; v(x, 1, t) = 0
u(0, y, t) = v(0, y, t) = 0
u(1, y, t) = v(1, y, t) = 0
u(x, 0, t) = v(x, 0, t) = 0
(1.68)
Las ecuaciones en derivadas parciales (1.65 a 1.67), representan las leyes fsicas de conservacin de
la masa y conservacin de los momentos en las direcciones x e y respectivamente, la ecuacin (1.65)
tambien suele llamarse condicin de incompresibilidad del fluido ya que la densidad es considerada
constante en todo el dominio espacial y temporal.
Ntese tambien que las condiciones de contorno (1.68) enmarcan el problema de tal modo que la
cavidad corresponde a un cuadrado unitario cuyo vertice inferior izquierdo coincide con el origen de
un sistema de coordenadas cartesianas en el plano bidimensional y para la garanta de la existencia de
soluciones numricas la variable temporal se fija en un intervalo cerrado finito con tiempo final t = T
Captulo
Resultados Numricos
tf
dt =
100
23
24
Las velocidades u en cada clula computacional en la direccin x con la cavidad girada 90 grados
sentido horario son:
2.5949e 07
2.8353e 06
3.3053e 06
6.4923e 06
6.9643e 06
5.3955e 06
2.5438e 06
1.0533e 06
5.6893e 08
0
0.02211
0.01258
0.02714 0.046367 0.070124 0.10404 0.13753 0.18409
0.0581
0.34707 1
Las velocidades v en cada clula computacional en la direccin y con la cavidad girada 90 grados
sentido horario son:
0 1.5689e 06
0.0020328
0.0059045
0.0069242
0.0030217
0.00017285
0.0049843
0.0063223
0.0063452
0.0028074
1.1005e 07
6.0394e 07
2.2634e 06
6.4734e 06
9.2002e 06
7.3004e 06
3.6148e 06
7.8241e 07
0.017167
0.036785
0.069292
0.10529
0.14852
0.1904
0.22435
3.8987e 07 0
0.20883
0.024846
0.056769
0.098467
0.15093
0.20259
0.2431
0.24881
0.15857
0
0.023034
0.049799
0.091449
0.13674
0.18506
0.2093
0.19049
0.10569
0
0.012741
0.025716
0.048183
0.084701
0.11589
0.13866
0.12069
0.064623
0
0.0033997
0.0052138
0.0084191
0.00090321
0.025678
0.042785
0.052984
0.03023
0
0.015872
0.036275
0.058744
0.086807
0.093763
0.059457
0.03072
0.0012219
0
0.024224
0.051425
0.093376
0.14106
0.19627
0.20488
0.12076
0.04904
0
0.022091
0.050672
0.09196
0.15306
0.22325
0.30601
0.29726
0.104
0
0.012579
0.029947
0.058946
0.10007
0.16298
0.23761
0.34707
0.29571
0
1.1066e 06 3.5215e 06 6.3216e 06 8.2976e 06 8.8457e 06 7.3841e 06 4.6114e 06 8.8925e 07 0
Las presiones p en cada clula computacional con la cavidad girada 90 grados sentido horario
son:
0.10244
0.16454
0.26684
0.45172
0.64711
0.14162
0.21502
0.33126
0.50955
0.77296
1.0234
0.10838
0.14955
0.10137
0.1358
0.1038
0.06164
0.02121
0.03323
0.04003
0.36052
0.55846
0.21713
0.3127
0.43225 0.54755 0.61886 0.31986
0.18888
0.25968 0.33217
0.3897
0.41115 0.20398
0.14153
0.19057 0.24182
0.279
0.28526 0.13956
0.0227
0.090402
0.13459
0.13215
0.12516
0.062608
0.032386
0.13574
0.31347
0.47356
0.49427
0.23167
0.012508
0.12196
0.36185
0.8626
1.3624
0.7642
0.0034784 0.053039
0.16933
0.51364
0.98934
0.59789
25
26
Bibliografa
[1] Charles-Henri Bruneau, Mazen Saad, The 2D lid-driven cavity problem revisited, Universit
Bordeaux 1, 2005, 331-334.
[2] C. Hirsch, Numerical computation of internal and external flows, second edition, vol
1, 2007, 337-370
[3] B. Seibold, A compact and fast Matlab code solving the incompressible Navier-Stokes
equations on rectangular domains, Massachusetts Institute of Technology, 2008
[4] A. de Oliveira Fortuna, Tcnicas Computacionais para Dinmica dos Fluidos. Conceitos Bsicos e Aplicaes, edusp.
[5] D. Sperandio, J. Texeira, L. Monken, Caractersticas Matemticas e computacionais
dos mtodos Numricos, Pearson Education, so Paulo, 2003, 335-337
[6] C.H. Marchi, The Lid-Driven Square Cavity Flow: Numerical Solution with a 1024 x
1024 Grid, Federal University of Paran UFPR, 81531-980 Curitiba, PR, Brazil, July-September
2009
[7] A. Lozano, Notas sobre dinmica de fluidos computacional, Rev.0.5.2, 24 de marzo
de 2014 versin online
27
28
BIBLIOGRAFA
[8] J.D. Hoffman, Numerical Methods for Engineers and Scientists, second edition Revised and Expanded, MacGraw Hill, New York 1992. 548-549
[9] R.A. Bofim de Queiroz, Desenvolvimento e teste de esquemas "upwind" de alta resoluo e suas aplicaes em escoamentos incompressveis com superfcies livres, Tesis de
Maestra, USP - So Carlos, marzo de 2009
[10] L. Correa, Um novo esquema upwind de alta resoluo para equaes de conservao no
estacionrias dominadas por conveco, Tesis de Maestra, USP- So Carlos, Abril/2011
[11] P.L. Roe. Some contribuitions to the modeling of discontinuous flows. Lect. Notes
Appl. Math, 22:163-193, 1985
[12] R.A.B. de Queiroz and V.G. Ferreira. Development and testing of hight-resolution
upwind schemes - Upwind schemes for incompressible free surface flows. VDM Verlag Dr M
uller,
2010