Robotron1715

5 Reputation

One Badge

1 years, 42 days

MaplePrimes Activity


These are questions asked by Robotron1715

numeric.mws

 

Dear Colleagues,
Please give advice - How correct in Maple can be solve the heat nonlinear equation system with perfect thermal contact conditions?

Using procedure dsolve (numeric) Maple cannot determine type of problem (initial or boundary).

****
initial/boundary conditions specified at too many points 1
Error, (in dsolve/numeric/type_check) cannot determine type of problem (initial or boundary), as conditions are given at 4 points
****
 

In attachment Maple file - numeric.

Maybe anyone can get a simple exaple of the maple code for numeric solution of this  task?

 

 


 

restart;

k1:=5; k2:=7;k3:=8; Q:=4;p:=0.001;l:=0.1; d:=0.6;

k1 := 5

k2 := 7

k3 := 8

Q := 4

p := 0.1e-2

l := .1

d := .6

(1)

Eq1:=(k1^2)*diff(T1(x),x$2)=0;

 

Eq1 := 25*(diff(T1(x), `$`(x, 2))) = 0

(2)

Eq2:=(k2^2)*diff(T2(x),x$2)=-Q*Dirac(x)+T2(x)^4;

Eq2 := 49*(diff(T2(x), `$`(x, 2))) = -4*Dirac(x)+T2(x)^4

(3)

Eq3:=(k3^2)*diff(T3(x),x$2)=0;

Eq3 := 64*(diff(T3(x), `$`(x, 2))) = 0

(4)

con1:=T1(0)=T2(0);con2:=T3(p)=T2(p); con3:=k1*(D(T1)(0))=Q-k2*(D(T2)(0)); con4:=k2*(D(T2)(p))=k3*(D(T3)(p)); con5:= T3(l)=300; con6:=T1(-d)=300;

con1 := T1(0) = T2(0)

con2 := T3(0.1e-2) = T2(0.1e-2)

con3 := 5*(D(T1))(0) = 4-7*(D(T2))(0)

con4 := 7*(D(T2))(0.1e-2) = 8*(D(T3))(0.1e-2)

con5 := T3(.1) = 300

con6 := T1(-.6) = 300

(5)

nsol:=dsolve({Eq1, Eq2, Eq3,con1,con2,con3,con4,con5,con6},{T1(x),T2(x),T3(x)}, numeric,bvp);

 

 

initial/boundary conditions specified at too many points 1

Error, (in dsolve/numeric/type_check) cannot determine type of problem (initial or boundary), as conditions are given at 4 points

 

 

 

 


 

Download numeric.mws

Page 1 of 1