El vórtice de Rankine (g.34)

De MateWiki
Saltar a: navegación, buscar
Trabajo realizado por estudiantes
Título El Vórtice de Rankine (Grupo 34)
Asignatura Teoría de Campos
Curso 2025-26
Autores Miguel Gómez-Hidalgo Rivas
Haytam Imhah Chatoual
Darío Pérez
Pablo Ramírez Serrano
Jorge Machín Menés
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción

El Vórtice de Rankine es un modelo matemático, ideado por William John Macquorn Rankine, ingeniero y físico escocés. Su diseño estuvo sujeto a la imperiosa necesidad de explicar de manera simplificada los fluidos rotatorios. Este modelo, aplicado a la vida cotidiana permite la descripción de la estructura básica de fenómenos meteorológicos como tornados y huracanes o en ciertos casos puede explicar ciertos aspectos de la ingeniería como la aerodinámica, ayudando a la creación de sistemas como las turbinas o ventiladores. En este trabajo, haremos algunos cálculos interesantes para la comprensión de este modelo. Además, utilizaremos códigos de Matlab para la representación de funciones y campos vectoriales de manera gráfica.

2 Circulación

Dada la función que representa la velocidad del vórtice [math]v_\theta(\rho) = \begin{cases} \dfrac{\Gamma}{2 \pi R^2} \, \rho & \text{si } \rho \le R \\[2mm] \dfrac{\Gamma}{2 \pi \rho} & \text{si } \rho \gt R \end{cases} [/math], para la situación [math]\rho = \text{R}[/math], tenemos la expresión [math]v_{\theta} = \frac{\Gamma}{2\pi R^{2}}\, \rho \;=\; \frac{\Gamma}{2\pi R} [/math]. Según los datos que nos proporcionan [math](R=250 m; v_{\theta}= 90 m/s )[/math], nos daría un resultado de: [math]\Gamma = 1.4137 \times 10^{5} \ \mathrm{m^2/s} [/math]

Función circulación
% Parámetros
R = 250;           % Radio del núcleo (m)
vR = 90;           % Velocidad tangencial en rho = R (m/s)
Gamma = 2*pi*R*vR; % Circulación total

% Dom
rho = 0:1:1000;    % Distancia radial desde 0 a 1000 m

%Campo de velocidad tangencial según Rankine
vtheta = zeros(size(rho));

inside = (rho <= R);
vtheta(inside) = (Gamma/(2*pi*R^2)) * rho(inside);
outside = (rho > R);
vtheta(outside) = Gamma ./ (2*pi*rho(outside));

vtheta(1) = 0;  % Velocidad nula en el centro

%% Gráfica del perfil de velocidad
figure('Color', 'w', 'Position', [200 200 700 450]);

plot(rho, vtheta, 'b-', 'linewidth', 2);
hold on

plot(R, vR, 'ro', 'MarkerFaceColor', 'r', 'Markersize', 8, ...
    'DisplayName', sprintf('(R, v_R) = (%d m, %d m/s)', R, vR));

xlabel('\rho (m)', 'FontSize', 12);
ylabel('v_\theta(\rho) (m/s)', 'FontSize', 12);
title(sprintf('Vórtice de Rankine - R = %d m, v_\\theta(R) = %d m/s, \\Gamma = %.2e m²/s', ...
    R, vR, Gamma), 'FontSize', 12);

% Leyenda
legend('v_\theta(\rho)', '\rho = R', sprintf('(R, v_R) = (%d, %d)', R, vR), ...
    'Location', 'NorthEast');

% Configuración de ejes
grid on
xlim([0 1000]);
ylim([0 max(vtheta)*1.15]);


% Información en 
fprintf('\n=== VÓRTICE DE RANKINE ===\n');
fprintf('Radio del núcleo: R = %.0f m\n', R);
fprintf('Velocidad en R: v_R = %.0f m/s\n', vR);
fprintf('Circulación: Γ = %.4e m²/s\n', Gamma);
fprintf('Velocidad máxima: v_max = %.2f m/s (en ρ = R)\n', vR);
fprintf('\nComportamiento:\n');
fprintf('- Interior (ρ ≤ R): v_θ ∝ ρ\n');
fprintf('- Exterior (ρ > R): v_θ ∝ 1/ρ\n');


3 Representación de [math]\vec{v}[/math]

Como podemos observar por la expresión del campo de la velocidad, [math]\mathbf{v} = v_\theta(\rho) \, \mathbf{e}_\theta[/math], la velocidad depende de la coordenada radial y no de [math]z[/math]. Además en la ecuación principal de la velocidad, que recordemos:[math]v_\theta(\rho) =\begin{cases} \dfrac{\Gamma}{2 \pi R^2} \, \rho & \text{si } \rho \le R \\[2mm]\dfrac{\Gamma}{2 \pi \rho} & \text{si } \rho \gt R\end{cases}[/math] , se puede apreciar cómo en sí misma no se puede apreciar ninguna [math]z[/math]. Es por ello que el campo vectorial de v es estrictamente horizontal, y está “conviviendo” en un plano bidimensional. Además al no tener la z, la estructura del campo vectorial es la misma a cualquier altura z.

Campo vectorial de [math]\vec{v}[/math]
R = 250;         % r del núcleo (m)
vR = 90;         % v tang en rho=R (m/s)
Gamma = 2*pi*R*vR;
% Dom
[x, y] = meshgrid(-800:40:800, -800:40:800);
rho = sqrt(x.^2 + y.^2);        % dist rad
theta = atan2(y, x);            % áng polar
% Campo de v_theta
vtheta = zeros(size(rho));
% Región interna flujo como sol
inside = rho <= R & rho ~= 0;
vtheta(inside) = (Gamma./(2*pi*R^2)) .* rho(inside);
% Región externa irrotacional
outside = rho > R;
vtheta(outside) = Gamma ./ (2*pi*rho(outside));
%%Componentes cartesianas de la velocidad
% Direcciones tangenciales
u = -vtheta .* sin(theta);   % componente en x
v =  vtheta .* cos(theta);   % componente en y
% Gráfica quiver
figure('Color','w','Position',[200 200 700 550])
hold on
% Parte interna del vórtice (color azul)
quiver(x(inside), y(inside), u(inside), v(inside), ...
    'Color',[0 0.3 1], 'LineWidth',1.2)
% Parte externa del vórtice (color rojo)
quiver(x(outside), y(outside), u(outside), v(outside), ...
    'Color',[1 0 0], 'LineWidth',1.2)
% Círculo marcando el ojo
th = linspace(0,2*pi,300);
plot(R*cos(th), R*sin(th), 'k--', 'LineWidth',1.4)

axis equal
xlabel('x (m)')
ylabel('y (m)')
title('Campo vectorial del vórtice de Rankine (plano horizontal)')
legend('Interior del ojo','Exterior del vórtice','Radio del núcleo')

grid on


4 Vórtices atmosféricos y de Burger-Rott

Uno de los vórtices más conocidos que existen son los tornados, estos son columnas de aire que giran violentamente y se extienden desde las nubes hasta el suelo. Este tipo de vórtice meteorológico se mide mediante la Escala Fujita Mejorada que mide desde EF0 hasta EF5 en función del daño que el tornado causa. Un tornado de intensidad EF5 puede alcanzar vientos de hasta 322 km/h. La formación de un tornado ocurre cuando el aire cálido y húmedo asciende y choca con el frío y seco, creando así inestabilidad atmosférica, acompañada de una tormenta conocida como supercélula. El cambio de dirección del viento provoca que en la parte superior de el tornado tenga corrientes que giran horizontalmente. A medida que se acerca al suelo, estas corrientes se giran e inclinan creando un efecto embudo que arrastra tierra y escombros.

Otro de los vórtices más conocidos son los huracanes, o ciclones tropicales. La intensad de este tipo de vórtices se mide mediante la Escala Saffir-Simpson, que clasifica los huracanes según el daño civil que pueden provocar, pudiendo sobrepasar los 252 km/h. Para que el huracán se forme, necesita agua de mar cálida, aire cálido y húmedo. El aire húmedo, asciende, se enfría y libera calor, creando una zona de baja presión que atrae más aire aumentando los vientos. Debido a la rotación de la tierra, este movimiento se transforma en un patrón espiral formándose los ciclones tropicales.

Otro ejemplo es la tromba marina, esta es un tornado que se forma sobre el agua. Estos tornados de agua se clasifican por la antes mencionada Escala Fujita Mejorada, aunque su intensidad no suele superar los 138 km/h (conocidas como no tornádicas), aunque podrían llegar a superar los 225 km/h. Estas se forman debido a la interacción entre la convección del aire mediante un movimiento ascendente y la rotación sostenida de aire.

Para terminar, los dust devils como su nombre indica, son remolinos de polvo o arena, ligeramente semejantes a los tornados. A este vórtice no se le asemeja como tal a una escala, aunque sin embargo, gracias a su semejanza con los tornados, se les puede asignar como al vórtice de antes, la Escala Fujita Mejorada. Pese a la similitud mencionada, los dust devils son mucho más débiles y menos intensos, alcanzando rara vez un EF1 según su escala. Estos remolinos de polvo se forman cuando una corriente de aire caliente se eleva y comienza a girar. Además, el terreno ha de ser seco y con pequeñas piedras o arena.

La alternativa más realista del vórtice de Rankine es la del vórtice de Burger-Rott. Ambos nos describen vórtices con un núcleo central dónde la velocidad es proporcional a r y que además tiene una región exterior donde disminuye la velocidad. Sin embargo, el modelo de Burger-Rott describe una transición entre el paso del núcleo a la región exterior que es suave, mientras que el de Rankine nos describe una transición bastante brusca. En ambos modelos, el fluido que se presenta tiende a ir hacia dentro, como podemos apreciar en el lavabo cuando taponamos el desagüe y al abrirlo se va moviendo hacia dentro. A parte, este fluido en los dos modelos va hacia arriba o hacia abajo en función de el eje, que es la línea imaginaria ubicada en el centro del vórtice.

5 Divergencia de [math]\vec{v}[/math]

Dada la ecuación de la velocidad (la podemos encontrar en el apartado 2 y 3), se puede apreciar que nos describe la función de movimiento mediante coordenadas cilíndricas. Para calcular la divergencia del campo, debemos aplicar la fórmula [math]\nabla \cdot \vec{v} = \frac{1}{\rho}\frac{\partial}{\partial \rho}(\rho v_{\rho}) + \frac{1}{\rho}\frac{\partial}{\partial \theta}(v_{\theta})+ \frac{\partial}{\partial z}(v_{z})[/math].Además, podemos observar cómo [math]v_{\rho} = 0; v_{z} = 0; v_{\theta} = v_{\theta}(\rho)[/math]. Por lo tanto; sustituyendo en la fórmula, nos quedaría:[math]\nabla \cdot \vec{v} = \left( \frac{1}{\rho}\frac{\partial}{\partial \rho}(\rho \cdot 0)+ \frac{1}{\rho}\frac{\partial}{\partial \theta}\big(v_{\theta}(\rho)\big) + \frac{\partial}{\partial z}(0) \right)= 0 + \frac{1}{\rho}\cdot 0 + 0 = 0.[/math].

Con esto concluimos que la divergencia del campo de velocidad es nulo. Además como ya sabemos, la divergencia de un campo nos indica si el campo se extiende o se contrae en su dominio. O dicho en otras palabras la tasa de cambio volumétrica por unidad de volumen. Por eso, podemos decir que si la divergencia de nuestro campo velocidad es 0 es porque no existe cambio de volumen. Además, por definición, el flujo es incompresible, esto es no hay fuentes (punto donde “sale” el fluido) ni sumideros (punto donde “entra” el fluido).

6 Calcula analíticamente el rotacional y representarlo gráficamente

El rotacional es un operador diferencial que mide la rotación de un campo vectorial alrededor de un punto. Antes de proceder al cálculo rotacional también hay que hacer hincapié en el significado de vorticidad. La vorticidad es el rotacional de un campo vectorial de velocidad v:[math]\vec{\zeta}=\nabla \times \vec{v}=\left(\frac{1}{r}\frac{\partial v_z}{\partial \theta}-\frac{\partial v_\theta}{\partial z}\right)\hat{e}_r+\left(\frac{\partial v_r}{\partial z}-\frac{\partial v_z}{\partial r}\right)\hat{e}_\theta+\frac{1}{r}\left(\frac{\partial (r v_\theta)}{\partial r}-\frac{\partial v_r}{\partial \theta}\right)\hat{e}_z[/math]

En este caso el rotacional solo tiene una componente[math]v_z[/math] . Por lo tanto, se reduce a: [math]\nabla \times \vec{v}=\frac{\Gamma}{2\pi r^2}\hat{e}_z[/math]

Además, la vorticidad es inversamente proporcional a la distancia, es decir, cuanto más lejos nos encontremos del vórtice menor será la intensidad. Para hacer el cálculo rotacional se hace por separado y para facilitar el cálculo lo hacemos todo en coordenadas cilíndricas:

1.Núcleo (r<R): Principalmente sacamos la velocidad angular en la que en el núcleo se determina a partir de la siguiente fórmula: [math]v_\theta(r)=\omega r[/math] Despejando la velocidad angular y partiendo de los datos proporcionados por el profesor la velocidad angular sería: [math]\omega=\frac{90\,\text{m/s}}{250\,\text{m}}=0,36\,\text{s}^{-1}[/math]. En el núcleo el rotacional es constante por tanto aplicamos la siguiente fórmula: [math]\nabla\times\vec{v}=2\omega\hat{e}_z=2\times 0,36\,\text{s}^{-1}\hat{e}_z=0,72\,\text{s}^{-1}\hat{e}_z[/math]. Por tanto, la vorticidad máxima y uniforme en el núcleo del tornado es [math]0,72\,\text{s}^{-1}[/math].

2.Región exterior (r>R): En la región exterior el vórtice es irrotacional por tanto la vorticidad fuera del núcleo es: [math]\nabla\times\vec{v}=\vec{0}[/math]

Representación gráfica: Perfil de Vortizidad del Vórtice de Rankine (Tornado EF4).
% 1. Definición de Parámetros y Constantes
R = 250;           % Radio del núcleo (m)
V_R = 90;          % Velocidad tangencial máxima en r = R (m/s)

% Cálculo de la velocidad angular (omega) y la vorticidad del núcleo (zeta_core)
omega = V_R / R;
zeta_core = 2 * omega; % Magnitud del Rotacional en el núcleo (s^-1)

% 2. Dominio Radial (Eje x del gráfico)
% Creamos un vector de distancias radiales (r)
R_max = 4 * R;
N = 500;           % Número de puntos
r = linspace(0, R_max, N); % Vector de distancias radiales desde el centro

% 3. Cálculo del Rotacional (Vorticidad)
% Inicializar el vector de vorticidad (Rotacional) a cero (para r > R)
Zeta = zeros(1, N);
Indices_core = (r <= R);
Zeta(Indices_core) = zeta_core;

% 4. Representación Gráfica del Perfil
figure;
plot(r, Zeta, 'r', 'LineWidth', 3); % Trazar el perfil de vorticidad

% Marcar el límite del núcleo
hold on;
plot([R R], [0 zeta_core] * 1.1, 'k--', 'LineWidth', 1.5); % Línea vertical en r=R

% Añadir etiquetas y título
xlabel('Distancia Radial r (m)');
ylabel(['Magnitud del Rotacional (Vorticidad) ||\nabla \times \vec{v}|| (s^{-1})']);
title('Perfil de Vortizidad del Vórtice de Rankine (Tornado EF4)');

% Añadir grid y leyenda
grid on;
legend('Magnitud del Rotacional', ['Límite del Núcleo (r = ' num2str(R) ' m)'], 'Location', 'NorthEast');

% Ajustar el eje y para ver claramente el salto
ylim([0 zeta_core * 1.1]);
hold off;



7 Mapa de colores del campo escalar [math]\vert \nabla \times \vec{v} \vert[/math]

La vorticidad está concentrada en el núcleo, es decir la circunferencia de radio [math]250\,\text{m}[/math] representada anteriormente por Matlab con color rojo. Y en la región exterior no existe vorticidad por lo tanto esta coloreada por azul, siendo el rojo la máxima vorticidad y el azul la mínima es decir [math]0[/math].

En matlab sería:

Mapa de colores de la Magnitud de la Vorticidad, donde el rojo representa el valor máximo ($0,72\text{s}^{-1}$) y el azul representa la vorticidad nula.
% 1. Parámetros Requeridos
R_ = 250;           % Radio del núcleo (m)
V_R = 90;           % Velocidad tangencial máxima en r = R (m/s)

% Cálculo de la velocidad angular (omega) y la vorticidad del núcleo
omega = V_R / R;
zeta_core = 2 * omega; % 0.72 s^-1

% 2. Creación del Dominio (Plano Horizontal)
R_max = 5 * R; % Rango hasta 1250 m del centro
N = 500;       % Número de puntos
x = linspace(-R_max, R_max, N);
y = linspace(-R_max, R_max, N);
[X, Y] = meshgrid(x, y);

% Calcular la distancia radial r para cada punto
R_grid = sqrt(X.^2 + Y.^2);

% 3. Cálculo del Campo Escalar (Magnitud de la Vorticidad)
% Inicializar la matriz de vorticidad (Zeta) a cero (para la región exterior)
Zeta = zeros(N, N);
% Encontrar los puntos que están dentro del núcleo (r <= R)
indices_core = (R_grid <= R);
% Asignar la vorticidad constante en el núcleo (el valor máximo)
Zeta(indices_core) = zeta_core;

% 4. Representación Gráfica (Mapa de Colores)
figure;
imagesc(x, y, Zeta); % Crea el mapa de colores
axis equal tight;    % Asegura que el círculo no se vea deformado
colorbar;            % Muestra la barra de color con la escala de vorticidad
colormap jet;        % Uso de un mapa de colores común (jet o otro de su elección)
xlabel('Coordenada x (m)');
ylabel('Coordenada y (m)');
title('Mapa de Colores de la Magnitud de la Vorticidad ||\nabla \times \vec{v}||');
set(gca, 'YDir', 'normal'); % Asegura que Y aumente hacia arriba

% Opcional: Marcar el límite del núcleo (R = 250 m)
hold on;
theta = linspace(0, 2*pi, 100);
plot(R*cos(theta), R*sin(theta), 'k--', 'LineWidth', 2);
legend('Límite del Núcleo');
hold off;


8 Comportamiento de una barca respecto al vórtice según su posición

Dependiendo de la zona en que se encuentre un objeto tendrá diferente comportamiento. En este caso comentaremos sobre el comportamiento de una barca según su posición.

1.En el núcleo del vórtice el flujo se comporta de una manera parecida al de un sólido en rotación ya que hay una velocidad tangencial que provoca una rotación uniforme, y dicha velocidad tangencial disminuye linealmente con la distancia al núcleo. Por lo tanto, si la barca está en el núcleo rota sobre si misma ya que existe una vorticidad o rotación que hace que todo el fluido situado en el núcleo rote como si fuera una sola pieza.

2.En cambio, en la región exterior el flujo se comporta de manera irrotacional, es decir, la vorticidad es nula por lo tanto la barca en el exterior no estaría sometida a ningún tipo de rotación.

9 Campo de presión p(ρ, z) y su representación

El campo de presión relaciona la posición de un punto en el espacio con el valor de la presión en ese punto. Se define en el espacio p(\rho, z) tal que:

  • [math]\rho \in [0, 1000][/math]
  • [math]z \in [0, z_0][/math]

Siendo la formula de presión: [math] p(\rho, z) = \begin{cases} P_0 + \frac{1}{2} \rho_{aire} \cdot v_{\theta}^2(\rho) - \rho_{aire} \cdot g \cdot z & \text{si } \rho \le R \\ P_{\infty} - \frac{1}{2} \rho_{aire} \cdot v_{\theta}^2(\rho) - \rho_{aire} \cdot g \cdot z & \text{si } \rho \gt R \end{cases} [/math] Donde la velocidad tangencial [math]v_{\theta}[/math] se define como: [math] v_{\theta} = \begin{cases} \frac{\Gamma}{2\pi R^2} \rho & \text{si } \rho \le R \\ \frac{\Gamma}{2\pi \rho} & \text{si } \rho \gt R \end{cases} [/math]

9.1 Código MATLAB para el Campo de Presión del Vórtice

Mapa de colores del campo de presión p(\rho, z) en un vórtice de Rankine.
% Datos
R = 250;              % Radio (m)
v_R = 90;             % Velocidad tangencial en el Radio R (m/s)
p_0 = 920* 100;       % Presión mínima en el centro (Pa)
z_0 = 2.8 * 1000;     % Altura máxima del dominio (m)
rho_aire = 1.225;     % Densidad aire (kg/m^3)
g = 9.81;             % Aceleración gravitacional (m/s^2)
p_inf = 101325;       % Presión atmosférica estándar (Pa)
Gamma = 2 * pi * R * v_R;   % Circulación (m^2/s)

% Función para la Velocidad Tangencial (v_theta)
v_theta = @(rho) arrayfun(@(r) ...
    (r <= R) * ((Gamma / (2 * pi * R^2)) * r) + ... % Núcleo
    (r > R) * ((Gamma / (2 * pi * r))), rho);       % Exterior

% Definición del Dominio
rho_max = 1000;     % Radio máximo del mapa (m)
N_rho = 101; 
N_z = 51;   
rho_vec = linspace(0, rho_max, N_rho);
z_vec = linspace(0, z_0, N_z);
[RHO, Z] = meshgrid(rho_vec, z_vec); % Crear la malla 2D

% Cálculo del Campo de Presión
V_THETA = v_theta(RHO);
P = zeros(size(RHO));

% 1. Región del Núcleo (rho <= R)
mascara_nucleo = (RHO <= R);
V_NUCLEO = V_THETA(mascara_nucleo);
% Ecuación para el núcleo (usando p_0):
P(mascara_nucleo) = p_0 + 0.5 * rho_aire * V_NUCLEO.^2 - rho_aire * g * Z(mascara_nucleo);

% 2. Región Exterior (rho > R)
mascara_ext = (RHO > R);
V_EXTERIOR = V_THETA(mascara_ext);
% Ecuación para el exterior (usando p_inf):
P(mascara_ext) = p_inf - 0.5 * rho_aire * V_EXTERIOR.^2 - rho_aire * g * Z(mascara_ext);

% 3. Representación Gráfica (Mapa de colores)
figure;
contourf(RHO, Z, P / 100, 20, 'LineStyle', 'none'); % Conversión a mbar (P/100)
colorbar;
colormap(jet);
xlabel('Radio \rho (m)');
ylabel('Altura z (m)');
title(['Presión p(\rho, z) Vórtice de Rankine (Tornado EF4)']);
c = colorbar;
c.Label.String = 'Presión (mbar)';
axis equal;
ylim([0, z_0]);
grid on;

10 Cálculo de la Caída de Presión (Δ𝑝)

Se calcula la caída de presión entre el exterior y el centro del ojo del vórtice al nivel del suelo ([math]z=0[/math]) y se compara con el valor dado [math]P_{\infty} - P_0[/math].

La caída de presión pedida es: [math]\Delta p = p(R^{+}, 0) - p(0, 0)[/math] . Donde [math]R^{+}[/math] es el radio inmediatamente fuera del núcleo del vórtice.

Presión en el centro del ojo: [math]p(0, 0)[/math]

El punto está en el núcleo, donde [math]\rho = 0[/math] y [math]z = 0[/math].

Como [math]\rho=0[/math], entonces [math]v_{\theta}(\rho)=0[/math] y la presión queda:

[math]p(0, 0) = P_0 + \frac{1}{2} \rho_{aire} \cdot v_{\theta}^2(0) - \rho_{aire} \cdot g \cdot 0[/math] [math]p(0, 0) = P_0[/math]

Presión justo fuera del núcleo: [math]p(R^{+}, 0)[/math]

El punto está en la región exterior, donde [math]\rho = R[/math] y [math]z = 0[/math].

Como el campo de velocidades es continuo en [math]\rho = R[/math] para el vórtice de Rankine, entonces [math]v_{\theta}(R^{+}) = v_{\theta}(R)[/math], y conocemos que: [math]v_{\theta}(R) = \frac{\Gamma}{2 \pi R}[/math]

Entonces, sustituyendo en la fórmula exterior: [math]p(R^{+}, 0) = P_{\infty} - \frac{1}{2} \rho_{aire} \cdot v_{\theta}^2(R) - \rho_{aire} \cdot g \cdot 0[/math]

[math]p(R^{+}, 0) = P_{\infty} - \frac{1}{2} \rho_{aire} \left( \frac{\Gamma}{2 \pi R} \right)^2[/math] Sustituyendo lo obtenido: [math]\Delta p = p(R^{+}, 0) - p(0, 0)[/math] [math]\Delta p = \left[ P_{\infty} - \frac{1}{2} \rho_{aire} \left( \frac{\Gamma}{2 \pi R} \right)^2 \right] - P_0[/math]


Simplificando: [math]\Delta p = (P_{\infty} - P_0) - \frac{1}{2} \rho_{aire} V_R^2[/math]

Usando los datos [math]\frac{1}{2} \rho_{aire} V_R^2 = \frac{1}{2} (1.225) (90)^2 = \frac{1}{2} (1.225)(8100) = 4961.25 \text{ Pa}[/math]

Sustituyendo los valores conocidos ([math]P_{\infty} = 101325 \text{ Pa}[/math], [math]P_0 = 92000 \text{ Pa}[/math]): [math]\Delta p = 101325 - 92000 - 4961.25 = 4363.75 \text{ Pa}[/math]

Comparado con [math]P_{\infty} - P_0 = 9325 \text{ Pa}[/math] no coinciden.

El valor [math]P_{\infty} - P_0[/math] representa la caída total de presión desde la atmósfera en reposo ([math]v_{\theta}=0[/math]) hasta el centro del ojo donde la velocidad es cero. En cambio, la caída calculada ([math]\Delta p[/math]) representa solo la caída de presión entre el punto de velocidad máxima (justo fuera del núcleo) y el centro del ojo.

11 Diferencia de presión entre la presión atmosférica estándar y en el centro del ojo. Limitaciones del modelo

Piden la diferencia de presión total: [math]P_{\infty} - p(0, 0)[/math]

Asumiendo que la presión es una función continua, la caída de presión total [math]P_{\infty} - p(0, 0)[/math] se obtiene al aplicar la ecuación de Bernoulli simplificada entre un punto muy lejano ([math]\rho \to \infty[/math]) y el centro del ojo ([math]\rho=0[/math]), ambos en el suelo ([math]z=0[/math]).

Presión en el Centro del Ojo [math]p(0, 0)[/math]

Sustituyendo en la fórmula de presión para el núcleo ([math]\rho \le R[/math]) con [math]\rho=0[/math] y [math]z=0[/math]: [math]p(0, 0) = P_0[/math]

Presión en el Exterior [math]p(\infty, 0)[/math]

El punto está en la región exterior ([math]\rho \gt R[/math]). Cuando [math]\rho \to \infty[/math], entonces [math]v_{\theta}(\rho) \to 0[/math]. Sustituyendo en la fórmula exterior, queda: [math]p(\infty, 0) = P_{\infty}[/math]

Derivación Teórica de la Caída de Presión

Para que el modelo sea consistente y [math]P_0[/math] sea la presión mínima real, la caída de presión teórica debe ser la suma de la caída en la región exterior y la caída en la región del núcleo. La mayor parte de la caída de presión ocurre en el límite [math]\rho=R[/math].

Sabemos que la presión debe ser continua en [math]\rho = R[/math]: [math]p(R, 0) = P_{\infty} - \frac{1}{2} \rho_{aire} v_{\theta}^2(R)[/math] [math]p(R, 0) = P_0 + \frac{1}{2} \rho_{aire} v_{\theta}^2(R)[/math]

Igualando estas dos expresiones y despejando la diferencia de presión teórica: [math]P_{\infty} - \frac{1}{2} \rho_{aire} v_{\theta}^2(R) = P_0 + \frac{1}{2} \rho_{aire} v_{\theta}^2(R)[/math]

Despejamos la diferencia de presión teórica y queda: [math]\Delta P_{\text{Teórico}} = P_{\infty} - P_0 = \rho_{aire} v_{\theta}^2(R)[/math]

Usamos los parámetros:

  • [math]\rho_{aire} = 1.225 \text{ kg/m}^3[/math]
  • [math]v_{\theta}(R) = 90 \text{ m/s}[/math]

El valor teórico de la caída de presión es: [math]\Delta P_{\text{Teórico}} = 1.225 \cdot (90)^2 = 9922.5 \text{ Pa}[/math]

[math]\Delta P_{\text{Teórico}} = 99.23 \text{ mbar}[/math]

El valor dado de la caída de presión total es: [math]P_{\infty} - P_0 = 9325 \text{ Pa}[/math] [math]P_{\infty} - P_0 = 93.25 \text{ mbar}[/math]

La diferencia numérica entre el valor teórico y el valor dado es: [math]\text{Diferencia} = \Delta P_{\text{Teórico}} - (P_{\infty} - P_0)_{\text{Dado}}[/math] [math]\text{Diferencia} = 99.23 \text{ mbar} - 93.25 \text{ mbar} = 5.98 \text{ mbar}[/math]

La diferencia de [math]5.98 \text{ mbar}[/math] entre el valor teórico y el valor dado es aceptable.

11.1 Limitaciones del Modelo de Rankine

La magnitud de un tornado de categoría EF4 es extrema, y el modelo de Rankine es una idealización simplificada. El error en la estimación de la caída de presión máxima es razonable, ya que los valores de [math]P_0[/math] y [math]V_R[/math] vienen de mediciones sujetas a incertidumbre.

Las limitaciones de este modelo son:

  • Ignora por completo los efectos de la compresibilidad y las turbulencias reales, asumiendo un flujo perfectamente simétrico.
  • Transición Discontinua: La suposición más fuerte es la transición brusca en [math]\rho=R[/math]. En un tornado real, la viscosidad y la turbulencia crean una región de transición suave.
  • Ignora la fricción con el suelo y con el aire circundante, así como el flujo radial y el flujo vertical que son esenciales para la alimentación de energía en un tornado real.
  • Los fenómenos reales están impulsados por la convección y los gradientes de temperatura, que el modelo de Rankine ignora.

12 Calculo del gradiente de presión

Para el cálculo del gradiente de presión utilizaremos la expresión:

[math]\frac{dp}{dr} = \rho_{\text{aire}} \frac{v_\theta^2}{r}[/math]

que es la fórmula del movimiento radial.

Se definen los parámetros, la malla y la distancia radial \(r\) para cada punto de dicha malla. Se crea la matriz base y se aplican las ecuaciones de la velocidad tangencial:

En el núcleo:

[math]v_\theta(r) = \tfrac{1}{2} \rho_{\text{aire}} \Omega^2 r^2[/math]

En el exterior del núcleo:

[math]v_\theta(r) = \frac{\Gamma}{2\pi r}[/math]

A continuación, se calcula el módulo del gradiente de presión y se realiza la descomposición en sus componentes \(X\) y \(Z\). 

ro = 1.2;      % densidad del aire (kg/m^3)
R = 200;        % radio del núcleo (m)
Omega = 0.3;    % velocidad angular (rad/s)
Gamma = Omega*R^2;  % circulación

% Mallado para la sección vertical (x,z)
[x,z] = meshgrid(-400:20:400, -500:20:0);
r = sqrt(x.^2 + z.^2);
r(r==0) = 1e-6;  % evita poder dividir por cero

% Velocidad tangencial
veltheta = zeros(size(r));
veltheta(r <= R) = Omega .* r(r <= R);
veltheta(r > R)  = Gamma ./ (2*pi*r(r > R));

% Gradiente de presión
Grad = ro .* (veltheta.^2 ./ r);

% Componentes de la seccion(x,z)
Cx = Grad .* (x ./ r);
Cz = Grad .* (z ./ r);

% Creacióon de la figura
figure;
quiver(x, z, Cx, Cz, 'b');
xlabel('x (m)')
ylabel('z (m)')
title(' \nabla p en una sección vertical del vórtice de Rankine')
axis equal
grid on


12.1 Representación de superficies isobáricas

Nos dan los valores de presión: 950, 970, 990 y 1000 mbar.

Al estar trabajando con un vórtice de Rankine, la presión solo depende del radio, de modo que el radio es la única variable. Por tanto, las isobaras son círculos concéntricos alrededor del eje del vórtice.

El vórtice se divide en dos regiones: presión dentro del núcleo y presión fuera del núcleo.

1. Análisis dentro del núcleo (r ≤ R):

En esta región, la velocidad angular \(\Omega\) es constante. Aplicando la ecuación de movimiento radial obtenemos:

[math]p(r) = p(0) + \tfrac{1}{2} \rho_{\text{aire}} \Omega^2 r^2[/math]

siendo \(p(0)\) la presión en el centro del vórtice.

Dado que \(1\ \text{mbar} = 100\ \text{Pa}\), la conversión deja la expresión final:

[math]p(r \le R) = p(0) + \tfrac{1}{2} \rho_{\text{aire}} \Omega^2 r^2 \cdot 0.01[/math]

2. Análisis fuera del núcleo (r > R):

En esta región el flujo es irrotacional y la circulación \(\Gamma\) es constante. La velocidad tangencial viene dada por:

[math]v_\theta(r) = \frac{\Gamma}{2\pi r}[/math]

Al aplicar la ecuación de movimiento radial y realizar la conversión de unidades obtenemos:

[math]p(r \gt R) = p(0) + \tfrac{1}{2} \rho_{\text{aire}} \left( \frac{\Gamma}{2\pi r} \right)^2 \cdot 0.01[/math]

Archivo:Campovectarea2
Campo vectorial de [math]\vec{v}[/math]
% Parámetros del vórtice de Rankine
ro = 1.2;    % kg/m^3
R = 200;     % m
Omega = 0.3; % rad/s
Gamma = Omega * R^2; % <-- Parámetro auxiliar usado en la fórmula A

% Presión ambiental
p0 = 1010;   % valor usual exterior a un tornado

% Mallado horizontal (x,y)
[x, y] = meshgrid(-500:5:500, -500:5:500);
r = sqrt(x.^2 + y.^2);
r(r==0) = 1e-6; % Evitar división por cero

%Presión en el núcleo (r<=R)
p = zeros(size(r));
p(r<= R) = p0 - 0.5 * ro * (Omega^2) * (r(r<= R)).^2 * 0.01; 
%Presión fuera del núcleo (r > R) 
pR = p0 - 0.5 * ro * Omega^2 * R^2 * 0.01; % Presión en el borde r=R
A = ro * Gamma^2 / (8*pi^2) * 0.01;        % A usa el Gamma auxiliar
p(r > R) = pR - A * ( (1./(r(r > R)).^2) - (1/R^2) ) * 0.01; 

% Representación de isóbaras
figure;
contour(x, y, p, [950 970 990 1000],'LineColor', 'r', 'LineWidth', 2);
clabel(contour(x,y,p,[950 970 990 1000]),'FontSize', 10);

xlabel('x (m)');
ylabel('y (m)');
title('Isóbaras del vórtice de Rankine');
axis equal;
grid on;


12.2 Estimación de la fuerza neta hacia el centro de la fachada de una pequeña edificación

La depresión en el núcleo del tornado genera fuerzas destructivas. Dicha edificación presenta un área \(A = 50\,\text{m}^2\), situada a una distancia \(\rho = \tfrac{3}{4}R\) del eje. La fachada está orientada perpendicularmente al flujo radial del tornado. Expresa el resultado en toneladas-fuerza (1 tf ≈ 10 kN).

Consideramos un vórtice de Rankine con núcleo sólido para \(r \le R\), cuya velocidad tangencial es:

[math]v_\theta(r) = \Omega r[/math]

La ecuación de equilibrio radial queda:

[math]\frac{dp}{dr} = \rho_{\text{aire}} \frac{v_\theta^2}{r}[/math]

Sustituyendo la velocidad tangencial:

[math]\frac{dp}{dr} = \rho_{\text{aire}} \frac{(\Omega r)^2}{r} = \rho_{\text{aire}} \Omega^2 r[/math]

Para hallar la depresión \(\Delta p\), integramos desde el exterior \(p(0)\) hasta un radio \(r\):

[math] p(0) - p(r) = \Delta p = \int_0^{r} \rho_{\text{aire}} \Omega^2 r' \, dr' = \tfrac{1}{2} \rho_{\text{aire}} \Omega^2 r^2 [/math]

La fuerza sobre la fachada, al estar orientada perpendicularmente al flujo radial, es:

[math]F = \Delta p \, A[/math]

Sustituimos los parámetros proporcionados:

[math]\rho_{\text{aire}} = 1.2\ \text{kg/m}^3,\qquad \Omega = 0.3\ \text{rad/s},\qquad R = 200\ \text{m},\qquad r = 0.75\,R = 150\ \text{m},\qquad A = 50\ \text{m}^2 [/math]

Cálculos intermedios:

[math]r^2 = 150^2 = 22500\ \text{m}^2[/math]

[math]\Omega^2 = (0.3)^2 = 0.09[/math]

Sustituyendo en la expresión de la depresión:

[math] \Delta p = \tfrac{1}{2} \cdot 1.2 \cdot 0.09 \cdot 22500 = 1215\ \text{Pa} [/math]

Relación: \(1\ \text{Pa} = 1\ \text{N/m}^2\)

[math]F = 1215 \times 50 = 60750\ \text{N}[/math]

Conversión a toneladas-fuerza (1 tf ≈ 10000 N):

[math]F_{\text{tf}} = \frac{60750}{10000} = 6.075\ \text{tf} \approx 6.08\ \text{tf}[/math]

Por tanto, la fuerza neta hacia el centro sobre la fachada es aproximadamente:

[math]F \approx 6.08\ \text{toneladas-fuerza}\ \ (\approx 60.8\ \text{kN})[/math]

Obtenido de «https://mat.caminos.upm.es/w/index.php?title=El_vórtice_de_Rankine_(g.34)&oldid=96985»