Ice formation on a smooth or rough cold surface due to the impact of a supercooled water droplet

Ice accretion is considered in the impact of a supercooled water droplet on a smooth or rough solid surface, the roughness accounting for earlier icing. In this theoretical investigation the emphasis and novelty lie in the full nonlinear interplay of the droplet motion and the growth of the ice surface being addressed for relatively small times, over a realistic range of Reynolds numbers, Froude numbers, Weber numbers, Stefan numbers and capillary underheating parameters. The Prandtl number and the kinetic under-heating parameter are taken to be order unity. The ice accretion brings inner layers into play forcibly, affecting the outer flow. (The work includes viscous effects in an isothermal impact without phase change, as a special case, and the differences between impact with and without freezing.) There are four main findings. First, the icing dynamically can accelerate or decelerate the spreading of the droplet whereas roughness on its own tends to decelerate spreading. The interaction between the two and the implications for successive freezings are found to be subtle. Second, a focus on the dominant physical effects reveals a multi-structure within which restricted regions of turbulence are implied. The third main finding is an essentially parabolic shape for a single droplet freezing under certain conditions. Fourth is a connection with a body of experimental and engineering work and with practical findings to the extent that the explicit predictions here for ice-accretion rates are found to agree with the experimental range.

Non-classical thermal transport and phase change at the nanoscale

Premi Extraordinari de Doctorat, promoció 2018-2019. Àmbit de CiènciesFor 200 years, Fourier’s law has been used to describe heat transfer with excellent results. However, as technology advances, more and more situations arise where heat conduction is not well described by the classical equations. Examples are applications with extremely short time scales such as ultra fast laser heating, or very small length scales such as the heat conduction through nanowires or nanostructures in general. In this thesis we investigate alternative models which aim to correctly describe the non-classical effects that appear in extreme situations and which Fourier’s law fails to describe. A popular approach is the Guyer-Krumhansl equation and the framework of phonon hydrodynamics. This formalism is particularly appealing from a mathematical point of view since it is analogous to the Navier-Stokes equations of fluid mechanics, and from a physical point of view, since it is able to describe the physics in a simple and elegant way. In the first part of the thesis we use phonon hydrodynamics to predict the size-dependent thermal conductivity observed experimentally in nanostructures such as nanowires or thin films. In particular, we show that the Guyer-Krumhansl equation is suitable to capture the dependence of the thermal conductivity on the size of the physical system under consider- ation. During the modelling process we use the analogy with fluids to incorporate a slip boundary condition with a slip coefficient that depends on the ratio of the phonon mean free path to the characteristic size of the system. With only one fitting parameter we are able to accurately reproduce experimental observations corresponding to nanowires and nanorods of different sizes. The second part of the thesis consists of studying the effect of the non-classical fea- tures on melting and solidification processes. We consider different extensions and in- corporate them into the mathematical description of a solidification process in a simple, one-dimensional geometry. In chapter 5 we employ an effective Fourier law which replaces the original thermal conductivity by a size-dependent expression that accounts for non-local effects. In chapter 6 we use the Maxwell-Cattaneo and the Guyer-Krumhansl equations to formulate the Maxwell-Cattaneo-Stefan and the Guyer-Krumhansl-Stefan problems respec- tively. After performing a detailed asymptotic analysis we are able to reduce both models to a system of two ordinary differential equations and obtain excellent agreement with the cor- responding numerical solutions. In situations near Fourier resonance, which is a particular case where non-classical effects in the Guyer-Krumhansl model cancel each other out, the solidification kinetics are very similar to those described by the classical model. However, in this case we see that non-classical effects are still observable in the evolution of the heat flux through the solid, which suggests that this is a quantity which is more convenient to determine the presence of these effects in phase change processes.La llei de Fourier ha estat una peça clau per a descriure la conducció de calor des de que fou proposada fa gairabé 200 anys. No obstant, a mesura que avança la tecnologia ens hi trobem més sovint amb situacions on les equacions clàssiques perden la seva validesa. En aquesta tesi investiguem alguns models alternatius que tenen com a objectiu descriure la conducció de calor en situacions on la llei de Fourier no és aplicable. Un model que s'ha aconseguit establir com un extensió vàlida de la llei de Fourier és l'equació de Guyer i Krumhansl i el marc de la hidrodinàmica de fonons derivat d'aquesta. Es tracta d'un model particularment interessant, ja que les equacions són anàlogues a les equacions per a fluids dins de la hidrodinàmica clàssica. A la primera part de la tesi considerem aquesta equació per a descriure la conducció de calor estàtica per nanofibres de seccions transversals circulars i rectangulars. En particular, calculem una conductivitat tèrmica efectiva i trobem que és possible reproduïr els resultats experimentals amb un sol paràmetre d'adjust. En el cas de nanofibres cilíndriques, no és necessari cap paràmetre d'adjust si es consideren unes certes condicions de vora per al flux. Una conseqüència d'haver de considerar extensions per a la llei de Fourier és que s'ha d'estudiar l'efecte que tenen aquests canvis en la descripció de processos de canvi de fase. En la segona part de la tesi investiguem els efectes que tenen diversos models sobre la solidificació d'un líquid unidimensional. Al capítol 5 estudiem el cas en el que considerem la conductivitat tèrmica com a una funció de la mida del sòlid i que incorpora característiques que són importants quan el tamany del sòlid és comparable a les longituds característiques dels fonons, mentres que al capítol 6 considerem l'equació de Guyer i Krumhansl dels capítols anteriors. En ambdós casos, un anàlisi asimptòtic ens permet reduïr la complexitat del problema i proposar models reduïts formats per un parell de'equacions diferencials ordinàries.Award-winningPostprint (published version

Interfacial dynamics and pinch-off singularities for axially symmetric Darcy flow

We study a model for the evolution of an axially symmetric bubble of inviscid fluid in a homogeneous porous medium otherwise saturated with a viscous fluid. The model is a moving boundary problem that is a higher-dimensional analogue of Hele-Shaw flow. Here we are concerned with the development of pinch-off singularities characterised by a blow-up of the interface curvature and the bubble subsequently breaking up into two; these singularities do not occur in the corresponding two-dimensional Hele-Shaw problem. By applying a novel numerical scheme based on the level set method, we show that solutions to our problem can undergo pinch-off in various geometries. A similarity analysis suggests that the minimum radius behaves as a power law in time with exponent $\alpha = 1/3$ just before and after pinch-off has occurred, regardless of the initial conditions; our numerical results support this prediction. Further, we apply our numerical scheme to simulate the time-dependent development and translation of axially symmetric Saffman-Taylor fingers and Taylor-Saffman bubbles in a cylindrical tube, highlighting key similarities and differences with the well-studied two-dimensional cases.Comment: 16 pages, 16 figure

Mathematics of moving boundary problems in diffusion

This thesis is concerned with the development, generahzation and apphcation of a formal series technique for classical one-dimensional moving boundary diffusion problems. The solution procedure consists of two major steps. Firstly, the introduction of a boundary fixing transformation, which fixes the moving boundary and simplifies the transformed equations. Secondly, the assumption of a formal series solution which leads to a system of ordinary differential equations for the unknown coefficients in the series. The method generalizes to multi-phase and heterogeneous moving boundary problems for both constant temperature and Newton\u27s radiation conditions and yields simple and highly accurate estimates for both the temperature and boundary motion