Analytical and numerical solutions of inward spherical solidification of a superheated melt with radiative-convective heat transfer and density jump at freezing front

Analytical and numerical solutions of a general problem related to the radially symmetric inward spherical solidification of a superheated melt have been studied in this paper. In the radiation-convection type boundary conditions, the heat transfer coefficient has been taken as time dependent which could be infinite, at time,t=0. This is necessary, for the initiation of instantaneous solidification of superheated melt, over its surface. The analytical solution consists of employing suitable fictitious initial temperatures and fictitious extensions of the original region occupied by the melt. The numerical solution consists of finite difference scheme in which the grid points move with the freezing front. The numerical scheme can handle with ease the density changes in the solid and liquid states and the shrinkage or expansions of volumes due to density changes. In the numerical results, obtained for the moving boundary and temperatures, the effects of several parameters such as latent heat, Boltzmann constant, density ratios, heat transfer coefficients, etc. have been shown. The correctness of numerical results has also been checked by satisfying the integral heat balance at every timestep.

Glass Transition in the Density Functional Theory of Freezing

A mean-field description of the glass transition in the hard-sphere system is obtained by numerically locating "glassy" minima of a model free-energy functional. These minima, characterized by inhomogeneous but aperiodic density distributions, appear as the average density is increased above the value at which equilibrium crystallization takes place. Investigations of the density distribution and local bond-orientational order at these minima yield results similar to those obtained from simulations.

Effect of fluctuations on the freezing of a colloidal suspension in an external periodic potential

We incorporate the effects of fluctuations in a density functional analysis of the freezing of a colloidal liquid in the presence of an external potential generated by interfering laser beams. A mean-field treatment, using a density functional theory, predicts that with the increase in the strength of the modulating potential, the freezing transition changes from a first order to a continuous one via a tricritical point for a suitable choice of the modulating wavevectors. We demonstrate here that the continuous nature of the freezing transition at large values of the external potential V-e survives the presence of fluctuations. We also show that fluctuations tend to stabilize the liquid phase in the large V-e regime.

Freezing of classical fluids in a quenched random potential

The phase diagram of a hard-sphere fluid in the presence of a random pinning potential is studied analytically and numerically. In the analytic work, replicas are introduced for averaging over the quenched disorder, and the hypernetted chain approximation is used to calculate density correlations in the replicated liquid. The freezing transition of the liquid into a nearly crystalline state is studied using a density-functional approach, and the liquid to glass transition is studied using a phenomenological replica symmetry breaking approach. In the numerical work, local minima of a discretized version of the Ramakrishnan-Yussouff free-energy functional are located and the phase diagram in the density-disorder plane is obtained from an analysis of the relative stability of these minima. Both approaches lead to similar results for the phase diagram. The first-order liquid to crystalline solid transition is found to change to a continuous liquid to glass transition as the strength of the disorder is increased above a threshold value.

Melting-freezing cycles in a relatively sheared pair of crystalline monolayers

The nonequilibrium dynamical behaviour that arises when two ordered two-dimensional monolayers of particles are sheared over each other is studied in Brownian dynamics simulations. A curious sequence of nonequilibrium states is observed as the driving rate is increased, the most striking of which is a sliding state with irregular alternation between disordered and ordered states. We comment on possible mechanisms underlying these cycles, and experiments that could observe them.

Analytical solutions for heat transfer during cyclic melting and freezing of a phase change material used in electronic or electrical packaging

In this paper we develop an analytical heat transfer model, which is capable of analyzing cyclic melting and solidification processes of a phase change material used in the context of electronics cooling systems. The model is essentially based on conduction heat transfer, with treatments for convection and radiation embedded inside. The whole solution domain is first divided into two main sub-domains, namely, the melting sub-domain and the solidification sub-domain. Each sub-domain is then analyzed for a number of temporal regimes. Accordingly, analytical solutions for temperature distribution within each subdomain are formulated either using a semi-infinity consideration, or employing a method of quasi-steady state, depending on the applicability. The solution modules are subsequently united, leading to a closed-form solution for the entire problem. The analytical solutions are then compared with experimental and numerical solutions for a benchmark problem quoted in the literature, and excellent agreements can be observed.

Spin Probe Signature of Freezing in Water Slow Cooling Experiment: Our Observation is Global or Local?

First systematic spin probe ESR study of water freezing has been conducted using TEMPOL and TEMPO as the probes. The spin probe signature of the water freezing has been described in terms of the collapse of narrow triplet spectrum into a single broad line. This spin probe signature of freezing has been observed at an anomalously low temperature when a milimoler solution of TEMPOL is slowly cooled from room temperature. A systematic observation has revealed a spin probe concentration dependence of these freezing and respective melting points. These results can be explained in terms of localization of spin probe and liquid water,most probably in the interstices of ice grains, in an ice matrix. The lowering of spin probe freezing point, along with the secondary evidences, like spin probe concentration dependence of peak-to-peak width in frozen limit signal, indicates a possible size dependence of these localizations/entrapments with spin probe concentration. A weak concentration dependence of spin probe assisted freezing and melting points, which has been observed for TEMPO in comparison to TEMPOL, indicates different natures of interactions with water of these two probes. This view is also supported by the relaxation behavior of the two probes.

Self-diffusion of colloids at freezing

This is an introduction to the theory of interacting Brownian particles, as applied to charge-stabilised colloidal suspensions near their equilibrium liquid-solid transition. The density functional approach to the statics of the transition is reviewed briefly, and the generalised Langevin equation method for the dynamics presented in detail. Work with A.V. Indrani [1] on a self-consistent approach for calculating the excess single-particle friction is presented, which explains the observed [2] ''universal'' suppression of self-diffusion at freezing as a consequence of the universal structure-factor height at this transition. Criticisms, open questions, and challenges to theory are discussed.

Freezing of the octahedral tilt near ferromagnetic transition and appearance of a glassy phase at low temperature driven by the tilt instabilities in SrRuO3

The emergence of low temperature glassy phase in widely known itinerant ferromagnet SrRuO3 is remotely understood. In order to understand this aspect, we have undertaken a detailed temperature dependent (5-250 K) neutron diffraction study. We observe a freezing of the octahedral tilt near the ferromagnetic transition and an unusual deviation in the octahedral tilt near the onset of low temperature spin glass like phase. A reduction of the ordered magnetic moment and a decline in the total integrated magnetic intensity is observed around the temperature where the glassy behaviour starts to appear. The magnetotransport study also reveals the possibility for an additional magnetic ordering by demonstrating a peak in magnetoresistance at the low temperature side as well. The neutron diffraction study presented here provides useful information to understand the observed unusual low temperature magnetic phenomena in SrRuO3.

Reversible shear-induced crystallization above equilibrium freezing temperature in a lyotropic surfactant system

We demonstrate a unique shear-induced crystallization phenomenon above the equilibrium freezing temperature (T-K(o)) in weakly swollen isotropic (L-i) and lamellar (L-alpha) mesophases with bilayers formed in a cationic-anionic mixed surfactant system. Synchrotron rheological X-ray diffraction study reveals the crystallization transition to be reversible under shear (i.e., on stopping the shear, the nonequilibrium crystalline phase L-c melts back to the equilibrium mesophase). This is different from the shear-driven crystallization below T-K(o), which is irreversible. Rheological optical observations show that the growth of the crystalline phase occurs through a preordering of the L-i phase to an L-alpha phase induced by shear flow, before the nucleation of the Lc phase. Shear diagram of the L-i phase constructed in the parameter space of shear rate ((gamma)) over dot vs. temperature exhibits L-i -> L-c and L-i -> L-alpha transitions above the equilibrium crystallization temperature (T-K(o)), in addition to the irreversible shear-driven nucleation of L-c in the L-i phase below T-K(o). In addition to revealing a unique class of nonequilibrium phase transition, the present study urges a unique approach toward understanding shear-induced phenomena in concentrated mesophases of mixed amphiphilic systems.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.