UNIVERSIDAD NACIONAL DE INGENIERIA
Facultad de Ingeniería Mecánica
ASME International
Sección estudiantil
Universidad Nacional de Ingeniería
“Aplicaciones de Mecánica de Fluidos
Computacional (CFD) en la Industria.”
Ing. Luis Tay Wo Chong Hilares
Temario
• Introducción sobre CFD
• ¿Cómo es un proceso de análisis con
CFD?
• Aplicaciones de CFD en la industria
• Referencias
• Apreciaciones finales
Mecánica de fluidos computacional
Definición
z La mecánica de fluidos computacional (CFD) es el análisis de
sistemas que envuelven flujos de fluidos, transferencia de
calor y otros fenómenos como reacciones químicas por medio
de simulaciones numéricas.
z El problema físico a ser resuelto puede ser descrito
matemáticamente por las ecuaciones de conservación de
masa, momentum y energía en su forma diferencial parcial.
z El CFD reemplaza las derivadas parciales en estas
ecuaciones con formas algebraicas discretizadas
Mecánica de fluidos computacional
Ventajas
z Reducción de costos y tiempo en el desarrollo de nuevos
diseños.
z Permite probar diversas variaciones hasta alcanzar el
resultado optimo, todo esto antes de realizar el prototipo.
z Reducción en prototipos permite ahorros en la creación
de plantas piloto, pruebas de testeo y rediseño de
productos.
z Capacidad de diseñar a grandes escalas y a condiciones
difíciles donde es difícil hacer mediciones (altas
presiones, o altas temperaturas por ejemplo).
z Permite comprender lo que sucede dentro del sistema a
analizar. Por ejemplo, permite mostrar al comprador
como trabajará su equipo al mostrar lo que esta pasando
dentro de él.
z Llegada del producto más rápido al mercado
Mecánica de fluidos computacional
Ecuaciones de transporte
z Conservación de Masa
∂ρ r
+ ∇ ⋅ (ρ u ) = 0
∂t
z Conservación de momentum
r r r
∂ρ u rr ∂u ∂u j 2
+ ∇ ⋅ (ρ uu ) = −∇ p + ∇ ⋅ τ ij + ρ ∑ y k f k + g τ ij = i +
∂x
− µ δ ij ∇ ⋅ ur
3
∂t j ∂xi
z Conservación de energía
[ ]
∂ρ ht r ∂p r r r
∂t
+ ∇ ⋅ ( ρ ht u ) =
∂t
+ ∇ ⋅ J k + u τij + ρ f k + Rh ht = h +
2
(
1 2 2
ux + u y + uz
2
)
z Ecuación de estado
p = ρ M k RT
Metodo numérico de volumenes finitos
z Reducir el sistema de ecuaciones diferenciales
parciales a un sistema algebraico mediante un
proceso de discretización e integración sobre un
volumen.
La ecuación general de transporte
∂ (ρφ ) r
+ ∇ ⋅ (ρφu ) = ∇ ⋅ (Γ ∇φ ) + Sφ
∂t
Término Término Término Término
transitivo Convectivo Difusivo fuente N
Esquema de discretización W
n
P E
w e
a PφP = ∑ a nbφnb + b s
S
v
u
Modelos de turbulencia
z Ecuaciones de Navier Stokes promediadas de
Reynolds (RANS)
r r r
τ = ρ u'u'
( )
t
∂ρ U rr ij
+ ∇ ⋅ ρ UU = −∇ p + ∇ ⋅ τ ij + ∇ ⋅τ t ij
∂t Transferencia de
momentum por la
turbulencia
z Aproximación de Boussinesq
Relacionar los esfuerzos de Reynolds con los
gradientes de velocidad y los modela usando una
viscosidad turbulenta análoga a la viscosidad
molecular.
τ t
ij
r r
( 2 r
) ∂U i ∂U j
= ρ u ' u ' = ρ ui ' u j ' = ρk + µt ∇ ⋅ U δ i j − µt
+
∂xi
3 ∂x j
k: Energía cinética turbulenta
Modelos de turbulencia
z Ecuaciones del modelo de turbulencia k-ε
- Energía cinética turbulenta (k)
k=
1
2
1
(
ui ' ui ' = u x ' 2 + u y ' 2 + u z ' 2
2
)
(ρk ) + ∇ ⋅ ( ρkU ) = ∇ ⋅ (µ∇k ) + ∇ ⋅ µt
∂ r
(∇k ) + Gk − ρε + S k
∂t σk
- Disipación (ε)
∂ui ' ∂u i '
ε =ν
∂x j ∂x j
(ρε ) + ∇ ⋅ ( ρε U ) = ∇ ⋅ (µ∇ε ) + ∇ ⋅ µ t ε2
∂ r
(∇ε ) + ρC1 S ε − C 2ε ρ + Sε
∂t σε k + νε
Metodo númerico de volumenes finitos
z Cálculo de la presión
Inicio
∂p
a Pui , P = ∑ anb ui , nb − + S P Asumir p*
∂xi P
Resolver las ecuaciones
S P + ∑ anbui , nb 1 ∂p discretizadas de momentum
ui , P = −
aP a P ∂xi P ui*
1 ∂p
ui , P = u~ *i , P − Resolver la ecuación de
a P ∂xi P corrección de la presión
p’
∂ (ρu~ *i , P )
ρ ∂p p*=p, ui*=ui
∂
=
Ec. Continuidad
∂xi aP ∂xi P ∂xi Corregir las velocidades y la
presión p añadiendo p’ a p*
z El algorítmo SIMPLE p, ui
p = p * + p' ui = ui * +u' Resolver las otras ecuaciones
discretizadas de transporte (temperatura,
especies, turbulencia, etc.)
Presión Presión Corrección
correcta Asumida de la Presión
Convergencia
u 'i , P =−
∑a nb u 'i ,nb 1 ∂p '
−
aP a P ∂xi P
Parar
1 ∂p '
u'i , P = −
a P ∂xi P
¿Cómo es un proceso de análisis con
CFD?
Paso 1
• Identificar el problema físico
• Identificar los modelos físicos a ser usados:
Turbulencia, combustión, multifásico, etc.
• ¿Qué precisión queremos?
¿Cómo es un proceso de análisis con
CFD?
Paso 2
• Crear la geometría
• Generar la malla
¿Cómo es un proceso de análisis
con CFD?
Paso 3
• Importar y revisar la malla en el programa de
CFD
¿Cómo es un proceso de análisis
con CFD?
Paso 4
• Seleccionar los modelos físicos apropiados
(turbulencia, combustión, transferencia de calor,
multifásico, etc).
¿Cómo es un proceso de análisis
con CFD?
Paso 5
• Definir los materiales
¿Cómo es un proceso de análisis con
CFD?
Paso 6
• Definir las condiciones de frontera
¿Cómo es un proceso de análisis
con CFD?
Paso 7
• Definir los esquemas de solución e iniciar el
caso
¿Cómo es un proceso de análisis con
CFD?
Paso 8
• Realizar el post-procesamiento
¿Cómo es un proceso de análisis
con CFD?
Paso 8
• Validación de los resultados
Sloan D, Smith P, Smoot D. “Modeling of swirl in turbulent flow systems”. Prog. Energy Combust.
Sci. Vol 12, pp 163-250, 1986 (No pertenece al caso presentado)
APLICACIONES
Áreas de trabajo en CFD
- Automotriz
- Aeroespacial
- Procesos químicos
- Medicina
- Aplicaciones ambientales
- Turbomáquinas
- Flujos bifásicos
- Transferencia de calor
- Reología, etc.
Sector automotriz
• Estudio de flujos
sobre el automóvil y
dentro del automóvil.
• Análisis de
catalizadores y flujos
de escape
Imagenes de Tecplot Inc., Fluent
Turbomáquinas
• Estudio de cargas térmicas en alabes de turbinas
• Interacción estator/rotor
Imagenes de Tecplot Inc., Fluent, CD-Adapco, cfdesign
Sector aeroespacial
Imagenes de Tecplot Inc., Fluent, CD-Adapco, Fieldview
Coeficiente de presión
Accesorios y tuberias
Distribución
de presión
Accesorios y tuberias
• Prediccion de erosion en una tuberia
Soldadura
Electrode -0.04
-0.08
-0.12
-0.16
-0.2
-0.24
Constant voltage -0.28
Wire feedrate -0.32
-0.36
-0.4
-0.44
10 ms 30 ms 35 ms 40 ms 45 ms 50 ms 100 ms 150 ms -0.48
r
-0.52
-0.56
-0.6
Jsij
200 ms 300 ms 345 ms 355 ms 359 ms 362 ms 372 ms 385 ms
Calculation domain 393 ms 399 ms 450 ms 500 ms 600 ms 650 ms 680 ms 688 ms
-0.4 -0.2 0 0.2 0.4
Wall
Symmetrical axis
Fig. 1 -- Schematic sketch of metal transfer
process in GMAW with boundary conditions
Imagenes de http://www.engr.uky.edu/~cfd/welding.html y fluent
Intercambiadores de calor y aire
acondicionado
Sin a/c Con a/c
Imagenes de Tecplot Inc., Fluent, CD-Adapco, Fieldview
Aplicaciones en Ingenieria civil y
arquitectura
• Para edificios
- Cargas generadas por
vientos en edificios y
puentes.
- Optimizar el confort y
seguridad en espacios
urbanos (generación de
vientos entre edificios,
vientos en calles, etc.).
• Para comprender como
actuaría en una habitación
un sistema de calefacción,
ventilación y aire
acondicionado.
• Para optimizar el control
térmico en un edificio
residencial
Referencias
• www.fluent.com
• www.ansys.com/products/cfx.asp
• www.cd-adapco.com
• www.cfd-online.com
• An introduction to Computational Fluid Dynamics –
Versteeg; Malalasekera.
• Numerical Heat Transfer and Fluid Flow – Patankar.
• Computational Methods for Fluid Dynamics –
Ferziger; Peric.
• Turbulence Modeling for CFD – Wilcox.
Apreciaciones finales
• El uso de CFD en la industria es una necesidad
debido a los diversos beneficios que ofrece en
el desarrollo de productos.
• El desarrollo de CFD en las universidades es
una necesidad para continuar las
investigaciones en los modelos con CFD.
• El uso de CFD no es la solución final, ya que los
resultados deben ser validados. Sin embargo, la
ayuda que proporciona es de gran utilidad como
herramienta de diseño.
Apreciaciones finales
• La aplicación del CFD debe trabajar de la
mano con las mediciones experimentales,
las cuales tendrán la función de “calibrar”
los modelos de CFD para el desarrollo de
nuevos diseños.
• La industria requiere de productos más
baratos, de más rápida producción y
mejores.
ASME International
Sección estudiantil
Universidad Nacional de Ingeniería
Ing. Luis Tay Wo Chong
[email protected] GRACIAS
