An Optimal Method for Diffusion Parameters of Nonlinear Diffusion Problem of Drug Releasing in 2 D-Disc Device by Separate Variable Method

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.


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 [1][2][3][4][5].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 [6][7][8][9][10].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 [6-8, 11, 12].However, in the practice application, many diffusion processes are subjected to the nonlinear partial differential equations [6][7][8][9][10][11][12].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.

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: where the drug concentration is uniform in the device and zero in liquid for the initial condition at  = 0; that is, 0 is the total mass of the diffusion material,   is the volume of Ω 1 , the diffusion coefficient () is the function of the concentration, (, , ) 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.
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.
Using the linear diffusion processes to replace the nonlinear diffusion process in the time, it is easy to obtain the following linear equation: Using the polar coordinate system to replace (3), the above equations can be changed as follows: From the above deduction, the nonlinear diffusion equation of the drug releasing in the section of [  ,  +1 ] and [  ,  +1 ] can be obtained as follows: Because the length of the time and the range of space by the discretion are very small, ( −1, ) can be considered as In the section of [  ,  +1 ], the linear equations can be represented as follows: where ] can be obtained and get If  = 0 and  = 1,  01 (, , 0) = (, , 0).
To solve ( − 1) ×  linear equations, by the separation variable method, we can get the solutions of the nonlinear diffusion equation as follows: Equation ( 9) can be changed into the following formula: where  (1)   ( = 1, 2, . ..) are the positive roots of  1 () and Γ( + 1) = !,  is the positive integer and

Least Square Method Based on Separate Variable Method for Solving Optimal Control Problem
Let where  = ( 1 ( 1 , . . .,   ),  2 ( 1 , . . .,   ), . ..,   ( 1 , . . .,   ))  , and For an initial diffusion parameter point  1 ( where (, , ) ,  ∈ (  ,  +1 ) ,  = 1, 2, . . ., , (16) equation ( 15) can be written in the following form: To solve this problem, suppose  () is the th approximation and let the function   () be the Taylor expansion function at  () ; the minimal point  (+1) and the ( + 1)th approximation can be computed by the iterative method.The iterative formula is deduced in detail as follows. Set We use () to replace () and compute the minimal point of () to estimate the function ().The least square problem: min () can be solved as follows.Set Equation ( 19) can be written as In order to search the stable point of (), set Taking  and  into the above formula, there is Moving the right  Τ  () to the left hand in the above equation, the following equation can be obtained: Obviously, this is a linear algebraic equation about the function value and the first order partial derivative   ( () ) at the point  () .If matrix  is the full column rank,    is a symmetry positive matrix.Therefore, there exists (  ) −1 .
We can get the stable point of () by ( 19) Set  (+1) as the ( + 1)th approximation of stable point ().

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  = 0, the inner concentration is 138.1553g/dm 2 and the outer is 0 g/dm 2 .After some diffusion process in a period of time, the inner and outer concentration will be equal and be 3.9468 g/dm 2 ; it is In the numerical example, suppose () =  +  where  and  are the unknown constants.We firstly suppose the given diffusion parameter (, ) 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: 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 ,  and the optimal values (, ) of the diffusion parameter are depicted in Figures 2, 3 and 4.  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  converges very fast because their values become from −3.0152 × 10 −4 to 3.1408 × 10 −6 and from Table 2 and Figure 4  varies from 1.1501 × 10 −3 to −6.1503 × 10 −6 .From the last column in Table 2, the optimized values ,  also become very fast from (−2.0152 × 10 −4 , 1.2502 × 10 −3 ) to (2.9856 × 10 −4 , 3.0283 × 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 ,  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 [9] 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 [9] 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 ,  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.

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.

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

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

1 Figure 3 : 5 Figure 4 :
Figure 3: The error value of  on the initial value for (0.0001 and 0.0001).

Figure 6 :
Figure 6: Computed data in each iteration and experimental data.

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

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