ANÁLISIS DE FLUJO Y TERMODINÁMICO POR MEDIO DE ANSYS A UN

MODELO DE REFRIGERACIÓN TERMOELÉCTRICO COMPACTO

PROYECTO FINAL

CARLOS ARTURO BOHORQUEZ

CARLOS ARTURO PEREZ BERNAL

SERGIO ESTEBAB CULMA ORTÍZ

UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS

FACULTAD TECNOLÓGICA
INGENIERÍA MECÁNICA
2021
BOGOTÁ D.C.
DELIMITACIONES DE ESTUDIOS

A la hora de poder observar el comportamiento de un fluido de índole gaseosa es importante el

tener en cuenta variables que puedan afectar su comportamiento a través de un elemento, así
mismo, propiedades tales como presión velocidad y funciones de posición (tanto en análisis 2D
como 3D), pueden ofrecer parámetros importantes a la hora de realizar estudios, tal es el ejemplo
de un estudio termodinámico que se abordará más adelante con mayor detalle, tal es el caso de
nuestro radiador de aluminio a una temperatura de 5° C, el cual será objeto de estudio debido a
que es un elemento y pieza esencial en un proyecto de refrigeración llevado a cabo en una
compañía, como una parte de un sistema de refrigeración termoeléctrico.

Ilustración No.1 (Conjunto de elementos de refrigeración termoeléctrico)

Las características del elemento ensamblado se ilustran en la imagen previamente expuesta, como
se puede observar existe una caja de almacenamiento del aire para deshumedecer en una
proporción importante el aire, previamente al ingreso del fluido gaseoso al radiador de baja
temperatura:

Ilustración No.2 (Conjunto de elemento de refrigeración termoeléctrico)

Como podemos observar (ilustración No.2), el objeto posee unas características particulares a la
hora de usarse como conductor térmico de las bajas temperaturas del efecto termoeléctrico
desarrollado en el proyecto con el fin de reducir la energía térmica del aire atmosférico que pase a
través de los múltiples canales del elemento mecánico de aluminio, en éste caso se realizará con
propiedades estándar del aluminio para los cuales se verá reflejado el uso de los siguientes datos:

Tabla No.1 (algunas propiedades térmicas del Aluminio)

Es con esto que la parametrización de análisis de fluidos, térmico y de presión del elemento y del
fluido (aire) que recorre los canales previamente mencionados, toma un papel importante a la
hora de determinar la utilidad y la viabilidad del desarrollo del proyecto.

ANÁLISIS DE FLUIDOS:

Para dar inicio a los análisis y resultados obtenidos en el análisis de fluidos se dispondrán primero
los ajustes iniciales generales para cada estudio en específico, esto con el fin de otorgar un
panorama de estandarización acerca del estudio de las interacciones del aire con la pieza de
aluminio:

Se tuvieron en cuenta las siguientes propiedades del aire:

Tabla No.2 (algunas propiedades del aire)

1. Análisis de movimiento y velocidad: Se realizaron pruebas de velocidad y así mismo de
direccionamiento de las partículas a través del dispositivo, lo cual nos indicó un valor
m
máximo de 2,0156 , al ser parametrizado como un gas ideal y configurarse de manera
s
que la velocidad a la que el fluido llega al dispositivo a través de la región ¨inlet¨ tiene una
m
magnitud vectorial paralela al eje y Y perpendicular al plano XZ de 1 . Se efectúa una
s
alteración en la velocidad en sentido de aumento a más del doble debido al cambio de
sección que se efectúa al entrar en el radiador a baja temperatura, este análisis se vuelve
aún mas evidente al aplicar la ecuación de continuidad de Bernoulli:

Ecuación No.1 (algunas propiedades del aire)

Despreciando los efectos de la gravedad y teniendo a la velocidad únicamente como

función de la presión y la densidad en el caso de un gas ideal, tenemos que:

2
V 2=(k −P) , donde ρ = P/ RT
ρ
Ecuación No.2 (algunas propiedades del aire)

Una de las ventajas de cálculo que ofrece el programa es la de condición de temperatura

estática, puesto que tanto sabemos existe una variedad en los cálculos cuando los cuerpos
restauran su contenido energético debido a su propia temperatura, por lo que la única
variable en la Ecuación No.2 es la presión que se manifiesta de la siguiente manera:

Ilustración No.3 (Velocidades al ingreso y en medio del recorrido del radiador a baja temperatura)

Ahora bien, para que un fluido pueda obtener un aumento de velocidad a través de una
sección deben haber dos condiciones esenciales para que se de ésta condición, primero
debe haber una condición de caída de la presión confirme su velocidad aumente y
segundo debe haber un estrangulamiento o reducción de sección que permita transformar
parte de esa alta presión en velocidad para el fluido, tal como lo evidencia ANSYS, vamos a
evaluar cada caso en particular y respectivamente en el mismo orden en el que fueron
enunciadas con el fin de respaldar lo que en el programa se desarrolla.
 Ilustración No.4 (Estrangulamiento del flujo de aire)

Como podemos apreciar en la ilustración No.4, las líneas cuyo corregimiento de color
m
tiende hacia el azul, denota la velocidad de entrada bajo las condiciones de 1 , a medida
s
que el fluido gaseoso (aire) va ingresando en los canales del radiador a baja temperatura
su velocidad comienza a sufrir un cambio en cuanto al incremento de su magnitud neta en
dirección a la cara parametrizada como ¨outlet¨, siendo esta misma la salida del aire bajo
las condiciones de salida del radiador, mas no las condiciones iniciales.

2. Análisis de presión: Como ya habíamos mencionado anteriormente los modelos de

velocidad muestran un incremento en la velocidad del fluido al poder aplicar y entender la
ecuación de continuidad de Bernoulli expresada la velocidad como función de la presión
del sistema (Ecuación No.2). es por ello que se hace necesaria una observación detallada
sobre un mapa de contorno con gama de presiones distribuidas sobre el área superficial
del elemento:
 Ilustración No.5 (Distribución de presión)
Se logra apreciar en la ilustración No. 5 la distribución de presión ejercida por el fluido
sobre las paredes del elemento a nivel interno de los canales, dicha distribución de presión
se distribuye a lo largo del canal, formando un visible decrecimiento del nivel de presión
ejercida, pues tal como se abordó en la ecuación No.2, así como el cuadrado de la
velocidad puede ser función del cuadrado de la presión, la misma ecuación puede
escribirse de la siguiente forma:

2 KRT 2kRT
V 2 +2 RT = ≤¿ 2 =P
P V +2 RT
Ecuación No.3 (algunas propiedades del aire)

La expresión nos explica (Ecuación No.3) que efectivamente existe una relación
inversamente proporcional a la cantidad de presión ejercida a través del recorrido por los
múltiples canales del dispositivo respecto de su velocidad, puesto que como habíamos
dicho antes (véase ilustración 4), los efectos de estrangulación y aumento de velocidad se
traducen en caídas de presión; para este caso hablamos que la presión máxima esta dada
en 21.56 Pa, sin tener encuenta el área roja, ya que por temas de geometría es razonable
que la presión ejercida directamente al contacto del muro con respecto a la entrada
¨inlet¨ sea la de mayor categoría al estar tan cerca a la superficie solida.

3. Análisis térmico: Por parte del análisis térmico tenemos que la variación dada por las
condiciones de temperatura estática se aproximan a una interacción donde las
propiedades se mantienen, sin embargo al ser dinámica la forma en que se presenta la
interacción por convección térmica entre el aire y el radiador a baja temperatura,
podemos inferir que la transferencia es exitosa y que la temperatura del aire se ve
reducida en un intervalo entre 4° y 3° C, por lo que el proyecto resulta bastante
interesante al evitar los temas de condensación de agua al ser regulado y calculado para
evitar dicha humidificación directa en el elemento mecánico sin necesidad de usar la caja
deshumidificadora en la que se hizo alusión en la ilustración número 1:

Ilustración No.6 (Distribución de temperatura estado inicial)

 Ilustración No.7 (Distribución de temperatura estado final)

Se puede concluir que en cuanto a iteraciones funciona bastante bien con pocos time
steps para las iteraciones de categorías como presión y velocidad, ya que se lograron
adquirir unos resultados de gran volumen y continuidad para los cuales existen valores
que representan una tendencia de datos hacia características de precisión y exactitud con
los cuales se desarrolla mayor confianza a la hora de evaluar los datos obtenidos con las
simulaciones y cálculos de ANSYS.

Mientras que con los térmicos es necesario lograr una muy buena cantidad de iteraciones
en el momento de alcanzar un estándar bueno en cuanto a visualización del fenómeno en
cuestión, y no solo time steps, si no también el intervalo entre cada uno de ellos, es decir,
debe haber un gran volumen de time steps con una cantidad moderada de iteraciones
para unos pocos segundos de visualización de interacción e intercambio energético, es por
ello que se debe tener presente los temas de enmallado para poder refinar los cálculos sin
necesidad de saturar el sistema siempre y cuando se obtengan la mejor tendencia de
datos por las características previamente mencionadas de precisión y exactitud, puesto
que estas pautas de datos nos pueden aproximar y dar mejor idea acerca de los sucesos
que ocurren entorno a un fenómeno como las interacciones térmicas por medio de
convección.