Question: PDE system of sixth order

hi

how pdsolve 2D couple  non linear differential equations which attached below?

thanks

2D.mw

PD1 := 12.6000000000000*(diff(U(x, theta), x, x, x, x))-7500*(diff(U(x, theta), x, x))+10.2112755389544*(diff(U(x, theta), x, x, theta, theta))+7.22165554279476*(diff(V(x, theta), x, x, x, theta))-3730.19397871630*(diff(V(x, theta), x, theta))+4.57316679658261*(diff(V(x, theta), x, theta, theta, theta))+36.0000000000000*(diff(W(x, theta), x, x, x))-9375.00*(diff(W(x, theta), x))-15.6731205946742*(diff(W(x, theta), x, theta, theta))-1947.26649812618*(diff(U(x, theta), theta, theta))+1.41357763467822*(diff(U(x, theta), theta, theta, theta, theta)):

PD2 := 3.52500000000000*(diff(V(x, theta), x, x, x, x))-3150.00*(diff(V(x, theta), x, x))+5.05278814097746*(diff(V(x, theta), theta, theta, theta, theta))-4972.65371594662*(diff(V(x, theta), theta, theta))+10.2112755389544*(diff(V(x, theta), x, x, theta, theta))+7.22165554279476*(diff(U(x, theta), x, x, x, theta))+4.57316679658261*(diff(U(x, theta), x, theta, theta, theta))-2536.53190552708*(diff(U(x, theta), x, theta))+97.2834589649212*(diff(W(x, theta), x, x, theta))+30.9917088694028*(diff(W(x, theta), theta, theta, theta))-31273.9463175575*(diff(W(x, theta), theta)):

PD3 := -0.168000000000000e-2*(diff(W(x, theta), x, x, x, x, x, x))+16.9000000000000*(diff(W(x, theta), x, x, x, x))-0.319161728473364e-2*(diff(W(x, theta), x, x, x, x, theta, theta))-37.5000000000000*(diff(W(x, theta), x, x))+796.276013656496*(diff(W(x, theta), x, x, theta, theta))-0.202111525639098e-2*(diff(W(x, theta), x, x, theta, theta, theta, theta))-36.0000000000000*(diff(U(x, theta), x, x, x))+9375.00*(diff(U(x, theta), x))+15.6731205946742*(diff(U(x, theta), x, theta, theta))-97.2834589649212*(diff(V(x, theta), x, x, theta))+31273.9463175575*(diff(V(x, theta), theta))-30.9917088694028*(diff(V(x, theta), theta, theta, theta))+2.06625000000000*10^5*W(x, theta)+313.462411893484*(diff(W(x, theta), theta, theta))+6.77715234781900*(diff(W(x, theta), theta, theta, theta, theta))-0.426628729281504e-3*(diff(W(x, theta), theta, theta, theta, theta, theta, theta))-400/((1-h3*f3(x))*ln(10-10*h3*f3(x))^2)-500/(1-h3*f3(x))^4:

BC := {D[1](U)*(0, theta) = 0, D[1](U)*(1, theta) = 0, D[1](U)*(x, 0) = 0, D[1](U)*(x, 1) = 0, D[1](V)*(0, theta) = 0, D[1](V)*(1, theta) = 0, D[1](V)*(x, 0) = 0, D[1](V)*(x, 1) = 0, D[1](W)*(0, theta) = 0, D[1](W)*(1, theta) = 0, D[1](W)*(x, 0) = 0, D[1](W)*(x, 1) = 0, D[2](W)*(0, theta) = 0, D[2](W)*(1, theta) = 0, D[2](W)*(x, 0) = 0, D[2](W)*(x, 1) = 0, U(0, theta) = 0, U(1, theta) = 0, U(x, 0) = 0, U(x, 1) = 0, V(0, theta) = 0, V(1, theta) = 0, V(x, 0) = 0, V(x, 1) = 0, W(0, theta) = 0, W(1, theta) = 0, W(x, 0) = 0, W(x, 1) = 0}

``

 

Download 2D.mw

Please Wait...