Diferencia entre revisiones de «Coordenadas cilíndricas parabólicas Grupo 6B»

Revisión del 13:46 24 nov 2024

Trabajo F: Coordenadas Cilíndricas Parabólicas

1 Introducción

En este trabajo vamos a estudiar y aplicar las denominadas coordenadas cilíndricas parabólicas, que se denotan por (u, v, z). Estas tienen la siguiente relación con las coordenadas cartesianas (x₁, x₂, x₃):

[math] \begin{aligned} x_1 &= \frac{u^2 - v^2}{2}, \\ x_2 &= uv, \\ x_3 &= z, \end{aligned} [/math]

donde u > 0.

Las coordenadas cilíndricas parabólicas son una generalización de las coordenadas cilíndricas estándar y extienden un cambio de coordenadas en R² a todo el espacio R³. A continuación, se presentan los cálculos, representaciones y aplicaciones.

---

2 Parametrizaciones de las líneas coordenadas \(\gamma_u\), \(\gamma_v\) y \(\gamma_z\)

Líneas coordenadas en cartesianas:

• \(\gamma_u\): [math](x_1, x_2, x_3) = \left( \frac{t^2 - v^2}{2}, tv, z \right)[/math], con t variable y v, z constantes.
• \(\gamma_v\): [math](x_1, x_2, x_3) = \left( \frac{u^2 - t^2}{2}, ut, z \right)[/math], con t variable y u, z constantes.
• \(\gamma_z\): [math](x_1, x_2, x_3) = \left( \frac{u^2 - v^2}{2}, uv, t \right)[/math], con t variable y u, v constantes.

Figura 1: Líneas coordendas.

clear,clc
%Parametrizaciones de las lineas coordenadas
u = linspace(0.1, 2, 50); % Valores de u
v = linspace(0.1, 2, 50); % Valores de v
%Dibujo de las lineas coordenadas
figure;
hold on;

% Curvas gamma_u (variando u, con v fijo)
v_fixed = 1;
x1_u = (u.^2 - v_fixed^2) / 2;
x2_u = u .* v_fixed;
plot(x1_u, x2_u, 'r', 'LineWidth', 1.5);

% Curvas gamma_v (variando v, con u fijo)
u_fixed = 1;
x1_v = (u_fixed^2 - v.^2) / 2;
x2_v = (u_fixed) .* v;
plot(x1_v, x2_v, 'b', 'LineWidth', 1.5);

% Estilo del gráfico
title('Líneas coordenadas en el plano z = 0');
xlabel('x_1');
ylabel('x_2');
legend({'Líneas gamma_u (u varía)', 'Líneas gamma_v (v varía)'});
grid on;
axis equal;
hold off;

3 Velocidades de las líneas coordenadas \(\gamma_u\), \(\gamma_v\) y \(\gamma_z\)

Cálculos: Los campos velocidad de las líneas coordenadas son:

• Para γₐ:

[math]\gamma'_u = \left( u, v, 0 \right)[/math].

• Para γᵥ:

[math]\gamma'_v = \left( -v, u, 0 \right)[/math].

• Para γ_z:

[math]\gamma'_z = \left( 0, 0, 1 \right)[/math].

Factores de escala: Los factores de escala h_u, hᵥ, h_z son los módulos de los campos velocidad:

[math] h_u = |\gamma_u'(u)| = \sqrt{u^2 + v^2}, \quad h_v = |\gamma_v'(v)| = \sqrt{u^2 + v^2}, \quad h_z = |\gamma_z'(z)| = 1. [/math]

Vectores tangentes: Los vectores tangentes unitarios son:

• [math]\mathbf{e}_u = \frac{\gamma'_u}{h_u} = \frac{1}{\sqrt{u^2 + v^2}} \left( u, v, 0 \right)[/math],
• [math]\mathbf{e}_v = \frac{\gamma'_v}{h_v} = \frac{1}{\sqrt{u^2 + v^2}} \left( -v, u, 0 \right)[/math],
• [math]\mathbf{e}_z = \left( 0, 0, 1 \right)[/math].

Figura 1: Líneas coordendas.

clear,clc,clf
% Punto de interés
u = 1;
v = 1;

% Vectores unitarios en ese punto
h = sqrt(u^2 + v^2);
eu = [u/h, v/h, 0]; % Vector e_u
ev = [-v/h, u/h, 0]; % Vector e_v

% Coordenadas cartesianas de las líneas coordenadas en el plano z = 0
x1_u = (u.^2 - v.^2) / 2;
x2_u = u .* v;

% Gráfico
figure;
hold on;
quiver(x1_u, x2_u, eu(1), eu(2), 'r', 'LineWidth', 1.5);
quiver(x1_u, x2_u, ev(1), ev(2), 'b', 'LineWidth', 1.5);
plot(x1_u, x2_u, 'k--', 'LineWidth', 1); % Línea coordenada
title('Vectores unitarios en el plano z = 0');
xlabel('x_1');
ylabel('x_2');
legend({'e_u', 'e_v'});
grid on;
axis equal;
hold off;

---

4 2. Campos velocidad y factores de escala

Cálculos: Los campos velocidad de las líneas coordenadas son:

• Para γₐ:

[math]\gamma'_u = \left( u, v, 0 \right)[/math].

• Para γᵥ:

[math]\gamma'_v = \left( -v, u, 0 \right)[/math].

• Para γ_z:

[math]\gamma'_z = \left( 0, 0, 1 \right)[/math].

Factores de escala: Los factores de escala hₐ, hᵥ, h_z son los módulos de los campos velocidad: [math] h_u = \sqrt{u^2 + v^2}, \quad h_v = \sqrt{u^2 + v^2}, \quad h_z = 1. [/math]

Vectores tangentes: Los vectores tangentes unitarios son:

• [math]\mathbf{e}_u = \frac{\gamma'_u}{h_u} = \frac{1}{\sqrt{u^2 + v^2}} \left( u, v, 0 \right)[/math],
• [math]\mathbf{e}_v = \frac{\gamma'_v}{h_v} = \frac{1}{\sqrt{u^2 + v^2}} \left( -v, u, 0 \right)[/math],
• [math]\mathbf{e}_z = \left( 0, 0, 1 \right)[/math].

Visualización: Dibujar con MATLAB las líneas coordenadas γₐ y γᵥ junto con los vectores tangentes en un punto arbitrario.

---

5 3. Matrices de cambio de base

Matriz de cambio de base: [math] Q = \begin{bmatrix} \mathbf{e}_u & \mathbf{e}_v & \mathbf{e}_z \end{bmatrix}. [/math]

Matriz inversa: [math] Q^{-1} = Q^\top, \quad \text{ya que la base es ortonormal}. [/math]

---

6 4. Campo posición en coordenadas cilíndricas parabólicas

El campo posición se expresa como: [math] \mathbf{r}(u, v, z) = \frac{u^2 - v^2}{2} \mathbf{i} + uv \mathbf{j} + z \mathbf{k}. [/math]

---

7 5. Gradiente de un campo escalar

El gradiente de un campo escalar en coordenadas cilíndricas parabólicas es: [math] \nabla f = \frac{\partial f}{\partial u} \frac{\mathbf{e}_u}{h_u} + \frac{\partial f}{\partial v} \frac{\mathbf{e}_v}{h_v} + \frac{\partial f}{\partial z} \frac{\mathbf{e}_z}{h_z}. [/math]

Aplicación: calcular el gradiente del campo f(x₁, x₂, x₃) = x₂. 1. Convertir f a coordenadas cilíndricas parabólicas: [math]f(u, v, z) = uv[/math]. 2. Evaluar el gradiente en el punto cartesiano (x₁, x₂, x₃) = (0, 1, 1).

---

8 6. Divergencia

La divergencia en este sistema es: [math] \nabla \cdot \mathbf{F} = \frac{1}{h_u h_v h_z} \left[ \frac{\partial}{\partial u}(h_v h_z F_u) + \frac{\partial}{\partial v}(h_u h_z F_v) + \frac{\partial}{\partial z}(h_u h_v F_z) \right]. [/math]

Ejemplo: calcular la divergencia del campo posición.

---

9 7. Rotacional

La expresión del rotacional en este sistema es: [math] \nabla \times \mathbf{F} = \frac{1}{h_u h_v h_z} \begin{vmatrix} \mathbf{e}_u & \mathbf{e}_v & \mathbf{e}_z \\ \frac{\partial}{\partial u} & \frac{\partial}{\partial v} & \frac{\partial}{\partial z} \\ h_u F_u & h_v F_v & h_z F_z \end{vmatrix}. [/math]

---

10 8. Superficies de nivel

Las superficies de nivel para los campos escalares son:

• f₁(u, v, z) = u: Superficie parabólica.
• f₂(u, v, z) = v: Superficie parabólica.
• f₃(u, v, z) = z: Plano horizontal.

Visualización: Dibujar cada superficie de nivel en MATLAB y analizar si son superficies regladas.

---

11 9. Curvatura de una parábola

La parábola es: [math] y = -2x^2 + 2. [/math]

Curvatura: La curvatura es: [math] \kappa(x) = \frac{|y''(x)|}{(1 + (y'(x))^2)^{3/2}}. [/math] Evaluar y graficar κ(x) en MATLAB para x ∈ [-1, 1].

---

12 10. Uso de la parábola en ingeniería

La parábola tiene múltiples aplicaciones en ingeniería, como:

• Diseño de puentes (arcos parabólicos).
• Antenas parabólicas (reflectores).
• Elementos arquitectónicos.

Añadir imágenes de ejemplos y explicar su uso.