Advanced Search

Question: How do I plot multiple pde numerical solutions on one plot using maple?

Question:How do I plot multiple pde numerical solutions on one plot using maple?

tobsam 15 Maple 17
plot pde numeric
December 30 2013
1

Please, I solved a pde system of equation problem numerically, using maple 17.

But I dont know how to plot multiple solutions on one graph.

I want to vary one of the parameters....

e.g Pr=0.71, Pr=7, Pr=10 where other parameters are kept constant

 

My working is attachedtobi_msc_solution.mw

restart

M := 1:

pde1 := diff(u(y, t), t)+Typesetting:-delayDotProduct(S, diff(u(y, t), y))-2*k^2*u(y, t) = diff(u(y, t), y, y)+theta(y, t)+Typesetting:-delayDotProduct(N, C(y, t))+Typesetting:-delayDotProduct(M, u(y, t))+u(y, t)/K:

                pde2 := theta(y, t)+t*(diff(theta(y, t), t))+S*(diff(theta(y, t), y)) = (diff(theta(y, t), y, y))/Pr-Typesetting:-delayDotProduct(alpha, theta(y, t)):

pde3 := C(y, t)+t*(diff(C(y, t), t))+S*(diff(C(y, t), y)) = (diff(C(y, t), y, y))/Sh-Typesetting:-delayDotProduct(R, C(y, t)):

PDE := {pde1, pde2, pde3}:

IBC := {C(0, t) = 1, C(1, t) = 0, C(y, 0) = 0, u(0, t) = 0, u(1, t) = 0, u(y, 0) = 0, theta(0, t) = 1, theta(1, t) = 0, theta(y, 0) = 0}:

pds := pdsolve(PDE, IBC, numeric)

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

 (1)

pds:-plot[display](u(y, t), t = .5, linestyle = "solid", colour = "blue", legend = "Pr=0.71", title = "Velocity Profile", labels = ["y", "theta"])

  

``


Download tobi_msc_solution.mw

 

Please, Any help will be gracefully appreciated

 
Please Wait...