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

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.

1 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.

Mallado que representa los puntos de la región ocupada por un fluido

% 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);   % obstáculo circular

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;

2 Velocidad de las partículas del fluido

En este apartado analizaremos la velocidad de las partículas dada por el gradiente de la siguiente función potencial:

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

2.1 .-Representación de la Función Potencial

Primero representaremos la función potencial que describe el flujo asociado al movimiento de un fluido incompresible alrededor de un obstáculo circular.

2.2 .-Representación del campo de velocidades

Representaremos el campo de velocidades derivado de la función potencial, 𝑢 ⃗ =∇φ. 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.