Diferencia entre revisiones de «La espiral de Ekman (Grupo 16)»

De MateWiki
Saltar a: navegación, buscar

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434
(→‎Campo vectorial V en varios planos horizontales)
(→‎Campo vectorial V en varios planos horizontales)
Línea 83: Línea 83:
  
 
==Campo vectorial V en varios planos horizontales==
 
==Campo vectorial V en varios planos horizontales==
−
{{matlab|codigo=%%
+
{{matlab|codigo=  
 
Apartado 5
 
Apartado 5
 +
 
clc;close all;clear;
 
clc;close all;clear;
−
omega = 7.2921e-5;  % Velocidad angular de la Tierra (rad/s)
+
omega = 7.2921e-5;  % Velocidad angular (rad/s)
−
phi = 30.10242;  % Latitud en grados
+
phi = 30.10242;  % Latitud (grados)
 
V0 = 0.15;  % Velocidad superficial inducida por el viento (m/s)
 
V0 = 0.15;  % Velocidad superficial inducida por el viento (m/s)
 
nu = 0.05;  % Viscosidad turbulenta (m^2/s)
 
nu = 0.05;  % Viscosidad turbulenta (m^2/s)
−
theta0 = 3*pi/4;  % Ángulo inicial en radianes
+
theta0 = 3*pi/4;  % Ángulo inicial (radianes)
−
% Calcular el parámetro de Coriolis f
+
 
 +
% Cálculo del parámetro de Coriolis f
 
f = 2 * omega * sin(phi);
 
f = 2 * omega * sin(phi);
−
% Calcular la profundidad de Ekman
+
 
 +
% Profundidad de Ekman
 
delta = sqrt(2 * nu / abs(f));
 
delta = sqrt(2 * nu / abs(f));
 +
 
z_max = 3 * delta;  % Profundidad máxima para la animación
 
z_max = 3 * delta;  % Profundidad máxima para la animación
 
n_frames = 40;  % Número de frames de la animación
 
n_frames = 40;  % Número de frames de la animación
 
z_vals = linspace(0, z_max, n_frames);  % Valores de profundidad
 
z_vals = linspace(0, z_max, n_frames);  % Valores de profundidad
−
[X, Y] = meshgrid(linspace(-1, 1, 10), linspace(-1, 1, 10));  % Malla para visualizar el campo vectorial
+
[X, Y] = meshgrid(linspace(-1, 1, 10), linspace(-1, 1, 10));   
 +
 
 
figure(1);
 
figure(1);
 
view(3)
 
view(3)
 
%axis equal;
 
%axis equal;
−
% Calcular las componentes de la velocidad
+
 
 +
% Calculo velocidad
 
u_m = (f/abs(f))*V0 * exp(-z_vals / delta)  .* cos(-z_vals / delta+theta0);  % Componente u(z)
 
u_m = (f/abs(f))*V0 * exp(-z_vals / delta)  .* cos(-z_vals / delta+theta0);  % Componente u(z)
−
v_m = V0 * exp(-z_vals / delta)  .* sin(-z_vals / delta+theta0);  % Componente v(z)
+
v_m = V0 * exp(-z_vals / delta)  .* sin(-z_vals / delta+theta0);   
 +
 
 +
% Componente v(z)
 
plot3(u_m,v_m,-z_vals,'Color','g')
 
plot3(u_m,v_m,-z_vals,'Color','g')
 
hold on
 
hold on
−
% Inicializar la animación
+
 
 +
% Animación
 
for k = 1:n_frames
 
for k = 1:n_frames
 
   z = z_vals(k);  % Profundidad actual
 
   z = z_vals(k);  % Profundidad actual
 
   u = u_m(k);  % Componente u(z)
 
   u = u_m(k);  % Componente u(z)
 
   v = v_m(k);  % Componente v(z)
 
   v = v_m(k);  % Componente v(z)
−
   % Calcular los campos vectoriales para este valor de z
+
 
 +
   % Calculo de campos vectoriales para este valor de z
 
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
 
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
 
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
 
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
Línea 118: Línea 128:
 
   ylim([-1.2 1.2]);
 
   ylim([-1.2 1.2]);
 
   zlim([-z_max 0]);
 
   zlim([-z_max 0]);
−
  % Pausa para actualizar la animación
+
 
 
   pause(0.1);
 
   pause(0.1);
−
  %Reinicio grafica
+
 
 
   cla
 
   cla
 
   plot3(u_m,v_m,-z_vals,'Color','g')
 
   plot3(u_m,v_m,-z_vals,'Color','g')
Línea 129: Línea 139:
 
end
 
end
 
hold off
 
hold off
−
%%
+
 
 +
 
 
figure(2);
 
figure(2);
 
view(3)
 
view(3)
Línea 139: Línea 150:
 
   u = u_m(k);  % Componente u(z)
 
   u = u_m(k);  % Componente u(z)
 
   v = v_m(k);  % Componente v(z)
 
   v = v_m(k);  % Componente v(z)
−
   % Calcular los campos vectoriales para este valor de z
+
 
 +
   % Calculo campos vectoriales para este valor de z
 
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
 
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
 
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
 
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
Línea 145: Línea 157:
 
   ylim([-1.2 1.2]);
 
   ylim([-1.2 1.2]);
 
   zlim([-z_max 0]);
 
   zlim([-z_max 0]);
−
  % Pausa para actualizar la animación
+
 
 
   pause(0.1);
 
   pause(0.1);
−
  %Reinicio grafica
+
 
   cla
 
   cla
 
   plot3(u_m,v_m,-z_vals,'Color','g')
 
   plot3(u_m,v_m,-z_vals,'Color','g')

Revisión del 16:18 3 dic 2025

Trabajo realizado por estudiantes
Título La espiral de Ekman. Grupo 16.
Asignatura Teoría de Campos
Curso 2025-26
Autores Angela María Clarambo Nassarre
Alejandra García-Agullo Canle
María Cabranes Gonzalez
Alvaro Román Aguilera
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


1 Introducción

2 Parametros de coriolis f y profundidad de Ekman dE

El parámetro de Coriolis cuantifica la influencia local de la rotación terrestre en el movimiento de fluidos y objetos sobre la Tierra. Se define como:

[math] f= 2 Ω \sin ( \phi ) [/math]

Sustituyendo la velocidad angular Ω y la latitud expresada en radianes [math] \phi = \pi/6 [/math], nos da la siguiente expresión:

[math] f = 2 \cdot 7.2921 \cdot 10 ^ { - 5 } \cdot \sin ( \pi / 6 ) = 7.2921 \cdot 10 ^ { - 5 } [/math]

Para identificar en qué lugares de la Tierra 𝑓 es positivo, negativo o nulo, debemos fijarnos en [math] \phi [/math], es decir, en la latitud.

  • 𝑓 será nulo en el Ecuador (latitud 0°). No hay efecto de Coriolis horizontal.
  • 𝑓 será positivo en todo el Hemisferio Norte. La fuerza de Coriolis desvía los movimientos hacia la derecha.
  • 𝑓 será negativo en todo el Hemisferio Sur. La fuerza de Coriolis desvía los movimientos hacia la izquierda.

La profundidad de Ekman 𝑑𝐸, define la capa de agua directamente influenciada por la fricción del viento en la superficie y desviada por la fuerza de Coriolis.

Para calcularla, es necesario saber |𝑓|, la viscosidad turbulenta y aplicar la siguiente fórmula:

[math] d _ { E} = \sqrt { \frac { 2 v_ { e } } { | f | }} = 37.0317 [/math] metros

3 Determinar el valor de ϑ

θ es una fase inicial que depende de la dirección del viento, que a su vez está relacionada con la fuerza de coriolis. En la superficie, la dirección del flujo de agua se desvía θ con respecto la dirección del viento.

Como hemos visto antes nuestro valor de f es positivo por lo cual nos encontramos en el hemisferio norte y por ende [math] sgn(f) =1 [/math]

[math]\color{white} ‘f= 2 Ω \sin ( \phi ) [/math]

[math] \color{white} “f = 2 \cdot 7.2921 \cdot 10 ^ { - 5 } \cdot \sin ( \pi / 6 ) = 7.2921 \cdot 10 ^ { - 5 } \gt 0 [/math]

[math] u ( z ) = V _ { 0 } e ^ { \frac { z } {d _ { E }} } \cos ( \frac { z } { d _ { E } } + ϑ)[/math]

[math]v ( z ) = V _ { 0 } e ^ { \frac { z } {d _ { E }} } s e n ( \frac { z } { d _ { E } } + ϑ )[/math]

[math] \rightarrow z = 0 \rightarrow [/math][math] u( 0 ) = v _ { 0 } . \cos ( ϑ )[/math][math] \color{white} "v( 0 ) = v _ { 0 } . \sin ( ϑ )[/math]

Como el viento va de norte a sur, la fuerza de coriolis genera una corriente hacia el este entonces el flujo en la superficie se desviará π/4 rad en esa dirección .

Por lo tanto, [math] ϑ = \frac { 3 \pi} { 4 } [/math]

4 Soluciones de las ecuaciones diferenciales de Ekman

Las soluciones propuestas para las componentes horizontales de la velocidad en la capa de Ekman son:

[math] u ( z ) =V_ { 0 } e ^ { \frac { z } {d _ { E }} } . \cos ( \frac { z } { d _ { E } } + ϑ ) [/math] [math] v ( z ) = V_ { 0 } e ^ { \frac { z } {d _ { E }} } s e n ( \frac { z } { d _ { E } } + ϑ) [/math]

Para verificar que estas expresiones son solución del sistema, deben satisfacer las ecuaciones diferenciales de Ekman:

[math] \frac { d ^ { 2 } u } { d z ^ { 2 } } = - \frac { f } { v _ { e } } v [/math] [math] \frac { d ^ { 2 } v } { d z ^ { 2 } } = \frac { f } { v _ { e } } u [/math]

Cálculo de primeras derivadas:

[math] \frac { d u } { d z } = V _ { 0 } \frac { 1 } { d_{E} } e ^ { \frac { z } { d_{E} } }( \cos ( \frac { z} { d_{E}} + ϑ ) - \sin ( \frac { z } { d_{E} } + ϑ ) ) [/math]


[math] \frac { d v} { d z } = V _ { 0 } \frac { 1 } { d_{E} } e ^ { \frac { z } { d_{E} } }( \sin ( \frac { z } { d_{E} } + ϑ ) + \cos ( \frac { z} { d_{E}} +ϑ )) [/math]

Cálculo de segundas derivadas:

[math] \frac { d ^ { 2 } u } { d z ^ { 2 } } = \frac { -2 } { d {E}^ { 2 } } V _ { 0 } e ^ { \frac { z } { d{E} } } \sin ( \frac { z } { d _{E}} + ϑ) [/math]


[math] \frac { d ^ { 2 } v } { d z ^ { 2 } } = v _ { 0} \frac { 1 } { d_{E} }( \frac { 1 } { d_{E} } e ^ { \frac { z} { d_{E}} } ( 2 \cos ( \frac { z} { d_{E}} + ϑ ) )) = V_ { 0} \frac { 2 } { d_{E}^2 }e ^ { \frac { z} { d_{E}} } ( \cos ( \frac { z} { d_{E}} + ϑ )) [/math]

Igualamos cada derivada a su derivada de Ekman:

[math] \frac { d ^ { 2 } v } { d z ^ { 2 } } = v _ { 0} \frac { 1 } { d_{E} }( \frac { 1 } { d_{E} } e ^ { \frac { z} { d_{E}} } ( 2 \cos ( \frac { z} { d_{E}} + ϑ) ) = V_ { 0} \frac { 2 } { d_{E}^2 }e ^ { \frac { z} { d_{E}} } ( \cos ( \frac { z} { d_{E}} + ϑ ) )[/math]

Sustituiyendo [math] d_{E} [/math] en la ecuación, nos queda:

[math] \frac { d ^ { 2 } u } { d z ^ { 2 } } = - \frac { f } { v e } v \rightarrow \frac { f } { v e } = \frac { 2 } { d_ { E } ^ { 2 } } [/math]

Ahora se sustituye [math] d_{E} [/math] y así se confirma que se verifica [math]f = |f |[/math]

5 Campo vectorial V en varios planos horizontales

Apartado 5

clc;close all;clear;
omega = 7.2921e-5;  % Velocidad angular (rad/s)
phi = 30.10242;  % Latitud (grados)
V0 = 0.15;  % Velocidad superficial inducida por el viento (m/s)
nu = 0.05;  % Viscosidad turbulenta (m^2/s)
theta0 = 3*pi/4;  % Ángulo inicial (radianes)

% Cálculo del parámetro de Coriolis f
f = 2 * omega * sin(phi);

% Profundidad de Ekman
delta = sqrt(2 * nu / abs(f));

z_max = 3 * delta;  % Profundidad máxima para la animación
n_frames = 40;  % Número de frames de la animación
z_vals = linspace(0, z_max, n_frames);  % Valores de profundidad
[X, Y] = meshgrid(linspace(-1, 1, 10), linspace(-1, 1, 10));  

figure(1);
view(3)
%axis equal;

% Calculo velocidad
u_m = (f/abs(f))*V0 * exp(-z_vals / delta)  .* cos(-z_vals / delta+theta0);  % Componente u(z)
v_m = V0 * exp(-z_vals / delta)  .* sin(-z_vals / delta+theta0);  

% Componente v(z)
plot3(u_m,v_m,-z_vals,'Color','g')
hold on

% Animación
for k = 1:n_frames
   z = z_vals(k);  % Profundidad actual
   u = u_m(k);  % Componente u(z)
   v = v_m(k);  % Componente v(z)

   % Calculo de campos vectoriales para este valor de z
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
   xlim([-1.2 1.2]);
   ylim([-1.2 1.2]);
   zlim([-z_max 0]);

   pause(0.1);
  
   cla
   plot3(u_m,v_m,-z_vals,'Color','g')
   xlabel('Oeste - Este (m)');
   ylabel('Norte - Sur (m)');
   zlabel('Profundidad (m)');
   grid on
end
hold off


figure(2);
view(3)
plot3(u_m,v_m,-z_vals,'Color','g')
hold on
% Inicializar la animación
for k = 1:n_frames
   z = z_vals(k);  % Profundidad actual
   u = u_m(k);  % Componente u(z)
   v = v_m(k);  % Componente v(z)

   % Calculo campos vectoriales para este valor de z
   quiver3(X, Y,-z*ones(size(X)), u * ones(size(X)), v * ones(size(Y)),zeros(size(X)), 'AutoScale', 'off', 'LineWidth', 1.5,'Color','b');
   title(sprintf('Campo Vectorial de Ekman en z/d_E = %.2f', -z / delta));
   xlim([-1.2 1.2]);
   ylim([-1.2 1.2]);
   zlim([-z_max 0]);

   pause(0.1);
 
   cla
   plot3(u_m,v_m,-z_vals,'Color','g')
   xlabel('Oeste - Este (m)');
   ylabel('Norte - Sur (m)');
   grid on
   view([0 90])
end
hold off


6 Representación de los vectores V evaluados en puntos del eje vertical

7 Divergencia de V

El campo está definido por: [math]\overrightarrow {V} = u(z) \overrightarrow { i } + v(z) \overrightarrow { j }[/math]

Donde las componentes para el modelo clásico de Ekman son:

[math] u ( z ) = s g n ( f ) \cdot v _ { 0 } e ^ { z / d _ { E } } . \cos ( \frac { z } { d _ { E } } + ϑ)[/math]

[math]v ( z ) = v _ { 0 } e ^ { \frac { z } {d _ { E }} } s e n ( \frac { z } { d _ { E } } + ϑ )[/math]

Un aspecto fundamental del flujo de Ekman es que las componentes de velocidad u y v dependen exclusivamente de la coordenada vertical z, no varían en las direcciones horizontales.

La divergencia horizontal del campo vectorial se define como: [math]\nabla\overline { V } = \frac { d u } { d x } + \frac { d v } { d y } = 0[/math]

Esto significa que el flujo de Ekman es incompresible a escala horizontal. No existe convergencia ni divergencia neta de masa de agua en el plano horizontal dentro de la capa de Ekman. Esto implica que el transporte de masa asociado a la espiral se produce sin que se acumule o se vacíe agua en ningún punto, reflejando la continuidad del fluido en el modelo ideal.

8 Rotacional de V

El rotacional en la espiral de Ekman describe el giro del flujo dentro de la capa de Ekman, resultado de la combinación entre la fricción y la fuerza de Coriolis. El rotacional cuantifica la circulación del flujo por unidad de área y es clave para entender cómo el transporte en esta capa afecta procesos de mayor escala, como la formación de vórtices o el desarrollo de corrientes oceánicas.

En la espiral de Ekman, el flujo horizontal tiene componentes en las direcciones [math] x [/math] e [math] y [/math]:

[math]\overrightarrow { U} = u(z) \overrightarrow { i } + v(z) \overrightarrow { j }[/math],

El transporte neto de masa dentro de la capa de Ekman queda orientado de forma perpendicular al viento en la superficie, como consecuencia del equilibrio entre la fuerza de Coriolis y la fricción.

El rotacional de un campo de velocidades tridimensional se define como:

[math]\nabla×\vec u(x,y) = {\begin{vmatrix} \vec i & \vec j & \vec k \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ \vec u_x & \vec u_y & \vec u_z \end{vmatrix}} [/math]


En el caso de la espiral de Ekman, donde asumimos un flujo horizontal [math](\vec u_z = 0)[/math] y homogéneo en las direcciones horizontales. Esto significa que:

1. La variación del componente [math] v [/math] con la profundidad genera una rotación en la dirección x.

2. La variación del componente [math] u [/math] con la profundidad genera una rotación en la dirección y.


El rotacional global será tridimensional, pero su componente dominante estará vinculado al gradiente de velocidad en la profundidad.


A continuación, se adjunta un programa de Matlab que representa el rotacional en varios planos paralelos, desde la superficie del mar hasta la profundidad de Ekman:

9 Reflexión del transporte de Ekman (por unidad de ancho)

10 Calculo flujo neto de V a través de una pared vertical

11 La espiral de Ekman en coordenadas cilíndricas

12 Curvatura y Torsión de la esprial de Ekman

La curvatura y la torsión ayudan a describir cómo se comporta la espiral de Ekman. La curvatura indica cuánto cambia la dirección de la trayectoria en cada punto, mientras que la torsión mide cómo va girando el plano en el que se curva la trayectoria, reflejando así su carácter tridimensional. Las expresiones para la curvatura y la torsión se definen como: :[math] \kappa (z) y \tau (z)[/math]

13 Calculo y representación del tiedro de Frenet

14 Calculo de la longitud del arco de la espiral de Ekman

15 La proyección de la espiral de Ekman sobre el plano horizontal

15.1 Representación logaritmica en el plano XY

15.2 Verificación de p=p0ebtheta

15.3 Calculo del ángulo entre vector posición y vector tangente. Verificación de que es constante

15.4 Aplicaciones de la espiral logaritmica en ingeniería

La espiral logarítmica (o espiral equiangular) es una de las curvas más estudiadas en matemáticas aplicadas debido a sus propiedades geométricas únicas y su frecuente aparición en sistemas naturales y artificiales.

Aplicaciones en ingeniería:

Diseño aeronáutico y de turbomáquinas:

Álabes de turbinas y compresores: El perfil de los álabes en turbinas eólicas, turbinas de gas y compresores centrífugos a menudo sigue una espiral logarítmica. Esto se debe a que mantiene un ángulo de ataque constante a lo largo del radio, lo que minimiza las pérdidas por choque y mejora la eficiencia aerodinámica o hidrodinámica.

Ductos de difusores y toberas: La forma espiral logarítmica permite una expansión o compresión gradual del flujo, reduciendo la formación de turbulencias y la separación de la capa límite.


Antenas y electromagnetismo:

Antenas espirales logarítmicas: Son ampliamente utilizadas en telecomunicaciones y radioastronomía por su característica de banda ancha. Pueden operar en un amplio rango de frecuencias sin cambios significativos en su impedancia y patrón de radiación.

Diseño de guías de onda y acopladores direccionales: La curva permite una transición suave de la energía electromagnética, minimizando reflexiones.


Ingeniería mecánica y diseño de resortes:

Resortes espirales no cilíndricos: En aplicaciones donde se requiere una relación no lineal entre fuerza y desplazamiento, los resortes con forma de espiral logarítmica ofrecen una respuesta mecánica más controlada.

Mecanismos de seguimiento solar: Algunos diseños de concentradores solares utilizan superficies reflectantes en forma de espiral logarítmica para enfocar la radiación de manera uniforme a lo largo del día.


Ingeniería civil y arquitectura:

Diseño estructural de escaleras y rampas: La espiral logarítmica proporciona una transición ergonómica y estéticamente agradable entre niveles.

Distribución de tensiones en cúpulas y arcos: Algunas estructuras de cubierta se inspiran en esta curva para redistribuir cargas de manera óptima.

Procesamiento de señales e imagen:

Algoritmos de compresión y reconocimiento de patrones: La espiral logarítmica sirve como base para transformadas matemáticas (como la transformada espiral) utilizadas en el análisis de texturas y formas naturales.

Obtenido de «https://mat.caminos.upm.es/w/index.php?title=La_espiral_de_Ekman_(Grupo_16)&oldid=96343»