Flujo alrededor de un obstáculo circular (Grupo 29)
De MateWiki
Revisión del 22:29 24 nov 2025 de Rodrigo de Pelayo (Discusión | contribuciones)
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Flujo alrededor de un obstáculo circular (Grupo 29) |
| Asignatura | Teoría de Campos |
| Curso | 2025-26 |
| Autores | Rodrigo de Pelayo García García
Sergio Resino Velayo Cayetano Gilabert Castejón Hugo Moreno Peregrina |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
El estudio del flujo alrededor de un cilindro circular es uno de los problemas más clásicos y analizados en la mecánica de fluidos teórica y aplicada. Su importancia en ingeniería civil es notable, pues muchas estructuras presentan formas aproximadamente cilíndricas: pilas de puentes, torres, tirantes, pilotes o tuberías expuestas a corrientes o viento.
Comprender el comportamiento del flujo permite predecir fenómenos como:
-Distribución de presiones alrededor del obstáculo
-Fuerzas de arrastre y sustentación,
-Desprendimiento periódico de vórtices
-Fenómenos de resonancia estructural
-Erosión local en estructuras hidráulicas
Contenido
- 1 Región del fluido y mallado en coordenadas polares
- 2 Función potencial y campo de velocidades
- 3 Divergencia y Rotacional
- 4 Lineas de Corriente
- 5 Puntos de Remanso
- 6 Presión según Bernouilli
- 7 Trayectorias de Particulas
- 8 Circulaón y Paradoja D' Alembert
- 9 Flujo con circulación impuesta
- 10 Véase también
- 11 Referencias
1 Región del fluido y mallado en coordenadas polares
% Parámetros del anillo
rho_inner = 1; % radio interior (círculo unidad)
rho_outer = 5; % radio exterior
n_rho = 40; % número de curvas rho = cte (radiales concéntricas)
n_theta = 40; % número de curvas theta = cte (radiales)
% Generar malla en coordenadas polares
rho = linspace(rho_inner, rho_outer, n_rho);
theta = linspace(0, 2*pi, n_theta);
[U,V] = meshgrid(rho,theta);
%Prarametrizar la superficie
hold on
X=U.*cos(V);
Y=U.*sin(V);
Z=0.*U;
%Representación del mallado y dibujo del círculo
mesh(X,Y,Z);
plot(1*cos(theta), 1*sin(theta), 'LineWidth',2);
%Delimitar los ejer
axis([-4,4,-4,4]);
% Marcas de eje
box on
axis equal
grid off
% Etiquetas y título
xlabel('x')
ylabel('y')
title('Mallado del anillo')
hold off