Análisis Computacional de Flujo en Canal Abierto con Módulos Generadores de Calor Profesor: Álvaro Valencia Alumnos: David Azocar Paredes Rodrigo Reyes Tapia DEPARTAMENTO DE INGENIERIA MECANICA UNIVERSIDAD DE CHILE Universidad de Chile Facultad de Ciencias Físicas y Matemáticas Departamento de Ingeniería Mecánica
Literatura
COMPUTATIONAL ANALYSIS OF TURBULENT FLUID FLOW AND HEAT TRANSFER OVER AN ARRAY OF HEATED MODULES USING TURBULENCE MODELS
 
Resultados
Formulación del Problema
Geometría del Problema
Geometría del Problema Parametros del canal: Altura del canal: H=10 cm Largo del canal: Lt=80 cm Largo de entrada: Le=8 cm Separación de módulos: S1=7 S2=12 S3=12 S4= 7 cm Alto del módulo: L/B=1, 0.75, 3 (B=3 fijo) Calor generado: qj= 80
Ecuaciones Ecuación de Continuidad: Ecuación de N-S según X: Ecuación de N-S según Y: Ecuación  de Energía:
Variables Adimensionales
Definición del Caso
Casos Estudiados Caso 0: v=0.1 y L/B=1.5 Caso 1: v=1 y L/B=1.5 Caso 2: v=3 y L/B=1.5 Caso 3: v=10 y L/B=1.5 Caso 4: v=10 y L/B=0.75 Caso 5: v=10 y L/B=3
Método Numérico
Volúmenes Finitos El dominio de cálculo es dividido en volúmenes de control disjuntos. Cada uno rodea un nodo de la red.  Para obtener la ecuación discretizada se integra la ecuación diferencial sobre cada VC. De esta manera la ED cumple el principio de conservación en un volumen finito, tal como la ec. diferencial lo cumple en un volumen infinitesimal.
Volúmenes Finitos P: punto genérico E, W: puntos vecinos este y oeste. Dx= ancho del VC. e, w, caras del volumen de control , que rodea al punto genérico
Discretización de Ecuación La ecuación:  se integra, obteniendo una ecuación de la forma: Luego se discretiza, para obtener:
Interpolación Diferencias Centrales Aguas Arriba Esquema Exponencial Híbrido Ley de Potencia
Ley de Potencia Este esquema se basa en interpolar “  T   ”, mediante una   aproximación del esquema exponencial.  Dicha aproximación se logra reemplazando la curva   solución del esquema exponencial, por dos rectas y dos polinomios de quinto orden
Algoritmo SIMPLEC 1.-Considerando el campo inicial de velocidades, se calculan los coeficientes de la ecuación de movimiento y luego se calculan las seudo-velocidades a partir de las ecuaciones de movimiento sin gradiente de presión. 2.-Con las seudo-velocidades se resuelve el campo de presiones usando la ecuación de continuidad. Con el campo de presiones obtenido se resuelven las ecuaciones de movimiento para determinar un nuevo campo de velocidad. 3.-Con el nuevo campo de velocidades, se resuelve nuevamente la ecuación de continuidad apara determinar el campo de correcciones de presión. 4.-A partir de las correcciones de presión, se corrige el campo de velocidades, obteniendo parámetros finales para ocupar en la ecuación de energía. 5.-Con estos datos se resuelve la ecuación de energía y se obtiene el campo de temperatura. 6.-Luego se determina si el número de Nusselt converge o no. Si el resultado de la iteración no converge a una cota requerida, el ciclo se repite desde el cálculo de los coeficientes de la ecuación de movimiento. 7.-Luego de lograr la convergencia, se incrementa el tiempo y se aplican de nuevo los pasos anteriores para obtener los campos de velocidad y temperatura en un nuevo instante de tiempo.
Resultados
Velocidades
Caso 0
Caso 1
Caso 2
Caso 3
Caso 4
Caso 5
Isotermas
Caso 0
Caso 1
Caso 2
Caso 3
Caso 4
Caso 5
Nusselt
Caso 0
Caso 1
Caso 2
Caso 3
Caso 4
Caso 5
Comparación con Literatura
Líneas de flujo obtenidas en la bibliografía
Coeficiente de caída de presión
Coeficientes de caída de presión obtenidos en la bibliografía
Conclusiones
El caso laminar presenta la mayor variación del Nusselt del primer modulo con respecto a los siguientes. Los casos turbulentos presentan Nusselt casi constante en los módulos. Al aumentar la altura de los módulos se generan mas turbulencias Con esto, la transferencia de calor aumenta.

Más contenido relacionado

PPT
Tratatimiento numerico de ecuaciones diferenciales (2)
DOCX
Descarga de condensadores
PPTX
Conducción bidimensional
PPTX
Conduccion bidireccional
PPTX
Trasnferencia de calor
DOCX
Conducción bidimensional de calor en estado estacionario
PPTX
Diaps 2do parcial (1)
PPTX
Asignación #4, samantha arrambide
Tratatimiento numerico de ecuaciones diferenciales (2)
Descarga de condensadores
Conducción bidimensional
Conduccion bidireccional
Trasnferencia de calor
Conducción bidimensional de calor en estado estacionario
Diaps 2do parcial (1)
Asignación #4, samantha arrambide

La actualidad más candente (19)

PDF
Ilovepdf merged (35)
PPTX
Transferencia de calor bidireccional
PPTX
Ley de gauss
PPTX
Diaps 2do parcial
PPTX
Mecanica de fluidos Bernoulli
PPTX
Flujo horizontal (1)
PPT
Ley de gauss
PDF
Presentacion de circuitos trifacicos.
PPTX
Presentación1
PPTX
Fluidos
PPTX
Electromagnetismo: Ley de Gauss
PDF
Tarea teoria e..
PPSX
Ley de gauss clase 5 ok
PPSX
Ley de Gaus
PPTX
Fisica 2 explicacion
PDF
Ley de gauss
PPTX
Eg035 conduccion superficies extendidas 2
PPTX
Eg051 lmtd
PPT
Ley De Gauss
Ilovepdf merged (35)
Transferencia de calor bidireccional
Ley de gauss
Diaps 2do parcial
Mecanica de fluidos Bernoulli
Flujo horizontal (1)
Ley de gauss
Presentacion de circuitos trifacicos.
Presentación1
Fluidos
Electromagnetismo: Ley de Gauss
Tarea teoria e..
Ley de gauss clase 5 ok
Ley de Gaus
Fisica 2 explicacion
Ley de gauss
Eg035 conduccion superficies extendidas 2
Eg051 lmtd
Ley De Gauss
Publicidad

Destacado (6)

DOCX
Transición de un canal
PDF
hidrologia-estudiantes-ing-civil
PDF
Problemas selectos de fenomenos de transporte
PPTX
PARALISIS CEREBRAL INFANTIL - PCI
PDF
HIDRAULICA DE CANALES - PEDRO RODRIGUEZ
PDF
Solucionario De Fenomenos De Transporte R Byron Bird
Transición de un canal
hidrologia-estudiantes-ing-civil
Problemas selectos de fenomenos de transporte
PARALISIS CEREBRAL INFANTIL - PCI
HIDRAULICA DE CANALES - PEDRO RODRIGUEZ
Solucionario De Fenomenos De Transporte R Byron Bird
Publicidad

Similar a Canal con generadores de vortices (20)

DOC
212art rogd1
PDF
Transferencia de calor jose manrique valadez - segunda edicion
PDF
Transferencia de calor, segunda edición [manrique] decrypted
PDF
1512 7159-1-pb
PDF
5106 22307-1-pb
PDF
Tema 4. primer principio sa
DOCX
Aplicación de análisis numérico en alabes
PPTX
PPT exposición CFD Dinámica de fluidos computacional.pptx
PDF
Intro to CFD II.pdf
PDF
Simulación numérica de una llama de difusión con un coflujo axial mediante AN...
PDF
Ejercicios
PDF
10 intercambiadores de-calor
PDF
Mecánica de fluidos II, Francisco Ugarte
PDF
358847270 mecanica-de-fluidos-2-ugarte-pdf
PDF
Ecuación de la energía
PPTX
Trabajo de Termodinámica.pptx
DOCX
Parctica #6 termo
PPTX
Transferenciadecalorbidireccional
212art rogd1
Transferencia de calor jose manrique valadez - segunda edicion
Transferencia de calor, segunda edición [manrique] decrypted
1512 7159-1-pb
5106 22307-1-pb
Tema 4. primer principio sa
Aplicación de análisis numérico en alabes
PPT exposición CFD Dinámica de fluidos computacional.pptx
Intro to CFD II.pdf
Simulación numérica de una llama de difusión con un coflujo axial mediante AN...
Ejercicios
10 intercambiadores de-calor
Mecánica de fluidos II, Francisco Ugarte
358847270 mecanica-de-fluidos-2-ugarte-pdf
Ecuación de la energía
Trabajo de Termodinámica.pptx
Parctica #6 termo
Transferenciadecalorbidireccional

Más de Alvaro Valencia (7)

PPT
CFD en modelo de Aneurisma Basilar
PPT
Flujo turbulento
PPT
PPT
Laminar Flow
PPT
Flujo 2 Cilindros
PPT
Burbujas ideales
PPT
Aneurismas cerebrales
CFD en modelo de Aneurisma Basilar
Flujo turbulento
Laminar Flow
Flujo 2 Cilindros
Burbujas ideales
Aneurismas cerebrales

Canal con generadores de vortices

  • 1. Análisis Computacional de Flujo en Canal Abierto con Módulos Generadores de Calor Profesor: Álvaro Valencia Alumnos: David Azocar Paredes Rodrigo Reyes Tapia DEPARTAMENTO DE INGENIERIA MECANICA UNIVERSIDAD DE CHILE Universidad de Chile Facultad de Ciencias Físicas y Matemáticas Departamento de Ingeniería Mecánica
  • 3. COMPUTATIONAL ANALYSIS OF TURBULENT FLUID FLOW AND HEAT TRANSFER OVER AN ARRAY OF HEATED MODULES USING TURBULENCE MODELS
  • 4.  
  • 8. Geometría del Problema Parametros del canal: Altura del canal: H=10 cm Largo del canal: Lt=80 cm Largo de entrada: Le=8 cm Separación de módulos: S1=7 S2=12 S3=12 S4= 7 cm Alto del módulo: L/B=1, 0.75, 3 (B=3 fijo) Calor generado: qj= 80
  • 9. Ecuaciones Ecuación de Continuidad: Ecuación de N-S según X: Ecuación de N-S según Y: Ecuación de Energía:
  • 12. Casos Estudiados Caso 0: v=0.1 y L/B=1.5 Caso 1: v=1 y L/B=1.5 Caso 2: v=3 y L/B=1.5 Caso 3: v=10 y L/B=1.5 Caso 4: v=10 y L/B=0.75 Caso 5: v=10 y L/B=3
  • 14. Volúmenes Finitos El dominio de cálculo es dividido en volúmenes de control disjuntos. Cada uno rodea un nodo de la red.  Para obtener la ecuación discretizada se integra la ecuación diferencial sobre cada VC. De esta manera la ED cumple el principio de conservación en un volumen finito, tal como la ec. diferencial lo cumple en un volumen infinitesimal.
  • 15. Volúmenes Finitos P: punto genérico E, W: puntos vecinos este y oeste. Dx= ancho del VC. e, w, caras del volumen de control , que rodea al punto genérico
  • 16. Discretización de Ecuación La ecuación: se integra, obteniendo una ecuación de la forma: Luego se discretiza, para obtener:
  • 17. Interpolación Diferencias Centrales Aguas Arriba Esquema Exponencial Híbrido Ley de Potencia
  • 18. Ley de Potencia Este esquema se basa en interpolar “ T ”, mediante una aproximación del esquema exponencial. Dicha aproximación se logra reemplazando la curva solución del esquema exponencial, por dos rectas y dos polinomios de quinto orden
  • 19. Algoritmo SIMPLEC 1.-Considerando el campo inicial de velocidades, se calculan los coeficientes de la ecuación de movimiento y luego se calculan las seudo-velocidades a partir de las ecuaciones de movimiento sin gradiente de presión. 2.-Con las seudo-velocidades se resuelve el campo de presiones usando la ecuación de continuidad. Con el campo de presiones obtenido se resuelven las ecuaciones de movimiento para determinar un nuevo campo de velocidad. 3.-Con el nuevo campo de velocidades, se resuelve nuevamente la ecuación de continuidad apara determinar el campo de correcciones de presión. 4.-A partir de las correcciones de presión, se corrige el campo de velocidades, obteniendo parámetros finales para ocupar en la ecuación de energía. 5.-Con estos datos se resuelve la ecuación de energía y se obtiene el campo de temperatura. 6.-Luego se determina si el número de Nusselt converge o no. Si el resultado de la iteración no converge a una cota requerida, el ciclo se repite desde el cálculo de los coeficientes de la ecuación de movimiento. 7.-Luego de lograr la convergencia, se incrementa el tiempo y se aplican de nuevo los pasos anteriores para obtener los campos de velocidad y temperatura en un nuevo instante de tiempo.
  • 43. Líneas de flujo obtenidas en la bibliografía
  • 44. Coeficiente de caída de presión
  • 45. Coeficientes de caída de presión obtenidos en la bibliografía
  • 47. El caso laminar presenta la mayor variación del Nusselt del primer modulo con respecto a los siguientes. Los casos turbulentos presentan Nusselt casi constante en los módulos. Al aumentar la altura de los módulos se generan mas turbulencias Con esto, la transferencia de calor aumenta.