Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Atomic hydrogen diffusion in novel magnesium nanostructures : the impact of incorporated subsurface carbon atoms

Ab initio Density Functional Theory (DFT) calculations are performed to study the diffusion of atomic hydrogen on a Mg(0001) surface and their migration into the subsurface layers. A carbon atom located initially on a Mg(0001) surface can migrate into the sub-surface layer and occupy a fcc site, with charge transfer to the C atom from neighboring Mg atoms. The cluster of postively charged Mg atoms surrounding a sub-surface C is then shown to facilitate the dissociative chemisorption of molecular hydrogen on the Mg(0001) surface, and the surface migration and subsequent diffusion into the subsurface of atomic hydrogen. This helps rationalize the experimentally-observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

A study of the diffusion characteristics of normal, delipidized and relipidized articular cartilage using magnetic resonance imaging

This paper assesses the capacity to provide semipermeability of the synthetic layer of surface-active phospholipids created to replace the depleted surface amorphous layer of articular cartilage. The surfaces of articular cartilage specimens in normal, delipidized, and relipidized conditions following incubation in dipalmitoyl-phosphatidylcholine and palmitoyl-oleoyl-phosphatidylcholine components of the joint lipid mixture were characterized nanoscopically with the atomic force microscope and also imaged as deuterium oxide (D2O) diffused transiently through these surfaces in a magnetic resonance imaging enclosure. The MR images were then used to determine the apparent diffusion coefficients in a purpose-built MATLAB®-based algorithm. Our results revealed that all surfaces were permeable to D2O, but that there was a significant difference in the semipermeability of the surfaces under the different conditions, relative to the apparent diffusion coefficients. Based on the results and observations, it can be concluded that the synthetic lipid that is deposited to replace the depleted SAL of articular cartilage is capable of inducing some level of semipermeability.

On the growth process of Cu2ZnSn(S,Se)4 absorber layer formed by selenizing Cu–ZnS–SnS precursors and its photovoltaic performance

Sputtering and subsequent sulfurization(orselenization)is one of the methods that have been extensively employed to fabricate Cu2ZnSn(S,Se)4 (CZTSSe) thin films. However, there are limited reports on the effect of precursor stacking order of the sputtered source materials on the properties of the synthesized CZTSSe films. In this work,the morphology and crystallization process of the CZTSSe films which were prepared by selenizing Cu–ZnS–SnS precursor layers with different stacking sequences and the adhesion property between the as-synthesized CZTSSe layer and Mosubstrate have been thoroughly investigated. It has been found that the growth of CZTSSe material and the morphology of the film strongly depend on the location of Culayer in the precursor film. The formation of CZTSSe starts from the diffusion of Cu–Se to Sn(S,Se)layert o form Cu–Sn–(S,Se) compound,followed by the reaction with Zn(S,Se). The investigation of themorphology of the CZTSSe films has shown that large grains are formed in the film with the precursor stacking order of Mo/SnS/ZnS/Cu,which is attributed to a bottom-to-top growth mechanism. In contrast, the film made from a precursor with a stacking sequence of Mo/ZnS/ SnS/Cu is mainly consisted of small grains due to a top-to-bottom growth mechanism. The best CZTSSe solar cell with energy conversion efficiency of3.35%has been achieved with the selenized Mo/ZnS/ SnS/Cu film, which is attributed to a good contact between the absorber layer and the Mosubstrate.

The effect of thermal radiation on the natural convection boundary layer flow over a wavy horizontal surface

In this article, natural convection boundary layer flow is investigated over a semi-infinite horizontal wavy surface. Such an irregular (wavy) surface is used to exchange heat with an external radiating fluid which obeys Rosseland diffusion approximation. The boundary layer equations are cast into dimensionless form by introducing appropriate scaling. Primitive variable formulations (PVF) and stream function formulations (SFF) are independently used to transform the boundary layer equations into convenient form. The equations obtained from the former formulations are integrated numerically via implicit finite difference iterative scheme whereas equations obtained from lateral formulations are simulated through Keller-box scheme. To validate the results, solutions produced by above two methods are compared graphically. The main parameters: thermal radiation parameter and amplitude of the wavy surface are discussed categorically in terms of shear stress and rate of heat transfer. It is found that wavy surface increases heat transfer rate compared to the smooth wall. Thus optimum heat transfer is accomplished when irregular surface is considered. It is also established that high amplitude of the wavy surface in the boundary layer leads to separation of fluid from the plate.

Interdiffusion and Growth of the Superconductor Nb3Sn in Nb/Cu(Sn) Diffusion Couples

The confusion over the growth rate of the Nb3Sn superconductor compound following the bronze technique is addressed. Furthermore, a possible explanation for the corrugated structure of the product phase in the multifilamentary structure is discussed. Kirkendall marker experiments are conducted to study the relative mobilities of the species, which also explains the reason for finding pores in the product phase layer. The movement of the markers after interdiffusion reflects that Sn is the faster diffusing species. Furthermore, different concentrations of Sn in the bronze alloy are considered to study the effect of Sn content on the growth rate. Based on the parabolic growth constant at different temperatures, the activation energy for the growth is determined.

nvestigations on the Stability and Extinction of a Laminar Diffusion Flame Over a Porous Flat Plate

The limits of stability and extinction of a laminar diffusion flame have been experimentally studied in a two-dimensional laminar boundary layer over a porous flat plate through which n-pentane vapour was uniformly injected. The stability and extinction boundaries are mapped on a plot of free stream oxidant velocity versus fuel injection velocity. Effects of free stream temperature and of dilution of fuel and oxidant on these boundaries have been examined. The results show that there exists a limiting oxidant flux beyond which the diffusion flame cannot be sustained. This limiting oxidant flux has been found to depend_on the free stream oxygen concentration, fuel concentration and injection'velocity of the fuel.

A semi-analytical solution for multilayer diffusion in a composite medium consisting of a large number of layers

Diffusion in a composite slab consisting of a large number of layers provides an ideal prototype problem for developing and analysing two-scale modelling approaches for heterogeneous media. Numerous analytical techniques have been proposed for solving the transient diffusion equation in a one-dimensional composite slab consisting of an arbitrary number of layers. Most of these approaches, however, require the solution of a complex transcendental equation arising from a matrix determinant for the eigenvalues that is difficult to solve numerically for a large number of layers. To overcome this issue, in this paper, we present a semi-analytical method based on the Laplace transform and an orthogonal eigenfunction expansion. The proposed approach uses eigenvalues local to each layer that can be obtained either explicitly, or by solving simple transcendental equations. The semi-analytical solution is applicable to both perfect and imperfect contact at the interfaces between adjacent layers and either Dirichlet, Neumann or Robin boundary conditions at the ends of the slab. The solution approach is verified for several test cases and is shown to work well for a large number of layers. The work is concluded with an application to macroscopic modelling where the solution of a fine-scale multilayered medium consisting of two hundred layers is compared against an “up-scaled” variant of the same problem involving only ten layers.

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

Study on the Growth of Nb3Sn Superconductor in Cu(Sn)/Nb Diffusion Couple

Nb3Sn growth following the bronze technique, (i.e. by interdiffusion between Cu(Sn) alloy (bronze) and Nb) is one of the important methodologies to produce this superconductor. In this study, we have addressed the confusion over the growth rate of the Nb3Sn phase. Furthermore, a possible explanation for the corrugated layer in the multifilamentary structure is discussed. Kirkendall marker experiments were conducted to study the relative mobilities of the species, which also explained the reason for finding pores in the product phase layer. Based on the parabolic growth constant at different temperatures, the activation energy for the growth is determined. We have further explained the dramatic increase in the growth rate of the prod

Flow structure and heat transfer characteristics behind a diaphragm in the presence of a diffusion flame

The characteristics of the separated flow behind a diaphragm over a burning surface are investigated experimentally. This complex problem of practical significance involving recirculation, blowing and combustion reactions is studied in a two-dimensional combustion tunnel. The flame structure, recirculation patterns and heat transfer to the surface are presented for a range of values of free stream and fuel injection velocities as well as for different heights of the diaphragm. The trends of heat transfer vs axial distance are shown to be similar to those resulting from a non-reactive heated stream with a diaphragm. Treating the case of a boundary layer diffusion flame as that corresponding to the zero height of the diaphragm, the heat transfer augmentation due to recirculation is estimated. It is found that at considerable downstream distances (xfh > 3), the heat transfer rates with diaphragm overtake the rates from a developing boundary layer case. Flow visualization studies with particle track photography show that there are many similarities between the reactive and the non-reactive cases.

Diffusion path during internal displacement reactions in multi-component oxides: Reaction between Fe and (Co,Mg)TiO3 solid solution at 1273 K

The reaction between Fe foil and a disc of ilmenite solid solution (Co-0.48 Ni-0.52) TiO3 was studied at 1273 K. At the metal/oxide interface, the displacement reaction, Fe + (Co,Mg)TiO3 = Co + (Fe,Mg)TiO3 occurs, resulting in an ilmenite solid solution containing three divalent cations. Ferrous ions diffuse into the oxide solid solution and cause the precipitation of Co-Fe alloy as discrete particles inside the oxide matrix. The morphology of the product layer was characterized by SEM. Only two phases, alloy and ilmenite, were detected in the reaction zone. This suggests that the local flux condition imposed by ilmenite stoichiometry (Co + Fe + Mg):Ti = 1:1] was satisfied during the reactive diffusion: (J(Co) + J(Fe) + J(Mg)) = J(Ti). The composition of the alloy and the oxide was determined using EPMA as a function of distance in the direction of diffusion. Although Mg does not participate in the displacement reaction, its composition in the ilmenite phase was found to be position dependent inside the reaction zone. The up-hill diffusion of inert Mg is caused by the development of chemical potential gradients as a result of displacement reaction. The evolution of composition gradients inside the reaction zone and the diffusion path in a ternary composition diagram of the system CoTiO3-FeTiO3-MgTiO3 are discussed. (C) 2010 Elsevier B.V. All rights reserved.

Reactive Diffusion Between Vanadium and Silicon

Interdiffusion study is conducted in the V-Si system to determine integrated diffusion coefficients of the phases. Activation energy values are calculated from the experiments conducted at different temperatures. The average values are found to be 208, 240 and 141 kJ/mol, respectively, for the V(3)Si, V(5)Si(3) and VSi(2) phases. The low activation energy for the VSi(2) phase indicates very high concentration of defects or the significant contribution from the grain boundary diffusion. The error in calculation of diffusion parameters from a very thin phase layer in a multiphase diffusion couple is discussed. Further the data available in the literature in this system is compared and the problems in the indirect methodology followed previously to calculate the diffusion parameters are discussed.

Growth and consumption rates of the phase layers during interdiffusion in a diffusion couple with finite end member

The relations for the growth and consumption rates of a layer with finite thickness as an end member and the product phases in the interdiffusion zone are developed. We have used two different methodologies, the diffusion based and the physico-chemical approach to develop the same relations. We have shown that the diffusion based approach is rather straightforward; however, the physico-chemical approach is much more versatile than the other method. It was found that the position of the marker plane becomes vague in the second stage of the interdiffusion process in pure A thin layer/B couple, where two phases grow simultaneously.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.