Dear all,

Now, mathematics really got me struck during my graduation project on drying of milk particles in a fluidized bed. For this project I am trying to calculate the profile for surface evaporation according to Crank (1955). The analytical solution is quite diffucult and displays lots of numerical errors in Excel.

Being practically desperate I hope someone could help me out with this one.

I will first explain the total problem and thereafter the practical issue.

We consider a moist particle of 200 mu within a an air stream. The general formula is for the concentration on time = t and position = r:

The 's are the roots of the following equation:

In the formula the letters have the following meaning:

concentration at time = t and place r in particle

initial concentration within the particle at t = 0

concentration around the particle at

1/2 particle diameter

Diffusion coefficient

mass transition coefficient between particle and air

Euler number

Now the real problem. The formula will be used to describe the profile of C over r and t. Resturcturing the formula clearly shows that there are several limits.

Some of these limits are not a big problem, like: or YET...

Since the diffusion coefficient is quite low to , the 's go to practically .(not really pi) And this makes the go to zero which means that goes to infinity.

The question to you now is whether you could help me out solving the following limit analytically:

Please not that the l'hopital rule can probably not be applied since only the denominator is 0. And note that is not really ! Has someone got another option?

Or even better to solve analytically:

with the as input for the 's

Could please someone help me out? I'm struggelling for weeks now. Excel turns mad when I try to solve this equation, since Excel uses taylor approximations for the sinus. Even Maple 11 could not help me with the rootfinding?

Wish you all the best

jurjen