MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 796812 10.1155/2014/796812 796812 Research Article An Optimal Method for Diffusion Parameters of Nonlinear Diffusion Problem of Drug Releasing in 2D-Disc Device by Separate Variable Method Li Youyun 1, 2 Ouyang Jinhui 1 Qin Jiaohua 3 Gao Yingli 2 Loukopoulos Vassilios C. 1 State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin Beijing 100038 China 2 College of Highway Engineering Changsha University of Science and Technology Hunan, 410004 China csust.edu.cn 3 College of Computer Central South University of Forestry and Technology Changsha 410004 China csuft.edu.cn 2014 2332014 2014 14 11 2013 16 02 2014 16 02 2014 23 3 2014 2014 Copyright © 2014 Youyun Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

An optimization control model and the corresponding computational method drawing the diffusion parameters of the nonlinear problem for the drug releasing in the 2D-disc device were given in this paper. Firstly, based on the nonlinear diffusion equation of the drug releasing in the 2D-disc device, we used the linear diffusion problem to discrete the nonlinear diffusion problem with the discrete space and the discrete time. Then, by the separate variable method, the solution of the linear problem was given. Next, the least square method based on the separate variable idea (LSMSV) was used to estimate the nonlinear appropriate diffusion parameters. Finally, a numerical example was presented to show that the control model and the numerical method were valid for computing the diffusion coefficient of the nonlinear problem for the drug releasing in the 2D-disc device.

1. Introduction

In engineering fields, there exist many diffusion processes in many fields such as geomechanics engineering, biomedical science, civil engineering, water pollution, and soil engineering . In order to simulate the diffusion processes to obtain their merits, it is important to draw the effective diffusiveness. There are many models for the simulation of the diffusion processes. Most of them are the nonlinear or linear models. For the linear models, most optimal control problems governed by the diffusion equations arose in many scientific and engineering applications such as the water pollution problems and the drug releasing fields . There are many various techniques for the identification for the effective diffusiveness based on the linear models. These techniques are based on either empirical or semiempirical models from drug delivery mechanisms or on analytic solutions of the diffusion equation in 2D or in the special cases [68, 11, 12].

However, in the practice application, many diffusion processes are subjected to the nonlinear partial differential equations . In order to illustrate the nonlinear diffusion processes in many fields, many nonlinear models are applied to estimate the properties. For many nonlinear diffusion fields, the diffusion coefficients are the functions of the diffusion concentration. Therefore, computing the diffusion parameters is mainly to determinate the parameters of the coefficient function called the diffusion parameters function. The diffusion parameter function of the concentration is considered as the main element to control the diffusion processes. Therefore, many researches were given to determine the parameter function to illustrate the diffusion processes. Many nonlinear optimal models drawing the diffusion parameters depended on both the lab technology and the shape of the container [9, 10, 12, 13]. In order to compute the diffusion parameters, many scientists and mathematicians provided some optimal methods to compute the diffusion parameters [9, 10, 12]. Most of them cost a lot of computing time and computing memory. Even some of them such as the different method or the finite element method cost more than one week [9, 14]. Therefore, in order to save the computing time and the computing memory, in our published papers, we had provided some numerical methods based on the separate variable method to compute the diffusion parameters of the linear process in the sphere device. For these reasons, in this paper, we will also propose a new numerical optimal method (the least square method based on the separate variable idea) to extract the diffusion parameters from the nonlinear diffusion problems.

The next two sections will give the nonlinear diffusion problem in the 2D-disc device based on the drug releasing property and the discrete method. Section 3 is devoted to providing the least square method based on the separate variable idea for the optimal control model of the nonlinear diffusion equation system governing the drug releasing process. In Section 4, the numerical example is presented to demonstrate the feasibility and the validity, the convergence of the model, and computing method. Finally, we will discuss the model and computing method and give the future work for the different control model.

2. The Nonlinear Diffusion Problem in the 2D-Disc Device and the Discrete Linear Problem

The nonlinear diffusion process in 2D-disc device is governed by the following partial differential equation: (1) C ( x , y , t ) t - · ( D ( C ) C ) = 0 , t > 0 , ( x , y ) Ω , C ( x , y , t ) n = 0 , t > 0 , ( x , y ) Ω , C ( x , y , 0 ) = H ( x , y ) , ( x , y ) Ω , where the drug concentration is uniform in the device and zero in liquid for the initial condition at t = 0 ; that is, (2) H ( x , y , 0 ) = { M 0 V d , ( x , y ) Ω 1 , 0 , ( x , y ) Ω Ω 1 ; M 0 is the total mass of the diffusion material, V d is the volume of Ω 1 , the diffusion coefficient D ( C ) is the function of the concentration, C ( x , y , t ) is the drug releasing concentration, Ω 1 is the small disc containing the drug, and Ω is the big disc containing the liquid shown in Figure 1.

Large disc Ω containing liquid and small disc Ω 1 containing drug.

In order to solve the nonlinear diffusion process by the separate variable method, it is necessary to use some linear equations to replace the nonlinear equation (1). The discretion process is given as follows.

Time Discretion. Setting Δ t i = t i + 1 - t i , i = 1,2 , , N , the time [ T 0 , T 1 ] of the diffusion process is divided into many sections [ t i , t i + 1 ] and the concentration C i is computed in [ t i , t i + 1 ] .

Space Discretion. Setting Δ r j = r j + 1 - r j , j = 1,2 , , K - 1 , R 1 = r 1 < r 2 < < r k = R 2 , the radius [ R 1 , R 2 ] is divided into many intervals and the diffusion concentration in the section of [ r 1 , r 2 ] , [ r 2 , r 3 ] , , [ r k - 1 , r k ] is obtained. So the nonlinear equation system (1) can be changed into ( k - 1 ) * N linear equations.

Using the linear diffusion processes to replace the nonlinear diffusion process in the time, it is easy to obtain the following linear equation: (3) C i t = D ( C i - 1 ) ( 2 C i x 2 + 2 C i y 2 ) .

Using the polar coordinate system to replace (3), the above equations can be changed as follows: (4) C i t = D ( C i - 1 ) ( 2 C i r 2 + 1 r C i r ) .

From the above deduction, the nonlinear diffusion equation of the drug releasing in the section of [ t i , t i + 1 ] and [ r j , r j + 1 ] can be obtained as follows: (5) C i j t = D ( C ( i - 1 ) j ) ( 2 C i j r 2 + 1 r C i j r ) , i = 1,2 , , N , j = 1,2 , , K - 1 .

Because the length of the time and the range of space by the discretion are very small, D ( C i - 1 , j ) can be considered as a constant in [ t i , t i + 1 ] × [ r j , r j + 1 ] .

In the section of [ t i , t i + 1 ] , the linear equations can be represented as follows: (6) C i j t = D ( C i - 1 , j ) ( 2 C i j r 2 + 1 r C i j r ) , t i t t i + 1 , m m l l i = 1 , , N , 0 < r < r j + 1 , j = 1 , , K - 1 , C i j r | r = r j + 1 = 0 , C i j ( r , t i ) = { M i - 1 , j V d j , 0 < r < r j , 0 , r j < r < r j + 1 , where V d j = π ( r j + 1 2 - r j 2 ) . Therefore, the concentrations C i j of this equation in [ t i , t i + 1 ] × [ r j , r j + 1 ] can be obtained and get (7) M i j = 0 2 π r j r j + 1 C i ( r , θ , t ) d r d θ - 0 2 π r j - 1 r j C i - 1 ( r , θ , t ) d r d θ , m m m m m m m m m m m m m l l l i = 1,2 , , N , j = 1,2 , , K . If i = 0 and j = 1 , C 01 ( r , θ , 0 ) = H ( r , θ , 0 ) .

To solve ( K - 1 ) × N linear equations, by the separation variable method, we can get the solutions of the nonlinear diffusion equation as follows: (8) C i j = M i j V d j + n = 1 2 M i j r j V d j r j + 1 μ n J 0 2 ( μ n ) J 1 ( r j μ n r j + 1 ) × exp ( - D ( C ( i - 1 ) j ) ( μ n r j + 1 ) 2 t i ) J 0 ( μ n r j + 1 r ) , m m m m m m l l i = 1,2 , , N , j = 1,2 , , K , (9) C i = j = 1 N l = 1 i C l j ( r , θ , t ) , r ( r j , r j + 1 ) , i = 1,2 , , K .

Equation (9) can be changed into the following formula: (10) C i = j = 1 N ( l = 1 i ( M l j V d + n = 1 2 M l j r j V d r j + 1 μ n ( 1 ) J 0 2 ( μ n ( 1 ) ) J 1 ( r j μ n r j + 1 ) × exp ( - D ( C ( l - 1 ) j ) ( μ n r j + 1 ) 2 t l ) × J 0 ( μ n r j + 1 r ) M l j V d + n = 1 2 M l j r j V d r j + 1 μ n ( 1 ) J 0 2 ( μ n ( 1 ) ) J 1 ( r j μ n r j + 1 ) ) l = 1 i ( M l j V d + n = 1 2 M l j r j V d r j + 1 μ n ( 1 ) J 0 2 ( μ n ( 1 ) ) J 1 ( r j μ n r j + 1 ) ) , r ( r 2 , r K ) , i = 1,2 , , K - 1 , where μ n ( 1 ) ( n = 1,2 , ) are the positive roots of J 1 ( x ) and Γ ( m + 1 ) = m ! , m is the positive integer and (11) J 1 ( μ n ) = m = 0 ( - 1 ) m 1 Γ ( m + 1 ) Γ ( m + 2 ) ( μ n 2 ) 2 m + 1 , J 0 ( μ n R 2 r ) = m = 0 ( - 1 ) m 1 Γ ( m + 1 ) Γ ( m + 1 ) ( μ n r 2 R 2 ) 2 m .

3. Control Problem Drawing Nonlinear Diffusion Parameters Problem 1.

Search the coefficient functions D ( a 1 , , a L , C ) to satisfy (12) min D > 0 { ( M T 1 - M T 1 0 ) 2 + ( M T 2 - M T 2 0 ) 2 + + ( M T e - M T e 0 ) 2 } , where M T 1 0 , M T 2 0 , , M T e 0 are the given experimental data and M T 1 , M T 2 , , M T e are the computed data by the following equations: (13) M T i = t k T i , j = 2 K M k j .

4. Least Square Method Based on Separate Variable Method for Solving Optimal Control Problem

Let (14) E ( D ) = ( M T 1 - M T 1 0 ) 2 + ( M T 2 - M T 2 0 ) 2 + + ( M T e - M T e 0 ) 2 = ( M - M * ) T ( M - M * ) , where M = ( M T 1 ( a 1 , , a L ) , M T 2 ( a 1 , , a L ) , , M T e ( a 1 , , a L ) ) T , and M * = ( M T 1 0 , M T 2 0 , , M T e 0 ) T .

For an initial diffusion parameter point D 1 ( a 1 , , a L ) , Problem 1 can be solved iteratively by the following deduction. If a 1 i , , a L i are the i th approximation and δ a 1 i , , δ a L i are the i th increment of a 1 , , a L , respectively, in each step, an increment δ D i ( δ a 1 i , , δ a L i ) will be computed as follows. To minimize E ( D i + δ D i ) with δ D i ( δ a 1 i , , δ a L i ) , let X = ( x 1 , , x L ) = ( a 1 , , a L ) , f T i ( X ) = M T i - M T i 0 , F ( X ) = E ( D ) , so F ( X ) = E ( X ) and δ D ( δ X ) = δ D i ( δ a 1 i , , δ a L i ) . If set (15) f i ( X ) = M T i - M T i 0 = 0 2 π r 2 r N C T i ( r , θ , T i ) d r d θ - M i 0 i = 1,2 , , e , where (16) C T i = j = 2 N l = 1 max { l T i } C l j ( r , θ , t ) , r ( r j , r j + 1 ) , i = 1,2 , , e , equation (15) can be written in the following form: (17) F ( x ) = i = 1 e f i 2 ( X ) .

To solve this problem, suppose X ( k ) is the k th approximation and let the function f i ( X ) be the Taylor expansion function at X ( k ) ; the minimal point X ( k + 1 ) and the ( k + 1 ) th approximation can be computed by the iterative method. The iterative formula is deduced in detail as follows.

Set (18) ϕ i ( X ) = f i ( X ( k ) ) + f i ( X ( k ) ) T ( X - X ( k ) ) = f i ( X ( k ) ) T X - [ f i ( X ( k ) ) T X ( k ) - f i ( X ( k ) ) ] , m m m m m m m m m m m m m m m m m i = 1,2 , , e , (19) ϕ ( X ) = i = 1 e ϕ i 2 ( X ) .

We use ϕ ( X ) to replace F ( X ) and compute the minimal point of ϕ ( X ) to estimate the function F ( X ) . The least square problem: min ϕ ( X ) can be solved as follows. Set (20) A = [ f 1 ( X ( k ) ) T f e ( X ( k ) ) T ] = [ f 1 ( X ( k ) ) x 1 f 1 ( X ( k ) ) x L f e ( X ( k ) ) x 1 f e ( X ( k ) ) x L ] , (21) B = [ f 1 ( X ( k ) ) T X ( k ) - f 1 ( X ( k ) ) f e ( X ( k ) ) T X ( k ) - f e ( X ( k ) ) ] = A X ( k ) - f ( k ) , f ( k ) = [ f 1 ( X ( k ) ) f 2 ( X ( k ) ) f e ( X ( k ) ) ] , f i ( X ( k ) ) x l = r 2 r N 0 2 π j = 1 N l = 1 l T i ( M l j V d + n = 1 2 M l j r j V d r j + 1 μ n ( 1 ) J 0 2 ( μ n ( 1 ) ) m m m m m m m m m m × J 1 ( r j μ n r j + 1 ) m m m m m m m m m m × exp ( ( μ n r j + 1 ) 2 t l - D ( X ( k ) , C ( l - 1 ) , j ) m m m m m m m m m m m m m m l × ( μ n r j + 1 ) 2 t l ) m m m m m m m m m m × J 0 ( μ n r j + 1 r ) D ( X ( k ) , C l - 1 , j ) x l m m m m m m m m m m × ( μ n r j + 1 ) 2 t l ) d θ d r .

Equation (19) can be written as (22) ϕ ( X ) = ( A X - B ) T ( A X - B ) = X T A T A X - 2 B T A X + B T B .

In order to search the stable point of ϕ ( X ) , set (23) ϕ ( X ) = 2 A T A X - 2 A T B = 0 .

Taking A and B into the above formula, there is (24) A T A X = A T ( A X ( k ) - f ( k ) ) .

Moving the right A Τ A X ( k ) to the left hand in the above equation, the following equation can be obtained: (25) A T A ( X - X ( k ) ) = - A T f ( k ) .

Obviously, this is a linear algebraic equation about the function value and the first order partial derivative f i ( X ( k ) ) at the point X ( k ) . If matrix A is the full column rank, A T A is a symmetry positive matrix. Therefore, there exists ( A T A ) - 1 . We can get the stable point of ϕ ( X ) by (19) (26) X ( k + 1 ) = X ( k ) - ( A T A ) - 1 A T f ( k ) .

Set X ( k + 1 ) as the ( k + 1 )th approximation of stable point F ( X ) .

Algorithm 2 (diffusion parameter of nonlinear process).

Step  1. Give the initial point X ( 1 ) = ( a 1 , ( 1 ) , , a n , ( 1 ) ) and set t 1 = 0 , M 1 = 0 ,   k = 1 , f 1 ( X ( k ) ) = 0 and the control error ε > 0 , go to Step  2.

Step  2. Compute the drug concentration M i j ( X ( k ) ) and C i j ( X ( k ) ) at t i in [ r j , r j + 1 ] by the formula (7) and (8), obtain the mass in all sections [ r j , r j + 1 ] at the time t i , and go to Step  3.

Step  3. Compute the mass M T i ( X ( K ) ) , i = 1,2 , , e . at the outer container Ω Ω 1 by the formula (13); go to Step  4.

Step  4. According to the formula (15), compute f i ( X ( k ) ) = M T i ( X ( k ) ) - M T i 0 ; get the vector: f ( k ) = [ f 1 ( X ( k ) ) f 2 ( X ( k ) ) f e ( X ( k ) ) ] . Compute the first order partial derivative a i j = f i ( X ( k ) ) / x j , i = 1,2 , , e , j = 1,2 , , L by formula (21). Get the matrix A k = ( a i j ) e × L by the formula (20); go to Step  5.

Step  5. Compute X ( k + 1 ) = X ( k ) - ( A T A k ) - 1 A T f ( k ) ; go to Step  6.

Step  6. If X ( k + 1 ) - X ( k ) L 2 ε , set X * = X ( k + 1 ) and get the optimum coefficient ( a 1 , , a L ) = X * , go to Step  7; otherwise, let k = k + 1 , go to Step  2.

Step  7. Output the optimal diffusion ( a 1 , , a L ) and stop the algorithm.

5. The Numerical Example

To investigate the feasibility and the validity of the proposed scheme, a numerical example is given in this section. The cylinder device for the drug releasing is divided into the large and small cylinder structures where the small container includes the drug and the large container is filled with the liquid. The radius of small and large cylinder devices are 0.4800  dm and 2.8399  dm, respectively, and there is 100 g of drug in small disc vessel. When t = 0 , the inner concentration is 138.1553  g / d m 2 and the outer is 0  g / d m 2 . After some diffusion process in a period of time, the inner and outer concentration will be equal and be 3.9468  g / d m 2 ; it is (27) C ( r , 0 ) = { M 0 π R 1 2 = 138.1553 , 0 < r < R 1 , 0 , R 1 < r < R 2 , C ( r , T e ) = 3.9468 0 < r < R 2 .

In the numerical example, suppose D ( C ) = a C + b where a and b are the unknown constants. We firstly suppose the given diffusion parameter ( a , b ) to be ( 0.0003 , 0.0003 ) then compute the diffusion process to obtain the computed data as the experiment data of the concentration shown in Table 1. Next, the above optimal model and the optimal method to estimate the diffusion parameter based on the optimal control model with the separation variable method of the drug releasing in the 2D-disc device are used to illustrate the feasibility and the validity of the LMSV for the nonlinear diffusion process. The following formula is given to compute the total error in order to estimate the convergence rate: (28) total  error   ( k ) = j = 1 e ( M T j - M T j k 0 ) 2 , δ a = a k - a k - 1 , δ b = b k - b k - 1 .

Diffusion qualities at the different times with diffusion parameters (0.0003 and 0.0003).

Time (second) M Time (second) M Time (second) M
0 0 160 50.0660 300 65.0085
20 42.1695 180 52.2352 320 67.0014
40 42.2755 200 54.4820 340 68.9612
60 42.7675 220 56.6823 360 70.8841
80 43.7504 240 58.8197 1000 95.6678
100 45.5094 260 60.9205 2000 97.1432
120 47.8335 280 62.9829 3000 97.1432

Supposing the space interval to be 0.012 dm and the time interval 20 seconds; getting 25 terms in the Bessel function and by the least square method algorithm, the optimal points and their error values in each optimal step are given in the Table 2. In order to illustrate the convergence of the least square algorithm, the optimal increment δ a , δ b and the optimal values ( a , b ) of the diffusion parameter are depicted in Figures 2, 3 and 4.

The different errors and optimization value based on initial value (0.0001 and 0.0001).

Iterated number Total error δ a δ b Optimization value ( a , b )
1 2.9804 × 104 −3.0152 × 10−4 1.150 1  × 10−3 (−2.0152 × 10−4, 1.2502 × 10−3)
2 8.2626 × 103 1.4436 × 10−5 1.7564 × 10−4 (−1.8708 × 10−4, 1.4258 × 10−3)
3 2.9995 × 103 1.3370 × 10−4 −2.3741 × 10−4 (−5.3387 × 10−4, 1.1884 × 10−3)
4 9.03742 × 102 1.5343 × 10−4 −3.6483 × 10−4 (1.0004 × 10−4, 8.2360  × 10−4)
5 2.3343 × 102 1.1775 × 10−4 −3.0877 × 10−4 (2.1779 × 10−4, 5.1482  × 10−4)
6 4.6960 × 10 6.0218 × 10−5 −1.6193 × 10−4 (2.7801 × 10−4, 3.5290 × 10−4)
7 5.0935 1.7408 × 10−5 −4.3911 × 10−5 ( 2.9542  × 10−4, 3.0898 × 10−4)
8 4.0114 3.1408 × 10−6 −6.1503 × 10−6 (2.9856 × 10−4, 3.0283  × 10−4)

The overall error with iterative times on the initial value for (0.0001 and 0.0001).

The error value of a on the initial value for (0.0001 and 0.0001).

The error value of b on the initial value for (0.0001 and 0.0001).

In order to illustrate the convergence of the algorithm for diffusion parameter of nonlinear process by the least square method by separate variables (LSMSV), we discuss the convergence data as follows: from the second column in Table 2 and Figure 3, the increments δ a converges very fast because their values become from −3.0152 × 10−4 to 3.1408 × 10−6 and from Table 2 and Figure 4   δ b varies from 1.1 501 × 10 −3 to −6.1503 × 10−6. From the last column in Table 2, the optimized values a , b also become very fast from (−2.0152 × 10−4, 1.2502 × 10−3) to (2.9856 × 10−4, 3.02 83 × 10 −4) by only eight iterated steps. Therefore, from Table 2 and the error value in Figures 2, 3 and 4, the data show the error and the increment δ a , δ b convergent stately. It is easy to illustrate the convergence of the algorithm for diffusion parameter of nonlinear process by least square method by separate variables idea.

In order to test the convergent velocity, we obtain the computing time for the optimal parameters by the algorithm in the numerical examples. The computing time is 1 minute and 56 seconds by the algorithm for diffusion parameter of nonlinear process. In order to test the merits of the algorithm for diffusion parameter of nonlinear process by LSMSV, it is hard for us to use the algorithm in the paper  to compute the parameter values because the computed time is very long. In addition, from the computing processes, it is also easy to understand why the computing velocity becomes very high. Because the algorithm in this paper only computes some of the polynomial functions in each iteration to cost only less than 2 minutes, however, the algorithm in  needs to solve millions linear algebra equations in each iteration for the nonlinear diffusion process to cost more than one day. Therefore, from the computing convergent time and the computational theory, we can obtain the conclusion that the convergent velocity of algorithm for diffusion parameter of nonlinear process by LSMSV is very fast.

In order to estimate the validity of the algorithm for diffusion parameter of nonlinear process by LSMSV, the total error values between the experiment values and the computed values in this paper are shown in Figure 5 and the first column in Table 2. The optimal computed drug mass depending on the different optimal diffusion parameters and the experiment mass are depicted in Figure 6. From the first column in Table 2, the errors of the optimized values a , b and the parameters are only 0.0156 × 10−4 and 0.213 × 10−4, the error between the optimal computed drug mass 100 g and the experiment mass data 97.1432 g is only 2.8568 g, and the relative error of the mass is only 2.9%. The error result shows the algorithm for diffusion parameter of nonlinear process by LSMSV being valid to extract the diffusion parameter of the drug releasing in the disc devices for the nonlinear diffusion process. From Figures 5 and 6 and the last optimal diffusion parameter value ( 2.9856 × 10 - 4 , 3.0283 × 10 - 4 ) , the algorithm is valid to extract the diffusion parameter of the drug releasing for the nonlinear diffusion process in the 2D-disc devices. From the numerical example, comparing the data of all iterative steps of the optimization values and the experimental data, it is easy to get the conclusion that the algorithm is the convergence and effectiveness to extract the diffusion parameters of the nonlinear drug releasing in the 2D-disc devices.

The optimization increments of ( a , b ) .

Computed data in each iteration and experimental data.

6. Conclusion

In this paper, we propose an optimal method to extract the diffusion parameters of the nonlinear diffusion process of the drug releasing in the 2D-disc device based on the separation variable method with discrete time and discrete space. The numerical result given in the previous section demonstrates the feasibility and validity of this algorithm for diffusion parameter of nonlinear process by LSMSV. The effectiveness of this optimal control model to estimate the diffusion parameters for the nonlinear drug releasing in the 2D-disc device is also discussed. How to establish the theorem of the algorithm is our future work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (NSFC) Grants (no. 11072041 and 61202496), by State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin (IWHRSKL201205), by State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology (GZ1005), by China Postdoctoral Science Foundation (20100480944 and 2012T50692), and by Hunan Provincial Natural Science Foundation of China (13JJ3070, 13JJ2031).

Siepmann J. Ainaoui A. Vergnaud J. M. Bodmeier R. Calculation of the dimensions of drug-polymer devices based on diffusion parameters Journal of Pharmaceutical Sciences 1998 87 7 827 832 2-s2.0-0031825269 10.1021/js980006a Crawford G. J. Hicks C. R. Lou X. Vijayasekaran S. Tan D. Mulholland B. Chirila T. V. Constable I. J. The Chirila Keratoprosthesis: phase I human clinical trial Ophthalmology 2002 109 5 883 889 2-s2.0-0036237911 10.1016/S0161-6420(02)00958-2 Wang S. Lou X. Novel mathematicsl models for extracting effective drug diffusivity of porous PHEMA hydrogels: from a planar matrix into a finite external volume Proceeding of the ICEAM May 2011 Changsha, China Grassi M. Grassi G. Mathematical modelling and controlled drug delivery: matrix systems Current Drug Delivery 2005 2 1 97 116 2-s2.0-17844390397 10.2174/1567201052772906 Lou X. Vijayasekaran S. Sugiharti R. Robertson T. Morphological and topographic effects on calcification tendency of pHEMA hydrogels Biomaterials 2005 26 29 5808 5871 2-s2.0-20444386223 Mellott M. B. Searcy K. Pishko M. V. Release of protein from highly cross-linked hydrogels of poly(ethylene glycol) diacrylate fabricated by UV polymerization Biomaterials 2001 22 9 929 941 2-s2.0-0035342614 10.1016/S0142-9612(00)00258-1 Zhu J. Zeng Q. A mathematical theoretical frame for control of air pollution Science in China. Series D 2002 32 864 870 Price P. E. Jr. Wang S. Romdhane I. H. Extracting effective diffusion parameters from drying experiments AIChE Journal 1997 43 8 1925 1934 2-s2.0-0031208917 Li Y. Xiang Z. Wang S. Identifying time-dependent drug diffusion parameters in the Cylindrical tube by the finite difference method Proceeding of the ICEAM May 2011 Changsha, China Parra-Guevara D. Skiba Y. N. Elements of the mathematical modeling in the control of pollutants emissions Ecological Modelling 2003 167 3 263 275 2-s2.0-0348103731 10.1016/S0304-3800(03)00191-1 Fu J. C. Hagemeir C. Moyer D. L. An unified mathematical model for diffusion from drug polymer composite tablets Journal of Biomedical Materials Research 1976 10 5 743 758 2-s2.0-0016994959 Hukka A. The effective diffusion coefficient and mass transfer coefficient of nordic softwoods as calculated from direct drying experiments Holzforschung 1999 53 5 534 540 2-s2.0-0033284507 Hicks C. R. Crawford G. J. Lou X. Tan D. T. Cornea replacement using a synthetic hydrogel cornea, AlphaCor: device, preliminary outcomes and complications Eye 2003 17 385 392 10.1038/sj.eye.6700333 Li Y. Xiang Z. Xiang X. Wang S. A computer algorithm for optimizing to extract effective diffusion coefficients of drug delivery from cylinders Information Technology Journal 2010 9 8 1647 1652 2-s2.0-78651360139