Asymptotic, numerical and approximate techniques for a free boundary problem arising in the diffusion of glassy polymers

作者：Mitchell S L^*; O' Brien S B G

来源：Applied Mathematics and Computation, 2012, 219(1): 376-388.

DOI：10.1016/j.amc.2012.06.026

• Controlled drug release
• Solvent penetration
• Asymptotics
• Boundary immobilisation
• Stefan problem
• Keller box scheme
• Heat balance integral method
• Artificial parameter

This paper considers approximate solution methods for a one dimensional Stefan problem describing solvent diffusion in glassy polymers. Similar to the classic Stefan problem, the region initially has zero thickness and so must be analysed carefully before performing a numerical computation. A small-time analysis gives the correct starting solution which is then incorporated into the second order accurate Keller box finite difference scheme. We also consider a detailed analysis of small and large time expansions, as well as the large control parameter limit, and show that our generalised approach enables us to obtain higher order terms than given in previous studies. Finally, we apply the combined integral method (CIM) to this problem, which is a refinement of the popular heat balance integral method (HBIM), and compare both the CIM and asymptotic solutions to the numerical results.

• 出版日期2012-9-15