Flujo alrededor de un obstáculo circular (Grupo 26)

Trabajo realizado por estudiantes
Título	Flujo alrededor de un obstáculo circular (Grupo 26)
Asignatura	Teoría de Campos
Curso	2025-26
Autores	Gonzalo Gallego Fulgencio Andrea García Carrasco Aaron García Martin Miryam Sanchez-Ferragut Samalea Guillermo Rodriguez Navadijos
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción

Vamos a estudiar el flujo de un fluido incompresible alrededor de un obstáculo circular, trabajando en el plano y utilizando coordenadas cilíndricas (polares) para describir el campo de velocidades y las condiciones en la superficie del cilindro. Este enfoque permite formular de manera directa las ecuaciones del flujo potencial y analizar cómo la presencia del obstáculo modifica la distribución de velocidades y presiones. A partir de este planteamiento se desarrollarán las cuestiones que se piden a continuación.

2 Dibujar mallado

En este primer apartado representaremos la región ocupada por el fluido, que corresponde al exterior del círculo unidad. Para ello construiremos un mallado en coordenadas polares que cubra el anillo comprendido entre los radios 1 y 5, con centro en el origen. Este mallado permitirá visualizar los puntos interiores de la zona de estudio y establecer la geometría sobre la que se formulará posteriormente el problema del flujo. Para completar la representación, dibujaremos también los ejes cartesianos en el dominio [ − 4 , 4 ] × [ − 4 , 4 ] [−4,4]×[−4,4], lo que facilitará interpretar la posición del obstáculo circular y la extensión del fluido respecto al sistema de referencia.

% Trabajo P - Apartado (1)
% Mallado del anillo 1 <= r <= 5 en coordenadas polares

clear; clc; close all;

R1  = 1;      % radio interior (obstáculo)
R2  = 5;      % radio exterior del fluido
Nr  = 25;     % número de divisiones radiales
Nth = 80;     % número de divisiones angulares

rho   = linspace(R1, R2, Nr);
theta = linspace(0, 2*pi, Nth);

[RHO, TH] = meshgrid(rho, theta);

X = RHO .* cos(TH);
Y = RHO .* sin(TH);
Z = 0.*RHO;
figure; hold on;

% Líneas radiales (theta = constante)
for i = 1:Nth
    plot(X(i,:), Y(i,:), 'g');
end

% Circunferencias (r = constante)
for j = 1:Nr
    plot(X(:,j), Y(:,j), 'g');
end

% Obstáculo circular (r = 1) representado solo con contorno
th_circ = linspace(0, 2*pi, 400);
x_circ  = R1 * cos(th_circ);
y_circ  = R1 * sin(th_circ);
plot(x_circ, y_circ, 'k', 'LineWidth', 2);   % círculo negro grueso

axis equal;
xlim([-4 4]);
ylim([-4 4]);
xlabel('x');
ylabel('y');
title('Mallado en el anillo 1 \leq r \leq 5 (flujo alrededor de un cilindro)');
grid off;
hold off;

3 Velocidad particular fluido

En este apartado analizaremos la función potencial:

[math] \varphi (\rho ,\theta)=(\rho +\frac{1}{\rho})\cos (\theta ) [/math]

que describe el flujo potencial asociado al movimiento de un fluido incompresible alrededor de un obstáculo circular. Dado que la velocidad de las partículas viene dada por el gradiente de esta función, 𝑢 ⃗ = ∇ 𝜑 u =∇φ, representaremos tanto la función potencial como el campo de velocidades derivado de ella. A partir de estas gráficas podremos visualizar la estructura del flujo y comprobar, mediante un zoom local, que el campo de velocidades es ortogonal a las curvas de nivel de 𝜑 φ, una propiedad característica de los flujos potenciales.