Klockjj

Ask a Question

Create a Post

One Badge

Member for: 8 years, 255 days

Contact Klockjj

MaplePrimes Activity

These are replies submitted by Klockjj

@Rouben Rostamian ah thank y...

Answered: Klockjj 10

February 05 2016

0 0

@Rouben Rostamian ah thank you very much for spotting that big mistake for me that equation was suppose to be using an exp() instead of ln().

That solved most of my problems however now i am having issues with creating an array for output. I am getting this signularity error this time, would you happen to have any ideas as to what may be causing this?

I am trying to generate a list of values for Cprop(z), T(z) and P(z) from 0 to 5 in .5 increments.

>	kinetic parameter

(24)

>	radial disperssion Coeff

Dr := Vs*dp/(9*(1+19.4*(dp/dt)^2))

(25)

>	thermal conductivity calculations,w/m-k

Kbs := kfluid*(`εb`+(1-`εb`)/(2/3*(1+kfluid/ksolid)))

(26)

Kbd := `εb`*Cpair*Gs*dp/(9*(1+19.4*(dp/dt)^2))

(27)

(28)

>	Heat of reaction (deltaH), Find heat of formation for each reactant and product--> dHrxn = heat of formation(product)-heat of formation(reactant)

Error, missing operator or `;`

>	heat of formation = A + BT + CT^2

(29)

(30)

(31)

(32)

>	Solving nodified Ergun Equation

f := (1-`εb`)*(1.75+(150*(1-`εb`))*`μ_air`/(dp*Gs))/`εb`

(33)

>	Determining U of the wall and Ueffective

Uwall := 1/(1/hi+ln(dto/dt)/(2*kpipe)+1/ho)

(34)

Ueff := 1/(1/Uwall+(1/2)*dt/(4*KB))

(35)

0.4117046714e-2*(T(z)+273.15)*(diff(Cprop(z), z))/P(z) = -exp(19.837-13636/(T(z)+273.15))*Cprop(z)

(36)

1.4645*(diff(T(z), z)) = -146.7382702*exp(19.837-13636/(T(z)+273.15))*Cprop(z)-666.3198268*T(z)+415217.2000

(37)

ode3 := diff(P(z), z) = -f*Gs^2/(rho*dp)

(38)

ode := Gs*(diff(Cprop(z), z))/rho = -K*Cprop(z), Gs*Cpair*(diff(T(z), z)) = -K*Cprop(z)*`ΔHrxn`-4*Ueff*(T(z)-(350+273.15))/dt, diff(P(z), z) = -f*Gs^2/(rho*dp)

0.4117046714e-2*(T(z)+273.15)*(diff(Cprop(z), z))/P(z) = -exp(19.837-13636/(T(z)+273.15))*Cprop(z), 1.4645*(diff(T(z), z)) = -146.7382702*exp(19.837-13636/(T(z)+273.15))*Cprop(z)-666.3198268*T(z)+415217.2000, diff(P(z), z) = -4.361346783*(T(z)+273.15)/P(z)

(39)

(40)

(41)

(42)

(43)

(44)

Warning, cannot evaluate the solution further right of .10547158e-2, probably a singularity

array( 1 .. 2, 1 .. 1, [( 2, 1 ) = (array( 1 .. 11, 1 .. 4, [( 1, 4 ) = (HFloat(350.0)), ( 1, 1 ) = (HFloat(0.0)), ( 1, 3 ) = (HFloat(2.5)), ( 1, 2 ) = (HFloat(0.03)) ] )), ( 1, 1 ) = (array( 1 .. 4, [( 1 ) = (z), ( 2 ) = (Cprop(z)), ( 3 ) = (P(z)), ( 4 ) = (T(z)) ] )) ] )