## Diffusion From Spheres

### Associated Hands-On Teaching Module: http://www.che.utah.edu/module/?p_id=35

This simulation is meant to demonstrate aspects of mass transfer in radial coordinates A summary of this system and its governing equations may be found below the simulation. Set up your initial conditions at the bottom of the applet and press play to run the reaction simulation in "real-time". Press fast forward to run as quickly as your computer is capable. Mouse-over each input to see a description of that option or variable.

**1-D Mass Transfer in a Sphere / Drug Delivery Simulation**

University of Utah - Department of Chemical Engineering

by Anthony Butterfield

May Not Work if Your Browser Isn't Up to Date

**Simplifying to Diffusion Out of a Single Bead into a Semi-Infinite Media**

While this simulation models a more complex system with diffusion into a finite media with changing concentration, the traditional 1-D radial diffusion problem may be simulated by making the volume of water to beads very large. For example, one could reduce the number of beads to 1 and raise the liquid volume to a kiloliter. When this is done and dt is small enough to avoid numerical error, the concentration profile across the bead should be the same as that given by the analytical solution, *C ^{a}*(

*r*), which may be plotted using the appropriate check-box below the plot.

**Comparing With Your Own Data**

- Running your alginate bead experiment in the lab and collect data from your spectrophotometer over time until the signal levels off.
- There will be a lot of data points that need to be converted from voltage to concentration but you can’t go directly from V to C.
- You will need to convert your voltage to transmittance and then to absorbance, using the methods described in your experiments with the spectrophotometers.
- Transmittance to Absorbance

A = -log_{10}(T) - Because A is proportional to
*C*, to get the percentage released,*Y*, you need only divide*A*by the maximum*A*measured,*if the experiment ran until the signal leveled off*. If you did not let the experiment run until the maximum was found, you can convert to*C*using your calibration for fuschin, and use a calculated value of*C*to determine_{inf}*Y*. - You will not want to use every data point, as it will slow the simulation down too much. Instead, you can take, for example, every 50th data point in a Matlab vector of time,
*t*, using the command

>> tnew=t(1:50:end) - Convert your data to a comma delimited form appropriate for the simulation. You can use the command

>> sprintf('%f,' , tnew) - Paste the comma delimited time and Y into the text boxes at the bottom of the simulation.
- Put the appropriate bead radius and other initial conditions in for your particular beads.
- Run the simulation and find a diffusivity that creates a release profile, Y, that best matches your data.
- Press eject to get that data for plotting later.