Advanced Search

Question: how I can Pdsolve this equation?

Question:how I can Pdsolve this equation?

torabi 85 Maple 2018
differential-equation pde pdsolve
February 20 2019
0

is possible to solvethis equation via maple?

hank you

EQUATION2.mw
 

> 

restartNULL
> 

alpha := 1.2*10^(-4); Betaa := 4.0*log(2); J := 13.4; delta := 15.3*10^(-9); tp := 10^(-13); tq := 8.5*10^(-12); tu := 90.0*10^(-12); kapa := 315; r0 := 2.0*10^(-7); Lx := 5.0*10^(-7); Ly := 5.0*10^(-7); Lz := 1.0*10^(-7); a := 0.7e-1*(Betaa/Pi)^.5*J/(15.3*10^(-22)); bb := exp(-((10^(-7)*x-(1/2)*Lx)^2+(10^(-7)*y-(1/2)*Ly)^2)/(2*r0^2)); print(aa = a); Q := a*exp(-z*10^(-7)/delta)*exp(-1.88*abs(t-2*tp)/tp)*bb

0.1200000000e-3

  

4.0*ln(2)

  

13.4

  

0.1530000000e-7

  

1/10000000000000

  

0.8500000000e-11

  

0.9000000000e-10

  

315

  

0.2000000000e-6

  

0.5000000000e-6

  

0.5000000000e-6

  

0.1000000000e-6

  

0.6917775548e21*ln(2)^.5

  

exp(-0.1250000000e14*((1/10000000)*x-0.2500000000e-6)^2-0.1250000000e14*((1/10000000)*y-0.2500000000e-6)^2)

  

aa = 0.6917775548e21*ln(2)^.5

  

0.6917775548e21*ln(2)^.5*exp(-6.535947712*z)*exp(-0.1880000000e14*abs(t-1/5000000000000))*exp(-0.1250000000e14*((1/10000000)*x-0.2500000000e-6)^2-0.1250000000e14*((1/10000000)*y-0.2500000000e-6)^2)

 (1)
> 

(diff(U(x, y, z, t), t)+tq*(diff(U(x, y, z, t), t, t)))/alpha = diff(U(x, y, z, t), x, x)+diff(U(x, y, z, t), y, y)+tu*(diff(U(x, y, z, t), x, x, t)+diff(U(x, y, z, t), y, y, t))+tu*(diff(U(x, y, z, t), z, z, t))+(Q+tq*(diff(Q, t)))/kapa

8333.333333*(diff(U(x, y, z, t), t))+0.7083333333e-7*(diff(diff(U(x, y, z, t), t), t)) = diff(diff(U(x, y, z, t), x), x)+diff(diff(U(x, y, z, t), y), y)+0.9000000000e-10*(diff(diff(diff(U(x, y, z, t), t), x), x))+0.9000000000e-10*(diff(diff(diff(U(x, y, z, t), t), y), y))+0.9000000000e-10*(diff(diff(diff(U(x, y, z, t), t), z), z))+0.2196119222e19*ln(2)^.5*exp(-6.535947712*z)*exp(-0.1880000000e14*abs(t-1/5000000000000))*exp(-0.1250000000e14*((1/10000000)*x-0.2500000000e-6)^2-0.1250000000e14*((1/10000000)*y-0.2500000000e-6)^2)-0.3509398517e21*ln(2)^.5*exp(-6.535947712*z)*abs(1, t-1/5000000000000)*exp(-0.1880000000e14*abs(t-1/5000000000000))*exp(-0.1250000000e14*((1/10000000)*x-0.2500000000e-6)^2-0.1250000000e14*((1/10000000)*y-0.2500000000e-6)^2)

 (2)
> 

``

 

Boundary condition:

> 

U(0, y, z, t) = 300; U(Lx, y, z, t) = 300; U(x, 0, z, t) = 300; U(x, Ly, z, t) = 300; U(x, y, 0, t) = 300; U(x, y, Lz, t) = 300

#####################################

INITIAL CONDITIONS:

 

> 

U(x, y, z, 0) = 300; (D[1](U))(x, y, z, 0) = 0

(D[1](U))(x, y, z, 0) = 0

 (3)

NULL

 

``


 

Download EQUATION2.mw

 

 
Please Wait...