Observar el comportamiento de las componentes de potencial utilizando paquetes de software estudiantil

OBJETIVOS

• Observar el comportamiento de las componentes de potencial utilizando paquetes de software estudiantil.

ACTIVIDADES

• Analice cuidadosamente el script correspondiente a la ecuación de Laplace y verifique que este consiente en definir la geometría, las fronteras y el tipo de ecuación diferencial parcial a resolver; finalmente dibujar los resultados. Verificar si los resultados son los mismos a los expuestos en la guía.

EQUIPOS NECESARIOS

• Computador .

• Software FlexPDE.

• Matlab

• Pdetool

DESARROLLO

• Analice cuidadosamente el script correspondiente a la ecuación de Laplace y verifique que este consiente en definir la geometría, las fronteras y el tipo de ecuación diferencial parcial a resolver; finalmente dibujar los resultados. Verificar si los resultados son los mismos a los expuestos en la guía.

Punto 1

TITLE

'Temperatura'

VARIABLES T

EQUATIONS

div(grad(T)) = 0

BOUNDARIES

region 1

start (0,0)

natural(T) = 0 line to (0,4)

value(T) = 32 line to (2,2)

value(T) = 32 line to (4,4)

natural(T)=0 line to (4,0)

value(T)=212 line to close

PLOTS

grid(x,y)

contour(T) as 'Lineas de igual temperatura'

vector(-dx(Y),-dy(T)) as 'Flujo de Calor'

END

Flujo de calor

Líneas de igual temperatura

Disipador de calor

Código Flexpde:

TITLE 'Disipador De Calor'

VARIABLES T

DEFINITIONS

coe = 1

longitud=3

! INITIAL VALUES

EQUATIONS

div(coe*grad(T))=0

BOUNDARIES

REGION 1

START(0,7)

natural (T)=30 line to (7,7)

natural (T)=30 line to (7,0)

natural (T)=30 line to (0,0)

natural(T)=30 line to close

REGION 2

coe=35

START(1,3)

value(T)=100 line to (1,3.5)

to (1+longitud,3.5) to (1+longitud,3)

line to close

monitors

contour(T)

plots

grid(x,y)

vector(dx(T),dy(T))

contour(T)

end

Mallado Disipador

Flujo De calor en el Disipador.

Líneas de temperatura.

En Pdetool

Canaleta:

Disipador

En matlab

Canaleta

dm=60;

t1=32;t2=212;

M=zeros(dm);

n=0;

while n<=1000

for j=2:dm-1

for i=2:dm-1

M(j,i)=(M(j+1,i)+M(j-1,i)+M(j,i+1)+M(j,i-1))/4;

end

end

for i=1:dm

M(i,1)=M(i,2);

M(i,dm)=M(i,dm-1);

end

for i=1:dm;

for j=1:dm

if j==1

M(j,i)=t2;

end

if j>dm/2

if i==j

M(j,i)=32;

end

if dm-i+1==j

M(j,i)=32;

end

end

end

end

n=n+1;

for i=2:dm/2

for j=dm:-1:dm-i+2

M(j,i)=0;

M(j,dm-i+1)=0;

end

end

end

disp(M);

x=linspace(0,4,dm);

...