The solution of the classical one-dimensional Stefan problem predicts that in time t the melt front goes as s(t)∼t12. In the presence of heterogeneity, however, anomalous behavior can be observed where the time exponent n12. In such a case, it may be appropriate to write down the governing equations of the Stefan problem in terms of fractional order time (1≥β>0) and space (1≥α>0) derivatives. Here, we present sharp and diffuse interface models of fractional Stefan problems and discuss available analytical solutions. We illustrate that in the fractional time case (β<1), a solution of the diffuse interface model in the sharp interface limit will not coincide with the solution of the sharp interface counterpart; negating a well know result of integer derivative Stefan problems. The paper concludes with the development of an implicit time stepping numerical solution for the diffuse interface fractional Stefan model. Results from this solution are verified with available analytical solutions.
|Original language||English (US)|
|Number of pages||9|
|Journal||International Journal of Heat and Mass Transfer|
|State||Published - Jul 2014|
Bibliographical noteFunding Information:
The author acknowledges support from the National Science Foundation, United Sates, through Grant EAR-1318593 , Generalized Transport Models in Earth Surface Dynamics. The author is also very greatful to Fede Falcini (Consiglio Nazionale delle Ricerche, Rome, Italy), Roberto Garra (La Sapienza, Rome, Italy) and Piotr Rybka (Instytut Matematyki Stosowanej i Mechaniki, Uniwersytet Warszawski, Poland) for many insights and fruitful discussions.
- Fractional derivative
- Numerical solution
- Stefan problem