Advanced Search

Question: How to obtain animated solution of PDE

Question:How to obtain animated solution of PDE

wzelik 70 Maple
animate pdsolve
March 21 2014
0

Hi

I need a temperature distribution inside a barrier during a heating process.
I will be appreciated for any help.

 

wz

How to animete BC using varying temperature in time?  How to obtain animated solution?

> 

restart

Diffusivity coefficent...

> 

a := 0.1e-5:

Thickness of barrier...

> 

L := .2:

Heating curve:
Time in heating curve (in hours form exmaple)...

> 

Time := seq(i, i = 1 .. 10):

Varying temperature in time [K]....

> 

Temp_in_Time := [433.15, 568.15, 703.15, 838.15, 973.15, 1108.15, 1243.15, 1378.15, 1513.15, 1616.15]:

Initial temperature [K]

> 

Tot := 298:
> 

PDE := diff(T(x, t), t) = a*(diff(T(x, t), x, x)):

--->>>

> 

BC1 := {T(0, t) = Temp_in_Time[2], T(L, t) = Temp_in_Time[2], T(x, 0) = Tot}:
> 

sol := pdsolve(PDE, BC1, numeric, timestep = 50):
> 

sol:-plot(t = 3*3600, thickness = 3, colour = red);

 
> 

``



Download heating.mw
Please Wait...