Drug Diffusion in Solid Square Bar

Given:

$-$ Initial concentration of Drug A in water $($C$\backslash$infty $=8\mathrm{.}0\backslash$times $10^(-3)$ kg mol$/$m$^3)$

$-$ Dimensions of the solid square bar: $0\mathrm{.}2$ m $\backslash$times $0\mathrm{.}2$ m $\backslash$times $10$ m

$-$ Initial concentration in the square bar $($Ci $=2\mathrm{.}0\backslash$times $10^(-3)$ kg mol$/$m$^3)$

$-$ Effective diffusivity $($Deff $=4\mathrm{.}5\backslash$times $10^(-3)$ cm$^2/$s$)$

$-$ Convection mass transfer coefficient $($kc $=9\backslash$times $10^(-5)$ cm$/$s$)$

$-$ Grid size $(\backslash$Delta x $=\backslash$Delta y $=5$ cm$)$

c$)$ Finding Concentrations at t $=10$ hours

Use the finite difference equations to find the concentrations C$,$ C$,...,$ C$$ at t $=10$ hours. Utilize Excel or MATLAB for computation.

d$)$ Estimating Time for Saturation

Estimate the time required for the bar to reach $95\%$ saturation of Drug A $($Cave $=7\mathrm{.}6\backslash$times $10^(-3)$ kg mol$/$m$^3).$

$$