Beyond the classical Stefan problem

Abstract

In this thesis we develop and analyse mathematical models describing phase change phenomena linked with novel technological applications. The models are based on modifications to standard phase change theory. The mathematical tools used to analyse such models include asymptotic analysis, similarity solutions, the Heat Balance Integral Method and standard numerical techniques such as finite differences.

In chapters 2 and 3 we study the melting of nanoparticles. The Gibbs-Thomson relation, accounting for melting point depression, is coupled to the heat equations for the solid and liquid and the associated Stefan condition. A perturbation approach, valid for large Stefan numbers, is used to reduce the governing system of partial differential equations to a less complex one involving two ordinary differential equations. Comparison between the reduced system and the numerical solution shows good agreement. Our results reproduce interesting behaviour observed experimentally such as the abrupt melting of nanoparticles. Standard analyses of the Stefan problem impose constant physical properties, such as density or specific heat. We formulate the Stefan problem to allow for variation at the phase change and show that this can lead to significantly different melting times when compared to the standard formulation.

In chapter 4 we study a mathematical model describing the solidification of supercooled liquids. For Stefan numbers larger than unity the classical Neumann solution provides an analytical expression to describe the solidification. For Stefan numbers equal or smaller than unity, the Neumann solution is no longer valid. Instead, a linear relationship between the phase change temperature and the front velocity is often used. This allows solutions for all values of the Stefan number. However, the linear relation is an approximation to a more complex, nonlinear relationship, valid for small amounts of supercooling. We look for solutions using the nonlinear relation and demonstrate the inaccuracy of the linear relation for large supercooling. Further, we show how the classical Neumann solution significantly over-predicts the solidification rate for values of the Stefan number near unity.

The Stefan problem is often reduced to a 'one-phase' problem (where one of the phases is neglected) in order to simplify the analysis. When the phase change temperature is variable it has been claimed that the standard reduction loses energy. In chapters 5 and 6, we examine the one-phase reduction of the Stefan problem when the phase change temperature is time-dependent. In chapter 5 we derive a one-phase reduction of the supercooled Stefan problem, and test its performance against the solution of the two-phase model. Our model conserves energy and is based on consistent physical assumptions, unlike one-phase reductions from previous studies. In chapter 6 we study the problem from a general perspective, and identify the main erroneous assumptions of previous studies leading to one-phase reductions that do not conserve energy or, alternatively, are based on non-physical assumptions. We also provide a general one-phase model of the Stefan problem with a generic variable phase change temperature, valid for spherical, cylindrical and planar geometries.

En aquesta tesi desenvolupem i analitzem models matemàtics que descriuen processos de transició de fase vinculats a noves tecnologies. Els models es basen en modificacions de la teoria estàndard de canvis de fase. Les tècniques matemàtiques per resoldre els models es basen en l'anàlisi asimptòtic, solucions autosimilars, el mètode de la integral del balanç de la calor i mètodes numèrics estàndards com ara el mètode de les diferències finites. En els capítols 2 i 3 estudiarem la transició solid-liquid d'una nanopartícula, acoblant la relació de Gibbs-Thomson, que descriu la depressió de la temperatura de fusió en una superfície corba, amb l'equació de la calor per la fase sòlida i líquida, i la condició de Stefan. Mitjançant el mètode de pertorbacions, per a valors grans del nombre de Stefan, el problema d'equacions en derivades parcials inicial és reduït a un sistema més senzill de dues equacions diferencials ordinàries. La solució del sistema reduït concorda perfectament amb la solució numèrica del sistema inicial d'equacions en derivades parcials. Els resultats confirmen la transició ultra ràpida de sòlid a líquid observada en experiments amb nanopartícules. Els anàlisis estàndards del problema de Stefan consideren propietats físiques com la densitat i el calor específic constants en la fase sòlida i líquida. En aquesta tesi, formularem el problema de Stefan relaxant la condició de densitat constant, el que portarà a diferències molt significatives en els temps totals de fusió al comparar-los amb els temps obtinguts mitjançant la formulació habitual del problema de Stefan. En el capítol 4 estudiarem un model matemàtic que descriu la solidificació de líquids sota-refredats. Per nombres de Stefan més grans que la unitat la solució clàssica de Neumann dona una expressió analítica que descriu el procés. Per valors del nombre de Stefan més petits o iguals que la unitat, la solució de Neumann no es vàlida i, habitualment, per tal de trobar solucions en aquest règim, s'estableix una relació lineal entre la temperatura de canvi de fase i la velocitat del front de solidificació. Aquesta relació lineal però, és una aproximació per sota-refredaments moderats d'una relació no lineal més complexa. En aquest capítol, buscarem solucions del problema incorporant la relació no lineal al model, i demostrarem la poca precisió a l'utilitzar l'aproximació lineal. A més a més, veurem com la solució de Neumann sobreestima de manera significativa la velocitat del procés per valors propers a la unitat. Habitualment, el problema de Stefan és simplificat a un problema d'una fase (on una de les fases es desestimada) per tal de reduir la dificultat de l'anàlisi. En casos on la temperatura de transició de fase és variable s'ha reclamat que la simplificació d'una fase del problema no conserva l'energia. En els capítols 5 i 6, examinarem les diferents reduccions d'una fase del problema de Stefan en el cas on la temperatura de transició depèn del temps. En el capítol 5 derivarem el problema de Stefan d'una fase associat a la solidificació de líquids sota-refredats, i compararem la solució del sistema resultant amb la solució del problema de Stefan estàndard de dues fases. A diferència dels models d'una fase descrits en estudis previs, el nostre model reduït conserva l'energia i està basat en suposicions físiques consistents. En el capítol 6 estudiarem el problema des d'una perspectiva més general i identificarem les suposicions errònies d'estudis previs que porten a la no conservació de l'energia o que, alternativament, estan basades en suposicions físiques poc consistents. A més a més, derivarem un model d'una fase amb temperatura de canvi de fase variable, vàlid per geometries esfèriques, cilíndriques i planes.