Diferencia entre revisiones de «La espiral de Ekman (grupo 5)»

Revisión del 12:34 3 dic 2024

En este artículo se analizará la influencia de la espiral de Ekman sobre las corrientes oceánicas. Cuando el viento sopla, se crea una corriente superficial sobre el mar en la dirección del viento, sin embargo, debido a la fuerza de Coriolis, esta dirección no permanece constante. La dirección de la velocidad del flujo se desvía gradualmente formando la espiral de Ekman. En consecuencia, el transporte total del agua seguirá un desplazamiento perpendicular a la dirección del viento en la superficie.

La velocidad del fluido en cada capa de profundidad se describe con [math] \vec{v}= u\vec{i}+v\vec{j} [/math]

Siendo:

• [math] V_{0} [/math] es la intensidad de la velocidad superficial inducida por el viento;
  
• [math] d_{E}=\sqrt{\frac{2\times V_{e}}{|f|}} [/math] es la profundidad de Ekman, la distancia vertical en la que la influencia del viento y la fuerza de Coriolis afectan significativamente el movimiento del agua
  
• [math] \vartheta [/math] es una fase inicial determinada por la dirección del viento respecto a la fuerza de Coriolis;
  
• sgn es la función signo:
  

1 Parámetro de Coriolis

El parámetro de Coriolis está definido por la siguiente fórmula: [math] f=2\Omega sen(\phi ) [/math], donde [math] \Omega [/math] es la velocidad angular de la Tierra y [math] \phi [/math] la latitud. Aproximadamente, la velocidad angular de la tierra es [math] \Omega =7.2921\cdot 10^{-5} rad/s [/math].

Para una latitud de 45°N, es decir, de [math] \phi =\frac{\sqrt{2}}{2}rad [/math], se determina sustituyendo en la ecuación anterior que el parámetro de Coriolis es de:

Dependiendo de la latitud, el valor del parámetro de Coriolis puede ser negativo o positivo. Sustituyendo en la fórmula dada, se demuestra que:

• Para el hemisferio norte [math] (0^{\circ}\lt \phi \lt 90^{\circ}) [/math], el parámetro [math] f\gt0 [/math]
• Para el Ecuador [math] (\phi = 0^{\circ}) [/math], el parámetro [math] f=0 [/math]
• Para el hemisferio sur [math] (-90^{\circ}\lt \phi \lt 0^{\circ}) [/math], el parámetro [math] f\lt0 [/math]

2 Valor de fase inicial θ

En el apartado anterior se ha obtenido que [math] f=1.0312\cdot 10^{-4} \frac{1}{s} [/math]. Para obtener la fase inicial, utilizamos las dos siguientes ecuaciones:

• [math]u(z)=sgn(f)V_{o}e^{z/d_{E}}cos\left ( \frac{z}{d_{E}}+\vartheta \right )[/math]
  
• [math]v(z)=V_{o}e^{z/d_{E}}sin\left ( \frac{z}{d_{E}}+\vartheta \right ) [/math]
  

Al ser [math] f\gt0 [/math], concluimos que [math] sgn(f)=1 [/math], y por lo tanto [math]u(z)=V_{o}e^{z/d_{E}}cos\left ( \frac{z}{d_{E}}+\vartheta \right )[/math].
Para [math] z=0 [/math]:

• [math]u(0)=V_{o}cos\left ( \vartheta \right )[/math]
  
• [math]v(0)=V_{o}sin\left ( \vartheta \right ) [/math]
  

La ángulo inicial del flujo de Ekman está desviado en superficie aproximadamente [math] 45^{\circ} [/math], es decir, [math] \pi /4 [/math] rad.
Resolviendo las ecuaciones anteriores, se obtiene que la fase inicial es [math] \vartheta= -\frac{\pi }{4}=\frac{3\pi }{4}[/math]

3 Soluciones de las ecuaciones diferenciales de Ekman

Las ecuaciones diferenciales de Ekman vienen definidas de la siguiente manera:

Cuyas soluciones se muestran:

Por lo tanto, para verificar la igualdad, tan solo hay que derivar dos veces respecto de z las soluciones propuestas:

4 Animación de campo vectorial \( \vec{v}\)

Animación del campo vectorial V

%Parametros
V0=0.2;                     
visc=0.1;                   
phi=pi/4;                   
omega=7.2921e-5;            
f=2*omega*sin(phi);         
dE=sqrt(2*visc/abs(f));     
theta=-pi/4;                

% Profundidad de ekman
zz=linspace(0,-dE,30);  % Desde la superficie (z=0) hasta -d_E

% Creación de la figura para la animación
figure;
hold on;
axis equal;
xlabel('U');
ylabel('V');
title('Campo vectorial de la espiral de Ekman');
for i=1:length(zz)
    z=zz(i);

    u=sign(f)*V0*exp(z/dE)*cos((z/dE)+theta);
    v=V0*exp(z/dE)*sin((z/dE)+theta);
   
    quiver(0,0,u,v);
    pause(0.4);
    
    hold on;
    axis equal;
    xlabel('Este (x)');
    ylabel('Norte (y)');
    title('Campo vectorial de la espiral de Ekman');
    grid on;
end

5 Representación de campo vectorial \( \vec{v}\)

Representación del campo vectorial

% Parámetros iniciales
f=1.03e-4;          
V0=0.2;             
visc=0.1;           
dE=sqrt(2 * visc / abs(f));  
theta=-pi/4;         

% Rango de profundidades
z_max =5*dE;      
z_vals=linspace(0,-z_max,30); 

% Crear la figura
figure;
hold on;
axis equal;
xlabel('Este (m)');
ylabel('Norte (m)');

% Recorrer los valores de z
for i = 1:length(z_vals)
    z = z_vals(i);
   
    u = sign(f) * V0 * exp(z/dE) * cos(z/dE + theta);
    v = V0 * exp(z/dE) * sin(z/dE + theta);

    quiver(0, 0, u, v);
end

6 Cálculo y representación de divergencia de \( \vec{v}\)

La divergencia del campo vectorial \( \vec{v}\) en coordenadas cartesianas se define como [math] \nabla\cdot \vec{v}(x,y,z)=\frac{\partial {u}}{\partial x}+\frac{\partial {v}}{\partial y}+\frac{\partial {w}}{\partial z} [/math]. Calculando, se determina que la divergencia de \( \vec{v}\) es:

Que la divergencia sea nula significa que no hay creación ni destrucción de flujo, es decir el volumen que entra es el mismo que el que sale, se trata de un fluido incompresible.

7 Cálculo y animación de rotacional de \( \vec{v}\)

El rotacional del campo vectorial \( \vec{v}\) en coordenadas cartesianas se define como [math] \nabla\times \vec{v}(x,y,z)=\begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ u_{x} & u_{y} & u_{z} \\ \end{vmatrix} [/math]. Sustituyendo, se halla el rotacional de \( \vec{v}\):

Animación del rotacional de V

%Parametros iniciales
V0=0.2;
theta=-pi/4;
dE=sqrt(2*0.1/abs(1.03e-4));

%Definicion de los rotacionales
z=linspace(0,-dE,50);
rot_x=((-V0*exp(z/dE))/dE).*(sin((z/dE)+theta)+cos((z/dE)+theta));
rot_y=((V0*exp(z/dE))/dE).*(-sin((z/dE)+theta)+cos((z/dE)+theta));
rot_z=zeros(size(z));


% Crear la animación
for k = 1:length(z)
    quiver3(0, 0, z(k), rot_x(k), rot_y(k), 0);  
    pause(0.4); 
    hold on;
end

8 Cálculo del flujo a través de una pared perpendicular

9 Parametrización de la espiral de Ekman en coordenadas cilíndricas

10 Cálculo de la curvatura y torsión de la espiral de Ekman

11 Calculo y animación del triedro de Frenet

12 Espiral logarítmica