Fluido alrededor de un obstáculo circular (15-B)

Trabajo realizado por estudiantes
Título	Fluido alrededor de un obstáculo circular. Grupo 15-B
Asignatura	Teoría de Campos
Curso	2022-23
Autores	Julian Lavandero Pacheco Javier Sesmero Zamarrón Ana Sarró Redel Sergio Navarro Czyz
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

2 Mallado

Para realizar el mallado vamos a utilizar Matlab. Vamos a dibujar la gráfica en los ejes [math](x,y) ∈ [-5,5]\times[-5,5] [/math].

Mallado del anillo

rho=linspace(1,5,30);                        %definimos rho y theta 
th=linspace(0,2*pi,30);

[U,V]=meshgrid(rho,th);                      %Creamos la malla

hold on
X=U.*cos(V);                                 %Parametrizamos la superficie
Y=U.*sin(V);
Z=0.*U;
mesh(X,Y,0*Z);                               %Dibujamos la placa 
plot(1*cos(th),1*sin(th),'k','lineWidth',1); %Representación del obstáculo
axis([-5,5,-5,5]);                           %Establecemos los ejes 
view(2);
title ('Placa');
xlabel 'EJE X'
ylabel 'EJE Y'
axis equal
hold off

3 Función potencial y campo de velocidades

Función potencial

rho=linspace(1,5,30);                             %Definimos rho y theta
th=linspace(0,2*pi,30);

[U,V]=meshgrid(rho,th);                           %Creamos la malla

hold on
X=U.*cos(V);                                      %Parametrización de la superficie
Y=U.*sin(V);

f=@(rho,th)(rho+(1./rho)).*cos(th) +sqrt(2).*th;  %Definimos la función potencial 
Z=f(U,V);                                         %Aplicamos la función potencial
surf(X,Y,Z);                                      %Dibujamos la función
plot(1*cos(th),1*sin(th),'k','lineWidth',3);      %Representamos el obstáculo

view(2);  
axis([-5,5,-5,5]);
title ('Función potencial');
xlabel ('EJE X');
ylabel ('EJE Y');
axis equal
hold off