Saha

70 Reputation

4 Badges

1 years, 313 days

MaplePrimes Activity


These are questions asked by Saha

Please help! How to solve this?

pp-velocity.mw

Could you Please Help me,the Maple code for Plot this equations in any Numerical Method

restart

with(PDETools, declare):

with(LinearAlgebra):

with(PolynomialIdeals):

with(plots):

declare(f(x));

f(x)*`will now be displayed as`*f

(1)

eq1 := diff(f(x), x, x, x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(f(x), x, x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta = 0;

diff(diff(diff(f(x), x), x), x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(diff(f(x), x), x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta = 0

(2)

eq2 := K[nf]*(diff(theta(x), x, x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x)) = 0;

K[nf]*(diff(diff(theta(x), x), x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x)) = 0

(3)

bcs := F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0;

F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0

(4)

a1 := [M = 1, alpha = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, `ρβ`[fl] = 20939.1, K[nf] = .6842, eta = 0];

[M = 1, alpha = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, `ρβ`[fl] = 20939.1, K[nf] = .6842, eta = 0]

 

[Pr = 6.2, M = 1, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, `ρβ`[fl] = 20939.1, K[nf] = .6842, K[f] = .613, eta = 0]

(5)

b1 := subs(a1, eq1);

diff(diff(diff(f(x), x), x), x)+.5325197465*f(x)*(diff(diff(f(x), x), x)) = 0

(6)

c1 := subs(a2, eq2);

1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x)) = 0

(7)

d1 := dsolve({b1, bcs}, {bcs, c1}, numeric); d1(0)

Error, (in dsolve/numeric/process_input) invalid argument: {1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x)) = 0, F(0) = 0, F(1) = 0, F(10) = 1, Theta(0) = 1, Theta(10) = 0}

 

d1(0)

(8)

``

Download  Numerical Method.mw

5 6 7 8 Page 7 of 8