Diferencia entre revisiones de «El vórtice de Rankine (Grupo 4)»

De MateWiki
Saltar a: navegación, buscar
(Divergencia y gradiente del campo de velocidades)
(Campo de velocidades)
Línea 36: Línea 36:
 
v_theta_R = 320; % Velocidad tangencial en el límite del núcleo en km/h
 
v_theta_R = 320; % Velocidad tangencial en el límite del núcleo en km/h
 
Gamma = 2 * pi * R * v_theta_R; % Circulación
 
Gamma = 2 * pi * R * v_theta_R; % Circulación
 
 
n = 100; % Número de puntos
 
n = 100; % Número de puntos
 
rho = linspace(0.1, 2*R, n); % Distancia radial en km
 
rho = linspace(0.1, 2*R, n); % Distancia radial en km
Línea 43: Línea 42:
 
x = Mrho .* cos(Mtht); % Coordenadas x
 
x = Mrho .* cos(Mtht); % Coordenadas x
 
y = Mrho .* sin(Mtht); % Coordenadas y
 
y = Mrho .* sin(Mtht); % Coordenadas y
 
 
% Velocidad tangencial
 
% Velocidad tangencial
 
Vtheta = Gamma ./ (2 * pi * R^2) .* Mrho; % General
 
Vtheta = Gamma ./ (2 * pi * R^2) .* Mrho; % General
 
Vtheta(Mrho > R) = Gamma ./ (2 * pi * Mrho(Mrho > R)); % Exterior
 
Vtheta(Mrho > R) = Gamma ./ (2 * pi * Mrho(Mrho > R)); % Exterior
 
 
% Componentes cartesianas
 
% Componentes cartesianas
 
Vx = -Vtheta .* sin(Mtht);
 
Vx = -Vtheta .* sin(Mtht);
 
Vy = Vtheta .* cos(Mtht);
 
Vy = Vtheta .* cos(Mtht);
 
 
% Graficar
 
% Graficar
 
figure;
 
figure;
Línea 58: Línea 54:
 
quiver(x(Mrho > R), y(Mrho > R), Vx(Mrho > R), Vy(Mrho > R), 1, 'b'); % Exterior en azul
 
quiver(x(Mrho > R), y(Mrho > R), Vx(Mrho > R), Vy(Mrho > R), 1, 'b'); % Exterior en azul
 
hold off;
 
hold off;
 
 
axis equal;
 
axis equal;
 
title('Campo de velocidades del vórtice de Rankine (Huracán Camille)');
 
title('Campo de velocidades del vórtice de Rankine (Huracán Camille)');

Revisión del 13:31 27 nov 2024

Trabajo realizado por estudiantes
Título El vórtice de Rankine (Grupo 4)
Asignatura Teoría de Campos
Curso 2024-25
Autores Nacira Faraji Bahja
Paula Gómez Pinilla
Beatriz Matía Esteban
Daniel Portincasa Navarro
Natasha del Carmen Vidal
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción

El modelo de vórtice de Rankine es una representación simplificada de los fluidos rotatorios, utilizada para estudiar fenómenos como huracanes y tornados. Este modelo divide el vórtice en dos regiones: un núcleo central, donde la rotación es uniforme, y una región exterior irrotacional. Aunque idealizado, el modelo es útil para analizar sistemas con núcleos bien definidos, como ciertos huracanes. En este trabajo, se analizarán las propiedades del vórtice de Rankine mediante representaciones gráficas y cálculos matemáticos, usando datos de eventos reales como el huracán Camille (1969). Esto permitirá evaluar su utilidad y limitaciones en la representación de fenómenos atmosféricos complejos.

2 Campo de velocidades

Representar el campo de velocidad en un plano paralelo al suelo en los vórtices de Rankine tiene varias ventajas:
1.)Claridad y simplicidad: Un plano paralelo al suelo permite visualizar claramente las líneas de flujo y la distribución de velocidades, facilitando la comprensión del vórtice.
2.)Consistencia con la geometría: Representar el flujo en este plano se alinea con la naturaleza predominantemente horizontal del vórtice, haciendo los resultados más intuitivos.
3.)Facilidad de análisis: En este plano es más sencillo calcular propiedades como la circulación y la vorticidad, ya que las componentes de velocidad están directamente relacionadas con las ecuaciones del flujo.

En el modelo de Rankine, el campo de velocidad se describe en términos de dos regiones: el “ojo” del vórtice y la región exterior. Para un vórtice con ojo de radio [math]\text{R}[/math] y circulación máxima [math]\Gamma[/math], el campo de velocidad se define en coordenadas cilíndricas [math] \left ( \rho ,\theta ,z \right ) \vec{V}=\upsilon _{\rho} \vec{e}_{\rho }+ \upsilon_{\theta }\vec{e}_{\theta }+\upsilon _{z}\vec{e}_{z}[/math], donde:
[math]v_{r}=0,[/math]

[math]v_{\theta}(r)=\begin{Bmatrix} \frac{\Gamma }{2\pi R^{2}}r & r\leq R \\\frac{\Gamma }{2\pi r} & r\gt R \\ \end{Bmatrix} \\[/math]
[math]v_{z}=0[/math]


Mallado
% Parámetros del huracán Camille
R = 46.3; % Radio del núcleo en km
v_theta_R = 320; % Velocidad tangencial en el límite del núcleo en km/h
Gamma = 2 * pi * R * v_theta_R; % Circulación
n = 100; % Número de puntos
rho = linspace(0.1, 2*R, n); % Distancia radial en km
tht = linspace(0, 2*pi, n); % Ángulo
[Mrho, Mtht] = meshgrid(rho, tht); % Mallado en coordenadas polares
x = Mrho .* cos(Mtht); % Coordenadas x
y = Mrho .* sin(Mtht); % Coordenadas y
% Velocidad tangencial
Vtheta = Gamma ./ (2 * pi * R^2) .* Mrho; % General
Vtheta(Mrho > R) = Gamma ./ (2 * pi * Mrho(Mrho > R)); % Exterior
% Componentes cartesianas
Vx = -Vtheta .* sin(Mtht);
Vy = Vtheta .* cos(Mtht);
% Graficar
figure;
hold on;
quiver(x(Mrho <= R), y(Mrho <= R), Vx(Mrho <= R), Vy(Mrho <= R), 1, 'r'); % Núcleo en rojo
quiver(x(Mrho > R), y(Mrho > R), Vx(Mrho > R), Vy(Mrho > R), 1, 'b'); % Exterior en azul
hold off;
axis equal;
title('Campo de velocidades del vórtice de Rankine (Huracán Camille)');
xlabel('x (km)');
ylabel('y (km)');
legend('Ojo del vórtice (r \leq R)', 'Exterior del vórtice (r > R)');


3 Divergencia y gradiente del campo de velocidades

La divergencia de un campo vectorial mide la tasa neta de flujo que sale de un punto en el espacio. Matemáticamente, la divergencia de un campo de velocidades