Diferencia entre revisiones de «Torres de enfriamiento hiperbólicas (Grupo 06)»
(→Fuerza total y por unidad de superficie del viento sobre la mitad expuesta) |
(→Fuerza total y por unidad de superficie del viento sobre la mitad expuesta) |
||
| Línea 339: | Línea 339: | ||
<center><math>A_{\text{half}}=\pi \sqrt{1 + R'^2}\left[R_{\min} H+\frac{\Delta R}{2}H\right].</math></center> | <center><math>A_{\text{half}}=\pi \sqrt{1 + R'^2}\left[R_{\min} H+\frac{\Delta R}{2}H\right].</math></center> | ||
<br> | <br> | ||
| − | <center><math>\boxed{A_{\text{half}}=\pi H \sqrt{1 + \left( \frac{\Delta R}{H} \right)^2}\left( R_{\min} + \frac{\Delta R}{2} \right)}. | + | <center><math>\boxed{A_{\text{half}}=\pi H \sqrt{1 + \left( \frac{\Delta R}{H} \right)^2}\left( R_{\min} + \frac{\Delta R}{2} \right)}.section*{Evaluación numérica}A_{\text{half}} \approx 2.0303909708635903 \times 10^{4}\ \text{m}^2,\qquad(\approx 2.03\times 10^{4}\ \text{m}^2).</math></center> |
== Aplicación practica del hiperboloide de revolución == | == Aplicación practica del hiperboloide de revolución == | ||
Revisión del 12:59 28 nov 2025
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Torres de enfriamiento hiperbólicas (Grupo 06) |
| Asignatura | Teoría de Campos |
| Curso | 2025-26 |
| Autores |
|
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Introducción
Las torres de enfriamiento hiperbólicas, son estructuras emblemáticas del paisaje energético. Estas se pueden encontrar en instalaciones como la Central Nuclear de Cofrentes (España) o en Niederaussem (Alemania). Son un elemento clave para la disipación térmica en el ciclo de producción de energía. Su diseño de hiperboloide de revolución de una hoja acelera el flujo de aire ascendente y enfría el agua de manera natural sin necesidad de usar ventiladores. La doble curvatura proporcionada por esta estructura hiperbólica otorga a la torre una enorme resistencia frente al viento y su propio peso.
En el desarrollo de este articulo vamos a analizar una torre de enfriamiento estándar de altura máxima ([math]H[/math]) con radio máximo ([math]Rmáx[/math]) y radio mínimo ([math]Rmín[/math]) el cual se encuentra a [math]\dfrac{\scriptsize 2}{\scriptsize 3}H[/math].La superficie de la torre se describe matemáticamente mediante la siguiente ecuación en coordenadas cartesianas:
Sean [math]a, c, z_0\gt0[/math] unos valores a determinar
Sean en nuestro modelo [math]Rmáx=55m[/math] ; [math]Rmín=30m[/math] ; [math]H=150m[/math]
2 La geometría hiperbólica en las torres de enfriamiento
Como se menciona anteriormente, estas estructuras siguen una morfología hiperbólica. La coordenada [math]z[/math] esta comprendida desde [math]z=0[/math] y [math]z=150[/math] al ser la altura máxima ([math]H=150m[/math]). Al ser la figura un solido de revolución, generado al girar una hipérbola alrededor del eje Z , tiene simetría radial.
2.1 Obtención de los parámetros
Como se ha visto anteriormente el parámetro [math]z_0[/math] es la posición vertical donde ocurre el estrechamiento, el cual está localizado a una altura de [math]\dfrac{\scriptsize 2}{\scriptsize 3}H[/math].
El parámetro [math]a[/math] representa el radio de la estrechez del hiperboloide, el cual es su radio mas pequeño [math]Rmín=30m[/math].
El parámetro [math]c[/math] es el factor de estiramiento vertical. Esta variable controla la curvatura de las paredes. Determina que tan rápido se ensancha la torre a medida que se alejas de la garganta.
Se necesita otro punto conocida para la obtención de [math]c[/math]. Se conoce que la base se encuentra a una altura [math]z=0[/math] y que tiene un radio [math]R=55[/math].
Cualquier punto del borde de la base cumple que [math]x^2+y^2=55^2[/math]. Para facilitar los cálculos se toma el siguiente punto de la base [math](x,y)=(55,0)[/math]. Se sustituye con este punto en la ecuación inicial.
Simplificando la ecuación
Elevando al cuadrado:
Despejando [math]c[/math]:
Los parámetros obtenidos son [math]a=30,z_0=100,c=65,07914[/math]
2.2 Ecuación de la superficie hiperbólica en coordenadas cilíndricas
Para obtener la ecuación en coordenadas cilíndricas se parte de la expresión anterior en coordenadas cartesianas:
Sean [math]a, c, z_0\gt0[/math] los valores determinados en el apartado anterior.
Con el objetivo de pasar la ecuación a coordenadas cilíndricas se definen los parámetros de las coordenadas cartesianas [math](x,y,z)[/math] en coordenadas cilíndricas [math](ρ,θ,z).[/math]
[math]y=ρ\cdot\sin(θ)[/math]
Cambiando [math]x^2+y^2[/math] a coordenadas cilíndricas:
[math]x^2+y^2=ρ^2\cdot\cos^2(θ)+ρ^2\cdot\sin^2(θ)[/math]
Se saca factor común [math]ρ^2[/math]:
Aplicando la siguiente identidad trigonométrica:
Sustituyendo en la ecuación
Al ser los demás datos constantes, no hace falta cambiar nada más. La ecuación en coordenadas cilíndricas queda:
3 Campo de temperaturas y transferencia de calor
El campo de temperaturas del aire en el interior de la torre se modeliza mediante la siguiente expresión:
Siendo:
- [math]T_{\text{base}}[/math] la temperatura en el centro de la base, en este caso [math]T_{\text{base}}=65ºC[/math];
- [math]\Delta T_z[/math] la caída de temperatura de base a tope, en este caso [math]\Delta T_z=38ºC[/math];
- [math]\Delta T_{\rho}[/math] la variación radial de temperatura, en este caso [math]\Delta T_{\rho}=8ºC[/math];
- [math]n[/math] el exponente de convección, en este caso [math]n=1,8[/math].
Sustituyendo estos valores en el campo y simplificando se llega a la siguiente expresión:
3.1 Cálculo de la masa del aire contenido en la torre
La densidad del aire varía con la temperatura, por ello y dadas las dimensiones de la torre de enfriamiento, hay que tener en cuenta esta variación para el correcto cálculo de la masa. La expresión que relaciona la densidad con la temperatura es la siguiente:
Donde [math]\rho_{\text{aire}}=1,225kg/m^3[/math] es la densidad del aire estándar a temperatura [math]T_0=15ºC[/math].
De este modo, se calcula la masa del aire dentro de la torre mediante la siguiente expresión:
Si se desarrolla la expresión se llega a lo siguiente:
Esta integral no es de fácil resolución, luego utilizando una herramienta de internet para realizar el cálculo numérico se llega al siguiente resultado:
Por otro lado, si se considera que todo el aire está a una temperatura [math]T=15ºC[/math] su densidad es constante y de valor [math]\rho_{\text{aire}}=1,225kg/m^3[/math]. Por lo tanto, la masa se calcula multiplicando dicha densidad por el volumen de aire contenido en la torre. Por lo que el problema pasa a ser calcular dicho volumen.
El volumen de la torre se calcula integrando el jacobiano en coordenadas cilíndricas en la región requerida:
En la sección 2.2 se obtiene la expresión del hiperboloide de revolución en coordenadas cilíndricas:
Despejando se obtiene una relación entre [math]\rho[/math] y [math]z[/math]
Volviendo a la integral, integrando primero [math]\theta[/math] y después [math]\rho[/math] se obtiene lo siguiente:
Sustituyendo [math]\rho^2(z) = a^2 \left[ 1 + \frac{(z - z_0)^2}{c^2} \right][/math] con [math]z_0=100[/math]:
Resolviendo ambas integrales:
Para la realización de [math]\int_0^{150} (z - 100)^2 dz[/math] usamos el cambio de variable [math]u = z - 100, du=dz[/math] con límites [math]u(0)=-100, u(150)=50[/math]
Sustituyendo en la integral, con [math]a=30[/math] y [math]c^2 =\frac{10000}{2,361}[/math]:
Por tanto, suponiendo que la masa de aire estuviera toda ella a 15ºC, esta tendría la siguiente masa:
Comparando ambos valores: [math]M_{\text{real}} \approx 781420 \text{ kg}[/math] y [math]M_{\text{15ºC}} \approx 825377 \text{ kg}[/math] se observa que el aire tiene una masa aproximadamente un 5% menor que la masa que tendría si estuviera a temperatura constante [math]T=15ºC[/math].
3.2 Cálculo de la temperatura media y la potencia térmica evacuada
La temperatura media se define como:
Se definirá [math]I_T=\iiint_V T(\rho, z) dV = \int_0^{2\pi} \int_0^{H} \int_0^{\rho(z)} T(\rho, z) \cdot \rho \, d\rho dz d\theta[/math]
Sustituyendo y definiendo los límites de integración se llega a la siguiente integral:
Debido al término exponencial esta integral no es de fácil resolución, luego utilizando una herramienta de internet se llega al siguiente resultado:
Luego la temperatura media será:
La potencia térmica evacuada puede calcularse mediante la siguiente expresión:
Siendo:
- [math]\rho_{\text{aire}} = 1,225 \text{ kg/m}^3[/math] la densidad del aire a temperatura [math]T_{\text{ambiente}}=15ºC[/math];
- [math]c_p = 1005 \text{ J/(kg K)}[/math] el calor específico a presión constante del aire y [math]T_{\text{ambiente}}=15ºC[/math];
- [math]Q = 1500 \text{ m}^3/\text{s}[/math] el caudal de aire;
- [math]T_{\text{media}} = 28,97 \text{ °C}[/math] la temperatura media calculada anteriormente;
- [math]T_{\text{ambiente}} = 15 \text{ °C}[/math] la temperatura ambiente.
Sustituyendo:
4 Análisis de la presión del viento
El viento es un factor que ejerce presión lateral sobre la torre de enfriamiento que varían a lo largo de la superficie en función de la altura. La velocidad escalar se puede representar de la siguiente forma:
[math]V_0[/math] es la velocidad de referencia a la altura [math]Zref[/math] que en nuestro modelo es [math]18m/s[/math].
[math]Zref[/math] es la altura de referencia sobre el suelo, que es [math]10m[/math].
[math]α[/math] es un exponente que depende del terreno, en nuestro modelo tendrá un valor de [math]0,1429[/math].
4.1 .-Representación del campo escalar de presiones.
Utilizando la velocidad del viento previamente mencionada, la presión del viento sobre la superficie de la torre se puede representar como:
Sea ρ la densidad del aire estándar, que tiene un valor de [math]1,225Kg/m^3[/math].
Sustituyendo en la fórmula obtenemos que la presión en función de la altura es:
Haremos el supuesto de que la torre es golpeada por un viento paralelo al vector [math]-\frac{1}{\sqrt{2}}(\vec{i} + \vec{j})[/math] para la mitad expuesta de la torre, se adjunta a continuación el código usado en Matlab y la representación obtenida:
IMAGENPRESIONES
a = 20; z0 = 80; c = 34.91;
V0 = 18;
zref = 10;
alpha = 0.1429;
rho_air = 1.225; %[kg/m3], (dato extraido de worldpress.com)
v = linspace(0, 150, 100);
u = linspace(0, pi, 100);
[U, V] = meshgrid(u, v);
Rho = a * sqrt(1 + ((V - z0) / c).^2);
X = Rho .* cos(U);
Y = Rho .* sin(U);
P_z = 102.79 * V.^0.286;
figure;
surf(X, Y, V, P_z, 'EdgeColor', 'none');
colormap(jet);
colorbar;axis equal;
title('Mapa de Presión del Viento en la Mitad de la Torre Expuesta');
xlabel('X (m)');ylabel('Y (m)');zlabel('Z (m)');
4.2 Campo de fuerza en la superficie expuesta
La fuerza ejercida por el viento sobre la superficie actúa en la dirección normal a ella y el campo vectorial se representa de la siguiente forma:
A continuación se adjunta la imagen con la representación que la fuerza creada por la presión del viento sobre la superficie de la torre, y el código de MATLAB con el que se he representado:
a = 20; c = 34.91; z0 = 80;
H = 150;zref=10;
theta = linspace(0, pi, 18);
z = linspace(0, H, 18);
[Theta, Z] = meshgrid(theta, z);
R = a * sqrt(1 + ((Z - z0).^2 / c^2));
X = R .* cos(Theta);
Y = R .* sin(Theta);
P_z = 102.79 * Z.^0.286
dX_dtheta = -R .* sin(Theta);% Derivadas parciales de la parametrización
dY_dtheta = R .* cos(Theta);
dZ_dtheta = zeros(size(Z));
dX_dz = a * (Z - z0) ./ (c^2 * sqrt(1 + ((Z - z0).^2 / c^2))) .* cos(Theta);
dY_dz = a * (Z - z0) ./ (c^2 * sqrt(1 + ((Z - z0).^2 / c^2))) .* sin(Theta);
dZ_dz = ones(size(Z)); % Derivada de Z respecto de Z
NX = dY_dtheta .* dZ_dz - dZ_dtheta .* dY_dz;% Producto vectorial para obtener el vector normal
NY = dZ_dtheta .* dX_dz - dX_dtheta .* dZ_dz;
NZ = dX_dtheta .* dY_dz - dY_dtheta .* dX_dz;
Norm = sqrt(NX.^2 + NY.^2 + NZ.^2);% Normalización de los vectores normales
NX = NX ./ Norm; NY = NY ./ Norm; NZ = NZ ./ Norm;
% Fuerza descompuesta generada por la presión del aire
FX = -P_z .* NX; FY = -P_z .* NY; FZ = -P_z .* NZ;
%representación del campo vectorial de fuerzas sobre dicha superfie:
quiver3(X, Y, Z, FX, FY, FZ, 1.5, 'r');
xlabel('X (m)'); ylabel('Y (m)'); zlabel('Z (m)');
title('Campo vectorial de fuerza sobre media torre');
axis equal;4.3 Fuerza total y por unidad de superficie del viento sobre la mitad expuesta
Teniendo en cuenta la presión dinámica total calculada previamente P(z), se calcula la fuerza total sobre la mitad frontal partiendo de la expresión (integrando en θ la mitad frontal da un factor 2):
Sustituyendo el perfil de viento y el radio lineal:
queda
Separando integrales:
Las integrales elementales resueltas son:
Por tanto:
Finalmente
Es necesario calcular el área de la mitad de la torre de enfriamiento, la superficie lateral de revolución correspondiente a la mitad frontal es:
Dado que
es constante, el factor geométrico puede sacarse fuera de la integral:
La integral del radio es:
Por tanto:
5 Aplicación practica del hiperboloide de revolución
Tras la investigación hecha anteriormente sobre las ventajas que tiene este tipo de geometría, estructurales y económicas, se presentaran algunos ejemplos del uso de esta geometría en las torres de enfriamiento y de la aplicación en otras estructuras mas del ámbito arquitectónico.
5.1 Ejemplo de torres de enfriamiento
Las torres de enfriamiento de las centrales energéticas son la máxima expresión del hiperboloide de revolución en la construcción. Ejemplos muy representativos son la central nuclear de Cofrentes en España, Didcot en Reino Unido y Niederaussem en Alemania.
5.1.1 Central nuclear de Cofrentes
Esta central se sitúa en la provincia de Valencia. Cuenta con dos torres de enfriamiento de 129 metros de altura. Tiene un diámetro de base de 90 metros y gracias a su funcionamiento ayudan a reducir el consumo de agua del río Júcar.
