Advanced Search

Question: ode bvp fluid problem

Question:ode bvp fluid problem

faisal 15 Maple 17
dsolve numeric differential-equations bvp
November 01 2017
0

how may i encounter this problem,

Please help me on this ode solve (dsolve) wit
 

> 

restart
> 

with(plots):
> 

ps := 8933:
> 

Eq1 := 8*(11-10*d)*(eta*(diff(f(eta), `$`(eta, 3)))+diff(f(eta), `$`(eta, 2)))/((1-phi)^2.5*(1-phi+phi*ps/pf))+4*f(eta)*(diff(f(eta), `$`(eta, 2)))+m*(1-4*(diff(f(eta), eta))*(diff(f(eta), eta))) = 0;

8*(11-10*d)*(eta*(diff(diff(diff(f(eta), eta), eta), eta))+diff(diff(f(eta), eta), eta))/((1-phi)^2.5*(1+7.958981045*phi))+4*f(eta)*(diff(diff(f(eta), eta), eta)) = 0

 (1)
> 

Eq2 := 2*knf*(eta*(diff(theta(eta), `$`(eta, 2)))+diff(theta(eta), eta))/(Pr*kf*(1-phi(eta)+phi*ps*cps/(pf*cpf)))+f(eta)*(diff(theta(eta), eta))-n*(diff(f(eta), eta))*theta(eta) = 0;

186.4367280*(eta*(diff(diff(theta(eta), eta), eta))+diff(theta(eta), eta))/(1-phi(eta)+.8253667629*phi)+f(eta)*(diff(theta(eta), eta)) = 0

 (2)
> 

Va := [.1, .1, .1]:
> 

`Vφ` := [0.5e-1, 0.5e-1, 0.5e-1]:
> 

etainf := 8:
> 

bcs := f(a) = 0, (D(f))(a) = 0, theta(a) = 1, (D(f))(etainf) = 1/2, theta(etainf) = 0;

f(a) = 0, (D(f))(a) = 0, theta(a) = 1, (D(f))(8) = 1/2, theta(8) = 0

 (3)
> 

dsys := {Eq1, Eq2, bcs}:
> 

for i to 3 do for h to 3 do phi := `Vφ`[i]; a := Va[h]; dsol[h][i] := dsolve(dsys, numeric, continuation = d); print(a); print(phi); print(dsol[h][i](0)) end do end do

Error, (in fproc) unable to store 'YP[3]' when datatype=float[8]

  
> 

``


 

Download Grosan.mw

h bcs
Please Wait...