Campos y temperatura en anillo plano (Grupo 14A)
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Campos y temperatura en anillo plano (Grupo 14A) |
| Asignatura | Teoría de Campos |
| Curso | 2014-15 |
| Autores | Jose Manuel Alonso de Caso, Raquel Álvarez,Alejandro del Rey, Raphael Gide,Miguel López, Alejandro Mendizabal. |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 . Mallado del anillo plano
Partimos de una superficie plana circular centrada en el origen, delimitada por los radios de valores 1 y 2 unidades, que expresaremos en coordenadas cilíndricas para el estudio en este trabajo.
%vectores de las coordenadas
rho=linspace(1,5,50);
phi=linspace(0,2*pi,50);
%malla
h=0.1; % Paso de muestreo
%coordenadas polares
r=1:h:2;
t=0:h:2*pi+h;
[rr,tt]=meshgrid(r,t); %mallado
%Parametrización
x=rr.*cos(tt);
y=rr.*sin(tt);
mesh(x,y,0*x)
axis([-3,3,-3,3])
2 . Distribución de temperaturas en el anillo
La temperatura se distribuye dentro del anillo siguiendo la función [math]T=(8-y^2+2y)exp(-x^2+2x-1)[/math]. La temperatura en cada punto del disco debido a esta distribución se puede intuir mediante las curvas de nivel de la representación gráfica adjunta abajo. Se trata de una distribución cuyo foco de calor se encuentra desplazado con respecto al origen de coordenadas. De esta forma se observan curvas concéntricas que están formadas por puntos que se encuentran a la misma temperatura (curvas de nivel).
Utilizaremos la expresión [math] max(max(T)) [/math] para obtener el valor máximo de la temperatura (T) para que matlab nos lo muestre. El máximo de temperatura en el dominio es [math] 8.9945 [/math], aproximadamente [math]9[/math]
h=0.1;
r=1:h:2;
t=0:h:2*pi+h;
[rr,tt]=meshgrid(r,t); %mallado
x=rr.*cos(tt);
y=rr.*sin(tt);
mesh(x,y,0*x);
axis([-3,3,-3,3])
T=(8-y.^2+2.*y).*exp(-x.^2+2.*x-1); % Expresión del campo escalar térmico
subplot(2,1,1)
surf(x,y,T)
axis([-3,3,-3,3])
subplot(2,1,2)
contour(x,y,T,25);
axis([-3,3,-3,3])
3 . Gradiente de la temperatura
El gradiente de la temperatura [math]∇T=(T_1,T_2)= \frac{\partial T}{\partial x} \vec i + \frac{\partial T}{\partial y} \vec j = (-2x+2).exp(-x^2+2x-1) \vec i + (-2*y+2).exp(-x^2+2*x-1) \vec j [/math] se define como ortongonal a las lineas de nivel, como podemos comprobar visualmente en la figura. Se observa también que el gradiente en cada punto converge hacia el foco de calor.
h=0.1;
r=1:h:2;
t=0:h:2*pi+h;
[r,t]=meshgrid(r,t);
x=r.*cos(tt);
y=r.*sin(tt);
%Relación polares-cartesianas
u=sqrt(x.^2+y.^2);
v=atan(y./x);
f=(8-y.^2+2.*y).*exp(-x.^2+2.*x-1);
fx=(-2.*x+2).*exp(-x.^2+2.*x-1); % Derivada parcial en x
fy=(-2.*y+2).*exp(-x.^2+2.*x-1); % Derivada parcial en y
hold on
quiver(x,y,fx,fy)
axis([-3,3,-3,3])
contour(x,y,f,25) %Curvas de nivel
axis([-3,3,-3,3])
hold off
4 . Desplazamiento provocado por [math]\vec u[/math]
A continuación aplicamos el desplazamiento provocado por el vector [math]\vec u = (1-ρ)^2 g_ρ [/math] en el mallado y lo representamos gráficamente. Cada vector representa la fuerza aplicada en cada punto del sólido y nos indica la dirección de desplazamiento que lleva cada punto.
h=0.1;
r=1:h:2;
t=0:h:2*pi+h;
[rr,tt]=meshgrid(r,t);
x=rr.*cos(tt);
y=rr.*sin(tt);
u=sqrt(x.^2+y.^2);
v=atan(y./x);
ux=((1-u).^2).*x./u;
uy=((1-u).^2).*y./u;
quiver(x,y,ux,uy)
axis([-3,3,-3,3])
5 . Sólido antes y despues del desplazamiento
%malla
h=0.1; %paso de muestreo
%coordenadas polares
r=1:h:2;
t=0:h:2*pi+h;
[rr,tt]=meshgrid(r,t); %mallado
x=rr.*cos(tt);
y=rr.*sin(tt);
u=sqrt(x.^2+y.^2);
v=atan(y./x);
%Desplazamientos originados por la vibración
ux=((1-u).^2).*x./u;
uy=((1-u).^2).*y./u;
%Posición final de los puntos de la placa
dx=x+ux;
dy=y+uy;
subplot(1,2,1) %antes
mesh(x,y,0.*x)
axis([-3,3,-3,3])
subplot(1,2,2) %despues
mesh(dx,dy,0.*dx)
axis([-3,3,-3,3])
