Abstract
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.
Original language | English |
---|---|
Pages (from-to) | 376-388 |
Number of pages | 13 |
Journal | Applied Mathematics and Computation |
Volume | 219 |
Issue number | 1 |
DOIs | |
Publication status | Published - 15 Sep 2012 |
Keywords
- Artificial parameter
- Asymptotics
- Boundary immobilisation
- Controlled drug release
- Heat balance integral method
- Keller box scheme
- Solvent penetration
- Stefan problem