Mathematical modelling of LiFePO4 cathodes

LiFePO4 is a commercially available battery material with good theoretical discharge capacity, excellent cycle life and increased safety compared with competing Li-ion chemistries. It has been the focus of considerable experimental and theoretical scrutiny in the past decade, resulting in LiFePO4 cathodes that perform well at high discharge rates. This scrutiny has raised several questions about the behaviour of LiFePO4 material during charge and discharge. In contrast to many other battery chemistries that intercalate homogeneously, LiFePO4 can phase-separate into highly and lowly lithiated phases, with intercalation proceeding by advancing an interface between these two phases. The main objective of this thesis is to construct mathematical models of LiFePO4 cathodes that can be validated against experimental discharge curves. This is in an attempt to understand some of the multi-scale dynamics of LiFePO4 cathodes that can be difficult to determine experimentally. The first section of this thesis constructs a three-scale mathematical model of LiFePO4 cathodes that uses a simple Stefan problem (which has been used previously in the literature) to describe the assumed phase-change. LiFePO4 crystals have been observed agglomerating in cathodes to form a porous collection of crystals and this morphology motivates the use of three size-scales in the model. The multi-scale model developed validates well against experimental data and this validated model is then used to examine the role of manufacturing parameters (including the agglomerate radius) on battery performance. The remainder of the thesis is concerned with investigating phase-field models as a replacement for the aforementioned Stefan problem. Phase-field models have recently been used in LiFePO4 and are a far more accurate representation of experimentally observed crystal-scale behaviour. They are based around the Cahn-Hilliard-reaction (CHR) IBVP, a fourth-order PDE with electrochemical (flux) boundary conditions that is very stiff and possesses multiple time and space scales. Numerical solutions to the CHR IBVP can be difficult to compute and hence a least-squares based Finite Volume Method (FVM) is developed for discretising both the full CHR IBVP and the more traditional Cahn-Hilliard IBVP. Phase-field models are subject to two main physicality constraints and the numerical scheme presented performs well under these constraints. This least-squares based FVM is then used to simulate the discharge of individual crystals of LiFePO4 in two dimensions. This discharge is subject to isotropic Li+ diffusion, based on experimental evidence that suggests the normally orthotropic transport of Li+ in LiFePO4 may become more isotropic in the presence of lattice defects. Numerical investigation shows that two-dimensional Li+ transport results in crystals that phase-separate, even at very high discharge rates. This is very different from results shown in the literature, where phase-separation in LiFePO4 crystals is suppressed during discharge with orthotropic Li+ transport. Finally, the three-scale cathodic model used at the beginning of the thesis is modified to simulate modern, high-rate LiFePO4 cathodes. High-rate cathodes typically do not contain (large) agglomerates and therefore a two-scale model is developed. The Stefan problem used previously is also replaced with the phase-field models examined in earlier chapters. The results from this model are then compared with experimental data and fit poorly, though a significant parameter regime could not be investigated numerically. Many-particle effects however, are evident in the simulated discharges, which match the conclusions of recent literature. These effects result in crystals that are subject to local currents very different from the discharge rate applied to the cathode, which impacts the phase-separating behaviour of the crystals and raises questions about the validity of using cathodic-scale experimental measurements in order to determine crystal-scale behaviour.

Identification and characterization of flowering genes in kiwifruit : sequence conservation and role in kiwifruit flower development

Background Flower development in kiwifruit (Actinidia spp.) is initiated in the first growing season, when undifferentiated primordia are established in latent shoot buds. These primordia can differentiate into flowers in the second growing season, after the winter dormancy period and upon accumulation of adequate winter chilling. Kiwifruit is an important horticultural crop, yet little is known about the molecular regulation of flower development. Results To study kiwifruit flower development, nine MADS-box genes were identified and functionally characterized. Protein sequence alignment, phenotypes obtained upon overexpression in Arabidopsis and expression patterns suggest that the identified genes are required for floral meristem and floral organ specification. Their role during budbreak and flower development was studied. A spontaneous kiwifruit mutant was utilized to correlate the extended expression domains of these flowering genes with abnormal floral development. Conclusions This study provides a description of flower development in kiwifruit at the molecular level. It has identified markers for flower development, and candidates for manipulation of kiwifruit growth, phase change and time of flowering. The expression in normal and aberrant flowers provided a model for kiwifruit flower development.

Azimuthal surface waves in a coaxial semiconductor structure with metal walls

The dispersion properties and topography of the fields of azimuthal surface wave (ASW) in a coaxial semiconductor structure with metal walls, placed in an external magnetic field, are investigated analytically and numerically. It is shown that an ASW phase-shifting device can be realized in the proposed structure. The conditions are indicated for which wave perturbations exist having frequencies that depend on the direction of phase change.

Multiphase porous media model for heat and mass transfer during drying of agricultural products

Modelling of food processing is complex because it involves sophisticated material and transport phenomena. Most of the agricultural products such fruits and vegetables are hygroscopic porous media containing free water, bound water, gas and solid matrix. Considering all phase in modelling is still not developed. In this article, a comprehensive porous media model for drying has been developed considering bound water, free water separately, as well as water vapour and air. Free water transport was considered as diffusion, pressure driven and evaporation. Bound water assumed to be converted to free water due to concentration difference and also can diffuse. Binary diffusion between water vapour and air was considered. Since, the model is fundamental physics based it can be applied to any drying applications and other food processing where heat and mass transfer takes place in porous media with significant evaporation and other phase change.

In situ investigation of explosive crystallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy

The crystallization of amorphous semiconductors is a strongly exothermic process. Once initiated the release of latent heat can be sufficient to drive a self-sustaining crystallization front through the material in a manner that has been described as explosive. Here, we perform a quantitative in situ study of explosive crystallization in amorphous germanium using dynamic transmission electron microscopy. Direct observations of the speed of the explosive crystallization front as it evolves along a laser-imprinted temperature gradient are used to experimentally determine the complete interface response function (i.e., the temperature-dependent front propagation speed) for this process, which reaches a peak of 16 m/s. Fitting to the Frenkel-Wilson kinetic law demonstrates that the diffusivity of the material locally/immediately in advance of the explosive crystallization front is inconsistent with those of a liquid phase. This result suggests a modification to the liquid-mediated mechanism commonly used to describe this process that replaces the phase change at the leading amorphous-liquid interface with a change in bonding character (from covalent to metallic) occurring in the hot amorphous material.

Effects of the filling process on the evolution of the mushy zone and macrosegregation in alloy casting

A numerical model of the entire casting process starting from the mould filling stage to complete solidification is presented. The model takes into consideration any phase change taking place during the filling process. A volume of fluid method is used for tracking the metal–air interface during filling and an enthalpy based macro-scale solidification model is used for the phase change process. The model is demonstrated for the case of filling and solidification of Pb–15 wt%Sn alloy in a side-cooled two-dimensional rectangular cavity, and the resulting evolution of a mushy region and macrosegregation are studied. The effects of process parameters related to filling, namely degree of melt superheat and filling velocity on macrosegregation in the cavity, are also investigated. Results show significant differences in the progress of the mushy zone and macrosegregation pattern between this analysis and conventional analysis without the filling effect.

Spin-state equilibria in the LCoO3 system from 59Co NMR

Spin-state equilibria in the whole set of LCoO3 (where L stands for a rare-earth metal or Y) have been investigated with the use of 59Co NMR as a probe for the polycrystalline samples (except Ce) in the temperature interval 110-550 K and frequency range 3- 11.6 MHz. Besides confirming the coexistence of the high-spin—low-spin state in this temperature range, a quadrupolar interaction of ∼0.1 -0.5 MHz has been detected for the first time from 59Co NMR. The NMR line shape is found to depend strongly on the relative magnitude of the magnetic and quadrupolar interactions present. Analysis of the powder pattern reveals two basically different types of transferred hyperfine interaction between the lighter and heavier members of the rare-earth series. The first three members of the lighter rare-earth metals La, Pr (rhombohedral), and Nd (tetragonal), exhibit second-order quadrupolar interaction with a zero-asymmetry parameter at lower temperatures. Above a critical temperature TS (dependent on the size of the rare-earth ion), the quadrupolar interaction becomes temperature dependent and eventually gives rise to a first-order interaction thus indicating a possible second-order phase change. Sm and Eu (orthorhombic) exhibit also a second-order quadrupolar interaction with a nonzero asymmetry parameter ((η∼0.47)) at 300 K, while the orthorhombic second-half members (Dy,..., Lu and Y) exhibit first-order quadrupolar interaction at all temperatures. Normal paramagnetic behavior, i.e., a linear variation of Kiso with T-1, has been observed in the heavier rare-earth cobaltites (Er,..., Lu and Y), whereas an anomalous variation has been observed in (La,..., Nd)CoO3. Thus, Kiso increases with increasing temperature in PrCoO3 and NdCoO3. These observations corroborate the model of the spin-state equilibria in LCoO3 originally proposed by Raccah and Goodenough. A high-spin—low-spin ratio, r=1, can be stabilized in the perovskite structure by a cooperative displacement of the oxygen atoms from the high-spin towards the low-spin cation. Where this ordering into high- and low-spin sublattices occurs at r=1, one can anticipate equivalent displacement of all near-neighbor oxygen atoms towards a low-spin cobalt ion. Thus the heavier LCoO3 exhibits a small temperature-independent first-order quadrupolar interaction. Where r<1, the high- and low-spin states are disordered, giving rise to a temperature-dependent second-order quadrupolar interaction with an anomalous Kiso for the lighter LCoO3.

Structural origin of set-reset processes in Ge15Te83Si2 glass investigated using in situ Raman scattering and transmission electron microscopy

We report here that the structural origin of an easily reversible Ge15Te83Si2 glass can be a promising candidate for phase change random access memories. In situ Raman scattering studies on Ge15Te83Si2 sample, undertaken during the amorphous set and reset processes, indicate that the degree of disorder in the glass is reduced from off to set state. It is also found that the local structure of the sample under reset condition is similar to that in the amorphous off state. Electron microscopic studies on switched samples indicate the formation of nanometric sized particles of c-SiTe2 structure. ©2009 American Institute of Physics

Studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring

Numerical and experimental studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring are performed. The alloy is electromagnetically stirred to produce semisolid slurry in a cylindrical graphite mould placed in the annulus of a linear electromagnetic stirrer. The mould is cooled at the bottom, such that solidification progresses from the bottom to the top of the cylindrical mould. A numerical model is developed for simulating the transport phenomena associated with the solidification process using a set of single-phase governing equations of mass. momentum, energy. and species conservation. The viscosity variation of the slurry, used in the model, is determined experimentally using a rotary viscometer. The set of governing equations is solved using a pressure-based finite volume technique, along with an enthalpy based phase change algorithm. The numerical study involves prediction of temperature, velocity, species and solid fraction distribution in the mould. Corresponding solidification experiments are performed, with time-temperature history recorded at key locations. The microstructures at various temperature measurement locations in the solidified billet are analyzed. The numerical predictions of temperature variations are in good agreement with experiments, and the predicted flow field evolution correlates well with the microstructures observed at various locations.

Effects of mould filling on evolution of the solid-liquid interface during solidification

This work presents a numerical analysis of simultaneous mould filling and phase change for solidification in a two-dimensional rectangular cavity. The role of residual flow strength and temperature gradients within the solidifying domain, caused by the filling process, on the evolution of solidification interface are investigated. An implicit volume of fluid (VOF)-based algorithm has been employed for simulating the free surface flows during the filling process, while the model for solidification is based on a fixed-grid enthalpy-based control volume approach. Solidification modeling is coupled with VOF through User Defined Functions developed in the commercial computational fluid dynamics (CFD) code FLUENT 6.3.26. Comparison between results of the conventional analysis without filling effect and those of the present analysis shows that the residual flow resulting from the filling process significantly influences the progress of the solidification interface. A parametric study is also performed with variables such as cooling rate, filling velocity and filling configuration, in order to investigate the coupled effects of the buoyancy-driven flow and the residual flow on the solidification behavior.

Design and analysis of ARROW-based racetrack resonator vibration sensor for inertial navigation application

In this paper we propose and analyze a novel racetrack resonator based vibration sensor for inertial grade application. The resonator is formed with an Anti Resonance Reflecting Optical Waveguide (ARROW) structure which offers the advantage of low loss and single mode propagation. The waveguide is designed to operate at 1310nm and TM mode of propagation since the Photo-elastic co-efficient is larger than TE mode in a SiO2/ Si3N4/ SiO2. The longer side of the resonator is placed over a cantilever beam with a proof mass. A single bus waveguide is coupled to the resonator structure. When the beam vibrates the resonator arm at the foot of the cantilever experiences maximum stress. Due to opto-mechanical coupling the effective refractive index of the resonator changes hence the resonance wavelength shifts. The non uniform cantilever beam has a dimension of 1.75mm X 0.45mm X 0.020mm and the proof mass has a dimension of 3mm X 3mm X 0.380mm. The proof mass lowers the natural frequency of vibration to 410Hz, hence designed for inertial navigation application. The operating band of frequency is from DC to 100Hz and acceleration of less than 1g. The resonator has a Free Spectral Range (FSR) of 893pm and produces a phase change of 22.4mrad/g.

Effect of process parameters on transport phenomena during solidification in the presence of electromagnetic stirring

In the present work, a numerical study is performed to predict the effect of process parameters on transport phenomena during solidification of aluminium alloy A356 in the presence of electromagnetic stirring. A set of single-phase governing equations of mass, momentum, energy and species conservation is used to represent the solidification process and the associated fluid flow, heat and mass transfer. In the model, the electromagnetic forces are incorporated using an analytical solution of Maxwell equation in the momentum conservation equations and the slurry rheology during solidification is represented using an experimentally determined variable viscosity function. Finally, the set of governing equations is solved for various process conditions using a pressure based finite volume technique, along with an enthalpy based phase change algorithm. In present work, the effect of stirring intensity and cooling rate are considered. It is found that increasing stirring intensity results in increase of slurry velocity and corresponding increase in the fraction of solid in the slurry. In addition, the increasing stirring intensity results uniform distribution of species and fraction of solid in the slurry. It is also found from the simulation that the distribution of solid fraction and species is dependent on cooling rate conditions. At low cooling rate, the fragmentation of dendrites from the solid/liquid interface is more.

Experimental approach for studying structural changes in axonal membrane upon nerve excitation

The Hodgkin and Huxley (HH) model of action potential has become a central paradigm of neuroscience. Despite its ability to predict action potentials with remarkable accuracy, it fails to explain several biophysical findings related to the initiation and propagation of the nerve impulse. The isentropic heat release and optical phenomena demonstrated by various experiments suggest that action potential is accompanied by a transient phase change in the axonal membrane. In this study a method was developed for preparing a giant axon from the crayfish abdominal cord for studying the molecular mechanisms of action potential simultaneously by electrophysiological and optical methods. Also an alternative setup using a single-cell culture of an Aplysia sensory neuron is presented. In addition to the description of the method, the preliminary results on the effect of phloretin, a dipole potential lowering compound, on the excitability of a crayfish giant axon are presented.

Heat transfer correlations for PCM-based heat sinks with plate fins

Characterization of melting process in a Phase Change Material (PCM)-based heat sink with plate fin type thermal conductivity enhancers (TCEs) is numerically studied in this paper. Detailed parametric investigations are performed to find the effect of aspect ratio of enclosure and the applied heat flux on the thermal performance of the heat sinks. Various non-dimensional numbers, such as Nusselt number (Nu), Rayleigh number (Ra), Stefan number (Ste) and Fourier number (Fo) based on a characteristic length scale, are identified as important parameters. The half fin thickness and the fin height are varied to obtain a wide range of aspect ratios of an enclosure. It is found that a single correlation of Nu with Ra is not applicable for all aspect ratios of enclosure with melt convection taken into account. To find appropriate length scales, enclosures with different aspect ratios are divided into three categories, viz. (a) shallow enclosure, (b) rectangular enclosure and (c) tall enclosure. Accordingly, an appropriate characteristic length scale is identified for each type of enclosure and correlation of Nu with Ra based on that characteristic length scale is developed. (C) 2010 Elsevier Ltd. All rights reserved.

Thin film amorphous semiconductor microwave switches

Microwave switches operating in the X band were designed and fabricated using amorphous chalcogenide semiconductors of composition GexTeyAsz. Threshold devices were shown to operate as microwave modulators at modulation frequencies of up to 100 MHz. No delay time was observed at the highest frequency although the modulation efficiency decreased above 10 MHz owing to the finite recovery time which was approximately 0.3 × 10−8s. The devices can also be used as variolossers, the insertion loss being 0.5 dB in the OFF state and increasing on switching from 5 dB at 1 mA device current to 18 dB at 100 mA.The behaviour of the threshold switches can be explained in terms of the formation of a conducting filament in the ON state with a constant current density of 2 × 104Acm−2 that is shunted by the device capacitance. The OFF state conductivity σ varies as ωn (0.5 < n < 1) which is characteristic of hopping in localized states. However, there was evidence of a decrease in n or a saturation of the conductivity at high frequencies.As a result of phase separation memory switches require no holding current in the ON state and may be used as novel latching semiconductor phase-shifters.