1 - Introducción Al CFD Con ANSYS
1 - Introducción Al CFD Con ANSYS
1 - Introducción Al CFD Con ANSYS
análisis fluido
dinámico con ANSYS
MULTIPHYSICS
1. ¿Que es CFD?
2. ¿Dónde se Usa?
3. Conceptos Básicos del CFD
4. Introducción al FLOTRAN
5. Aplicaciones del FLOTRAN
6. Pasos involucraron en el CFD & FLOTRAN
7. Potencialidades del Solucionador (Solver) del
FLOTRAN
8. Interfaz para el Usuario del FLOTRAN
¿Que es CFD?
Combustión de carbón en un
horno de quemadores múltiples
¿Dónde se Usa?
•Arquitectura
•Medioambiente
Superficie de
Control
Volumen
de
Control
Superficie de
Control
Volumen
de
Si N = m entonces η= 1 Control
xx xy xz
. g V x V x V x V y V x V z V x
x y z x
t x y z
yy yx yz V y V y V y V y
. g V x V y V z
y x z z
t x y
y
zz zx zy
. g V z V x V z V y V z V z V z
z x y z
t x y z
Ecuaciones Constitutivas
Muchos fluidos exhiben una relación lineal entre las componentes del esfuerzo y
las velocidades de deformación. Estos fluidos denominados NEWTONIANOS por
lo general son asociados con las características isotrópicas, entonces es posible
relacionar las componentes del esfuerzo y las velocidades de deformación
empleando solo dos propiedades del flujo, la viscosidad μ y el segundo
coeficiente de viscosidad λ, que en conjunto forman las seis ecuaciones
constitutivas que relacionan los esfuerzos normales y cortantes, con la presión, y
los gradientes de velocidades, planteadas por Navier (1785-1836) y Stokes (1819-
1903), estas son las siguientes:
z
2
x
yy p
x
y
z y
V x V y V x
2 V x
zz p
x
y
z z
2
Según la hipótesis de Stokes se cumple que:
3
Ecuaciones de Navier-Stokes
x
2 V V V z
y t x y
y 2 2 x y z
y z
Vz V z V z V z V z
2 2 2
p V V z
g z
2
V V V
x z
z y z t x y
z 2 2 x y z
Estas ecuaciones son conocidas como las de Navier-Stokes y representan el
flujo de un fluido viscoso en su forma general. Su solución analítica no ha
sido lograda, siendo considerada uno de los grandes retos mundiales para la
comunidad científica.
D V
p
g V
Dt
Ecuaciones de la Energía Flujo Compresible
La conservación de energía puede expresarse en términos de la temperatura de
estancamiento total, (a menudo es útil en los flujos muy comprimibles, o
temperaturas estáticas, y apropiado para velocidades bajas en los análisis de
fluidos incompresibles. Dicha ecuación es:
t
.C p .T 0
x
.V x C p .T 0
y
.V y C p .T 0
z
.V z C p .T 0
T 0 T 0 T 0 P
K. K. W E Q
V K
K.
x x y y z z V t
Donde:
Cp = Calor específico
To = temperatura total (o estancamiento)
K = conductividad térmica
Wv = termino de trabajo viscoso
Qv = la fuente de calor volumétrica
Φ = termino de generación de calor viscosa
Ek = energía cinética (definida posteriormente)
Ecuaciones de la Energía Flujo Compresible+
La temperatura estática es calculada de la temperatura total de la energía cinética:
T T0 V
Donde: 2.C p
T = temperatura estática
v = magnitud del vector de velocidad del fluido
Cp = calor especifico del fluido
u
u j .. . uk
V j
W .
xi xi xk x j
K 1 2 K 1 2 K 1 2
. .V . .V . .V
k
E x C p x 2 y C p y 2 z C p z 2
2 . 2
x
y z y x z x z y
2 2
Ecuaciones de la Energía Flujo Incompresible
T 0 T 0 T 0
K. K. K. Q
x x y y z z V
t
C p .T o V C p T o K T QV
El Fenómeno de La Turbulencia
Si los efectos inerciales son bastante grandes con respecto a los efectos
viscosos, entonces el flujo se vuelve turbulento. La Turbulencia se mide en
función la velocidad instantánea, esta fluctúa en cada punto del campo de flujo.
La velocidad se expresa en términos de un valor medio y un componente
fluctuante Tal como se aprecia en la figura siguiente:
V’
V V
V x V x V x
'
De la figura se deduce que:
Donde: V x
Componente x de la velocidad instantánea
V x
Componente x de la velocidad media
'
V x
Componente x de la fluctuación
El Fenómeno de La Turbulencia +
Después de la substitución de la ecuación anterior en las ecuaciones
constitutivas, las velocidades en las ecuaciones de momento son valores
promediados, y se representan colocando la barra en la expresión subsiguiente
de las ecuaciones de momento. Mientras que las ecuaciones constitutivas toman
la forma:
' '
'
.V x .V z
' ' '
.V x .V x .V x .V y
R
x
x y z
' ' '
z x .V z .V x y .V z .V y z .V z .V z
R ' ' '
Donde:
σR = Termino de esfuerzo de Reynolds
El Fenómeno de La Turbulencia +
Una aproximación de la viscosidad de remolino al modelo turbulento se agrega al
término de esfuerzo viscoso con un coeficiente desconocido, la viscosidad
turbulenta queda así:
V x
.V x .V y .
t y
La principal ventaja de este cambio, es que la representación del σR es de la
misma forma que las condiciones de difusión en las ecuaciones originales. Las
dos condiciones pueden combinarse en una viscosidad eficaz que se define
como la suma de la viscosidad del flujo laminar y la viscosidad turbulenta, y se
expresa así:
e
t
t
. Ls
2
Donde:
μt = Viscosidad Turbulenta
Φ = viscosidad de disipación L si Lx 0
x
Lx = Escala de longitud
Ln = la distancia más corta del nodo a la pared más cercana
Lc = Escala de longitud característica (el valor más grande de Ln encontrado)
Modelo k-ε
Para el Modelo k-ε la viscosidad turbulenta es calculada como una función de los
parámetros de la turbulencia la energía cinética k y su rapidez de dispersión ε que
usa la Ecuación anterior y queda expresada así:
2
C . . k
Donde:
t
k =energía cinética turbulenta
ε = rapidez de disipación de la energía cinética turbulenta
t k t k t k
.k
t
.V x .k
x
y V y
z
. .k .V z .k . . .
x k x y k y z k z
C . . T T T
. .
4
.t
g . g . g .
t
t x x y y z z
t t t
.
t
.V x .
x
y V y
. . .V z .
z
. . .
x x y y z z
2
C .1 C 3. . .k T T T
C1, t k C 2 k
. . . .
. g .
x x g . g .
z z
t y y
Modelo k-ε Normal+
Se pueden expresar de forma compacta así:
.k C . .
.V .k t
.k . .
4 t
. g .T
k t
t t
.
2
C .1 C 3. . .k
. C 1, . . C 2 . .
.V . t
. g .T
t
t t k k
C2 1,92
Cμ 0,09
Coeficientes del Modelo Normal σk 1,0
σε 1,3
σt 0,9
C3 1,0
C4 0,0
β 0,0
Modelo k-ε Normal+
La solución a las ecuaciones de turbulencia se usa para calcular la viscosidad
eficaz y la conductibilidad térmica eficaz:
2
.C
C .. K K K t p
e e
t
Donde:
μe = viscosidad efectiva
Ke = conductividad efectiva
σt = Numero Prandtl Turbulento (Schmidt)
Modelo Turbulento RNG
En el modelo de RNG, La constante C1,ε en la Ecuación anterior de dispersión, se
reemplaza por la siguiente función:
. 1
C 1,42
1 .
1, 3
En el modelo RNG Cμ, C2 son constantes. Como se mostró en la tabla anterior, los
multiplicadores de difusión tienen valores diferentes que el modelo predefinido (k-ε
Normal). El valor de la constante rotacional Cr en el modelo de RNG es 0. Los
valores contemplados en la Ecuación 2.56 que no son especificados en la tabla
anterior son cubiertos por la tabla siguiente
Variable Valor
β 0.12
C2 1.68
Cμ 0.085
σk 0.72
σε 0.72
4.38
Inicio Algoritmo de Solución
Fin
CONDICIONES DE CONTORNO
Nota:
Los esfuerzos de esquila de Pared están dados por la siguiente expresión:
V tan 1 E
ln
w
w
k
Donde:
μ = es la viscosidad absoluta
δ = distancia del nodo de la pared al nodo cercano a la pared
vtan = velocidad de un nodo paralelo a otro cercano a la pared
E = constante del modelo de turbulencia (por defecto es 9.0)
κ = constante del modelo de turbulencia (por defecto es 0.4)
τw = esfuerzos de esquila de Pared
También es usado: y
w
Donde: y+ es una distancia adimensional de la pared